Penetrance and expressivity of mitochondrial variants in a large clinically unselected population.

Whole genome sequencing (WGS) from large clinically unselected cohorts provides a unique opportunity to assess the penetrance and expressivity of rare and/or known pathogenic mitochondrial variants in population. Using WGS from 179 862 clinically unselected individuals from the UK Biobank, we performed extensive single and rare variant aggregation association analyses of 15 881 mtDNA variants and 73 known pathogenic variants with 15 mitochondrial disease-relevant phenotypes. We identified 12 homoplasmic and one heteroplasmic variant (m.3243A>G) with genome-wide significant associations in our clinically unselected cohort. Heteroplasmic m.3243A>G (MAF = 0.0002, a known pathogenic variant) was associated with diabetes, deafness and heart failure and 12 homoplasmic variants increased aspartate aminotransferase levels including three low-frequency variants (MAF ~0.002 and beta~0.3 SD). Most pathogenic mitochondrial disease variants (n = 66/74) were rare in the population (<1:9000). Aggregated or single variant analysis of pathogenic variants showed low penetrance in unselected settings for the relevant phenotypes, except m.3243A>G. Multi-system disease risk and penetrance of diabetes, deafness and heart failure greatly increased with m.3243A>G level ≥ 10%. The odds ratio of these traits increased from 5.61, 12.3 and 10.1 to 25.1, 55.0 and 39.5, respectively. Diabetes risk with m.3243A>G was further influenced by type 2 diabetes genetic risk. Our study of mitochondrial variation in a large-unselected population identified novel associations and demonstrated that pathogenic mitochondrial variants have lower penetrance in clinically unselected settings. m.3243A>G was an exception at higher heteroplasmy showing a significant impact on health making it a good candidate for incidental reporting.

Identification and analysis of individuals who deviate from their genetically-predicted phenotype.

Findings from genome-wide association studies have facilitated the generation of genetic predictors for many common human phenotypes. Stratifying individuals misaligned to a genetic predictor based on common variants may be important for follow-up studies that aim to identify alternative causal factors. Using genome-wide imputed genetic data, we aimed to classify 158,951 unrelated individuals from the UK Biobank as either concordant or deviating from two well-measured phenotypes. We first applied our methods to standing height: our primary analysis classified 244 individuals (0.15%) as misaligned to their genetically predicted height. We show that these individuals are enriched for self-reporting being shorter or taller than average at age 10, diagnosed congenital malformations, and rare loss-of-function variants in genes previously catalogued as causal for growth disorders. Secondly, we apply our methods to LDL cholesterol (LDL-C). We classified 156 (0.12%) individuals as misaligned to their genetically predicted LDL-C and show that these individuals were enriched for both clinically actionable cardiovascular risk factors and rare genetic variants in genes previously shown to be involved in metabolic processes. Individuals whose LDL-C was higher than expected based on the genetic predictor were also at higher risk of developing coronary artery disease and type-two diabetes, even after adjustment for measured LDL-C, BMI and age, suggesting upward deviation from genetically predicted LDL-C is indicative of generally poor health. Our results remained broadly consistent when performing sensitivity analysis based on a variety of parametric and non-parametric methods to define individuals deviating from polygenic expectation. Our analyses demonstrate the potential importance of quantitatively identifying individuals for further follow-up based on deviation from genetic predictions.

Rare genetic variants in genes and loci linked to dominant monogenic developmental disorders cause milder related phenotypes in the general population.

