The genetic landscape of Scotland and the Isles.

Britain and Ireland are known to show population genetic structure; however, large swathes of Scotland, in particular, have yet to be described. Delineating the structure and ancestry of these populations will allow variant discovery efforts to focus efficiently on areas not represented in existing cohorts. Thus, we assembled genotype data for 2,554 individuals from across the entire archipelago with geographically restricted ancestry, and performed population structure analyses and comparisons to ancient DNA. Extensive geographic structuring is revealed, from broad scales such as a NE to SW divide in mainland Scotland, through to the finest scale observed to date: across 3 km in the Northern Isles. Many genetic boundaries are consistent with Dark Age kingdoms of Gaels, Picts, Britons, and Norse. Populations in the Hebrides, the Highlands, Argyll, Donegal, and the Isle of Man show characteristics of isolation. We document a pole of Norwegian ancestry in the north of the archipelago (reaching 23 to 28% in Shetland) which complements previously described poles of Germanic ancestry in the east, and "Celtic" to the west. This modern genetic structure suggests a northwestern British or Irish source population for the ancient Gaels that contributed to the founding of Iceland. As rarer variants, often with larger effect sizes, become the focus of complex trait genetics, more diverse rural cohorts may be required to optimize discoveries in British and Irish populations and their considerable global diaspora.

Genome-wide survey of parent-of-origin effects on DNA methylation identifies candidate imprinted loci in humans.

