A single nucleotide polymorphism genetic risk score to aid diagnosis of coeliac disease: a pilot study in clinical care

BACKGROUND

Single nucleotide polymorphism-based genetic risk scores (GRS) model genetic risk as a continuum and can discriminate coeliac disease but have not been validated in clinic. Human leukocyte antigen (HLA) DQ gene testing is available in clinic but does not include non-HLA attributed risk and is limited by discrete risk stratification.

AIMS

To accurately characterise both HLA and non-HLA coeliac disease genetic risk as a single nucleotide polymorphism-based GRS and evaluate diagnostic utility.

METHODS

We developed a 42 single nucleotide polymorphism coeliac disease GRS from a European case-control study (12 041 cases vs 12 228 controls) using HLA-DQ imputation and published genome-wide association studies. We validated the GRS in UK Biobank (1237 cases) and developed direct genotyping assays. We tested the coeliac disease GRS in a pilot clinical cohort of 128 children presenting with suspected coeliac disease.

RESULTS

The GRS was more discriminative of coeliac disease than HLA-DQ stratification in UK Biobank (receiver operating characteristic area under the curve [ROC-AUC] = 0.88 [95% CIs: 0.87-0.89] vs 0.82 [95% CIs: 0.80-0.83]). We demonstrated similar discrimination in the pilot clinical cohort (114 cases vs 40 controls, ROC-AUC = 0.84 [95% CIs: 0.76-0.91]). As a rule-out test, no children with coeliac disease in the clinical cohort had a GRS below 38th population centile.

CONCLUSIONS

A single nucleotide polymorphism-based GRS may offer more effective and cost-efficient testing of coeliac disease genetic risk in comparison to HLA-DQ stratification. As a comparatively inexpensive test it could facilitate non-invasive coeliac disease diagnosis but needs detailed assessment in the context of other diagnostic tests and against current diagnostic algorithms.

Large Copy-Number Variants in UK Biobank Caused by Clonal Hematopoiesis May Confound Penetrance Estimates

Large copy-number variants (CNVs) are strongly associated with both developmental delay and cancer, but the type of disease depends strongly on when and where the mutation occurred, i.e., germline versus somatic. We used microarray data from UK Biobank to investigate the prevalence and penetrance of large autosomal CNVs and chromosomal aneuploidies using a standard CNV detection algorithm not designed for detecting mosaic variants. We found 160 individuals that carry >10 Mb copy number changes, including 56 with whole chromosome aneuploidies. Nineteen (12%) individuals had a diagnosis of Down syndrome or other developmental disorder, while 84 (52.5%) individuals had a diagnosis of hematological malignancies or chronic myeloproliferative disorders. Notably, there was no evidence of mosaicism in the blood for many of these large CNVs, so they could easily be mistaken for germline alleles even when caused by somatic mutations. We therefore suggest that somatic mutations associated with blood cancers may result in false estimates of rare variant penetrance from population biobanks.

Refinement of the critical genomic region for congenital hyperinsulinism in the Chromosome 9p deletion syndrome

Background: Large contiguous gene deletions at the distal end of the short arm of chromosome 9 result in the complex multi-organ condition chromosome 9p deletion syndrome.  A range of clinical features can result from these deletions with the most common being facial dysmorphisms and neurological impairment. Congenital hyperinsulinism is a rarely reported feature of the syndrome with the genetic mechanism for the dysregulated insulin secretion being unknown.  Methods: We studied the clinical and genetic characteristics of 12 individuals with chromosome 9p deletions who had a history of neonatal hypoglycaemia. Using off-target reads generated from targeted next-generation sequencing of the genes known to cause hyperinsulinaemic hypoglycaemia (n=9), or microarray analysis (n=3), we mapped the minimal shared deleted region on chromosome 9 in this cohort. Targeted sequencing was performed in three patients to search for a recessive mutation unmasked by the deletion. Results: In 10/12 patients with hypoglycaemia, hyperinsulinism was confirmed biochemically. A range of extra-pancreatic features were also reported in these patients consistent with the diagnosis of the Chromosome 9p deletion syndrome. The minimal deleted region was mapped to 7.2 Mb, encompassing 38 protein-coding genes. In silico analysis of these genes highlighted SMARCA2 and RFX3 as potential candidates for the hypoglycaemia. Targeted sequencing performed on three of the patients did not identify a second disease-causing variant within the minimal deleted region. Conclusions: This study identifies 9p deletions as an important cause of hyperinsulinaemic hypoglycaemia and increases the number of cases reported with 9p deletions and hypoglycaemia to 15 making this a more common feature of the syndrome than previously appreciated.  Whilst the precise genetic mechanism of the dysregulated insulin secretion could not be determined in these patients, mapping the deletion breakpoints highlighted potential candidate genes for hypoglycaemia within the deleted region.

A combined risk score enhances prediction of type 1 diabetes among susceptible children

