Transcriptome variation in human tissues revealed by long-read sequencing.

Regulation of transcript structure generates transcript diversity and plays an important role in human disease1-7. The advent of long-read sequencing technologies offers the opportunity to study the role of genetic variation in transcript structure8-16. In this Article, we present a large human long-read RNA-seq dataset using the Oxford Nanopore Technologies platform from 88 samples from Genotype-Tissue Expression (GTEx) tissues and cell lines, complementing the GTEx resource. We identified just over 70,000 novel transcripts for annotated genes, and validated the protein expression of 10% of novel transcripts. We developed a new computational package, LORALS, to analyse the genetic effects of rare and common variants on the transcriptome by allele-specific analysis of long reads. We characterized allele-specific expression and transcript structure events, providing new insights into the specific transcript alterations caused by common and rare genetic variants and highlighting the resolution gained from long-read data. We were able to perturb the transcript structure upon knockdown of PTBP1, an RNA binding protein that mediates splicing, thereby finding genetic regulatory effects that are modified by the cellular environment. Finally, we used this dataset to enhance variant interpretation and study rare variants leading to aberrant splicing patterns.

TDP-43 represses cryptic exon inclusion in the FTD-ALS gene UNC13A.

A hallmark pathological feature of the neurodegenerative diseases amyotrophic lateral sclerosis (ALS) and frontotemporal dementia (FTD) is the depletion of RNA-binding protein TDP-43 from the nucleus of neurons in the brain and spinal cord1. A major function of TDP-43 is as a repressor of cryptic exon inclusion during RNA splicing2-4. Single nucleotide polymorphisms in UNC13A are among the strongest hits associated with FTD and ALS in human genome-wide association studies5,6, but how those variants increase risk for disease is unknown. Here we show that TDP-43 represses a cryptic exon-splicing event in UNC13A. Loss of TDP-43 from the nucleus in human brain, neuronal cell lines and motor neurons derived from induced pluripotent stem cells resulted in the inclusion of a cryptic exon in UNC13A mRNA and reduced UNC13A protein expression. The top variants associated with FTD or ALS risk in humans are located in the intron harbouring the cryptic exon, and we show that they increase UNC13A cryptic exon splicing in the face of TDP-43 dysfunction. Together, our data provide a direct functional link between one of the strongest genetic risk factors for FTD and ALS (UNC13A genetic variants), and loss of TDP-43 function.

Evaluating drug targets through human loss-of-function genetic variation.

Naturally occurring human genetic variants that are predicted to inactivate protein-coding genes provide an in vivo model of human gene inactivation that complements knockout studies in cells and model organisms. Here we report three key findings regarding the assessment of candidate drug targets using human loss-of-function variants. First, even essential genes, in which loss-of-function variants are not tolerated, can be highly successful as targets of inhibitory drugs. Second, in most genes, loss-of-function variants are sufficiently rare that genotype-based ascertainment of homozygous or compound heterozygous 'knockout' humans will await sample sizes that are approximately 1,000 times those presently available, unless recruitment focuses on consanguineous individuals. Third, automated variant annotation and filtering are powerful, but manual curation remains crucial for removing artefacts, and is a prerequisite for recall-by-genotype efforts. Our results provide a roadmap for human knockout studies and should guide the interpretation of loss-of-function variants in drug development.

Transcript expression-aware annotation improves rare variant interpretation.

The acceleration of DNA sequencing in samples from patients and population studies has resulted in extensive catalogues of human genetic variation, but the interpretation of rare genetic variants remains problematic. A notable example of this challenge is the existence of disruptive variants in dosage-sensitive disease genes, even in apparently healthy individuals. Here, by manual curation of putative loss-of-function (pLoF) variants in haploinsufficient disease genes in the Genome Aggregation Database (gnomAD)1, we show that one explanation for this paradox involves alternative splicing of mRNA, which allows exons of a gene to be expressed at varying levels across different cell types. Currently, no existing annotation tool systematically incorporates information about exon expression into the interpretation of variants. We develop a transcript-level annotation metric known as the 'proportion expressed across transcripts', which quantifies isoform expression for variants. We calculate this metric using 11,706 tissue samples from the Genotype Tissue Expression (GTEx) project2 and show that it can differentiate between weakly and highly evolutionarily conserved exons, a proxy for functional importance. We demonstrate that expression-based annotation selectively filters 22.8% of falsely annotated pLoF variants found in haploinsufficient disease genes in gnomAD, while removing less than 4% of high-confidence pathogenic variants in the same genes. Finally, we apply our expression filter to the analysis of de novo variants in patients with autism spectrum disorder and intellectual disability or developmental disorders to show that pLoF variants in weakly expressed regions have similar effect sizes to those of synonymous variants, whereas pLoF variants in highly expressed exons are most strongly enriched among cases. Our annotation is fast, flexible and generalizable, making it possible for any variant file to be annotated with any isoform expression dataset, and will be valuable for the genetic diagnosis of rare diseases, the analysis of rare variant burden in complex disorders, and the curation and prioritization of variants in recall-by-genotype studies.