Genomic imprinting is an epigenetic mechanism leading to parent-of-origin silencing of alleles. So far, the precise number of imprinted regions in humans is uncertain. In this study, we leveraged genome-wide DNA methylation in whole blood measured longitudinally at three time points (birth, childhood and adolescence) and genome-wide association studies (GWAS) data in 740 mother-child duos from the Avon Longitudinal Study of parents and children to identify candidate imprinted loci. We reasoned that cis-meQTLs at genomic regions that were imprinted would show strong evidence of parent-of-origin associations with DNA methylation, enabling the detection of imprinted regions. Using this approach, we identified genome-wide significant cis-meQTLs that exhibited parent-of-origin effects (POEs) at 82 loci, 34 novel and 48 regions previously implicated in imprinting (3.7-10

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.

A genome-wide association study of mitochondrial DNA copy number in two population-based cohorts.

BACKGROUND: Mitochondrial DNA copy number (mtDNA CN) exhibits interindividual and intercellular variation, but few genome-wide association studies (GWAS) of directly assayed mtDNA CN exist. We undertook a GWAS of qPCR-assayed mtDNA CN in the Avon Longitudinal Study of Parents and Children (ALSPAC) and the UK Blood Service (UKBS) cohort. After validating and harmonising data, 5461 ALSPAC mothers (16-43 years at mtDNA CN assay) and 1338 UKBS females (17-69 years) were included in a meta-analysis. Sensitivity analyses restricted to females with white cell-extracted DNA and adjusted for estimated or assayed cell proportions. Associations were also explored in ALSPAC children and UKBS males. RESULTS: A neutrophil-associated locus approached genome-wide significance (rs709591 [MED24], ß (change in SD units of mtDNA CN per allele) [SE] - 0.084 [0.016], p = 1.54e-07) in the main meta-analysis of adult females. This association was concordant in magnitude and direction in UKBS males and ALSPAC neonates. SNPs in and around ABHD8 were associated with mtDNA CN in ALSPAC neonates (rs10424198, ß [SE] 0.262 [0.034], p = 1.40e-14), but not other study groups. In a meta-analysis of unrelated individuals (N = 11,253), we replicated a published association in TFAM (ß [SE] 0.046 [0.017], p = 0.006), with an effect size much smaller than that observed in the replication analysis of a previous in silico GWAS. CONCLUSIONS: In a hypothesis-generating GWAS, we confirm an association between TFAM and mtDNA CN and present putative loci requiring replication in much larger samples. We discuss the limitations of our work, in terms of measurement error and cellular heterogeneity, and highlight the need for larger studies to better understand nuclear genomic control of mtDNA copy number.

Association of maternal circulating 25(OH)D and calcium with birth weight: A mendelian randomisation analysis.

BACKGROUND: Systematic reviews of randomised controlled trials (RCTs) have suggested that maternal vitamin D (25[OH]D) and calcium supplementation increase birth weight. However, limitations of many trials were highlighted in the reviews. Our aim was to combine genetic and RCT data to estimate causal effects of these two maternal traits on offspring birth weight. METHODS AND FINDINGS: We performed two-sample mendelian randomisation (MR) using genetic instrumental variables associated with 25(OH)D and calcium that had been identified in genome-wide association studies (GWAS; sample 1; N = 122,123 for 25[OH]D and N = 61,275 for calcium). Associations between these maternal genetic variants and offspring birth weight were calculated in the UK Biobank (UKB) (sample 2; N = 190,406). We used data on mother-child pairs from two United Kingdom birth cohorts (combined N = 5,223) in sensitivity analyses to check whether results were influenced by fetal genotype, which is correlated with the maternal genotype (r ≈ 0.5). Further sensitivity analyses to test the reliability of the results included MR-Egger, weighted-median estimator, 'leave-one-out', and multivariable MR analyses. We triangulated MR results with those from RCTs, in which we used randomisation to supplementation with vitamin D (24 RCTs, combined N = 5,276) and calcium (6 RCTs, combined N = 543) as an instrumental variable to determine the effects of 25(OH)D and calcium on birth weight. In the main MR analysis, there was no strong evidence of an effect of maternal 25(OH)D on birth weight (difference in mean birth weight -0.03 g [95% CI -2.48 to 2.42 g, p = 0.981] per 10% higher maternal 25[OH]D). The effect estimate was consistent across our MR sensitivity analyses. Instrumental variable analyses applied to RCTs suggested a weak positive causal effect (5.94 g [95% CI 2.15-9.73, p = 0.002] per 10% higher maternal 25[OH]D), but this result may be exaggerated because of risk of bias in the included RCTs. The main MR analysis for maternal calcium also suggested no strong evidence of an effect on birth weight (-20 g [95% CI -44 to 5 g, p = 0.116] per 1 SD higher maternal calcium level). Some sensitivity analyses suggested that the genetic instrument for calcium was associated with birth weight via exposures that are independent of calcium levels (horizontal pleiotropy). Application of instrumental variable analyses to RCTs suggested that calcium has a substantial effect on birth weight (178 g [95% CI 121-236 g, p = 1.43 × 10-9] per 1 SD higher maternal calcium level) that was not consistent with any of the MR results. However, the RCT instrumental variable estimate may have been exaggerated because of risk of bias in the included RCTs. Other study limitations include the low response rate of UK Biobank, which may bias MR estimates, and the lack of suitable data to test whether the effects of genetic instruments on maternal calcium levels during pregnancy were the same as those outside of pregnancy. CONCLUSIONS: Our results suggest that maternal circulating 25(OH)D does not influence birth weight in otherwise healthy newborns. However, the effect of maternal circulating calcium on birth weight is unclear and requires further exploration with more research including RCT and/or MR analyses with more valid instruments.

Better governance, better access: practising responsible data sharing in the METADAC governance infrastructure.

BACKGROUND: Genomic and biosocial research data about individuals is rapidly proliferating, bringing the potential for novel opportunities for data integration and use. The scale, pace and novelty of these applications raise a number of urgent sociotechnical, ethical and legal questions, including optimal methods of data storage, management and access. Although the open science movement advocates unfettered access to research data, many of the UK's longitudinal cohort studies operate systems of managed data access, in which access is governed by legal and ethical agreements between stewards of research datasets and researchers wishing to make use of them. Amongst other things, these agreements aim to respect the reasonable expectations of the research participants who provided data and samples, as expressed in the consent process. Arguably, responsible data management and governance of data and sample use are foundational to the consent process in longitudinal studies and are an important source of trustworthiness in the eyes of those who contribute data to genomic and biosocial research. METHODS: This paper presents an ethnographic case study exploring the foundational principles of a governance infrastructure for Managing Ethico-social, Technical and Administrative issues in Data ACcess (METADAC), which are operationalised through a committee known as the METADAC Access Committee. METADAC governs access to phenotype, genotype and 'omic' data and samples from five UK longitudinal studies. FINDINGS: Using the example of METADAC, we argue that three key structural features are foundational for practising responsible data sharing: independence and transparency; interdisciplinarity; and participant-centric decision-making. We observe that the international research community is proactively working towards optimising the use of research data, integrating/linking these data with routine data generated by health and social care services and other administrative data services to improve the analysis, interpretation and utility of these data. The governance of these new complex data assemblages will require a range of expertise from across a number of domains and disciplines, including that of study participants. Human-mediated decision-making bodies will be central to ensuring achievable, reasoned and responsible decisions about the use of these data; the METADAC model described in this paper provides an example of how this could be realised.

Genome-wide association study identifies common and low-frequency variants at the AMH gene locus that strongly predict serum AMH levels in males.

Anti-Müllerian hormone (AMH) is an essential messenger of sexual differentiation in the foetus and is an emerging biomarker of postnatal reproductive function in females. Due to a paucity of adequately sized studies, the genetic determinants of circulating AMH levels are poorly characterized. In samples from 2815 adolescents aged 15 from the ALSPAC study, we performed the first genome-wide association study of serum AMH levels across a set of ∼9 m '1000 Genomes Reference Panel' imputed genetic variants. Genetic variants at the AMH protein-coding gene showed considerable allelic heterogeneity, with both common variants [rs4807216 (P(Male) = 2 × 10(-49), Beta: ∼0.9 SDs per allele), rs8112524 (P(Male) = 3 × 10(-8), Beta: ∼0.25)] and low-frequency variants [rs2385821 (P(Male) = 6 × 10(-31), Beta: ∼1.2, frequency 3.6%)] independently associated with apparently large effect sizes in males, but not females. For all three SNPs, we highlight mechanistic links to AMH gene function and demonstrate highly significant sex interactions (P(Het) 0.0003-6.3 × 10(-12)), culminating in contrasting estimates of trait variance explained (24.5% in males versus 0.8% in females). Using these SNPs as a genetic proxy for AMH levels, we found no evidence in additional datasets to support a biological role for AMH in complex traits and diseases in men.

Collapsed methylation quantitative trait loci analysis for low frequency and rare variants.

BACKGROUND: Single variant approaches have been successful in identifying DNA methylation quantitative trait loci (mQTL), although as with complex traits they lack the statistical power to identify the effects from rare genetic variants. We have undertaken extensive analyses to identify regions of low frequency and rare variants that are associated with DNA methylation levels. METHODS: We used repeated measurements of DNA methylation from five different life stages in human blood, taken from the Avon Longitudinal Study of Parents and Children (ALSPAC) cohort. Variants were collapsed across CpG islands and their flanking regions to identify variants collectively associated with methylation, where no single variant was individually responsible for the observed signal. All analyses were undertaken using the sequence kernel association test. RESULTS: For loci where no individual variant mQTL was observed based on a single variant analysis, we identified 95 unique regions where the combined effect of low frequency variants (MAF ≤ 5%) provided strong evidence of association with methylation. For loci where there was previous evidence of an individual variant mQTL, a further 3 regions provided evidence of association between multiple low frequency variants and methylation levels. Effects were observed consistently across 5 different time points in the lifecourse and evidence of replication in the TwinsUK and Exeter cohorts was also identified. CONCLUSION: We have demonstrated the potential of this novel approach to mQTL analysis by analysing the combined effect of multiple low frequency or rare variants. Future studies should benefit from applying this approach as a complementary follow up to single variant analyses.

Prenatal and early life influences on epigenetic age in children: a study of mother-offspring pairs from two cohort studies.

DNA methylation-based biomarkers of aging are highly correlated with actual age. Departures of methylation-estimated age from actual age can be used to define epigenetic measures of child development or age acceleration (AA) in adults. Very little is known about genetic or environmental determinants of these epigenetic measures of aging. We obtained DNA methylation profiles using Infinium HumanMethylation450 BeadChips across five time-points in 1018 mother-child pairs from the Avon Longitudinal Study of Parents and Children. Using the Horvath age estimation method, we calculated epigenetic age for these samples. AA was defined as the residuals from regressing epigenetic age on actual age. AA was tested for associations with cross-sectional clinical variables in children. We identified associations between AA and sex, birth weight, birth by caesarean section and several maternal characteristics in pregnancy, namely smoking, weight, BMI, selenium and cholesterol level. Offspring of non-drinkers had higher AA on average but this difference appeared to resolve during childhood. The associations between sex, birth weight and AA found in ARIES were replicated in an independent cohort (GOYA). In children, epigenetic AA measures are associated with several clinically relevant variables, and early life exposures appear to be associated with changes in AA during adolescence. Further research into epigenetic aging, including the use of causal inference methods, is required to better our understanding of aging.

Genome-wide association study identifies SNPs in the MHC class II loci that are associated with self-reported history of whooping cough.

Whooping cough is currently seeing resurgence in countries despite high vaccine coverage. There is considerable variation in subject-specific response to infection and vaccine efficacy, but little is known about the role of human genetics. We carried out a case-control genome-wide association study of adult or parent-reported history of whooping cough in two cohorts from the UK: the ALSPAC cohort and the 1958 British Birth Cohort (815/758 cases and 6341/4308 controls, respectively). We also imputed HLA alleles using dense SNP data in the MHC region and carried out gene-based and gene-set tests of association and estimated the amount of additive genetic variation explained by common SNPs. We observed a novel association at SNPs in the MHC class II region in both cohorts [lead SNP rs9271768 after meta-analysis, odds ratio [95% confidence intervals (CIs)] 1.47 (1.35, 1.6), P-value 1.21E - 18]. Multiple strong associations were also observed at alleles at the HLA class II loci. The majority of these associations were explained by the lead SNP rs9271768. Gene-based and gene-set tests and estimates of explainable common genetic variation could not establish the presence of additional associations in our sample. Genetic variation at the MHC class II region plays a role in susceptibility to whooping cough. These findings provide additional perspective on mechanisms of whooping cough infection and vaccine efficacy.

Longitudinal analysis of DNA methylation associated with birth weight and gestational age.

Gestational age (GA) and birth weight have been implicated in the determination of long-term health. It has been hypothesized that changes in DNA methylation may mediate these long-term effects. We obtained DNA methylation profiles from cord blood and peripheral blood at ages 7 and 17 in the same children from the Avon Longitudinal Study of Parents and Children. Repeated-measures data were used to investigate changes in birth-related methylation during childhood and adolescence. Ten developmental phenotypes (e.g. height) were analysed to identify possible mediation of health effects by DNA methylation. In cord blood, methylation at 224 CpG sites was found to be associated with GA and 23 CpG sites with birth weight. Methylation changed in the majority of these sites over time, but neither birth characteristic was strongly associated with methylation at age 7 or 17 (using a conservative correction for multiple testing of P < 1.03 × 10(-7)), suggesting resolution of differential methylation by early childhood. Associations were observed between birth weight-associated CpG sites and phenotypic characteristics in childhood. One strong association involved birth weight, methylation of a CpG site proximal to the NFIX locus and bone mineral density at age 17. Analysis of serial methylation from birth to adolescence provided evidence for a lack of persistence of methylation differences beyond early childhood. Sites associated with birth weight were linked to developmental genes and have methylation levels which are associated with developmental phenotypes. Replication and interrogation of causal relationships are needed to substantiate whether methylation differences at birth influence the association between birth weight and development.

Prenatal exposure to maternal smoking and offspring DNA methylation across the lifecourse: findings from the Avon Longitudinal Study of Parents and Children (ALSPAC).

Maternal smoking during pregnancy has been found to influence newborn DNA methylation in genes involved in fundamental developmental processes. It is pertinent to understand the degree to which the offspring methylome is sensitive to the intensity and duration of prenatal smoking. An investigation of the persistence of offspring methylation associated with maternal smoking and the relative roles of the intrauterine and postnatal environment is also warranted. In the Avon Longitudinal Study of Parents and Children, we investigated associations between prenatal exposure to maternal smoking and offspring DNA methylation at multiple time points in approximately 800 mother-offspring pairs. In cord blood, methylation at 15 CpG sites in seven gene regions (AHRR, MYO1G, GFI1, CYP1A1, CNTNAP2, KLF13 and ATP9A) was associated with maternal smoking, and a dose-dependent response was observed in relation to smoking duration and intensity. Longitudinal analysis of blood DNA methylation in serial samples at birth, age 7 and 17 years demonstrated that some CpG sites showed reversibility of methylation (GFI1, KLF13 and ATP9A), whereas others showed persistently perturbed patterns (AHRR, MYO1G, CYP1A1 and CNTNAP2). Of those showing persistence, we explored the effect of postnatal smoke exposure and found that the major contribution to altered methylation was attributed to a critical window of in utero exposure. A comparison of paternal and maternal smoking and offspring methylation showed consistently stronger maternal associations, providing further evidence for causal intrauterine mechanisms. These findings emphasize the sensitivity of the methylome to maternal smoking during early development and the long-term impact of such exposure.

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.

TRPA1 gene polymorphisms and childhood asthma.

BACKGROUND: Animal data have suggested that the transient receptor potential ankyrin-1 (TRPA1) ion channel plays a key role in promoting airway inflammation in asthma and may mediate effects of paracetamol on asthma, yet confirmatory human data are lacking. To study associations of TRPA1 gene variants with childhood asthma and total IgE concentration, and interactions between TRPA1 and prenatal paracetamol exposure on these outcomes. METHODS: We analysed associations between 31 TRPA1 single nucleotide polymorphisms (SNPs) and current doctor-diagnosed asthma and total IgE concentration at 7.5 years in the Avon Longitudinal Study of Parents and Children (ALSPAC) birth cohort. We sought to confirm the most significant associations with comparable outcomes in the Prevention and Incidence of Asthma and Mite Allergy (PIAMA) and Generation R birth cohorts. In ALSPAC, we explored interactions with prenatal paracetamol exposure. RESULTS: In ALSPAC, there was strong evidence for association between six SNPs and asthma: rs959974 and rs1384001 (per-allele odds ratio for both: 1.30 (95% CI: 1.15-1.47), p = 0.00001), rs7010969 (OR 1.28 (1.13-1.46), p = 0.00004), rs3735945 (OR 1.30 (1.09-1.55), p = 0.003), rs920829 (OR 1.30 (1.09-1.54), p = 0.004) and rs4738202 (OR 1.22 (1.07-1.39), p = 0.004). In a meta-analysis across the three cohorts, the pooled effect estimates confirmed that all six SNPs were significantly associated with asthma. In ALSPAC, TRPA1 associations with asthma were not modified by prenatal paracetamol, although associations with IgE concentration were. CONCLUSION: This study suggests that TRPA1 may play a role in the development of childhood asthma. (249 words).

Age at menarche and lung function: a Mendelian randomization study.

A trend towards earlier menarche in women has been associated with childhood factors (e.g. obesity) and hypothesised environmental exposures (e.g. endocrine disruptors present in household products). Observational evidence has shown detrimental effects of early menarche on various health outcomes including adult lung function, but these might represent spurious associations due to confounding. To address this we used Mendelian randomization where genetic variants are used as proxies for age at menarche, since genetic associations are not affected by classical confounding. We estimated the effects of age at menarche on forced vital capacity (FVC), a proxy for restrictive lung impairment, and ratio of forced expiratory volume in one second to FVC (FEV1/FVC), a measure of airway obstruction, in both adulthood and adolescence. We derived SNP-age at menarche association estimates for 122 variants from a published genome-wide meta-analysis (N = 182,416), with SNP-lung function estimates obtained by meta-analysing three studies of adult women (N = 46,944) and two of adolescent girls (N = 3025). We investigated the impact of departures from the assumption of no pleiotropy through sensitivity analyses. In adult women, in line with previous evidence, we found an effect on restrictive lung impairment with a 24.8 mL increase in FVC per year increase in age at menarche (95% CI 1.8-47.9; p = 0.035); evidence was stronger after excluding potential pleiotropic variants (43.6 mL; 17.2-69.9; p = 0.001). In adolescent girls we found an opposite effect (-56.5 mL; -108.3 to -4.7; p = 0.033), suggesting that the detrimental effect in adulthood may be preceded by a short-term post-pubertal benefit. Our secondary analyses showing results in the same direction in men and boys, in whom age at menarche SNPs have also shown association with sexual development, suggest a role for pubertal timing in general rather than menarche specifically. We found no effect on airway obstruction (FEV1/FVC).

The regulatory element READ1 epistatically influences reading and language, with both deleterious and protective alleles.

BACKGROUND: Reading disability (RD) and language impairment (LI) are heritable learning disabilities that obstruct acquisition and use of written and spoken language, respectively. We previously reported that two risk haplotypes, each in strong linkage disequilibrium (LD) with an allele of READ1, a polymorphic compound short tandem repeat within intron 2 of risk gene DCDC2, are associated with RD and LI. Additionally, we showed a non-additive genetic interaction between READ1 and KIAHap, a previously reported risk haplotype in risk gene KIAA0319, and that READ1 binds the transcriptional regulator ETV6. OBJECTIVE: To examine the hypothesis that READ1 is a transcriptional regulator of KIAA0319. METHODS: We characterised associations between READ1 alleles and RD and LI in a large European cohort, and also assessed interactions between READ1 and KIAHap and their effect on performance on measures of reading, language and IQ. We also used family-based data to characterise the genetic interaction, and chromatin conformation capture (3C) to investigate the possibility of a physical interaction between READ1 and KIAHap. RESULTS AND CONCLUSIONS: READ1 and KIAHap show interdependence--READ1 risk alleles synergise with KIAHap, whereas READ1 protective alleles act epistatically to negate the effects of KIAHap. The family data suggest that these variants interact in trans genetically, while the 3C results show that a region of DCDC2 containing READ1 interacts physically with the region upstream of KIAA0319. These data support a model in which READ1 regulates KIAA0319 expression through KIAHap and in which the additive effects of READ1 and KIAHap alleles are responsible for the trans genetic interaction.

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 identifies common variants at the SERPINA6/SERPINA1 locus influencing plasma cortisol and corticosteroid binding globulin.

Variation in plasma levels of cortisol, an essential hormone in the stress response, is associated in population-based studies with cardio-metabolic, inflammatory and neuro-cognitive traits and diseases. Heritability of plasma cortisol is estimated at 30-60% but no common genetic contribution has been identified. The CORtisol NETwork (CORNET) consortium undertook genome wide association meta-analysis for plasma cortisol in 12,597 Caucasian participants, replicated in 2,795 participants. The results indicate that <1% of variance in plasma cortisol is accounted for by genetic variation in a single region of chromosome 14. This locus spans SERPINA6, encoding corticosteroid binding globulin (CBG, the major cortisol-binding protein in plasma), and SERPINA1, encoding α1-antitrypsin (which inhibits cleavage of the reactive centre loop that releases cortisol from CBG). Three partially independent signals were identified within the region, represented by common SNPs; detailed biochemical investigation in a nested sub-cohort showed all these SNPs were associated with variation in total cortisol binding activity in plasma, but some variants influenced total CBG concentrations while the top hit (rs12589136) influenced the immunoreactivity of the reactive centre loop of CBG. Exome chip and 1000 Genomes imputation analysis of this locus in the CROATIA-Korcula cohort identified missense mutations in SERPINA6 and SERPINA1 that did not account for the effects of common variants. These findings reveal a novel common genetic source of variation in binding of cortisol by CBG, and reinforce the key role of CBG in determining plasma cortisol levels. In turn this genetic variation may contribute to cortisol-associated degenerative diseases.

Alleles of a polymorphic ETV6 binding site in DCDC2 confer risk of reading and language impairment.

Reading disability (RD) and language impairment (LI) are common learning disabilities that make acquisition and utilization of reading and verbal language skills, respectively, difficult for affected individuals. Both disorders have a substantial genetic component with complex inheritance. Despite decades of study, reading and language, like many other complex traits, consistently evade identification of causative and functional variants. We previously identified a putative functional risk variant, named BV677278 for its GenBank accession number, for RD in DCDC2. This variant consists of an intronic microdeletion and a highly polymorphic short tandem repeat (STR) within its breakpoints. We have also shown this STR to bind to an unknown nuclear protein with high specificity. Here, we replicate BV677278's association with RD, expand its association to LI, identify the BV677278-binding protein as the transcription factor ETV6, and provide compelling genetic evidence that BV677278 is a regulatory element that influences reading and language skills. We also provide evidence that BV677278 interacts nonadditively with KIAA0319, an RD-associated gene, to adversely affect several reading and cognitive phenotypes. On the basis of these data, we propose a new name for BV677278: "READ1" or "regulatory element associated with dyslexia 1."

Common variants in left/right asymmetry genes and pathways are associated with relative hand skill.

Humans display structural and functional asymmetries in brain organization, strikingly with respect to language and handedness. The molecular basis of these asymmetries is unknown. We report a genome-wide association study meta-analysis for a quantitative measure of relative hand skill in individuals with dyslexia [reading disability (RD)] (n = 728). The most strongly associated variant, rs7182874 (P = 8.68 × 10(-9)), is located in PCSK6, further supporting an association we previously reported. We also confirmed the specificity of this association in individuals with RD; the same locus was not associated with relative hand skill in a general population cohort (n = 2,666). As PCSK6 is known to regulate NODAL in the development of left/right (LR) asymmetry in mice, we developed a novel approach to GWAS pathway analysis, using gene-set enrichment to test for an over-representation of highly associated variants within the orthologs of genes whose disruption in mice yields LR asymmetry phenotypes. Four out of 15 LR asymmetry phenotypes showed an over-representation (FDR ≤ 5%). We replicated three of these phenotypes; situs inversus, heterotaxia, and double outlet right ventricle, in the general population cohort (FDR ≤ 5%). Our findings lead us to propose that handedness is a polygenic trait controlled in part by the molecular mechanisms that establish LR body asymmetry early in development.