Natural selection of immune and metabolic genes associated with health in two lowland Bolivian populations.

A growing body of work has addressed human adaptations to diverse environments using genomic data, but few studies have connected putatively selected alleles to phenotypes, much less among underrepresented populations such as Amerindians. Studies of natural selection and genotype-phenotype relationships in underrepresented populations hold potential to uncover previously undescribed loci underlying evolutionarily and biomedically relevant traits. Here, we worked with the Tsimane and the Moseten, two Amerindian populations inhabiting the Bolivian lowlands. We focused most intensively on the Tsimane, because long-term anthropological work with this group has shown that they have a high burden of both macro and microparasites, as well as minimal cardiometabolic disease or dementia. We therefore generated genome-wide genotype data for Tsimane individuals to study natural selection, and paired this with blood mRNA-seq as well as cardiometabolic and immune biomarker data generated from a larger sample that included both populations. In the Tsimane, we identified 21 regions that are candidates for selective sweeps, as well as 5 immune traits that show evidence for polygenic selection (e.g., C-reactive protein levels and the response to coronaviruses). Genes overlapping candidate regions were strongly enriched for known involvement in immune-related traits, such as abundance of lymphocytes and eosinophils. Importantly, we were also able to draw on extensive phenotype information for the Tsimane and Moseten and link five regions (containing PSD4, MUC21 and MUC22, TOX2, ANXA6, and ABCA1) with biomarkers of immune and metabolic function. Together, our work highlights the utility of pairing evolutionary analyses with anthropological and biomedical data to gain insight into the genetic basis of health-related traits.

Infant Growth Trajectories and Lipid Levels in Adolescence: Evidence From a Chilean Infancy Cohort.

Growth in early infancy is hypothesized to affect chronic disease risk factors later in life. To date, most reports draw on European-ancestry cohorts with few repeated observations in early infancy. We investigated the association between infant growth before 6 months and lipid levels in adolescents in a Hispanic/Latino cohort. We characterized infant growth from birth to 5 months in male (n = 311) and female (n = 285) infants from the Santiago Longitudinal Study (1991-1996) using 3 metrics: weight (kg), length (cm), and weight-for-length (g/cm). Superimposition by translation and rotation (SITAR) and latent growth mixture models (LGMMs) were used to estimate the association between infant growth characteristics and lipid levels at age 17 years. We found a positive relationship between the SITAR length velocity parameter before 6 months of age and high-density lipoprotein cholesterol levels in adolescence (11.5, 95% confidence interval; 3.4, 19.5), indicating higher high-density lipoprotein cholesterol levels occurring with faster length growth. The strongest associations from the LGMMs were between higher low-density lipoprotein cholesterol and slower weight-for-length growth, following a pattern of associations between slower growth and adverse lipid profiles. Further research in this window of time can confirm the association between early infant growth as an exposure and adolescent cardiovascular disease risk factors.

Genome-wide association study identifying novel variant for fasting insulin and allelic heterogeneity in known glycemic loci in Chilean adolescents: The Santiago Longitudinal Study.

BACKGROUND: The genetic underpinnings of glycemic traits have been understudied in adolescent and Hispanic/Latino (H/L) populations in comparison to adults and populations of European ancestry. OBJECTIVE: To identify genetic factors underlying glycemic traits in an adolescent H/L population. METHODS: We conducted a genome-wide association study (GWAS) of fasting glucose (FG) and fasting insulin (FI) in H/L adolescents from the Santiago Longitudinal Study. RESULTS: We identified one novel variant positioned in the CSMD1 gene on chromosome 8 (rs77465890, effect allele frequency = 0.10) that was associated with FI (ß = -0.299, SE = 0.054, p = 2.72×10-8 ) and was only slightly attenuated after adjusting for body mass index z-scores (ß = -0.252, SE = 0.047, p = 1.03×10-7 ). We demonstrated directionally consistent, but not statistically significant results in African and Hispanic adults of the Population Architecture Using Genomics and Epidemiology Consortium. We also identified secondary signals for two FG loci after conditioning on known variants, which demonstrate allelic heterogeneity in well-known glucose loci. CONCLUSION: Our results exemplify the importance of including populations with diverse ancestral origin and adolescent participants in GWAS of glycemic traits to uncover novel risk loci and expand our understanding of disease aetiology.

Sociodemographic predictors of early postnatal growth: evidence from a Chilean infancy cohort.

