REVEL Is Better at Predicting Pathogenicity of Loss-of-Function than Gain-of-Function Variants

This is the final version. Available on open access from Hindawi via the DOI in this recordData Availability: The list of variants used in this study are included in Supplementary Table 1.In silico predictive tools can help determine the pathogenicity of variants. The 2015 American College of Medical Genetics and Genomics (ACMG) guidelines recommended that scores from these tools can be used as supporting evidence of pathogenicity. A subsequent publication by the ClinGen Sequence Variant Interpretation Working Group suggested that high scores from some tools were sufficiently predictive to be used as moderate or strong evidence of pathogenicity. REVEL is a widely used metapredictor that uses the scores of 13 individual in silico tools to calculate the pathogenicity of missense variants. Its ability to predict missense pathogenicity has been assessed extensively; however, no study has previously tested whether its performance is affected by whether the missense variant acts via a loss-of-function (LoF) or gain-of-function (GoF) mechanism. We used a highly curated dataset of 66 confirmed LoF and 65 confirmed GoF variants to evaluate whether this affected the performance of REVEL. 98% of LoF and 100% of GoF variants met the author-recommended REVEL threshold of 0.5 for pathogenicity, while 89% of LoF and 88% of GoF variants exceeded the 0.75 threshold. However, while 55% of LoF variants met the threshold recommended for a REVEL score to count as strong evidence of pathogenicity from the ACMG guidelines (0.932), only 35% of GoF variants met this threshold (). GoF variants are therefore less likely to receive the highest REVEL scores which would enable the REVEL score to be used as strong evidence of pathogenicity. This has implications for classification with the ACMG guidelines as GoF variants are less likely to meet the criteria for pathogenicity. P = 0.0352 ). GoF variants are therefore less likely to receive the highest REVEL scores which would enable the REVEL score to be used as strong evidence of pathogenicity. This has implications for classification with the ACMG guidelines as GoF variants are less likely to meet the criteria for pathogenicity.Wellcome TrustResearch EnglandNational Institute for Health and Care Research (NIHR

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

This is the uncorrected proof. The final version is available on open access from Public Library of Science via the DOI in this recordData Availability: The SavvyCNV tool and the code used to run the benchmarking comparisons are freely available on github. The tool is available at https://github.com/rdemolgen/SavvySuite. The code used to run the benchmarking comparisons is available at: https://github.com/exeter-matthew-wakeling/SavvyCNV_benchmarking. Our study uses the ICR96 data set for benchmarking, which is publicly available and can be accessed through the European-Genome phenome Archive (EGA) under the accession number EGAS00001002428. The dataset of 2591 samples referred to the molecular genetics department at the Royal Devon and Exeter Hospital for genetic testing cannot be shared due to patient confidentiality issues, as the genotype data could be used to identify individuals and so cannot be made openly available. Requests for access to the anonymised data by researchers will be considered following an application to the Genetic Beta Cell Research Bank (https://www.diabetesgenes.org/current-research/genetic-beta-cell-research-bank/) with proposals reviewed by the Genetic Data Access Committee.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.Medical Research Council (MRC)Research EnglandDiabetes U

Evaluation of Evidence for Pathogenicity Demonstrates that BLK, KLF11 and PAX4 Should not be Included in Diagnostic Testing for MODY

This is the author accepted manuscript. The final version is available from the American Diabetes Association via the DOI in this recordData and Resource Availability: UK Biobank data is accessible via application: https://www.ukbiobank.ac.uk/enable-yourresearch. GnomAD data is publically available: https://gnomad.broadinstitute.org/. The MODY cohort data is not publicly available due the limitations of the current ethics and to protect patient confidentiality but is available from the corresponding authors on reasonable request. No applicable resources were generated or analyzed during the current study.Maturity Onset Diabetes of the Young (MODY) is an autosomal dominant form of monogenic diabetes, reported to be caused by variants in 16 genes. Concern has been raised about whether variants in BLK (MODY11), KLF11 (MODY7) and PAX4 (MODY9) cause MODY. We examined variant-level genetic evidence (co-segregation with diabetes and frequency in population) for published putative pathogenic variants in these genes and used burden testing to test gene-level evidence in a MODY cohort (n=1227) compared to population control (UK Biobank, n=185,898). For comparison we analysed well-established causes of MODY, HNF1A and HNF4A. The published variants in BLK, KLF11 and PAX4 showed poor co-segregation with diabetes (combined LOD scores ≤1.2), compared to HNF1A and HNF4A (LOD scores >9), and are all too common to cause MODY (minor allele frequency >4.95x10-5). Ultra-rare missense and protein-truncating variants (PTVs) were not enriched in a MODY cohort compared to the UK Biobank (PTVs P>0.05, missense P>0.1 for all three genes) while HNF1A and HNF4A were enriched (P<10-6). Sensitivity analyses using different population cohorts supported our results. Variant and gene-level genetic evidence does not support BLK, KLF11 or PAX4 as causes of MODY. They should not be included in MODY diagnostic genetic testing.Medical Research Council (MRC)Diabetes UKResearch EnglandWellcome Trus

The Role of ONECUT1 Variants in Monogenic and Type 2 Diabetes Mellitus

This is the author accepted manuscript. The final version is available from the American Diabetes Association via the DOI in this recordData and Resource Availability: The monogenic diabetes patient data used in this study are available from the corresponding author upon reasonable request. The type 2 diabetes data used in this study are available via UK Biobank (https://www.ukbiobank.ac.uk/), subject to necessary approvals.ONECUT1 (also known as HNF6) is a transcription factor involved in pancreatic development and beta-cell function. Recently, biallelic variants in ONECUT1 were reported as a cause of neonatal diabetes mellitus (NDM) in 2 subjects and missense monoallelic variants were associated with type 2 diabetes and possibly maturity-onset diabetes of the young (MODY). Here we examine the role of ONECUT1 variants in NDM, MODY and Type 2 diabetes in large international cohorts of subjects with monogenic diabetes and >400,000 subjects from UK Biobank. We identified a biallelic frameshift ONECUT1 variant as the cause of NDM in one individual. However, we found no enrichment of missense or null ONECUT1 variants among 484 individuals clinically suspected of MODY, in whom all known genes had been excluded. Finally, using a rare variant burden test in the UK Biobank European cohort, we identified a significant association between heterozygous ONECUT1 null variants and type 2 diabetes (P=0.006) but did not find association between missense variants and type 2 diabetes. Our results confirm biallelic ONECUT1 variants as a cause of NDM and highlight monoallelic null variants as a risk factor for type 2 diabetes. These findings confirm the critical role of ONECUT1 in human beta-cell function.Diabetes UKMedical Research Council (MRC)Wellcome TrustNational Institute for Health Research (NIHR

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

Version 2; peer review: 3 approved. Available from F1000 Research via the DOI in this recordBackground: 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.Wellcome TrustRoyal Societ

Type 1 diabetes genetic risk score discriminates between monogenic and Type 1 diabetes in children diagnosed at the age of < 5 years in the Iranian population

This is the final version. Available on open access from Wiley via the DOI in this recordAim To examine the extent to which discriminatory testing using antibodies and Type 1 diabetes genetic risk score, validated in European populations, is applicable in a non‐European population. Methods We recruited 127 unrelated children with diabetes diagnosed between 9 months and 5 years from two centres in Iran. All children underwent targeted next‐generation sequencing of 35 monogenic diabetes genes. We measured three islet autoantibodies (islet antigen 2, glutamic acid decarboxylase and zinc transporter 8) and generated a Type 1 diabetes genetic risk score in all children. Results We identified six children with monogenic diabetes, including four novel mutations: homozygous mutations in WFS1 (n=3), SLC19A2 and SLC29A3, and a heterozygous mutation in GCK. All clinical features were similar in children with monogenic diabetes (n=6) and in the rest of the cohort (n=121). The Type 1 diabetes genetic risk score discriminated children with monogenic from Type 1 diabetes [area under the receiver‐operating characteristic curve 0.90 (95% CI 0.83–0.97)]. All children with monogenic diabetes were autoantibody‐negative. In children with no mutation, 59 were positive to glutamic acid decarboxylase, 39 to islet antigen 2 and 31 to zinc transporter 8. Measuring zinc transporter 8 increased the number of autoantibody‐positive individuals by eight. Conclusions The present study provides the first evidence that Type 1 diabetes genetic risk score can be used to distinguish monogenic from Type 1 diabetes in an Iranian population with a large number of consanguineous unions. This test can be used to identify children with a higher probability of having monogenic diabetes who could then undergo genetic testing. Identification of these individuals would reduce the cost of treatment and improve the management of their clinical course.Wellcome TrustDiabetes U

Chromosome 20p11.2 deletions cause congenital hyperinsulinism via the loss of FOXA2 or its regulatory elements

This is the final version. Available on open access from Springer Nature via the DOI in this recordData availability: All non-clinical data analyzed during this study are included in this article (and its Supplementary Information). The 20p11.2 variants reported in this study were uploaded to ClinVar (SUB14235415). Clinical and genotype data can be used to identify individuals and are therefore available only through collaboration to experienced teams working on approved studies examining the mechanisms, cause, diagnosis and treatment of diabetes and other beta cell disorders. Requests for collaboration will be considered by a steering committee following an application to the Genetic Beta Cell Research Bank (https://www.diabetesgenes.org/current-research/genetic-beta-cell-research-bank/). Contact by email should be directed to S. Flanagan ([email protected]). All requests for access to data will be responded to within 14 d. Accession codes and DOI numbers for all ChIP-seq, ATAC-seq, RNA-seq and scRNA-seq datasets are provided in Supplementary Table 2. We used the Genome Reference Consortium Human Build 37 (GRCh37) to annotate genetic data (accession number GCF_000001405.13). Details of this assembly are provided at https://www.ncbi.nlm.nih.gov/assembly/GCF_000001405.13/.Persistent congenital hyperinsulinism (HI) is a rare genetically heterogeneous condition characterised by dysregulated insulin secretion leading to life-threatening hypoglycaemia. For up to 50% of affected individuals screening of the known HI genes does not identify a disease-causing variant. Large deletions have previously been used to identify novel regulatory regions causing HI. Here, we used genome sequencing to search for novel large (>1 Mb) deletions in 180 probands with HI of unknown cause and replicated our findings in a large cohort of 883 genetically unsolved individuals with HI using off-target copy number variant calling from targeted gene panels. We identified overlapping heterozygous deletions in five individuals (range 3-8 Mb) spanning chromosome 20p11.2. The pancreatic beta-cell transcription factor gene, FOXA2, a known cause of HI was deleted in two of the five individuals. In the remaining three, we found a minimal deleted region of 2.4 Mb adjacent to FOXA2 that encompasses multiple non-coding regulatory elements that are in conformational contact with FOXA2. Our data suggests that the deletions in these three children may cause disease through the dysregulation of FOXA2 expression. These findings provide new insights into the regulation of FOXA2 in the beta-cell and confirm an aetiological role for chromosome 20p11.2 deletions in syndromic HI.Wellcome Trus

Non-coding variants disrupting a tissue-specific regulatory element in HK1 cause congenital hyperinsulinism.

This is the author accepted manuscript. The final version is available from Nature Research via the DOI in this recordData availability statement: All non‐clinical data analysed during this study are included in this published article (and its supplementary information files). Clinical and genotype data is available only through collaboration as this can be used to identify individuals and so cannot be made openly available. Requests for collaboration will be considered following an application to the Genetic Beta Cell Research Bank (https://www.diabetesgenes.org/current‐research/genetic‐ beta‐cell‐research‐bank/). Contact by email should be directed to the Corresponding author.Code availability statement: All code and software versions used specified in Methods.Gene expression is tightly regulated, with many genes exhibiting cell-specific silencing when their protein product would disrupt normal cellular function1. This silencing is largely controlled by non-coding elements, and their disruption might cause human disease2. We performed gene-agnostic screening of the non-coding regions to discover new molecular causes of congenital hyperinsulinism. This identified 14 non-coding de novo variants affecting a 42-bp conserved region encompassed by a regulatory element in intron 2 of the hexokinase 1 gene (HK1). HK1 is widely expressed across all tissues except in the liver and pancreatic beta cells and is thus termed a 'disallowed gene' in these specific tissues. We demonstrated that the variants result in a loss of repression of HK1 in pancreatic beta cells, thereby causing insulin secretion and congenital hyperinsulinism. Using epigenomic data accessed from public repositories, we demonstrated that these variants reside within a regulatory region that we determine to be critical for cell-specific silencing. Importantly, this has revealed a disease mechanism for non-coding variants that cause inappropriate expression of a disallowed gene.Wellcome Trus

Viral Bcl-2-Mediated Evasion of Autophagy Aids Chronic Infection of γHerpesvirus 68

γ-herpesviruses (γHVs) have developed an interaction with their hosts wherein they establish a life-long persistent infection and are associated with the onset of various malignancies. One critical virulence factor involved in the persistency of murine γ-herpesvirus 68 (γHV68) is the viral homolog of the Bcl-2 protein (vBcl-2), which has been implicated to counteract both host apoptotic responses and autophagy pathway. However, the relative significance of the two activities of vBcl-2 in viral persistent infection has yet to be elucidated. Here, by characterizing a series of loss-of-function mutants of vBcl-2, we have distinguished the vBcl-2-mediated antagonism of autophagy from the vBcl-2-mediated inhibition of apoptosis in vitro and in vivo. A mutant γHV68 virus lacking the anti-autophagic activity of vBcl-2 demonstrates an impaired ability to maintain chronic infections in mice, whereas a mutant virus lacking the anti-apoptotic activity of vBcl-2 establishes chronic infections as efficiently as the wild-type virus but displays a compromised ability for ex vivo reactivation. Thus, the vBcl-2-mediated antagonism of host autophagy constitutes a novel mechanism by which γHVs confer persistent infections, further underscoring the importance of autophagy as a critical host determinant in the in vivo latency of γ-herpesviruses