Genome-Wide Analyses of Vocabulary Size in Infancy and Toddlerhood: Associations With Attention-Deficit/Hyperactivity Disorder, Literacy, and Cognition-Related Traits

BACKGROUND

The number of words children produce (expressive vocabulary) and understand (receptive vocabulary) changes rapidly during early development, partially due to genetic factors. Here, we performed a meta-genome-wide association study of vocabulary acquisition and investigated polygenic overlap with literacy, cognition, developmental phenotypes, and neurodevelopmental conditions, including attention-deficit/hyperactivity disorder (ADHD).

METHODS

We studied 37,913 parent-reported vocabulary size measures (English, Dutch, Danish) for 17,298 children of European descent. Meta-analyses were performed for early-phase expressive (infancy, 15-18 months), late-phase expressive (toddlerhood, 24-38 months), and late-phase receptive (toddlerhood, 24-38 months) vocabulary. Subsequently, we estimated single nucleotide polymorphism-based heritability (SNP-h2) and genetic correlations (rg) and modeled underlying factor structures with multivariate models.

RESULTS

Early-life vocabulary size was modestly heritable (SNP-h2 = 0.08-0.24). Genetic overlap between infant expressive and toddler receptive vocabulary was negligible (rg = 0.07), although each measure was moderately related to toddler expressive vocabulary (rg = 0.69 and rg = 0.67, respectively), suggesting a multifactorial genetic architecture. Both infant and toddler expressive vocabulary were genetically linked to literacy (e.g., spelling: rg = 0.58 and rg = 0.79, respectively), underlining genetic similarity. However, a genetic association of early-life vocabulary with educational attainment and intelligence emerged only during toddlerhood (e.g., receptive vocabulary and intelligence: rg = 0.36). Increased ADHD risk was genetically associated with larger infant expressive vocabulary (rg = 0.23). Multivariate genetic models in the ALSPAC (Avon Longitudinal Study of Parents and Children) cohort confirmed this finding for ADHD symptoms (e.g., at age 13; rg = 0.54) but showed that the association effect reversed for toddler receptive vocabulary (rg = -0.74), highlighting developmental heterogeneity.

CONCLUSIONS

The genetic architecture of early-life vocabulary changes during development, shaping polygenic association patterns with later-life ADHD, literacy, and cognition-related traits.

Genome-Wide Interaction Analysis With DASH Diet Score Identified Novel Loci for Systolic Blood Pressure

BACKGROUND

The Dietary Approaches to Stop Hypertension (DASH) diet score lowers blood pressure (BP). We examined interactions between genotype and the DASH diet score in relation to systolic BP.

METHODS

We analyzed up to 9 420 585 single nucleotide polymorphisms in up to 127 282 individuals of 6 population groups (91% of European population) from the Cohorts for Heart and Aging Research in Genomic Epidemiology consortium (n=35 660) and UK Biobank (n=91 622) and performed European population-specific and cross-population meta-analyses.

RESULTS

We identified 3 loci in European-specific analyses and an additional 4 loci in cross-population analyses at Pinteraction<5e-8. We observed a consistent interaction between rs117878928 at 15q25.1 (minor allele frequency, 0.03) and the DASH diet score (Pinteraction=4e-8; P for heterogeneity, 0.35) in European population, where the interaction effect size was 0.42±0.09 mm Hg (Pinteraction=9.4e-7) and 0.20±0.06 mm Hg (Pinteraction=0.001) in Cohorts for Heart and Aging Research in Genomic Epidemiology and the UK Biobank, respectively. The 1 Mb region surrounding rs117878928 was enriched with cis-expression quantitative trait loci (eQTL) variants (P=4e-273) and cis-DNA methylation quantitative trait loci variants (P=1e-300). Although the closest gene for rs117878928 is MTHFS, the highest narrow sense heritability accounted by single nucleotide polymorphisms potentially interacting with the DASH diet score in this locus was for gene ST20 at 15q25.1.

CONCLUSIONS

We demonstrated gene-DASH diet score interaction effects on systolic BP in several loci. Studies with larger diverse populations are needed to validate our findings.

Trans-ancestral genome-wide association study of longitudinal pubertal height growth and shared heritability with adult health outcomes

BACKGROUND

Pubertal growth patterns correlate with future health outcomes. However, the genetic mechanisms mediating growth trajectories remain largely unknown. Here, we modeled longitudinal height growth with Super-Imposition by Translation And Rotation (SITAR) growth curve analysis on ~ 56,000 trans-ancestry samples with repeated height measurements from age 5 years to adulthood. We performed genetic analysis on six phenotypes representing the magnitude, timing, and intensity of the pubertal growth spurt. To investigate the lifelong impact of genetic variants associated with pubertal growth trajectories, we performed genetic correlation analyses and phenome-wide association studies in the Penn Medicine BioBank and the UK Biobank.

RESULTS

Large-scale growth modeling enables an unprecedented view of adolescent growth across contemporary and 20th-century pediatric cohorts. We identify 26 genome-wide significant loci and leverage trans-ancestry data to perform fine-mapping. Our data reveals genetic relationships between pediatric height growth and health across the life course, with different growth trajectories correlated with different outcomes. For instance, a faster tempo of pubertal growth correlates with higher bone mineral density, HOMA-IR, fasting insulin, type 2 diabetes, and lung cancer, whereas being taller at early puberty, taller across puberty, and having quicker pubertal growth were associated with higher risk for atrial fibrillation.

CONCLUSION

We report novel genetic associations with the tempo of pubertal growth and find that genetic determinants of growth are correlated with reproductive, glycemic, respiratory, and cardiac traits in adulthood. These results aid in identifying specific growth trajectories impacting lifelong health and show that there may not be a single "optimal" pubertal growth pattern.

Dietary and supplemental intake of vitamins C and E is associated with altered DNA methylation in an epigenome-wide association study meta-analysis

BACKGROUND

Dietary intake of antioxidants such as vitamins C and E protect against oxidative stress, and may also be associated with altered DNA methylation patterns.

METHODS

We meta-analysed epigenome-wide association study (EWAS) results from 11,866 participants across eight population-based cohorts to evaluate the association between self-reported dietary and supplemental intake of vitamins C and E with DNA methylation. EWAS were adjusted for age, sex, BMI, caloric intake, blood cell type proportion, smoking status, alcohol consumption, and technical covariates. Significant results of the meta-analysis were subsequently evaluated in gene set enrichment analysis (GSEA) and expression quantitative trait methylation (eQTM) analysis.

RESULTS

In meta-analysis, methylation at 4,656 CpG sites was significantly associated with vitamin C intake at FDR ≤ 0.05. The most significant CpG sites associated with vitamin C (at FDR ≤ 0.01) were enriched for pathways associated with systems development and cell signalling in GSEA, and were associated with downstream expression of genes enriched in the immune response in eQTM analysis. Furthermore, methylation at 160 CpG sites was significantly associated with vitamin E intake at FDR ≤ 0.05, but GSEA and eQTM analysis of the top most significant CpG sites associated with vitamin E did not identify significant enrichment of any biological pathways investigated.

CONCLUSIONS

We identified significant associations of many CpG sites with vitamin C and E intake, and our results suggest that vitamin C intake may be associated with systems development and the immune response.

The Role of Genetics in Preterm Birth

Preterm birth (PTB), defined as the birth of a child before 37 completed weeks gestation, affects approximately 11% of live births and is the leading cause of death in children under 5 years. PTB is a complex disease with multiple risk factors including genetic variation. Much research has aimed to establish the biological mechanisms underlying PTB often through identification of genetic markers for PTB risk. The objective of this review is to present a comprehensive and updated summary of the published data relating to the field of PTB genetics. A literature search in PubMed was conducted and English studies related to PTB genetics were included. Genetic studies have identified genes within inflammatory, immunological, tissue remodeling, endocrine, metabolic, and vascular pathways that may be involved in PTB. However, a substantial proportion of published data have been largely inconclusive and multiple studies had limited power to detect associations. On the contrary, a few large hypothesis-free approaches have identified and replicated multiple novel variants associated with PTB in different cohorts. Overall, attempts to predict PTB using single "-omics" datasets including genomic, transcriptomic, and epigenomic biomarkers have been mostly unsuccessful and have failed to translate to the clinical setting. Integration of data from multiple "-omics" datasets has yielded the most promising results.

