Novel loci and Mapuche genetic ancestry are associated with pubertal growth traits in Chilean boys

Joint models reveal genetic architecture of pubertal stage transitions and their association with BMI in admixed Chilean population

Early or late pubertal onset can lead to disease in adulthood, including cancer, obesity, type 2 diabetes, metabolic disorders, bone fractures, and psychopathologies. Thus, knowing the age at which puberty is attained is crucial as it can serve as a risk factor for future diseases. Pubertal development is divided into five stages of sexual maturation in boys and girls according to the standardized Tanner scale. We performed genome-wide association studies (GWAS) on the "Growth and Obesity Chilean Cohort Study" cohort composed of admixed children with mainly European and Native American ancestry. Using joint models that integrate time-to-event data with longitudinal trajectories of body mass index (BMI), we identified genetic variants associated with phenotypic transitions between pairs of Tanner stages. We identified $42$ novel significant associations, most of them in boys. The GWAS on Tanner $3\rightarrow 4$ transition in boys captured an association peak around the growth-related genes LARS2 and LIMD1 genes, the former of which causes ovarian dysfunction when mutated. The associated variants are expression and splicing Quantitative Trait Loci regulating gene expression and alternative splicing in multiple tissues. Further, higher individual Native American genetic ancestry proportions predicted a significantly earlier puberty onset in boys but not in girls. Finally, the joint models identified a longitudinal BMI parameter significantly associated with several Tanner stages' transitions, confirming the association of BMI with pubertal timing.

Genetic associations with disease in populations with Indigenous American ancestries

The genetic architecture of complex diseases affecting populations with Indigenous American ancestries is poorly understood due to their underrepresentation in genomics studies. While most of the genetic diversity associated with disease trait variation is shared among worldwide populations, a fraction of this component is expected to be unique to each continental group, including Indigenous Americans. Here, I describe the current state of knowledge from genome-wide association studies on Indigenous populations, as well as non-Indigenous populations with partial Indigenous ancestries from the American continent, focusing on disease susceptibility and anthropometric traits. While some studies identified risk alleles unique to Indigenous populations, their effects on trait variation are mostly small. I suggest that the associations rendered by many inter-population studies are probably inflated due to the absence of socio-cultural-economic covariates in the association models. I encourage the inclusion of admixed individuals in future GWAS studies to control for inter-ancestry differences in environmental factors. I suggest that some complex diseases might have arisen as trade-off costs of adaptations to past evolutionary selective pressures. Finally, I discuss how expanding panels with Indigenous ancestries in GWAS studies is key to accurately assess genetic risk in populations from the American continent, thus decreasing global health disparities.

Novel signals and polygenic score for height are associated with pubertal growth traits in Southwestern American Indians

Most genetic variants associated with adult height have been identified through large genome-wide association studies (GWASs) in European-ancestry cohorts. However, it is unclear how these variants influence linear growth during adolescence. This study uses anthropometric and genotypic data from a longitudinal study conducted in an American Indian community in Arizona between 1965-2007. Growth parameters (i.e. height, velocity, and timing of growth spurt) were derived from the Preece-Baines growth model, a parametric growth curve fitted to longitudinal height data, in 787 participants with height measurements spanning the whole period of growth. Heritability estimates suggested that genetic factors could explain 25% to 71% of the variance of pubertal growth traits. We performed a GWAS of growth parameters, testing their associations with 5 077 595 imputed or directly genotyped variants. Six variants associated with height at peak velocity (P < 5 × 10-8, adjusted for sex, birth year and principal components). Implicated genes include NUDT3, previously associated with adult height, and PACSIN1. Two novel variants associated with duration of growth spurt (P < 5 × 10-8) in LOC105375344, an uncharacterized gene with unknown function. We finally examined the association of growth parameters with a polygenic score for height derived from 9557 single nucleotide polymorphisms (SNPs) identified in the GIANT meta-analysis for which genotypic data were available for the American Indian study population. Height polygenic score was correlated with the magnitude and velocity of height growth that occurred before and at the peak of the adolescent growth spurt, indicating overlapping genetic architecture, with no influence on the timing of adolescent growth.

Genetic determinants of serum bilirubin using inferred native American gene variants in Chilean adolescents