Many rare monogenic diseases are known to be caused by deleterious variants in thousands of genes, however the same variants can also be found in people without the associated clinical phenotypes. The penetrance of these monogenic variants is generally unknown in the wider population, as they are typically identified in small clinical cohorts of affected individuals and families with highly penetrant variants. Here, we investigated the phenotypic effect of rare, potentially deleterious variants in genes and loci where similar variants are known to cause monogenic developmental disorders (DDs) in a large population cohort. We used UK Biobank to investigate phenotypes associated with rare protein-truncating and missense variants in 599 monoallelic DDG2P genes by using whole-exome-sequencing data from ∼200,000 individuals and rare copy-number variants overlapping known DD loci by using SNP-array data from ∼500,000 individuals. We found that individuals with these likely deleterious variants had a mild DD-related phenotype, including lower fluid intelligence, slower reaction times, lower numeric memory scores, and longer pairs matching times compared to the rest of the UK Biobank cohort. They were also shorter, had a higher BMI, and had significant socioeconomic disadvantages: they were less likely to be employed or be able to work and had a lower income and higher deprivation index. Our findings suggest that many genes routinely tested within pediatric genetics have deleterious variants with intermediate penetrance that may cause lifelong sub-clinical phenotypes in the general adult population.

Reduced penetrance of MODY-associated HNF1A/HNF4A variants but not GCK variants in clinically unselected cohorts.

The true prevalence and penetrance of monogenic disease variants are often not known because of clinical-referral ascertainment bias. We comprehensively assess the penetrance and prevalence of pathogenic variants in HNF1A, HNF4A, and GCK that account for >80% of monogenic diabetes. We analyzed clinical and genetic data from 1,742 clinically referred probands, 2,194 family members, clinically unselected individuals from a US health system-based cohort (n = 132,194), and a UK population-based cohort (n = 198,748). We show that one in 1,500 individuals harbor a pathogenic variant in one of these genes. The penetrance of diabetes for HNF1A and HNF4A pathogenic variants was substantially lower in the clinically unselected individuals compared to clinically referred probands and was dependent on the setting (32% in the population, 49% in the health system cohort, 86% in a family member, and 98% in probands for HNF1A). The relative risk of diabetes was similar across the clinically unselected cohorts highlighting the role of environment/other genetic factors. Surprisingly, the penetrance of pathogenic GCK variants was similar across all cohorts (89%-97%). We highlight that pathogenic variants in HNF1A, HNF4A, and GCK are not ultra-rare in the population. For HNF1A and HNF4A, we need to tailor genetic interpretation and counseling based on the setting in which a pathogenic monogenic variant was identified. GCK is an exception with near-complete penetrance in all settings. This along with the clinical implication of diagnosis makes it an excellent candidate for the American College of Medical Genetics secondary gene list.

The impact of Mendelian sleep and circadian genetic variants in a population setting.

Rare variants in ten genes have been reported to cause Mendelian sleep conditions characterised by extreme sleep duration or timing. These include familial natural short sleep (ADRB1, DEC2/BHLHE41, GRM1 and NPSR1), advanced sleep phase (PER2, PER3, CRY2, CSNK1D and TIMELESS) and delayed sleep phase (CRY1). The association of variants in these genes with extreme sleep conditions were usually based on clinically ascertained families, and their effects when identified in the population are unknown. We aimed to determine the effects of these variants on sleep traits in large population-based cohorts. We performed genetic association analysis of variants previously reported to be causal for Mendelian sleep and circadian conditions. Analyses were performed using 191,929 individuals with data on sleep and whole-exome or genome-sequence data from 4 population-based studies: UK Biobank, FINRISK, Health-2000-2001, and the Multi-Ethnic Study of Atherosclerosis (MESA). We identified sleep disorders from self-report, hospital and primary care data. We estimated sleep duration and timing measures from self-report and accelerometery data. We identified carriers for 10 out of 12 previously reported pathogenic variants for 8 of the 10 genes. They ranged in frequency from 1 individual with the variant in CSNK1D to 1,574 individuals with a reported variant in the PER3 gene in the UK Biobank. No carriers for variants reported in NPSR1 or PER2 were identified. We found no association between variants analyzed and extreme sleep or circadian phenotypes. Using sleep timing as a proxy measure for sleep phase, only PER3 and CRY1 variants demonstrated association with earlier and later sleep timing, respectively; however, the magnitude of effect was smaller than previously reported (sleep midpoint ~7 mins earlier and ~5 mins later, respectively). We also performed burden tests of protein truncating (PTVs) or rare missense variants for the 10 genes. Only PTVs in PER2 and PER3 were associated with a relevant trait (for example, 64 individuals with a PTV in PER2 had an odds ratio of 4.4 for being "definitely a morning person", P = 4x10-8; and had a 57-minute earlier midpoint sleep, P = 5x10-7). Our results indicate that previously reported variants for Mendelian sleep and circadian conditions are often not highly penetrant when ascertained incidentally from the general population.

