Skip to main content

Genetic distance and ancestry proportion modify the association between maternal genetic risk score of type 2 diabetes and fetal growth

Abstract

Background

Maternal genetic risk of type 2 diabetes (T2D) has been associated with fetal growth, but the influence of genetic ancestry is not yet fully understood. We aimed to investigate the influence of genetic distance (GD) and genetic ancestry proportion (GAP) on the association of maternal genetic risk score of T2D (GRST2D) with fetal weight and birthweight.

Methods

Multi-ancestral pregnant women (n = 1,837) from the NICHD Fetal Growth Studies – Singletons cohort were included in the current analyses. Fetal weight (in grams, g) was estimated from ultrasound measurements of fetal biometry, and birthweight (g) was measured at delivery. GRST2D was calculated using T2D-associated variants identified in the latest trans-ancestral genome-wide association study and was categorized into quartiles. GD and GAP were estimated using genotype data of four reference populations. GD was categorized into closest, middle, and farthest tertiles, and GAP was categorized as highest, medium, and lowest. Linear regression analyses were performed to test the association of GRST2D with fetal weight and birthweight, adjusted for covariates, in each GD and GAP category.

Results

Among women with the closest GD from African and Amerindigenous ancestries, the fourth and third GRST2D quartile was significantly associated with 5.18 to 7.48 g (weeks 17–20) and 6.83 to 25.44 g (weeks 19–27) larger fetal weight compared to the first quartile, respectively. Among women with middle GD from European ancestry, the fourth GRST2D quartile was significantly associated with 5.73 to 21.21 g (weeks 18–26) larger fetal weight. Furthermore, among women with middle GD from European and African ancestries, the fourth and second GRST2D quartiles were significantly associated with 117.04 g (95% CI = 23.88–210.20, p = 0.014) and 95.05 g (95% CI = 4.73–185.36, p = 0.039) larger birthweight compared to the first quartile, respectively. The absence of significant association among women with the closest GD from East Asian ancestry was complemented by a positive significant association among women with the highest East Asian GAP.

Conclusions

The association between maternal GRST2D and fetal growth began in early-second trimester and was influenced by GD and GAP. The results suggest the use of genetic GD and GAP could improve the generalizability of GRS.

Background

Observational studies have shown an association between extremes of birthweight and an increased risk of perinatal morbidity and mortality, as well as long-term cardiovascular and metabolic complications such as type 2 diabetes (T2D) [1,2,3,4]. Genetic and non-genetic factors contribute to variations in fetal growth [5,6,7,8,9]. The maternal genome influences fetal growth through a direct impact on the intrauterine environment [8, 9], accounting for 7.6% of the genome-wide heritability of birthweight [9]. Previous studies, largely involving European ancestry individuals, showed that maternal genetic susceptibility to T2D contributes to increased fetal growth and birthweight [9,10,11,12,13]. It is known that the distribution of genetic variants, lifestyle factors, and their interactions can vary by ancestry [14, 15]. However, there is a knowledge gap in the link between maternal genetic susceptibility to T2D and fetal growth across different genetic ancestries.

The association between individual or aggregate effect (i.e., genetic risk score (GRS)) of genetic variants and health outcomes is poorly transferable across diverse ancestry groups [16]. GRST2D based on summary statistics from genome-wide association studies (GWAS) in African Americans was associated with increased birthweight in populations of African ancestry [13]. Likewise, GRST2D based on summary statistics from GWAS involving individuals of European ancestry showed a positive association with fetal weight beginning from 24 gestational week among European American women, but not among non-European American women [12]. The difference in GRST2D effects across the genetic ancestry continuum and the gestational weeks in which GRST2D begins to influence fetal growth are, therefore, unclear for most ancestral populations.

Addressing these research gaps through the use of self-reported race/ethnic groups or discrete genetic ancestry classifications has shortcomings [17]. For example, individuals categorized into the same ancestry group may exhibit variations along the genetic ancestry continuum. Consequently, it is challenging to establish clear boundaries between clusters, and some individuals may be excluded because of difficulty in assigning them to a specific ancestry group [17]. Moreover, the definitions of socially constructed race/ethnic groups and genetic ancestry are intricate and continuously evolving [18].

To address the aforementioned challenges, recent studies suggest using genetic distance (GD) [17, 19]. GD is a continuous metric that quantifies the similarity of an individual’s genome to the center of a standard reference population [17]. GD is considered a less biased estimator of genetic similarity as it takes into account inter-individual genetic variation within a given sample [17]. It has been shown that the accuracy and transferability of GRS derived from associations discovered in European GWAS decrease as the GD from a European ancestry reference population increases [17]. Additionally, the genetic ancestry proportion (GAP) of admixed populations affects the strength of the association between GRS and disease traits in Hispanic Americans [20]. Based on these recent developments, we posit that understanding the relationship between maternal genetic risk for T2D and fetal growth, accounting for GD and GAP, can provide insights about targeted intervention across different population groups. Therefore, in a cohort of multi-ancestral US pregnant women, we investigated the influence of maternal GD and GAP in the association between maternal GRST2D and fetal weight trajectories throughout gestational weeks 10–40.

Methods

