Unique Capabilities of Genome Sequencing for Rare Disease Diagnosis

Background

Causal variants underlying rare disorders may remain elusive even after expansive gene panels or exome sequencing (ES). Clinicians and researchers may then turn to genome sequencing (GS), though the added value of this technique and its optimal use remain poorly defined. We therefore investigated the advantages of GS within a phenotypically diverse cohort.

Methods

GS was performed for 744 individuals with rare disease who were genetically undiagnosed. Analysis included review of single nucleotide, indel, structural, and mitochondrial variants.

Results

We successfully solved 218/744 (29.3%) cases using GS, with most solves involving established disease genes (157/218, 72.0%). Of all solved cases, 148 (67.9%) had previously had non-diagnostic ES. We systematically evaluated the 218 causal variants for features requiring GS to identify and 61/218 (28.0%) met these criteria, representing 8.2% of the entire cohort. These included small structural variants (13), copy neutral inversions and complex rearrangements (8), tandem repeat expansions (6), deep intronic variants (15), and coding variants that may be more easily found using GS related to uniformity of coverage (19).

Conclusion

We describe the diagnostic yield of GS in a large and diverse cohort, illustrating several types of pathogenic variation eluding ES or other techniques. Our results reveal a higher diagnostic yield of GS, supporting the utility of a genome-first approach, with consideration of GS as a secondary or tertiary test when higher-resolution structural variant analysis is needed or there is a strong clinical suspicion for a condition and prior targeted genetic testing has been negative.

Biallelic and monoallelic variants in PLXNA1 are implicated in a novel neurodevelopmental disorder with variable cerebral and eye anomalies

PURPOSE

To investigate the effect of PLXNA1 variants on the phenotype of patients with autosomal dominant and recessive inheritance patterns and to functionally characterize the zebrafish homologs plxna1a and plxna1b during development.

METHODS

We assembled ten patients from seven families with biallelic or de novo PLXNA1 variants. We describe genotype-phenotype correlations, investigated the variants by structural modeling, and used Morpholino knockdown experiments in zebrafish to characterize the embryonic role of plxna1a and plxna1b.

RESULTS

Shared phenotypic features among patients include global developmental delay (9/10), brain anomalies (6/10), and eye anomalies (7/10). Notably, seizures were predominantly reported in patients with monoallelic variants. Structural modeling of missense variants in PLXNA1 suggests distortion in the native protein. Our zebrafish studies enforce an embryonic role of plxna1a and plxna1b in the development of the central nervous system and the eye.

CONCLUSION

We propose that different biallelic and monoallelic variants in PLXNA1 result in a novel neurodevelopmental syndrome mainly comprising developmental delay, brain, and eye anomalies. We hypothesize that biallelic variants in the extracellular Plexin-A1 domains lead to impaired dimerization or lack of receptor molecules, whereas monoallelic variants in the intracellular Plexin-A1 domains might impair downstream signaling through a dominant-negative effect.

Bi-allelic Loss-of-Function CACNA1B Mutations in Progressive Epilepsy-Dyskinesia

The occurrence of non-epileptic hyperkinetic movements in the context of developmental epileptic encephalopathies is an increasingly recognized phenomenon. Identification of causative mutations provides an important insight into common pathogenic mechanisms that cause both seizures and abnormal motor control. We report bi-allelic loss-of-function CACNA1B variants in six children from three unrelated families whose affected members present with a complex and progressive neurological syndrome. All affected individuals presented with epileptic encephalopathy, severe neurodevelopmental delay (often with regression), and a hyperkinetic movement disorder. Additional neurological features included postnatal microcephaly and hypotonia. Five children died in childhood or adolescence (mean age of death: 9 years), mainly as a result of secondary respiratory complications. CACNA1B encodes the pore-forming subunit of the pre-synaptic neuronal voltage-gated calcium channel Ca_v2.2/N-type, crucial for SNARE-mediated neurotransmission, particularly in the early postnatal period. Bi-allelic loss-of-function variants in CACNA1B are predicted to cause disruption of Ca²⁺ influx, leading to impaired synaptic neurotransmission. The resultant effect on neuronal function is likely to be important in the development of involuntary movements and epilepsy. Overall, our findings provide further evidence for the key role of Ca_v2.2 in normal human neurodevelopment.

