Interaction Between Apolipoprotein M Gene Single-Nucleotide Polymorphisms and Obesity and its Effect on Type 2 Diabetes Mellitus Susceptibility

This study investigated the correlation of four single nucleotide polymorphisms (SNPs) in Apolipoprotein M (ApoM) with the risk of type 2 diabetes mellitus (T2DM) and effects of the interactions of this gene and obesity. The effects of SNP and obesity interaction on T2DM was examined by generalized multifactor dimensionality reduction (GMDR) combined with the logistic regression model. T2DM patient-control haplotype was analyzed in silico using the haplotype analysis algorithm SHEsis. The rs805296-C allele or 724-del allele indicted high risk of T2DM. The incidence of T2DM in individuals with rs805296-C allele polymorphism (TC + CC) was higher than those without (TT), adjusted OR (95%CI) = 1.29 (1.10–1.66) (p < 0.001). Moreover, the individuals with 724-delallele have a higher risk of T2DM compared to those with 724-ins variants, adjusted OR (95%CI) = 1.66 (1.40–2.06), p < 0.001. GMDR analysis suggested that the interaction model composed of the two factors, rs805296 and obesity, was the best model with statistical significance (P value from sign test [Psign]=0.0107). The T2DM risk in obese individuals having TC or CC genotype was higher than non-obese individuals with TT genotype (OR = 2.38, 95% CI = 1.58–3.53). Haplotype analysis suggests that rs805297-C and rs9404941-C alleles haplotype indicate high risk of T2DM, OR (95%CI) = 1.62 (1.29–2.16), p < 0.001. Our results suggested that rs805296 and 724-del minor allele of ApoM gene, interaction of rs805296 and obesity, rs805297-C and rs9404941-C alleles haplotype were indicators of high T2DM risk.