Genome-Wide Interaction Analysis with DASH Diet Score Identified Novel Loci for Systolic Blood Pressure

Objective

We examined interactions between genotype and a Dietary Approaches to Stop Hypertension (DASH) diet score in relation to systolic blood pressure (SBP).

Methods

We analyzed up to 9,420,585 biallelic imputed single nucleotide polymorphisms (SNPs) in up to 127,282 individuals of six population groups (91% of European population) from the Cohorts for Heart and Aging Research in Genomic Epidemiology consortium (CHARGE; n=35,660) and UK Biobank (n=91,622) and performed European population-specific and cross-population meta-analyses.

Results

We identified three loci in European-specific analyses and an additional four loci in cross-population analyses at P for interaction < 5e-8. We observed a consistent interaction between rs117878928 at 15q25.1 (minor allele frequency = 0.03) and the DASH diet score (P for interaction = 4e-8; P for heterogeneity = 0.35) in European population, where the interaction effect size was 0.42±0.09 mm Hg (P for interaction = 9.4e-7) and 0.20±0.06 mm Hg (P for interaction = 0.001) in CHARGE and the UK Biobank, respectively. The 1 Mb region surrounding rs117878928 was enriched with cis-expression quantitative trait loci (eQTL) variants (P = 4e-273) and cis-DNA methylation quantitative trait loci (mQTL) variants (P = 1e-300). While the closest gene for rs117878928 is MTHFS, the highest narrow sense heritability accounted by SNPs potentially interacting with the DASH diet score in this locus was for gene ST20 at 15q25.1.

Conclusion

We demonstrated gene-DASH diet score interaction effects on SBP in several loci. Studies with larger diverse populations are needed to validate our findings.

Genome-wide association study of placental weight identifies distinct and shared genetic influences between placental and fetal growth

A well-functioning placenta is essential for fetal and maternal health throughout pregnancy. Using placental weight as a proxy for placental growth, we report genome-wide association analyses in the fetal (n = 65,405), maternal (n = 61,228) and paternal (n = 52,392) genomes, yielding 40 independent association signals. Twenty-six signals are classified as fetal, four maternal and three fetal and maternal. A maternal parent-of-origin effect is seen near KCNQ1. Genetic correlation and colocalization analyses reveal overlap with birth weight genetics, but 12 loci are classified as predominantly or only affecting placental weight, with connections to placental development and morphology, and transport of antibodies and amino acids. Mendelian randomization analyses indicate that fetal genetically mediated higher placental weight is causally associated with preeclampsia risk and shorter gestational duration. Moreover, these analyses support the role of fetal insulin in regulating placental weight, providing a key link between fetal and placental growth.

Changes in Refractive Error During Young Adulthood: The Effects of Longitudinal Screen Time, Ocular Sun Exposure, and Genetic Predisposition

Purpose

Changes in refractive error during young adulthood is common yet risk factors at this age are largely unexplored. This study explored risk factors for these changes, including gene-environmental interactions.

Methods

Spherical equivalent refraction (SER) and axial length (AL) for 624 community-based adults were measured at 20 (baseline) and 28 years old. Participants were genotyped and their polygenic scores (PGS) for refractive error calculated. Self-reported screen time (computer, television, and mobile devices) from 20 to 28 years old were collected prospectively and longitudinal trajectories were generated. Past sun exposure was quantified using conjunctival ultraviolet autofluorescence (CUVAF) area.

Results

Median change in SER and AL were -0.023 diopters (D)/year (interquartile range [IQR] = -0.062 to -0.008) and +0.01 mm/year (IQR = 0.000 to 0.026), respectively. Sex, baseline myopia, parental myopia, screen time, CUVAF, and PGS were significantly associated with myopic shift. Collectively, these factors accounted for approximately 20% of the variance in refractive error change, with screen time, CUVAF, and PGS each explaining approximately 1% of the variance. Four trajectories for total screen time were found: "consistently low" (n = 148), "consistently high" (n = 250), "consistently very high" (n = 76), and "increasing" (n = 150). Myopic shift was faster in those with "consistently high" or "consistently very high" screen time compared to "consistently-low" (P ≤ 0.031). For each z-score increase in PGS, changes in SER and AL increased by -0.005 D/year and 0.002 mm/year (P ≤ 0.045). Of the three types of screen time, only computer time was associated with myopic shift (P ≤ 0.040). There was no two- or three-way interaction effect between PGS, CUVAF, or screen time (P ≥ 0.26).

Conclusions

Higher total or computer screen time, less sun exposure, and genetic predisposition are each independently associated with greater myopic shifts during young adulthood. Given that these factors explained only a small amount of the variance, there are likely other factors driving refractive error change during young adulthood.

European and multi-ancestry genome-wide association meta-analysis of atopic dermatitis highlights importance of systemic immune regulation

Atopic dermatitis (AD) is a common inflammatory skin condition and prior genome-wide association studies (GWAS) have identified 71 associated loci. In the current study we conducted the largest AD GWAS to date (discovery N = 1,086,394, replication N = 3,604,027), combining previously reported cohorts with additional available data. We identified 81 loci (29 novel) in the European-only analysis (which all replicated in a separate European analysis) and 10 additional loci in the multi-ancestry analysis (3 novel). Eight variants from the multi-ancestry analysis replicated in at least one of the populations tested (European, Latino or African), while two may be specific to individuals of Japanese ancestry. AD loci showed enrichment for DNAse I hypersensitivity and eQTL associations in blood. At each locus we prioritised candidate genes by integrating multi-omic data. The implicated genes are predominantly in immune pathways of relevance to atopic inflammation and some offer drug repurposing opportunities.

Combined vaginal progesterone and cervical cerclage in the prevention of preterm birth: a systematic review and meta-analysis

OBJECTIVE

Vaginal progesterone and cervical cerclage are both effective interventions for reducing preterm birth. It is currently unclear whether combined therapy offers superior effectiveness than single therapy. This study aimed to determine the efficacy of combining cervical cerclage and vaginal progesterone in the prevention of preterm birth.

DATA SOURCES

We searched Medline (Ovid), EMBASE (Ovid), PsycINFO (Ovid), CINAHL (EBSCOhost), Cochrane Library (Wiley), and Scopus (from their inception to 2020).

STUDY ELIGIBILITY CRITERIA

The review accepted randomized and pseudorandomized control trials, nonrandomized experimental control trials, and cohort studies. High risk patients (shortened cervical length <25mm or previous preterm birth) who were assigned cervical cerclage, vaginal progesterone, or both for the prevention of preterm birth were included. Only singleton pregnancies were assessed.

METHODS

The primary outcome was birth <37 weeks. Secondary outcomes included birth <28 weeks, <32 weeks and <34 weeks, gestational age at delivery, days between intervention and delivery, preterm premature rupture of membranes, cesarean delivery, neonatal mortality, neonatal intensive care unit admission, intubation, and birthweight. Following title and full-text screening, 11 studies were included in the final analysis. Risk of bias was assessed using the Cochrane Collaboration tool for assessing the risk of bias (ROBINS-I and RoB-2). Quality of evidence was assessed using the GRADE (Grading of Recommendations, Assessment, Development, and Evaluations) tool.

RESULTS