The mutational constraint spectrum quantified from variation in 141,456 humans.

Genetic variants that inactivate protein-coding genes are a powerful source of information about the phenotypic consequences of gene disruption: genes that are crucial for the function of an organism will be depleted of such variants in natural populations, whereas non-essential genes will tolerate their accumulation. However, predicted loss-of-function variants are enriched for annotation errors, and tend to be found at extremely low frequencies, so their analysis requires careful variant annotation and very large sample sizes1. Here we describe the aggregation of 125,748 exomes and 15,708 genomes from human sequencing studies into the Genome Aggregation Database (gnomAD). We identify 443,769 high-confidence predicted loss-of-function variants in this cohort after filtering for artefacts caused by sequencing and annotation errors. Using an improved model of human mutation rates, we classify human protein-coding genes along a spectrum that represents tolerance to inactivation, validate this classification using data from model organisms and engineered human cells, and show that it can be used to improve the power of gene discovery for both common and rare diseases.

Corrigendum: Landscape of X chromosome inactivation across human tissues.

This corrects the article DOI: 10.1038/nature24265.

Landscape of X chromosome inactivation across human tissues.

X chromosome inactivation (XCI) silences transcription from one of the two X chromosomes in female mammalian cells to balance expression dosage between XX females and XY males. XCI is, however, incomplete in humans: up to one-third of X-chromosomal genes are expressed from both the active and inactive X chromosomes (Xa and Xi, respectively) in female cells, with the degree of 'escape' from inactivation varying between genes and individuals. The extent to which XCI is shared between cells and tissues remains poorly characterized, as does the degree to which incomplete XCI manifests as detectable sex differences in gene expression and phenotypic traits. Here we describe a systematic survey of XCI, integrating over 5,500 transcriptomes from 449 individuals spanning 29 tissues from GTEx (v6p release) and 940 single-cell transcriptomes, combined with genomic sequence data. We show that XCI at 683 X-chromosomal genes is generally uniform across human tissues, but identify examples of heterogeneity between tissues, individuals and cells. We show that incomplete XCI affects at least 23% of X-chromosomal genes, identify seven genes that escape XCI with support from multiple lines of evidence and demonstrate that escape from XCI results in sex biases in gene expression, establishing incomplete XCI as a mechanism that is likely to introduce phenotypic diversity. Overall, this updated catalogue of XCI across human tissues helps to increase our understanding of the extent and impact of the incompleteness in the maintenance of XCI.

A quantitative framework for characterizing the evolutionary history of mammalian gene expression.

The evolutionary history of a gene helps predict its function and relationship to phenotypic traits. Although sequence conservation is commonly used to decipher gene function and assess medical relevance, methods for functional inference from comparative expression data are lacking. Here, we use RNA-seq across seven tissues from 17 mammalian species to show that expression evolution across mammals is accurately modeled by the Ornstein-Uhlenbeck process, a commonly proposed model of continuous trait evolution. We apply this model to identify expression pathways under neutral, stabilizing, and directional selection. We further demonstrate novel applications of this model to quantify the extent of stabilizing selection on a gene's expression, parameterize the distribution of each gene's optimal expression level, and detect deleterious expression levels in expression data from individual patients. Our work provides a statistical framework for interpreting expression data across species and in disease.

Analysis of protein-coding genetic variation in 60,706 humans.

