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.

Development and Standardization of an Improved Type 1 Diabetes Genetic Risk Score for Use in Newborn Screening and Incident Diagnosis

OBJECTIVE

Previously generated genetic risk scores (GRSs) for type 1 diabetes (T1D) have not captured all known information at non-HLA loci or, particularly, at HLA risk loci. We aimed to more completely incorporate HLA alleles, their interactions, and recently discovered non-HLA loci into an improved T1D GRS (termed the "T1D GRS2") to better discriminate diabetes subtypes and to predict T1D in newborn screening studies.

RESEARCH DESIGN AND METHODS

In 6,481 case and 9,247 control subjects from the Type 1 Diabetes Genetics Consortium, we analyzed variants associated with T1D both in the HLA region and across the genome. We modeled interactions between variants marking strongly associated HLA haplotypes and generated odds ratios to create the improved GRS, the T1D GRS2. We validated our findings in UK Biobank. We assessed the impact of the T1D GRS2 in newborn screening and diabetes classification and sought to provide a framework for comparison with previous scores.

RESULTS

The T1D GRS2 used 67 single nucleotide polymorphisms (SNPs) and accounted for interactions between 18 HLA DR-DQ haplotype combinations. The T1D GRS2 was highly discriminative for all T1D (area under the curve [AUC] 0.92; P < 0.0001 vs. older scores) and even more discriminative for early-onset T1D (AUC 0.96). In simulated newborn screening, the T1D GRS2 was nearly twice as efficient as HLA genotyping alone and 50% better than current genetic scores in general population T1D prediction.

CONCLUSIONS

An improved T1D GRS, the T1D GRS2, is highly useful for classifying adult incident diabetes type and improving newborn screening. Given the cost-effectiveness of SNP genotyping, this approach has great clinical and research potential in T1D.

Disentangling the genetics of lean mass

Background

Lean body mass (LM) plays an important role in mobility and metabolic function. We previously identified five loci associated with LM adjusted for fat mass in kilograms. Such an adjustment may reduce the power to identify genetic signals having an association with both lean mass and fat mass.

Objectives

To determine the impact of different fat mass adjustments on genetic architecture of LM and identify additional LM loci.

Methods

We performed genome-wide association analyses for whole-body LM (20 cohorts of European ancestry with n = 38,292) measured using dual-energy X-ray absorptiometry) or bioelectrical impedance analysis, adjusted for sex, age, age2, and height with or without fat mass adjustments (Model 1 no fat adjustment; Model 2 adjustment for fat mass as a percentage of body mass; Model 3 adjustment for fat mass in kilograms).

Results

Seven single-nucleotide polymorphisms (SNPs) in separate loci, including one novel LM locus (TNRC6B), were successfully replicated in an additional 47,227 individuals from 29 cohorts. Based on the strengths of the associations in Model 1 vs Model 3, we divided the LM loci into those with an effect on both lean mass and fat mass in the same direction and refer to those as "sumo wrestler" loci (FTO and MC4R). In contrast, loci with an impact specifically on LM were termed "body builder" loci (VCAN and ADAMTSL3). Using existing available genome-wide association study databases, LM increasing alleles of SNPs in sumo wrestler loci were associated with an adverse metabolic profile, whereas LM increasing alleles of SNPs in "body builder" loci were associated with metabolic protection.

Conclusions

In conclusion, we identified one novel LM locus (TNRC6B). Our results suggest that a genetically determined increase in lean mass might exert either harmful or protective effects on metabolic traits, depending on its relation to fat mass.

Genome-wide association analyses of chronotype in 697,828 individuals provides insights into circadian rhythms

Being a morning person is a behavioural indicator of a person’s underlying circadian rhythm. Using genome-wide data from 697,828 UK Biobank and 23andMe participants we increase the number of genetic loci associated with being a morning person from 24 to 351. Using data from 85,760 individuals with activity-monitor derived measures of sleep timing we find that the chronotype loci associate with sleep timing: the mean sleep timing of the 5% of individuals carrying the most morningness alleles is 25 min earlier than the 5% carrying the fewest. The loci are enriched for genes involved in circadian regulation, cAMP, glutamate and insulin signalling pathways, and those expressed in the retina, hindbrain, hypothalamus, and pituitary. Using Mendelian Randomisation, we show that being a morning person is causally associated with better mental health but does not affect BMI or risk of Type 2 diabetes. This study offers insights into circadian biology and its links to disease in humans.