Body measurements. General clinicopathological information of enrolled participants were recorded by trained staff. Body weight in kilograms divided by the square of the height in meters was used to measure body mass index (BMI). Individuals with cigarette smoking history for one year or above and smoke at least once per day were classified as cigarette smokers. Alcohol amount consumed was the sum alcohol consumed per week from all kinds of wine. After 8 or more hours of fasting, blood samples were obtainedfrom the individuals in the next morning. Samples were stored at -80 °C until use. Oxidase enzymatic method was used to measure plasma glucose concentration. High-density lipoprotein (HDL)-cholesterol and triglyceride (TG) concentrations were measured using a Hitachi biochemistry.
Genomic DNA extraction and genotyping. dbSNP algorithm (http://www.ncbi.nlm.nih.gov/projects/ SNP)was used to select the SNPs to be investigated (criteria: MAF > 5% in the dbSNP database). Four SNPs (rs805296,−724 ins/del, rs805297 and rs9404941) were selected for further investigations in this study. The genomic DNA was extracted from the collected patients or healthy individuals blood samples using Genomic DNA extraction kit (Roche, USA) and stored at −80 °C. The genotype of the selected SNPs in the samples was analyzed by Restriction Fragment Length Polymorphism (RFLP) on the basis of polymerase Chain Reaction (PCR). The primer sequences and descriptions of 4 SNPs are shown in Table 1. The PCR reaction system included: Taq DNA polymerase, dNTPs, PCR buffer, and MgCl 2 . The measurement of PCR detection reagent is as follows: less than 0.1 μg genomic DNA template, 12.5 μl 2 × Taq PCR Mastermix, 10 μmol of each primer and add ddH 2 O to a final reaction volume of 25 μl. PCR was carried out at an Applied Biosystems PCR equipment using the following procedures: 1 cycle of 94 °C denaturation for 3 min, 30 cycles of 95 °C denaturation for 30 s, 60 °C annealing for 30 s, and 72 °C extension for 30 s. The resulted products were sequenced using an automatic sequencer (Model 3730, BGI, Shanghai, China). Diagnostic criteria. Individuals with a fasting glucose ≥ 126 mg/dl (7.0 mmol/l) or having undergone hypoglycemic therapy were diagnosed as diabetics. Individuals with a fasting glucose ≥ 126 mg/dl (7.0 mmol/l), or blood glucoselevels2 h postprandial ≥ 200 mg/dl (11.0 mmol/l), or having undergone hypoglycemic therapy in the interim were classified into T2DM group 18 . Individuals with a BMI ≥ 28 kg/m 2 were classified as obese 19 .
Statistical analysis. All data analysis was performedon SPSS 22.0 software (Chicago, IL). Mean and standard deviation (SD) were measured for continuous variables and the differences were analyzed using Students' t-test; percentages were measured for categorical variables and the differences were analyzed using chi-square test. The genotype distribution differences among individuals with T2DM and healthy controls were analyzed usingChi-square test. In silico analysis algorithm SHEsis was used toanalyze T2DM patient-control haplotype (http://analysis.bio-x.cn/myAnalysis.php). Generalized multifactor dimensionality reduction (GMDR) was performed to investigate all the interactions. Effects of SNPs and obesity interaction on the risk of T2DM were measured by logistic regression model. The collected clinicopathologicalinformation was used to adjust odds.

Results
A total of 1371 participants, consisting of 681 in the T2DM group and 690 as healthy controls, were enrolled (606 males and 765 females) with a mean age at 61.1 ± 13.8. The general characteristics of these enrolled T2DM patients and healthy controls are shown in Table 2. We observed a significant distribution difference between T2DM patients andhealthy controls in BMI, TG, TC, and HDL. However, no close association was observed for males, smoking, alcohol consumption, and mean age between T2DM patients and controls. Genotype distribution was analyzed using the Hardy-Weinberg equilibrium. The frequenciesof C allele of rs805296and 724-delwere higher in individuals with T2DM compared to healthy controls (30.4% of T2DM patients and 22.2% of controls, p < 0.001 for C allele of rs805296; 28.9% of T2DM patients and 21.2% of controls, p < 0.001 for 724-del) ( Table 3). Logistic regression model revealedthe incidence of T2DM in individuals with rs805296-C allele or 724-del allelewas higher compared to those with TT variants or 724-ins variants respectively (Table 3). However, no significant correlationswere found when the association of C-1065Ars805297 and T-855Crs9404941 with T2DM risk (Table 3) were analyzed.
After covariates adjustment, GMDR analysis was performed to analyze the correlation of ApoM gene and obesity interaction with the risk of T2DM (Table 4). The results revealed that the interaction model composed of the two factors, rs805296 and obesity, which was the best model with statistical significance (P value from sign test [P sign ]=0.0107). Meanwhile, the cross-validation consistency and testing accuracy for this two-locus were 10/10 and 62.17%, respectively. After adjustment for the collected clinicopathological information, we found the incidence of T2DM in individuals with TC or CC genotype and high BMI was higher than those with TT genotype and normal BMI (Table 5).
D' values of4 ApoM gene SNPs were analyzed using Pairwise LD statistics. The results presented in Table 6 showed that value calculated for rs9404941 and rs805297was 0.817. Therefore, analysis for rs9404941 and rs805297 was performed within silico haplotype analysis softwareSHEsis. As a result, the frequency of T-T haplotype was the highest in both populations (47.01% for individuals with T2DM, 54.67% for healthy controls). Also, our results demonstrated that rs805297-C and rs9404941-C alleles were indicators for a high T2DM risk (Table 7).

Discussion
Our research indicated that T2DM risk was positively correlated with rs805296-C or 724-del allele but did not have any association with the other two SNPs. Although different studies have identified several potential candidate genes associated with the risk of T2DM risk factors, including SNPs polymorphism within ApoM gene [12][13][14][15] . However, the conclusions drawn from these reports regarding the relationship between ApoM SNPs and incidence of T2DMare inconsistent. For example, two Chinese studies indicated that ApoM rs805296 (T-778C) was an indicator forthe risk of both type 1 and type 2 diabetes 12,13 . However, nocorrelationbetween this specific gene polymorphism and the risk of T2DM was found in an independent study 14 . Zhang et al. 20 analyzed the associations between four SNPs investigated in this study and the risk of T2DM in a total of 335 eastern Han Chinese participants. The resultsillustrated that C-724del polymorphism predicts a high risk of T2DM but rs805296 (T-778C) polymorphism did not correlate with the risk of T2DM. Moreover, the genotype frequency and distribution difference of rs805297 (C-1065A) and rs9404941 (T-855C) inindividuals with T2DM and healthy controls was not significant. Xu et al 11 . suggested that ApoM rs9404941 (T-855C) predicts highincidenceof CAD 8 . But another study 21 indicated that rs805297 SNP in ApoM gene predicts high risk of dyslipidemia but did not have any influence on the incidence of CAD. A previous study 22 also suggested that plasma ApoMexpressionwasupregulated in individuals with hyperlipidemia butdownregulated in individuals with T2DM compared with that in healthy controls. In addition, previous studiesimplied that ApoM expression was inversely correlated withglycemia 23 . Moreover, overexpression of ApoM in Goto-Kakizaki rats enhanced the effects of insulin 24 , indicating that ApoMhas the potential to be used a therapeutic target for T2DM.
It was reported that genetic background is a key modulator for the human reactions to environmental risk factors. The importance of environmental risk factors including obesity in the pathogenesis of T2DM has been widely recognized 15,16 . However, no study was performed to analyze the gene-environment interaction, especially the interaction ofApoM and obesity on the incidence of T2DM. Therefore, GMDR model was used to analyze the interaction of ApoMand obesity onthe incidence of T2DM in Chinese population;it revealed that the interaction www.nature.com/scientificreports www.nature.com/scientificreports/ of rs805296 and obesity was significant. We also found that TC or CC genotype and high BMI increased the incidence of T2DM in comparisonto TT genotype and normal BMI. In this study, we showed a string chain reaction between rs805297 and rs9404941 since the D' value was above 0.8. We also found that individuals with rs805297-C and rs9404941-C allelestend to have high risk of T2DM using haplotype analysis.
There were several shortcomings in this study. Firstly, limited numbers ofApoM SNPswere investigated and may result in the genetic information of ApoM was not sufficientlyfactored into the analysis. Secondly, the

Genotypes and Alleles
Frequencies n (%)     Table 7. ApoM gene and T2DM risk association measured by haplotype analysis. *Adjusted for gender, age, smoking and BMI.