Large-scale reference data sets of human genetic variation are critical for the medical and functional interpretation of DNA sequence changes. Here we describe the aggregation and analysis of high-quality exome (protein-coding region) DNA sequence data for 60,706 individuals of diverse ancestries generated as part of the Exome Aggregation Consortium (ExAC). This catalogue of human genetic diversity contains an average of one variant every eight bases of the exome, and provides direct evidence for the presence of widespread mutational recurrence. We have used this catalogue to calculate objective metrics of pathogenicity for sequence variants, and to identify genes subject to strong selection against various classes of mutation; identifying 3,230 genes with near-complete depletion of predicted protein-truncating variants, with 72% of these genes having no currently established human disease phenotype. Finally, we demonstrate that these data can be used for the efficient filtering of candidate disease-causing variants, and for the discovery of human 'knockout' variants in protein-coding genes.

The Genetic Landscape of Diamond-Blackfan Anemia.

Diamond-Blackfan anemia (DBA) is a rare bone marrow failure disorder that affects 7 out of 1,000,000 live births and has been associated with mutations in components of the ribosome. In order to characterize the genetic landscape of this heterogeneous disorder, we recruited a cohort of 472 individuals with a clinical diagnosis of DBA and performed whole-exome sequencing (WES). We identified relevant rare and predicted damaging mutations for 78% of individuals. The majority of mutations were singletons, absent from population databases, predicted to cause loss of function, and located in 1 of 19 previously reported ribosomal protein (RP)-encoding genes. Using exon coverage estimates, we identified and validated 31 deletions in RP genes. We also observed an enrichment for extended splice site mutations and validated their diverse effects using RNA sequencing in cell lines obtained from individuals with DBA. Leveraging the size of our cohort, we observed robust genotype-phenotype associations with congenital abnormalities and treatment outcomes. We further identified rare mutations in seven previously unreported RP genes that may cause DBA, as well as several distinct disorders that appear to phenocopy DBA, including nine individuals with biallelic CECR1 mutations that result in deficiency of ADA2. However, no new genes were identified at exome-wide significance, suggesting that there are no unidentified genes containing mutations readily identified by WES that explain >5% of DBA-affected case subjects. Overall, this report should inform not only clinical practice for DBA-affected individuals, but also the design and analysis of rare variant studies for heterogeneous Mendelian disorders.

Variants in the Oxidoreductase PYROXD1 Cause Early-Onset Myopathy with Internalized Nuclei and Myofibrillar Disorganization.

This study establishes PYROXD1 variants as a cause of early-onset myopathy and uses biospecimens and cell lines, yeast, and zebrafish models to elucidate the fundamental role of PYROXD1 in skeletal muscle. Exome sequencing identified recessive variants in PYROXD1 in nine probands from five families. Affected individuals presented in infancy or childhood with slowly progressive proximal and distal weakness, facial weakness, nasal speech, swallowing difficulties, and normal to moderately elevated creatine kinase. Distinctive histopathology showed abundant internalized nuclei, myofibrillar disorganization, desmin-positive inclusions, and thickened Z-bands. PYROXD1 is a nuclear-cytoplasmic pyridine nucleotide-disulphide reductase (PNDR). PNDRs are flavoproteins (FAD-binding) and catalyze pyridine-nucleotide-dependent (NAD/NADH) reduction of thiol residues in other proteins. Complementation experiments in yeast lacking glutathione reductase glr1 show that human PYROXD1 has reductase activity that is strongly impaired by the disease-associated missense mutations. Immunolocalization studies in human muscle and zebrafish myofibers demonstrate that PYROXD1 localizes to the nucleus and to striated sarcomeric compartments. Zebrafish with ryroxD1 knock-down recapitulate features of PYROXD1 myopathy with sarcomeric disorganization, myofibrillar aggregates, and marked swimming defect. We characterize variants in the oxidoreductase PYROXD1 as a cause of early-onset myopathy with distinctive histopathology and introduce altered redox regulation as a primary cause of congenital muscle disease.

Cryptic intronic NBAS variant reveals the genetic basis of recurrent liver failure in a child.