Combined therapy was associated with lower risk of preterm birth at <37 weeks than cerclage alone (risk ratio, 0.51; 95% confidence interval, 0.37-0.79) or progesterone alone (risk ratio, 0.75; 95% confidence interval, 0.58-0.96). Compared with cerclage only, combined therapy was associated with preterm birth at <34 weeks, <32 weeks, or <28 weeks, decreased neonatal mortality, increased birthweight, increased gestational age, and a longer interval between intervention and delivery. Compared with progesterone alone, combined therapy was associated with preterm birth at <32 weeks, <28 weeks, decreased neonatal mortality, increased birthweight, and increased gestational age. There were no differences in any other secondary outcomes.

CONCLUSION

Combined treatment of cervical cerclage and vaginal progesterone could potentially result in a greater reduction in preterm birth than in single therapy. Further, well-conducted and adequately powered randomized controlled trials are needed to assess these promising findings.

Language and reading impairments are associated with increased prevalence of non-right-handedness

Handedness has been studied for association with language-related disorders because of its link with language hemispheric dominance. No clear pattern has emerged, possibly because of small samples, publication bias, and heterogeneous criteria across studies. Non-right-handedness (NRH) frequency was assessed in N = 2503 cases with reading and/or language impairment and N = 4316 sex-matched controls identified from 10 distinct cohorts (age range 6-19 years old; European ethnicity) using a priori set criteria. A meta-analysis (Ncases  = 1994) showed elevated NRH % in individuals with language/reading impairment compared with controls (OR = 1.21, CI = 1.06-1.39, p = .01). The association between reading/language impairments and NRH could result from shared pathways underlying brain lateralization, handedness, and cognitive functions.

Understanding the genetic complexity of puberty timing across the allele frequency spectrum

Pubertal timing varies considerably and has been associated with a range of health outcomes in later life. To elucidate the underlying biological mechanisms, we performed multi-ancestry genetic analyses in ~800,000 women, identifying 1,080 independent signals associated with age at menarche. Collectively these loci explained 11% of the trait variance in an independent sample, with women at the top and bottom 1% of polygenic risk exhibiting a ~11 and ~14-fold higher risk of delayed and precocious pubertal development, respectively. These common variant analyses were supported by exome sequence analysis of ~220,000 women, identifying several genes, including rare loss of function variants in ZNF483 which abolished the impact of polygenic risk. Next, we implicated 660 genes in pubertal development using a combination of in silico variant-to-gene mapping approaches and integration with dynamic gene expression data from mouse embryonic GnRH neurons. This included an uncharacterized G-protein coupled receptor GPR83, which we demonstrate amplifies signaling of MC3R, a key sensor of nutritional status. Finally, we identified several genes, including ovary-expressed genes involved in DNA damage response that co-localize with signals associated with menopause timing, leading us to hypothesize that the ovarian reserve might signal centrally to trigger puberty. Collectively these findings extend our understanding of the biological complexity of puberty timing and highlight body size dependent and independent mechanisms that potentially link reproductive timing to later life disease.

Development and internal validation of a non-invasive clinical tool to predict sufficient omega-3 levels in early pregnancy

BACKGROUND

Complications from preterm birth (PTB) are the leading cause of death and disability in those under five years. Whilst the role of omega-3 (n-3) supplementation in reducing PTB is well-established, growing evidence suggests supplementation use in those replete may increase the risk of early PTB.

AIM

To develop a non-invasive tool to identify individuals with total n-3 serum levels above 4.3% of total fatty acids in early pregnancy.

METHODS

We conducted a prospective observational study recruiting 331 participants from three clinical sites in Newcastle, Australia. Eligible participants (n = 307) had a singleton pregnancy between 8 and 20 weeks' gestation at recruitment. Data on factors associated with n-3 serum levels were collected using an electronic questionnaire; these included estimated intake of n-3 (including food type, portion size, frequency of consumption), n-3 supplementation, and sociodemographic factors. The optimal cut-point of estimated n-3 intake that predicted mothers with total serum n-3 levels likely above 4.3% was developed using multivariate logistic regression, adjusting for maternal age, body mass index, socioeconomic status, and n-3 supplementation use. Total serum n-3 levels above 4.3% was selected as previous research has demonstrated that mothers with these levels are at increased risk of early PTB if they take additional n-3 supplementation during pregnancy. Models were evaluated using various performance metrics including sensitivity, specificity, area under receiver operator characteristic (AUROC) curve, true positive rate (TPR) at 10% false positive rate (FPR), Youden Index, Closest to (0,1) Criteria, Concordance Probability, and Index of Union. Internal validation was performed using 1000-bootstraps to generate 95% confidence intervals for performance metrics generated.

RESULTS

Of 307 eligible participants included for analysis, 58.6% had total n-3 serum levels above 4.3%. The optimal model had a moderate discriminative ability (AUROC 0.744, 95% CI 0.742-0.746) with 84.7% sensitivity, 54.7% specificity and 37.6% TPR at 10% FPR.

CONCLUSIONS

Our non-invasive tool was a moderate predictor of pregnant women with total serum n-3 levels above 4.3%; however, its performance is not yet adequate for clinical use.

TRIAL REGISTRATION

This trial was approved by the Hunter New England Human Research Ethics Committee of the Hunter New England Local Health District (Reference 2020/ETH00498 on 07/05/2020 and 2020/ETH02881 on 08/12/2020).

Genetic effects on the timing of parturition and links to fetal birth weight

The timing of parturition is crucial for neonatal survival and infant health. Yet, its genetic basis remains largely unresolved. We present a maternal genome-wide meta-analysis of gestational duration (n = 195,555), identifying 22 associated loci (24 independent variants) and an enrichment in genes differentially expressed during labor. A meta-analysis of preterm delivery (18,797 cases, 260,246 controls) revealed six associated loci and large genetic similarities with gestational duration. Analysis of the parental transmitted and nontransmitted alleles (n = 136,833) shows that 15 of the gestational duration genetic variants act through the maternal genome, whereas 7 act both through the maternal and fetal genomes and 2 act only via the fetal genome. Finally, the maternal effects on gestational duration show signs of antagonistic pleiotropy with the fetal effects on birth weight: maternal alleles that increase gestational duration have negative fetal effects on birth weight. The present study provides insights into the genetic effects on the timing of parturition and the complex maternal-fetal relationship between gestational duration and birth weight.

Multi-ancestry genome-wide association analyses improve resolution of genes and pathways influencing lung function and chronic obstructive pulmonary disease risk

Lung-function impairment underlies chronic obstructive pulmonary disease (COPD) and predicts mortality. In the largest multi-ancestry genome-wide association meta-analysis of lung function to date, comprising 580,869 participants, we identified 1,020 independent association signals implicating 559 genes supported by ≥2 criteria from a systematic variant-to-gene mapping framework. These genes were enriched in 29 pathways. Individual variants showed heterogeneity across ancestries, age and smoking groups, and collectively as a genetic risk score showed strong association with COPD across ancestry groups. We undertook phenome-wide association studies for selected associated variants as well as trait and pathway-specific genetic risk scores to infer possible consequences of intervening in pathways underlying lung function. We highlight new putative causal variants, genes, proteins and pathways, including those targeted by existing drugs. These findings bring us closer to understanding the mechanisms underlying lung function and COPD, and should inform functional genomics experiments and potentially future COPD therapies.

Implicating genes, pleiotropy, and sexual dimorphism at blood lipid loci through multi-ancestry meta-analysis

BACKGROUND

Genetic variants within nearly 1000 loci are known to contribute to modulation of blood lipid levels. However, the biological pathways underlying these associations are frequently unknown, limiting understanding of these findings and hindering downstream translational efforts such as drug target discovery.

RESULTS

To expand our understanding of the underlying biological pathways and mechanisms controlling blood lipid levels, we leverage a large multi-ancestry meta-analysis (N = 1,654,960) of blood lipids to prioritize putative causal genes for 2286 lipid associations using six gene prediction approaches. Using phenome-wide association (PheWAS) scans, we identify relationships of genetically predicted lipid levels to other diseases and conditions. We confirm known pleiotropic associations with cardiovascular phenotypes and determine novel associations, notably with cholelithiasis risk. We perform sex-stratified GWAS meta-analysis of lipid levels and show that 3-5% of autosomal lipid-associated loci demonstrate sex-biased effects. Finally, we report 21 novel lipid loci identified on the X chromosome. Many of the sex-biased autosomal and X chromosome lipid loci show pleiotropic associations with sex hormones, emphasizing the role of hormone regulation in lipid metabolism.

