Fetal alleles predisposing to metabolically favorable adiposity are associated with higher birth weight.

BACKGROUND: Higher birthweight is associated with higher adult body mass index (BMI). Alleles that predispose to greater adult adiposity might act in fetal life to increase fetal growth and birthweight. Whether there are fetal effects of recently identified adult metabolically favorable adiposity alleles on birthweight is unknown. AIM: We aimed to test the effect on birthweight of fetal genetic predisposition to higher metabolically favorable adult adiposity and compare that with the effect of fetal genetic predisposition to higher adult BMI. METHODS: We used published genome wide association study data (n = upto 406 063) to estimate fetal effects on birthweight (adjusting for maternal genotype) of alleles known to raise metabolically favorable adult adiposity or BMI. We combined summary data across single nucleotide polymorphisms (SNPs) with random effects meta-analyses. We performed weighted linear regression of SNP-birthweight effects against SNP-adult adiposity effects to test for a dose-dependent association. RESULTS: Fetal genetic predisposition to higher metabolically favorable adult adiposity and higher adult BMI were both associated with higher birthweight (3 g per effect allele (95% CI: 1-5) averaged over 14 SNPs; P = 0.002; 0.5 g per effect allele (95% CI: 0-1) averaged over 76 SNPs; P = 0.042, respectively). SNPs with greater effects on metabolically favorable adiposity tended to have greater effects on birthweight (R2 = 0.2912, P = 0.027). There was no dose-dependent association for BMI (R2 = -0.0019, P = 0.602). CONCLUSIONS: Fetal genetic predisposition to both higher adult metabolically favorable adiposity and BMI is associated with birthweight. Fetal effects of metabolically favorable adiposity-raising alleles on birthweight are modestly proportional to their effects on future adiposity, but those of BMI-raising alleles are not.

Meta-regression of genome-wide association studies to estimate age-varying genetic effects.

Fixed-effect meta-analysis has been used to summarize genetic effects on a phenotype across multiple Genome-Wide Association Studies (GWAS) assuming a common underlying genetic effect. Genetic effects may vary with age (or other characteristics), and not allowing for this in a GWAS might lead to bias. Meta-regression models between study heterogeneity and allows effect modification of the genetic effects to be explored. The aim of this study was to explore the use of meta-analysis and meta-regression for estimating age-varying genetic effects on phenotypes. With simulations we compared the performance of meta-regression to fixed-effect and random -effects meta-analyses in estimating (i) main genetic effects and (ii) age-varying genetic effects (SNP by age interactions) from multiple GWAS studies under a range of scenarios. We applied meta-regression on publicly available summary data to estimate the main and age-varying genetic effects of the FTO SNP rs9939609 on Body Mass Index (BMI). Fixed-effect and random-effects meta-analyses accurately estimated genetic effects when these did not change with age. Meta-regression accurately estimated both main genetic effects and age-varying genetic effects. When the number of studies or the age-diversity between studies was low, meta-regression had limited power. In the applied example, each additional minor allele (A) of rs9939609 was inversely associated with BMI at ages 0 to 3, and positively associated at ages 5.5 to 13. Our findings challenge the assumption that genetic effects are consistent across all ages and provide a method for exploring this. GWAS consortia should be encouraged to use meta-regression to explore age-varying genetic effects.

Using Genomic Structural Equation Modeling to Partition the Genetic Covariance Between Birthweight and Cardiometabolic Risk Factors into Maternal and Offspring Components in the Norwegian HUNT Study.

The Barker Hypothesis posits that adverse intrauterine environments result in fetal growth restriction and increased risk of cardiometabolic disease through developmental compensations. Here we introduce a new statistical model using the genomic SEM software that is capable of simultaneously partitioning the genetic covariation between birthweight and cardiometabolic traits into maternally mediated and offspring mediated contributions. We model the covariance between birthweight and later life outcomes, such as blood pressure, non-fasting glucose, blood lipids and body mass index in the Norwegian HUNT study, consisting of 15,261 mother-eldest offspring pairs with genetic and phenotypic data. Application of this model showed some evidence for maternally mediated effects of systolic blood pressure on offspring birthweight, and pleiotropy between birthweight and non-fasting glucose mediated through the offspring genome. This underscores the importance of genetic links between birthweight and cardiometabolic phenotypes and offer alternative explanations to environmentally based hypotheses for the phenotypic correlation between these variables.