Resolving the full spectrum of human genome variation using Linked-Reads

Large-scale population analyses coupled with advances in technology have demonstrated that the human genome is more diverse than originally thought. To date, this diversity has largely been uncovered using short-read whole-genome sequencing. However, these short-read approaches fail to give a complete picture of a genome. They struggle to identify structural events, cannot access repetitive regions, and fail to resolve the human genome into haplotypes. Here, we describe an approach that retains long range information while maintaining the advantages of short reads. Starting from ∼1 ng of high molecular weight DNA, we produce barcoded short-read libraries. Novel informatic approaches allow for the barcoded short reads to be associated with their original long molecules producing a novel data type known as "Linked-Reads". This approach allows for simultaneous detection of small and large variants from a single library. In this manuscript, we show the advantages of Linked-Reads over standard short-read approaches for reference-based analysis. Linked-Reads allow mapping to 38 Mb of sequence not accessible to short reads, adding sequence in 423 difficult-to-sequence genes including disease-relevant genes STRC, SMN1, and SMN2 Both Linked-Read whole-genome and whole-exome sequencing identify complex structural variations, including balanced events and single exon deletions and duplications. Further, Linked-Reads extend the region of high-confidence calls by 68.9 Mb. The data presented here show that Linked-Reads provide a scalable approach for comprehensive genome analysis that is not possible using short reads alone.

Contribution of noncoding pathogenic variants to RPGRIP1-mediated inherited retinal degeneration

Purpose:

With the advent of gene therapies for inherited retinal degenerations (IRDs), genetic diagnostics will have an increasing role in clinical decision-making. Yet the genetic cause of disease cannot be identified using exon-based sequencing for a significant portion of patients. We hypothesized that non-coding mutations contribute significantly to the genetic causality of IRDs and evaluated patients with single coding mutations in RPGRIP1 to test this hypothesis.

Methods:

IRD families underwent targeted panel sequencing. Unsolved cases were explored by whole exome and genome sequencing looking for additional mutations. Candidate mutations were then validated by Sanger sequencing, quantitative PCR, and in vitro splicing assays in two cell lines analyzed through amplicon sequencing.

Results:

Among 1722 families, three had biallelic loss of function mutations in RPGRIP1 while seven had a single disruptive coding mutation. Whole exome and genome sequencing revealed potential non-coding mutations in these seven families. In six, the non-coding mutations were shown to lead to loss of function in vitro.

Conclusion:

Non-coding mutations were identified in 6 of 7 families with single coding mutations in RPGRIP1. The results suggest that non-coding mutations contribute significantly to the genetic causality of IRDs and RPGRIP1–mediated IRDs are more common than previously thought.

A synthetic-diploid benchmark for accurate variant-calling evaluation

Existing benchmark datasets for use in evaluating variant-calling accuracy are constructed from a consensus of known short-variant callers, and they are thus biased toward easy regions that are accessible by these algorithms. We derived a new benchmark dataset from the de novo PacBio assemblies of two fully homozygous human cell lines, which provides a relatively more accurate and less biased estimate of small-variant-calling error rates in a realistic context.

Quantifying the Impact of Rare and Ultra-rare Coding Variation across the Phenotypic Spectrum