Study population

Data from the Eunice Kennedy Shriver National Institute of Child Health and Human Development (NICHD) Fetal Growth Studies – Singletons cohort study, which recruited 2,802 women from four self-identified racial/ethnic groups (i.e., non-Hispanic white, non-Hispanic black, Hispanic white and Asian and Pacific Islanders), was used [21, 22]. The NICHD Fetal Growth Studies – Singletons cohort study was designed to develop a fetal growth standard for the US by race/ethnicity. Details on the inclusion criteria, data collection and quality assurance, and ethical approval were published elsewhere [21, 22]. Women completed sociodemographic, reproductive and pregnancy history questionnaires at enrollment in addition to providing blood specimen [21]. Clinical data was also extracted from medical records. The current analyses included 1,837 women with genotype data, at least two measures of fetal weight, and without gestational diabetes or hyperglycemia based on glucose challenge test.

Fetal weight and birthweight measures

After ultrasonographic confirmation of gestational age at 8+ 0 and 13+ 6 weeks of gestation, participants were randomized to one of four groups and scheduled for five follow-up appointments to measure fetal growth biometrics using identical standardized obstetrical ultrasonography equipment and protocols [21, 22]. This randomization scheme was designed to capture all gestational week’s windows. Measurements between site sonographers and experts had a high correlation (> 0.99) and low coefficient of variation (< 3%) [23]. Fetal weight was estimated from head circumference, abdominal circumference, and femur length using Hadlock’s formula [24]. Fetal weight measures at gestational weeks 10–40 were estimated from the five measurements using a linear mixed model with a cubic spline mean structure and a cubic polynomial random effect [25]. Birthweight was measured in grams (g) using an electronic infant scale or beam balance scale [21].

Genotyping and quality control

Genomic DNA was extracted from stored maternal blood specimens obtained from 2,215 (i.e., 641 self-identified non-Hispanic white, 652 non-Hispanic black, 582 Hispanic white and 340 Asian and Pacific Islander) women and genotyped using the Infinium Multiethnic Global BeadChip microarray (Illumina) with > 1.7 million single nucleotide polymorphisms (SNPs). A total of 2,056 women passed genotype quality control. SNP genotypes were imputed using the entire 1000 Genomes Phase 3 reference sequence data in the Michigan Imputation Server, implementing Eagle2 for haplotype phasing followed by Minimac3 for imputing non-typed SNPs. Details on genotyping, quality control and imputation procedures were described in our previous study [26].

Genetic risk score for T2D (GRST2D)

Publicly available summary statistics of 302 of the 338 genome-wide significant SNPs associated with T2D in the latest T2D multi-ancestry GWAS were used to calculate GRST2D (Supplementary Table 1) [27]. GRS aggregates the effects of genetic variants into a single score [28]. Weighted GRST2D was calculated by multiplying the dosage (i.e., values range 0 to 2) of the T2D-increasing allele for each SNP by logOR and the resulting values summed up [27]. We categorized the GRST2D into quartiles, and the lowest quartile (quartile 1) was used as a reference group.

Genetic distance (GD) and ancestry proportion (GAP)

We used GD and ancestry proportion estimates to determine, for each woman’s genome, similarity with genotype of a reference population [29], and genome-wide average ancestry from each of four continental populations [30], respectively. First, genome-wide principal components analysis (PCA) of our cohort and the reference population (i.e., 1000 Genome and Human Genome Diversity Project population (HGDP)) genotype was performed using ‘flashpcaR’ R package. Next, our cohort genotype was projected on the PC space of the reference population genotype. Unlike the previous studies that estimated GD using only European reference population [17, 19], we performed four types of GD estimations using three reference populations from the 1000 Genomes project [31]: European - Utah residents with Northern and Western European ancestry (CEU); African - Yoruba in Ibadan, Nigeria (YRI); and East Asian - Han Chinese in Beijing, China (CHB); and one reference population from the 1000 Genomes and Human Genome Diversity Project (HGDP) [32] – Amerindigenous – Native American (NAM). Then, GD of each woman’s genotype from the center of the reference population genotype projected on the PC space was calculated [17]. Finally, women were categorized into tertiles based on their GD and conventionally named as ‘closest GD’ (tertile 1), ‘middle GD’ (tertile 2) or ‘farthest GD’ (tertile 3) (Supplementary Table 2).

To supplement our findings based on GD, GAP was estimated for each sample by performing linkage disequilibrium pruning of the genotype data and running unsupervised clustering in ADMIXTURE version 1.3 software [30] with the assumption of four ancestral components (k = 4). These clusters represent European, African, Amerindigenous and East Asian genetic ancestry proportions. Reference samples for African, European and East Asian ancestries were obtained from genotype data of the 1000 Genomes project [31], whereas reference samples for Amerindigenous ancestry were obtained by combining samples from the 1000 Genomes and HGDP [31, 32]. The analyses were described in detail in our previous study [33]. Women were categorized into tertile based on their GAP and conventionally named as ‘lowest GAP’ (tertile 1), ‘medium GAP’ (tertile 2), and ‘highest GAP’ (tertile 3) (Supplementary Table 3).

Data analyses