Associations between diabetes-related genetic risk scores and residual beta cell function in type 1 diabetes: the GUTDM1 study.

AIMS/HYPOTHESIS: Use of genetic risk scores (GRS) may help to distinguish between type 1 diabetes and type 2 diabetes, but less is known about whether GRS are associated with disease severity or progression after diagnosis. Therefore, we tested whether GRS are associated with residual beta cell function and glycaemic control in individuals with type 1 diabetes. METHODS: Immunochip arrays and TOPMed were used to genotype a cross-sectional cohort (n=479, age 41.7 ± 14.9 years, duration of diabetes 16.0 years [IQR 6.0-29.0], HbA1c 55.6 ± 12.2 mmol/mol). Several GRS, which were originally developed to assess genetic risk of type 1 diabetes (GRS-1, GRS-2) and type 2 diabetes (GRS-T2D), were calculated. GRS-C1 and GRS-C2 were based on SNPs that have previously been shown to be associated with residual beta cell function. Regression models were used to investigate the association between GRS and residual beta cell function, assessed using the urinary C-peptide/creatinine ratio, and the association between GRS and continuous glucose monitor metrics. RESULTS: Higher GRS-1 and higher GRS-2 both showed a significant association with undetectable UCPCR (OR 0.78; 95% CI 0.69, 0.89 and OR 0.84: 95% CI 0.75, 0.93, respectively), which were attenuated after correction for sex and age of onset (GRS-2) and disease duration (GRS-1). Higher GRS-C2 was associated with detectable urinary C-peptide/creatinine ratio (≥0.01 nmol/mmol) after correction for sex and age of onset (OR 6.95; 95% CI 1.19, 40.75). A higher GRS-T2D was associated with less time below range (TBR) (OR for TBR<4% 1.41; 95% CI 1.01 to 1.96) and lower glucose coefficient of variance (ß -1.53; 95% CI -2.76, -0.29). CONCLUSIONS/INTERPRETATION: Diabetes-related GRS are associated with residual beta cell function in individuals with type 1 diabetes. These findings suggest some genetic contribution to preservation of beta cell function.

Diabetes Study of Children of Diverse Ethnicity and Race: Study design.

AIMS: Determining diabetes type in children has become increasingly difficult due to an overlap in typical characteristics between type 1 diabetes (T1D) and type 2 diabetes (T2D). The Diabetes Study in Children of Diverse Ethnicity and Race (DISCOVER) programme is a National Institutes of Health (NIH)-supported multicenter, prospective, observational study that enrols children and adolescents with non-secondary diabetes. The primary aim of the study was to develop improved models to differentiate between T1D and T2D in diverse youth. MATERIALS AND METHODS: The proposed models will evaluate the utility of three existing T1D genetic risk scores in combination with data on islet autoantibodies and other parameters typically available at the time of diabetes onset. Low non-fasting serum C-peptide (<0.6 nmol/L) between 3 and 10 years after diabetes diagnosis will be considered a biomarker for T1D as it reflects the loss of insulin secretion ability. Participating centres are enrolling youth (<19 years old) either with established diabetes (duration 3-10 years) for a cross-sectional evaluation or with recent onset diabetes (duration 3 weeks-15 months) for the longitudinal observation with annual visits for 3 years. Cross-sectional data will be used to develop models. Longitudinal data will be used to externally validate the best-fitting model. RESULTS: The results are expected to improve the ability to classify diabetes type in a large and growing subset of children who have an unclear form of diabetes at diagnosis. CONCLUSIONS: Accurate and timely classification of diabetes type will help establish the correct clinical management early in the course of the disease.

A founder UMOD variant is a common cause of hereditary nephropathy in the British population.