Estimating indirect parental genetic effects on offspring phenotypes using virtual parental genotypes derived from sibling and half sibling pairs.

Indirect parental genetic effects may be defined as the influence of parental genotypes on offspring phenotypes over and above that which results from the transmission of genes from parents to their children. However, given the relative paucity of large-scale family-based cohorts around the world, it is difficult to demonstrate parental genetic effects on human traits, particularly at individual loci. In this manuscript, we illustrate how parental genetic effects on offspring phenotypes, including late onset conditions, can be estimated at individual loci in principle using large-scale genome-wide association study (GWAS) data, even in the absence of parental genotypes. Our strategy involves creating "virtual" mothers and fathers by estimating the genotypic dosages of parental genotypes using physically genotyped data from relative pairs. We then utilize the expected dosages of the parents, and the actual genotypes of the offspring relative pairs, to perform conditional genetic association analyses to obtain asymptotically unbiased estimates of maternal, paternal and offspring genetic effects. We apply our approach to 19066 sibling pairs from the UK Biobank and show that a polygenic score consisting of imputed parental educational attainment SNP dosages is strongly related to offspring educational attainment even after correcting for offspring genotype at the same loci. We develop a freely available web application that quantifies the power of our approach using closed form asymptotic solutions. We implement our methods in a user-friendly software package IMPISH (IMputing Parental genotypes In Siblings and Half Siblings) which allows users to quickly and efficiently impute parental genotypes across the genome in large genome-wide datasets, and then use these estimated dosages in downstream linear mixed model association analyses. We conclude that imputing parental genotypes from relative pairs may provide a useful adjunct to existing large-scale genetic studies of parents and their offspring.

Mendelian randomization analysis of factors related to ovulation and reproductive function and endometrial cancer risk.

BACKGROUND: Observational epidemiological studies suggest a link between several factors related to ovulation and reproductive function and endometrial cancer (EC) risk; however, it is not clear whether these relationships are causal, and whether the risk factors act independently of each other. The aim of this study was to investigate putative causal relationships between the number of live births, age at last live birth, and years ovulating and EC risk.  METHODS: We conducted a series of observational analyses to investigate various risk factors and EC risk in the UK Biobank (UKBB). Additionally, multivariate analysis was performed to elucidate the relationship between the number of live births, age at last live birth, and years ovulating and other related factors such as age at natural menopause, age at menarche, and body mass index (BMI). Secondly, we used Mendelian randomization (MR) to assess if these observed relationships were causal. Genome-wide significant single nucleotide polymorphisms (SNPs) were extracted from previous studies of woman's number of live births, age at menopause and menarche, and BMI. We conducted a genome-wide association analysis using the UKBB to identify SNPs associated with years ovulating, years using the contraceptive pill, and age at last live birth. RESULTS: We found evidence for a causal effect of the number of live births (inverse variance weighted (IVW) odds ratio (OR): 0.537, p = 0.006), the number of years ovulating (IVW OR: 1.051, p = 0.014), in addition to the known risk factors BMI, age at menarche, and age at menopause on EC risk in the univariate MR analyses. Due to the close relationships between these factors, we followed up with multivariable MR (MVMR) analysis. Results from the MVMR analysis showed that number of live births had a causal effect on EC risk (OR: 0.783, p = 0.036) independent of BMI, age at menarche and age at menopause. CONCLUSIONS: MVMR analysis showed that the number of live births causally reduced the risk of EC.

Genome-wide associations for birth weight and correlations with adult disease.

