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 (cosegregation 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 = 1,227) compared with a control population (UK Biobank [n = 185,898]). For comparison we analyzed well-established causes of MODY, HNF1A, and HNF4A. The published variants in BLK, KLF11, and PAX4 showed poor cosegregation with diabetes (combined logarithm of the odds [LOD] scores ≤1.2), compared with HNF1A and HNF4A (LOD scores >9), and are all too common to cause MODY (minor allele frequency >4.95 × 10-5). Ultra-rare missense and protein-truncating variants (PTV) were not enriched in a MODY cohort compared with the UK Biobank population (PTV P > 0.05, missense P > 0.1 for all three genes) while HNF1A and HNF4A were enriched (P < 10-6). Findings of sensitivity analyses with different population cohorts supported our results. Variant and gene-level genetic evidence does not support BLK, KLF11, or PAX4 as a cause of MODY. They should not be included in MODY diagnostic genetic testing.

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.

HNF4A mutations cause increased birth weight, transient neonatal hypoglycemia, and maturity onset diabetes of the young (MODY). The most frequently reported HNF4A mutation is p.R114W (previously p.R127W), but functional studies have shown inconsistent results; there is a lack of cosegregation in some pedigrees and an unexpectedly high frequency in public variant databases. We confirm that p.R114W is a pathogenic mutation with an odds ratio of 30.4 (95% CI 9.79-125, P = 2 × 10(-21)) for diabetes in our MODY cohort compared with control subjects. p.R114W heterozygotes did not have the increased birth weight of patients with other HNF4A mutations (3,476 g vs. 4,147 g, P = 0.0004), and fewer patients responded to sulfonylurea treatment (48% vs. 73%, P = 0.038). p.R114W has reduced penetrance; only 54% of heterozygotes developed diabetes by age 30 years compared with 71% for other HNF4A mutations. We redefine p.R114W as a pathogenic mutation that causes a distinct clinical subtype of HNF4A MODY with reduced penetrance, reduced sensitivity to sulfonylurea treatment, and no effect on birth weight. This has implications for diabetes treatment, management of pregnancy, and predictive testing of at-risk relatives. The increasing availability of large-scale sequence data is likely to reveal similar examples of rare, low-penetrance MODY mutations.

The Oxford Nanopore Technologies (ONT) MinION is a new sequencing technology that potentially offers read lengths of tens of kilobases (kb) limited only by the length of DNA molecules presented to it. The device has a low capital cost, is by far the most portable DNA sequencer available, and can produce data in real-time. It has numerous prospective applications including improving genome sequence assemblies and resolution of repeat-rich regions. Before such a technology is widely adopted, it is important to assess its performance and limitations in respect of throughput and accuracy. In this study we assessed the performance of the MinION by re-sequencing three bacterial genomes, with very different nucleotide compositions ranging from 28.6% to 70.7%; the high G. +. C strain was underrepresented in the sequencing reads. We estimate the error rate of the MinION (after base calling) to be 38.2%. Mean and median read lengths were 2. kb and 1. kb respectively, while the longest single read was 98. kb. The whole length of a 5. kb rRNA operon was covered by a single read. As the first nanopore-based single molecule sequencer available to researchers, the MinION is an exciting prospect; however, the current error rate limits its ability to compete with existing sequencing technologies, though we do show that MinION sequence reads can enhance contiguity of de novo assembly when used in conjunction with Illumina MiSeq data.

AbstractObjectiveHyperinsulinaemic hypoglycaemia (HH) can occur in isolation or more rarely feature as part of a syndrome. Screening for mutations in the ‘syndromic HH’ genes is guided by phenotype with genetic testing used to confirm the clinical diagnosis. As HH can be the presenting feature of a syndrome it is possible that mutations will be missed as these genes are not routinely screened in all newly diagnosed individuals. We investigated the frequency of pathogenic variants in syndromic genes in individuals with HH who had not been clinically diagnosed with a syndromic disorder at referral for genetic testing.DesignWe used genome sequencing data to assess the prevalence of mutations in syndromic HH genes in an international cohort of patients with HH of unknown genetic cause.MethodsWe undertook genome sequencing in 82 individuals with HH without a clinical diagnosis of a known syndrome at referral for genetic testing. Within this cohort we searched for the genetic aetiologies causing 20 different syndromes where HH had been reported as a feature.ResultsWe identified a pathogenicKMT2Dvariant in a patient with HH diagnosed at birth, confirming a genetic diagnosis of Kabuki syndrome. Clinical data received following the identification of the mutation highlighted additional features consistent with the genetic diagnosis. Pathogenic variants were not identified in the remainder of the cohort.ConclusionsPathogenic variants in the syndromic HH genes are rare but should be considered in newly diagnosed individuals as HH may be the presenting feature.