Maternal and newborn background characteristics were summarized using frequency (%) and mean (± SD). The correlation between GD and GAP was estimated using Spearman correlation method. Linear regression models were fitted to test the association of GRST2D with fetal weight at each gestation week (weeks 10–40) and birthweight, adjusted for years lived in the United States (continuous), marital status (categorical), maternal age (continuous), pre-pregnancy body mass index (continuous), parity (categorical), and fetal sex (categorical) in each GD and GAP category. The association between GRST2D and birthweight was further adjusted for gestational age at delivery (continuous). Additionally, interaction model was fitted by including interaction between GRST2D (continuous) and GD and GAP (continuous). Statistical significance was set to be p-value < 0.05. All analyses, unless specified otherwise, were implemented using the software package R version 4.2.2.

Results

Characteristics of study participants

The mean (± SD) of maternal age, number of years lived in the US, and birthweight were 28 (± 6) years, 23 (± 9) years, and 3,314 (± 529) grams, respectively. Approximately three-fourths of the women were married or living with a partner (73.54%). The mean (± SD) GRST2D was 16.19 (± 0.55). The mean (± SD) GD of the women from European, African, Amerindigenous, and East Asian reference panels was 0.03 (± 0.01), 0.03 (± 0.01), 0.13 (± 0.02), and 0.04 (± 0.01) unit, respectively. The mean (± SD) European, African, Amerindigenous, and East Asian GAPs in the maternal genome were 0.48 (± 0.37), 0.30 (± 0.36), 0.10 (± 0.20), and 0.12 (± 0.29), respectively (Table 1).

Table 1 Description of maternal and fetal background characteristics. (N = 1,837)

There was weak to moderate negative correlation between GD from an ancestry reference and GAP of the same ancestry, for instance, with the lowest correlation between African GD and European GD (r = -0.0.06, p < 0.05) and the highest correlation between European GD and European GAP (r = -0.56, p < 0.05) (Fig. 1). Nearly all GD and GAP tertile groups included each of the four self-reported race/ ethnicity. No GD or GAP tertile represents less than three race/ethnicity groups, demonstrating the continuum nature of human population variation (Supplementary Figs. 1 and 2).

Fig. 1
figure 1

Correlation between genetic distance (GD) and genetic ancestry proportion (GAP). Numbers inside boxes denote correlation coefficients and are only shown for those with p < 0.05

GD and association of GRST2D with fetal weight

GRST2D showed significant positive association with fetal weight beginning from week 17 among women with the closest African GD, medium European GD, and irrespective of Amerindigenous GD. Specifically, the fourth GRST2D quartile was significantly associated with 5.73 to 21.21 g (weeks 18–26) and 5.18 to 7.48 g (weeks 17–20) larger fetal weight compared to the first quartile among women with the middle GD from European ancestry and with the closest GD from African ancestry, respectively. Among women with the closest and farthest GD from Amerindigenous ancestry, the third and second GRST2D quartiles were significantly associated with 6.83 to 25.44 g (weeks 19–27) and 7.65 to 149.71 g (weeks 19–40) larger fetal weight compared to the first quartile, respectively. Among women with middle GD from Amerindigenous ancestry, the second GRST2D quartile was significantly associated with 8.75 to 69.75 g (weeks 21–35) lower fetal weight. No significant associations were observed regardless of GD from East Asian ancestry (Fig. 2, Supplementary Fig. 3, Supplementary Table 4). The interaction model showed that the interaction between European GD and GRST2D significantly associated with fetal weight (weeks 21–40, p < 0.05).

Fig. 2
figure 2

Weekly change in standardized fetal weight (in SD) associated with GRST2D based on GD. The dotted curved lines represent the average change in fetal weight, while the broken lines with colored zones indicate the 95% CI. The solid horizontal lines indicate the null hypothesis of no change in fetal weight. GRS (Q1) is the reference group

GD and association of GRST2D with birthweight

GRST2D and birthweight showed significant positive association among women with medium and farthest European GD, medium African GD, and farthest Amerindigenous GD. Specifically, among women with middle and farthest GD from European ancestry, the fourth and second GRST2D quartile was significantly associated with a 117.04 g (95% CI = 23.88–210.20, p = 0.014) and 110.02 g (95% CI = 17.09–202.95, p = 0.02) larger birthweight compared to the first quartile, respectively. Compared to the first quartile, the second GRST2D quartile was significantly associated with a 95.05 g (95% CI = 4.73–185.36, p = 0.039) and 95.28 g (95% CI = 4.54–186.03, p = 0.04) larger birthweight among women with middle GD from African ancestry and with farthest GD from Amerindigenous ancestry, respectively. The interaction model showed that the interaction between European GD and GRST2D significantly associated with birthweight (p = 0.04). No significant association was observed regardless of GD from East Asian ancestry (Fig. 3).

Fig. 3
figure 3

Change in birthweight (g) associated with GRST2D based on GD. The squares indicate the average change in birthweight, the horizontal lines represent the 95% CI, and the broken vertical lines indicate the null hypothesis of no change. GRS (Q1) is the reference group

GAP and association of GRST2D with fetal weight

