Nicole M. Warrington

Sequence variants affecting eosinophil numbers associate with asthma and myocardial infarction

Eosinophils are pleiotropic multifunctional leukocytes involved in initiation and propagation of inflammatory responses and thus have important roles in the pathogenesis of inflammatory diseases. Here we describe a genome-wide association scan for sequence variants affecting eosinophil counts in blood of 9,392 Icelanders. The most significant SNPs were studied further in 12,118 Europeans and 5,212 East Asians. SNPs at 2q12 (rs1420101), 2q13 (rs12619285), 3q21 (rs4857855), 5q31 (rs4143832) and 12q24 (rs3184504) reached genome-wide significance (P = 5.3 × 10−14, 5.4 × 10−10, 8.6 × 10−17, 1.2 × 10−10 and 6.5 × 10−19, respectively). A SNP at IL1RL1 associated with asthma (P = 5.5 × 10−12) in a collection of ten different populations (7,996 cases and 44,890 controls). SNPs at WDR36, IL33 and MYB that showed suggestive association with eosinophil counts were also associated with atopic asthma (P = 4.2 × 10−6, 2.2 × 10−5 and 2.4 × 10−4, respectively). We also found that a nonsynonymous SNP at 12q24, in SH2B3, associated significantly (P = 8.6 × 10−8) with myocardial infarction in six different populations (6,650 cases and 40,621 controls).

Meta-analysis of genome-wide association studies identifies ten loci influencing allergic sensitization

Allergen-specific immunoglobulin E (present in allergic sensitization) has a central role in the pathogenesis of allergic disease. We performed the first large-scale genome-wide association study (GWAS) of allergic sensitization in 5,789 affected individuals and 10,056 controls and followed up the top SNP at each of 26 loci in 6,114 affected individuals and 9,920 controls. We increased the number of susceptibility loci with genome-wide significant association with allergic sensitization from three to ten, including SNPs in or near TLR6, C11orf30, STAT6, SLC25A46, HLA-DQB1, IL1RL1, LPP, MYC, IL2 and HLA-B. All the top SNPs were associated with allergic symptoms in an independent study. Risk-associated variants at these ten loci were estimated to account for at least 25% of allergic sensitization and allergic rhinitis. Understanding the molecular mechanisms underlying these associations may provide new insights into the etiology of allergic disease.

Association between common variation at the FTO locus and changes in body mass index from infancy to late childhood: the complex nature of genetic association through growth and development.

An age-dependent association between variation at the FTO locus and BMI in children has been suggested. We meta-analyzed associations between the FTO locus (rs9939609) and BMI in samples, aged from early infancy to 13 years, from 8 cohorts of European ancestry. We found a positive association between additional minor (A) alleles and BMI from 5.5 years onwards, but an inverse association below age 2.5 years. Modelling median BMI curves for each genotype using the LMS method, we found that carriers of minor alleles showed lower BMI in infancy, earlier adiposity rebound (AR), and higher BMI later in childhood. Differences by allele were consistent with two independent processes: earlier AR equivalent to accelerating developmental age by 2.37% (95% CI 1.87, 2.87, p = 10−20) per A allele and a positive age by genotype interaction such that BMI increased faster with age (p = 10−23). We also fitted a linear mixed effects model to relate genotype to the BMI curve inflection points adiposity peak (AP) in infancy and AR. Carriage of two minor alleles at rs9939609 was associated with lower BMI at AP (−0.40% (95% CI: −0.74, −0.06), p = 0.02), higher BMI at AR (0.93% (95% CI: 0.22, 1.64), p = 0.01), and earlier AR (−4.72% (−5.81, −3.63), p = 10−17), supporting cross-sectional results. Overall, we confirm the expected association between variation at rs9939609 and BMI in childhood, but only after an inverse association between the same variant and BMI in infancy. Patterns are consistent with a shift on the developmental scale, which is reflected in association with the timing of AR rather than just a global increase in BMI. Results provide important information about longitudinal gene effects and about the role of FTO in adiposity. The associated shifts in developmental timing have clinical importance with respect to known relationships between AR and both later-life BMI and metabolic disease risk.

Intestinal dysbiosis in ankylosing spondylitis.

Ankylosing spondylitis (AS) is a common, highly heritable immune‐mediated arthropathy that occurs in genetically susceptible individuals exposed to an unknown but likely ubiquitous environmental trigger. There is a close relationship between the gut and spondyloarthritis, as exemplified in patients with reactive arthritis, in whom a typically self‐limiting arthropathy follows either a gastrointestinal or urogenital infection. Microbial involvement in AS has been suggested; however, no definitive link has been established. The aim of this study was to determine whether the gut in patients with AS carries a distinct microbial signature compared with that in the gut of healthy control subjects.

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.uk Supplementary information: Supplementary data are available at Bioinformatics online.