Birth weight (BW) has been shown to be influenced by both fetal and maternal factors and in observational studies is reproducibly associated with future risk of adult metabolic diseases including type 2 diabetes (T2D) and cardiovascular disease. These life-course associations have often been attributed to the impact of an adverse early life environment. Here, we performed a multi-ancestry genome-wide association study (GWAS) meta-analysis of BW in 153,781 individuals, identifying 60 loci where fetal genotype was associated with BW (P < 5 × 10-8). Overall, approximately 15% of variance in BW was captured by assays of fetal genetic variation. Using genetic association alone, we found strong inverse genetic correlations between BW and systolic blood pressure (Rg = -0.22, P = 5.5 × 10-13), T2D (Rg = -0.27, P = 1.1 × 10-6) and coronary artery disease (Rg = -0.30, P = 6.5 × 10-9). In addition, using large -cohort datasets, we demonstrated that genetic factors were the major contributor to the negative covariance between BW and future cardiometabolic risk. Pathway analyses indicated that the protein products of genes within BW-associated regions were enriched for diverse processes including insulin signalling, glucose homeostasis, glycogen biosynthesis and chromatin remodelling. There was also enrichment of associations with BW in known imprinted regions (P = 1.9 × 10-4). We demonstrate that life-course associations between early growth phenotypes and adult cardiometabolic disease are in part the result of shared genetic effects and identify some of the pathways through which these causal genetic effects are mediated.

Higher maternal adiposity reduces offspring birthweight if associated with a metabolically favourable profile.

AIMS/HYPOTHESIS: Higher maternal BMI during pregnancy is associated with higher offspring birthweight, but it is not known whether this is solely the result of adverse metabolic consequences of higher maternal adiposity, such as maternal insulin resistance and fetal exposure to higher glucose levels, or whether there is any effect of raised adiposity through non-metabolic (e.g. mechanical) factors. We aimed to use genetic variants known to predispose to higher adiposity, coupled with a favourable metabolic profile, in a Mendelian randomisation (MR) study comparing the effect of maternal 'metabolically favourable adiposity' on offspring birthweight with the effect of maternal general adiposity (as indexed by BMI). METHODS: To test the causal effects of maternal metabolically favourable adiposity or general adiposity on offspring birthweight, we performed two-sample MR. We used variants identified in large, published genetic-association studies as being associated with either higher adiposity and a favourable metabolic profile, or higher BMI (n = 442,278 and n = 322,154 for metabolically favourable adiposity and BMI, respectively). We then extracted data on the metabolically favourable adiposity and BMI variants from a large, published genetic-association study of maternal genotype and offspring birthweight controlling for fetal genetic effects (n = 406,063 with maternal and/or fetal genotype effect estimates). We used several sensitivity analyses to test the reliability of the results. As secondary analyses, we used data from four cohorts (total n = 9323 mother-child pairs) to test the effects of maternal metabolically favourable adiposity or BMI on maternal gestational glucose, anthropometric components of birthweight and cord-blood biomarkers. RESULTS: Higher maternal adiposity with a favourable metabolic profile was associated with lower offspring birthweight (-94 [95% CI -150, -38] g per 1 SD [6.5%] higher maternal metabolically favourable adiposity, p = 0.001). By contrast, higher maternal BMI was associated with higher offspring birthweight (35 [95% CI 16, 53] g per 1 SD [4 kg/m2] higher maternal BMI, p = 0.0002). Sensitivity analyses were broadly consistent with the main results. There was evidence of outlier SNPs for both exposures; their removal slightly strengthened the metabolically favourable adiposity estimate and made no difference to the BMI estimate. Our secondary analyses found evidence to suggest that a higher maternal metabolically favourable adiposity decreases pregnancy fasting glucose levels while a higher maternal BMI increases them. The effects on neonatal anthropometric traits were consistent with the overall effect on birthweight but the smaller sample sizes for these analyses meant that the effects were imprecisely estimated. We also found evidence to suggest that higher maternal metabolically favourable adiposity decreases cord-blood leptin while higher maternal BMI increases it. CONCLUSIONS/INTERPRETATION: Our results show that higher adiposity in mothers does not necessarily lead to higher offspring birthweight. Higher maternal adiposity can lead to lower offspring birthweight if accompanied by a favourable metabolic profile. DATA AVAILABILITY: The data for the genome-wide association studies (GWAS) of BMI are available at https://portals.broadinstitute.org/collaboration/giant/index.php/GIANT_consortium_data_files . The data for the GWAS of body fat percentage are available at https://walker05.u.hpc.mssm.edu .

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.

Genome-wide association study of offspring birth weight in 86 577 women identifies five novel loci and highlights maternal genetic effects that are independent of fetal genetics.