GRST2D showed significant positive associations with fetal weight among women with medium European GAP, highest Amerindigenous GAP, and highest East Asian GAP. Specifically, among women with medium European GAP, the second GRST2D quartile was significantly associated with 7.14 to 21.63 g (weeks 19–26) larger fetal weight compared to the lowest quartile. Among women with the highest Amerindigenous GAP, the fourth GRST2D quartile was significantly associated with 3.48 to 3.65 g (weeks 10–12) smaller fetal weight and 25.35 to 117.53 g (weeks 27–40) larger fetal weight. Among women with the highest East Asian GAP, the fourth GRST2D quartile was significantly associated with 5.93 to 13.26 g (weeks 18–23) larger fetal weight. No significant association was found regardless of African GAP (Fig. 4, Supplementary Fig. 4, Supplementary Table 5). The interaction model showed that the interaction between Amerindigenous GAP and GRST2D significantly associated with fetal weight (weeks 25–40, p < 0.05).

Fig. 4
figure 4

Weekly change in fetal weight (in SD) associated with GRST2D based on GAP. The dotted curved lines represent the average change in fetal weight, the broken lines with colored zones represent the 95% CI, and the solid horizontal lines represent the null hypothesis of no change in fetal weight. GRS (Q1) is the reference group

GAP and association of GRST2D with birthweight

GRST2D and birthweight showed significant positive association among women with medium European GAP, medium African GAP, highest East Asian GAP, and highest Amerindigenous GAP. Specifically, among women with medium European GAP, the second, third, and fourth GRST2D quartiles were significantly associated with 95.25 g (95% CI = 4.04–186.46, p = 0.041), 131.11 g (95% CI = 40.39–221.83, p = 0.005), and 102.59 g (95% CI = 11.47–193.71, p = 0.027) larger birthweight compared to the lowest quartile, respectively. Compared to the lowest quartile, the fourth GRST2D quartile was significantly associated with 111.88 g (95% CI = 18.03–205.73, p = 0.02) among women with medium African GAP; and the third GRST2D quartile was significantly associated with 106.45 g (95% CI = 15.48–197.42, p = 0.022) larger birthweight among women with the highest East Asian GAP. Among women with the highest Amerindigenous GAP, the third and fourth GRST2D quartile was significantly associated with 126.54 g (95% CI = 35.04–218.04, p = 0.007) and 152.17 g (95% CI = 59.42–244.92, p = 0.001) larger birthweight, respectively (Fig. 5). The interaction model showed significant interaction between Amerindigenous GAP and GRST2D (p = 0.004).

Fig. 5
figure 5

Change in birthweight (g) associated with GRST2D based on GAP. Squares represents the average change in birthweight, the horizontal lines represent the 95% CI, and the broken vertical lines represent the null hypothesis of no change. GRS (Q1) is the reference group

Discussion

By leveraging T2D-associated loci from a large multi-ancestral GWAS and incorporating genetic distance and ancestry proportion estimates in a multi-ancestral pregnancy cohort, we found that maternal GRST2D begins to influence fetal weight in early second trimester of pregnancy. Presence and gestational timing of associations between GRST2D and fetal weight varied by maternal GD and GAP. In many instances, the association between maternal GRST2D and fetal weight was influenced by both GD and GAP. For example, GRST2D was consistently associated with higher birthweight in women with middle or farthest GD from European, African, and Amerindigenous ancestries, as well as in those with the highest GAP. In other instances, the use of GAP yielded an association undetected when using GD. For example, while GRST2D was not associated with fetal weight or birthweight irrespective of GD from East Asian ancestry, significant association was detected with East Asian GAP. These findings and the weak to moderate correlation between GD and GA suggest the importance of considering both GD and GAP as alternative measures in identifying the timing and impact of GRST2D on fetal weight.

A key finding of our study is the detection of associations between maternal GRST2D and fetal weight in early second trimester (gestational weeks 17–19) based on GD from African, European, and Amerindigenous ancestry. In a previous study, maternal GRST2D was associated with increased fetal weight starting at week 24 among European American women, but not among non-European women [12]. In addition, maternal glucose levels measured at gestational weeks 10–14 have been significantly associated with fetal weight starting at week 27 [34]. The use of genetic distance estimates and implementation of GRS derived from a multi-ancestral GWAS may have facilitated our ability to identify a change in fetal weight at gestational weeks much earlier than previous studies.

There may be several factors that underlie the observed differences in the association of GRST2D with fetal growth by GD and GAP, especially after adjustment for maternal pre-pregnancy body mass index. Other cardiometabolic factors may have trajectories during pregnancy that differ across GD and GAP group, hence partly explaining these differences. For example, total gestational weight gain significantly differed across the continuum of Amerindigenous GD, systolic blood pressure differed across European and African GD, and diastolic blood pressure differed across European and Amerindigenous GD (Supplementary Table 2). Furthermore, out of the 302 SNPs included in the GRS calculation, 40% (121 out of 302) were associated with at least one cardiometabolic or anthropometric trait such as lipid profiles, blood pressure, adiposity and body composition, over-weight or obesity, and birth weight (Supplementary Table 1). Interactions among these factors and shared genetic effects could set the intrauterine metabolic, endocrine, and inflammation environment, potentially explaining differential effects of GRST2D on fetal growth by GD and GAP groups. Further investigation of the role of T2D-related SNPs in pregnancy-related cardiometabolic traits such as gestational weight gain, gestational diabetes, and gestational hypertension could provide additional insights into the association between GRST2D and fetal growth.