There is a limited understanding about the impact of rare protein-truncating variants across multiple phenotypes. We explore the impact of this class of variants on 13 quantitative traits and 10 diseases using whole-exome sequencing data from 100,296 individuals. Protein-truncating variants in genes intolerant to this class of mutations increased risk of autism, schizophrenia, bipolar disorder, intellectual disability, and ADHD. In individuals without these disorders, there was an association with shorter height, lower education, increased hospitalization, and reduced age at enrollment. Gene sets implicated from GWASs did not show a significant protein-truncating variants burden beyond what was captured by established Mendelian genes. In conclusion, we provide a thorough investigation of the impact of rare deleterious coding variants on complex traits, suggesting widespread pleiotropic risk.

A Loss-of-Function Splice Acceptor Variant in IGF2 Is Protective for Type 2 Diabetes

Type 2 diabetes (T2D) affects more than 415 million people worldwide, and its costs to the health care system continue to rise. To identify common or rare genetic variation with potential therapeutic implications for T2D, we analyzed and replicated genome-wide protein coding variation in a total of 8,227 individuals with T2D and 12,966 individuals without T2D of Latino descent. We identified a novel genetic variant in the IGF2 gene associated with ∼20% reduced risk for T2D. This variant, which has an allele frequency of 17% in the Mexican population but is rare in Europe, prevents splicing between IGF2 exons 1 and 2. We show in vitro and in human liver and adipose tissue that the variant is associated with a specific, allele-dosage-dependent reduction in the expression of IGF2 isoform 2. In individuals who do not carry the protective allele, expression of IGF2 isoform 2 in adipose is positively correlated with both incidence of T2D and increased plasma glycated hemoglobin in individuals without T2D, providing support that the protective effects are mediated by reductions in IGF2 isoform 2. Broad phenotypic examination of carriers of the protective variant revealed no association with other disease states or impaired reproductive health. These findings suggest that reducing IGF2 isoform 2 expression in relevant tissues has potential as a new therapeutic strategy for T2D, even beyond the Latin American population, with no major adverse effects on health or reproduction.

Exome sequences versus sequential gene testing in the UK highly specialised Service for Limb Girdle Muscular Dystrophy

BACKGROUND

Limb girdle muscular dystrophies are a group of rare and genetically heterogeneous diseases that share proximal weakness as a common feature; however they are often lacking very specific phenotypic features to allow an accurate differential diagnosis based on the clinical signs only, limiting the diagnostic rate using phenotype driven genetic testing. Next generation sequencing provides an opportunity to obtain molecular diagnoses for undiagnosed patients, as well as identifying novel genetic causes of muscle diseases. We performed whole exome sequencing (WES) on 104 affected individuals from 75 families in who standard gene by gene testing had not yielded a diagnosis. For comparison we also evaluated the diagnostic rate using sequential gene by gene testing for 91 affected individuals from 84 families over a 2 year period.

RESULTS

Patients selected for WES had undergone more extensive prior testing than those undergoing standard genetic testing and on average had had 8 genes screened already. In this extensively investigated cohort WES identified the genetic diagnosis in 28 families (28/75, 37%), including the identification of the novel gene ZAK and two unpublished genes. WES of a single affected individual with sporadic disease yielded a diagnosis in 13/38 (34%) of cases. In comparison, conventional gene by gene testing provided a genetic diagnosis in 28/84 (33%) families. Titinopathies and collagen VI related dystrophy were the most frequent diagnoses made by WES. Reasons why mutations in known genes were not identified previously included atypical phenotypes, reassignment of pathogenicity of variants, and in one individual mosaicism for a COL6A1 mutation which was undetected by prior direct sequencing.

CONCLUSION

WES was able to overcome many limitations of standard testing and achieved a higher rate of diagnosis than standard testing even in this cohort of extensively investigated patients. Earlier application of WES is therefore likely to yield an even higher diagnostic rate. We obtained a high diagnosis rate in simplex cases and therefore such individuals should be included in exome or genome sequencing projects. Disease due to somatic mosaicism may be increasingly recognised due to the increased sensitivity of next generation sequencing techniques to detect low level mosaicism.

A 'second truncation' in TTN causes early onset recessive muscular dystrophy