BACKGROUND: Monogenic disorders are estimated to account for 10%-12% of patients with kidney failure. We report the unexpected finding of an unusual uromodulin (UMOD) variant in multiple pedigrees within the British population and demonstrate a shared haplotype indicative of an ancestral variant. METHODS: Probands from 12 apparently unrelated pedigrees with a family history of kidney failure within a geographically contiguous UK region were shown to be heterozygous for a pathogenic variant of UMOD c.278_289delTCTGCCCCGAAG insCCGCCTCCT. RESULTS: A total of 88 clinically affected individuals were identified, all born in the UK and of white British ethnicity. 20 other individuals with the variant were identified in the UK 100,000 Genomes (100K) Project and 9 from UK Biobank (UKBB). A common extended haplotype was present in 5 of the UKBB individuals who underwent genome sequencing which was only present in <1 in 5000 of UKBB controls. Significantly, rare variants (<1 in 250 general population) identified within 1 Mb of the UMOD variant by genome sequencing were detected in all of the 100K individuals, indicative of an extended shared haplotype. CONCLUSION: Our data confirm a likely founder UMOD variant with a wide geographical distribution within the UK. It should be suspected in cases of unexplained familial nephropathy presenting in patients of white British ancestry.

Prevalence of Fabry disease-causing variants in the UK Biobank.

BACKGROUND: Fabry disease is an X-linked lysosomal storage disorder resulting from deficiency of the alpha-galactosidase A enzyme leading to accumulation of globotriaosylceramide in multiple organ sites with prominent cardiovascular and renal involvement. Global prevalence estimates of Fabry disease based on clinical ascertainment range from 1 in 40 000 to 1 in 170 000. We aimed to determine the prevalence of Fabry disease-causing variants in UK Biobank. METHODS: We sought GLA gene variants in exome sequencing data from 200 643 individuals from UK Biobank. We used ACMG/AMP guidelines (American College of Medical Genetics/Association for Molecular Pathology) to classify pathogenicity and compared baseline biomarker data, hospital ICD-10 (International Classification of Diseases version-10) codes, general practitioner records and self-reported health data with those without pathogenic variants. RESULTS: We identified 81 GLA coding variants. We identified eight likely pathogenic variants on the basis of being rare (<1/10 000 individuals) and either previously reported to cause Fabry disease, or being protein-truncating variants. Thirty-six individuals carried one of these variants. In the UK Biobank, the prevalence of likely pathogenic Fabry disease-causing variants is 1/5732 for late-onset disease-causing variants and 1/200 643 for variants causing classic Fabry disease. CONCLUSION: Fabry disease-causing GLA variants are more prevalent in an unselected population sample than the reported prevalence of Fabry disease. These are overwhelmingly variants associated with later onset. It is possible the prevalence of later-onset Fabry disease exceeds current estimates.

Recurrent 17q12 microduplications contribute to renal disease but not diabetes.

BACKGROUND: 17q12 microdeletion and microduplication syndromes present as overlapping, multisystem disorders. We assessed the disease phenotypes of individuals with 17q12 CNV in a population-based cohort. METHODS: We investigated 17q12 CNV using microarray data from 450 993 individuals in the UK Biobank and calculated disease status associations for diabetes, liver and renal function, neurological and psychiatric traits. RESULTS: We identified 11 17q12 microdeletions and 106 microduplications. Microdeletions were strongly associated with diabetes (p=2×10-7) but microduplications were not. Estimated glomerular filtration rate (eGFR mL/min/1.73 m2) was consistently lower in individuals with microdeletions (p=3×10-12) and microduplications (p=6×10-25). Similarly, eGFR <60, including end-stage renal disease, was associated with microdeletions (p=2×10-9, p<0.003) and microduplications (p=1×10-9, p=0.009), respectively, highlighting sometimes substantially reduced renal function in each. Microduplications were associated with decreased fluid intelligence (p=3×10-4). SNP association analysis in the 17q12 region implicated changes to HNF1B as causing decreased eGFR (NC_000017.11:g.37741642T>G, rs12601991, p=4×10-21) and diabetes (NC_000017.11:g.37741165C>T, rs7501939, p=6×10-17). A second locus within the region was also associated with fluid intelligence (NC_000017.11:g.36593168T>C, rs1005552, p=6×10-9) and decreased eGFR (NC_000017.11:g.36558947T>C, rs12150665, p=4×10-15). CONCLUSION: We demonstrate 17q12 microdeletions but not microduplications are associated with diabetes in a population-based cohort, likely caused by HNF1B haploinsufficiency. We show that both 17q12 microdeletions and microduplications are associated with renal disease, and multiple genes within the region likely contribute to renal and neurocognitive phenotypes.