The mechanisms through which maternal T2D genetic risk induces fetal growth in early pregnancy are complex. A possible mechanism of the genetic effect is enhanced trans-placental glucose transfer which leads to increased secretion of fetal insulin, a growth hormone that promotes fetal growth [8, 9, 35]. Our finding is supported by this mechanism because fetal pancreatic release of insulin begins as early as 11 weeks of gestation [36]. Further studies should clarify the relationship between fetal insulin and glucose levels during early gestation. Moreover, non-glycemic pathways may underlie our observed associations because T2D shares genetic risk with metabolic and cardiovascular phenotypes, all of which are related to fetal growth [12]. Given the interdependent influence of the parental and fetal genomes on birthweight [8, 9], integrating fetal and parental genotype data will help elucidate the potential mechanisms underlying these associations [37,38,39,40].

Using multiple ancestry reference genotypes as anchors, we generated four GD estimates that allowed evaluation of GRS performance accounting for broad human population genetic variations. Previous studies tested the performance of GRS based on only GD estimates using a European ancestry reference genotype [17, 19], which could disregard the complex picture in human genetic variation and sociocultural differences. Using GD based on a single reference population could undermine predicting genetic influence in multi-ancestral populations that lead to increased heath disparity. For instance, although GRST2D was associated with fetal weight in the same gestational period for both the group with the closest GD from African and Amerindigenous references, the composition of those two groups by race/ethnicity is different. Even within a single ancestry group, the prediction accuracy of GRS may vary depending on factors such as socio-economic status, geographic distance, age, or sex of the individuals involved in the GWAS and prediction cohorts [15]. Moreover, race/ethnicity and genetic ancestry are complex constructs with evolving definitions due to the social nature of self-identified racial/ethnic constructs and phenotypic heterogeneity within broad racial/ethnic groups [12, 18]. Therefore, the use of multiple GD can partly overcome this limitation and lead to improved accuracy in prediction of genetic influence in non-European populations.

Although the GRST2D used in the present analysis was based on a T2D GWAS involving diverse populations, individuals of European ancestry form the vast majority of the study participants [27], which is a prevailing limitation of many genetic studies. Including ancestrally diverse individuals in genetic studies enhances insights into the genetic architecture of diseases [41]. Addressing the diversity gap in genomics is critical to ensure that all communities benefit from research innovations and precision medicine [42,43,44]. GRS derived from ancestry-matched and multi-ancestry studies will have better predictive accuracy, facilitating translation of genetic research into clinical care and public health policies that benefit a wide range of populations [45, 46].

Our study has several strengths. The GRS was derived from summary statistics of the largest T2D-GWAS in diverse ancestry groups [27]. The longitudinal ultrasound measures of fetal biometry in the NICHD Fetal Growth Studies allowed a comprehensive understanding of the influence of genetic variants on fetal growth at specific time points during pregnancy and at birth [21]. The study recruited low-risk pregnant women without previous history of major medical conditions and pregnancy complications, thereby minimizing potential confounding due to maternal morbidity [21].

We acknowledge certain limitations. The fetal effects of the SNPs used to calculate maternal GRS have not been accounted because fetal genotypes are unavailable in our dataset. Several studies consistently reported T2D-risk variants in the maternal and fetal genomes may have directionally opposite effects on insulin-mediated fetal growth [10, 47]. By design, women with type 2 diabetes were not enrolled in the study, limiting generalizability of the findings to women with glycemic dysregulations. Given the exploratory nature of our study, we did not apply multiple test correction. The Hadlock formula was used to calculate fetal weight, which may introduce inaccuracies particularly in the third trimester [24]. To address this issue, the average of three measurements was used, and the same sonographers were involved throughout the study period, minimizing intra-observer and inter-observer variability [23]. While GRS holds promise as a potential biomarker in medicine [48], it is crucial to validate its performance using studies with larger sample sizes and by integrating with non-genetic factors.

Conclusions

Our study showed that the association between GRST2D and fetal growth begins in early second trimester of pregnancy and may vary depending on the reference population used to infer genetic similarity. Additionally, the use of both GD and GAP yielded valuable and alternative insights that may not be obtained when using either alone. The research question, study population, non-genetic factors, and their relationship with the chosen reference should be considered when selecting a more appropriate measure of genetic similarity. Further investigation is needed to evaluate the extent of maternal genetic susceptibility to T2D, integrated with other clinical and social determinants of health, can facilitate early detection of fetal growth changes, thereby mitigating adverse pregnancy outcomes and long-term diseases.

Data availability

The data sets generated and/or analyzed during the current study are available from the corresponding author upon reasonable request.

Abbreviations

BMI:

Body mass index

GAP:

Genetic ancestry proportion

GD:

Genetic distance

GRS:

Genetic risk score

GRST2D :