GWAS have previously found 24 genomic loci associated with chronotype, an individual’s preference for early or late sleep timing. Here, the authors identify 327 additional loci in a sample of 697,828 individuals and further explore the relationships of chronotype with metabolic and psychiatric diseases.

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⁻⁸. 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.

GWAS Identifies Risk Locus for Erectile Dysfunction and Implicates Hypothalamic Neurobiology and Diabetes in Etiology

Erectile dysfunction (ED) is a common condition affecting more than 20% of men over 60 years, yet little is known about its genetic architecture. We performed a genome-wide association study of ED in 6,175 case subjects among 223,805 European men and identified one locus at 6q16.3 (lead variant rs57989773, OR 1.20 per C-allele; p = 5.71 × 10⁻¹⁴), located between MCHR2 and SIM1. In silico analysis suggests SIM1 to confer ED risk through hypothalamic dysregulation. Mendelian randomization provides evidence that genetic risk of type 2 diabetes mellitus is a cause of ED (OR 1.11 per 1-log unit higher risk of type 2 diabetes). These findings provide insights into the biological underpinnings and the causes of ED and may help prioritize the development of future therapies for this common disorder.

Genome-Wide and Abdominal MRI Data Provide Evidence That a Genetically Determined Favorable Adiposity Phenotype Is Characterized by Lower Ectopic Liver Fat and Lower Risk of Type 2 Diabetes, Heart Disease, and Hypertension

Recent genetic studies have identified alleles associated with opposite effects on adiposity and risk of type 2 diabetes. We aimed to identify more of these variants and test the hypothesis that such favorable adiposity alleles are associated with higher subcutaneous fat and lower ectopic fat. We combined MRI data with genome-wide association studies of body fat percentage (%) and metabolic traits. We report 14 alleles, including 7 newly characterized alleles, associated with higher adiposity but a favorable metabolic profile. Consistent with previous studies, individuals carrying more favorable adiposity alleles had higher body fat % and higher BMI but lower risk of type 2 diabetes, heart disease, and hypertension. These individuals also had higher subcutaneous fat but lower liver fat and a lower visceral-to-subcutaneous adipose tissue ratio. Individual alleles associated with higher body fat % but lower liver fat and lower risk of type 2 diabetes included those in PPARG, GRB14, and IRS1, whereas the allele in ANKRD55 was paradoxically associated with higher visceral fat but lower risk of type 2 diabetes. Most identified favorable adiposity alleles are associated with higher subcutaneous and lower liver fat, a mechanism consistent with the beneficial effects of storing excess triglycerides in metabolically low-risk depots.

The Common HNF1A Variant I27L Is a Modifier of Age at Diabetes Diagnosis in Individuals With HNF1A-MODY

There is wide variation in the age at diagnosis of diabetes in individuals with maturity-onset diabetes of the young (MODY) due to a mutation in the HNF1A gene. We hypothesized that common variants at the HNF1A locus (rs1169288 [I27L], rs1800574 [A98V]), which are associated with type 2 diabetes susceptibility, may modify age at diabetes diagnosis in individuals with HNF1A-MODY. Meta-analysis of two independent cohorts, comprising 781 individuals with HNF1A-MODY, found no significant associations between genotype and age at diagnosis. However after stratifying according to type of mutation (protein-truncating variant [PTV] or missense), we found each 27L allele to be associated with a 1.6-year decrease (95% CI -2.6, -0.7) in age at diagnosis, specifically in the subset (n = 444) of individuals with a PTV. The effect size was similar and significant across the two independent cohorts of individuals with HNF1A-MODY. We report a robust genetic modifier of HNF1A-MODY age at diagnosis that further illustrates the strong effect of genetic variation within HNF1A upon diabetes phenotype.

Sex and BMI Alter the Benefits and Risks of Sulfonylureas and Thiazolidinediones in Type 2 Diabetes: A Framework for Evaluating Stratification Using Routine Clinical and Individual Trial Data

OBJECTIVE