CONCLUSIONS

Taken together, our findings provide insights into the biological mechanisms through which associated variants lead to altered lipid levels and potentially cardiovascular disease risk.

A saturated map of common genetic variants associated with human height

Common single-nucleotide polymorphisms (SNPs) are predicted to collectively explain 40-50% of phenotypic variation in human height, but identifying the specific variants and associated regions requires huge sample sizes1. Here, using data from a genome-wide association study of 5.4 million individuals of diverse ancestries, we show that 12,111 independent SNPs that are significantly associated with height account for nearly all of the common SNP-based heritability. These SNPs are clustered within 7,209 non-overlapping genomic segments with a mean size of around 90 kb, covering about 21% of the genome. The density of independent associations varies across the genome and the regions of increased density are enriched for biologically relevant genes. In out-of-sample estimation and prediction, the 12,111 SNPs (or all SNPs in the HapMap 3 panel2) account for 40% (45%) of phenotypic variance in populations of European ancestry but only around 10-20% (14-24%) in populations of other ancestries. Effect sizes, associated regions and gene prioritization are similar across ancestries, indicating that reduced prediction accuracy is likely to be explained by linkage disequilibrium and differences in allele frequency within associated regions. Finally, we show that the relevant biological pathways are detectable with smaller sample sizes than are needed to implicate causal genes and variants. Overall, this study provides a comprehensive map of specific genomic regions that contain the vast majority of common height-associated variants. Although this map is saturated for populations of European ancestry, further research is needed to achieve equivalent saturation in other ancestries.

Genome-wide analyses of individual differences in quantitatively assessed reading- and language-related skills in up to 34,000 people

The use of spoken and written language is a fundamental human capacity. Individual differences in reading- and language-related skills are influenced by genetic variation, with twin-based heritability estimates of 30 to 80% depending on the trait. The genetic architecture is complex, heterogeneous, and multifactorial, but investigations of contributions of single-nucleotide polymorphisms (SNPs) were thus far underpowered. We present a multicohort genome-wide association study (GWAS) of five traits assessed individually using psychometric measures (word reading, nonword reading, spelling, phoneme awareness, and nonword repetition) in samples of 13,633 to 33,959 participants aged 5 to 26 y. We identified genome-wide significant association with word reading (rs11208009, P = 1.098 × 10-8) at a locus that has not been associated with intelligence or educational attainment. All five reading-/language-related traits showed robust SNP heritability, accounting for 13 to 26% of trait variability. Genomic structural equation modeling revealed a shared genetic factor explaining most of the variation in word/nonword reading, spelling, and phoneme awareness, which only partially overlapped with genetic variation contributing to nonword repetition, intelligence, and educational attainment. A multivariate GWAS of word/nonword reading, spelling, and phoneme awareness maximized power for follow-up investigation. Genetic correlation analysis with neuroimaging traits identified an association with the surface area of the banks of the left superior temporal sulcus, a brain region linked to the processing of spoken and written language. Heritability was enriched for genomic elements regulating gene expression in the fetal brain and in chromosomal regions that are depleted of Neanderthal variants. Together, these results provide avenues for deciphering the biological underpinnings of uniquely human traits.

A multi-layer functional genomic analysis to understand noncoding genetic variation in lipids

A major challenge of genome-wide association studies (GWASs) is to translate phenotypic associations into biological insights. Here, we integrate a large GWAS on blood lipids involving 1.6 million individuals from five ancestries with a wide array of functional genomic datasets to discover regulatory mechanisms underlying lipid associations. We first prioritize lipid-associated genes with expression quantitative trait locus (eQTL) colocalizations and then add chromatin interaction data to narrow the search for functional genes. Polygenic enrichment analysis across 697 annotations from a host of tissues and cell types confirms the central role of the liver in lipid levels and highlights the selective enrichment of adipose-specific chromatin marks in high-density lipoprotein cholesterol and triglycerides. Overlapping transcription factor (TF) binding sites with lipid-associated loci identifies TFs relevant in lipid biology. In addition, we present an integrative framework to prioritize causal variants at GWAS loci, producing a comprehensive list of candidate causal genes and variants with multiple layers of functional evidence. We highlight two of the prioritized genes, CREBRF and RRBP1, which show convergent evidence across functional datasets supporting their roles in lipid biology.

The power of genetic diversity in genome-wide association studies of lipids

Increased blood lipid levels are heritable risk factors of cardiovascular disease with varied prevalence worldwide owing to different dietary patterns and medication use1. Despite advances in prevention and treatment, in particular through reducing low-density lipoprotein cholesterol levels2, heart disease remains the leading cause of death worldwide3. Genome-wideassociation studies (GWAS) of blood lipid levels have led to important biological and clinical insights, as well as new drug targets, for cardiovascular disease. However, most previous GWAS4-23 have been conducted in European ancestry populations and may have missed genetic variants that contribute to lipid-level variation in other ancestry groups. These include differences in allele frequencies, effect sizes and linkage-disequilibrium patterns24. Here we conduct a multi-ancestry, genome-wide genetic discovery meta-analysis of lipid levels in approximately 1.65 million individuals, including 350,000 of non-European ancestries. We quantify the gain in studying non-European ancestries and provide evidence to support the expansion of recruitment of additional ancestries, even with relatively small sample sizes. We find that increasing diversity rather than studying additional individuals of European ancestry results in substantial improvements in fine-mapping functional variants and portability of polygenic prediction (evaluated in approximately 295,000 individuals from 7 ancestry groupings). Modest gains in the number of discovered loci and ancestry-specific variants were also achieved. As GWAS expand emphasis beyond the identification of genes and fundamental biology towards the use of genetic variants for preventive and precision medicine25, we anticipate that increased diversity of participants will lead to more accurate and equitable26 application of polygenic scores in clinical practice.

Rare variant analysis in eczema identifies exonic variants in DUSP1, NOTCH4 and SLC9A4

Previous genome-wide association studies revealed multiple common variants involved in eczema but the role of rare variants remains to be elucidated. Here, we investigate the role of rare variants in eczema susceptibility. We meta-analyze 21 study populations including 20,016 eczema cases and 380,433 controls. Rare variants are imputed with high accuracy using large population-based reference panels. We identify rare exonic variants in DUSP1, NOTCH4, and SLC9A4 to be associated with eczema. In DUSP1 and NOTCH4 missense variants are predicted to impact conserved functional domains. In addition, five novel common variants at SATB1-AS1/KCNH8, TRIB1/LINC00861, ZBTB1, TBX21/OSBPL7, and CSF2RB are discovered. While genes prioritized based on rare variants are significantly up-regulated in the skin, common variants point to immune cell function. Over 20% of the single nucleotide variant-based heritability is attributable to rare and low-frequency variants. The identified rare/low-frequency variants located in functional protein domains point to promising targets for novel therapeutic approaches to eczema.

Genetic studies of eczema to date have mostly explored common genetic variation. Here, the authors perform a large meta-analysis for common and rare variants and discover 8 loci associated with eczema. Over 20% of the heritability of the condition is attributable to rare variants.

Genome-wide analysis of thyroid function in Australian adolescents highlights SERPINA7 and NCOA3

OBJECTIVE

Genetic factors underpin the narrow intraindividual variability of thyroid function, although precise contributions of environmental vs genetic factors remain uncertain. We sought to clarify the heritability of thyroid function traits and thyroid peroxidase antibody (TPOAb) positivity and identify single nucleotide polymorphisms (SNPs) contributing to the trait variance.

METHODS