Genetic risk score of type 2 diabetes

GWAS:

Genome-wide association study

HGDP:

Human genome diversity project

NICHD:

National Institute of Child Health and Human Development

PCA:

Principal component analysis

SD:

Standard deviation

SNPs:

Single nucleotide polymorphisms

T2D:

Type 2 diabetes

References

  1. Risnes KR, Vatten LJ, Baker JL, Jameson K, Sovio U, Kajantie E, et al. Birthweight and mortality in adulthood: a systematic review and meta-analysis. Int J Epidemiol. 2011;40(3):647–61.

    Article  PubMed  Google Scholar 

  2. Mi D, Fang H, Zhao Y, Zhong L. Birth weight and type 2 diabetes: a meta-analysis. Exp Ther Med. 2017;14(6):5313–20.

    PubMed  PubMed Central  Google Scholar 

  3. Knop MR, Geng TT, Gorny AW, Ding R, Li C, Ley SH, et al. Birth weight and risk of type 2 diabetes Mellitus, Cardiovascular Disease, and hypertension in adults: a Meta-analysis of 7 646 267 participants from 135 studies. J Am Heart Association. 2018;7(23):e008870.

    Article  Google Scholar 

  4. Harder T, Rodekamp E, Schellong K, Dudenhausen JW, Plagemann A. Birth weight and subsequent risk of type 2 diabetes: a meta-analysis. Am J Epidemiol. 2007;165(8):849–57.

    Article  PubMed  Google Scholar 

  5. Zheng T, Zhang J, Sommer K, Bassig BA, Zhang X, Braun J, et al. Effects of Environmental exposures on fetal and Childhood Growth trajectories. Ann Glob Health. 2016;82(1):41–99.

    Article  PubMed  Google Scholar 

  6. Hoffman DJ, Reynolds RM, Hardy DB. Developmental origins of health and disease: current knowledge and potential mechanisms. Nutr Rev. 2017;75(12):951–70.

    Article  PubMed  Google Scholar 

  7. Gluckman PD, Buklijas T, Hanson MA. Chapter 1 - the Developmental origins of Health and Disease (DOHaD) Concept: past, Present, and Future. In: Rosenfeld CS, editor. The Epigenome and Developmental origins of Health and Disease. Boston: Academic; 2016. pp. 1–15.

    Google Scholar 

  8. Juliusdottir T, Steinthorsdottir V, Stefansdottir L, Sveinbjornsson G, Ivarsdottir EV, Thorolfsdottir RB, et al. Distinction between the effects of parental and fetal genomes on fetal growth. Nat Genet. 2021;53(8):1135–42.

    Article  CAS  PubMed  Google Scholar 

  9. Warrington NM, Beaumont RN, Horikoshi M, Day FR, Helgeland Ø, Laurin C, et al. Maternal and fetal genetic effects on birth weight and their relevance to cardio-metabolic risk factors. Nat Genet. 2019;51(5):804–14.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  10. Beaumont RN, Warrington NM, Cavadino A, Tyrrell J, Nodzenski M, Horikoshi M, et al. Genome-wide association study of offspring birth weight in 86 577 women identifies five novel loci and highlights maternal genetic effects that are independent of fetal genetics. Hum Mol Genet. 2018;27(4):742–56.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  11. Freathy RM, Weedon MN, Bennett A, Hypponen E, Relton CL, Knight B, et al. Type 2 diabetes TCF7L2 risk genotypes alter birth weight: a study of 24,053 individuals. Am J Hum Genet. 2007;80(6):1150–61.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  12. Ouidir M, Zeng X, Chatterjee S, Zhang C, Tekola-Ayele F. Ancestry-matched and Cross-ancestry Genetic Risk scores of type 2 diabetes in pregnant women and fetal growth: a study in an Ancestrally Diverse Cohort. Diabetes. 2022;71(2):340–9.

    Article  CAS  PubMed  Google Scholar 

  13. Rahman ML, Shrestha D, Workalemahu T, Wu J, Zhu C, Zhang C, et al. Maternal and offspring genetic risk of type 2 diabetes and offspring birthweight among African ancestry populations. J Clin Endocrinol Metab. 2019;104(11):5032–42.

    Article  PubMed  PubMed Central  Google Scholar 

  14. McVean GA, Altshuler DM, Durbin RM, Abecasis GR, Bentley DR, Chakravarti A, et al. An integrated map of genetic variation from 1,092 human genomes. Nature. 2012;491(7422):56–65.

    Article  CAS  PubMed  Google Scholar 

  15. Mostafavi H, Harpak A, Agarwal I, Conley D, Pritchard JK, Przeworski M. Variable prediction accuracy of polygenic scores within an ancestry group. eLife. 2020;9:e48376.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  16. Wang Y, Tsuo K, Kanai M, Neale BM, Martin AR. Challenges and opportunities for developing more generalizable polygenic risk scores. Annu Rev Biomed Data Sci. 2022;5:293–320.

    Article  PubMed  PubMed Central  Google Scholar 

  17. Ding Y, Hou K, Xu Z, Pimplaskar A, Petter E, Boulier K, et al. Polygenic scoring accuracy varies across the genetic ancestry continuum. Nature. 2023;618(7966):774–81.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  18. Lewis ACF, Molina SJ, Appelbaum PS, Dauda B, Di Rienzo A, Fuentes A, et al. Getting genetic ancestry right for science and society. Science. 2022;376(6590):250–2.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  19. Privé F, Aschard H, Carmi S, Folkersen L, Hoggart C, O’Reilly PF, et al. Portability of 245 polygenic scores when derived from the UK Biobank and applied to 9 ancestry groups from the same cohort. Am J Hum Genet. 2022;109(1):12–23.

    Article  PubMed  PubMed Central  Google Scholar 

  20. Sharma J. Investigating the influence of genetic ancestry on gene-environment interactions of polygenic risk score and acculturation in Hispanic/Latino groups. Johns Hopkins University; 2022.

  21. Grewal J, Grantz KL, Zhang C, Sciscione A, Wing DA, Grobman WA, et al. Cohort Profile: NICHD fetal growth studies-singletons and twins. Int J Epidemiol. 2018;47(1):25–l.

    Article  PubMed  Google Scholar 

  22. Buck Louis GM, Grewal J, Albert PS, Sciscione A, Wing DA, Grobman WA, et al. Racial/ethnic standards for fetal growth: the NICHD fetal growth studies. Am J Obstet Gynecol. 2015;213(4):449. e1-.e41.

    Article  PubMed Central  Google Scholar 

  23. Hediger ML, Fuchs KM, Grantz KL, Grewal J, Kim S, Gore-Langton RE, et al. Ultrasound Quality Assurance for Singletons in the National Institute of Child Health and Human Development fetal growth studies. J Ultrasound Med. 2016;35(8):1725–33.

    Article  PubMed  PubMed Central  Google Scholar 

  24. Hadlock FP, Harrist RB, Sharman RS, Deter RL, Park SK. Estimation of fetal weight with the use of head, body, and femur measurements–a prospective study. Am J Obstet Gynecol. 1985;151(3):333–7.

    Article  CAS  PubMed  Google Scholar 

  25. Pinheiro J, Bates D. Mixed-effects models in S and S-PLUS. Springer science & business media; 2006.

  26. Tekola-Ayele F, Zhang C, Wu J, Grantz KL, Rahman ML, Shrestha D, et al. Trans-ethnic meta-analysis of genome-wide association studies identifies maternal ITPR1 as a novel locus influencing fetal growth during sensitive periods in pregnancy. PLoS Genet. 2020;16(5):e1008747.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  27. Mahajan A, Spracklen CN, Zhang W, Ng MCY, Petty LE, Kitajima H, et al. Multi-ancestry genetic study of type 2 diabetes highlights the power of diverse populations for discovery and translation. Nat Genet. 2022;54(5):560–72.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  28. Igo RP Jr., Kinzy TG, Cooke Bailey JN. Genetic risk scores. Curr Protoc Hum Genet. 2019;104(1):e95.

    Article  PubMed  PubMed Central  Google Scholar 

  29. National Academies of Sciences E. Medicine. Using population descriptors in genetics and genomics research: A new framework for an evolving field2023.

  30. Alexander DH, Novembre J, Lange K. Fast model-based estimation of ancestry in unrelated individuals. Genome Res. 2009;19(9):1655–64.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  31. Siva N. 1000 genomes project. Nat Biotechnol. 2008;26(3):256–7.

    Article  PubMed  Google Scholar 

  32. Cann HM. Human genome diversity. C R Acad Sci III. 1998;321(6):443–6.

    Article  CAS  PubMed  Google Scholar 

  33. Tekola-Ayele F, Ouidir M, Shrestha D, Workalemahu T, Rahman ML, Mendola P, et al. Admixture mapping identifies African and amerindigenous local ancestry loci associated with fetal growth. Hum Genet. 2021;140(7):985–97.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  34. Li M, Hinkle SN, Grantz KL, Kim S, Grewal J, Grobman WA, et al. Glycaemic status during pregnancy and longitudinal measures of fetal growth in a multi-racial US population: a prospective cohort study. Lancet Diabetes Endocrinol. 2020;8(4):292–300.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  35. Tyrrell J, Richmond RC, Palmer TM, Feenstra B, Rangarajan J, Metrustry S, et al. Genetic evidence for Causal relationships between maternal obesity-related traits and Birth Weight. JAMA. 2016;315(11):1129–40.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  36. Reiher H, Fuhrmann K, Noack S, Woltanski KP, Jutzi E, Hahn von Dorsche H, et al. Age-dependent insulin secretion of the endocrine pancreas in vitro from fetuses of diabetic and nondiabetic patients. Diabetes Care. 1983;6(5):446–51.

    Article  CAS  PubMed  Google Scholar 

  37. Vaessen N, Janssen JA, Heutink P, Hofman A, Lamberts SW, Oostra BA, et al. Association between genetic variation in the gene for insulin-like growth factor-I and low birthweight. Lancet. 2002;359(9311):1036–7.

    Article  CAS  PubMed  Google Scholar 

  38. Lindsay RS, Hanson RL, Wiedrich C, Knowler WC, Bennett PH, Baier LJ. The insulin gene variable number tandem repeat class I/III polymorphism is in linkage disequilibrium with birth weight but not type 2 diabetes in the Pima population. Diabetes. 2003;52(1):187–93.

    Article  CAS  PubMed  Google Scholar 

  39. Chen J, Bacelis J, Sole-Navais P, Srivastava A, Juodakis J, Rouse A, et al. Dissecting maternal and fetal genetic effects underlying the associations between maternal phenotypes, birth outcomes, and adult phenotypes: a mendelian-randomization and haplotype-based genetic score analysis in 10,734 mother-infant pairs. PLoS Med. 2020;17(8):e1003305.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  40. Tyrrell JS, Yaghootkar H, Freathy RM, Hattersley AT, Frayling TM. Parental diabetes and birthweight in 236 030 individuals in the UK biobank study. Int J Epidemiol. 2013;42(6):1714–23.

    Article  PubMed  PubMed Central  Google Scholar 

  41. Caliebe A, Tekola-Ayele F, Darst BF, Wang X, Song YE, Gui J, et al. Including diverse and admixed populations in genetic epidemiology research. Genet Epidemiol. 2022;46(7):347–71.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  42. Landry LG, Ali N, Williams DR, Rehm HL, Bonham VL. Lack of diversity in genomic databases is a barrier to translating precision medicine research into practice. Health Aff. 2018;37(5):780–5.

    Article  Google Scholar 

  43. Fine MJ, Ibrahim SA, Thomas SB. The role of race and genetics in health disparities research. Am J Public Health. 2005;95(12):2125–8.

    Article  PubMed  PubMed Central  Google Scholar 

  44. Konuma T, Okada Y. Statistical genetics and polygenic risk score for precision medicine. Inflamm Regen. 2021;41(1):18.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  45. Martin AR, Kanai M, Kamatani Y, Okada Y, Neale BM, Daly MJ. Clinical use of current polygenic risk scores may exacerbate health disparities. Nat Genet. 2019;51(4):584–91.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  46. Sirugo G, Williams SM, Tishkoff SA. The missing diversity in Human Genetic studies. Cell. 2019;177(1):26–31.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  47. Hattersley AT, Beards F, Ballantyne E, Appleton M, Harvey R, Ellard S. Mutations in the glucokinase gene of the fetus result in reduced birth weight. Nat Genet. 1998;19(3):268–70.

    Article  CAS  PubMed  Google Scholar 

  48. Wand H, Lambert SA, Tamburro C, Iacocca MA, O’Sullivan JW, Sillari C, et al. Improving reporting standards for polygenic scores in risk prediction studies. Nature. 2021;591(7849):211–9.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