The choice of therapy for type 2 diabetes after metformin is guided by overall estimates of glycemic response and side effects seen in large cohorts. A stratified approach to therapy would aim to improve on this by identifying subgroups of patients whose glycemic response or risk of side effects differs markedly. We assessed whether simple clinical characteristics could identify patients with differing glycemic response and side effects with sulfonylureas and thiazolidinediones.

RESEARCH DESIGN AND METHODS

We studied 22,379 patients starting sulfonylurea or thiazolidinedione therapy in the U.K. Clinical Practice Research Datalink (CPRD) to identify features associated with increased 1-year HbA_1c fall with one therapy class and reduced fall with the second. We then assessed whether prespecified patient subgroups defined by the differential clinical factors showed differing 5-year glycemic response and side effects with sulfonylureas and thiazolidinediones using individual randomized trial data from ADOPT (A Diabetes Outcome Progression Trial) (first-line therapy, n = 2,725) and RECORD (Rosiglitazone Evaluated for Cardiovascular Outcomes and Regulation of Glycemia in Diabetes) (second-line therapy, n = 2,222). Further replication was conducted using routine clinical data from GoDARTS (Genetics of Diabetes Audit and Research in Tayside Scotland) (n = 1,977).

RESULTS

In CPRD, male sex and lower BMI were associated with greater glycemic response with sulfonylureas and a lesser response with thiazolidinediones (both P < 0.001). In ADOPT and RECORD, nonobese males had a greater overall HbA_1c reduction with sulfonylureas than with thiazolidinediones (P < 0.001); in contrast, obese females had a greater HbA_1c reduction with thiazolidinediones than with sulfonylureas (P < 0.001). Weight gain and edema risk with thiazolidinediones were greatest in obese females; however, hypoglycemia risk with sulfonylureas was similar across all subgroups.

CONCLUSIONS

Patient subgroups defined by sex and BMI have different patterns of benefits and risks on thiazolidinedione and sulfonylurea therapy. Subgroup-specific estimates can inform discussion about the choice of therapy after metformin for an individual patient. Our approach using routine and shared trial data provides a framework for future stratification research in type 2 diabetes.

A Type 1 Diabetes Genetic Risk Score Predicts Progression of Islet Autoimmunity and Development of Type 1 Diabetes in Individuals at Risk

OBJECTIVE

We tested the ability of a type 1 diabetes (T1D) genetic risk score (GRS) to predict progression of islet autoimmunity and T1D in at-risk individuals.

RESEARCH DESIGN AND METHODS