Genome-wide association studies of birth weight have focused on fetal genetics, whereas relatively little is known about the role of maternal genetic variation. We aimed to identify maternal genetic variants associated with birth weight that could highlight potentially relevant maternal determinants of fetal growth. We meta-analysed data on up to 8.7 million SNPs in up to 86 577 women of European descent from the Early Growth Genetics (EGG) Consortium and the UK Biobank. We used structural equation modelling (SEM) and analyses of mother-child pairs to quantify the separate maternal and fetal genetic effects. Maternal SNPs at 10 loci (MTNR1B, HMGA2, SH2B3, KCNAB1, L3MBTL3, GCK, EBF1, TCF7L2, ACTL9, CYP3A7) were associated with offspring birth weight at P < 5 × 10-8. In SEM analyses, at least 7 of the 10 associations were consistent with effects of the maternal genotype acting via the intrauterine environment, rather than via effects of shared alleles with the fetus. Variants, or correlated proxies, at many of the loci had been previously associated with adult traits, including fasting glucose (MTNR1B, GCK and TCF7L2) and sex hormone levels (CYP3A7), and one (EBF1) with gestational duration. The identified associations indicate that genetic effects on maternal glucose, cytochrome P450 activity and gestational duration, and potentially on maternal blood pressure and immune function, are relevant for fetal growth. Further characterization of these associations in mechanistic and causal analyses will enhance understanding of the potentially modifiable maternal determinants of fetal growth, with the goal of reducing the morbidity and mortality associated with low and high birth weights.

Introducing M-GCTA a Software Package to Estimate Maternal (or Paternal) Genetic Effects on Offspring Phenotypes.

There is increasing interest within the genetics community in estimating the relative contribution of parental genetic effects on offspring phenotypes. Here we describe the user-friendly M-GCTA software package used to estimate the proportion of phenotypic variance explained by maternal (or alternatively paternal) and offspring genotypes on offspring phenotypes. The tool requires large studies where genome-wide genotype data are available on mother- (or alternatively father-) offspring pairs. The software includes several options for data cleaning and quality control, including the ability to detect and automatically remove cryptically related pairs of individuals. It also allows users to construct genetic relationship matrices indexing genetic similarity across the genome between parents and offspring, enabling the estimation of variance explained by maternal (or alternatively paternal) and offspring genetic effects. We evaluated the performance of the software using a range of data simulations and estimated the computing time and memory requirements. We demonstrate the use of M-GCTA on previously analyzed birth weight data from two large population based birth cohorts, the Avon Longitudinal Study of Parents and Children (ALSPAC) and the Norwegian Mother and Child Cohort Study (MoBa). We show how genetic variation in birth weight is predominantly explained by fetal genetic rather than maternal genetic sources of variation.

Calculating Power to Detect Maternal and Offspring Genetic Effects in Genetic Association Studies.

Offspring outcomes are a function of maternal genetics operating on the intrauterine and postnatal environment, offspring genetics and environmental factors. Partitioning genetic effects into maternal and offspring components requires genotyped mother-offspring pairs or genotyped individuals with phenotypic information on themselves and their offspring. We performed asymptotic power calculations and data simulations to estimate power to detect maternal and offspring genetic effects under a range of different study designs and models. We also developed the "Maternal and offspring Genetic effects Power Calculator" (M-GPC), an online utility which allows users to estimate the power to detect maternal and offspring genetic effects in their own studies. We find that approximately 50,000 genotyped mother-offspring pairs will be required to detect realistically sized maternal or offspring genetic effects (> 0.1% variance explained) with appreciable power (power > 90%, α = 5 × 10-8, two degree of freedom test), whereas greater than 10,000 pairs will be required to determine whether known genetic loci have maternal and/or offspring genetic effects (power > 78%, α = 0.05). The structural equation modelling framework espoused in this manuscript provides a natural method of combining different data structures including those present in large scale biobanks in order to maximize power to detect maternal and offspring genetic effects. We conclude that the sample sizes required to detect maternal or offspring genetic effects that explain realistic proportions of the trait variance with appreciable power are achievable and within the range of current research consortia.

LD Hub: a centralized database and web interface to perform LD score regression that maximizes the potential of summary level GWAS data for SNP heritability and genetic correlation analysis.