Download references

Acknowledgements

We acknowledge the study participants of the NICHD Fetal Growth Studies. We thank research teams at all participating clinical centers (which include Christina Care Health Systems, Columbia University, Fountain Valley Hospital, California, Long Beach Memorial Medical Center, New York Hospital, Queens, Northwestern University, University of Alabama at Birmingham, University of California, Irvine, Medical University of South Carolina, Saint Peters University Hospital, Tufts University, and Women and Infants Hospital of Rhode Island). The authors also acknowledge the Wadsworth Center, C-TASC and The EMMES Corporations in providing data and imaging support. This work utilized the computational resources of the NIH HPC Biowulf cluster (http://hpc.nih.gov).

Funding

This research was supported by the Intramural Research Program of the Eunice Kennedy Shriver National Institute of Child Health and Human Development, National Institutes of Health including American Recovery and Reinvestment Act funding via contract numbers HHSN275200800013C; HHSN275200800002I; HHSN27500006; HHSN275200800003IC; HHSN275200800014C; HHSN275200800012C; HHSN275200800028C; HHSN275201000009C and HHSN27500008. Additional support was obtained from the NIH Office of the Director, the National Institute on Minority Health and Health Disparities and the National Institute of Diabetes and Digestive and Kidney Diseases.

Open access funding provided by the National Institutes of Health

Author information

Authors and Affiliations

Authors

Contributions

F.T-A. conceived and designed the study, supervised the analyses, and critically reviewed and edited the manuscript. R.J.B and T.D. performed data analyses and interpretation. T.D. and P.W. wrote the manuscript draft. All authors provided critical intellectual content and approved the final manuscript.

Corresponding author

Correspondence to Fasil Tekola-Ayele.

Ethics declarations

Ethics approval and consent to participate

The study was approved by the Institutional Review Boards of NICHD and all participating clinical sites. Informed consent was obtained from each of the study participants. The study has been registered at ClinicalTrials.gov (Trial registration: NCT00912132).

Consent for publication

Not applicable.

Competing interests

The authors declare no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Electronic supplementary material

Below is the link to the electronic supplementary material.

Supplementary Material 1

Supplementary Material 2

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Habtewold, T.D., Wijesiriwardhana, P., Biedrzycki, R.J. et al. Genetic distance and ancestry proportion modify the association between maternal genetic risk score of type 2 diabetes and fetal growth. Hum Genomics 18, 81 (2024). https://doi.org/10.1186/s40246-024-00645-1

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s40246-024-00645-1

Keywords