A genome-wide association study identifies 5 loci associated with frozen shoulder and implicates diabetes as a causal risk factor.

Frozen shoulder is a painful condition that often requires surgery and affects up to 5% of individuals aged 40-60 years. Little is known about the causes of the condition, but diabetes is a strong risk factor. To begin to understand the biological mechanisms involved, we aimed to identify genetic variants associated with frozen shoulder and to use Mendelian randomization to test the causal role of diabetes. We performed a genome-wide association study (GWAS) of frozen shoulder in the UK Biobank using data from 10,104 cases identified from inpatient, surgical and primary care codes. We used data from FinnGen for replication and meta-analysis. We used one-sample and two-sample Mendelian randomization approaches to test for a causal association of diabetes with frozen shoulder. We identified five genome-wide significant loci. The most significant locus (lead SNP rs28971325; OR = 1.20, [95% CI: 1.16-1.24], p = 5x10-29) contained WNT7B. This variant was also associated with Dupuytren's disease (OR = 2.31 [2.24, 2.39], p<1x10-300) as were a further two of the frozen shoulder associated variants. The Mendelian randomization results provided evidence that type 1 diabetes is a causal risk factor for frozen shoulder (OR = 1.03 [1.02-1.05], p = 3x10-6). There was no evidence that obesity was causally associated with frozen shoulder, suggesting that diabetes influences risk of the condition through glycemic rather than mechanical effects. We have identified genetic loci associated with frozen shoulder. There is a large overlap with Dupuytren's disease associated loci. Diabetes is a likely causal risk factor. Our results provide evidence of biological mechanisms involved in this common painful condition.

Utility of genetic risk scores in type 1 diabetes.

Iterative advances in understanding of the genetics of type 1 diabetes have identified >70 genetic regions associated with risk of the disease, including strong associations across the HLA class II region that account for >50% of heritability. The increased availability of genetic data combined with the decreased costs of generating these data, have facilitated the development of polygenic scores that aggregate risk variants from associated loci into a single number: either a genetic risk score (GRS) or a polygenic risk score (PRS). PRSs incorporate the risk of many possibly correlated variants from across the genome, even if they do not reach genome-wide significance, whereas GRSs estimate the cumulative contribution of a smaller subset of genetic variants that reach genome-wide significance. Type 1 diabetes GRSs have utility in diabetes classification, aiding discrimination between type 1 diabetes, type 2 diabetes and MODY. Type 1 diabetes GRSs are also being used in newborn screening studies to identify infants at risk of future presentation of the disease. Most early studies of type 1 diabetes genetics have been conducted in European ancestry populations, but, to develop accurate GRSs across diverse ancestries, large case-control cohorts from non-European populations are still needed. The current barriers to GRS implementation within healthcare are mainly related to a lack of guidance and knowledge on integration with other biomarkers and clinical variables. Once these limitations are addressed, there is huge potential for 'test and treat' approaches to be used to tailor care for individuals with type 1 diabetes.

The relationship between islet autoantibody status and the genetic risk of type 1 diabetes in adult-onset type 1 diabetes.