Type 1 diabetes (T1D)-an autoimmune disease that destroys the pancreatic islets, resulting in insulin deficiency-often begins early in life when islet autoantibody appearance signals high risk1. However, clinical diabetes can follow in weeks or only after decades, and is very difficult to predict. Ketoacidosis at onset remains common2,3 and is most severe in the very young4,5, in whom it can be life threatening and difficult to treat6-9. Autoantibody surveillance programs effectively prevent most ketoacidosis10-12 but require frequent evaluations whose expense limits public health adoption13. Prevention therapies applied before onset, when greater islet mass remains, have rarely been feasible14 because individuals at greatest risk of impending T1D are difficult to identify. To remedy this, we sought accurate, cost-effective estimation of future T1D risk by developing a combined risk score incorporating both fixed and variable factors (genetic, clinical and immunological) in 7,798 high-risk children followed closely from birth for 9.3 years. Compared with autoantibodies alone, the combined model dramatically improves T1D prediction at ≥2 years of age over horizons up to 8 years of age (area under the receiver operating characteristic curve ≥ 0.9), doubles the estimated efficiency of population-based newborn screening to prevent ketoacidosis, and enables individualized risk estimates for better prevention trial selection.

IgA Nephropathy Genetic Risk Score to Estimate the Prevalence of IgA Nephropathy in UK Biobank

Background

IgA nephropathy (IgAN) is the commonest glomerulonephritis worldwide. Its prevalence is difficult to estimate, as people with mild disease do not commonly receive a biopsy diagnosis. We aimed to generate an IgA nephropathy genetic risk score (IgAN-GRS) and estimate the proportion of people with hematuria who had IgAN in the UK Biobank (UKBB).

Methods

We calculated an IgAN-GRS using 14 single-nucleotide polymorphisms (SNPs) drawn from the largest European Genome-Wide Association Study (GWAS) and validated the IgAN-GRS in 464 biopsy-proven IgAN European cases from the UK Glomerulonephritis DNA Bank (UKGDB) and in 379,767 Europeans in the UKBB. We used the mean of IgAN-GRS to calculate the proportion of potential IgAN in 14,181 with hematuria and other nonspecific renal phenotypes from 379,767 Europeans in the UKBB.

Results