Gene variants in the UGT1A1 gene are strongly associated with circulating bilirubin levels in several populations, as well as other variants of modest effect across the genome. However, the effects of such variants are unknown regarding the Native American ancestry of the admixed Latino population. Our objective was to assess the Native American genetic determinants of serum bilirubin in Chilean admixed adolescents using the local ancestry deconvolution approach. We measured total serum bilirubin levels in 707 adolescents of the Chilean Growth and Obesity Cohort Study (GOCS) and performed high-density genotyping using the Illumina-MEGA array (>1.7 million genotypes). We constructed a local ancestry reference panel with participants from the 1000 Genomes Project, the Human Genome Diversity Project, and our GOCS cohort. Then, we inferred and isolated haplotype tracts of Native American, European, or African origin to perform genome-wide association studies. In the whole cohort, the rs887829 variant and others near UGT1A1 were the unique signals achieving genome-wide statistical significance (b = 0.30; p = 3.34 × 10-57). After applying deconvolution methods, we found that significance is also maintained in Native American (b = 0.35; p = 3.29 × 10-17) and European (b = 0.28; p = 1.14 × 10-23) ancestry components. The rs887829 variant explained a higher percentage of the variance of bilirubin in the Native American (37.6%) compared to European ancestry (28.4%). In Native American ancestry, carriers of the TT genotype of this variant averaged 4-fold higher bilirubinemia compared to the CC genotype (p = 2.82 × 10-12). We showed for the first time that UGT1A1 variants are the primary determinant of bilirubin levels in Native American ancestry, confirming its pan-ethnic relevance. Our study illustrates the general value of the local ancestry deconvolution approach to assessing isolated ancestry effects in admixed populations.

New insights from GWAS on BMI-related growth traits in a longitudinal cohort of admixed children with Native American and European ancestry

Body-mass index (BMI) is a hallmark of adiposity. In contrast with adulthood, the genetic architecture of BMI during childhood is poorly understood. The few genome-wide association studies (GWAS) on children have been performed almost exclusively in Europeans and at single ages. We performed cross-sectional and longitudinal GWAS for BMI-related traits on 904 admixed children with mostly Mapuche Native American and European ancestries. We found regulatory variants of the immune gene HLA-DQB3 strongly associated with BMI at 1.5 - 2.5 years old. A variant in the sex-determining gene DMRT1 was associated with the age at adiposity rebound (Age-AR) in girls (P = 9.8 × 10 - 9 ). BMI was significantly higher in Mapuche than in Europeans between 5.5 and 16.5 years old. Finally, Age-AR was significantly lower (P = 0.004 ) by 1.94 years and BMI at AR was significantly higher (P = 0.04 ) by 1.2 kg/m 2 , in Mapuche children compared with Europeans.

Genome-Wide Association Study and Polygenic Risk Scores of Serum DHEAS Levels in a Chilean Children Cohort

CONTEXT

Adrenarche reflects the developmental growth of the adrenal zona reticularis, which produces increasing adrenal androgen secretion (eg, dehydroepiandrosterone [DHEA]/dehydroepiandrosterone sulfate [DHEAS]) from approximately age 5 to 15 years.

OBJECTIVE

We hypothesized that the study of the genetic determinants associated with variations in serum DHEAS during adrenarche might detect genetic variants influencing the rate or timing of this process.

METHODS

Genome-wide genotyping was performed in participants of the Chilean pediatric Growth and Obesity Chilean Cohort Study (GOCS) cohort (n = 788). We evaluated the genetic determinants of DHEAS levels at the genome-wide level and in targeted genes associated with steroidogenesis. To corroborate our findings, we evaluated a polygenic risk score (PRS) for age at pubarche, based on the discovered variants, in children from the same cohort.

RESULTS

We identified one significant variant at the genome-wide level in the full cohort, close to the GALR1 gene (P = 3.81 × 10-8). In addition, variants suggestive of association (P < 1 × 10-5) were observed in PRLR, PITX1, PTPRD, NR1H4, and BCL11B. Stratifying by sex, we found variants suggestive of association in SERBP1 and CAMTA1/VAMP3 for boys and near ZNF98, TRPC6, and SULT2A1 for girls. We also found significant reductions in age at pubarche in those children with higher PRS for greater DHEAS based on these newly identified variants.

CONCLUSION

Our results disclose one variant associated with DHEAS concentrations at the level of genome-wide association study significance, and several variants with a suggestive association that may be involved in the genetic regulation of adrenarche.