We studied the 1,244 TrialNet Pathway to Prevention study participants (T1D patients' relatives without diabetes and with one or more positive autoantibodies) who were genotyped with Illumina ImmunoChip (median [range] age at initial autoantibody determination 11.1 years [1.2-51.8], 48% male, 80.5% non-Hispanic white, median follow-up 5.4 years). Of 291 participants with a single positive autoantibody at screening, 157 converted to multiple autoantibody positivity and 55 developed diabetes. Of 953 participants with multiple positive autoantibodies at screening, 419 developed diabetes. We calculated the T1D GRS from 30 T1D-associated single nucleotide polymorphisms. We used multivariable Cox regression models, time-dependent receiver operating characteristic curves, and area under the curve (AUC) measures to evaluate prognostic utility of T1D GRS, age, sex, Diabetes Prevention Trial-Type 1 (DPT-1) Risk Score, positive autoantibody number or type, HLA DR3/DR4-DQ8 status, and race/ethnicity. We used recursive partitioning analyses to identify cut points in continuous variables.

RESULTS

Higher T1D GRS significantly increased the rate of progression to T1D adjusting for DPT-1 Risk Score, age, number of positive autoantibodies, sex, and ethnicity (hazard ratio [HR] 1.29 for a 0.05 increase, 95% CI 1.06-1.6; P = 0.011). Progression to T1D was best predicted by a combined model with GRS, number of positive autoantibodies, DPT-1 Risk Score, and age (7-year time-integrated AUC = 0.79, 5-year AUC = 0.73). Higher GRS was significantly associated with increased progression rate from single to multiple positive autoantibodies after adjusting for age, autoantibody type, ethnicity, and sex (HR 2.27 for GRS >0.295, 95% CI 1.47-3.51; P = 0.0002).

CONCLUSIONS

The T1D GRS independently predicts progression to T1D and improves prediction along T1D stages in autoantibody-positive relatives.

PLIN1 Haploinsufficiency Is Not Associated With Lipodystrophy

CONTEXT

Monogenic partial lipodystrophy is a genetically heterogeneous disease where only variants with specific genetic mechanisms are causative. Three heterozygous protein extending frameshift variants in PLIN1 have been reported to cause a phenotype of partial lipodystrophy and insulin resistance.

OBJECTIVE

We investigated if null variants in PLIN1 cause lipodystrophy.

METHODS

As part of a targeted sequencing panel test, we sequenced PLIN1 in 2208 individuals. We also investigated the frequency of PLIN1 variants in the gnomAD database, and the type 2 diabetes knowledge portal.

RESULTS

We identified 6/2208 (1 in 368) individuals with a PLIN1 null variant. None of these individuals had clinical or biochemical evidence of overt lipodystrophy. Additionally, 14/17,000 (1 in 1214) individuals with PLIN1 null variants in the type 2 diabetes knowledge portal showed no association with biomarkers of lipodystrophy. PLIN1 null variants occur too frequently in gnomAD (126/138,632; 1 in 1100) to be a cause of rare overt monogenic partial lipodystrophy.

CONCLUSIONS

Our study suggests that heterozygous variants that are predicted to result in PLIN1 haploinsufficiency are not a cause of familial partial lipodystrophy and should not be reported as disease-causing variants by diagnostic genetic testing laboratories. This finding is in keeping with other known monogenic causes of lipodystrophy, such as PPARG and LMNA, where only variants with specific genetic mechanisms cause lipodystrophy.

Estimating sleep parameters using an accelerometer without sleep diary

Wrist worn raw-data accelerometers are used increasingly in large-scale population research. We examined whether sleep parameters can be estimated from these data in the absence of sleep diaries. Our heuristic algorithm uses the variance in estimated z-axis angle and makes basic assumptions about sleep interruptions. Detected sleep period time window (SPT-window) was compared against sleep diary in 3752 participants (range = 60-82 years) and polysomnography in sleep clinic patients (N = 28) and in healthy good sleepers (N = 22). The SPT-window derived from the algorithm was 10.9 and 2.9 minutes longer compared with sleep diary in men and women, respectively. Mean C-statistic to detect the SPT-window compared to polysomnography was 0.86 and 0.83 in clinic-based and healthy sleepers, respectively. We demonstrated the accuracy of our algorithm to detect the SPT-window. The value of this algorithm lies in studies such as UK Biobank where a sleep diary was not used.

Meta-genome-wide association studies identify a locus on chromosome 1 and multiple variants in the MHC region for serum C-peptide in type 1 diabetes

AIMS/HYPOTHESIS

The aim of this study was to identify genetic variants associated with beta cell function in type 1 diabetes, as measured by serum C-peptide levels, through meta-genome-wide association studies (meta-GWAS).

METHODS

We performed a meta-GWAS to combine the results from five studies in type 1 diabetes with cross-sectionally measured stimulated, fasting or random C-peptide levels, including 3479 European participants. The p values across studies were combined, taking into account sample size and direction of effect. We also performed separate meta-GWAS for stimulated (n = 1303), fasting (n = 2019) and random (n = 1497) C-peptide levels.

RESULTS

In the meta-GWAS for stimulated/fasting/random C-peptide levels, a SNP on chromosome 1, rs559047 (Chr1:238753916, T>A, minor allele frequency [MAF] 0.24-0.26), was associated with C-peptide (p = 4.13 × 10^-8), meeting the genome-wide significance threshold (p < 5 × 10^-8). In the same meta-GWAS, a locus in the MHC region (rs9260151) was close to the genome-wide significance threshold (Chr6:29911030, C>T, MAF 0.07-0.10, p = 8.43 × 10^-8). In the stimulated C-peptide meta-GWAS, rs61211515 (Chr6:30100975, T/-, MAF 0.17-0.19) in the MHC region was associated with stimulated C-peptide (β [SE] = - 0.39 [0.07], p = 9.72 × 10^-8). rs61211515 was also associated with the rate of stimulated C-peptide decline over time in a subset of individuals (n = 258) with annual repeated measures for up to 6 years (p = 0.02). In the meta-GWAS of random C-peptide, another MHC region, SNP rs3135002 (Chr6:32668439, C>A, MAF 0.02-0.06), was associated with C-peptide (p = 3.49 × 10^-8). Conditional analyses suggested that the three identified variants in the MHC region were independent of each other. rs9260151 and rs3135002 have been associated with type 1 diabetes, whereas rs559047 and rs61211515 have not been associated with a risk of developing type 1 diabetes.

CONCLUSIONS/INTERPRETATION

We identified a locus on chromosome 1 and multiple variants in the MHC region, at least some of which were distinct from type 1 diabetes risk loci, that were associated with C-peptide, suggesting partly non-overlapping mechanisms for the development and progression of type 1 diabetes. These associations need to be validated in independent populations. Further investigations could provide insights into mechanisms of beta cell loss and opportunities to preserve beta cell function.

Clinical and research uses of genetic risk scores in type 1 diabetes

Type 1 diabetes (T1D) is a chronic disease of high blood glucose caused by autoimmune destruction of pancreatic beta cells eventually resulting in severe insulin deficiency. T1D has a significant heritable risk. Genetic associations found are particularly strong in the HLA class II region but T1D is a polygenic disease associated with over 60 loci across the genome. Polygenic risk scores are one method of summing these genetic risk elements as a single continuous variable. This review discusses the clinical and research utility of genetic risk scores in T1D particularly in disease prediction and progression. We also explore creative uses of genetic risk scores in big data and the limitations of using a genetic risk score. The increase in publically available genetic data and rapid fall in costs of genotyping mean that a T1D genetic risk score (T1D GRS) is likely to prove useful for disease prediction, discrimination, investigation of unusual cohorts, and investigation of biology in large datasets where genetic data are available.

A Common Allele in FGF21 Associated with Sugar Intake Is Associated with Body Shape, Lower Total Body-Fat Percentage, and Higher Blood Pressure

Fibroblast growth factor 21 (FGF21) is a hormone that has insulin-sensitizing properties. Some trials of FGF21 analogs show weight loss and lipid-lowering effects. Recent studies have shown that a common allele in the FGF21 gene alters the balance of macronutrients consumed, but there was little evidence of an effect on metabolic traits. We studied a common FGF21 allele (A:rs838133) in 451,099 people from the UK Biobank study, aiming to use the human allele to inform potential adverse and beneficial effects of targeting FGF21. We replicated the association between the A allele and higher percentage carbohydrate intake. We then showed that this allele is more strongly associated with higher blood pressure and waist-hip ratio, despite an association with lower total body-fat percentage, than it is with BMI or type 2 diabetes. These human phenotypes of variation in the FGF21 gene will inform research into FGF21's mechanisms and therapeutic potential.

Precision Medicine in Type 2 Diabetes: Clinical Markers of Insulin Resistance Are Associated With Altered Short- and Long-term Glycemic Response to DPP-4 Inhibitor Therapy

OBJECTIVE

A precision approach to type 2 diabetes therapy would aim to target treatment according to patient characteristics. We examined if measures of insulin resistance and secretion were associated with glycemic response to dipeptidyl peptidase 4 (DPP-4) inhibitor therapy.

RESEARCH DESIGN AND METHODS

We evaluated whether markers of insulin resistance and insulin secretion were associated with 6-month glycemic response in a prospective study of noninsulin-treated participants starting DPP-4 inhibitor therapy (Predicting Response to Incretin Based Agents [PRIBA] study; n = 254), with replication for routinely available markers in U.K. electronic health care records (Clinical Practice Research Datalink [CPRD]; n = 23,001). In CPRD, we evaluated associations between baseline markers and 3-year durability of response. To test the specificity of findings, we repeated analyses for glucagon-like peptide 1 (GLP-1) receptor agonists (PRIBA, n = 339; CPRD, n = 4,464).

RESULTS

In PRIBA, markers of higher insulin resistance (higher fasting C-peptide [P = 0.03], HOMA2 insulin resistance [P = 0.01], and triglycerides [P < 0.01]) were associated with reduced 6-month HbA_1c response to DPP-4 inhibitors. In CPRD, higher triglycerides and BMI were associated with reduced HbA_1c response (both P < 0.01). A subgroup defined by obesity (BMI ≥30 kg/m²) and high triglycerides (≥2.3 mmol/L) had reduced 6-month response in both data sets (PRIBA HbA_1c reduction 5.3 [95% CI 1.8, 8.6] mmol/mol [0.5%] [obese and high triglycerides] vs. 11.3 [8.4, 14.1] mmol/mol [1.0%] [nonobese and normal triglycerides]; P = 0.01). In CPRD, the obese, high- triglycerides subgroup also had less durable response (hazard ratio 1.28 [1.16, 1.41]; P < 0.001). There was no association between markers of insulin resistance and response to GLP-1 receptor agonists.

CONCLUSIONS

Markers of higher insulin resistance are consistently associated with reduced glycemic response to DPP-4 inhibitors. This finding provides a starting point for the application of a precision diabetes approach to DPP-4 inhibitor therapy.

Genetic scores to stratify risk of developing multiple islet autoantibodies and type 1 diabetes: A prospective study in children

BACKGROUND

Around 0.3% of newborns will develop autoimmunity to pancreatic beta cells in childhood and subsequently develop type 1 diabetes before adulthood. Primary prevention of type 1 diabetes will require early intervention in genetically at-risk infants. The objective of this study was to determine to what extent genetic scores (two previous genetic scores and a merged genetic score) can improve the prediction of type 1 diabetes.

METHODS AND FINDINGS

The Environmental Determinants of Diabetes in the Young (TEDDY) study followed genetically at-risk children at 3- to 6-monthly intervals from birth for the development of islet autoantibodies and type 1 diabetes. Infants were enrolled between 1 September 2004 and 28 February 2010 and monitored until 31 May 2016. The risk (positive predictive value) for developing multiple islet autoantibodies (pre-symptomatic type 1 diabetes) and type 1 diabetes was determined in 4,543 children who had no first-degree relatives with type 1 diabetes and either a heterozygous HLA DR3 and DR4-DQ8 risk genotype or a homozygous DR4-DQ8 genotype, and in 3,498 of these children in whom genetic scores were calculated from 41 single nucleotide polymorphisms. In the children with the HLA risk genotypes, risk for developing multiple islet autoantibodies was 5.8% (95% CI 5.0%-6.6%) by age 6 years, and risk for diabetes by age 10 years was 3.7% (95% CI 3.0%-4.4%). Risk for developing multiple islet autoantibodies was 11.0% (95% CI 8.7%-13.3%) in children with a merged genetic score of >14.4 (upper quartile; n = 907) compared to 4.1% (95% CI 3.3%-4.9%, P < 0.001) in children with a genetic score of ≤14.4 (n = 2,591). Risk for developing diabetes by age 10 years was 7.6% (95% CI 5.3%-9.9%) in children with a merged score of >14.4 compared with 2.7% (95% CI 1.9%-3.6%) in children with a score of ≤14.4 (P < 0.001). Of 173 children with multiple islet autoantibodies by age 6 years and 107 children with diabetes by age 10 years, 82 (sensitivity, 47.4%; 95% CI 40.1%-54.8%) and 52 (sensitivity, 48.6%, 95% CI 39.3%-60.0%), respectively, had a score >14.4. Scores were higher in European versus US children (P = 0.003). In children with a merged score of >14.4, risk for multiple islet autoantibodies was similar and consistently >10% in Europe and in the US; risk was greater in males than in females (P = 0.01). Limitations of the study include that the genetic scores were originally developed from case-control studies of clinical diabetes in individuals of mainly European decent. It is, therefore, possible that it may not be suitable to all populations.

CONCLUSIONS

A type 1 diabetes genetic score identified infants without family history of type 1 diabetes who had a greater than 10% risk for pre-symptomatic type 1 diabetes, and a nearly 2-fold higher risk than children identified by high-risk HLA genotypes alone. This finding extends the possibilities for enrolling children into type 1 diabetes primary prevention trials.

Frequency and phenotype of type 1 diabetes in the first six decades of life: a cross-sectional, genetically stratified survival analysis from UK Biobank

BACKGROUND

Type 1 diabetes is typically considered a disease of children and young adults. Genetic susceptibility to young-onset type 1 diabetes is well defined and does not predispose to type 2 diabetes. It is not known how frequently genetic susceptibility to type 1 diabetes leads to a diagnosis of diabetes after age 30 years. We aimed to investigate the frequency and phenotype of type 1 diabetes resulting from high genetic susceptibility in the first six decades of life.

METHODS

In this cross-sectional analysis, we used a type 1 diabetes genetic risk score based on 29 common variants to identify individuals of white European descent in UK Biobank in the half of the population with high or low genetic susceptibility to type 1 diabetes. We used Kaplan-Meier analysis to evaluate the number of cases of diabetes in both groups in the first six decades of life. We genetically defined type 1 diabetes as the additional cases of diabetes that occurred in the high genetic susceptibility group compared with the low genetic susceptibility group. All remaining cases were defined as type 2 diabetes. We assessed the clinical characteristics of the groups with genetically defined type 1 or type 2 diabetes.

FINDINGS

13 250 (3·5%) of 379 511 white European individuals in UK Biobank had developed diabetes in the first six decades of life. 1286 more cases of diabetes were in the half of the population with high genetic susceptibility to type 1 diabetes than in the half of the population with low genetic susceptibility. These genetically defined cases of type 1 diabetes were distributed across all ages of diagnosis; 537 (42%) were in individuals diagnosed when aged 31-60 years, representing 4% (537/12 233) of all diabetes cases diagnosed after age 30 years. The clinical characteristics of the group diagnosed with type 1 diabetes when aged 31-60 years were similar to the clinical characteristics of the group diagnosed with type 1 diabetes when aged 30 years or younger. For individuals diagnosed with diabetes when aged 31-60 years, the clinical characteristics of type 1 diabetes differed from those of type 2 diabetes: they had a lower BMI (27·4 kg/m2 [95% CI 26·7-28·0] vs 32·4 kg/m2 [32·2-32·5]; p<0·0001), were more likely to use insulin in the first year after diagnosis (89% [476/537] vs 6% [648/11 696]; p<0·0001), and were more likely to have diabetic ketoacidosis (11% [61/537] vs 0·3% [30/11 696]; p<0·0001).

INTERPRETATION

Genetic susceptibility to type 1 diabetes results in non-obesity-related, insulin-dependent diabetes, which presents throughout the first six decades of life. Our results highlight the difficulty of identifying type 1 diabetes after age 30 years because of the increasing background prevalence of type 2 diabetes. Failure to diagnose late-onset type 1 diabetes can have serious consequences because these patients rapidly develop insulin dependency.

FUNDING

Wellcome Trust and Diabetes UK.

Cohort profile for the MASTERMIND study: using the Clinical Practice Research Datalink (CPRD) to investigate stratification of response to treatment in patients with type 2 diabetes

PURPOSE

This is a retrospective cohort study using observational data from anonymised primary care records. We identify and extract all patients with type 2 diabetes and associated clinical data from the Clinical Practice Research Datalink (CPRD) to inform models of disease progression and stratification of treatment.

PARTICIPANTS

Data were extracted from CPRD on 8 August 2016. The initial data set contained all patients (n=313 485) in the database who had received a type 2 diabetes medication. Criteria were applied to identify and exclude those with type 1 diabetes, polycystic ovarian syndrome or other forms of diabetes (n=40 204), and for data quality control (n=12). We identified 251 338 patients for inclusion in future analyses of diabetes progression and treatment response.

FINDINGS TO DATE

For 6-month response to treatment, measured by change in glycated haemoglobin (HbA1c), we have 91 765 patients with 119 785 treatment response episodes. The greatest impact on reduction of HbA1c occurs with first-line and second-line treatments, metformin and sulfonylurea. Patients moving to third-line treatments tend to have greater weights and higher body mass index. We have investigated the impact of non-adherence to commonly used glucose-lowering medications on HbA1c. For baseline-adjusted HbA1c change over 1 year, non-adherent patients had lower HbA1c reductions than adherent patients, with mean and 95% CI of -4.4 (-4.7 to -4.0) mmol/mol (-0.40 (-0.43 to -0.37) %).

FUTURE PLANS

Findings from studies using these data will help inform future treatment plans and guidelines. Additional data are added with updates from CPRD. This will increase the numbers of patients on newer medications and add more data on those already receiving treatment. There are several ongoing studies investigating different hypotheses regarding differential response to treatment and progression of diabetes. For side effects, links to Hospital Episode Statistics data, where severe events such as hypoglycaemia will be recorded, will also be explored.