MOTIVATION: LD score regression is a reliable and efficient method of using genome-wide association study (GWAS) summary-level results data to estimate the SNP heritability of complex traits and diseases, partition this heritability into functional categories, and estimate the genetic correlation between different phenotypes. Because the method relies on summary level results data, LD score regression is computationally tractable even for very large sample sizes. However, publicly available GWAS summary-level data are typically stored in different databases and have different formats, making it difficult to apply LD score regression to estimate genetic correlations across many different traits simultaneously. RESULTS: In this manuscript, we describe LD Hub - a centralized database of summary-level GWAS results for 173 diseases/traits from different publicly available resources/consortia and a web interface that automates the LD score regression analysis pipeline. To demonstrate functionality and validate our software, we replicated previously reported LD score regression analyses of 49 traits/diseases using LD Hub; and estimated SNP heritability and the genetic correlation across the different phenotypes. We also present new results obtained by uploading a recent atopic dermatitis GWAS meta-analysis to examine the genetic correlation between the condition and other potentially related traits. In response to the growing availability of publicly accessible GWAS summary-level results data, our database and the accompanying web interface will ensure maximal uptake of the LD score regression methodology, provide a useful database for the public dissemination of GWAS results, and provide a method for easily screening hundreds of traits for overlapping genetic aetiologies. AVAILABILITY AND IMPLEMENTATION: The web interface and instructions for using LD Hub are available at http://ldsc.broadinstitute.org/ CONTACT: jie.zheng@bristol.ac.ukSupplementary information: Supplementary data are available at Bioinformatics online.

Shared genetic variants suggest common pathways in allergy and autoimmune diseases.

BACKGROUND: The relationship between allergy and autoimmune disorders is complex and poorly understood. OBJECTIVE: We sought to investigate commonalities in genetic loci and pathways between allergy and autoimmune diseases to elucidate shared disease mechanisms. METHODS: We meta-analyzed 2 genome-wide association studies on self-reported allergy and sensitization comprising a total of 62,330 subjects. These results were used to calculate enrichment for single nucleotide polymorphisms (SNPs) previously associated with autoimmune diseases. Furthermore, we probed for enrichment within genetic pathways and of transcription factor binding sites and characterized commonalities in variant burden on tissue-specific regulatory sites by calculating the enrichment of allergy SNPs falling in gene regulatory regions in various cells using Encode Roadmap DNase-hypersensitive site data. Finally, we compared the allergy data with those of all known diseases. RESULTS: Among 290 loci previously associated with 16 autoimmune diseases, we found a significant enrichment of loci also associated with allergy (P = 1.4e-17) encompassing 29 loci at a false discovery rate of less than 0.05. Such enrichment seemed to be a general characteristic for autoimmune diseases. Among the common loci, 48% had the same direction of effect for allergy and autoimmune diseases. Additionally, we observed an enrichment of allergy SNPs falling within immune pathways and regions of chromatin accessible in immune cells that was also represented in patients with autoimmune diseases but not those with other diseases. CONCLUSION: We identified shared susceptibility loci and commonalities in pathways between allergy and autoimmune diseases, suggesting shared disease mechanisms. Further studies of these shared genetic mechanisms might help in understanding the complex relationship between these diseases, including the parallel increase in disease prevalence.

Genetic variants in adult bone mineral density and fracture risk genes are associated with the rate of bone mineral density acquisition in adolescence.

Previous studies have identified 63 single-nucleotide polymorphisms (SNPs) associated with bone mineral density (BMD) in adults. These SNPs are thought to reflect variants that influence bone maintenance and/or loss in adults. It is unclear whether they affect the rate of bone acquisition during adolescence. Bone measurements and genetic data were available on 6397 individuals from the Avon Longitudinal Study of Parents and Children at up to five follow-up clinics. Linear mixed effects models with smoothing splines were used for longitudinal modelling of BMD and its components bone mineral content (BMC) and bone area (BA), from 9 to 17 years. Genotype data from the 63 adult BMD associated SNPs were investigated individually and as a genetic risk score in the longitudinal model. Each additional BMD lowering allele of the genetic risk score was associated with lower BMD at age 13 [per allele effect size, 0.002 g/cm(2) (SE = 0.0001, P = 1.24 × 10(-38))] and decreased BMD acquisition from 9 to 17 years (P = 9.17 × 10(-7)). This association was driven by changes in BMC rather than BA. The genetic risk score explained ∼2% of the variation in BMD at 9 and 17 years, a third of that explained in adults (6%). Genetic variants that putatively affect bone maintenance and/or loss in adults appear to have a small influence on the rate of bone acquisition through adolescence.