Heritability of thyroid-stimulating hormone (TSH), free T4 (fT4), free T3 (fT3) and TPOAb in a cohort of 2854 euthyroid, dizygous and monozygous twins (age range 11.9-16.9 years) from the Brisbane Longitudinal Twin Study (BLTS) was assessed using structural equation modelling. A genome-wide analysis was conducted on 2832 of these individuals across 7 522 526 SNPs as well as gene-based association analyses. Replication analysis of the association results was performed in the Raine Study (n = 1115) followed by meta-analysis to maximise power for discovery.

RESULTS

Heritability of thyroid function parameters in the BLTS was 70.8% (95% CI: 66.7-74.9%) for TSH, 67.5% (59.8-75.3%) for fT4, 59.7% (54.4-65.0%) for fT3 and 48.8% (40.6-56.9%) for TPOAb. The genome-wide association study (GWAS) in the discovery cohort identified a novel association between rs2026401 upstream of NCOA3 and TPOAb. GWAS meta-analysis found associations between TPOAb and rs445219, also near NCOA3, and fT3 and rs12687280 near SERPINA7. Gene-based association analysis highlighted SERPINA7 for fT3 and NPAS3 for fT4.

CONCLUSION

Our findings resolve former contention regarding heritability estimates of thyroid function traits and TPOAb positivity. GWAS and gene-based association analysis identified variants accounting for a component of this heritability.

A Life Course Approach to the Relationship Between Fetal Growth and Hypothalamic-Pituitary-Adrenal Axis Function

CONTEXT

Human and animal studies suggest that hypothalamic-pituitary-adrenal axis (HPA-A) function may be programmed in utero; however, these findings are inconsistent. Given the powerful metabolic actions of cortisol, it is important to clarify the influence of early life on adult HPA-A function.

OBJECTIVE

To determine the relationship between fetal growth and HPA-A stress response to a psychosocial stressor in young adults.

DESIGN

Multigenerational, prospective cohort study (the Raine Study) conducted between 1989 and 1991.

SETTING

King Edward Memorial Hospital, Perth, Western Australia, Australia.

PARTICIPANTS

A total of 917 participants aged 18 years from Gen2 of the Raine Study.

MAIN OUTCOME MEASURES

Measures of hypothalamic-pituitary-adrenal axis function before and after exposure to the Trier Social Stress Test.

RESULTS

In fully adjusted models, an inverse linear relationship was observed between birthweight and plasma measures of (1) baseline cortisol (β = -0.90%, 95% CI: -1.73 to -0.07; P = 0.03); (2) peak cortisol (β = -0.78%, 95% CI -1.51 to -0.06; P = 0.03); (3) area under the curve with respect to ground (β = -0.89%, 95% CI -1.60 to -0.18; P = 0.01); and (4) adrenal sensitivity (β = -1.02, 95% CI: -1.85 to -0.18; P = 0.02). Similar results were demonstrated for percent optimal birthweight. No consistent quadratic relationships were identified. No associations were found between measures of fetal adiposity and HPA-A function at age 18 years, or fetal growth and HPA-A response pattern. Removal of anticipatory responders from the models substantially attenuated the observed relationships.

CONCLUSION

We observed an inverse linear relationship between fetal growth and HPA-A function at age 18 years. This differs from the inverse parabolic relationship (inverted U curve) reported in adults of advanced age. Altered adrenal sensitivity may underlie this relationship.

Genetic association study of childhood aggression across raters, instruments, and age

Childhood aggressive behavior (AGG) has a substantial heritability of around 50%. Here we present a genome-wide association meta-analysis (GWAMA) of childhood AGG, in which all phenotype measures across childhood ages from multiple assessors were included. We analyzed phenotype assessments for a total of 328 935 observations from 87 485 children aged between 1.5 and 18 years, while accounting for sample overlap. We also meta-analyzed within subsets of the data, i.e., within rater, instrument and age. SNP-heritability for the overall meta-analysis (AGGoverall) was 3.31% (SE = 0.0038). We found no genome-wide significant SNPs for AGGoverall. The gene-based analysis returned three significant genes: ST3GAL3 (P = 1.6E-06), PCDH7 (P = 2.0E-06), and IPO13 (P = 2.5E-06). All three genes have previously been associated with educational traits. Polygenic scores based on our GWAMA significantly predicted aggression in a holdout sample of children (variance explained = 0.44%) and in retrospectively assessed childhood aggression (variance explained = 0.20%). Genetic correlations (rg) among rater-specific assessment of AGG ranged from rg = 0.46 between self- and teacher-assessment to rg = 0.81 between mother- and teacher-assessment. We obtained moderate-to-strong rgs with selected phenotypes from multiple domains, but hardly with any of the classical biomarkers thought to be associated with AGG. Significant genetic correlations were observed with most psychiatric and psychological traits (range [Formula: see text]: 0.19-1.00), except for obsessive-compulsive disorder. Aggression had a negative genetic correlation (rg = ~-0.5) with cognitive traits and age at first birth. Aggression was strongly genetically correlated with smoking phenotypes (range [Formula: see text]: 0.46-0.60). The genetic correlations between aggression and psychiatric disorders were weaker for teacher-reported AGG than for mother- and self-reported AGG. The current GWAMA of childhood aggression provides a powerful tool to interrogate the rater-specific genetic etiology of AGG.

The interactions between genetics and early childhood nutrition influence adult cardiometabolic risk factors

It is well established that genetics, environment, and interplay between them play a crucial role in adult disease. We aimed to evaluate the role of genetics, early life nutrition, and the interaction between them, on optimal adult health. As part of a large international consortium (n ~ 154,000), we identified 60 SNPs associated with both birthweight and adult disease. Utilising the Raine Study, we developed a birthweight polygenic score (BW-PGS) based on the 60 SNPs and examined relationships between BW-PGS and adulthood cardiovascular risk factors, specifically evaluating interactions with early life nutrition. Healthy nutrition was beneficial for all individuals; longer duration of any breastfeeding was particularly associated with lower BMI and lower Systolic Blood Pressure in those with higher BW-PGS. Optimal breastfeeding offers the greatest benefit to reduce adult obesity and hypertension in those genetically predisposed to high birthweight. This provides an example of how precision medicine in early life can improve adult health.

Sugar-Sweetened Beverage Consumption May Modify Associations Between Genetic Variants in the CHREBP (Carbohydrate Responsive Element Binding Protein) Locus and HDL-C (High-Density Lipoprotein Cholesterol) and Triglyceride Concentrations

Supplemental Digital Content is available in the text.

Background:

ChREBP (carbohydrate responsive element binding protein) is a transcription factor that responds to sugar consumption. Sugar-sweetened beverage (SSB) consumption and genetic variants in the CHREBP locus have separately been linked to HDL-C (high-density lipoprotein cholesterol) and triglyceride concentrations. We hypothesized that SSB consumption would modify the association between genetic variants in the CHREBP locus and dyslipidemia.

Methods:

Data from 11 cohorts from the Cohorts for Heart and Aging Research in Genomic Epidemiology consortium (N=63 599) and the UK Biobank (N=59 220) were used to quantify associations of SSB consumption, genetic variants, and their interaction on HDL-C and triglyceride concentrations using linear regression models. A total of 1606 single nucleotide polymorphisms within or near CHREBP were considered. SSB consumption was estimated from validated questionnaires, and participants were grouped by their estimated intake.

Results:

In a meta-analysis, rs71556729 was significantly associated with higher HDL-C concentrations only among the highest SSB consumers (β, 2.12 [95% CI, 1.16–3.07] mg/dL per allele; P<0.0001), but not significantly among the lowest SSB consumers (P=0.81; P_Diff <0.0001). Similar results were observed for 2 additional variants (rs35709627 and rs71556736). For triglyceride, rs55673514 was positively associated with triglyceride concentrations only among the highest SSB consumers (β, 0.06 [95% CI, 0.02–0.09] ln-mg/dL per allele, P=0.001) but not the lowest SSB consumers (P=0.84; P_Diff=0.0005).