Mutations in the gene encoding the giant skeletal muscle protein titin are associated with a variety of muscle disorders, including recessive congenital myopathies ±cardiomyopathy, limb girdle muscular dystrophy (LGMD) and late onset dominant distal myopathy. Heterozygous truncating mutations have also been linked to dilated cardiomyopathy. The phenotypic spectrum of titinopathies is emerging and expanding, as next generation sequencing techniques make this large gene amenable to sequencing. We undertook whole exome sequencing in four individuals with LGMD. An essential splice site mutation, previously reported in dilated cardiomyopathy, was identified in all families in combination with a second truncating mutation. Affected individuals presented with childhood onset proximal weakness associated with joint contractures and elevated CK. Cardiac dysfunction was present in two individuals. Muscle biopsy showed increased internal nuclei and immunoblotting identified reduction or absence of calpain-3 and demonstrated a marked reduction of C-terminal titin fragments. We confirm the co-occurrence of cardiac and skeletal myopathies associated with recessive truncating titin mutations. Compound heterozygosity of a truncating mutation previously associated with dilated cardiomyopathy and a 'second truncation' in TTN was identified as causative in our skeletal myopathy patients. These findings add to the complexity of interpretation and genetic counselling for titin mutations.

Using high-resolution variant frequencies to empower clinical genome interpretation

Purpose

Whole-exome and whole-genome sequencing have transformed the discovery of genetic variants that cause human Mendelian disease, but discriminating pathogenic from benign variants remains a daunting challenge. Rarity is recognized as a necessary, although not sufficient, criterion for pathogenicity, but frequency cutoffs used in Mendelian analysis are often arbitrary and overly lenient. Recent very large reference datasets, such as the Exome Aggregation Consortium (ExAC), provide an unprecedented opportunity to obtain robust frequency estimates even for very rare variants.

Methods

We present a statistical framework for the frequency-based filtering of candidate disease-causing variants, accounting for disease prevalence, genetic and allelic heterogeneity, inheritance mode, penetrance, and sampling variance in reference datasets.

Results

Using the example of cardiomyopathy, we show that our approach reduces by two-thirds the number of candidate variants under consideration in the average exome, without removing true pathogenic variants (false-positive rate<0.001).

Conclusion

We outline a statistically robust framework for assessing whether a variant is “too common” to be causative for a Mendelian disorder of interest. We present precomputed allele frequency cutoffs for all variants in the ExAC dataset.

Evaluation of HFC 245ca and HFC 236ea as Foam Blowing Agents

Hydrochlorofluorocarbons (HCFCs) have been widely adopted as an interim blowing agent in urethane insulations. Due to their scheduled phase-out, research efforts are being devoted to the identification and development of alternatives with zero ozone depletion potential. Physical blowing agents identified have included hydrocarbons, fluorocarbons, hydrofluoroethers, and more predominantly, hydrofluorocarbons (HFCs). The majority of the HFC evaluations have focused on the more readily available, low boiling candidates such as HFC 134a. Higher boiling HFC candidates that could be handled at ambient conditions and use current processing equipment would be more desirable.

This paper will describe results from a research program of two such candidate HFCs performed as a cooperative effort between the Environmental Protection Agency, Oak Ridge National Laboratory and Martin Marietta Manned Space Systems. The purpose of this effort was to perform a limited evaluation of the developmental HFCs 245ca and 236ea as blowing agents in urethane based insulations. These two materials were selected from screening tests of 37 C2, C3 and C4 isomers based on physical properties, atmospheric lifetime, flammability, estimated toxicity, difficulty of synthesis, suitability for dual use as a refrigerant and other factors. Solubility of the two materials in typical foam components was tested, pour foaming trials were performed and preliminary data were gathered regarding foam insulation performance.

Temple-Baraitser Syndrome and Zimmermann-Laband Syndrome: one clinical entity?

BACKGROUND