Genome-wide association study of blood lead shows multiple associations near ALAD.

Exposure to high levels of environmental lead, or biomarker evidence of high body lead content, is associated with anaemia, developmental and neurological deficits in children, and increased mortality in adults. Adverse effects of lead still occur despite substantial reduction in environmental exposure. There is genetic variation between individuals in blood lead concentration but the polymorphisms contributing to this have not been defined. We measured blood or erythrocyte lead content, and carried out genome-wide association analysis, on population-based cohorts of adult volunteers from Australia and UK (N = 5433). Samples from Australia were collected in two studies, in 1993-1996 and 2002-2005 and from UK in 1991-1992. One locus, at ALAD on chromosome 9, showed consistent association with blood lead across countries and evidence for multiple independent allelic effects. The most significant single nucleotide polymorphism (SNP), rs1805313 (P = 3.91 × 10(-14) for lead concentration in a meta-analysis of all data), is known to have effects on ALAD expression in blood cells but other SNPs affecting ALAD expression did not affect blood lead. Variants at 12 other loci, including ABO, showed suggestive associations (5 × 10(-6) > P > 5 × 10(-8)). Identification of genetic polymorphisms affecting blood lead reinforces the view that genetic factors, as well as environmental ones, are important in determining blood lead levels. The ways in which ALAD variation affects lead uptake or distribution are still to be determined.

A novel common variant in DCST2 is associated with length in early life and height in adulthood.

Common genetic variants have been identified for adult height, but not much is known about the genetics of skeletal growth in early life. To identify common genetic variants that influence fetal skeletal growth, we meta-analyzed 22 genome-wide association studies (Stage 1; N = 28 459). We identified seven independent top single nucleotide polymorphisms (SNPs) (P < 1 × 10(-6)) for birth length, of which three were novel and four were in or near loci known to be associated with adult height (LCORL, PTCH1, GPR126 and HMGA2). The three novel SNPs were followed-up in nine replication studies (Stage 2; N = 11 995), with rs905938 in DC-STAMP domain containing 2 (DCST2) genome-wide significantly associated with birth length in a joint analysis (Stages 1 + 2; ß = 0.046, SE = 0.008, P = 2.46 × 10(-8), explained variance = 0.05%). Rs905938 was also associated with infant length (N = 28 228; P = 5.54 × 10(-4)) and adult height (N = 127 513; P = 1.45 × 10(-5)). DCST2 is a DC-STAMP-like protein family member and DC-STAMP is an osteoclast cell-fusion regulator. Polygenic scores based on 180 SNPs previously associated with human adult stature explained 0.13% of variance in birth length. The same SNPs explained 2.95% of the variance of infant length. Of the 180 known adult height loci, 11 were genome-wide significantly associated with infant length (SF3B4, LCORL, SPAG17, C6orf173, PTCH1, GDF5, ZNFX1, HHIP, ACAN, HLA locus and HMGA2). This study highlights that common variation in DCST2 influences variation in early growth and adult height.

Phenotypic dissection of bone mineral density reveals skeletal site specificity and facilitates the identification of novel loci in the genetic regulation of bone mass attainment.