Conclusions:

Our results identified genetic variants in the CHREBP locus that may protect against SSB-associated reductions in HDL-C and other variants that may exacerbate SSB-associated increases in triglyceride concentrations.

Registration:

URL: https://www.clinicaltrials.gov; Unique identifier: NCT00005133, NCT00005121, NCT00005487, and NCT00000479.

Effectiveness of combined vaginal progesterone and cervical cerclage in preventing preterm birth: a systematic review and meta-analysis protocol

INTRODUCTION

Preterm birth (PTB) is the leading cause of death in children under 5 years. Preventive therapies targeted towards women with risk factors such as a prior PTB or a short cervix reduce the rate of PTB. Cervical cerclage, vaginal progesterone and a combination of the two have been used with no consensus as to whether combined treatment is more effective than any single treatment alone. The objective of this review is to determine the efficacy of combined treatment compared with cerclage alone and combined treatment compared with progesterone alone.

METHODS AND ANALYSIS

Studies will be sourced from the electronic databases Medline (Ovid), EMBASE (Ovid), PsycINFO (Ovid), Scopus, CINAHL (EBSCOhost) and Cochrane Library (Wiley) and reference lists. We will not exclude any papers due to publication date. Randomised control trials (RCTs), non-RCTs and cohort studies assessing single therapy (either progesterone or cerclage) versus combined therapy in women with a singleton pregnancy will be included. Two independent reviewers will conduct study screening (at abstract and full-text level), data extraction and risk of bias assessment with disagreements resolved by an experienced researcher. Random or fixed effects models will be used depending on data heterogeneity and data will be presented as risk ratio for dichotomous data or mean difference for continuous data with a CI of 95% used for all outcomes.

ETHICS AND DISSEMINATION

Not applicable due to nature of the study type.

PROSPERO REGISTRATION NUMBER

CRD42020195975.

Variation in the SERPINA6/SERPINA1 locus alters morning plasma cortisol, hepatic corticosteroid binding globulin expression, gene expression in peripheral tissues, and risk of cardiovascular disease

The stress hormone cortisol modulates fuel metabolism, cardiovascular homoeostasis, mood, inflammation and cognition. The CORtisol NETwork (CORNET) consortium previously identified a single locus associated with morning plasma cortisol. Identifying additional genetic variants that explain more of the variance in cortisol could provide new insights into cortisol biology and provide statistical power to test the causative role of cortisol in common diseases. The CORNET consortium extended its genome-wide association meta-analysis for morning plasma cortisol from 12,597 to 25,314 subjects and from ~2.2 M to ~7 M SNPs, in 17 population-based cohorts of European ancestries. We confirmed the genetic association with SERPINA6/SERPINA1. This locus contains genes encoding corticosteroid binding globulin (CBG) and α1-antitrypsin. Expression quantitative trait loci (eQTL) analyses undertaken in the STARNET cohort of 600 individuals showed that specific genetic variants within the SERPINA6/SERPINA1 locus influence expression of SERPINA6 rather than SERPINA1 in the liver. Moreover, trans-eQTL analysis demonstrated effects on adipose tissue gene expression, suggesting that variations in CBG levels have an effect on delivery of cortisol to peripheral tissues. Two-sample Mendelian randomisation analyses provided evidence that each genetically-determined standard deviation (SD) increase in morning plasma cortisol was associated with increased odds of chronic ischaemic heart disease (0.32, 95% CI 0.06-0.59) and myocardial infarction (0.21, 95% CI 0.00-0.43) in UK Biobank and similarly in CARDIoGRAMplusC4D. These findings reveal a causative pathway for CBG in determining cortisol action in peripheral tissues and thereby contributing to the aetiology of cardiovascular disease.

The trans-ancestral genomic architecture of glycemic traits

Glycemic traits are used to diagnose and monitor type 2 diabetes and cardiometabolic health. To date, most genetic studies of glycemic traits have focused on individuals of European ancestry. Here we aggregated genome-wide association studies comprising up to 281,416 individuals without diabetes (30% non-European ancestry) for whom fasting glucose, 2-h glucose after an oral glucose challenge, glycated hemoglobin and fasting insulin data were available. Trans-ancestry and single-ancestry meta-analyses identified 242 loci (99 novel; P < 5 × 10-8), 80% of which had no significant evidence of between-ancestry heterogeneity. Analyses restricted to individuals of European ancestry with equivalent sample size would have led to 24 fewer new loci. Compared with single-ancestry analyses, equivalent-sized trans-ancestry fine-mapping reduced the number of estimated variants in 99% credible sets by a median of 37.5%. Genomic-feature, gene-expression and gene-set analyses revealed distinct biological signatures for each trait, highlighting different underlying biological pathways. Our results increase our understanding of diabetes pathophysiology by using trans-ancestry studies for improved power and resolution.

Variants associated with HHIP expression have sex-differential effects on lung function

Background: Lung function is highly heritable and differs between the sexes throughout life. However, little is known about sex-differential genetic effects on lung function. We aimed to conduct the first genome-wide genotype-by-sex interaction study on lung function to identify genetic effects that differ between males and females. Methods: We tested for interactions between 7,745,864 variants and sex on spirometry-based measures of lung function in UK Biobank (N=303,612), and sought replication in 75,696 independent individuals from the SpiroMeta consortium. Results: Five independent single-nucleotide polymorphisms (SNPs) showed genome-wide significant (P<5x10 -8) interactions with sex on lung function, and 21 showed suggestive interactions (P<1x10 -6). The strongest signal, from rs7697189 (chr4:145436894) on forced expiratory volume in 1 second (FEV 1) (P=3.15x10 -15), was replicated (P=0.016) in SpiroMeta. The C allele increased FEV 1 more in males (untransformed FEV 1 β=0.028 [SE 0.0022] litres) than females (β=0.009 [SE 0.0014] litres), and this effect was not accounted for by differential effects on height, smoking or pubertal age. rs7697189 resides upstream of the hedgehog-interacting protein ( HHIP) gene and was previously associated with lung function and HHIP lung expression. We found HHIP expression was significantly different between the sexes (P=6.90x10 -6), but we could not detect sex differential effects of rs7697189 on expression. Conclusions: We identified a novel genotype-by-sex interaction at a putative enhancer region upstream of the HHIP gene. Establishing the mechanism by which HHIP SNPs have different effects on lung function in males and females will be important for our understanding of lung health and diseases in both sexes.

Genetic Studies of Leptin Concentrations Implicate Leptin in the Regulation of Early Adiposity

Leptin influences food intake by informing the brain about the status of body fat stores. Rare LEP mutations associated with congenital leptin deficiency cause severe early-onset obesity that can be mitigated by administering leptin. However, the role of genetic regulation of leptin in polygenic obesity remains poorly understood. We performed an exome-based analysis in up to 57,232 individuals of diverse ancestries to identify genetic variants that influence adiposity-adjusted leptin concentrations. We identify five novel variants, including four missense variants, in LEP, ZNF800, KLHL31, and ACTL9, and one intergenic variant near KLF14. The missense variant Val94Met (rs17151919) in LEP was common in individuals of African ancestry only, and its association with lower leptin concentrations was specific to this ancestry (P = 2 × 10^-16, n = 3,901). Using in vitro analyses, we show that the Met94 allele decreases leptin secretion. We also show that the Met94 allele is associated with higher BMI in young African-ancestry children but not in adults, suggesting that leptin regulates early adiposity.

Novel loci for childhood body mass index and shared heritability with adult cardiometabolic traits