KCNH1 encodes a voltage-gated potassium channel that is predominantly expressed in the central nervous system. Mutations in this gene were recently found to be responsible for Temple-Baraitser Syndrome (TMBTS) and Zimmermann-Laband syndrome (ZLS).

METHODS

Here, we report a new case of TMBTS diagnosed in a Lebanese child. Whole genome sequencing was carried out on DNA samples of the proband and his parents to identify mutations associated with this disease. Sanger sequencing was performed to confirm the presence of detected variants.

RESULTS

Whole genome sequencing revealed three missense mutations in TMBTS patient: c.1042G > A in KCNH1, c.2131 T > C in STK36, and c.726C > A in ZNF517. According to all predictors, mutation in KCNH1 is damaging de novo mutation that results in substitution of Glycine by Arginine, i.e., p.(Gly348Arg). This mutation was already reported in a patient with ZLS that could affect the connecting loop between helices S4-S5 of KCNH1 with a gain of function effect.

CONCLUSIONS

Our findings demonstrate that KCNH1 mutations cause TMBTS and expand the mutational spectrum of KCNH1 in TMBTS. In addition, all cases of TMBTS were reviewed and compared to ZLS. We suggest that the two syndromes are a continuum and that the variability in the phenotypes is the result of the involvement of genetic modifiers.

RNAseq analysis for the diagnosis of muscular dystrophy

The precise genetic cause remains elusive in nearly 50% of patients with presumed neurogenetic disease, representing a significant barrier for clinical care. This is despite significant advances in clinical genetic diagnostics, including the application of whole‐exome sequencing and next‐generation sequencing‐based gene panels. In this study, we identify a deep intronic mutation in the DMD gene in a patient with muscular dystrophy using both conventional and RNAseq‐based transcriptome analyses. The implications of our data are that noncoding mutations likely comprise an important source of unresolved genetic disease and that RNAseq is a powerful platform for detecting such mutations.

Advancing the promise of genomic medicine

Genomic technologies present powerful tools with the potential to diagnose, treat, and elucidate the underpinnings of human disease. Technological advances have allowed researchers and clinicians to apply novel approaches and gain statistical power to delve more deeply and broadly into our genomes than ever before, enabling them to discover associations between genomic variants and disease phenotypes that were previously beyond our reach. In this webinar, we will explore how genomics is currently being applied to characterize rare disease and cancer, highlight recent advances in these respective areas, and look toward the future to predict what is needed to positively impact human health.

The Challenge of Next Generation Sequencing in the Context of Neuromuscular Diseases

New genomic technologies, such as exome, whole-genome and transcriptome sequencing, are transforming the genetic diagnosis of neuromuscular diseases and dramatically accelerating the discovery of new disease-associated genes. The increasingly widespread availability of these technologies creates both opportunities and challenges for neuromuscular disease researchers. Here we survey the current literature on the application of new genomic technologies to the diagnosis of severe muscle diseases, with a focus on assessment of the approaches used for data processing, analysis and interpretation. We also highlight several key areas requiring improvement.

HI: haplotype improver using paired-end short reads

Summary: We present a program to improve haplotype reconstruction by incorporating information from paired-end reads, and demonstrate its utility on simulated data. We find that given a fixed coverage, longer reads (implying fewer of them) are preferable.

Availability: The executable and user manual can be freely downloaded from ftp://ftp.sanger.ac.uk/pub/zn1/HI.

Contact:ql2@sanger.ac.uk

Intrathecal chemotherapy delivered by a lumbar-thecal catheter in metastatic medulloblastoma: a case illustration

BACKGROUND

Medulloblastoma is the most common malignant brain tumour in children. Despite recent advances, the prognosis in high risk patients remains poor. Further improvement in survival is dependent upon the development of strategies to attack the tumour more effectively, but with less toxicity. Intrathecal chemotherapy, is an ideal but currently underused method of directly targeting residual tumour within the area of resection and the leptomeningeal disease commonly associated with this tumour.