Heritability of bone mineral density (BMD) varies across skeletal sites, reflecting different relative contributions of genetic and environmental influences. To quantify the degree to which common genetic variants tag and environmental factors influence BMD, at different sites, we estimated the genetic (rg) and residual (re) correlations between BMD measured at the upper limbs (UL-BMD), lower limbs (LL-BMD) and skull (SK-BMD), using total-body DXA scans of ∼ 4,890 participants recruited by the Avon Longitudinal Study of Parents and their Children (ALSPAC). Point estimates of rg indicated that appendicular sites have a greater proportion of shared genetic architecture (LL-/UL-BMD rg = 0.78) between them, than with the skull (UL-/SK-BMD rg = 0.58 and LL-/SK-BMD rg = 0.43). Likewise, the residual correlation between BMD at appendicular sites (r(e) = 0.55) was higher than the residual correlation between SK-BMD and BMD at appendicular sites (r(e) = 0.20-0.24). To explore the basis for the observed differences in rg and re, genome-wide association meta-analyses were performed (n ∼ 9,395), combining data from ALSPAC and the Generation R Study identifying 15 independent signals from 13 loci associated at genome-wide significant level across different skeletal regions. Results suggested that previously identified BMD-associated variants may exert site-specific effects (i.e. differ in the strength of their association and magnitude of effect across different skeletal sites). In particular, variants at CPED1 exerted a larger influence on SK-BMD and UL-BMD when compared to LL-BMD (P = 2.01 × 10(-37)), whilst variants at WNT16 influenced UL-BMD to a greater degree when compared to SK- and LL-BMD (P = 2.31 × 10(-14)). In addition, we report a novel association between RIN3 (previously associated with Paget's disease) and LL-BMD (rs754388: ß = 0.13, SE = 0.02, P = 1.4 × 10(-10)). Our results suggest that BMD at different skeletal sites is under a mixture of shared and specific genetic and environmental influences. Allowing for these differences by performing genome-wide association at different skeletal sites may help uncover new genetic influences on BMD.

Genome-wide association and longitudinal analyses reveal genetic loci linking pubertal height growth, pubertal timing and childhood adiposity.

The pubertal height growth spurt is a distinctive feature of childhood growth reflecting both the central onset of puberty and local growth factors. Although little is known about the underlying genetics, growth variability during puberty correlates with adult risks for hormone-dependent cancer and adverse cardiometabolic health. The only gene so far associated with pubertal height growth, LIN28B, pleiotropically influences childhood growth, puberty and cancer progression, pointing to shared underlying mechanisms. To discover genetic loci influencing pubertal height and growth and to place them in context of overall growth and maturation, we performed genome-wide association meta-analyses in 18 737 European samples utilizing longitudinally collected height measurements. We found significant associations (P < 1.67 × 10(-8)) at 10 loci, including LIN28B. Five loci associated with pubertal timing, all impacting multiple aspects of growth. In particular, a novel variant correlated with expression of MAPK3, and associated both with increased prepubertal growth and earlier menarche. Another variant near ADCY3-POMC associated with increased body mass index, reduced pubertal growth and earlier puberty. Whereas epidemiological correlations suggest that early puberty marks a pathway from rapid prepubertal growth to reduced final height and adult obesity, our study shows that individual loci associating with pubertal growth have variable longitudinal growth patterns that may differ from epidemiological observations. Overall, this study uncovers part of the complex genetic architecture linking pubertal height growth, the timing of puberty and childhood obesity and provides new information to pinpoint processes linking these traits.

Robustness of the linear mixed effects model to error distribution assumptions and the consequences for genome-wide association studies.

Genome-wide association studies have been successful in uncovering novel genetic variants that are associated with disease status or cross-sectional phenotypic traits. Researchers are beginning to investigate how genes play a role in the development of a trait over time. Linear mixed effects models (LMM) are commonly used to model longitudinal data; however, it is unclear if the failure to meet the models distributional assumptions will affect the conclusions when conducting a genome-wide association study. In an extensive simulation study, we compare coverage probabilities, bias, type 1 error rates and statistical power when the error of the LMM is either heteroscedastic or has a non-Gaussian distribution. We conclude that the model is robust to misspecification if the same function of age is included in the fixed and random effects. However, type 1 error of the genetic effect over time is inflated, regardless of the model misspecification, if the polynomial function for age in the fixed and random effects differs. In situations where the model will not converge with a high order polynomial function in the random effects, a reduced function can be used but a robust standard error needs to be calculated to avoid inflation of the type 1 error. As an illustration, a LMM was applied to longitudinal body mass index (BMI) data over childhood in the ALSPAC cohort; the results emphasised the need for the robust standard error to ensure correct inference of associations of longitudinal BMI with chromosome 16 single nucleotide polymorphisms.