The genetic background of childhood body mass index (BMI), and the extent to which the well-known associations of childhood BMI with adult diseases are explained by shared genetic factors, are largely unknown. We performed a genome-wide association study meta-analysis of BMI in 61,111 children aged between 2 and 10 years. Twenty-five independent loci reached genome-wide significance in the combined discovery and replication analyses. Two of these, located near NEDD4L and SLC45A3, have not previously been reported in relation to either childhood or adult BMI. Positive genetic correlations of childhood BMI with birth weight and adult BMI, waist-to-hip ratio, diastolic blood pressure and type 2 diabetes were detected (Rg ranging from 0.11 to 0.76, P-values <0.002). A negative genetic correlation of childhood BMI with age at menarche was observed. Our results suggest that the biological processes underlying childhood BMI largely, but not completely, overlap with those underlying adult BMI. The well-known observational associations of BMI in childhood with cardio-metabolic diseases in adulthood may reflect partial genetic overlap, but in light of previous evidence, it is also likely that they are explained through phenotypic continuity of BMI from childhood into adulthood.

Variants associated with HHIP expression have sex-differential effects on lung function

Background: Lung function is highly heritable and differs between the sexes throughout life. However, little is known about sex-differential genetic effects on lung function. We aimed to conduct the first genome-wide genotype-by-sex interaction study on lung function to identify genetic effects that differ between males and females. Methods: We tested for interactions between 7,745,864 variants and sex on spirometry-based measures of lung function in UK Biobank (N=303,612), and sought replication in 75,696 independent individuals from the SpiroMeta consortium. Results: Five independent single-nucleotide polymorphisms (SNPs) showed genome-wide significant (P<5x10 -8) interactions with sex on lung function, and 21 showed suggestive interactions (P<1x10 -6). The strongest signal, from rs7697189 (chr4:145436894) on forced expiratory volume in 1 second (FEV 1) (P=3.15x10 -15), was replicated (P=0.016) in SpiroMeta. The C allele increased FEV 1 more in males (untransformed FEV 1 β=0.028 [SE 0.0022] litres) than females (β=0.009 [SE 0.0014] litres), and this effect was not accounted for by differential effects on height, smoking or pubertal age. rs7697189 resides upstream of the hedgehog-interacting protein ( HHIP) gene and was previously associated with lung function and HHIP lung expression. We found HHIP expression was significantly different between the sexes (P=6.90x10 -6), but we could not detect sex differential effects of rs7697189 on expression. Conclusions: We identified a novel genotype-by-sex interaction at a putative enhancer region upstream of the HHIP gene. Establishing the mechanism by which HHIP SNPs have different effects on lung function in males and females will be important for our understanding of lung health and diseases in both sexes.

Mendelian randomization analysis does not support causal associations of birth weight with hypertension risk and blood pressure in adulthood

Epidemiology studies suggested that low birthweight was associated with a higher risk of hypertension in later life. However, little is known about the causality of such associations. In our study, we evaluated the causal association of low birthweight with adulthood hypertension following a standard analytic protocol using the study-level data of 183,433 participants from 60 studies (CHARGE-BIG consortium), as well as that with blood pressure using publicly available summary-level genome-wide association data from EGG consortium of 153,781 participants, ICBP consortium and UK Biobank cohort together of 757,601 participants. We used seven SNPs as the instrumental variable in the study-level analysis and 47 SNPs in the summary-level analysis. In the study-level analyses, decreased birthweight was associated with a higher risk of hypertension in adults (the odds ratio per 1 standard deviation (SD) lower birthweight, 1.22; 95% CI 1.16 to 1.28), while no association was found between genetically instrumented birthweight and hypertension risk (instrumental odds ratio for causal effect per 1 SD lower birthweight, 0.97; 95% CI 0.68 to 1.41). Such results were consistent with that from the summary-level analyses, where the genetically determined low birthweight was not associated with blood pressure measurements either. One SD lower genetically determined birthweight was not associated with systolic blood pressure (β = - 0.76, 95% CI - 2.45 to 1.08 mmHg), 0.06 mmHg lower diastolic blood pressure (β = - 0.06, 95% CI - 0.93 to 0.87 mmHg), or pulse pressure (β = - 0.65, 95% CI - 1.38 to 0.69 mmHg, all p > 0.05). Our findings suggest that the inverse association of birthweight with hypertension risk from observational studies was not supported by large Mendelian randomization analyses.

Associations of autozygosity with a broad range of human phenotypes

In many species, the offspring of related parents suffer reduced reproductive success, a phenomenon known as inbreeding depression. In humans, the importance of this effect has remained unclear, partly because reproduction between close relatives is both rare and frequently associated with confounding social factors. Here, using genomic inbreeding coefficients (FROH) for >1.4 million individuals, we show that FROH is significantly associated (p < 0.0005) with apparently deleterious changes in 32 out of 100 traits analysed. These changes are associated with runs of homozygosity (ROH), but not with common variant homozygosity, suggesting that genetic variants associated with inbreeding depression are predominantly rare. The effect on fertility is striking: FROH equivalent to the offspring of first cousins is associated with a 55% decrease [95% CI 44-66%] in the odds of having children. Finally, the effects of FROH are confirmed within full-sibling pairs, where the variation in FROH is independent of all environmental confounding.

A trans-ancestral meta-analysis of genome-wide association studies reveals loci associated with childhood obesity

Although hundreds of genome-wide association studies-implicated loci have been reported for adult obesity-related traits, less is known about the genetics specific for early-onset obesity and with only a few studies conducted in non-European populations to date. Searching for additional genetic variants associated with childhood obesity, we performed a trans-ancestral meta-analysis of 30 studies consisting of up to 13 005 cases (≥95th percentile of body mass index (BMI) achieved 2-18 years old) and 15 599 controls (consistently <50th percentile of BMI) of European, African, North/South American and East Asian ancestry. Suggestive loci were taken forward for replication in a sample of 1888 cases and 4689 controls from seven cohorts of European and North/South American ancestry. In addition to observing 18 previously implicated BMI or obesity loci, for both early and late onset, we uncovered one completely novel locus in this trans-ancestral analysis (nearest gene, METTL15). The variant was nominally associated with only the European subgroup analysis but had a consistent direction of effect in other ethnicities. We then utilized trans-ancestral Bayesian analysis to narrow down the location of the probable causal variant at each genome-wide significant signal. Of all the fine-mapped loci, we were able to narrow down the causative variant at four known loci to fewer than 10 single nucleotide polymorphisms (SNPs) (FAIM2, GNPDA2, MC4R and SEC16B loci). In conclusion, an ethnically diverse setting has enabled us to both identify an additional pediatric obesity locus and further fine-map existing loci.

Association of Birth Weight With Type 2 Diabetes and Glycemic Traits: A Mendelian Randomization Study

IMPORTANCE

Observational studies have shown associations of birth weight with type 2 diabetes (T2D) and glycemic traits, but it remains unclear whether these associations represent causal associations.

OBJECTIVE

To test the association of birth weight with T2D and glycemic traits using a mendelian randomization analysis.

DESIGN, SETTING, AND PARTICIPANTS

This mendelian randomization study used a genetic risk score for birth weight that was constructed with 7 genome-wide significant single-nucleotide polymorphisms. The associations of this score with birth weight and T2D were tested in a mendelian randomization analysis using study-level data. The association of birth weight with T2D was tested using both study-level data (7 single-nucleotide polymorphisms were used as an instrumental variable) and summary-level data from the consortia (43 single-nucleotide polymorphisms were used as an instrumental variable). Data from 180 056 participants from 49 studies were included.

MAIN OUTCOMES AND MEASURES

Type 2 diabetes and glycemic traits.

RESULTS

This mendelian randomization analysis included 49 studies with 41 155 patients with T2D and 80 008 control participants from study-level data and 34 840 patients with T2D and 114 981 control participants from summary-level data. Study-level data showed that a 1-SD decrease in birth weight due to the genetic risk score was associated with higher risk of T2D among all participants (odds ratio [OR], 2.10; 95% CI, 1.69-2.61; P = 4.03 × 10-5), among European participants (OR, 1.96; 95% CI, 1.42-2.71; P = .04), and among East Asian participants (OR, 1.39; 95% CI, 1.18-1.62; P = .04). Similar results were observed from summary-level analyses. In addition, each 1-SD lower birth weight was associated with 0.189 SD higher fasting glucose concentration (β = 0.189; SE = 0.060; P = .002), but not with fasting insulin, 2-hour glucose, or hemoglobin A1c concentration.

