- 1Laboratory of Genetic Disease and Perinatal Medicine and Key Laboratory of Birth Defects and Related Diseases of Women and Children of the Ministry of Education, West China Second University Hospital, Sichuan University, Chengdu, China
- 2West China School of Nursing, Sichuan University, Chengdu, China
- 3Department of Obstetrics and Gynecology, West China Second University Hospital, Sichuan University, Chengdu, China
- 4Division of Peptides Related With Human Disease, West China Hospital, Sichuan University, Chengdu, China
- 5Department of Biochemistry and Molecular Biology, West China School of Preclinical and Forensic Medicine, Sichuan University, Chengdu, China
Background: ATP-binding cassette transporter A1 (ABCA1) has important roles in high-density lipoprotein (HDL) metabolism and reverse cholesterol transport, and is implicated in lipid-related disorders. Genetic variants are involved in the pathogenesis of gestational diabetes mellitus (GDM). The objective of this study was to investigate the association of rs2230806 (R219K), a single nucleotide polymorphism (SNP) in the lipid-related gene, with the risk of GDM and related traits.
Methods: The SNP, rs2230806, was genotyped, and clinical and metabolic parameters were determined in 660 GDM patients and 1,097 control subjects. Genetic associations with related traits were also analyzed.
Results: The genotype distributions were similar in GDM patients and normal controls. However, significant differences in the variables examined in the study subjects were noted across the three genotypes. The genotype at the rs2230806 polymorphism was significantly associated with HDL-cholesterol (HDL-C) levels and atherogenic index (AI) values in GDM patients and total cholesterol (TC) and low-density lipoprotein-cholesterol (LDL-C) levels in control subjects. Subgroup analysis showed that the polymorphism was associated with diastolic blood pressure, in addition to HDL-C levels and AI, in overweight/obese GDM patients, while it was associated with TC levels, AI, pre-pregnancy body mass index (BMI), and BMI at delivery in non-obese GDM patients. In addition, this polymorphism was associated with TC, LDL-C, and apoB levels in overweight/obese control subjects.
Conclusions: The rs2230806 polymorphism in the ABCA1 gene was associated with variations in atherometabolic traits in GDM patients, with characteristics of BMI dependency, but not with GDM. Our findings highlight a link between related phenotypes in women with GDM and genetic factors.
Introduction
Gestational diabetes mellitus (GDM) is the most common metabolic disorder that occurs during pregnancy. It is defined as glucose intolerance beginning in the second or third trimester of pregnancy (1). The prevalence of GDM varies between 1% and 28% globally (2), and is especially prevalent in developed countries (3). The reported prevalence of GDM in pregnant Chinese women is also high, with a rate of 14.8% (4). The increased incidence of risk factors, such as obesity, poor dietary habits, physical inactivity, and older maternal age, has led to an increased occurrence of the disease globally in recent years (4, 5). GDM can be dangerous to mothers and infants in the perinatal period (5), and it is associated with an increased risk of the development of metabolic-related disorders, such as type 2 diabetes mellitus, metabolic syndrome, and cardiovascular diseases in the years following pregnancy (5, 6). The etiology of GDM may involve genetic factors (7), lipid abnormalities (8), oxidative stress (9), and inflammatory reactions (10).
The ABCA1 gene is located on the human chromosome 9q31 and it encodes an integral cell membrane protein that facilitates the efflux of cholesterol and phospholipids. ABCA1 is highly expressed in various human tissues, including the placenta, uterus, liver, adrenal glands, small intestine, lungs, and heart (11), and it is believed to be functionally relevant in these tissues. In addition, activated pancreatic beta cell ABCA1 would stimulate insulin production. One study demonstrated that GDM affects different key players of cholesterol transfer including ABCA1, in placenta depending on the maternal BMI (12). Another study reported that up-regulation of ABCA1 and ABCG1 mRNA and protein levels were observed in GDM versus control human fetoplacental endothelial cells (13). In addition, recent study evaluated effects of GDM and obesity on the intrauterine metabolic milieu in the fetus, showing exposure to maternal obesity in utero leads to sex-specific changes in miRNA and expression of several target genes, including ABCA1 in human fetal liver (14). Furthermore, one study showed that serum ABCA1 could be detected in women with normal pregnancy, and the patients with preeclampsia had significantly lower serum ABCA1 levels than those in the control subjects (15). A substantial number of studies have shown a new role for ABCA1 in the liver and a metabolic relationship between increased TG levels and decreased high-density lipoprotein (HDL) levels (16, 17). Moreover, recent studies have investigated the role of ABCA1 in other features of metabolic syndrome and demonstrated its relationship with body weight, decreased insulin secretion and sensitivity, and increased blood glucose levels (18–20). Changes in these metabolic components are also frequently encountered in pregnant women, especially in pregnant women with GDM. In addition, ABCA1 has been reported to be expressed at high levels in various compartments of the placenta, including the villous syncytiotrophoblasts and fetal endothelial cells (21). ABCA1 is also linked to the feto-maternal transport of cholesterol, indicating its crucial role in normal fetal development in humans. Several studies of the human placenta have shown that ABCA1 exerts important beneficial effects against cardiovascular and metabolic diseases, and thus, is still of major interest as a therapeutic target for the prevention of such disorders. Approaches to targeting ABCA1 include the upregulation of ABCA1 expression by LXR activators (22) and recently discovered specific microRNAs (23).
Since ABCA1 has demonstrated roles in various GDM-related traits, such as lipid levels, obesity-related phenotypes, and insulin secretion and sensitivity, it may be a strong candidate gene for predisposition to GDM. The ABCA1 gene is highly polymorphic. One of the most common single nucleotide polymorphisms (SNPs) is the rs2230806 SNP (also known as R219K or G1051A) in exon 7. This SNP results in the substitution of an arginine amino acid (G allele) with a lysine (A allele) amino acid at position 219 of the polypeptide. This amino acid substitution is functionally relevant to the role of ABCA1 in lipid metabolism. For example, the A allele of the SNP is associated with higher levels of HDL-cholesterol (HDL-C) in various populations, such as Turkish (24), European (25), Japanese (26), and Chinese (27) populations. In addition, the A allele of the SNP is reported to be significantly associated with lower triglyceride (TG) (24, 28), total cholesterol (TC) (29–31), and low-density lipoprotein cholesterol (LDL-C) (28, 32, 33) levels. Clinically, the allelic variant of this SNP has been found to be associated with a reduced incidence of coronary artery disease (CAD) events (25, 34), premature CAD in patients with familial hypercholesterolemia (35), and type 2 diabetes mellitus (36). We speculated that there may be a relationship between this ABCA1 SNP and the risk of GDM and/or related traits or phenotypes.
We tested the hypothesis that ABCA1 rs2230806 SNP genotype is associated with the risk of GDM and/or atherometabolic traits in GDM patients. To date, there is little information on the possible relationship between ABCA1 polymorphisms and GDM. In this study, we evaluated the distribution of the ABCA1 rs2230806 SNP genotypes in a relatively small sample size of pregnant women with (n = 660) or without (n = 1,097) GDM from a well-characterized Chinese population. We also investigated whether this SNP was associated with the risk of GDM and related traits in our study subjects.
Materials and Methods
Study Subjects
Blood samples were taken from 1,757 Chinese women from 2013–2020, including 660 patients (average age, 35.57 ± 4.03) with GDM and 1,097 control subjects (average age, 34.54 ± 4.99). All study subjects were of Chinese Han ethnicity. This study was conducted in accordance with the principles of the Declaration of Helsinki. Written informed consent was obtained from all participants. The study was approved by the Institutional Review Board of West China Second University Hospital, Sichuan University (No. 2017–033).
All patients with GDM were recruited from the Obstetrics Department of West China Second University Hospital at Sichuan University in Chengdu. Each patient with GDM met the diagnostic criteria recommended by the International Association of Diabetes and Pregnancy Study Groups (24-28 weeks of gestation, fasting blood glucose concentration ≥ 5.1 mmol/L by oral 75 g glucose tolerance test, blood glucose concentration ≥ 10.0 mmol/L 1 h after and ≥ 8.5 mmol/L 2 h after glucose consumption) (37). Control subjects with uncomplicated pregnancies were recruited from the same department of the hospital during the same period. Patients were included if they underwent elective cesarean section or planned vaginal delivery at term, and had a singleton pregnancy. None of the study subjects had clinically evident chronic or acute diseases, such as infections, tumors, thyroid dysfunction, or cardiovascular disease.
Body weight in light clothing and barefoot height were measured, and body mass index (BMI, kg/m2) was determined. Normal weight and overweight/obesity were defined as BMI < 25 and ≥ 25 kg/m2, respectively. Sitting blood pressure was measured on the right arm after a 10-min rest with a semi-automatic device (Omron HBP-1300). The infant’s length and weight were measured after birth. Blood samples of the study subjects were taken in a stage of late pregnancy after an overnight fast. Plasma aliquots were stored at -80°C.
DNA Extraction and Genotyping
Genomic DNA was isolated from 500 μL of peripheral blood according to the method described by Erlich (38). PCR was performed in a final volume of 25 μL containing 12.5 μL of 2× Taq PCR Master Mix, 10.5 μL of ddH2O, and 0.5 mM each upper- and lower-strand primer. We used 1.0 μL of DNA sample per PCR. The investigated DNA sequences were amplified using the forward primer, 5′-GTA TTT TTG CAA GGC TAC CAG TTA CAT TTG ACA A-3′ and the reverse primer, 5′-GAT TGG CTT CAG GAT GTC CAT GTT GGA A-3′, which produced a 174 bp product. Amplifications were performed using a MyCycler™ thermal cycler system (Bio-Rad, Hercules, CA, USA) under the following conditions: initial denaturation at 94°C for 4 min, followed by 32 amplification cycles of denaturation at 94°C for 50s, primer binding (annealing) at 58°C for 50s, and chain elongation at 72°C for 1 min. The PCR program ended with a 10-min chain elongation step at 72°C. PCR products were digested with 10 U/μL Styl-HF (New England Biolabs, Ipswich, MA, USA), analyzed by electrophoresis on a 2.5% agarose gel, and visualized by staining with Genecolour fluorescent dye. Enzyme digestion resulted in 109- and 65-bp fragments for the GG genotype; a non-digested 174-bp fragment for the AA genotype; and 174-,109-, and 65-bp fragments for the GA genotype. To document the results, gel images were captured under ultraviolet light.
Analysis of Metabolic Parameters
Plasma insulin levels were measured using chemiluminescence assays (IMMULITE 2000; Diagnostic Products Corporation, Los Angeles, CA, USA). The concentrations of TC, HDL-C, LDL-C, and TGs were measured in plasma samples using an enzyme assay (Boehringer, Mannheim, Germany), and serum apoA1 and apoB levels were determined by polyethylene glycol-enhanced immunoturbidimetric assay (Siemens Healthcare Diagnostics, Inc., Munich, Germany) using a Hitachi 7600-010 automatic analyzer (Hitachi, Tokyo, Japan). Plasma glucose levels were measured using the glucose oxidase technique (Roche Diagnostics, Mannheim, Germany). The intra- and inter-assay coefficients of variation for all measurements were less than 5% and 10%, respectively. The homeostatic model assessment of insulin resistance (HOMA-IR) was calculated as [fasting glucose concentration (mmol/L) × fasting insulin concentration (μU/mL)]/22.5. The atherogenic index (AI) was calculated using the following equation: AI = [(TC)-(HDL-C)]/[HDL-C] (39).
Statistical Analysis
Data are presented as the mean ± standard deviation. Allele frequencies of ABCA1 polymorphisms were estimated by gene counting. Hardy-Weinberg equilibrium was tested in cases and controls using a chi-square test. Allele and genotype frequencies were compared between cases and controls using chi-square analysis. Differences in variables between the GDM and control subjects were evaluated using independent-sample Student t-tests. To assess the effects of ABCA1 gene polymorphism and changes in quantitative variables, we performed an analysis of variance. Adjustments for age and BMI were performed using analysis of covariance. The effects of other parameters, including the rs2230806 G/A genotypes of ABCA1, glycemic status (GDM = 1, control = 0), pre-pregnancy BMI, and BMI at delivery, on insulin, TG, TC, HDL-C, and LDL-C levels were evaluated by multivariate stepwise regression analyses. Statistical significance was set at p < 0.05. All statistical analyses were performed using the Statistical Program for Social Sciences (SPSS) 21.0 (IBM, Armonk, NY, USA).
Results
Clinical and Biochemical Characteristics
The comparison of clinical characteristics and biochemical parameters between the patients with GDM and controls is shown in Table 1. Compared to the controls, prepregnancy BMI, basic systolic blood pressure (SBP), and diastolic blood pressure (DBP) were found to be significantly higher, and weight gain during pregnancy and neonatal birth height were significantly lower in the patients with GDM. Women with GDM also had significantly higher fasting insulin and glucose concentrations, HOMA-IR scores, and TG levels, and lower LDL-C levels, AI values, and apo A1 levels compared with the control group after adjusting for differences in age and BMI at delivery. Of the 660 women with GDM, 71 women required insulin treatment and 589 women had dietary and exercise-based treatment only (Table 1). These results were generally consistent with the clinical and metabolic profiles frequently observed in GDM subjects, as observed in our recent study (40).
ABCA1 Genotype and Allele Frequencies
Table 2 displayed the genotypic or allelic frequencies of SNP rs2230806 in both GDM patients and controls. All 1757 samples were successfully genotyped, and the frequencies of the SNP were in Hardy-Weinberg equilibrium in both the GDM patient and control groups. The G and A allele frequencies of the rs2230806 locus were 58.3% and 41.7%, respectively, in the GDM group and 55.9% and 44.1%, respectively, in the control group. These genotype and allele frequencies were not significantly different between groups (p > 0.05).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/782453/fendo-12-782453-HTML-r1/image_m/fendo-12-782453-t002.jpg)
Table 2 Genotype and allele frequency for the rs2230806 SNP of the ABCA1 gene in women with GDM and control subjects.
In addition, when GDM patients and control subjects were further divided into overweight/obese and non-obese subgroups, no significant differences in genotype or allele frequencies were observed between these subgroups (Table 3).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/782453/fendo-12-782453-HTML-r1/image_m/fendo-12-782453-t003.jpg)
Table 3 Genotype and allele frequency for the rs2230806 SNP of the ABCA1 gene in women with GDM and control subjects divided into overweight/obese and nonobese subgroups.
Effect of ABCA1 SNP on Clinical and Metabolic Parameters of the GDM and Control Groups
Associations of the SNP rs2230806 and clinical characteristics and metabolic traits are shown in Table 4. In the GDM patient group, AA homozygotes at the rs2230806 locus had lower HDL-C levels and higher AI values than GG homozygotes and GA heterozygotes (p < 0.05). These differences remained significant after adjusting for age and BMI. There were no significant changes in other parameters according to the genotype at the rs2230806 locus in GDM patients.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/782453/fendo-12-782453-HTML-r1/image_m/fendo-12-782453-t004.jpg)
Table 4 Clinical characteristics and metabolic profile of ABCA1 rs2230806 genotypes in GDM patients and controls.
In addition, in the control group, TC and LDL-C levels were lower in AA homozygotes than GG homozygotes (all p < 0.05, Table 4). There is an insignificant trend of genotype-related effect on HDL-C variation similar as that in GDM patients. There were no significant changes according to genotype for the other tested parameters in control subjects.
Effect of ABCA1 SNP on Clinical and Metabolic Parameters in Overweight/Obese and Nonobese Subgroups
Since several polymorphisms of the ABCA1 gene affect HDL-C levels in the general population, and low HDL-C levels were associated with high BMI (obese phenotype), we further analyzed the effects of the ABCA1 gene polymorphism on clinical and metabolic parameters in GDM patients and control subjects stratified into overweight/obese (BMI ≥ 25 kg/m2) and non-obese (BMI < 25 kg/m2) subgroups. In the overweight/obese GDM group, we found that AA homozygotes at the rs2230806 polymorphic site had higher DBP and AI values than heterozygotes and GG homozygotes, lower TC levels than heterozygotes, and lower HDL-C levels than heterozygotes and GG homozygotes (p < 0.05, Table 5). In addition, in the non-obese GDM group, the AA homozygotes had a lower pre-pregnancy BMI and delivery BMI, and higher TC and AI than GG homozygotes or GG homozygotes and GA heterozygotes (Table 5).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/782453/fendo-12-782453-HTML-r1/image_m/fendo-12-782453-t005.jpg)
Table 5 Clinical characteristics and metabolic profile of ABCA1 rs2230806 genotypes in overweight/obese and non-obese GDM patients.
In the overweight/obese control group, TC, LDL-C, and apoB levels were lower in the AA homozygotes than the GG homozygotes (all p < 0.05). Genotype-related effects at this locus were not observed in non-obese control subjects (Table 6).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/782453/fendo-12-782453-HTML-r1/image_m/fendo-12-782453-t006.jpg)
Table 6 Clinical characteristics and metabolic profile of ABCA1 rs2230806 genotypes in overweight/obese and non-obese control subjects.
Association of Important Independent Variables With Insulin and Lipid Levels in Women With and Without GDM
Based on multivariate stepwise regression analysis, Table 7 showed that maternal BMI was a significant predictor of fasting insulin levels. In addition, prepregnancy BMI and glycemic status were significant predictors of TG and TC levels, prepregnancy BMI and allelic variant were significant predictors of HDL-C levels, and prepregnancy BMI, glycemic status, and allelic variant were significant predictors of LDL-C levels in women with and without GDM (Table 7).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/782453/fendo-12-782453-HTML-r1/image_m/fendo-12-782453-t007.jpg)
Table 7 Association of important independent predictors/confounders with biomarker levels by multivariate regression analysis in women with and without GDM.
Discussion
In this study we used candidate gene approach to explore the association of a lipid-related gene rs2230806 SNP (R219K) with the risk of GDM and related traits. Our results showed, for the first time, that the polymorphism was significantly associated with HDL-C levels and AI values in GDM patients. Further subgroup analysis showed that the polymorphism was associated with DBP in addition to HDL-C levels and AI values in overweight/obese GDM patients, while it was associated with TC levels, AI values, prepregnancy BMI, and BMI at delivery in non-obese GDM patients. These results provide evidence that the polymorphism is associated with variations in atherometabolic traits in GDM patients, with characteristics of BMI dependency. These findings highlight a link between related phenotypes in women with GDM and the genetic component.
The rs2230806 polymorphism controls the substitution of an arginine at residue 219 of the peptide with a lysine. This substitution may affect the function of ABCA1; that is, the first step in reverse cholesterol transport (41). Our results in a Chinese cohort showed that the rs2230806 polymorphism was associated with HDL-C concentration in GDM patients, especially those with higher BMI, supporting the notion that ABCA1 plays a central role in reverse cholesterol transport as a modulator of HDL-C. In addition, the same trend of the genotype-related effect on HDL-C levels in the control group reflexes that the difference of the HDL-C level among genotypes was relatively small in this group of people, which is consisting with the result in a Japanese population reported by Harada et al. (42).
It is unclear why overweight/obesity modulated the effect of the polymorphism on HDL-C levels. A previous study showed that, compared with small adipocytes, large adipocytes contain reduced membrane cholesterol concentrations, and that cholesterol traffic genes, including ABCA1, are involved in the response to cholesterol depletion in the membranes (43). Another study reported that, compared with lean men, overweight men have a lower ability to promote ABCA1-mediated cholesterol efflux (44). In addition, insulin may reduce ABCA1 expression levels in adipocytes, resulting in insulin resistance (45). Therefore, ABCA1 levels are downregulated to a greater extent in overweight/obese GDM patients, who are usually hyperinsulinemic, than in normal-weight subjects. Thus, obesity or overweight may induce effects directly linked to ABCA1 function that may interact with the ABCA1 polymorphism. However, the minor A allele was associated with decreased HDL-C levels in GDM patients, especially in those with overweight/obesity, which is in contrast to several previous reports in the general population (26, 27, 46). Similar results were also observed in obese non-pregnant subjects in the French D.E.S.I.R. study (47). The authors of the D.E.S.I.R. proposed that, when the minor A allele is associated with higher HDL levels, overweight/obesity cancels this effect, and when there is an additional minor allele associated with lower levels (as for the 378 C allele of rs1800978), overweight enhances this effect, leading to a lower HDL-C concentration outcome. We speculate that there might be a similar explanation for the association of the minor A allele with lower HDL-C levels in overweight/obese GDM patients in the present study.
While we observed an association between the ABCA1 polymorphism and HDL-C levels in GDM patients, we also found a BMI-dependent relationship between this polymorphism and other lipid-related phenotypes, such as TC and/or LDL-C levels. In overweight/obese GDM patients, significantly lower TC levels were observed between patients with AA and GA genotypes (p < 0.05). In contrast, non-obese GDM subjects with the AA genotype showed increased TC levels compared to those with GG or GA genotypes (all p < 0.05). This difference may be attributed to the different BMI statuses of GDM patients. Several studies have shown that the A allele of this polymorphism is significantly associated with lower TC levels (29–31). The result observed in the overweight/obese GDM group in this study is in line with previous reports. The result observed in non-obese GDM women may be attributed to smaller BMI-mediated effect on subjects with the AA genotype, since in the non-obese GDM group, there were clear genotype-related effects on prepregnancy BMI and BMI at delivery. The relatively lower BMI in subjects with the AA genotype compared to those with the GG or GA genotypes may be linked to changes in ABCA1 levels (and thus, perhaps protein function). Indeed, a previous study demonstrated that the expression of ABCA1 is downregulated to a greater extent in obese subjects (relatively higher BMI) than in normal-weight subjects (relatively lower BMI), suggesting that BMI-related effects on cholesterol levels may also exist among non-obese GDM patients with different genotypes.
In addition, in control subjects with normal pregnancies, the AA genotype was associated with decreased TC and LDL-C levels compared with the GG genotype (both p < 0.05). This is consistent with the results of previous studies (28–33). When overweight/obese and non-obese subgroups were further divided, genotype-related effects were only evident in overweight/obese control subjects. This indicates that the genotype-related effect on these two cholesterol fractions may also be BMI-dependent in normal pregnant women.
In the current study, overweight/obese GDM patients with the AA genotype showed significantly higher DBP (75.11 ± 7.99 mmHg) than those with the GA genotype (72.15 ± 8.59 mmHg). However, this difference was not observed in non-obese subjects with GDM. Clinical diabetes, including GDM, is considered an independent risk factor for cardiovascular complications. Women with GDM have a higher risk of developing hypertensive disorders or preeclampsia. One study using transcriptomics showed in the placenta a high level of pro-inflammatory factor expression with signs of major vascular dysfunction (48). Another study found increased levels of an array of pro-inflammatory molecules (49, 50), as well as a pro-inflammatory pattern of upregulated adipokines and decreased levels of the anti-inflammatory factor, adiponectin, in adipose tissue from women with GDM (51). This pro-inflammatory milieu, together with the dysregulated secretion of placental factors and/or β-cell injury, may trigger metabolic and cardiovascular abnormalities, including hypertension, in women with GDM and their offspring (51).
Several studies have established an association between lipid-related gene polymorphisms and blood pressure variation or hypertension in different populations and ethnic groups. For example, one study reported that the genotype frequencies of ABCA1 V825I were different between normotensive and hypertensive women, and DBP differed according to ABCA1 genotype (52). Yin et al. also reported that this polymorphism, as well as several other lipid-related gene polymorphisms, interact with overweight/obesity to modulate blood pressure (SBP, DBP and pulse pressure) levels (53). In addition, another report demonstrated that the -14 C→T polymorphism of the ABCA1 gene is a determinant of blood pressure and the development of hypertension in Japanese individuals (54), and nine lipid-related SNPs have interactions with the effect of cigarette smoking on blood pressure levels (55). These findings, together with our current study results, indicate the importance of lipid-related gene polymorphisms (including ABCA1 SNPs) in mediating blood pressure variation, and modifying the additional interaction of components such as overweight/obese and smoking.
In the present study, we also found an association between the ABCA1 rs2230806 SNP and AI value in all GDM patients and in individual overweight/obese and non-obese subgroups. AI values were higher in AA homozygotes in the GDM group and its separate subgroups than in GG homozygotes (Tables 4 and 5, all p < 0.05). These data strongly support the notion that patients with GDM have a higher risk of CAD and suggest that the genetic component is associated with atherogenic risk markers, including the AI values.
Several epidemiologic studies have shown positive associations between maternal gestational weight gain (GWG) and prepregnancy BMI and cardiometabolic risk factors in the offspring, such as fat mass (56), weight trajectory (57), waist circumference, BMI, and blood pressure (58). Recent epigenetic epidemiological studies suggest that low maternal prepregnancy BMI (59, 60) and GWG (61) are related to differential peripheral blood DNA methylation in offspring, indicating the importance of maternal prepregnancy BMI on the potential life-long health of offspring. In this study, we found an association between the ABCA1 rs2230806 SNP and prepregnancy BMI variation, which suggests that the maternal genetic component may have implications in determining the cardiometabolic health of the offspring by affecting maternal prepregnancy BMI.
Additionally, we have analyzed our data by regression analysis. The results demonstrated that prepregnancy BMI was an important predictor of TG, TC, HDL-C, and LDL-C levels, while the allelic variant, and both the allelic variant and glycemic status were also important predictors of HDL-C and LDL-C levels, respectively (Table 7). These results highlight the importance of prepregnancy BMI in determining these lipid phenotypes, and the genetic component in modulating HDL-C and LDL-C concentrations. Moreover, these results suggest that ABCA1 genetic variation was significantly associated with some lipid abnormalities in women with GDM and with metabolic adaption in normal pregnancies.
Our study has some strengths and limitations that should be noted. The strengths of this study are as follows. First, the study included both maternal and fetal/neonatal data. Second, the study adds to the limited body of evidence regarding the genetics of hyperglycemia during pregnancy. Third, we provided insights into the genetic regulator of lipid metabolism and its effect on biomarkers of insulin resistance in human pregnancy. The study limitations include a lack of data on ABCA1 expression levels according to different ABCA1 genotypes, which would have provided further insights into the mechanisms responsible for the genetic association. In addition, the sample size included in this study is relatively small. Moreover, our findings were limited to descriptive characteristics.
In conclusion, our study found no evidence that the ABCA1 rs2230806 SNP is associated with GDM in Southwest Chinese women. However, it demonstrated that the rs2230806 polymorphism in the ABCA1 gene is associated with HDL-C levels and AI values in GDM patients, and with TC and LDL-C levels in control subjects. In addition, a subgroup analysis showed that this SNP was also associated with DBP levels in addition to HDL-C levels and AI values in overweight/obese GDM patients, while it was associated with TC levels, AI values, pre-pregnancy BMI, and BMI at delivery in non-obese GDM women. Moreover, this SNP was associated with TC, LDL-C, and apoB levels in overweight/obese control subjects. Our results suggest that the rs2230806 polymorphism may be associated with increased maternal cardiovascular risk in women with GDM, and thus potentially, with an increased future risk of type 2 diabetes and CHD in mothers and their offspring. Screening for this variant in GDM patients, especially those with different BMI levels, and weight loss through changes in diet and lifestyle in women with GDM, may be beneficial and have preventive effects against complications in these patients.
Data Availability Statement
The datasets for this study can be found at the found here:ABCA1 (rs2230806): https://www.ncbi.nlm.nih.gov/snp/rs2230806.
Ethics Statement
The studies involving human participants were reviewed and approved by The Institutional Review Board of West China Second University Hospital, Sichuan University. The patients/participants provided their written informed consent to participate in this study.
Author Contributions
HB conceived and designed the experiments, analyzed the data and revised the paper. FT performed experiments and wrote the paper. LG and XL was responsible for patient screening. MZ, YW, RL, YL, and SL helped with the experiments. PF and DL helped with the experiments and revised the paper. All authors have read and approved the final manuscript.
Funding
This study was supported by grants from the National Natural Science Foundation of China 39870749, the National Key Research and Development Program of China 2016YFC1000400, Program for Changjiang Scholars and Innovative Research Team in University IRT0935, the Fundamental Research Funds for the Central Universities, and Research Seed Fund from West China Second University Hospital of Sichuan University (HB).
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s Note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Acknowledgments
The authors thank all the study subjects who donated blood samples. We are also thankful to Qian Gao, Xiaoli Yan, Zejun Li, and Fangyuan Luo for collecting the samples.
References
1. Classification and Diagnosis of Diabetes: American Diabetes Association (ADA), Standards of Medical Care in Diabetes-2018. Diabetes Care (2018) 41(Suppl 1):S13–s27. doi: 10.2337/dc18-S002
2. Hod M, Kapur A, Sacks DA, Hadar E, Agarwal M, Di Renzo GC, et al. The International Federation of Gynecology and Obstetrics (FIGO) Initiative on Gestational Diabetes Mellitus: A Pragmatic Guide for Diagnosis, Management, and Care. Int J Gynaecol Obstet (2015) 131 Suppl 3:S173–211. doi: 10.1016/s0020-7292(15)30033-3
3. Burlina S, Dalfrà MG, Lapolla A. Clinical and Biochemical Approach to Predicting Post-Pregnancy Metabolic Decompensation. Diabetes Res Clin Pract (2018) 145:178–83. doi: 10.1016/j.diabres.2018.02.035
4. Gao C, Sun X, Lu L, Liu F, Yuan J. Prevalence of Gestational Diabetes Mellitus in Mainland China: A Systematic Review and Meta-Analysis. J Diabetes Investig (2019) 10(1):154–62. doi: 10.1111/jdi.12854
5. Practice Bulletin No ACOG. 190: Gestational Diabetes Mellitus. Obstet Gynecol (2018) 131(2):e49–64. doi: 10.1097/aog.0000000000002501
6. Bellamy L, Casas JP, Hingorani AD, Williams D. Type 2 Diabetes Mellitus After Gestational Diabetes: A Systematic Review and Meta-Analysis. Lancet (2009) 373(9677):1773–9. doi: 10.1016/s0140-6736(09)60731-5
7. Zhang C, Bao W, Rong Y, Yang H, Bowers K, Yeung E, et al. Genetic Variants and the Risk of Gestational Diabetes Mellitus: A Systematic Review. Hum Reprod Update (2013) 19(4):376–90. doi: 10.1093/humupd/dmt013
8. Ryckman KK, Spracklen CN, Smith CJ, Robinson JG, Saftlas AF. Maternal Lipid Levels During Pregnancy and Gestational Diabetes: A Systematic Review and Meta-Analysis. Bjog (2015) 122(5):643–51. doi: 10.1111/1471-0528.13261
9. López-Tinoco C, Roca M, García-Valero A, Murri M, Tinahones FJ, Segundo C, et al. Oxidative Stress and Antioxidant Status in Patients With Late-Onset Gestational Diabetes Mellitus. Acta Diabetol (2013) 50(2):201–8. doi: 10.1007/s00592-011-0264-2
10. Mrizak I, Arfa A, Fekih M, Debbabi H, Bouslema A, Boumaiza I, et al. Inflammation and Impaired Endothelium-Dependant Vasodilatation in non Obese Women With Gestational Diabetes Mellitus: Preliminary Results. Lipids Health Dis (2013) 12:93. doi: 10.1186/1476-511x-12-93
11. Langmann T, Klucken J, Reil M, Liebisch G, Luciani MF, Chimini G, et al. Molecular Cloning of the Human ATP-Binding Cassette Transporter 1 (Habc1): Evidence for Sterol-Dependent Regulation in Macrophages. Biochem Biophys Res Commun (1999) 257(1):29–33. doi: 10.1006/bbrc.1999.0406
12. Dubé E, Ethier-Chiasson M, Lafond J. Modulation of Cholesterol Transport by Insulin-Treated Gestational Diabetes Mellitus in Human Full-Term Placenta. Biol Reprod (2013) 88(1):16. doi: 10.1095/biolreprod.112.105619
13. Sun Y, Kopp S, Strutz J, Gali CC, Zandl-Lang M, Fanaee-Danesh E, et al. Gestational Diabetes Mellitus Modulates Cholesterol Homeostasis in Human Fetoplacental Endothelium. Biochim Biophys Acta Mol Cell Biol Lipids (2018) 1863(9):968–79. doi: 10.1016/j.bbalip
14. Joshi A, Azuma R, Akumuo R, Goetzl L, Pinney SE. Gestational Diabetes and Maternal Obesity are Associated With Sex-Specific Changes in Mirna and Target Gene Expression in the Fetus. Int J Obes (Lond) (2020) 44(7):1497–507. doi: 10.1038/s41366-019-0485-y
15. Liu L, Zhang M, Li L, Li C, Min X, Zheng M. Expression and Clinical Significance of ATP-Binding Cassette Transporter 1 in Serum and Placental Tissue in Chinese Patients With Preeclampsia. Gynecol Obstet Invest (2014) 78(3):194–200. doi: 10.1159/000362551
16. Liu M, Chung S, Shelness GS, Parks JS. Hepatic ABCA1 Deficiency is Associated With Delayed Apolipoprotein B Secretory Trafficking and Augmented VLDL Triglyceride Secretion. Biochim Biophys Acta Mol Cell Biol Lipids (2017) 1862(10 Pt A):1035–43. doi: 10.1016/j.bbalip.2017.07.001
17. Chung S, Gebre AK, Seo J, Shelness GS, Parks JS. A Novel Role for ABCA1-Generated Large Pre-Beta Migrating Nascent HDL in the Regulation of Hepatic VLDL Triglyceride Secretion. J Lipid Res (2010) 51(4):729–42. doi: 10.1194/jlr.M900083
18. Brunham LR, Kruit JK, Pape TD, Timmins JM, Reuwer AQ, Vasanji Z, et al. Beta-Cell ABCA1 Influences Insulin Secretion, Glucose Homeostasis and Response to Thiazolidinedione Treatment. Nat Med (2007) 13(3):340–7. doi: 10.1038/nm1546
19. de Haan W, Bhattacharjee A, Ruddle P, Kang MH, Hayden MR. ABCA1 in Adipocytes Regulates Adipose Tissue Lipid Content, Glucose Tolerance, and Insulin Sensitivity. J Lipid Res (2014) 55(3):516–23. doi: 10.1194/jlr.M045294
20. Sánchez-Aguilera P, Diaz-Vegas A, Campos C, Quinteros-Waltemath O, Cerda-Kohler H, Barrientos G, et al. Role of ABCA1 on Membrane Cholesterol Content, Insulin-Dependent Akt Phosphorylation and Glucose Uptake in Adult Skeletal Muscle Fibers From Mice. Biochim Biophys Acta Mol Cell Biol Lipids (2018) 1863(12):1469–77. doi: 10.1016/j.bbalip.2018.09.005
21. Bhattacharjee J, Ietta F, Giacomello E, Bechi N, Romagnoli R, Fava A, et al. Expression and Localization of ATP Binding Cassette Transporter A1 (ABCA1) in First Trimester and Term Human Placenta. Placenta (2010) 31(5):423–30. doi: 10.1016/j.placenta.2010.02.015
22. Shrestha E, Hussein MA, Savas JN, Ouimet M, Barrett TJ, Leone S, et al. Poly(ADP-Ribose) Polymerase 1 Represses Liver X Receptor-Mediated ABCA1 Expression and Cholesterol Efflux in Macrophages. J Biol Chem (2016) 291(21):11172–84. doi: 10.1074/jbc.M116.726729
23. Sedgeman LR, Michell DL, Vickers KC. Integrative Roles of Micrornas in Lipid Metabolism and Dyslipidemia. Curr Opin Lipidol (2019) 30(3):165–71. doi: 10.1097/mol.0000000000000603
24. Hodoğlugil U, Williamson DW, Huang Y, Mahley RW. Common Polymorphisms of ATP Binding Cassette Transporter A1, Including a Functional Promoter Polymorphism, Associated With Plasma High Density Lipoprotein Cholesterol Levels in Turks. Atherosclerosis (2005) 183(2):199–212. doi: 10.1016/j.atherosclerosis.2005.03.004
25. Clee SM, Zwinderman AH, Engert JC, Zwarts KY, Molhuizen HO, Roomp K, et al. Common Genetic Variation in ABCA1 is Associated With Altered Lipoprotein Levels and a Modified Risk for Coronary Artery Disease. Circulation (2001) 103(9):1198–205. doi: 10.1161/01.cir.103.9.1198
26. Yamakawa-Kobayashi K, Yanagi H, Yu Y, Endo K, Arinami T, Hamaguchi H. Associations Between Serum High-Density Lipoprotein Cholesterol or Apolipoprotein AI Levels and Common Genetic Variants of the ABCA1 Gene in Japanese School-Aged Children. Metabolism (2004) 53(2):182–6. doi: 10.1016/j.metabol.2003.08.009
27. Huang Y, Wu Y, Liu R, Fan P, Zhang J, Wang F, et al. Differential Effect of ATP Binding Cassette Transporter A1 R219K and Cholesteryl Ester Transfer Protein Taqib Genotypes on HDL-C Levels in Overweight/Obese and Non-Obese Chinese Subjects. Acta Cardiol (2011) 66(2):231–7. doi: 10.1080/ac.66.2.2071256
28. Benton JL, Ding J, Tsai MY, Shea S, Rotter JI, Burke GL, et al. Associations Between Two Common Polymorphisms in the ABCA1 Gene and Subclinical Atherosclerosis: Multi-Ethnic Study of Atherosclerosis (MESA). Atherosclerosis (2007) 193(2):352–60. doi: 10.1016/j.atherosclerosis.2006.06.024
29. Fawzy MS, Alhadramy O, Hussein MH, Ismail HM, Ismail NM, Biomy NM, et al. Functional and Structural Impact of ATP-Binding Cassette Transporter A1 R219K and I883M Gene Polymorphisms in Obese Children and Adolescents. Mol Diagn Ther (2015) 19(4):221–34. doi: 10.1007/s40291-015-0150-7
30. Katzov H, Bennet AM, Höglund K, Wiman B, Lütjohann D, Brookes AJ, et al. Quantitative Trait Loci in ABCA1 Modify Cerebrospinal Fluid Amyloid-Beta 1-42 and Plasma Apolipoprotein Levels. J Hum Genet (2006) 51(3):171–9. doi: 10.1007/s10038-005-0341-x
31. Young Guk K, Eun Young C, Hyun Young P, Yang Soo J, Sook K, Jong Eun L. Association of R219K Polymorphism in the ABCA1 Gene With Plasma Lipid Levels and Coronary Artery Disease in Koreans. Korean Circ J (2003) 33:44–51. doi: 10.4070/kcj.2003.33.1.44
32. Çoban N, Onat A, Kömürcü Bayrak E, Güleç Ç, Can G, Erginel Ünaltuna N. Gender Specific Association of ABCA1 Gene R219K Variant in Coronary Disease Risk Through Interactions With Serum Triglyceride Elevation in Turkish Adults. Anadolu Kardiyol Derg (2014) 14(1):18–25. doi: 10.5152/akd.2013.234
33. Abd El-Aziz TA, Mohamed RH, Hagrass HA. Increased Risk of Premature Coronary Artery Disease in Egyptians With ABCA1 (R219K), CETP (Taqib), and LCAT (4886C/T) Genes Polymorphism. J Clin Lipidol (2014) 8(4):381–9. doi: 10.1016/j.jacl.2014.06.001
34. Ma XY, Liu JP, Song ZY. Associations of the ATP-Binding Cassette Transporter A1 R219K Polymorphism With HDL-C Level and Coronary Artery Disease Risk: A Meta-Analysis. Atherosclerosis (2011) 215(2):428–34. doi: 10.1016/j.atherosclerosis.2011.01.008
35. Cenarro A, Artieda M, Castillo S, Mozas P, Reyes G, Tejedor D, et al. A Common Variant in the ABCA1 Gene is Associated With a Lower Risk for Premature Coronary Heart Disease in Familial Hypercholesterolaemia. J Med Genet (2003) 40(3):163–8. doi: 10.1136/jmg.40.3.163
36. Jung D, Cao S, Liu M. Park s. A Meta-Analysis of the Associations Between the ATP-Binding Cassette Transporter ABCA1 R219K (Rs2230806) Polymorphism and the Risk of Type 2 Diabetes in Asians. Horm Metab Res (2018) 50(4):308–16. doi: 10.1055/a-0583-0201
37. Metzger BE, Gabbe SG, Persson B, Buchanan TA, Catalano PA, Damm P, et al. International Association of Diabetes and Pregnancy Study Groups Recommendations on the Diagnosis and Classification of Hyperglycemia in Pregnancy. Diabetes Care (2010) 33(3):676–82. doi: 10.2337/dc09-1848
38. Erlich HA. PCR Technology: Principles and Applications for DNA Amplification. 1st Ed. New York: Stockton (1989). p. 36.
39. Demirel F, Bideci A, Cinaz P, Camurdan MO, Biberoğlu G, Yesilkaya E, et al. Serum Leptin, Oxidized Low Density Lipoprotein and Plasma Asymmetric Dimethylarginine Levels and Their Relationship With Dyslipidaemia in Adolescent Girls With Polycystic Ovary Syndrome. Clin Endocrinol (Oxf) (2007) 67(1):129–34. doi: 10.1111/j.1365-2265.2007.02849.x
40. Guan L, Fan P, Liu X, Zhou M, Wu Y, Liu R, et al. Maternal GALNT2 Variations Affect Blood Pressure, Atherogenic Index, and Fetal Growth, Depending on BMI in Gestational Diabetes Mellitus. Front Endocrinol (Lausanne) (2021) 12:690229. doi: 10.3389/fendo.2021.690229
41. Phillips MC. Is ABCA1 a Lipid Transfer Protein? J Lipid Res (2018) 59(5):749–63. doi: 10.1194/jlr.R082313
42. Harada T, Imai Y, Nojiri T, Morita H, Hayashi D, Maemura K, et al. A Common Ile 823 Met Variant of ATP-Binding Cassette Transporter A1 Gene (ABCA1) Alters High Density Lipoprotein Cholesterol Level in Japanese Population. Atherosclerosis (2003) 169:105–12. doi: 10.1016/S0021-9150(03)00135-7
43. Le Lay S, Krief S, Farnier C, Lefrère I, Le Liepvre X, Bazin R, et al. Cholesterol, a Cell Size-Dependent Signal That Regulates Glucose Metabolism and Gene Expression in Adipocytes. J Biol Chem (2001) 276(20):16904–10. doi: 10.1074/jbc.M010955200
44. Beulens JW, Sierksma A, van Tol A, Fournier N, van Gent T, Paul JL, et al. Moderate Alcohol Consumption Increases Cholesterol Efflux Mediated by ABCA1. J Lipid Res (2004) 45(9):1716–23. doi: 10.1194/jlr.M400109-JLR200
45. Sartipy P, Loskutoff DJ. Expression Profiling Identifies Genes That Continue to Respond to Insulin in Adipocytes Made Insulin-Resistant by Treatment With Tumor Necrosis Factor-Alpha. J Biol Chem (2003) 278(52):52298–306. doi: 10.1074/jbc.M306922200
46. Woll PS, Hanson NQ, Arends VL, Tsai MY. Effect of Two Common Polymorphisms in the ATP Binding Cassette Transporter A1 Gene on HDL-Cholesterol Concentration. Clin Chem (2005) 51(5):907–9. doi: 10.1373/clinchem.2004.047126
47. Porchay I, Péan F, Bellili N, Royer B, Cogneau J, Chesnier MC, et al. ABCA1 Single Nucleotide Polymorphisms on High-Density Lipoprotein-Cholesterol and Overweight: The D.E.S.I.R. Study. Obes (Silver Spring) (2006) 14(11):1874–9. doi: 10.1038/oby.2006.217
48. Radaelli T, Varastehpour A, Catalano P, Hauguel-de Mouzon S. Gestational Diabetes Induces Placental Genes for Chronic Stress and Inflammatory Pathways. Diabetes (2003) 52(12):2951–8. doi: 10.2337/diabetes.52.12.2951
49. Rota S, Yildirim B, Kaleli B, Aybek H, Duman K. Kaptanoğlu B. C-Reactive Protein Levels in non-Obese Pregnant Women With Gestational Diabetes. Tohoku J Exp Med (2005) 206(4):341–5. doi: 10.1620/tjem.206.341
50. Jensen LA, Chik CL, Ryan EA. Review of Gestational Diabetes Mellitus Effects on Vascular Structure and Function. Diabetes Vasc Dis Res (2016) 13(3):170–82. doi: 10.1177/1479164115624681
51. Pantham P, Aye IL, Powell TL. Inflammation in Maternal Obesity and Gestational Diabetes Mellitus. Placenta (2015) 36(7):709–15. doi: 10.1016/j.placenta.2015.04.006
52. Yin RX, Wu JZ, Liu WY, Wu DF, Cao XL, Miao L, et al. Association of Several Lipid-Related Gene Polymorphisms and Blood Pressure Variation in the Bai Ku Yao Population. Am J Hypertens (2012) 25(8):927–36. doi: 10.1038/ajh.2012.55
53. Yin RX, Wu DF, Miao L, Aung LH, Cao XL, Yan TT, et al. Several Genetic Polymorphisms Interact With Overweight/Obesity to Influence Serum Lipid Levels. Cardiovasc Diabetol (2012) 11:123. doi: 10.1186/1475-2840-11-123
54. Yamada Y, Kato K, Yoshida T, Yokoi K, Matsuo H, Watanabe S, et al. Association of Polymorphisms of ABCA1 and ROS1 With Hypertension in Japanese Individuals. Int J Mol Med (2008) 21(1):83–9. doi: 10.3892/ijmm.21.1.83
55. Yin RX, Wu DF, Wu JZ, Cao XL, Aung LH, Miao L, et al. Interactions of Several Lipid-Related Gene Polymorphisms and Cigarette Smoking on Blood Pressure Levels. Int J Biol Sci (2012) 8(5):685–96. doi: 10.7150/ijbs.4401
56. Santos S, Severo M, Gaillard R, Santos AC, Barros H, Oliveira A. The Role of Prenatal Exposures on Body Fat Patterns at 7 Years: Intrauterine Programming or Birthweight Effects? Nutr Metab Cardiovasc Dis (2016) 26(11):1004–10. doi: 10.1016/j.numecd.2016.06.010
57. Lawrence GM, Shulman S, Friedlander Y, Sitlani CM, Burger A, Savitsky B, et al. Associations of Maternal Pre-Pregnancy and Gestational Body Size With Offspring Longitudinal Change in BMI. Obes (Silver Spring) (2014) 22(4):1165–71. doi: 10.1002/oby.20643
58. Hochner H, Friedlander Y, Calderon-Margalit R, Meiner V, Sagy Y, Avgil-Tsadok M, et al. Associations of Maternal Prepregnancy Body Mass Index and Gestational Weight Gain With Adult Offspring Cardiometabolic Risk Factors: The Jerusalem Perinatal Family Follow-Up Study. Circulation (2012) 125(11):1381–9. doi: 10.1161/circulationaha.111.070060
59. Sharp GC, Lawlor DA, Richmond RC, Fraser A, Simpkin A, Suderman M, et al. Maternal Pre-Pregnancy BMI and Gestational Weight Gain, Offspring DNA Methylation and Later Offspring Adiposity: Findings From the Avon Longitudinal Study of Parents and Children. Int J Epidemiol (2015) 44(4):1288–304. doi: 10.1093/ije/dyv042
60. Sharp GC, Salas LA, Monnereau C, Allard C, Yousefi P, Everson TM, et al. Maternal BMI at the Start of Pregnancy and Offspring Epigenome-Wide DNA Methylation: Findings From the Pregnancy and Childhood Epigenetics (PACE) Consortium. Hum Mol Genet (2017) 26(20):4067–85. doi: 10.1093/hmg/ddx290
Keywords: ABCA1 gene, single nucleotide polymorphism (SNP), gestational diabetes mellitus, obesity, atherometabolic traits
Citation: Tang F, Guan L, Liu X, Fan P, Zhou M, Wu Y, Liu R, Liu Y, Liu S, Li D and Bai H (2021) A Common R219K Variant of ATP-Binding Cassette Transporter A1 Gene Alters Atherometabolic Traits in Pregnant Women With Gestational Diabetes Mellitus. Front. Endocrinol. 12:782453. doi: 10.3389/fendo.2021.782453
Received: 24 September 2021; Accepted: 30 November 2021;
Published: 17 December 2021.
Edited by:
Ondřej Šeda, Charles University, CzechiaReviewed by:
Alpna Tyagi, University of Colorado Anschutz Medical Campus, United StatesRamkumar Mohan, University of Michigan, United States
Manel Ben-Aissa, University of Illinois at Chicago, United States
Copyright © 2021 Tang, Guan, Liu, Fan, Zhou, Wu, Liu, Liu, Liu, Li and Bai. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Huai Bai, baihuai@scu.edu.cn