AIMS/HYPOTHESIS: The reason for the observed lower rate of islet autoantibody positivity in clinician-diagnosed adult-onset vs childhood-onset type 1 diabetes is not known. We aimed to explore this by assessing the genetic risk of type 1 diabetes in autoantibody-negative and -positive children and adults. METHODS: We analysed GAD autoantibodies, insulinoma-2 antigen autoantibodies and zinc transporter-8 autoantibodies (ZnT8A) and measured type 1 diabetes genetic risk by genotyping 30 type 1 diabetes-associated variants at diagnosis in 1814 individuals with clinician-diagnosed type 1 diabetes (1112 adult-onset, 702 childhood-onset). We compared the overall type 1 diabetes genetic risk score (T1DGRS) and non-HLA and HLA (DR3-DQ2, DR4-DQ8 and DR15-DQ6) components with autoantibody status in those with adult-onset and childhood-onset diabetes. We also measured the T1DGRS in 1924 individuals with type 2 diabetes from the Wellcome Trust Case Control Consortium to represent non-autoimmune diabetes control participants. RESULTS: The T1DGRS was similar in autoantibody-negative and autoantibody-positive clinician-diagnosed childhood-onset type 1 diabetes (mean [SD] 0.274 [0.034] vs 0.277 [0.026], p=0.4). In contrast, the T1DGRS in autoantibody-negative adult-onset type 1 diabetes was lower than that in autoantibody-positive adult-onset type 1 diabetes (mean [SD] 0.243 [0.036] vs 0.271 [0.026], p<0.0001) but higher than that in type 2 diabetes (mean [SD] 0.229 [0.034], p<0.0001). Autoantibody-negative adults were more likely to have the more protective HLA DR15-DQ6 genotype (15% vs 3%, p<0.0001), were less likely to have the high-risk HLA DR3-DQ2/DR4-DQ8 genotype (6% vs 19%, p<0.0001) and had a lower non-HLA T1DGRS (p<0.0001) than autoantibody-positive adults. In contrast to children, autoantibody-negative adults were more likely to be male (75% vs 59%), had a higher BMI (27 vs 24 kg/m2) and were less likely to have other autoimmune conditions (2% vs 10%) than autoantibody-positive adults (all p<0.0001). In both adults and children, type 1 diabetes genetic risk was unaffected by the number of autoantibodies (p>0.3). These findings, along with the identification of seven misclassified adults with monogenic diabetes among autoantibody-negative adults and the results of a sensitivity analysis with and without measurement of ZnT8A, suggest that the intermediate type 1 diabetes genetic risk in autoantibody-negative adults is more likely to be explained by the inclusion of misclassified non-autoimmune diabetes (estimated to represent 67% of all antibody-negative adults, 95% CI 61%, 73%) than by the presence of unmeasured autoantibodies or by a discrete form of diabetes. When these estimated individuals with non-autoimmune diabetes were adjusted for, the prevalence of autoantibody positivity in adult-onset type 1 diabetes was similar to that in children (93% vs 91%, p=0.4). CONCLUSIONS/INTERPRETATION: The inclusion of non-autoimmune diabetes is the most likely explanation for the observed lower rate of autoantibody positivity in clinician-diagnosed adult-onset type 1 diabetes. Our data support the utility of islet autoantibody measurement in clinician-suspected adult-onset type 1 diabetes in routine clinical practice.

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.

SavvyCNV: Genome-wide CNV calling from off-target reads.

Identifying copy number variants (CNVs) can provide diagnoses to patients and provide important biological insights into human health and disease. Current exome and targeted sequencing approaches cannot detect clinically and biologically-relevant CNVs outside their target area. We present SavvyCNV, a tool which uses off-target read data from exome and targeted sequencing data to call germline CNVs genome-wide. Up to 70% of sequencing reads from exome and targeted sequencing fall outside the targeted regions. We have developed a new tool, SavvyCNV, to exploit this 'free data' to call CNVs across the genome. We benchmarked SavvyCNV against five state-of-the-art CNV callers using truth sets generated from genome sequencing data and Multiplex Ligation-dependent Probe Amplification assays. SavvyCNV called CNVs with high precision and recall, outperforming the five other tools at calling CNVs genome-wide, using off-target or on-target reads from targeted panel and exome sequencing. We then applied SavvyCNV to clinical samples sequenced using a targeted panel and were able to call previously undetected clinically-relevant CNVs, highlighting the utility of this tool within the diagnostic setting. SavvyCNV outperforms existing tools for calling CNVs from off-target reads. It can call CNVs genome-wide from targeted panel and exome data, increasing the utility and diagnostic yield of these tests. SavvyCNV is freely available at https://github.com/rdemolgen/SavvySuite.