AbstractMaturity 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 inBLK(MODY11),KLF11(MODY7) andPAX4(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,HNF1AandHNF4A. The published variants inBLK, KLF11andPAX4showed poor co-segregation with diabetes (combined LOD scores ≤1.2), compared toHNF1AandHNF4A(LOD scores &gt;9), and are all too common to cause MODY (minor allele frequency &gt;4.95×10−5). Ultra-rare missense and protein-truncating variants (PTVs) were not enriched in a MODY cohort compared to the UK Biobank (PTVsP&gt;0.05, missenseP&gt;0.1 for all three genes) whileHNF1AandHNF4Awere enriched (P&lt;10−6). Sensitivity analyses using different population cohorts supported our results. Variant and gene-level genetic evidence does not supportBLK, KLF11orPAX4as causes of MODY. They should not be included in MODY diagnostic genetic testing.

AbstractFinding new genetic causes of monogenic diabetes can help to understand development and function of the human pancreas. We aimed to find novel protein–truncating variants causing Maturity–Onset Diabetes of the Young (MODY), a subtype of monogenic diabetes. We used a combination of next–generation sequencing of MODY cases with unknown aetiology along with comparisons to the ExAC database to identify new MODY genes. In the discovery cohort of 36 European patients, we identified two probands with novel RFX6 heterozygous nonsense variants. RFX6 protein truncating variants were enriched in the MODY discovery cohort compared to the European control population within ExAC (odds ratio, OR=131, P=l×l0‐4). We found similar results in non–Finnish European (n=348, OR=43, P=5×l0‐5) and Finnish (n=80, OR=22, P=1×l0‐6) replication cohorts. The overall meta–analysis OR was 34 (P=l×l0‐16). RFX6 heterozygotes had reduced penetrance of diabetes compared to common HNF1A and HNF4A–MODY mutations (27%, 70% and 55% at 25 years of age, respectively). The hyperglycaemia resulted from beta–cell dysfunction and was associated with lower fasting and stimulated gastric inhibitory polypeptide (GIP) levels. Our study demonstrates that heterozygous RFX6 protein truncating variants are associated with MODY with reduced penetrance.

Background
Congenital hyperinsulinism (HI) is characterised by inappropriate insulin secretion despite low blood glucose. Persistent HI is often monogenic, with the majority of cases diagnosed in infancy. Less is known about the contribution of monogenic forms of disease in those presenting in childhood. We investigated the likelihood of finding a genetic cause in childhood-onset HI and explored potential factors leading to the later age at presentation of disease.

Methods
We screened known disease-causing genes in 1848 individuals with HI, referred for genetic testing as part of routine clinical care. Individuals were classified as infancy-onset (when diagnosed with HI

AbstractIdentifying 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 to call CNVs genome-wide. Up to 70% of sequencing reads from exome and targeted sequencing fall outside the targeted regions - SavvyCNV exploits this ‘free data’.We benchmarked SavvyCNV using truth sets generated from genome sequencing data and Multiplex Ligation-dependent Probe Amplification assays. SavvyCNV called CNVs with high precision and recall, outperforming five state-of-the-art CNV callers at calling CNVs genome-wide using off-target or on-target reads from targeted panel and exome sequencing. Furthermore SavvyCNV was able to call previously undetected clinically-relevant CNVs from targeted panel data highlighting the utility of this tool within the diagnostic setting. SavvyCNV is freely available.

AbstractBACKGROUNDAccurate penetrance of monogenic disorders is often unknown due to a phenotype-first approach to genetic testing. Here, we use a genotype-first approach in four large cohorts with different ascertainment contexts to accurately estimate penetrance of the three commonest causes of monogenic diabetes, Maturity Onset Diabetes of the Young (MODY). We contrast HNF1A-MODY / HNF4A-MODY which causes an age-related progressive diabetes and GCK-MODY, which causes life-long mild hyperglycaemia.METHODSWe analysed clinical and genetic sequencing data from four different cohorts: 1742 probands referred for clinical MODY testing; 2194 family members of the MODY probands; 132,194 individuals from an American hospital-based cohort; and 198,748 individuals from a UK population-based cohort.RESULTSAge-related penetrance of diabetes for pathogenic variants in HNF1A and HNF4A was substantially lower in the clinically unselected cohorts compared to clinically referred probands (ranging from 32% to 98% at age 40yrs for HNF1A, and 21% to 99% for HNF4A). The background rate of diabetes, but not clinical features or variant type, explained the reduced penetrance in the unselected cohorts. In contrast, penetrance of mild hyperglycaemia for pathogenic GCK variants was similarly high across cohorts (ranging from 89 to 97%) despite substantial variation in the background rates of diabetes.CONCLUSIONSAscertainment context is crucial when interpreting the consequences of monogenic variants for age-related variably penetrant disorders. This finding has important implications for opportunistic screening during genomic testing.

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 (cosegregation 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 = 1,227) compared with a control population (UK Biobank [n = 185,898]). For comparison we analyzed well-established causes of MODY, HNF1A, and HNF4A. The published variants in BLK, KLF11, and PAX4 showed poor cosegregation with diabetes (combined logarithm of the odds [LOD] scores ≤1.2), compared with HNF1A and HNF4A (LOD scores >9), and are all too common to cause MODY (minor allele frequency >4.95 × 10-5). Ultra-rare missense and protein-truncating variants (PTV) were not enriched in a MODY cohort compared with the UK Biobank population (PTV P > 0.05, missense P > 0.1 for all three genes) while HNF1A and HNF4A were enriched (P < 10-6). Findings of sensitivity analyses with different population cohorts supported our results. Variant and gene-level genetic evidence does not support BLK, KLF11, or PAX4 as a cause of MODY. They should not be included in MODY diagnostic genetic testing.

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.

CONTEXT: PLIN1 encodes perilipin-1, which coats lipid droplets in adipocytes and is involved in droplet formation, triglyceride storage, and lipolysis. Rare PLIN1 frameshift variants that extend the translated protein have been described to cause lipodystrophy. OBJECTIVE: This work aimed to test whether PLIN1 protein-truncating variants (PTVs) cause lipodystrophy in a large population-based cohort. METHODS: We identified individuals with PLIN1 PTVs in individuals with exome data in the UK Biobank. We performed gene-burden testing for individuals with PLIN1 PTVs. We replicated the associations using data from the T2D Knowledge portal. We performed a phenome-wide association study using publicly available association statistics. A total of 362 791 individuals in the UK Biobank, a population-based cohort, and 43 125 individuals in the T2D Knowledge portal, a type 2 diabetes (T2D) case-control study, were included in the analyses. Main outcome measures included 22 diseases and traits relevant to lipodystrophy. RESULTS: the 735 individuals with PLIN1 PTVs had a favorable metabolic profile. These individuals had increased high-density lipoprotein cholesterol (0.12 mmol/L; 95% CI, 0.09 to 0.14, P�
=�
2�
�
10-18), reduced triglycerides (-0.22 mmol/L; 95% CI, -0.29 to -0.14, P�
=�
3�
�
10-11), reduced waist-to-hip ratio (-0.02; 95% CI, -0.02 to -0.01, P�
=�
9�
�
10-12), and reduced systolic blood pressure (-1.67 mm Hg; 95% CI, -3.25 to -0.09, P�
=�
.05). These associations were consistent in the smaller T2D Knowledge portal cohort. In the UK Biobank, PLIN1 PTVs were associated with reduced risk of myocardial infarction (odds ratio [OR]�
=�
0.59; 95% CI, 0.35 to 0.93, P�
=�
.02) and hypertension (OR�
=�
0.85; 95% CI, 0.73 to 0.98, P�
=�
.03), but not T2D (OR�
=�
0.99; 95% CI, 0.63-1.51, P�
=�
.99). CONCLUSION: Our study suggests that PLIN1 haploinsufficiency causes a favorable metabolic profile and may protect against cardiovascular disease.

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.

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.

Gene expression is tightly regulated with many genes exhibiting cell-specific silencing when their protein product would disrupt normal cellular function. This silencing is largely controlled by non-coding elements and their disruption might cause human disease. 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 mutations affecting a 42bp conserved region encompassed by a regulatory element in intron 2 of Hexokinase 1 ( HK1 ), a pancreatic beta-cell disallowed gene. We demonstrated that these mutations resulted in expression of HK1 in the pancreatic beta-cells causing inappropriate insulin secretion and congenital hyperinsulinism. These mutations identify a regulatory region critical for cell-specific silencing. Importantly, this has revealed a new disease mechanism for non-coding mutations that cause inappropriate expression of a disallowed gene.

Objective
Mutations in the KATP channel genes, ABCC8 and KCNJ11, are the most common cause of congenital hyperinsulinism. The diagnosis of KATP-hyperinsulinism is important for the clinical management of the condition. We aimed to determine the clinical features that help to identify KATP-hyperinsulinism at diagnosis.


Design
We studied 761 individuals with KATP-hyperinsulinism and 862 probands with hyperinsulinism of unknown aetiology diagnosed before 6 months of age. All were referred as part of routine clinical care.


Methods
We compared the clinical features of KATP-hyperinsulinism and unknown hyperinsulinism cases. We performed logistic regression and receiver operator characteristic (ROC) analysis to identify the features that predict KATP-hyperinsulinism.


Results
Higher birth weight, diazoxide unresponsiveness and diagnosis in the first week of life were independently associated with KATP-hyperinsulinism (adjusted odds ratio: 4.5 (95% CI: 3.4–5.9), 0.09 (0.06–0.13) and 3.3 (2.0–5.0) respectively). Birth weight and diazoxide unresponsiveness were additive and highly discriminatory for identifying KATP-hyperinsulinism (ROC area under the curve for birth weight 0.80, diazoxide responsiveness 0.77, and together 0.88, 95% CI: 0.85–0.90). In this study, 86% born large for gestation and 78% born appropriate for gestation and who did not respond to diazoxide treatment had KATP-hyperinsulinism. In contrast, of those individuals born small for gestation, none who were diazoxide responsive and only 4% of those who were diazoxide unresponsive had KATP-hyperinsulinism.


Conclusions
Individuals with hyperinsulinism born appropriate or large for gestation and unresponsive to diazoxide treatment are most likely to have an ABCC8 or KCNJ11 mutation. These patients should be prioritised for genetic testing of KATP channel genes.

Multiple Nucleotide Variants (MNVs) are miscalled by the most widely utilised next generation sequencing analysis (NGS) pipelines, presenting the potential for missing diagnoses. These variants, which should be treated as a single insertion-deletion mutation event, are commonly called as separate single nucleotide variants. This can result in misannotation, incorrect amino acid predictions and potentially false positive and false negative diagnostic results. Using simulated data and re-analysis of sequencing data from a diagnostic targeted gene panel, we demonstrate that the widely adopted pipeline, GATK best practices, results in miscalling of MNVs and that alternative tools can call these variants correctly. The adoption of calling methods that annotate MNVs correctly would present a solution for individual laboratories, however GATK best practices are the basis for important public resources such as the gnomAD database. We suggest integrating a solution into these guidelines would be the optimal approach.

Multiple Nucleotide Variants (MNVs) are miscalled by the most widely utilised next generation sequencing analysis (NGS) pipelines, presenting the potential for missing diagnoses. These variants, which should be treated as a single insertion-deletion mutation event, are commonly called as separate single nucleotide variants. This can result in misannotation, incorrect amino acid predictions and potentially false positive and false negative diagnostic results. Using simulated data and re-analysis of sequencing data from a diagnostic targeted gene panel, we demonstrate that the widely adopted pipeline, GATK best practices, results in miscalling of MNVs and that alternative tools can call these variants correctly. The adoption of calling methods that annotate MNVs correctly would present a solution for individual laboratories, however GATK best practices are the basis for important public resources such as the gnomAD database. We suggest integrating a solution into these guidelines would be the optimal approach.

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.

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

PURPOSE: One of the greatest challenges currently facing those studying Mendelian disease is identifying the pathogenic variant from the long list produced by a next-generation sequencing test. We investigate the predictive ability of homozygosity mapping for identifying the regions likely to contain the causative variant. METHODS: We use 179 homozygous pathogenic variants from three independent cohorts to investigate the predictive power of homozygosity mapping. RESULTS: We demonstrate that homozygous pathogenic variants in our cohorts are disproportionately likely to be found within one of the largest regions of homozygosity: 80% of pathogenic variants are found in a homozygous region that is in the ten largest regions in a sample. The maximal predictive power is achieved in patients with 3 Mb from a telomere; this gives an area under the curve (AUC) of 0.735 and results in 92% of the causative variants being in one of the ten largest homozygous regions. CONCLUSION: This predictive power can be used to prioritize the list of candidate variants in gene discovery studies. When classifying a homozygous variant the size and rank of the region of homozygosity in which the candidate variant is located can also be considered as supporting evidence for pathogenicity.

Multiple Nucleotide Variants (MNVs) are miscalled by the most widely utilised next generation sequencing analysis (NGS) pipelines, presenting the potential for missing diagnoses that would previously have been made by standard Sanger (dideoxy) sequencing. These variants, which should be treated as a single insertion-deletion mutation event, are commonly called as separate single nucleotide variants. This can result in misannotation, incorrect amino acid predictions and potentially false positive and false negative diagnostic results. This risk will be increased as confirmatory Sanger sequencing of Single Nucleotide variants (SNVs) ceases to be standard practice. Using simulated data and re-analysis of sequencing data from a diagnostic targeted gene panel, we demonstrate that the widely adopted pipeline, GATK best practices, results in miscalling of MNVs and that alternative tools can call these variants correctly. The adoption of calling methods that annotate MNVs correctly would present a solution for individual laboratories, however GATK best practices are the basis for important public resources such as the gnomAD database. We suggest integrating a solution into these guidelines would be the optimal approach.

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.

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.

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.

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.

Congenital hyperinsulinaemic hypoglycaemia (HH) can occur in isolation or it may present as part of a wider syndrome. For approximately 40%-50% of individuals with this condition, sequence analysis of the known HH genes identifies a causative mutation. Identifying the underlying genetic aetiology in the remaining cases is important as a genetic diagnosis will inform on recurrence risk, may guide medical management and will provide valuable insights into β-cell physiology. We sequenced the exome of a child with persistent diazoxide-responsive HH, mild aortic insufficiency, severe hypotonia, and developmental delay as well as the unaffected parents. This analysis identified a de novo mutation, p.G403D, in the proband's CACNA1D gene. CACNA1D encodes the main L-type voltage-gated calcium channel in the pancreatic β-cell, a key component of the insulin secretion pathway. The p.G403D mutation had been reported previously as an activating mutation in an individual with primary hyper-aldosteronism, neuromuscular abnormalities, and transient hypoglycaemia. Sequence analysis of the CACNA1D gene in 60 further cases with HH did not identify a pathogenic mutation. Identification of an activating CACNA1D mutation in a second patient with congenital HH confirms the aetiological role of CACNA1D mutations in this disorder. A genetic diagnosis is important as treatment with a calcium channel blocker may be an option for the medical management of this patient.

Finding new causes of monogenic diabetes helps understand glycaemic regulation in humans. To find novel genetic causes of maturity-onset diabetes of the young (MODY), we sequenced MODY cases with unknown aetiology and compared variant frequencies to large public databases. From 36 European patients, we identify two probands with novel RFX6 heterozygous nonsense variants. RFX6 protein truncating variants are enriched in the MODY discovery cohort compared to the European control population within ExAC (odds ratio = 131, P = 1 × 10-4). We find similar results in non-Finnish European (n = 348, odds ratio = 43, P = 5 × 10-5) and Finnish (n = 80, odds ratio = 22, P = 1 × 10-6) replication cohorts. RFX6 heterozygotes have reduced penetrance of diabetes compared to common HNF1A and HNF4A-MODY mutations (27, 70 and 55% at 25 years of age, respectively). The hyperglycaemia results from beta-cell dysfunction and is associated with lower fasting and stimulated gastric inhibitory polypeptide (GIP) levels. Our study demonstrates that heterozygous RFX6 protein truncating variants are associated with MODY with reduced penetrance.Maturity-onset diabetes of the young (MODY) is the most common subtype of familial diabetes. Here, Patel et al. use targeted DNA sequencing of MODY patients and large-scale publically available data to show that RFX6 heterozygous protein truncating variants cause reduced penetrance MODY.

The long-read sequencers from Pacific Bioscience (PacBio) and Oxford Nanopore Technologies (ONT) offer the opportunity to phase mutations multiple kilobases apart directly from sequencing reads. In this study, we used long-range PCR with ONT and PacBio sequencing to phase two variants 9 kb apart in the RET gene. We also re-analysed data from a recent paper which had apparently successfully used ONT to phase clinically important haplotypes at the CYP2D6 and HLA loci. From these analyses, we demonstrate PCR-chimera formation during PCR amplification and reference alignment bias are pitfalls that need to be considered when attempting to phase variants using amplicon-based long-read sequencing technologies. These methodological pitfalls need to be avoided if the opportunities provided by long-read sequencers are to be fully exploited.

HNF4A mutations cause increased birth weight, transient neonatal hypoglycemia, and maturity onset diabetes of the young (MODY). The most frequently reported HNF4A mutation is p.R114W (previously p.R127W), but functional studies have shown inconsistent results; there is a lack of cosegregation in some pedigrees and an unexpectedly high frequency in public variant databases. We confirm that p.R114W is a pathogenic mutation with an odds ratio of 30.4 (95% CI 9.79-125, P = 2 × 10(-21)) for diabetes in our MODY cohort compared with control subjects. p.R114W heterozygotes did not have the increased birth weight of patients with other HNF4A mutations (3,476 g vs. 4,147 g, P = 0.0004), and fewer patients responded to sulfonylurea treatment (48% vs. 73%, P = 0.038). p.R114W has reduced penetrance; only 54% of heterozygotes developed diabetes by age 30 years compared with 71% for other HNF4A mutations. We redefine p.R114W as a pathogenic mutation that causes a distinct clinical subtype of HNF4A MODY with reduced penetrance, reduced sensitivity to sulfonylurea treatment, and no effect on birth weight. This has implications for diabetes treatment, management of pregnancy, and predictive testing of at-risk relatives. The increasing availability of large-scale sequence data is likely to reveal similar examples of rare, low-penetrance MODY mutations.

The Oxford Nanopore Technologies (ONT) MinION is a new sequencing technology that potentially offers read lengths of tens of kilobases (kb) limited only by the length of DNA molecules presented to it. The device has a low capital cost, is by far the most portable DNA sequencer available, and can produce data in real-time. It has numerous prospective applications including improving genome sequence assemblies and resolution of repeat-rich regions. Before such a technology is widely adopted, it is important to assess its performance and limitations in respect of throughput and accuracy. In this study we assessed the performance of the MinION by re-sequencing three bacterial genomes, with very different nucleotide compositions ranging from 28.6% to 70.7%; the high G. +. C strain was underrepresented in the sequencing reads. We estimate the error rate of the MinION (after base calling) to be 38.2%. Mean and median read lengths were 2. kb and 1. kb respectively, while the longest single read was 98. kb. The whole length of a 5. kb rRNA operon was covered by a single read. As the first nanopore-based single molecule sequencer available to researchers, the MinION is an exciting prospect; however, the current error rate limits its ability to compete with existing sequencing technologies, though we do show that MinION sequence reads can enhance contiguity of de novo assembly when used in conjunction with Illumina MiSeq data.

The application of high-throughput sequencing in profiling microbial communities is providing an unprecedented ability to investigate microbiomes. Such studies typically apply one of two methods: amplicon sequencing using PCR to target a conserved orthologous sequence (typically the 16S ribosomal RNA gene) or whole (meta)genome sequencing (WGS). Both methods have been used to catalog the microbial taxa present in a sample and quantify their respective abundances. However, a comparison of the inherent precision or bias of the different sequencing approaches has not been performed. We previously developed a metagenomic control material (MCM) to investigate error when performing different sequencing strategies. Amplicon sequencing using four different primer strategies and two 16S rRNA regions was examined (Roche 454 Junior) and compared to WGS (Illumina HiSeq). All sequencing methods generally performed comparably and in good agreement with organism specific digital PCR (dPCR); WGS notably demonstrated very high precision. Where discrepancies between relative abundances occurred they tended to differ by less than twofold. Our findings suggest that when alternative sequencing approaches are used for microbial molecular profiling they can perform with good reproducibility, but care should be taken when comparing small differences between distinct methods. This work provides a foundation for future work comparing relative differences between samples and the impact of extraction methods. We also highlight the value of control materials when conducting microbial profiling studies to benchmark methods and set appropriate thresholds.

Advances in DNA sequencing technology provide the possibility to analyse and characterize the genetic material from microbial populations (the microbiome) as a whole. Such comprehensive analysis of a microbiome using these 'metagenomic' approaches offers the potential to understand industrial, clinical and environmental microbiology to a level of detail that is unfeasible using conventional molecular or culture-based methods. However, the complexity offered by metagenomic analysis is also the weakness of this method and poses considerable challenges during analytical standardisation. In this manuscript, we discuss options for developing control materials for metagenomic analysis and describe our preliminary work investigating how such materials can be used to assist metagenomic measurements. The control materials we have developed demonstrate that, when performing 16S rDNA sequencing, different library preparation methods (incorporating adapters before and after the PCR) and small primer mismatches can alter the reported metagenomic profile. These findings illustrate that metagenomic analysis can be heavily biased by the choice of method and underpin the need for control materials that can provide a useful tool in informing choice of protocol for accurate analysis. © 2012 the Author(s).