The IgAN-GRS was higher in the IgAN cohort (4.30; 95% confidence interval [95% CI: 4.23–4.38) than in controls (3.98; 3.97–3.98; P < 0.0001). The mean GRS in UKBB participants with hematuria (n = 12,858) was higher (4.04; 4.02–4.06) than UKBB controls (3.98; 3.97–3.98; P < 0.0001) and higher in those with hematuria, hypertension, and microalbuminuria (n = 1323) (4.07; 4.02–4.13) versus (3.98; 3.97–3.98; P = 0.0003). Using the difference in these means, we estimated that IgAN accounted for 19% of noncancer hematuria and 28% with hematuria, hypertension, and microalbuminuria in UKBB.

Conclusions

We used an IgAN-GRS to estimate the prevalence of IgAN contributing to common phenotypes that are not always biopsied. The noninvasive use of polygenic risk in this setting may have further utility to identify likely etiology of nonspecific renal phenotypes in large population cohorts.

A genome-wide association study implicates multiple mechanisms influencing raised urinary albumin-creatinine ratio

Raised albumin–creatinine ratio (ACR) is an indicator of microvascular damage and renal disease. We aimed to identify genetic variants associated with raised ACR and study the implications of carrying multiple ACR-raising alleles with metabolic and vascular-related disease. We performed a genome-wide association study of ACR using 437 027 individuals from the UK Biobank in the discovery phase, 54 527 more than previous studies, and followed up our findings in independent studies. We identified 62 independent associations with ACR across 56 loci (P < 5 × 10^–8), of which 20 were not previously reported. Pathway analyses and the identification of 20 of the 62 variants (at r² > 0.8) coinciding with signals for at least 16 related metabolic and vascular traits, suggested multiple pathways leading to raised ACR levels. After excluding variants at the CUBN locus, known to alter ACR via effects on renal absorption, an ACR genetic risk score was associated with a higher risk of hypertension, and less strongly, type 2 diabetes and stroke. For some rare genotype combinations at the CUBN locus, most individuals had ACR levels above the microalbuminuria clinical threshold. Contrary to our hypothesis, individuals carrying more CUBN ACR-raising alleles, and above the clinical threshold, had a higher frequency of vascular disease. The CUBN allele effects on ACR were twice as strong in people with diabetes—a result robust to an optimization-algorithm approach to simulating interactions, validating previously reported gene–diabetes interactions (P ≤ 4 × 10^–5). In conclusion, a variety of genetic mechanisms and traits contribute to variation in ACR.

Type 1 diabetes genetic risk score is discriminative of diabetes in non-Europeans: evidence from a study in India

Type 1 diabetes (T1D) is a significant problem in Indians and misclassification of T1D and type 2 diabetes (T2D) is a particular problem in young adults in this population due to the high prevalence of early onset T2D at lower BMI. We have previously shown a genetic risk score (GRS) can be used to discriminate T1D from T2D in Europeans. We aimed to test the ability of a T1D GRS to discriminate T1D from T2D and controls in Indians. We studied subjects from Pune, India of Indo-European ancestry; T1D (n = 262 clinically defined, 200 autoantibody positive), T2D (n = 345) and controls (n = 324). We used the 9 SNP T1D GRS generated in Europeans and assessed its ability to discriminate T1D from T2D and controls in Indians. We compared Indians with Europeans from the Wellcome Trust Case Control Consortium study; T1D (n = 1963), T2D (n = 1924) and controls (n = 2938). The T1D GRS was discriminative of T1D from T2D in Indians but slightly less than in Europeans (ROC AUC 0.84 v 0.87, p < 0.0001). HLA SNPs contributed the majority of the discriminative power in Indians. A T1D GRS using SNPs defined in Europeans is discriminative of T1D from T2D and controls in Indians. As with Europeans, the T1D GRS may be useful for classifying diabetes in Indians.

The acute transcriptional response to resistance exercise: impact of age and contraction mode

Optimization of resistance exercise (RE) remains a hotbed of research for muscle building and maintenance. However, the interactions between the contractile components of RE (i.e. concentric (CON) and eccentric (ECC)) and age, are poorly defined. We used transcriptomics to compare age-related molecular responses to acute CON and ECC exercise. Eight young (21±1 y) and eight older (70±1 y) exercise-naïve male volunteers had vastus lateralis biopsies collected at baseline and 5 h post unilateral CON and contralateral ECC exercise. RNA was subjected to next-generation sequencing and differentially expressed (DE) genes tested for pathway enrichment using Gene Ontology (GO). The young transcriptional response to CON and ECC was highly similar and older adults displayed moderate contraction-specific profiles, with no GO enrichment. Age-specific responses to ECC revealed 104 DE genes unique to young, and 170 DE genes in older muscle, with no GO enrichment. Following CON, 15 DE genes were young muscle-specific, whereas older muscle uniquely expressed 147 up-regulated genes enriched for cell adhesion and blood vessel development, and 28 down-regulated genes involved in mitochondrial respiration, amino acid and lipid metabolism. Thus, older age is associated with contraction-specific regulation often without clear functional relevance, perhaps reflecting a degree of stochastic age-related dysregulation.

Comparison of Genetic Liability for Sleep Traits Among Individuals With Bipolar Disorder I or II and Control Participants

IMPORTANCE

Insomnia, hypersomnia, and an evening chronotype are common in individuals with bipolar disorder (BD), but whether this reflects shared genetic liability is unclear. Stratifying by BD subtypes could elucidate this association and inform sleep and BD research.

OBJECTIVE

To assess whether polygenic risk scores (PRSs) for sleep traits are associated with BD subtypes I and II.

DESIGN, SETTING, AND PARTICIPANTS

This case-control study was conducted in the United Kingdom and Sweden with participants with BD and control participants. Multinomial regression was used to assess whether PRSs for insomnia, daytime sleepiness, sleep duration, and chronotype are associated with BD subtypes compared with control participants. Affected individuals were recruited from the Bipolar Disorder Research Network. Control participants were recruited from the 1958 British Birth Cohort and the UK Blood Service. Analyses were repeated in an independent Swedish sample from August 2018 to July 2019. All participants were of European ancestry.

EXPOSURES

Standardized PRSs derived using alleles from genome-wide association studies of insomnia, sleep duration, daytime sleepiness, and chronotype. These were adjusted for the first 10 population principal components, genotyping platforms, and sex.

MAIN OUTCOMES AND MEASURES

Association of PRSs with BD subtypes, determined by semistructured psychiatric interview and case notes.

RESULTS

The main analysis included 4672 participants with BD (3132 female participants [67.0%]; 3404 with BD-I [72.9%]) and 5714 control participants (2812 female participants [49.2%]). Insomnia PRS was associated with increased risk of BD-II (relative risk [RR], 1.14 [95% CI, 1.07-1.21]; P = 8.26 × 10-5) but not BD-I (RR, 0.98 [95% CI, 0.94-1.03]; P = .409) relative to control participants. Sleep-duration PRS was associated with BD-I (RR, 1.10 [95% CI, 1.06-1.15]; P = 1.13 × 10-5) but not BD-II (RR, 0.99 [95% CI, 0.93-1.06]; P = .818). Associations between (1) insomnia PRS and BD-II and (2) sleep-duration PRS and BD-I were replicated in the Swedish sample of 4366 individuals with BD (2697 female participants [61.8%]; 2627 with BD-I [60.2%]) and 6091 control participants (3767 female participants [61.8%]). Chronotype and daytime-sleepiness PRS were not associated with BD subtypes.

CONCLUSIONS AND RELEVANCE

Per this analysis, BD subtypes differ in genetic liability to insomnia and hypersomnia, providing further evidence that the distinction between BD-I and BD-II has genetic validity. This distinction will be crucial in selecting participants for future research on the role of sleep disturbance in BD.

Assessment of MTNR1B Type 2 Diabetes Genetic Risk Modification by Shift Work and Morningness-Eveningness Preference in the UK Biobank

Night shift work, behavioral rhythms, and the common MTNR1B risk single nucleotide polymorphism (SNP), rs10830963, associate with type 2 diabetes; however, whether they exert joint effects to exacerbate type 2 diabetes risk is unknown. Among employed participants of European ancestry in the UK Biobank (N = 189,488), we aimed to test the cross-sectional independent associations and joint interaction effects of these risk factors on odds of type 2 diabetes (n = 5,042 cases) and HbA_1c levels (n = 175,156). Current shift work, definite morning or evening preference, and MTNR1B rs10830963 risk allele associated with type 2 diabetes and HbA_1c levels. The effect of rs10830963 was not modified by shift work schedules. While marginal evidence of interaction between self-reported morningness-eveningness preference and rs10830963 on risk of type 2 diabetes was seen, this interaction did not persist when analysis was expanded to include all participants regardless of employment status and when accelerometer-derived sleep midpoint was used as an objective measure of morningness-eveningness preference. Our findings suggest that MTNR1B risk allele carriers who carry out shift work or have more extreme morningness-eveningness preference may not have enhanced risk of type 2 diabetes.

Genome-Wide Association Study of Microscopic Colitis in the UK Biobank Confirms Immune-Related Pathogenesis

BACKGROUND AND AIMS

The causes of microscopic colitis are currently poorly understood. Previous reports have found clinical associations with coeliac disease and genetic associations at the human leukocyte antigen [HLA] locus on the ancestral 8.1 haplotype. We investigated pharmacological and genetic factors associated with microscopic colitis in the UK Biobank.

METHODS

In total, 483 European UK Biobank participants were identified by ICD10 coding, and a genome-wide association study was performed using BOLT-LMM, with a sensitivity analysis performed excluding potential confounders. The HLA*IMP:02 algorithm was used to estimate allele frequency at 11 classical HLA genes, and downstream analysis was performed using FUMA. Genetic overlap with inflammatory bowel disease [Crohn's disease and ulcerative colitis] was investigated using genetic risk scores.

RESULTS

We found significant phenotypic associations with smoking status, coeliac disease and the use of proton-pump inhibitors but not with other commonly reported pharmacological risk factors. Using the largest sample size to date, we confirmed a recently reported association with the MHC Ancestral 8.1 Haplotype. Downstream analysis suggests association with digestive tract morphogenesis. By calculating genetic risk scores, we also report suggestive evidence of shared genetic risk with Crohn's disease, but not with ulcerative colitis.

CONCLUSIONS

This report confirms the role of genetic determinants in the HLA in the pathogenesis of microscopic colitis. The genetic overlap with Crohn's disease suggests a common underlying mechanism of disease.

Refinement of the critical genomic region for hypoglycaemia in the Chromosome 9p deletion syndrome

Background: Large contiguous gene deletions at the distal end of the short arm of chromosome 9 result in the complex multi-organ condition chromosome 9p deletion syndrome.  A range of clinical features can result from these deletions with the most common being facial dysmorphisms and neurological impairment. Congenital hyperinsulinism is a rarely reported feature of the syndrome with the genetic mechanism for the dysregulated insulin secretion being unknown.  Methods: We studied the clinical and genetic characteristics of 12 individuals with chromosome 9p deletions who had a history of neonatal hypoglycaemia. Using off-target reads generated from targeted next-generation sequencing of the genes known to cause hyperinsulinaemic hypoglycaemia (n=9), or microarray analysis (n=3), we mapped the minimal shared deleted region on chromosome 9 in this cohort. Targeted sequencing was performed in three patients to search for a recessive mutation unmasked by the deletion. Results: In 10/12 patients with hypoglycaemia, hyperinsulinism was confirmed biochemically. A range of extra-pancreatic features were also reported in these patients consistent with the diagnosis of the Chromosome 9p deletion syndrome. The minimal deleted region was mapped to 7.2 Mb, encompassing 38 protein-coding genes. In silico analysis of these genes highlighted SMARCA2 and RFX3 as potential candidates for the hypoglycaemia. Targeted sequencing performed on three of the patients did not identify a second disease-causing variant within the minimal deleted region. Conclusions: This study identifies 9p deletions as an important cause of hyperinsulinaemic hypoglycaemia and increases the number of cases reported with 9p deletions and hypoglycaemia to 15 making this a more common feature of the syndrome than previously appreciated.  Whilst the precise genetic mechanism of the dysregulated insulin secretion could not be determined in these patients, mapping the deletion breakpoints highlighted potential candidate genes for hypoglycaemia within the deleted region.

Genome-wide association analysis of self-reported daytime sleepiness identifies 42 loci that suggest biological subtypes

Excessive daytime sleepiness (EDS) affects 10-20% of the population and is associated with substantial functional deficits. Here, we identify 42 loci for self-reported daytime sleepiness in GWAS of 452,071 individuals from the UK Biobank, with enrichment for genes expressed in brain tissues and in neuronal transmission pathways. We confirm the aggregate effect of a genetic risk score of 42 SNPs on daytime sleepiness in independent Scandinavian cohorts and on other sleep disorders (restless legs syndrome, insomnia) and sleep traits (duration, chronotype, accelerometer-derived sleep efficiency and daytime naps or inactivity). However, individual daytime sleepiness signals vary in their associations with objective short vs long sleep, and with markers of sleep continuity. The 42 sleepiness variants primarily cluster into two predominant composite biological subtypes - sleep propensity and sleep fragmentation. Shared genetic links are also seen with obesity, coronary heart disease, psychiatric diseases, cognitive traits and reproductive ageing.

Type 1 diabetes defined by severe insulin deficiency occurs after 30 years of age and is commonly treated as type 2 diabetes

AIMS/HYPOTHESIS

Late-onset type 1 diabetes can be difficult to identify. Measurement of endogenous insulin secretion using C-peptide provides a gold standard classification of diabetes type in longstanding diabetes that closely relates to treatment requirements. We aimed to determine the prevalence and characteristics of type 1 diabetes defined by severe endogenous insulin deficiency after age 30 and assess whether these individuals are identified and managed as having type 1 diabetes in clinical practice.

METHODS

We assessed the characteristics of type 1 diabetes defined by rapid insulin requirement (within 3 years of diagnosis) and severe endogenous insulin deficiency (non-fasting C-peptide <200 pmol/l) in 583 participants with insulin-treated diabetes, diagnosed after age 30, from the Diabetes Alliance for Research in England (DARE) population cohort. We compared characteristics with participants with retained endogenous insulin secretion (>600 pmol/l) and 220 participants with severe insulin deficiency who were diagnosed under age 30.

RESULTS

Twenty-one per cent of participants with insulin-treated diabetes who were diagnosed after age 30 met the study criteria for type 1 diabetes. Of these participants, 38% did not receive insulin at diagnosis, of whom 47% self-reported type 2 diabetes. Rapid insulin requirement was highly predictive of severe endogenous insulin deficiency: 85% required insulin within 1 year of diagnosis, and 47% of all those initially treated without insulin who progressed to insulin treatment within 3 years of diagnosis had severe endogenous insulin deficiency. Participants with late-onset type 1 diabetes defined by development of severe insulin deficiency had similar clinical characteristics to those with young-onset type 1 diabetes. However, those with later onset type 1 diabetes had a modestly lower type 1 diabetes genetic risk score (0.268 vs 0.279; p < 0.001 [expected type 2 diabetes population median, 0.231]), a higher islet autoantibody prevalence (GAD-, islet antigen 2 [IA2]- or zinc transporter protein 8 [ZnT8]-positive) of 78% at 13 years vs 62% at 26 years of diabetes duration; (p = 0.02), and were less likely to identify as having type 1 diabetes (79% vs 100%; p < 0.001) vs those with young-onset disease.

CONCLUSIONS/INTERPRETATION

Type 1 diabetes diagnosed over 30 years of age, defined by severe insulin deficiency, has similar clinical and biological characteristics to that occurring at younger ages, but is frequently not identified. Clinicians should be aware that patients progressing to insulin within 3 years of diagnosis have a high likelihood of type 1 diabetes, regardless of initial diagnosis.

Investigating causal relations between sleep traits and risk of breast cancer in women: mendelian randomisation study

OBJECTIVE

To examine whether sleep traits have a causal effect on risk of breast cancer.

DESIGN

Mendelian randomisation study.

SETTING

UK Biobank prospective cohort study and Breast Cancer Association Consortium (BCAC) case-control genome-wide association study.

PARTICIPANTS

156 848 women in the multivariable regression and one sample mendelian randomisation (MR) analysis in UK Biobank (7784 with a breast cancer diagnosis) and 122 977 breast cancer cases and 105 974 controls from BCAC in the two sample MR analysis.

EXPOSURES

Self reported chronotype (morning or evening preference), insomnia symptoms, and sleep duration in multivariable regression, and genetic variants robustly associated with these sleep traits.

MAIN OUTCOME MEASURE

Breast cancer diagnosis.

RESULTS

In multivariable regression analysis using UK Biobank data on breast cancer incidence, morning preference was inversely associated with breast cancer (hazard ratio 0.95, 95% confidence interval 0.93 to 0.98 per category increase), whereas there was little evidence for an association between sleep duration and insomnia symptoms. Using 341 single nucleotide polymorphisms (SNPs) associated with chronotype, 91 SNPs associated with sleep duration, and 57 SNPs associated with insomnia symptoms, one sample MR analysis in UK Biobank provided some supportive evidence for a protective effect of morning preference on breast cancer risk (0.85, 0.70, 1.03 per category increase) but imprecise estimates for sleep duration and insomnia symptoms. Two sample MR using data from BCAC supported findings for a protective effect of morning preference (inverse variance weighted odds ratio 0.88, 95% confidence interval 0.82 to 0.93 per category increase) and adverse effect of increased sleep duration (1.19, 1.02 to 1.39 per hour increase) on breast cancer risk (both oestrogen receptor positive and oestrogen receptor negative), whereas evidence for insomnia symptoms was inconsistent. Results were largely robust to sensitivity analyses accounting for horizontal pleiotropy.

CONCLUSIONS

Findings showed consistent evidence for a protective effect of morning preference and suggestive evidence for an adverse effect of increased sleep duration on breast cancer risk.

Using genetics to understand the causal influence of higher BMI on depression

BACKGROUND

Depression is more common in obese than non-obese individuals, especially in women, but the causal relationship between obesity and depression is complex and uncertain. Previous studies have used genetic variants associated with BMI to provide evidence that higher body mass index (BMI) causes depression, but have not tested whether this relationship is driven by the metabolic consequences of BMI nor for differences between men and women.

METHODS

We performed a Mendelian randomization study using 48 791 individuals with depression and 291 995 controls in the UK Biobank, to test for causal effects of higher BMI on depression (defined using self-report and Hospital Episode data). We used two genetic instruments, both representing higher BMI, but one with and one without its adverse metabolic consequences, in an attempt to 'uncouple' the psychological component of obesity from the metabolic consequences. We further tested causal relationships in men and women separately, and using subsets of BMI variants from known physiological pathways.

RESULTS

Higher BMI was strongly associated with higher odds of depression, especially in women. Mendelian randomization provided evidence that higher BMI partly causes depression. Using a 73-variant BMI genetic risk score, a genetically determined one standard deviation (1 SD) higher BMI (4.9 kg/m2) was associated with higher odds of depression in all individuals [odds ratio (OR): 1.18, 95% confidence interval (CI): 1.09, 1.28, P = 0.00007) and women only (OR: 1.24, 95% CI: 1.11, 1.39, P = 0.0001). Meta-analysis with 45 591 depression cases and 97 647 controls from the Psychiatric Genomics Consortium (PGC) strengthened the statistical confidence of the findings in all individuals. Similar effect size estimates were obtained using different Mendelian randomization methods, although not all reached P < 0.05. Using a metabolically favourable adiposity genetic risk score, and meta-analysing data from the UK biobank and PGC, a genetically determined 1 SD higher BMI (4.9 kg/m2) was associated with higher odds of depression in all individuals (OR: 1.26, 95% CI: 1.06, 1.50], P = 0.010), but with weaker statistical confidence.

CONCLUSIONS

Higher BMI, with and without its adverse metabolic consequences, is likely to have a causal role in determining the likelihood of an individual developing depression.

Meta-analysis of genome-wide association studies for height and body mass index in ∼700000 individuals of European ancestry

Recent genome-wide association studies (GWAS) of height and body mass index (BMI) in ∼250000 European participants have led to the discovery of ∼700 and ∼100 nearly independent single nucleotide polymorphisms (SNPs) associated with these traits, respectively. Here we combine summary statistics from those two studies with GWAS of height and BMI performed in ∼450000 UK Biobank participants of European ancestry. Overall, our combined GWAS meta-analysis reaches N ∼700000 individuals and substantially increases the number of GWAS signals associated with these traits. We identified 3290 and 941 near-independent SNPs associated with height and BMI, respectively (at a revised genome-wide significance threshold of P < 1 × 10-8), including 1185 height-associated SNPs and 751 BMI-associated SNPs located within loci not previously identified by these two GWAS. The near-independent genome-wide significant SNPs explain ∼24.6% of the variance of height and ∼6.0% of the variance of BMI in an independent sample from the Health and Retirement Study (HRS). Correlations between polygenic scores based upon these SNPs with actual height and BMI in HRS participants were ∼0.44 and ∼0.22, respectively. From analyses of integrating GWAS and expression quantitative trait loci (eQTL) data by summary-data-based Mendelian randomization, we identified an enrichment of eQTLs among lead height and BMI signals, prioritizing 610 and 138 genes, respectively. Our study demonstrates that, as previously predicted, increasing GWAS sample sizes continues to deliver, by the discovery of new loci, increasing prediction accuracy and providing additional data to achieve deeper insight into complex trait biology. All summary statistics are made available for follow-up studies.

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.

Genome-wide association analysis of diverticular disease points towards neuromuscular, connective tissue and epithelial pathomechanisms

OBJECTIVE

Diverticular disease is a common complex disorder characterised by mucosal outpouchings of the colonic wall that manifests through complications such as diverticulitis, perforation and bleeding. We report the to date largest genome-wide association study (GWAS) to identify genetic risk factors for diverticular disease.

DESIGN

Discovery GWAS analysis was performed on UK Biobank imputed genotypes using 31 964 cases and 419 135 controls of European descent. Associations were replicated in a European sample of 3893 cases and 2829 diverticula-free controls and evaluated for risk contribution to diverticulitis and uncomplicated diverticulosis. Transcripts at top 20 replicating loci were analysed by real-time quatitative PCR in preparations of the mucosal, submucosal and muscular layer of colon. The localisation of expressed protein at selected loci was investigated by immunohistochemistry.

RESULTS

We discovered 48 risk loci, of which 12 are novel, with genome-wide significance and consistent OR in the replication sample. Nominal replication (p<0.05) was observed for 27 loci, and additional 8 in meta-analysis with a population-based cohort. The most significant novel risk variant rs9960286 is located near CTAGE1 with a p value of 2.3×10^-10 and 0.002 (OR_allelic=1.14 (95% CI 1.05 to 1.24)) in the replication analysis. Four loci showed stronger effects for diverticulitis, PHGR1 (OR 1.32, 95% CI 1.12 to 1.56), FAM155A-2 (OR 1.21, 95% CI 1.04 to 1.42), CALCB (OR 1.17, 95% CI 1.03 to 1.33) and S100A10 (OR 1.17, 95% CI 1.03 to 1.33).

CONCLUSION

In silico analyses point to diverticulosis primarily as a disorder of intestinal neuromuscular function and of impaired connective fibre support, while an additional diverticulitis risk might be conferred by epithelial dysfunction.

Genetic studies of accelerometer-based sleep measures yield new insights into human sleep behaviour

Sleep is an essential human function but its regulation is poorly understood. Using accelerometer data from 85,670 UK Biobank participants, we perform a genome-wide association study of 8 derived sleep traits representing sleep quality, quantity and timing, and validate our findings in 5,819 individuals. We identify 47 genetic associations at P < 5 × 10-8, of which 20 reach a stricter threshold of P < 8 × 10-10. These include 26 novel associations with measures of sleep quality and 10 with nocturnal sleep duration. The majority of identified variants associate with a single sleep trait, except for variants previously associated with restless legs syndrome. For sleep duration we identify a missense variant (p.Tyr727Cys) in PDE11A as the likely causal variant. As a group, sleep quality loci are enriched for serotonin processing genes. Although accelerometer-derived measures of sleep are imperfect and may be affected by restless legs syndrome, these findings provide new biological insights into sleep compared to previous efforts based on self-report sleep measures.

Chronotype Genetic Variant in PER2 is Associated with Intrinsic Circadian Period in Humans

The PERIOD2 (PER2) gene is a core molecular component of the circadian clock and plays an important role in the generation and maintenance of daily rhythms. Rs35333999, a missense variant of PER2 common in European populations, has been shown to associate with later chronotype. Chronotype relates to the timing of biological and behavioral activities, including when we sleep, eat, and exercise, and later chronotype is associated with longer intrinsic circadian period (cycle length), a fundamental property of the circadian system. Thus, we tested whether this PER2 variant was associated with circadian period and found significant associations with longer intrinsic circadian period as measured under forced desynchrony protocols, the 'gold standard' for intrinsic circadian period assessment. Minor allele (T) carriers exhibited significantly longer circadian periods when determinations were based on either core body temperature or plasma melatonin measurements, as compared to non-carriers (by 12 and 11 min, respectively; accounting for ~7% of inter-individual variance). These findings provide a possible underlying biological mechanism for inter-individual differences in chronotype, and support the central role of PER2 in the human circadian timing system.

Genome-wide association study identifies genetic loci for self-reported habitual sleep duration supported by accelerometer-derived estimates

Sleep is an essential state of decreased activity and alertness but molecular factors regulating sleep duration remain unknown. Through genome-wide association analysis in 446,118 adults of European ancestry from the UK Biobank, we identify 78 loci for self-reported habitual sleep duration (p < 5 × 10⁻⁸; 43 loci at p < 6 × 10⁻⁹). Replication is observed for PAX8, VRK2, and FBXL12/UBL5/PIN1 loci in the CHARGE study (n = 47,180; p < 6.3 × 10⁻⁴), and 55 signals show sign-concordant effects. The 78 loci further associate with accelerometer-derived sleep duration, daytime inactivity, sleep efficiency and number of sleep bouts in secondary analysis (n = 85,499). Loci are enriched for pathways including striatum and subpallium development, mechanosensory response, dopamine binding, synaptic neurotransmission and plasticity, among others. Genetic correlation indicates shared links with anthropometric, cognitive, metabolic, and psychiatric traits and two-sample Mendelian randomization highlights a bidirectional causal link with schizophrenia. This work provides insights into the genetic basis for inter-individual variation in sleep duration implicating multiple biological pathways.

Sleep is essential for homeostasis and insufficient or excessive sleep are associated with adverse outcomes. Here, the authors perform GWAS for self-reported habitual sleep duration in adults, supported by accelerometer-derived measures, and identify genetic correlation with psychiatric and metabolic traits

Biological and clinical insights from genetics of insomnia symptoms

Insomnia is a common disorder linked with adverse long-term medical and psychiatric outcomes. The underlying pathophysiological processes and causal relationships of insomnia with disease are poorly understood. Here we identified 57 loci for self-reported insomnia symptoms in the UK Biobank (n = 453,379) and confirmed their effects on self-reported insomnia symptoms in the HUNT Study (n = 14,923 cases and 47,610 controls), physician-diagnosed insomnia in the Partners Biobank (n = 2,217 cases and 14,240 controls), and accelerometer-derived measures of sleep efficiency and sleep duration in the UK Biobank (n = 83,726). Our results suggest enrichment of genes involved in ubiquitin-mediated proteolysis and of genes expressed in multiple brain regions, skeletal muscle, and adrenal glands. Evidence of shared genetic factors was found between frequent insomnia symptoms and restless legs syndrome, aging, and cardiometabolic, behavioral, psychiatric, and reproductive traits. Evidence was found for a possible causal link between insomnia symptoms and coronary artery disease, depressive symptoms, and subjective well-being.

Association of Genetic Variants in NUDT15 With Thiopurine-Induced Myelosuppression in Patients With Inflammatory Bowel Disease

IMPORTANCE

Use of thiopurines may be limited by myelosuppression. TPMT pharmacogenetic testing identifies only 25% of at-risk patients of European ancestry. Among patients of East Asian ancestry, NUDT15 variants are associated with thiopurine-induced myelosuppression (TIM).

OBJECTIVE

To identify genetic variants associated with TIM among patients of European ancestry with inflammatory bowel disease (IBD).

DESIGN, SETTING, AND PARTICIPANTS

Case-control study of 491 patients affected by TIM and 679 thiopurine-tolerant unaffected patients who were recruited from 89 international sites between March 2012 and November 2015. Genome-wide association studies (GWAS) and exome-wide association studies (EWAS) were conducted in patients of European ancestry. The replication cohort comprised 73 patients affected by TIM and 840 thiopurine-tolerant unaffected patients.

EXPOSURES

Genetic variants associated with TIM.

MAIN OUTCOMES AND MEASURES

Thiopurine-induced myelosuppression, defined as a decline in absolute white blood cell count to 2.5 × 109/L or less or a decline in absolute neutrophil cell count to 1.0 × 109/L or less leading to a dose reduction or drug withdrawal.

RESULTS

Among 1077 patients (398 affected and 679 unaffected; median age at IBD diagnosis, 31.0 years [interquartile range, 21.2 to 44.1 years]; 540 [50%] women; 602 [56%] diagnosed as having Crohn disease), 919 (311 affected and 608 unaffected) were included in the GWAS analysis and 961 (328 affected and 633 unaffected) in the EWAS analysis. The GWAS analysis confirmed association of TPMT (chromosome 6, rs11969064) with TIM (30.5% [95/311] affected vs 16.4% [100/608] unaffected patients; odds ratio [OR], 2.3 [95% CI, 1.7 to 3.1], P = 5.2 × 10-9). The EWAS analysis demonstrated an association with an in-frame deletion in NUDT15 (chromosome 13, rs746071566) and TIM (5.8% [19/328] affected vs 0.2% [1/633] unaffected patients; OR, 38.2 [95% CI, 5.1 to 286.1], P = 1.3 × 10-8), which was replicated in a different cohort (2.7% [2/73] affected vs 0.2% [2/840] unaffected patients; OR, 11.8 [95% CI, 1.6 to 85.0], P = .03). Carriage of any of 3 coding NUDT15 variants was associated with an increased risk (OR, 27.3 [95% CI, 9.3 to 116.7], P = 1.1 × 10-7) of TIM, independent of TPMT genotype and thiopurine dose.

CONCLUSIONS AND RELEVANCE

Among patients of European ancestry with IBD, variants in NUDT15 were associated with increased risk of TIM. These findings suggest that NUDT15 genotyping may be considered prior to initiation of thiopurine therapy; however, further study including additional validation in independent cohorts is required.

Assessing the Pathogenicity, Penetrance, and Expressivity of Putative Disease-Causing Variants in a Population Setting

More than 100,000 genetic variants are classified as disease causing in public databases. However, the true penetrance of many of these rare alleles is uncertain and might be over-estimated by clinical ascertainment. Here, we use data from 379,768 UK Biobank (UKB) participants of European ancestry to assess the pathogenicity and penetrance of putatively clinically important rare variants. Although rare variants are harder to genotype accurately than common variants, we were able to classify as high quality 1,244 of 4,585 (27%) putatively clinically relevant rare (MAF < 1%) variants genotyped on the UKB microarray. We defined as "clinically relevant" variants that were classified as either pathogenic or likely pathogenic in ClinVar or are in genes known to cause two specific monogenic diseases: maturity-onset diabetes of the young (MODY) and severe developmental disorders (DDs). We assessed the penetrance and pathogenicity of these high-quality variants by testing their association with 401 clinically relevant traits. 27 of the variants were associated with a UKB trait, and we were able to refine the penetrance estimate for some of the variants. For example, the HNF4A c.340C>T (p.Arg114Trp) (GenBank: NM_175914.4) variant associated with diabetes is <10% penetrant by the time an individual is 40 years old. We also observed associations with relevant traits for heterozygous carriers of some rare recessive conditions, e.g., heterozygous carriers of the ERCC4 c.2395C>T (p.Arg799Trp) variant that causes Xeroderma pigmentosum were more susceptible to sunburn. Finally, we refute the previous disease association of RNF135 in developmental disorders. In conclusion, this study shows that very large population-based studies will help refine our understanding of the pathogenicity of rare genetic variants.