Brief Report: Intestinal Dysbiosis in Ankylosing Spondylitis

Ankylosing spondylitis (AS) is a common, highly heritable immune‐mediated arthropathy that occurs in genetically susceptible individuals exposed to an unknown but likely ubiquitous environmental trigger. There is a close relationship between the gut and spondyloarthritis, as exemplified in patients with reactive arthritis, in whom a typically self‐limiting arthropathy follows either a gastrointestinal or urogenital infection. Microbial involvement in AS has been suggested; however, no definitive link has been established. The aim of this study was to determine whether the gut in patients with AS carries a distinct microbial signature compared with that in the gut of healthy control subjects.

The Association of C-Reactive Protein and CRP Genotype with Coronary Heart Disease: Findings from Five Studies with 4,610 Cases amongst 18,637 Participants

Background It is unclear whether C-reactive protein (CRP) is causally related to coronary heart disease (CHD). Genetic variants that are known to be associated with CRP levels can be used to provide causal inference of the effect of CRP on CHD. Our objective was to examine the association between CRP genetic variant +1444C>T (rs1130864) and CHD risk in the largest study to date of this association. Methods and Results We estimated the association of CRP genetic variant +1444C>T (rs1130864) with CRP levels and with CHD in five studies and then pooled these analyses (N = 18,637 participants amongst whom there were 4,610 cases). CRP was associated with potential confounding factors (socioeconomic position, physical activity, smoking and body mass) whereas genotype (rs1130864) was not associated with these confounders. The pooled odds ratio of CHD per doubling of circulating CRP level after adjustment for age and sex was 1.13 (95%CI: 1.06, 1.21), and after further adjustment for confounding factors it was 1.07 (95%CI: 1.02, 1.13). Genotype (rs1130864) was associated with circulating CRP; the pooled ratio of geometric means of CRP level among individuals with the TT genotype compared to those with the CT/CC genotype was 1.21 (95%CI: 1.15, 1.28) and the pooled ratio of geometric means of CRP level per additional T allele was 1.14 (95%CI: 1.11, 1.18), with no strong evidence in either analyses of between study heterogeneity (I2 = 0%, p>0.9 for both analyses). There was no association of genotype (rs1130864) with CHD: pooled odds ratio 1.01 (95%CI: 0.88, 1.16) comparing individuals with TT genotype to those with CT/CC genotype and 0.96 (95%CI: 0.90, 1.03) per additional T allele (I2<7.5%, p>0.6 for both meta-analyses). An instrumental variables analysis (in which the proportion of CRP levels explained by rs1130864 was related to CHD) suggested that circulating CRP was not associated with CHD: the odds ratio for a doubling of CRP level was 1.04 (95%CI: 0.61, 1.80). Conclusions We found no association of a genetic variant, which is known to be related to CRP levels, (rs1130864) and having CHD. These findings do not support a causal association between circulating CRP and CHD risk, but very large, extended, genetic association studies would be required to rule this out.

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.

Association of Genetic Loci With Glucose Levels in Childhood and Adolescence: A Meta-Analysis of Over 6,000 Children

OBJECTIVE To investigate whether associations of common genetic variants recently identified for fasting glucose or insulin levels in nondiabetic adults are detectable in healthy children and adolescents. RESEARCH DESIGN AND METHODS A total of 16 single nucleotide polymorphisms (SNPs) associated with fasting glucose were genotyped in six studies of children and adolescents of European origin, including over 6,000 boys and girls aged 9–16 years. We performed meta-analyses to test associations of individual SNPs and a weighted risk score of the 16 loci with fasting glucose. RESULTS Nine loci were associated with glucose levels in healthy children and adolescents, with four of these associations reported in previous studies and five reported here for the first time (GLIS3, PROX1, SLC2A2, ADCY5, and CRY2). Effect sizes were similar to those in adults, suggesting age-independent effects of these fasting glucose loci. Children and adolescents carrying glucose-raising alleles of G6PC2, MTNR1B, GCK, and GLIS3 also showed reduced β-cell function, as indicated by homeostasis model assessment of β-cell function. Analysis using a weighted risk score showed an increase [β (95% CI)] in fasting glucose level of 0.026 mmol/L (0.021–0.031) for each unit increase in the score. CONCLUSIONS Novel fasting glucose loci identified in genome-wide association studies of adults are associated with altered fasting glucose levels in healthy children and adolescents with effect sizes comparable to adults. In nondiabetic adults, fasting glucose changes little over time, and our results suggest that age-independent effects of fasting glucose loci contribute to long-term interindividual differences in glucose levels from childhood onwards.

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 (re = 0.55) was higher than the residual correlation between SK-BMD and BMD at appendicular sites (re = 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 Pagets 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.