Systematic genetic testing for recessively inherited monogenic diabetes: a cross-sectional study in paediatric diabetes clinics.

AIMS/HYPOTHESIS: Current clinical guidelines for childhood-onset monogenic diabetes outside infancy are mainly focused on identifying and testing for dominantly inherited, predominantly MODY genes. There are no systematic studies of the recessively inherited causes of monogenic diabetes that are likely to be more common in populations with high rates of consanguinity. We aimed to determine the contribution of recessive causes of monogenic diabetes in paediatric diabetes clinics and to identify clinical criteria by which to select individuals for recessive monogenic diabetes testing. METHODS: We conducted a cross-sectional study of 1093 children from seven paediatric diabetes clinics across Turkey (a population with high rates of consanguinity). We undertook genetic testing of 50 known dominant and recessive causes of monogenic diabetes for 236 children at low risk of type 1 diabetes. As a comparison, we used monogenic diabetes cases from UK paediatric diabetes clinics (a population with low rates of consanguinity). RESULTS: Thirty-four children in the Turkish cohort had monogenic diabetes, equating to a minimal prevalence of 3.1%, similar to that in the UK cohort (p = 0.40). Forty-one per cent (14/34) had autosomal recessive causes in contrast to 1.6% (2/122) in the UK monogenic diabetes cohort (p < 0.0001). All conventional criteria for identifying monogenic diabetes (parental diabetes, not requiring insulin treatment, HbA1c ≤ 58 mmol/mol [≤7.5%] and a composite clinical probability of MODY >10%) assisted the identification of the dominant (all p ≤ 0.0003) but not recessive cases (all p ≥ 0.2) in Turkey. The presence of certain non-autoimmune extra-pancreatic features greatly assisted the identification of recessive (p < 0.0001, OR 66.9) but not dominant cases. CONCLUSIONS/INTERPRETATION: Recessively inherited mutations are a common cause of monogenic diabetes in populations with high rates of consanguinity. Present MODY-focused genetic testing strategies do not identify affected individuals. To detect all cases of monogenic paediatric diabetes, it is crucial that recessive genes are included in genetic panels and that children are selected for testing if they have certain non-autoimmune extra-pancreatic features in addition to current criteria.

Applying a genetic risk score for prostate cancer to men with lower urinary tract symptoms in primary care to predict prostate cancer diagnosis: a cohort study in the UK Biobank.

BACKGROUND: Prostate cancer is highly heritable, with >250 common variants associated in genome-wide association studies. It commonly presents with non-specific lower urinary tract symptoms that are frequently associated with benign conditions. METHODS: Cohort study using UK Biobank data linked to primary care records. Participants were men with a record showing a general practice consultation for a lower urinary tract symptom. The outcome measure was prostate cancer diagnosis within 2 years of consultation. The predictor was a genetic risk score of 269 genetic variants for prostate cancer. RESULTS: A genetic risk score (GRS) is associated with prostate cancer in symptomatic men (OR per SD increase = 2.12 [1.86-2.41] P = 3.5e-30). An integrated risk model including age and GRS applied to symptomatic men predicted prostate cancer (AUC 0.768 [0.739-0.796]). Prostate cancer incidence was 8.1% (6.7-9.7) in the highest risk quintile. In the lowest quintile, prostate cancer incidence was <1%. CONCLUSIONS: This study is the first to apply GRS in primary care to improve the triage of symptomatic patients. Men with the lowest genetic risk of developing prostate cancer could safely avoid invasive investigation, whilst those identified with the greatest risk could be fast-tracked for further investigation. These results show that a GRS has potential application to improve the diagnostic pathway of symptomatic patients in primary care.

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.

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-14), 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.