METHOD

We describe the case of a 12 yr old child with metastatic medulloblastoma, who received intrathecal topotecan via a spinal catheter.

CONCLUSION

This method represents a simple, safe and effective method of delivering an even and widespread distribution of drug within the cerebrospinal fluid (CSF) of the neuroaxis. With new agents being identified and others in the early stages of development, intrathecal chemotherapy may emerge as an important therapeutic option to consider when faced with such challenging cases.

Oxygen costs using a reciprocating gait orthosis in a paraplegic (T9) patient with a bilateral below-knee amputation: case report

The Reciprocating Gait Orthosis (RGO) is a useful aid to ambulation for patients with paraplegia. Its use has been described previously though not in conjunction with limb prostheses. We report here the energy costs of ambulation of a patient, disabled by paraplegia at T9 and bilateral below-knee amputations, walking at her preferred rate using an RGO while gas exchange was measured by the Douglas bag method. Oxygen uptake (VO2) rose from 0.198 1 min-1 at rest to 0.582 1 min-1 in the last minute of exercise, representing a VO2 of 14.3 ml kg-1 min-1. During the fourth minute of ambulation, energy consumption was 30.44 J kg-1 s-1 with an energy cost of 4.17 J kg-1 m-1 at a velocity of 0.13 m s-1, Ambulation with this combination of disability is possible with the aid of limb prostheses and an RGO though it is slow and the energy expenditure as consumption per second and cost per metre are high.

An artificial neural network system for diagnosis of acute myocardial infarction (AMI) in the accident and emergency department: evaluation and comparison with serum myoglobin measurements

Recent studies have confirmed that artificial neural networks (ANNs) are adept at recognising patterns in sets of clinical data. The diagnosis of acute myocardial infarction (AMI) in patients presenting with chest pain remains one of the greatest challenges in emergency medicine. The aim of this study was to evaluate the performance of an ANN trained to analyse clinical data from chest pain patients. The ANN was compared with serum myoglobin measurements--cardiac damage is associated with increased circulating myoglobin levels, and this is widely used as an early marker for evolving AMI. We used 39 items of clinical and ECG data from the time of presentation to derive 53 binary inputs to a back propagation network. On test data (200 cases), overall accuracy, sensitivity, specificity and positive predictive value (PPV) of the ANN were 91.8, 91.2, 90.2 and 84.9% respectively. Corresponding figures using linear discriminant analysis were 81.0, 77.9, 82.6 and 69.7% (P < 0.01). Using a further test set from a different centre (91 cases), the accuracy, sensitivity, specificity and PPV for the admitting physicians were 65.1, 28.5, 76.9 and 28.6% respectively compared with 73.6, 52.4, 80.0 and 44.0% for the ANN. Although myoglobin at presentation was highly specific, it was only 38.0% sensitive, compared with 85.7% at 3 h. Simple strategies to combine clinical opinion, ANN output and myoglobin at presentation could greatly improve sensitivity and specificity of AMI diagnosis. The ideal support for emergency room physicians may come from a combination of computer-aided analysis of clinical factors and biochemical markers such as myoglobin. This study demonstrates that the two approaches could be usefully combined, the major benefit of the decision support system being in the first 3 h before biochemical markers have become abnormal.

A Comparison of Alternative Techniques of Instruction

The present study examined some of the problems involved in communication and learning processes when inexperienced subjects were presented with information supplied through unfamiliar types of instruction media. In the study, 110 Xhosa-speaking novices to the gold mining industry were instructed through the media of video-tape or lecture, in either Fanakalo or Xhosa, with repeated or non-repeated presentations. Results showed that those who received repeated instruction recalled more information. A significant interaction also occurred between language and instruction medium. Level of education also affected comprehension. It appeared that the more the presentation methods were tailored to the specific needs of the individuals — in this case presentation through a familiar language, using an empathie lecturer — the more effective the communication and learning was likely to be.