OBJECTIVES: Infant anthropometric growth varies across socioeconomic factors, including maternal education and income, and may serve as an indicator of environmental influences in early life with long-term health consequences. Previous research has identified sociodemographic gradients in growth with a focus on the first year and beyond, but estimates are sparse for growth before 6 months. Thus, our objective was to examine the relationship between sociodemographic factors and infant growth patterns between birth and 5 months of age. DESIGN: Prospective cohort study. SETTINGS: Low-income to middle-income neighbourhoods in Santiago, Chile (1991-1996). PARTICIPANTS: 1412 participants from a randomised iron-deficiency anaemia preventive trial in healthy infants. MAIN OUTCOME MEASURES: Longitudinal anthropometrics including monthly weight (kg), length (cm) and weight-for-length (WFL) values. For each measure, we estimated three individual-level growth parameters (size, timing and velocity) from SuperImposition by Translation and Rotation models. Size and timing changes represent vertical and horizontal growth curve shifts, respectively, and velocity change represents growth rate shifts. We estimated the linear association between growth parameters and gestational age, maternal age, education and socioeconomic position (SEP). RESULTS: Lower SEP was associated with a slower linear (length) velocity growth parameter (-0.22, 95% CI -0.31 to -0.13)-outcome units are per cent change in velocity from the average growth curve. Lower SEP was associated with later WFL growth timing as demonstrated through the tempo growth parameter for females (0.25, 95% CI 0.05 to 0.42)-outcome units are shifts in days from the average growth curve. We found no evidence of associations between SEP and the weight size, timing or velocity growth rate parameters. CONCLUSION: Previous research on growth in older infants and children shows associations between lower SEP with slower length velocity. We found evidence supporting this association in the first 5 months of life, which may inform age-specific prevention efforts aimed at infant length growth.

A genome-wide association study of polycystic ovary syndrome identified from electronic health records.

BACKGROUND: Polycystic ovary syndrome is the most common endocrine disorder affecting women of reproductive age. A number of criteria have been developed for clinical diagnosis of polycystic ovary syndrome, with the Rotterdam criteria being the most inclusive. Evidence suggests that polycystic ovary syndrome is significantly heritable, and previous studies have identified genetic variants associated with polycystic ovary syndrome diagnosed using different criteria. The widely adopted electronic health record system provides an opportunity to identify patients with polycystic ovary syndrome using the Rotterdam criteria for genetic studies. OBJECTIVE: To identify novel associated genetic variants under the same phenotype definition, we extracted polycystic ovary syndrome cases and unaffected controls based on the Rotterdam criteria from the electronic health records and performed a discovery-validation genome-wide association study. STUDY DESIGN: We developed a polycystic ovary syndrome phenotyping algorithm on the basis of the Rotterdam criteria and applied it to 3 electronic health record-linked biobanks to identify cases and controls for genetic study. In the discovery phase, we performed an individual genome-wide association study using the Geisinger MyCode and the Electronic Medical Records and Genomics cohorts, which were then meta-analyzed. We attempted validation of the significant association loci (P<1×10-6) in the BioVU cohort. All association analyses used logistic regression, assuming an additive genetic model, and adjusted for principal components to control for population stratification. An inverse-variance fixed-effect model was adopted for meta-analysis. In addition, we examined the top variants to evaluate their associations with each criterion in the phenotyping algorithm. We used the STRING database to characterize protein-protein interaction network. RESULTS: Using the same algorithm based on the Rotterdam criteria, we identified 2995 patients with polycystic ovary syndrome and 53,599 population controls in total (2742 cases and 51,438 controls from the discovery phase; 253 cases and 2161 controls in the validation phase). We identified 1 novel genome-wide significant variant rs17186366 (odds ratio [OR]=1.37 [1.23, 1.54], P=2.8×10-8) located near SOD2. In addition, 2 loci with suggestive association were also identified: rs113168128 (OR=1.72 [1.42, 2.10], P=5.2×10-8), an intronic variant of ERBB4 that is independent from the previously published variants, and rs144248326 (OR=2.13 [1.52, 2.86], P=8.45×10-7), a novel intronic variant in WWTR1. In the further association tests of the top 3 single-nucleotide polymorphisms with each criterion in the polycystic ovary syndrome algorithm, we found that rs17186366 (SOD2) was associated with polycystic ovaries and hyperandrogenism, whereas rs11316812 (ERBB4) and rs144248326 (WWTR1) were mainly associated with oligomenorrhea or infertility. We also validated the previously reported association with DENND1A1. Using the STRING database to characterize protein-protein interactions, we found both ERBB4 and WWTR1 can interact with YAP1, which has been previously associated with polycystic ovary syndrome. CONCLUSION: Through a discovery-validation genome-wide association study on polycystic ovary syndrome identified from electronic health records using an algorithm based on Rotterdam criteria, we identified and validated a novel genome-wide significant association with a variant near SOD2. We also identified a novel independent variant within ERBB4 and a suggestive association with WWTR1. With previously identified polycystic ovary syndrome gene YAP1, the ERBB4-YAP1-WWTR1 network suggests involvement of the epidermal growth factor receptor and the Hippo pathway in the multifactorial etiology of polycystic ovary syndrome.