CONCLUSIONS AND RELEVANCE

In this study, a genetic predisposition to lower birth weight was associated with increased risk of T2D and higher fasting glucose concentration, suggesting genetic effects on retarded fetal growth and increased diabetes risk that either are independent of each other or operate through alterations of integrated biological mechanisms.

Variants in the fetal genome near pro-inflammatory cytokine genes on 2q13 associate with gestational duration

The duration of pregnancy is influenced by fetal and maternal genetic and non-genetic factors. Here we report a fetal genome-wide association meta-analysis of gestational duration, and early preterm, preterm, and postterm birth in 84,689 infants. One locus on chromosome 2q13 is associated with gestational duration; the association is replicated in 9,291 additional infants (combined P = 3.96 × 10⁻¹⁴). Analysis of 15,588 mother-child pairs shows that the association is driven by fetal rather than maternal genotype. Functional experiments show that the lead SNP, rs7594852, alters the binding of the HIC1 transcriptional repressor. Genes at the locus include several interleukin 1 family members with roles in pro-inflammatory pathways that are central to the process of parturition. Further understanding of the underlying mechanisms will be of great public health importance, since giving birth either before or after the window of term gestation is associated with increased morbidity and mortality.

Gestational duration depends on both maternal and fetal genetic influences. Here, the authors perform a fetal genome-wide association meta-analysis and find that a locus on 2q13 is associated with pregnancy duration and further show that the lead SNP rs7594852 changes the binding properties of transcriptional repressor HIC1.

Exome-Derived Adiponectin-Associated Variants Implicate Obesity and Lipid Biology

Circulating levels of adiponectin, an adipocyte-secreted protein associated with cardiovascular and metabolic risk, are highly heritable. To gain insights into the biology that regulates adiponectin levels, we performed an exome array meta-analysis of 265,780 genetic variants in 67,739 individuals of European, Hispanic, African American, and East Asian ancestry. We identified 20 loci associated with adiponectin, including 11 that had been reported previously (p < 2 × 10^-7). Comparison of exome array variants to regional linkage disequilibrium (LD) patterns and prior genome-wide association study (GWAS) results detected candidate variants (r² > .60) spanning as much as 900 kb. To identify potential genes and mechanisms through which the previously unreported association signals act to affect adiponectin levels, we assessed cross-trait associations, expression quantitative trait loci in subcutaneous adipose, and biological pathways of nearby genes. Eight of the nine loci were also associated (p < 1 × 10^-4) with at least one obesity or lipid trait. Candidate genes include PRKAR2A, PTH1R, and HDAC9, which have been suggested to play roles in adipocyte differentiation or bone marrow adipose tissue. Taken together, these findings provide further insights into the processes that influence circulating adiponectin levels.

Publisher Correction: Protein-altering variants associated with body mass index implicate pathways that control energy intake and expenditure in obesity

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

Maternal and fetal genetic effects on birth weight and their relevance to cardio-metabolic risk factors

Birth weight variation is influenced by fetal and maternal genetic and non-genetic factors, and has been reproducibly associated with future cardio-metabolic health outcomes. In expanded genome-wide association analyses of own birth weight (n = 321,223) and offspring birth weight (n = 230,069 mothers), we identified 190 independent association signals (129 of which are novel). We used structural equation modeling to decompose the contributions of direct fetal and indirect maternal genetic effects, then applied Mendelian randomization to illuminate causal pathways. For example, both indirect maternal and direct fetal genetic effects drive the observational relationship between lower birth weight and higher later blood pressure: maternal blood pressure-raising alleles reduce offspring birth weight, but only direct fetal effects of these alleles, once inherited, increase later offspring blood pressure. Using maternal birth weight-lowering genotypes to proxy for an adverse intrauterine environment provided no evidence that it causally raises offspring blood pressure, indicating that the inverse birth weight-blood pressure association is attributable to genetic effects, and not to intrauterine programming.

Protein-coding variants implicate novel genes related to lipid homeostasis contributing to body-fat distribution

Body-fat distribution is a risk factor for adverse cardiovascular health consequences. We analyzed the association of body-fat distribution, assessed by waist-to-hip ratio adjusted for body mass index, with 228,985 predicted coding and splice site variants available on exome arrays in up to 344,369 individuals from five major ancestries (discovery) and 132,177 European-ancestry individuals (validation). We identified 15 common (minor allele frequency, MAF ≥5%) and nine low-frequency or rare (MAF <5%) coding novel variants. Pathway/gene set enrichment analyses identified lipid particle, adiponectin, abnormal white adipose tissue physiology and bone development and morphology as important contributors to fat distribution, while cross-trait associations highlight cardiometabolic traits. In functional follow-up analyses, specifically in Drosophila RNAi-knockdowns, we observed a significant increase in the total body triglyceride levels for two genes (DNAH10 and PLXND1). We implicate novel genes in fat distribution, stressing the importance of interrogating low-frequency and protein-coding variants.

Author Correction: Genome-wide association and HLA fine-mapping studies identify risk loci and genetic pathways underlying allergic rhinitis

In the version of this article initially published, in Fig. 3, the y-axis numbering did not match the log scale indicated in the axis label. The error has been corrected in the HTML and PDF version of the article.

Genome-wide association and HLA fine-mapping studies identify risk loci and genetic pathways underlying allergic rhinitis

Allergic rhinitis is the most common clinical presentation of allergy, affecting 400 million people worldwide, with increasing incidence in westernized countries^1,2. To elucidate the genetic architecture and understand the underlying disease mechanisms, we carried out a meta-analysis of allergic rhinitis in 59,762 cases and 152,358 controls of European ancestry and identified a total of 41 risk loci for allergic rhinitis, including 20 loci not previously associated with allergic rhinitis, which were confirmed in a replication phase of 60,720 cases and 618,527 controls. Functional annotation implicated genes involved in various immune pathways, and fine mapping of the HLA region suggested amino acid variants important for antigen binding. We further performed genome-wide association study (GWAS) analyses of allergic sensitization against inhalant allergens and nonallergic rhinitis, which suggested shared genetic mechanisms across rhinitis-related traits. Future studies of the identified loci and genes might identify novel targets for treatment and prevention of allergic rhinitis.

Genome-wide association study identifies nine novel loci for 2D:4D finger ratio, a putative retrospective biomarker of testosterone exposure in utero

The ratio of the length of the index finger to that of the ring finger (2D:4D) is sexually dimorphic and is commonly used as a non-invasive biomarker of prenatal androgen exposure. Most association studies of 2D:4D ratio with a diverse range of sex-specific traits have typically involved small sample sizes and have been difficult to replicate, raising questions around the utility and precise meaning of the measure. In the largest genome-wide association meta-analysis of 2D:4D ratio to date (N = 15 661, with replication N = 75 821), we identified 11 loci (9 novel) explaining 3.8% of the variance in mean 2D:4D ratio. We also found weak evidence for association (β = 0.06; P = 0.02) between 2D:4D ratio and sensitivity to testosterone [length of the CAG microsatellite repeat in the androgen receptor (AR) gene] in females only. Furthermore, genetic variants associated with (adult) testosterone levels and/or sex hormone-binding globulin were not associated with 2D:4D ratio in our sample. Although we were unable to find strong evidence from our genetic study to support the hypothesis that 2D:4D ratio is a direct biomarker of prenatal exposure to androgens in healthy individuals, our findings do not explicitly exclude this possibility, and pathways involving testosterone may become apparent as the size of the discovery sample increases further. Our findings provide new insight into the underlying biology shaping 2D:4D variation in the general population.