BACKGROUND: In almost half of patients with acute liver failure the cause is unknown, making targeted treatment and decisions about liver transplantation a challenge. Monogenic disorders may contribute to a significant proportion of these undiagnosed patients, and so the incorporation of technologies such as next generation sequencing (NGS) in the clinic could aid in providing a definitive diagnosis. However, this technology may present a major challenge in interpretation of sequence variants, particularly those in non-coding regions. RESULTS: In this report we describe a case of Infantile liver failure syndrome 2 (ILFS2; MIM 616483) due to novel bi-allelic variants in the NBAS gene. A missense variant NM_015909.3(NBAS):c.2617C > T, NP_056993.2(NBAS):p.(Arg873Trp) was identified by whole genome sequencing (WGS). By combining WGS and reverse transcription-polymerase chain reaction (RT-PCR) we were able to identify a novel deep intronic variant, NM_015909.3(NBAS):c.2423 + 404G > C, leading to the inclusion of a pseudo-exon. This mechanism has not been described previously in this syndrome. CONCLUSIONS: This study highlights the utility of analyzing NGS data in conjunction with investigating complementary DNA (cDNA) using techniques such as RT-PCR for detection of variants that otherwise would be likely to be missed in common NGS bioinformatic analysis pipelines. Combining these approaches, particularly when the phenotype match is strong, could lead to an increase in the diagnostic yield in acute liver failure and thus aid in targeted treatment, accurate genetic counseling and restoration of reproductive confidence.

The ExAC browser: displaying reference data information from over 60 000 exomes.

Worldwide, hundreds of thousands of humans have had their genomes or exomes sequenced, and access to the resulting data sets can provide valuable information for variant interpretation and understanding gene function. Here, we present a lightweight, flexible browser framework to display large population datasets of genetic variation. We demonstrate its use for exome sequence data from 60 706 individuals in the Exome Aggregation Consortium (ExAC). The ExAC browser provides gene- and transcript-centric displays of variation, a critical view for clinical applications. Additionally, we provide a variant display, which includes population frequency and functional annotation data as well as short read support for the called variant. This browser is open-source, freely available at http://exac.broadinstitute.org, and has already been used extensively by clinical laboratories worldwide.

Nemaline myopathy and distal arthrogryposis associated with an autosomal recessive TNNT3 splice variant.

A male neonate presented with severe weakness, hypotonia, contractures and congenital scoliosis. Skeletal muscle specimens showed marked atrophy and degeneration of fast fibers with striking nemaline rods and hypertrophy of slow fibers that were ultrastructurally normal. A neuromuscular gene panel identified a homozygous essential splice variant in TNNT3 (chr11:1956150G > A, NM_006757.3:c.681+1G > A). TNNT3 encodes skeletal troponin-Tfast and is associated with autosomal dominant distal arthrogryposis. TNNT3 has not previously been associated with nemaline myopathy (NM), a rare congenital myopathy linked to defects in proteins associated with thin filament structure and regulation. cDNA studies confirmed pathogenic consequences of the splice variant, eliciting exon-skipping and intron retention events leading to a frameshift. Western blot showed deficiency of troponin-Tfast protein with secondary loss of troponin-Ifast . We establish a homozygous splice variant in TNNT3 as the likely cause of severe congenital NM with distal arthrogryposis, characterized by specific involvement of Type-2 fibers and deficiency of troponin-Tfast .

Pathogenic ASXL1 somatic variants in reference databases complicate germline variant interpretation for Bohring-Opitz Syndrome.

The clinical interpretation of genetic variants has come to rely heavily on reference population databases such as the Exome Aggregation Consortium (ExAC) database. Pathogenic variants in genes associated with severe, pediatric-onset, highly penetrant, autosomal dominant conditions are assumed to be absent or rare in these databases. Exome sequencing of a 6-year-old female patient with seizures, developmental delay, dysmorphic features, and failure to thrive identified an ASXL1 variant previously reported as causative of Bohring-Opitz syndrome (BOS). Surprisingly, the variant was observed seven times in the ExAC database, presumably in individuals without BOS. Although the BOS phenotype fit, the presence of the variant in reference population databases introduced ambiguity in result interpretation. Review of the literature revealed that acquired somatic mosaicism of ASXL1 variants (including pathogenic variants) during hematopoietic clonal expansion can occur with aging in healthy individuals. We examined all ASXL1 truncating variants in the ExAC database and determined most are likely somatic. Failure to consider somatic mosaicism may lead to the inaccurate assumption that conditions like BOS have reduced penetrance, or the misclassification of potentially pathogenic variants.