- Research article
- Open Access
The transcriptome-wide association search for genes and genetic variants which associate with BMI and gestational weight gain in women with type 1 diabetes
Molecular Medicine volume 27, Article number: 6 (2021)
Clinical data suggest that BMI and gestational weight gain (GWG) are strongly interconnected phenotypes; however, the genetic basis of the latter is rather unclear. Here we aim to find genes and genetic variants which influence BMI and/or GWG.
We have genotyped 316 type 1 diabetics using Illumina Infinium Omni Express Exome-8 v1.4 arrays. The GIANT, ARIC and T2D-GENES summary statistics were used for TWAS (performed with PrediXcan) in adipose tissue. Next, the analysis of association of imputed expression with BMI in the general and diabetic cohorts (Analysis 1 and 2) or GWG (Analysis 3 and 4) was performed, followed by variant association analysis (1 Mb around identified loci) with the mentioned phenotypes.
In Analysis 1 we have found 175 BMI associated genes and 19 variants (p < 10–4) which influenced GWG, with the strongest association for rs11465293 in CCL24 (p = 3.18E−06). Analysis 2, with diabetes included in the model, led to discovery of 1812 BMI associated loci and 207 variants (p < 10–4) influencing GWG, with the strongest association for rs9690213 in PODXL (p = 9.86E−07). In Analysis 3, among 648 GWG associated loci, 2091 variants were associated with BMI (FDR < 0.05). In Analysis 4, 7 variants in GWG associated loci influenced BMI in the ARIC cohort.
Here, we have shown that loci influencing BMI might have an impact on GWG and GWG associated loci might influence BMI, both in the general and T1DM cohorts. The results suggest that both phenotypes are related to insulin signaling, glucose homeostasis, mitochondrial metabolism, ubiquitinoylation and inflammatory responses.
In the recent years, significant attention was paid to transcriptomic-wide association studies (TWAS), which enable the prediction of gene-level associations between the 'imputed expression based on genetic data' and complex phenotypes (Gamazon et al. 2015). This is possible, since it was shown that a proportion of GWAS risk variants co-localize with genetic variants that regulate gene expression (i.e. expression quantitative trait loci, eQTL) (Hormozdiari et al. 2016). The prediction of gene expression based on genotype removes noise created by environmental factors as well as potential reverse causation (when the trait affects gene expression), therefore the TWAS analysis not only increases the statistical power, but also enables focusing on regions for which finding a functional interpretation of the associations is often easier (Mancuso et al. 2017; Li et al. 2018).
Gestational weight gain (GWG) has been extensively studied over the last several years as it was suggested that its inadequacy may lead to adverse both maternal (preeclampsia, hypertension, obesity later in life, cesarean section) and neonatal (preterm birth, stillbirth, inadequate neonatal weight—small or large for gestational age, obesity later in life) outcomes (Kominiarek and Peaceman 2017; Voerman et al. 2019). Simultaneously, clinicians tried to optimize pregnancy care for women with type 1 diabetes whose pregnancy outcomes (including developmental abnormalities, spontaneous abortions, neonatal hyperglycemia, neonatal hyperinsulinemia, maternal retinopathy, maternal nephropathy, preeclampsia) are far worse than in women from the general population (Celia et al. 2016). For a long time it was believed that it stemmed from maternal hyperglycemia; however, even though their glycemic goals have been achieved (HbA1c < 6.0%), large for gestational age (LGA) or macrosomic neonates births are more frequent among these women than in the general population (Scifres et al. 2014; Bashir et al. 2019; Dori-Dayan et al. 2020). Thus, there must be other potential contributors to adverse maternal and fetal outcomes (Secher et al. 2014; Mastella et al. 2018; Rys et al. 2018). Among them are also those which affect the general population—maternal body lipids, pre-pregnancy BMI and GWG (McWhorter et al. 2018).
Even though GWG can have an impact on the health of future generations, to date, most of the performed studies are retrospective or observational and their main goal was to assess the relationship between GWG and environmental factors (Nunnery et al. 2018; Siega-Riz et al. 2020), while the genetic risks for inadequate GWG have been scarcely analyzed. In a recent paper, the authors show that 43% and 26% of the variation in GWG can be explained by genetic factors in the first pregnancy and second pregnancy, respectively (Andersson et al. 2015). The first GWAS study on GWG in the general, multiethnic cohort was performed in 2018. The study has shown that 20% of the variability in GWG can be explained by maternal genetic variants. Unfortunately, it did not find significant associations of GWG with any genomic loci (Warrington et al. 2018).
Since both environmental and genetic data show that there is a correlation between pre-pregnancy BMI and GWG (Luecke et al. 2019) we took the TWAS approach to study the genetic correlation between BMI and GWG—in particular, to find loci and variants which are associated with GWG and BMI, or (on the contrary) only with GWG. We find this particularly important, as the reports on the overlap of genes and/or variants which influence both phenotypes appear to be conflicting (Lawlor et al. 2011; Kawai et al. 2019).
In this study, we aimed to investigate the genetic factors associated with GWG from the perspective of the genetics of obesity. To this aim, we performed TWAS analyses on the GIANT cohort representing the general population and two diabetic cohorts—T2D-GENES and ARIC—to find genes associated with BMI. The ARIC and T2D-GENES cohorts were used since diabetes mellitus itself might influence BMI and creates a special metabolic context, while no cohort with patients with type 1 diabetes was available. Only regions which were associated with BMI in these cohorts were subjected to the analysis of association with GWG in T1DM patients. Secondly, we also searched for reverse associations and checked whether genes in TWAS associated with GWG in the T1DM cohort might affect BMI in the ARIC or GIANT cohorts. If so, we decided to search for genes which imputed expression correlated with BMI in order to restrict the analysis of genetic association with GWG to these genes only, and conversely to check whether the genes which correlate with GWG (at the imputed expression level) do influence BMI. Thus, this work comprises 4 different analyses, as presented in Fig. 1.
Patients were recruited either in the Department of Metabolic Diseases University Hospital in Krakow or in the Division of Reproduction Department of Obstetrics, Gynecology and Gynecological Oncology, Poznan University of Medical Sciences, two important academic referral centers. Both cohorts included all consecutive pregnant women with pre-existing T1DM. The Kraków cohort included women who were referred to the clinic between the years 1998 and 2018, while the Poznań cohort of women who attended it between years 2006 and 2019. The data were collected at the time of clinic attendance. The women were either pregnant (no later than 1st trimester in Krakow, no later than 3rd trimester in Poznań) or planned pregnancy. Their clinical diagnosis of T1DM was established at least one year prior to conception. All women with diabetes who entered the pregnancy planning program received intensive diabetes management in the clinic. Forty percent of women in both cohorts planned their pregnancy. Women who did not plan their pregnancies entered the intensive diabetes care program after conception, at the first clinic visit. Pre-pregnancy care program includes education in the field of diet, self-monitoring of blood-glucose, glycemic targets, and self-adjustment of insulin doses, supplementation of folic acid and modification of therapy to improve glycemic control. The clinical characteristics in these women reflect the effects of routine diabetes management in non-pregnant patients. Two insulin regimens were used before pregnancy: multiple daily injections (MDI) or continuous subcutaneous insulin infusion (CSII) with a personal insulin pump. These patients were followed through the pregnancy. All women were registered in the clinic and medical data was created. Clinical characteristics of the women, information on the course of pregnancy and glycemic control markers were collected during regular clinic visits. Their pre-pregnancy weight was measured on the visit before pregnancy during preparation to conception. Women with unplanned pregnancy self-reported pre-pregnancy weight on first pregnancy visit in the clinic. The gestational weight gain was calculated as a difference between the weight collected on the last visit before delivery and pre-pregnancy weight. The partial clinical characteristics of these cohorts was published before (Cyganek et al. 2017; Gutaj et al. 2017). The prevalence of retinopathy was 25.6%. None of the patients showed evidence of chronic kidney disease of stage 3 or higher. Among the patients from the Poznań cohort the prevalence of retinopathy was 18%, nephropathy 13%. 16.8% of patients were hypertensive. The exclusion criteria were as follows: miscarriage, stillbirth, lack of data caused by not taking part in intensive diabetes management care or lack of regular out-patient visits in the clinic. We included only singleton pregnancies in the study. Whole blood samples were drawn and stored at − 80 °C. This study was approved by the Bioethical Committees of Jagiellonian University and Poznan University of Medical Sciences and performed according to the Helsinki Declaration. Written informed consent was collected from all patients.
DNA was extracted from whole blood with the use of automated nucleic acid extraction system Maxwell (Promega). Five hundred twenty-seven samples were genotyped on Illumina Infinium Omni Express Exome-8 v1.4 arrays. Only in-term live births with information regarding age, pre-pregnancy BMI, GWG, diabetes duration, treatment method and daily insulin dose available were included in the final analysis.
Data processing and imputation
The detailed protocol of the data processing, QC analysis and imputation is presented in (Ludwig-Słomczyńska et al. 2020).
Genotype data and the analysis of genetic variants
The GWAS analysis on the T1DM cohort was performed on a group of 316 females with complete data on: age, parity, insulin dose prior pregnancy, pre-pregnancy BMI and GWG. The mixed-effects model approach was used as implemented in the package GENESIS in R. The random effect was associated with the individual ID and the genetic relatedness matrix was estimated via the PCARelate method in package GENESIS in R. The 'null' mixed effects model was considered as:
where the OriginOfSample is an indicator of the sample being collected in the Division of Reproduction Department of Obstetrics, Gynecology and Gynecological Oncology, Poznan University of Medical Sciences. To the aim of testing significance of the (additive) effect of the genotype, the Wald's test statistics was used.
The summary statistics for BMI tested in the GIANT Consortium cohort were downloaded from https://portals.broadinstitute.org/collaboration/giant/images/1/15/SNP_gwas_mc_merge_nogc.tbl.uniq.gz.
The genotype and phenotype data for the ARIC cohort (GENEVA study) were accessed via dbGaP. As far as the genotype data is concerned, the imputed data for participants of European ancestry were used as deposited under the phg000248.v1 code. Per-chromosome genotype probability data were transformed to dosages. From these dosages, genotype relatedness matrices were estimated on a per-chromosome basis and merged with the aid of the SNPRelate package in R. For the GWAS analysis the GENESIS package was used. As far as the phenotype files are concerned, these were accessed under the phs000090.v3 study code. The 'null' mixed effects model was defined as:
As GENEVA is a longitudinal study, for diabetic participants, we used the BMI and age measurement at the first time when the status of the participant was 'diabetic', whereas for the non-diabetic participants we used the measurements at the study entry.
The T2Dgenes genotype files were accessed via dbGaP under the phg000573.v1 code. The raw vcf files were transformed to dosages prior analysis. The phenotype data were accessed from the phs000462.v2 study code. Similarly to the data analysis in the ARIC cohort, we recorded BMI and age at the first time when the diabetic status was positive or at the study entry, otherwise.
The TWAS analysis was performed via the PrediXcan and MetaXcan software with Adipose as the tissue of interest. For the GIANT Cohort data, the MetaXcan framework was implemented based on the GTEx.v7 models. For the ARIC, T2Dgenes and T1DM cohorts the PrediXcan software was used also with GTEx.v7 models. In the TWAS analyses in the ARIC and T1DM cohorts the same independent variables were used as with the GWAS models defined above—with the exception that the classical linear models were used as implemented in the limma package (with no random effects). For the T2Dgenes TWAS analysis, from the 'predicted expression' matrix the principal components were first estimated and the first component was added to the models as it did not correlate with any independent variable—and it differed between families (as tested with the Kruskal–Wallis test). As before, the target analysis in the T2Dgenes cohort was performed with the aid of linear models as implemented in the limma package.
The FUMA analysis was performed through the webserver at fuma.ctglab.nl.
The COJO and fast-BAT software were used as implemented by in the GTCA framework. The LD structure was estimated via the 1000 genomes data.
GO enrichment analysis was performed as implemented in the topGO package in R.
The flow of the subjects through the study as well as the accession numbers for the primary datasets used by us is present in Additional file 1: Data S1.
Gestational weight gain analysis in T1DM cohort
Our analysis comprised 316 women with T1DM for whom full phenotype data were available. The basic characterization of patients included in the analysis is presented in Table 1.
PrediXcan was used to impute gene expression in subcutaneous and visceral adipose tissue based on imputed genotype data, and the gene-based search for associations with GWG and BMI was performed using a linear model approach. Due to small cohort size, genes nominally significantly associated with GWG (442 genes in subcutaneous and 328 in visceral adipose tissue) and BMI (416 genes in subcutaneous and 326 in visceral adipose tissue) were considered for further analysis. It is worth noting that only sixty genes overlapped between the two phenotypes (Additional file 2: Table S1). We further explored this phenomenon by looking at the strength of the association between gene expression for these two traits. We found a negative correlation between logFCs of genes associated with GWG and BMI. This is consistent with clinical recommendations as patients with high pregestational BMI are advised to restrict their GWG; however, we believe that the genetic background might also play a role.
PrediXcan analysis on BMI in the general population (GIANT cohort) and variant associations with GWG in T1DM cohort—Analysis 1
Gene expression prediction in visceral and subcutaneous adipose tissue in the GIANT cohort (234,069 patients (Vinet and Zhedanov 2011)) was performed using PrediXcan software. One hundred seventy-five genes significantly associated (p < 1E−04) with BMI in the adipose (subcutaneous and visceral combined) tissue were found (Additional file 3: Table S2). The Venn diagram of genes which affected BMI, GWG or both phenotypes is presented in Additional file 4: Table S3. Only 15 genes influenced both phenotypes and 160 influenced only BMI while 633 had an impact on GWG only (Fig. 2).
We further analyzed the loci of interest and searched for variants (within the 1 Mb window around the gene) which influence BMI and/or GWG. We first mapped variants in 175 BMI associated genes (± 500,000 bp) and applied these restricted variant panels to the analysis of the association with GWG in the T1DM cohort. We filtered the results by MAF > 5% and p-value < 1E−04 and found an association with GWG for 19 variants in the T1DM cohort (Table 2). We also performed the analysis of the association with GWG only for 15 genes which showed an overlap in TWAS. These data are presented in Additional file 5: Data S2 (marked as subAnalysis 1).
The FUMA GWAS method further prioritized the variants obtained in the analysis. Twelve leading SNPs were detected. Three of them (rs7564856, rs11465293 and rs7961894) were associated with several traits in the GWAS Catalog, mainly blood parameters, but none significantly with metabolic outcome. At the same time, however, these were also associated at the level of significance p-value < 1E−03 with LDL cholesterol levels, nonalcoholic fatty liver disease, body mass index, mean arterial pressure or metabolite levels. These variants were eQTLs for 40 genes in several tissues, while rs9659938, rs11807240, rs2136682 (NEXN, NEXN-AS1), rs7564856 (SPEG), rs9352745 (SH3BGRL2), rs13340504 (CCL24), rs11465293 (CCL24), rs886131 (RPL29P25) influenced gene expression in both subcutaneous and visceral adipose. Thus, we conclude that there is a subset of genes which influence both BMI and GWG as well as genetic variants associated with GWG which co-localize to these loci.
Next, we tried to determine whether the signals for the two traits in the genes of interest are correlated. Thus, we checked whether the genetic variants associated with BMI and GWG are in linkage disequilibrium or they belong to different haplotypes. We performed the analysis of LD for all genes with variants significantly associated with GWG. However, below we present only the most interesting examples. In the gene GPN3, we have found 2 bins of eQTLs which are not in LD either with BMI or GWG variants. Variants associated with BMI also create a cluster, which is separate from a small cluster of 2 variants associated with GWG. One of GWG associated variants (rs876373) is in moderate LD with the BMI cluster (Additional file 6: Figure S1a). In another gene, PMS2P3, we found a cluster of 3 variants which are associated with GWG, all in LD with one eQTL (rs707395). At the same time, variants associated with BMI form a separate bin which is in moderate LD with a subset of eQTLs for this gene. This cluster is not in LD with the GWG bin or rs707395 (Additional file 6: Figure S1b). The third example is the STAG3L1 gene in which three clusters—GWG associated variants, BMI associated variants and eQTLs can be seen (r2), however, in each, few members of any cluster are in LD with members of other clusters (D’) making them rather dependent on each other (Additional file 6: Figure S1c).
PrediXcan analysis on BMI in the diabetic cohorts (ARIC and T2D-GENES) and variant associations with GWG in T1DM cohort—Analysis 2
Since it is known that diabetes and its treatment might impact BMI, we searched for variants that might affect BMI in the context of diabetes status. Since no T1DM cohort was available, we used T2D patients’ cohorts—T2D-GENES (590 patients) and ARIC (8746 patients) (being aware that the two phenotypes obesity and diabetes are interconnected and T2D patients are predisposed to obesity, while T1DM patients are not). After adjusting for diabetes, we found that 1812 genes were associated with BMI p < 0.05 (Additional file 7: Table S4). Among them, 135 influenced both phenotypes, 1677 influenced only BMI, and 513 influenced GWG (Additional file 8: Table S5, Fig. 3).
Again, we mapped variants to 1812 BMI associated genes (± 500000 bp) and used this panel to look for associations with GWG in T1DM cohort. When filtered by MAF > 5% and p-value < 1E−04 the analysis has shown that GWG correlates with 207 variants (Additional file 9: Table S6). Table 3 shows 20 variants which were significant at the level of p-value < 1E−05. We also performed the analysis of the association with GWG only for 135 genes which showed an overlap in TWAS. These data are presented in Additional file 5: Data 2 (marked as subAnalysis 2).
Most of the variants which influenced both GWG and BMI (p-value < 1E−05) were eQTLs associated with the expression of 13 genes (LRRFIP2, MLH1, GOLGA4, CCL26, RHBDD2, AC004980, POR, CCL26, GTF2IRD2, UPK3B, ART2P, LINC00974, KTP15, KRT17) in 14 tissues (Thyroid, Skeletal Muscle, Whole Blood, Artery, Cultured Fibroblasts, Tibial Nerve, Testis, Skin (Unexposed), Visceral Adipose Tissue, Minor Salivary Gland, Small Intestine, Colon, Brain (several regions).
FUMA analysis showed that among variants significant for association with GWG (p-value < 1E−04), 18 were listed for associations with multiple phenotypes in the GWAS Catalog. The KEGG analysis showed enrichment for the autoimmunological disease systemic lupus erythematosus (adj. p-value = 3.32E−07) and taste transduction (adj. p-value = 2.18E−02). In the molecular functions GO enrichment analysis we found these variants to be responsible, among others, for carbohydrate binding (adj. p-value = 3.15E−08), bitter taste receptor activity (adj. p-value = 4.51E−05), trace amine receptor activity (adj. p-value = 9.69E−05) and taste receptor activity (adj. p-value = 9.04E−04) (Additional file 10: Figure S2).
Finally, when comparing the lists of variants obtained in Analyses 1 and 2, we found a significant overlap, which supports our hypothesis that loci associated with BMI, both in the general population and diabetic cohorts, may also associate with GWG. Several SNPs—rs7564856 (SPEG), rs1541725 (AC009502.4), rs9347707 (PCARG), rs4329088 (PCARG), rs1978202 (CCL26), rs13340504 (CCL24), rs11465293 (CCL24), rs11037234 (RP11-111A24.2), rs876373 (RPL29P25), rs886131 (RPL29P25), rs7961894 (WDR66), rs4808209 (ZNF101) were found to influence GWG in T1DM, ARIC, T2D-GENES and GIANT cohorts. We conclude, that despite the relatively large overlap between traits at the gene level, the genetic variants which localize with these loci most often are linked to a single phenotype only.
Next, we aimed to study the “reverse association”—i.e. we asked whether genes associated with GWG influence BMI as well.
PrediXcan analysis on GWG in T1DM cohort and variant associations with BMI in the general population (GIANT cohort)—Analysis 3
Variants were mapped to 648 GWG associated genes (± 500,000 bp) (Additional file 1: Table S1) and a search for associations with BMI in the GIANT cohort was performed. Among 394,149 variants analyzed, 2091 were significantly associated with BMI (at FDR < 0.05). These significant SNPs encompassed 0.53% of the whole list of variants analyzed. In comparison, the significant SNPs in the analysis of BMI in the GIANT cohort make up 0.41%, which means that genes which influence GWG are significantly enriched in the association test for BMI.
We used the FUMA tool to prioritize these 2091 variants. Twelve lead SNPs and 12 genomic risk loci were found. Not surprisingly, the GWAS Catalog analysis has shown a significant enrichment for genes associated with BMI, body fat distribution, attendance to gym or sport groups, sleep duration, height, lipoprotein levels, low HDL-level cholesterol etc. as presented in Fig. 4 and Additional file 11: Table S7.
Thus, we find a significant enrichment of variants associated with BMI in the general population within loci associated with GWG.
PrediXcan analysis on GWG in T1DM cohort and variant associations with BMI in the diabetic cohort (ARIC)—Analysis 4
A GWAS analysis on BMI in the ARIC cohort was performed. Due to a large number of variants in the analysis we searched for lead SNPs using LD score regression. COJO analysis has returned 15 variants statistically significantly associated with the trait of interest (Table 4A).
Only 4 variants from this list were in GWG associated genes in PrediXcan, three of which overlapped with those associated with BMI (Table 4B). None of the variants were reported in GWAS Catalog or PheWas for the association with any trait. 6 variants were eQTLs—two for genes in which they were localized rs13359900 and rs73419031 for MCTP1 and SOX6, respectively; rs144768710 for ATP5G3, rs74036632 for ZFHX2-AS1, rs4131847 for HSPB8, while rs34427781 for several genes all localized 200000 bp upstream or downstream from the variant.
Since the COJO analysis was performed on loci associated with GWG (± 500000 bp around the gene), it was possible that the most important variants influencing BMI were located outside those windows. Thus, we used a different tool, fastBAT, to find genes associated with BMI. The analysis revealed 8 genes which nominally associate with the trait of interest (Table 5).
Again, none the of variants was previously mentioned to be significantly associated with any trait in GWAS Catalog or PheWas. Two were eQTLs—rs6563700 for LHFPL6, while rs12027661 for CAMK2N1.
In this work we performed TWAS as well as single variant associations to determine the loci related to GWG and/or BMI. About 15% of loci overlap between BMI and GWG at the TWAS level. In the general population (GIANT cohort) the pathway analysis has shown that loci which contribute to both phenotypes affect mainly energy metabolism. Among the enriched pathways are regulation of mitochondrion (GO:1903749, GO:1903747, GO:0010822, GO:0070585, GO:0010821, GO:0006839) and Golgi apparatus (GO:0051683, GO:0048313, GO:0051645). There are also few pathways which affect transcription, polyadenylation or methylation of RNA (GO:1903311, GO:0016071, GO:1900363, GO:0080009). The analysis of the loci which TWAS associates with GWG only leads to necrosis (GO:0070266, GO:0097300, GO:0070265), protein kinases signaling (GO:0046330, GO:0043507, GO:0032874, GO:0070304, GO:0043506, GO:0007256) and metabolism of sugars (GO:0034033, GO:0034030, GO:0033866, GO:0009226). The BMI associated loci enrich cell cycle progression pathways (GO:0010972, GO:1902750) host to pathogen signaling (GO:0043921, GO:0052312, GO:0052472) and macromolecule metabolism (GO:0010604, GO:0034641, GO:0006139, GO:0009308, GO:0046483) and insulin signaling (GO:0046626, GO:1900076). The results of the GO enrichment analysis are presented in Additional file 12: Table S8 and Additional file 13: Figures S3a, b.
The analysis of the overlap between BMI and GWG associated loci in the ARIC cohort and T1DM cohort respectively showed enrichment in the TGFβ signaling pathway (GO:0071559, GO:0071560), regulation of stem cell differentiation (GO:1901532, GO:1902036, GO:0060218, GO:2000736) and polyol pathway (GO:0019751, GO:0046173). Pathways associated with GWG were similar to those enriched in the general population, however the important contribution of genes which affect antigen presentation could be seen (GO:0002495, GO:0002504, GO:0019886). The pathways characteristic for BMI associated loci were very much involved into response to environmental stimuli i.e. temperature (GO:0009266, GO:0009408), oxygen levels (GO:0070482, GO:0001666, GO:0036293) and sensory perception (GO:0050954, GO:0007600). The results of the GO enrichment analysis are presented in Additional file 14: Table S9 and Additional file 15: Figures S4a, b.
At the genetic variant level, we detected SNPs in BMI associated loci that are related to GWG, and variants associated with both GWG and BMI (Analysis 1). At the same time, we performed a similar analysis with diabetes adjusted BMI and GWG, which led us to congruent findings (Analysis 2). These results point us towards inflammatory response, TGFβ signaling, ER stress and glucose homeostasis. Subsequently, we investigated the association of SNPs in loci which influenced GWG (based on TWAS) with BMI and found 2091 variants in the GIANT cohort (Analysis 3). A parallel analysis was done on diabetes adjusted BMI (in the ARIC cohort)—resulting in 15 variants localized to GWG associated loci (Analysis 4). The results of these analyses show an impact of lipid biosynthesis, appetite regulation, Ca2+ homeostasis (and ER stress) and inflammatory response on obesity. Our results point to the source of the genetic correlation between GWG and BMI and confirms the interconnection of the phenotypes. We compared the results of our gene focused analysis with GWAS on GWG published in 2018 (Warrington et al. 2018); even though variants from the discovery cohort did not replicate, their potential functional roles overlap with those found in our study. TMEM163 is one of the best known genes associated with obesity, LCORL was shown to associate with height, UGDH with TGFβ signaling, HLA-C with autoimmune response, HSD17B3 with fatty acids, GLRX3 with oxidative stress, RBM19 with ribosomal biogenesis, SYT4 with Ca2+ binding and pancreatic functioning, PSG5 with pregnancy development, and NTS4 with neutrophilin signaling pathway. Our results are also congruent with those obtained in the computational search for functional annotation of 445 loci associated with obesity (Cheng et al. 2018). Therefore, we hypothesize that a further search for clinical phenotypes which affect GWG might allow identification of other loci that associate with it. Below we shed more light on the biological interpretation of our results and links between them (Figs. 5, 6).
In analyses 1 and 2, the strongest association with GWG was observed for SNPs in chemokine receptor ligands—rs1978202—CCL26 and rs13340504, rs11465293—CCL24. Apart from variants located in these genes we identified several other (rs4742339, rs4796675, rs11235519) which serve as eQTLs for them. Chemokines are chemoattractant with the proinflammatory role. Their receptor expression was shown to correlate with BMI (being higher in obese individuals) (Ignacio et al. 2016) and BMI changes (decreased after bariatric surgery) (Gentili et al. 2016). CCL26 is also involved in adipose tissue beiging in response to cold (Finlin et al. 2017). At the same time, we find three variants rs9659938, rs11807240, rs2136682 associated with GWG in the NEXN gene (or serving as eQTL for it) also involved in the inflammatory response. Of note, the augmentation of NEXN antisense RNA (NEXN-AS1) inhibits TLR4 oligomerization and NFκB activity and leads to suppression of proinflammatory response (Hu et al. 2019). NEXN is also known for its abundant expression in striated muscles (Zhu et al. 2018) and cardiomyocytes; its expression is higher under high glucose conditions (Barbati et al. 2017) (Fig. 5).
Interestingly, we find GWG associated variants in genes involved in TGFβ signaling: SPEG, LINC000974, GIPC2. The GIPC2 gene (rs2136682 of which associates with GWG), belongs to the GIPC family known for regulating proliferation, cytokinesis or migration and involved in the trafficking of various transmembrane proteins. GIPC1 is necessary for cell-surface expression of transmembrane receptors, such as IGF1R and TGFβR3 (Katoh 2013; Song et al. 2016). The inhibition of LINC000974 leads to TGFβ secretion and repression of phosphorylated SMAD2 expression (Tang et al. 2014; Fang et al. 2019). Both TGFβ and SMADs are involved in insulin expression and signaling; SMAD4 deficiency was shown to improve glucose tolerance and glucose induced insulin release (Lin et al. 2009). Of note, three GWG associated variants (rs9690213, rs9393623, rs1978202) are eQTLs for LINC000974 and four (rs8080053, rs1978202, rs4796675, rs8080053) are located in the gene itself (Fig. 6).
The rs7564856, negatively associated with GWG is located in SPEG, and its expression is significantly greater in endurance than in power athletes (Kusić et al. 2020). Also, at the protein level, SPEG differs between muscles of high-response vs low response rat trainers. The coimmunoprecipitation experiment with SPEG in the former showed a significant number of proteins involved in JNK and TGFβ signaling (Kusić et al. 2020), which, as said before, is known to influence insulin signaling and ER stress response (Liu et al. 2019) (Figs. 4, 5).
Two other variants—rs9347707 and rs4329088—located in PACRG1 gene—are known to be involved in the insulin pathway. In C. elegans PACRG mutants have reduced insulin signaling. Also, PCARG plays a key role in the regulation of mitophagy (Stephenson et al. 2018) and ubiquitynylation (Loucks et al. 2016), its knockout leads to defective NFκB signaling (Meschede et al. 2020) and CHOP upregulation (Han et al. 2017), which in adipose tissue result in increased proinflammatory macrophage polarization (M1) and insulin resistance (Suzuki et al. 2017). Thus, we note that PACRG is also involved in ER stress. The RHBDD2 gene is also involved in this biological process, with six eQTLs (rs4742339, rs11235519, rs1978202, rs13340504, rs4796675 and rs13340504) associated with GWG. It is worth noting that silencing of RHBDD2 leads to increased expression of ATF6, IRE1, PERK, CRT, BiP, ATF4 and CHOP (Fig. 4).
The association of rs72849841 (in RNF213) with GWG draws our particular attention. Namely, in Akita mice (which develop diabetes spontaneously) the knockout of RNF213 lowers glucose tolerance by 20% and leads to increased insulin contents in pancreases (by 150%) compared to wild type animals. This is obtained via the inhibition of ubiquitinoylation (RNF213 is a ubiquitin ligase) as shown by the 30% lower percentage of CHOP positive B-cells (Kobayashi et al. 2013). The depletion of RNF213 protects cells from lipotoxicity (Piccolis et al. 2019), reduces palmitate-induced cell death and modifies non-mitochondrial oxygen consumption rate. Moreover, RNF213 is targeted to lipid droplets, which are specialized for neutral lipid storage and increases their abundance and stability in cells through the elimination of adipose triglyceride lipase (ATGL) (Sugihara et al. 2019). It is interesting to note that the latter is also negatively regulated by insulin (Kershaw et al. 2006). RNF213 variants in GWAS analyses suggestively associated (p-value ~ E−03) with diabetic polyneuropathy, portal hypertension, diabetic retinopathy and abnormal glucose (Fig. 4). Its expression differed between cyclic and pregnant heifers (Forde et al. 2012).
We identified GWG associated loci involved in glucose homeostasis or pancreatic functioning (Fig. 5)—e.g. rs11807240 in FUBP1, which upregulates the mRNA levels of the two hexokinase genes Hk1 and Hk2—the rate limiting enzymes of glycolysis. A positive correlation between FUBP1 mRNA and both of the hexokinases was found in several types of cancers (Kang et al. 2019). Variants in LCA5 gene were shown to modify glucose response in a clinical trial of insulin and potassium (GIK) infusion in acute coronary syndromes (Ellis et al. 2015). The rs9352745 associated with GWG is an eQTL for SH3BRGL2, whose expression is associated with diabetes (type 1, type 2 and gestational) in a meta-analysis study. An interaction between SH3BRGL2 SNP and total fat intake was found to affect LDL-PPD (Rudkowska et al. 2015). We identified GWG associated SNP rs7609526 in KANK1 related to fasting proinsulin (Huyghe et al. 2013) and fasting plasma glucose (Hwang et al. 2015) in GWAS. CADPS (with GWG associated SNP rs1398095) plays a role in insulin secretion (Speidel et al. 2008). Cg14451730 located near CADPS is differentially methylated with HbA1c levels in the Taiwan Biobank cohort (Hsiung et al. 2020). The suggestive association with fasting plasma glucose of rs487321 in the CADPS gene was found in an Arabian population of T2D patients (Hebbar et al. 2020). We also found a variant correlating with GWG in RPL29P25 gene. The presence of anti-60S ribosomal protein L29 (RPL29) antibody in human serum was shown to inhibit the proliferation of pancreatic cancer cells in various cancers and is believed to be a novel candidate for a prognostic marker for unrespectable pancreatic cancer (Muro et al. 2015). The knockdown of RPL29 leads to suppression of cell proliferation, induces cell arrest (at G0/G1phase), enhances cell apoptosis and decreases intracellular ROS generation (Li et al. 2012). The rs4808209 (in ZNF1 gene) is an eQTL for LPAR2—the receptor for a signaling lipid LPA. It is a part of the LPAR2/Gab1/PI3K/Akt pathway, which influences glucose uptake (Rodriguez-Araujo et al. 2013) and is differentially expressed in fatty vs normal liver in the extreme obese cohort (DiStefano et al. 2014). Its expression correlates with the invasiveness of pancreatic cancers (Gong et al. 2012) as well as gynecological disorders (Yu et al. 2016; Wasniewski and Woclawek-Potocka 2018; Fujii et al. 2019; Kowalczyk-Zieba et al. 2019).
In Analyses 3 and 4, we set to determine whether genes associated with GWG also impact BMI. Basing on the summary results from the GIANT cohort we show a number of such genes. Some of their SNPs were previously associated with BMI, body fat distribution, attendance to gym or sport groups, sleep duration, height, lipoprotein levels, and low HDL-level cholesterol. At the same time, in a much smaller ARIC cohort we find 15 variants associated with diabetes adjusted BMI. None of them is listed in GWAS Catalog or PheWas; however, other variants in these genes were already mentioned to be associated with BMI.
Our analysis shows association with rs12685009 in the SLC24A2 loci, which belongs to the family of K+ dependent Na+/Ca2+ exchanger (Schnetkamp 2013). A member of the family, SLC4 influences MC4R dependent satiety and the loss of MC4R function leads to prolonged obesity and reduced energy expenditure, while mice lacking SLC4 display anorexia (Li and Lytton 2014). Also, variants in SOX6, a transcription factor shown to be involved in the promotion of adipogenesis (Iguchi et al. 2005; Leow et al. 2016), were previously related to obesity (Liu et al. 2009; Correa-Rodríguez et al. 2018). The rs77346003 in the ADIPOR1 loci—a receptor for adiponectin known to regulate glucose metabolism and fatty acid oxidation—associates with BMI in our study. Its variants in GWAS are related to obesity (Beckers et al. 2013; Peters et al. 2013; Keustermans et al. 2017) as well as fetal weight (Fensterseifer et al. 2018; Muñoz-Muñoz et al. 2018). ADIPOR1 expression correlates with improved insulin sensitivity and PGC-1A expression—a master regulator of mitochondrial gene expression (Za’don et al. 2019). Moreover, our analysis identifies eQTLs for mitochondria associated genes ATP5G3 (rs144768710) and HSPB8 (rs4131847). ATP5G3 is involved in the proton pathway and acts as an energy-driving motor (Huang et al. 2013; He et al. 2017; Spataru et al. 2019). HSPB8 is known to prevent oxidative tissue damage and its expression in serum is used as a biomarker for virus induced type 1 diabetes (Karthik et al. 2012; Li et al. 2017; Yu et al. 2019). We also identified rs34427781—an eQTL for several genes (CYP2D8P, NDUFA6-AS1, CYP2D6, CCDC134, CENPM, CHADL, FAM109B, NAGA, NDUFA6, OGFRP1, OLA1P1, SEPT3, SHISA8, SLC25A5P1, SMDT1, SREBF2, TNFRSF13C, WBP2NL)—to be associated with BMI. It is worth noting that a subset of these genes is involved in lipid metabolism (NAGA, SREBP) or mitochondrial signaling (SLC25A5P1, NDUFA6, SMDT1). SREBP is a master regulator of lipid and sterol biosynthesis (Düvel et al. 2010); the gene–gene interactions between variants in the INSIG-SCAP-SREBP pathway are associated with risk of obesity in Chinese children (Liu et al. 2014). The SMDT1 gene builds the mitochondrial calcium uniporter subunit (mtCU) localized in the inner mitochondrial membrane and is responsible for the Ca2+ transport to the mitochondrial matrix (Pendin et al. 2014), while its expression was shown to associate with fetal development (Vishnyakova et al. 2019). SLC25 genes, members of SLC channels mentioned above, are known to be subunits of another complex essential for proper mitochondrial dynamics and energy production—ANC. ANC is responsible for mitochondrial/cytoplasmic ADP/ATP exchange (Clémençon et al. 2013) and its expression was shown to be consistently upregulated in obesity (Padilla et al. 2014). NDUFA6 is part of the first complex of the mitochondrial respiratory chain, which expression impacts oxidative stress and energy production efficacy (Fiedorczuk and Sazanov 2018). Lastly, the TNFRSF13C (BAFF) gene regulates insulin sensitivity (Kawasaki et al. 2013) and is associated with autoimmune diseases (Moisini and Davidson 2009).
The above results are congruent with the gene-based fBAT analysis that points us towards genes which correlate with BMI (GATAD2, LHFPL6) and immune response/autoimmunity (HLA-DQB, HLA-DRB, HLA-B). GATAD2A is involved in embryonic development (Wang et al. 2017) and associates with obesity (Saxena et al. 2012). A copy number variation of LHFPL6 associates with average daily gain in cattle (Xu et al. 2019), while its rs4073643 associates with systemic lupus erythematosus in the Chinese population (Zhang et al. 2016). While the association between HLA variants and T1DM is very well known (Zhao et al. 2016), HLA loci are linked to other phenotypes (Karnes et al. 2017): type 2 diabetes (Ng et al. 2014), fatty liver disease (Doganay et al. 2014) as well as BMI (Shen et al. 2018) or waist-to-hip circumference (Wen et al. 2016). More importantly, these variants associate with trends in fetal birth weight (Capittini et al. 2009), level of inflammation in visceral adipose tissue in pregnant women (Eyerahi et al. 2018) and different mRNA levels in visceral omental adipose tissue of pregnant women with gestational diabetes (Deng et al. 2018). The HLA-DQA1, HLADQB1 are differentially methylated in siblings born before vs after maternal bariatric surgery (Berglind et al. 2016).
In this study we identified several loci which contribute to the genetic correlation between BMI and GWG. Variants identified in those loci are associated with genes linked to insulin signaling, glucose homeostasis, mitochondrial metabolism, ubiquitinylation and inflammatory responses and placenta functioning, not only in the diabetic cohorts, but also in the general population. The genetic contribution to GWG is clearly connected with BMI associated loci.
Availability of data and materials
Sequence data has been deposited at the European Genome-phenome Archive (EGA), which is hosted by the EBI and the CRG, under accession number EGAS00001004408. Further information about EGA can be found on https://ega-archive.org. “The European Genome-phenome Archive of human data consented for biomedical research” (http://www.nature.com/ng/journal/v47/n7/full/ng.3312.html).
Transcriptomic-wide association study
Genome-wide association study
Body mass index
Gestational weight gain
Expression quantitative trait loci
Type 1 diabetes
Andersson ES, Silventoinen K, Tynelius P, Nohr EA, Sørensen TIA, Rasmussen F. Heritability of gestational weight gain - a swedish register-based twin study. Twin Res Hum Genet. 2015;18(4):410–8.
Barbati SA, Colussi C, Bacci L, Aiello A, Re A, Stigliano E, et al. Transcription factor CREM mediates high glucose response in cardiomyocytes and in a male mouse model of prolonged hyperglycemia. Endocrinology. 2017;158(7):2391–405.
Bashir M, Naem E, Taha F, Konje JC, Abou-Samra AB. Outcomes of type 1 diabetes mellitus in pregnancy; effect of excessive gestational weight gain and hyperglycaemia on fetal growth. Diabetes Metab Syndr Clin Res Rev. 2019;13(1):84–8.
Beckers S, De Freitas F, Zegers D, Mertens IL, Verrijken A, Van Camp JK, et al. No conclusive evidence for association of polymorphisms in the adiponectin receptor 1 gene, AdipoR1, with common obesity. Endocrine. 2013;43(1):120–6.
Berglind D, Müller P, Willmer M, Sinha I, Tynelius P, Näslund E, et al. Differential methylation in inflammation and type 2 diabetes genes in siblings born before and after maternal bariatric surgery. Obesity. 2016;24(1):250–61.
Capittini C, Pasi A, Bergamaschi P, Tinelli C, De Silvestri A, Mercati MP, et al. HLA haplotypes and birth weight variation: Is your future going to be light or heavy? Tissue Antigens. 2009;74(2):156–63.
Celia B, Richard H, Chineze O, Casey E, Elizabeth P, Macauley A. The St Vincent Declaration–25 years on: review of maternal deaths reported to the Confidential Enquiry in the UK and the lessons learnt. Eur J Obstet Gynecol Reprod Biol. 2016;206:e139–40.
Cheng M, Mei B, Zhou Q, Zhang M, Huang H, Han L, et al. Computational analyses of obesity associated loci generated by genome-wide association studies. PLoS ONE. 2018;13:7.
Clémençon B, Babot M, Trézéguet V. The mitochondrial ADP/ATP carrier (SLC25 family): Pathological implications of its dysfunction. Vol. 34, Molecular Aspects of Medicine. 2013. p. 485–93.
Correa-Rodríguez M, Schmidt-RioValle J, Rueda-Medina B. SOX6 rs7117858 polymorphism is associated with osteoporosis and obesity-related phenotypes. Eur J Clin Invest. 2018;48:10.
Cyganek K, Skupien J, Katra B, Hebda-Szydlo A, Janas I, Trznadel-Morawska I, et al. Risk of macrosomia remains glucose-dependent in a cohort of women with pregestational type 1 diabetes and good glycemic control. Endocrine. 2017;55(2):447–55.
Deng X, Yang Y, Sun H, Qi W, Duan Y, Qian Y. Analysis of whole genome-wide methylation and gene expression profiles in visceral omental adipose tissue of pregnancies with gestational diabetes mellitus. J Chinese Med Assoc. 2018;81(7):623–30.
DiStefano JK, Kingsley C, Craig Wood G, Chu X, Argyropoulos G, Still CD, et al. Genome-wide analysis of hepatic lipid content in extreme obesity. Acta Diabetol. 2014;52(2):373–82.
Doganay L, Katrinli S, Colak Y, Senates E, Zemheri E, Ozturk O, et al. HLA DQB1 alleles are related with nonalcoholic fatty liver disease. Mol Biol Rep. 2014;41(12):7937–43.
Dori-Dayan N, Cukierman-Yaffe T, Zemet R, Cohen O, Mazaki-Tovi S, Yoeli-Ullman R. 1013: Gestational weight gain does-not affect insulin requirement during pregnancy in women with Type 1 diabetes. Am J Obstet Gynecol. 2020;222(1):S630–1.
Düvel K, Yecies JL, Menon S, Raman P, Lipovsky AI, Souza AL, et al. Activation of a metabolic gene regulatory network downstream of mTOR complex 1. Mol Cell. 2010;39(2):171–83.
Ellis KL, Zhou Y, Beshansky JR, Ainehsazan E, Selker HP, Cupples LA, et al. Genetic modifiers of response to glucose-insulin-potassium (GIK) infusion in acute coronary syndromes and associations with clinical outcomes in the IMMEDIATE trial. Pharmacogenomics J. 2015;15(6):488–95.
Eyerahi BF, Mancilla-Herrera I, Espino Y Sosa S, Ortiz-Ramirez M, Flores-Rueda V, Ibargüengoitia-Ochoa F, et al. Macrophage populations in visceral adipose tissue from pregnant women: Potential role of obesity in maternal inflammation. Int J Mol Sci. 2018;19: 4.
Fang CY, Yu CC, Liao YW, Hsieh PL, Lu MY, Lin KC, et al. LncRNA LINC00974 activates TGF-β/Smad signaling to promote oral fibrogenesis. J Oral Pathol Med. 2019;48(2):151–8.
Fensterseifer SR, Austin KJ, Ford SP, Alexander BM. Effects of maternal obesity on maternal and fetal plasma concentrations of adiponectin and expression of adiponectin and its receptor genes in cotyledonary and adipose tissues at mid- and late-gestation in sheep. Anim Reprod Sci. 2018;197:231–9.
Fiedorczuk K, Sazanov LA. Mammalian Mitochondrial Complex I Structure and Disease-Causing Mutations. Vol. 28, Trends in Cell Biology. 2018. p. 835–67.
Finlin BS, Zhu B, Confides AL, Westgate PM, Harfmann BD, Dupont-Versteegden EE, et al. Mast cells promote seasonal white adipose beiging in humans. Diabetes. 2017;66(5):1237–46.
Forde N, Duffy GB, McGettigan PA, Browne JA, Mehta JP, Kelly AK, et al. Evidence for an early endometrial response to pregnancy in cattle: Both dependent upon and independent of interferon tau. Physiol Genomics. 2012;44(16):799–810.
Fujii T, Nagamatsu T, Schust DJ, Ichikawa M, Kumasawa K, Yabe S, et al. Placental expression of lysophosphatidic acid receptors in normal pregnancy and preeclampsia. Am J Reprod Immunol. 2019;82:5.
Gamazon ER, Wheeler HE, Shah KP, Mozaffari SV, Aquino-Michaels K, Carroll RJ, et al. A gene-based association method for mapping traits using reference transcriptome data. Nat Genet. 2015;47(9):1091–8.
Gentili A, Zaibi MS, Alomar SY, De Vuono S, Ricci MA, Alaeddin A, et al. Circulating Levels of the Adipokines Monocyte Chemotactic Protein-4 (MCP-4), Macrophage Inflammatory Protein-1β (MIP-1β), and Eotaxin-3 in Severe Obesity and Following Bariatric Surgery. Horm Metab Res. 2016;48(12):847–53.
Gong YL, Tao CJ, Hu M, Chen JF, Cao XF, Lv GM, et al. Expression of lysophosphatidic acid receptors and local invasiveness and metastasis in Chinese pancreatic cancers. Curr Oncol. 2012;19:1.
Gutaj P, Zawiejska A, Mantaj U, Wender-Ożegowska E. Determinants of preeclampsia in women with type 1 diabetes. Acta Diabetol. 2017;54(12):1115–21.
Han K, Hassanzadeh S, Singh K, Menazza S, Nguyen TT, Stevens MV, et al. Parkin regulation of CHOP modulates susceptibility to cardiac endoplasmic reticulum stress. Sci Rep. 2017;7:1.
He J, Ford HC, Carroll J, Ding S, Fearnley IM, Walker JE. Persistence of the mitochondrial permeability transition in the absence of subunit c of human ATP synthase. Proc Natl Acad Sci USA. 2017;114(13):3409–14.
Hebbar P, Abu-Farha M, Alkayal F, Nizam R, Elkum N, Melhem M, et al. Genome-wide association study identifies novel risk variants from RPS6KA1, CADPS, VARS, and DHX58 for fasting plasma glucose in Arab population. Sci Rep. 2020;10:1.
Hormozdiari F, van de Bunt M, Segrè AV, Li X, Joo JWJ, Bilow M, et al. Colocalization of GWAS and eQTL signals detects target genes. Am J Hum Genet. 2016;99(6):1245–60.
Hsiung CN, Chang YC, Lin CW, Chang CW, Chou WC, Chu HW, et al. The causal relationship of circulating triglyceride and glycated hemoglobin: a mendelian randomization study. J Clin Endocrinol Metab. 2020;105:3.
Hu YW, Guo FX, Xu YJ, Li P, Lu ZF, McVey DG, et al. Long noncoding RNA NEXN-AS1 mitigates atherosclerosis by regulating the actin-binding protein NEXN. J Clin Invest. 2019;129(3):1115–28.
Huang Y, Wang L, Bennett B, Williams RW, Wang YJ, Gu WK, et al. Potential role of Atp5g3 in epigenetic regulation of alcohol preference or obesity from a mouse genomic perspective. Genet Mol Res. 2013;12(3):3662–74.
Huyghe JR, Jackson AU, Fogarty MP, Buchkovich ML, Stančáková A, Stringham HM, et al. Exome array analysis identifies new loci and low-frequency variants influencing insulin processing and secretion. Nat Genet. 2013;45(2):197–201.
Hwang JY, Sim X, Wu Y, Liang J, Tabara Y, Hu C, et al. Genome-wide association meta-analysis identifies novel variants associated with fasting plasma glucose in East Asians. Diabetes. 2015;64(1):291–8.
Ignacio RMC, Gibbs CR, Lee ES, Son DS. Differential chemokine signature between human preadipocytes and adipocytes. Immune Netw. 2016;16(3):189–94.
Iguchi H, Ikeda Y, Okamura M, Tanaka T, Urashima Y, Ohguchi H, et al. SOX6 attenuates glucose-stimulated insulin secretion by repressing PDX1 transcriptional activity and is down-regulated in hyperinsulinemic obese mice. J Biol Chem. 2005;280(45):37669–80.
Kang M, Lee SM, Kim W, Lee KH, Kim DY. Fubp1 supports the lactate-Akt-mTOR axis through the upregulation of Hk1 and Hk2. Biochem Biophys Res Commun. 2019;512(1):93–9.
Karnes JH, Bastarache L, Shaffer CM, Gaudieri S, Xu Y, Glazer AM, et al. Phenome-wide scanning identifies multiple diseases and disease severity phenotypes associated with HLA variants. Sci Transl Med. 2017;9:389.
Karthik D, Ilavenil S, Kaleeswaran B, Sunil S, Ravikumar S. Proteomic analysis of plasma proteins in diabetic rats by 2D electrophoresis and MALDI-TOF-MS. Appl Biochem Biotechnol. 2012;166(6):1507–19.
Katoh M. Functional proteomics, human genetics and cancer biology of GIPC family members. Exp Mol Med. 2013;45:1.
Kawai VK, Nwosu SK, Kurnik D, Harrell FE, Stein CM. Variants in BMI-associated genes and adrenergic genes are not associated with gestational weight trajectory. Obesity. 2019;27(7):1184–9.
Kawasaki K, Abe M, Tada F, Tokumoto Y, Chen S, Miyake T, et al. Blockade of B-cell-activating factor signaling enhances hepatic steatosis induced by a high-fat diet and improves insulin sensitivity. Lab Investig. 2013;93(3):311–21.
Kershaw EE, Hamm JK, Verhagen LAW, Peroni O, Katic M, Flier JS. Adipose triglyceride lipase: Function, regulation by insulin, and comparison with adiponutrin. Diabetes. 2006;55(1):148–57.
Keustermans G, Van Der Heijden LB, Boer B, Scholman R, Nuboer R, Pasterkamp G, et al. Differential adipokine receptor expression on circulating leukocyte subsets in lean and obese children. PLoS ONE. 2017;12:10.
Kobayashi H, Yamazaki S, Takashima S, Liu W, Okuda H, Yan J, et al. Ablation of Rnf213 retards progression of diabetes in the Akita mouse. Biochem Biophys Res Commun. 2013;432(3):519–25.
Kominiarek MA, Peaceman AM. Gestational weight gain. Am J Obst Gynecol. 2017;217:642–51.
Kowalczyk-Zieba I, Woclawek-Potocka I, Wasniewski T, Boruszewska D, Grycmacher K, Sinderewicz E, et al. LPAR2 and LPAR4 are the main receptors responsible for LPA actions in ovarian endometriotic cysts. Reprod Sci. 2019;26(1):139–50.
Kusić D, Connolly J, Kainulainen H, Semenova EA, Borisov OV, Larin AK, et al. Striated muscle-specific serine/threonine-protein kinase beta segregates with high versus low responsiveness to endurance exercise training. Physiol Genomics. 2020;52(1):35–46.
Lawlor DA, Fraser A, Macdonald-Wallis C, Nelson SM, Palmer TM, Smith GD, et al. Maternal and offspring adiposity-related genetic variants and gestational weight gain. Am J Clin Nutr. 2011;94(1):149–55.
Leow SC, Poschmann J, Too PG, Yin J, Joseph R, McFarlane C, et al. The transcription factor SOX6 contributes to the developmental origins of obesity by promoting adipogenesis. Dev. 2016;143(6):950–61.
Li B, Verma SS, Veturi YC, Verma A, Bradford Y, Haas DW, et al. Evaluation of predixcan for prioritizing GWAS associations and predicting gene expression. In: Pacific Symposium on Biocomputing. 2018. p. 448–59.
Li C, Ge M, Yin Y, Luo M, Chen D. Silencing expression of ribosomal protein L26 and L29 by RNA interfering inhibits proliferation of human pancreatic cancer PANC-1 cells. Mol Cell Biochem. 2012;370(1–2):127–39.
Li XC, Hu QK, Chen L, Liu SY, Su S, Tao H, et al. HSPB8 promotes the fusion of autophagosome and lysosome during autophagy in diabetic neurons. Int J Med Sci. 2017;14(13):1335–41.
Li XF, Lytton J. An essential role for the K+-dependent Na+/Ca2+-exchanger, NCKX4, in melanocortin-4-receptor-dependent satiety. J Biol Chem. 2014;289(37):25445–59.
Lin HM, Lee JH, Yadav H, Kamaraju AK, Liu E, Zhigang D, et al. Transforming growth factor-β/Smad3 signaling regulates insulin gene transcription and pancreatic islet β-cell function. J Biol Chem. 2009;284(18):12246–57.
Liu FH, Song JY, Shang XR, Meng XR, Ma J, Wang HJ. The gene-gene interaction of INSIG-SCAP-SREBP pathway on the risk of obesity in Chinese children. Biomed Res Int. 2014;2014:1.
Liu YZ, Pei YF, Liu JF, Yang F, Guo Y, Zhang L, et al. Powerful bivariate Genome-wide association analyses suggest the SOX6 gene influencing both obesity and osteoporosis phenotypes in males. PLoS ONE. 2009;4:8.
Liu Z, Li C, Kang N, Malhi H, Shah VH, Maiers JL. Transforming growth factor (TGF) cross-talk with the unfolded protein response is critical for hepatic stellate cell activation. J Biol Chem. 2019;294(9):3137–51.
Loucks CM, Bialas NJ, Dekkers MPJ, Walker DS, Grundy LJ, Li C, et al. PACRG, a protein linked to ciliary motility, mediates cellular signaling. Mol Biol Cell. 2016;27(13):2133–44.
Ludwig-Słomczyńska AH, Seweryn MT, Kapusta P, Pitera E, Handelman SK, Mantaj U, et al. Mitochondrial GWAS and association of nuclear - Mitochondrial epistasis with BMI in T1DM patients. BMC Med Genomics. 2020;13:1.
Luecke E, Cohen AK, Brillante M, Rehkopf DH, Coyle J, Hendrick CE, et al. Similarities in maternal weight and birth weight across pregnancies and across sisters. Maternal Child Health J. 2019;1:138–47.
Mancuso N, Shi H, Goddard P, Kichaev G, Gusev A, Pasaniuc B. Integrating gene expression with summary association statistics to identify genes associated with 30 complex traits. Am J Hum Genet. 2017;100(3):473–87.
Mastella LS, Weinert LS, Gnielka V, Hirakata VN, Oppermann MLR, Silveiro SP, et al. Influence of maternal weight gain on birth weight: a gestational diabetes cohort. Arch Endocrinol Metab. 2018;62(1):48–56.
McWhorter KL, Bowers K, Dolan LM, Deka R, Jackson CL, Khoury JC. Impact of gestational weight gain and prepregnancy body mass index on the prevalence of large-for-gestational age infants in two cohorts of women with type 1 insulin-dependent diabetes: a cross-sectional population study. BMJ Open. 2018;8:3.
Meschede J, Šadić M, Furthmann N, Miedema T, Sehr DA, Müller-Rischart AK, et al. The Parkin-coregulated gene product PACRG promotes TNF signaling by stabilizing LUBAC. Sci Signal. 2020;13:617.
Moisini I, Davidson A. BAFF: A local and systemic target in autoimmune diseases. Clin Exp Immunol. 2009;158:155–63.
Muñoz-Muñoz E, Krause BJ, Uauy R, Casanello P. LGA-newborn from patients with pregestational obesity present reduced adiponectin-mediated vascular relaxation and endothelial dysfunction in fetoplacental arteries. J Cell Physiol. 2018;233(10):6723–33.
Muro S, Miyake Y, Kato H, Tsutsumi K, Yamamoto K. Serum anti-60s ribosomal protein L29 antibody as a novel prognostic marker for unresectable pancreatic cancer. Digestion. 2015;91(2):164–73.
Ng MCY, Shriner D, Chen BH, Li J, Chen WM, Guo X, et al. Meta-analysis of genome-wide association studies in African Americans provides insights into the genetic architecture of type 2 diabetes. PLoS Genet. 2014;10:8.
Nunnery D, Ammerman A, Dharod J. Predictors and outcomes of excess gestational weight gain among low-income pregnant women. Health Care Women Int. 2018;39(1):19–33.
Padilla J, Jenkins NT, Thorne PK, Martin JS, Scott Rector R, Wade Davis J, et al. Identification of genes whose expression is altered by obesity throughout the arterial tree. Physiol Genomics. 2014;46(22):821–32.
Pendin D, Greotti E, Pozzan T. The elusive importance of being a mitochondrial Ca2+ uniporter. Cell Calcium. 2014;55:139–45.
Peters KE, Beilby J, Cadby G, Warrington NM, Bruce DG, Davis WA, et al. A comprehensive investigation of variants in genes encoding adiponectin (ADIPOQ) and its receptors (ADIPOR1/R2), and their association with serum adiponectin, type 2 diabetes, insulin resistance and the metabolic syndrome. BMC Med Genet. 2013;14:1.
Piccolis M, Bond LM, Kampmann M, Pulimeno P, Chitraju C, Jayson CBK, et al. Probing the global cellular responses to lipotoxicity caused by saturated fatty acids. Mol Cell. 2019;74(1):32-44.e8.
Rodriguez-Araujo G, Nakagami H, Hayashi H, Mori M, Shiuchi T, Minokoshi Y, et al. Alpha-synuclein elicits glucose uptake and utilization in adipocytes through the Gab1/PI3K/Akt transduction pathway. Cell Mol Life Sci. 2013;70(6):1123–33.
Rudkowska I, Pérusse L, Bellis C, Blangero J, Després JP, Bouchard C, et al. Interaction between common genetic variants and total fat intake on low-density lipoprotein peak particle diameter: a genome-wide association study. J Nutrigenet Nutrigenomics. 2015;8(1):44–53.
Rys PM, Ludwig-Slomczynska AH, Cyganek K, Malecki MT. Continuous subcutaneous insulin infusion vs multiple daily injections in pregnant women with type 1 diabetes mellitus: a systematic review and meta-analysis of randomised controlled trials and observational studies. Eur J Endocrinol. 2018;178(5):545–63.
Saxena R, Elbers CC, Guo Y, Peter I, Gaunt TR, Mega JL, et al. Large-scale gene-centric meta-analysis across 39 studies identifies type 2 diabetes loci. Am J Hum Genet. 2012;90(3):410–25.
Schnetkamp PPM. The SLC24 gene family of Na+/Ca2+-K+ exchangers: From sight and smell to memory consolidation and skin pigmentation. Mol Aspects Med. 2013;34:455–64.
Scifres CM, Feghali MN, Althouse AD, Caritis SN, Catov JM. Effect of excess gestational weight gain on pregnancy outcomes in women with type 1 diabetes. Obstet Gynecol. 2014;123(6):1295–302.
Secher AL, Parellada CB, Ringholm L, Ásbjörnsdóttir B, Damm P, Mathiesen ER. Higher gestational weight gain is associated with increasing offspring birth weight independent of maternal glycemic control in women with type 1 diabetes. Diabetes Care. 2014;37(10):2677–84.
Shen J, Guo T, Wang T, Zhen Y, Ma X, Wang Y, et al. HLA-B∗07, HLA-DRB1∗07, HLA-DRB1∗12, and HLA-c∗03:02 strongly associate with BMI: Data from 1.3 million healthy Chinese adults. Diabetes. 2018;67(5):861–71.
Siega-Riz AM, Bodnar LM, Stotland NE, Stang J. The Current Understanding of Gestational Weight Gain Among Women with Obesity and the Need for Future Research. NAM Perspect. 2020;
Song NJ, Kim S, Jang BH, Chang SH, Yun UJ, Park KM, et al. Small molecule-induced complement factor D (Adipsin) promotes lipid accumulation and adipocyte differentiation. PLoS ONE. 2016;11:9.
Spataru A, Le Duc D, Zagrean L, Zagrean AM. Ethanol exposed maturing rat cerebellar granule cells show impaired energy metabolism and increased cell death after oxygen-glucose deprivation. Neural Regen Res. 2019;14(3):485–90.
Speidel D, Salehi A, Obermueller S, Lundquist I, Brose N, Renström E, et al. CAPS1 and CAPS2 regulate stability and recruitment of insulin granules in mouse pancreatic β cells. Cell Metab. 2008;7(1):57–67.
Stephenson SEM, Aumann TD, Taylor JM, Riseley JR, Li R, Mann JR, et al. Generation and characterisation of a parkin-Pacrg knockout mouse line and a Pacrg knockout mouse line. Sci Rep. 2018;8:1.
Sugihara M, Morito D, Ainuki S, Hirano Y, Ogino K, Kitamura A, et al. The AAA+ ATPase/ubiquitin ligase mysterin stabilizes cytoplasmic lipid droplets. J Cell Biol. 2019;218(3):949–60.
Suzuki T, Gao J, Ishigaki Y, Kondo K, Sawada S, Izumi T, et al. ER stress protein CHOP mediates insulin resistance by modulating adipose tissue macrophage polarity. Cell Rep. 2017;18(8):2045–57.
Tang J, Zhuo H, Zhang X, Jiang R, Ji J, Deng L, et al. A novel biomarker Linc00974 interacting with KRT19 promotes proliferation and metastasis in hepatocellular carcinoma. Cell Death Dis. 2014;5:12.
Vinet L, Zhedanov A. A “missing” family of classical orthogonal polynomials. J Phys A Math Theor. 2011;44(8):197–206.
Vishnyakova PA, Tarasova NV, Volodina MA, Tsvirkun DV, Sukhanova IA, Kurchakova TA, et al. Gestation age-associated dynamics of mitochondrial calcium uniporter subunits expression in feto-maternal complex at term and preterm delivery. Sci Rep. 2019;9:1.
Voerman E, Santos S, Inskip H, Amiano P, Barros H, Charles MA, et al. Association of gestational weight gain with adverse maternal and infant outcomes. JAMA. 2019;321(17):1702–15.
Wang Z, Kang J, Deng X, Guo B, Wu B, Fan Y. Knockdown of GATAD2A suppresses cell proliferation in thyroid cancer in vitro. Oncol Rep. 2017;37(4):2147–52.
Warrington NM, Richmond R, Fenstra B, Myhre R, Gaillard R, Paternoster L, et al. Maternal and fetal genetic contribution to gestational weight gain. Int J Obes. 2018;42(4):775–84.
Wasniewski T, Woclawek-Potocka I. Altered expression of lysophosphatidic acid receptors, in association with the synthesis of estrogens and androgens in type 1 endometrial cancer biology. Gynecol Endocrinol. 2018;34(5):422–7.
Wen W, Kato N, Hwang JY, Guo X, Tabara Y, Li H, et al. Genome-wide association studies in East Asians identify new loci for waist-hip ratio and waist circumference. Sci Rep. 2016;6:1.
Xu L, Yang L, Wang L, Zhu B, Chen Y, Gao H, et al. Probe-based association analysis identifies several deletions associated with average daily gain in beef cattle. BMC Genomics. 2019;20:1.
Yu L, Liang Q, Zhang W, Liao M, Wen M, Zhan B, et al. HSP22 suppresses diabetes-induced endothelial injury by inhibiting mitochondrial reactive oxygen species formation. Redox Biol. 2019;21:1.
Yu X, Zhang Y, Chen H. LPA receptor 1 mediates LPA-induced ovarian cancer metastasis: an in vitro and in vivo study. BMC Cancer. 2016;16:1.
Za’don NHA, Kamal AFM, Ismail F, Othman SIT, Appukutty M, Salim N, et al. High-intensity interval training induced PGC-1α and Adipor1 gene expressions and improved insulin sensitivity in obese individuals. Med J Malaysia. 2019;74(6):461–7.
Zhang Y, Yang J, Zhang J, Sun L, Hirankarn N, Pan HF, et al. Genome-wide search followed by replication reveals genetic interaction of CD80 and ALOX5AP associated with systemic lupus erythematosus in asian populations. Ann Rheum Dis. 2016;75(5):891–8.
Zhao LP, Alshiekh S, Zhao M, Carlsson A, Larsson HE, Forsander G, et al. Next-generation sequencing reveals that HLA-DRB3, -DRB4, and -DRB5 may be associated with islet autoantibodies and risk for childhood type 1 diabetes. Diabetes. 2016;65(3):710–8.
Zhu B, Rippe C, Holmberg J, Zeng S, Perisic L, Albinsson S, et al. Nexilin/NEXN controls actin polymerization in smooth muscle and is regulated by myocardin family coactivators and YAP. Sci Rep. 2018;8:1.
Influence of Pregnancy Weight on Maternal and Child Health. Influence of Pregnancy Weight on Maternal and Child Health. 2007.
The authors thank Prof. Katarzyna Cypryk and Dr. Monika Żurawska-Kliś for their help in recruiting patients. This research was supported in part by PLGrid Infrastructure. The bioinformatic analysis was performed using Prometheus (AGH, Krakow, Poland), Michigan Imputaton Server (Michigan, MI, USA) and The Ohio Supercomputer Center (Columbus, OH, USA). For manuscripts with ARIC GWAS data (including HapMap and 1000G imputed data): The Atherosclerosis Risk in Communities study has been funded in whole or in part with Federal funds from the National Heart, Lung, and Blood Institute, National Institutes of Health, Department of Health and Human Services (contract numbers HHSN268201700001I, HHSN268 201700002I, HHSN268201700003I, HHSN268201700004I and HHSN268201700005I), R0 1HL087641, R01HL059367 and R01HL086694; National Human Genome Research Institute contract U01HG004402; and National Institutes of Health contract HHSN268200625226C. The authors thank the staff and participants of the ARIC study for their important contributions. Infrastructure was partly supported by Grant Number UL1RR025005, a component of the National Institutes of Health and NIH Roadmap for Medical Research. The BioMe Biobank Program T2D-GENES Exome Sequencing Study was conducted by the BioMe Biobank Program T2D-GENES Exome Sequencing Study Investigators and supported by the National Institute of Diabetes and Digestive and Kidney Diseases (NIDDK). The data from the BioMe Biobank Program T2D-GENES Exome Sequencing Study reported here were supplied by the Broad Institute and Icahn School of Medicine at Mount Sinai. This manuscript was not prepared in collaboration with Investigators of the BioMe Biobank Program T2D-GENES Exome Sequencing study and does not necessarily reflect the opinions or views of the BioMe Biobank Program T2D-GENES Exome Sequencing study, or the NIDDK.
The study was funded by the National Science Centre in Poland through the Sonata Grant “Search for genetic variants influencing gestational weight gain in type 1 diabetes patients by genome wide association method” to ALS (Nr 2013/11/D/NZ5/03219).
Ethics approval and consent to participate
This study was approved by the Bioethical Committees of Jagiellonian University and Poznan University of Medical Sciences and performed according to the Helsinki Declaration. Written informed consent was collected from all patients.
Consent for publication
All authors contributed to critical revision of the manuscript and approved its publication.
The authors declare they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
The Table presenting the flow of the subject through the analysis process.
: Table S1a. The results of PrediXcan on BMI in T1D cohort in Subcutaneous Adipose Tissue. b. The results of PrediXcan on BMI in T1D cohort in Visceral Adipose Tissue. c. The results of PrediXcan on GWG in T1D cohort in Subcutaneous Adipose Tissue. d. The results of PrediXcan on GWG in T1D cohort in Visceral Adipose Tissue.
: Table S2. The overlap of PrediXcan results in Subcutaneous and Visceral Adipose Tissue in the Giant cohort.
: Table S3a. The results of PrediXcan on Giant cohort on BMI only. b. The results of PrediXcan on Giant cohort on BMI and GWG. c. The results of PrediXcan on Giant cohort on GWG only.
The subanalyses performed for genes in the overlap between GWG and BMI only.
: Figure S1a. The LD analysis for GPN3 gene. b. The LD analysis for PMS2P3 gene. c. The LD analysis for STAG3L1 gene.
: Table S4a.The results of PrediXcan on BMI in T2D cohort in Subcutaneous Adipose Tissue. b. The results of PrediXcan on BMI in T2D cohort in Visceral Adipose Tissue. c. The results of PrediXcan on BMI in ARIC cohort in Subcutaneous Adipose Tissue. d. The results of PrediXcan on BMI in ARIC cohort in Visceral Adipose Tissue.
: Table S5a. The results of PrediXcan on T2D&ARIC cohorts on BMI only. b. The results of PrediXcan on T2D&ARIC cohorts on BMI and GWG. c. The results of PrediXcan on T2D&ARIC cohorts on GWG only.
: Table S6. The list of variants associated with GWG in BMI associated genes in T2D&ARIC cohorts.
: Figure S2. The KEGG analysis on GWG associated variants in BMI associated genes in T2D&ARIC cohorts.
: Table S7. The FUMA analysis on BMI associated variants in GWG associated genes in T1DM cohort.
Go Enrichment analysis on the overlap between BMI and GWG genes in the Giant cohort. b. Go Enrichment analysis on BMI only associated genes in the Giant cohort. c. Go Enrichment analysis on GWG only associated genes in the Giant cohort.
: Figure S3a. Go Enrichment analysis on GWG only associated genes in the Giant cohort. b. Go Enrichment analysis on the overlap between BMI and GWG genes in the Giant cohort.
: Table S9a. Go Enrichment analysis on the overlap between BMI and GWG genes in T2D and ARIC cohorts.b. Go Enrichment analysis on BMI only associated genes in T2D and ARIC cohorts. c. Go Enrichment analysis on GWG only associated genes in T2D and ARIC cohorts.
: Figure S4a. Go Enrichment analysis on GWG only associated genes in T2D and ARIC cohorts. b. Go Enrichment analysis on the overlap between BMI and GWG genes in T2D and ARIC cohorts.
About this article
Cite this article
Ludwig-Słomczyńska, A.H., Seweryn, M.T., Kapusta, P. et al. The transcriptome-wide association search for genes and genetic variants which associate with BMI and gestational weight gain in women with type 1 diabetes. Mol Med 27, 6 (2021). https://doi.org/10.1186/s10020-020-00266-z
- Gestational weight gain (GWG)
- Transcriptomic-wide association study (TWAS)
- Association analysis