Death after High-Dose rAAV9 Gene Therapy in a Patient with Duchenne's Muscular Dystrophy.

We treated a 27-year-old patient with Duchenne's muscular dystrophy (DMD) with recombinant adeno-associated virus (rAAV) serotype 9 containing dSaCas9 (i.e., "dead" Staphylococcus aureus Cas9, in which the Cas9 nuclease activity has been inactivated) fused to VP64; this transgene was designed to up-regulate cortical dystrophin as a custom CRISPR-transactivator therapy. The dose of rAAV used was 1×1014 vector genomes per kilogram of body weight. Mild cardiac dysfunction and pericardial effusion developed, followed by acute respiratory distress syndrome (ARDS) and cardiac arrest 6 days after transgene treatment; the patient died 2 days later. A postmortem examination showed severe diffuse alveolar damage. Expression of transgene in the liver was minimal, and there was no evidence of AAV serotype 9 antibodies or effector T-cell reactivity in the organs. These findings indicate that an innate immune reaction caused ARDS in a patient with advanced DMD treated with high-dose rAAV gene therapy. (Funded by Cure Rare Disease.).

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.

Multi-omics identifies large mitoribosomal subunit instability caused by pathogenic MRPL39 variants as a cause of pediatric onset mitochondrial disease.

MRPL39 encodes one of 52 proteins comprising the large subunit of the mitochondrial ribosome (mitoribosome). In conjunction with 30 proteins in the small subunit, the mitoribosome synthesizes the 13 subunits of the mitochondrial oxidative phosphorylation (OXPHOS) system encoded by mitochondrial Deoxyribonucleic acid (DNA). We used multi-omics and gene matching to identify three unrelated individuals with biallelic variants in MRPL39 presenting with multisystem diseases with severity ranging from lethal, infantile-onset (Leigh syndrome spectrum) to milder with survival into adulthood. Clinical exome sequencing of known disease genes failed to diagnose these patients; however quantitative proteomics identified a specific decrease in the abundance of large but not small mitoribosomal subunits in fibroblasts from the two patients with severe phenotype. Re-analysis of exome sequencing led to the identification of candidate single heterozygous variants in mitoribosomal genes MRPL39 (both patients) and MRPL15. Genome sequencing identified a shared deep intronic MRPL39 variant predicted to generate a cryptic exon, with transcriptomics and targeted studies providing further functional evidence for causation. The patient with the milder disease was homozygous for a missense variant identified through trio exome sequencing. Our study highlights the utility of quantitative proteomics in detecting protein signatures and in characterizing gene-disease associations in exome-unsolved patients. We describe Relative Complex Abundance analysis of proteomics data, a sensitive method that can identify defects in OXPHOS disorders to a similar or greater sensitivity to the traditional enzymology. Relative Complex Abundance has potential utility for functional validation or prioritization in many hundreds of inherited rare diseases where protein complex assembly is disrupted.

Neuromuscular disorders: finding the missing genetic diagnoses.

Neuromuscular disorders (NMDs) are a wide-ranging group of diseases that seriously affect the quality of life of affected individuals. The development of next-generation sequencing revolutionized the diagnosis of NMD, enabling the discovery of hundreds of NMD genes and many more pathogenic variants. However, the diagnostic yield of genetic testing in NMD cohorts remains incomplete, indicating a large number of genetic diagnoses are not identified through current methods. Fortunately, recent advancements in sequencing technologies, analytical tools, and high-throughput functional screening provide an opportunity to circumvent current challenges. Here, we discuss reasons for missing genetic diagnoses in NMD, how emerging technologies and tools can overcome these hurdles, and examine future approaches to improving diagnostic yields in NMD.

Mitochondrial DNA variation across 56,434 individuals in gnomAD.

Genomic databases of allele frequency are extremely helpful for evaluating clinical variants of unknown significance; however, until now, databases such as the Genome Aggregation Database (gnomAD) have focused on nuclear DNA and have ignored the mitochondrial genome (mtDNA). Here, we present a pipeline to call mtDNA variants that addresses three technical challenges: (1) detecting homoplasmic and heteroplasmic variants, present, respectively, in all or a fraction of mtDNA molecules; (2) circular mtDNA genome; and (3) misalignment of nuclear sequences of mitochondrial origin (NUMTs). We observed that mtDNA copy number per cell varied across gnomAD cohorts and influenced the fraction of NUMT-derived false-positive variant calls, which can account for the majority of putative heteroplasmies. To avoid false positives, we excluded contaminated samples, cell lines, and samples prone to NUMT misalignment due to few mtDNA copies. Furthermore, we report variants with heteroplasmy ≥10%. We applied this pipeline to 56,434 whole-genome sequences in the gnomAD v3.1 database that includes individuals of European (58%), African (25%), Latino (10%), and Asian (5%) ancestry. Our gnomAD v3.1 release contains population frequencies for 10,850 unique mtDNA variants at more than half of all mtDNA bases. Importantly, we report frequencies within each nuclear ancestral population and mitochondrial haplogroup. Homoplasmic variants account for most variant calls (98%) and unique variants (85%). We observed that 1/250 individuals carry a pathogenic mtDNA variant with heteroplasmy above 10%. These mtDNA population allele frequencies are freely accessible and will aid in diagnostic interpretation and research studies.

The clinical and genetic spectrum of autosomal-recessive TOR1A-related disorders.

In the field of rare diseases, progress in molecular diagnostics led to the recognition that variants linked to autosomal-dominant neurodegenerative diseases of later onset can, in the context of biallelic inheritance, cause devastating neurodevelopmental disorders and infantile or childhood-onset neurodegeneration. TOR1A-associated arthrogryposis multiplex congenita 5 (AMC5) is a rare neurodevelopmental disorder arising from biallelic variants in TOR1A, a gene that in the heterozygous state is associated with torsion dystonia-1 (DYT1 or DYT-TOR1A), an early-onset dystonia with reduced penetrance. While 15 individuals with AMC5-TOR1A have been reported (less than 10 in detail), a systematic investigation of the full disease-associated spectrum has not been conducted. Here, we assess the clinical, radiological and molecular characteristics of 57 individuals from 40 families with biallelic variants in TOR1A. Median age at last follow-up was 3 years (0-24 years). Most individuals presented with severe congenital flexion contractures (95%) and variable developmental delay (79%). Motor symptoms were reported in 79% and included lower limb spasticity and pyramidal signs, as well as gait disturbances. Facial dysmorphism was an integral part of the phenotype, with key features being a broad/full nasal tip, narrowing of the forehead and full cheeks. Analysis of disease-associated manifestations delineated a phenotypic spectrum ranging from normal cognition and mild gait disturbance to congenital arthrogryposis, global developmental delay, intellectual disability, absent speech and inability to walk. In a subset, the presentation was consistent with foetal akinesia deformation sequence with severe intrauterine abnormalities. Survival was 71%, with higher mortality in males. Death occurred at a median age of 1.2 months (1 week-9 years), due to respiratory failure, cardiac arrest or sepsis. Analysis of brain MRI studies identified non-specific neuroimaging features, including a hypoplastic corpus callosum (72%), foci of signal abnormality in the subcortical and periventricular white matter (55%), diffuse white matter volume loss (45%), mega cisterna magna (36%) and arachnoid cysts (27%). The molecular spectrum included 22 distinct variants, defining a mutational hotspot in the C-terminal domain of the Torsin-1A protein. Genotype-phenotype analysis revealed an association of missense variants in the 3-helix bundle domain to an attenuated phenotype, while missense variants near the Walker A/B motif as well as biallelic truncating variants were linked to early death. In summary, this systematic cross-sectional analysis of a large cohort of individuals with biallelic TOR1A variants across a wide age-range delineates the clinical and genetic spectrum of TOR1A-related autosomal-recessive disease and highlights potential predictors for disease severity and survival.

MitoVisualize: a resource for analysis of variants in human mitochondrial RNAs and DNA.

SUMMARY: We present MitoVisualize, a new tool for analysis of the human mitochondrial DNA (mtDNA). MitoVisualize enables visualization of: (i) the position and effect of variants in mitochondrial transfer RNA and ribosomal RNA secondary structures alongside curated variant annotations, (ii) data across RNA structures, such as to show all positions with disease-associated variants or with post-transcriptional modifications and (iii) the position of a base, gene or region in the circular mtDNA map, such as to show the location of a large deletion. All visualizations can be easily downloaded as figures for reuse. MitoVisualize can be useful for anyone interested in exploring mtDNA variation, though is designed to facilitate mtDNA variant interpretation in particular. AVAILABILITY AND IMPLEMENTATION: MitoVisualize can be accessed via https://www.mitovisualize.org/. The source code is available at https://github.com/leklab/mito_visualize/. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Corrigendum: High-throughput discovery of novel developmental phenotypes.

This corrects the article DOI: 10.1038/nature19356.

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.

Centers for Mendelian Genomics: A decade of facilitating gene discovery.

PURPOSE: Mendelian disease genomic research has undergone a massive transformation over the past decade. With increasing availability of exome and genome sequencing, the role of Mendelian research has expanded beyond data collection, sequencing, and analysis to worldwide data sharing and collaboration. METHODS: Over the past 10 years, the National Institutes of Health-supported Centers for Mendelian Genomics (CMGs) have played a major role in this research and clinical evolution. RESULTS: We highlight the cumulative gene discoveries facilitated by the program, biomedical research leveraged by the approach, and the larger impact on the research community. Beyond generating a list of gene-phenotype relationships and participating in widespread data sharing, the CMGs have created resources, tools, and training for the larger community to foster understanding of genes and genome variation. The CMGs have participated in a wide range of data sharing activities, including deposition of all eligible CMG data into the Analysis, Visualization, and Informatics Lab-space (AnVIL), sharing candidate genes through the Matchmaker Exchange and the CMG website, and sharing variants in Genotypes to Mendelian Phenotypes (Geno2MP) and VariantMatcher. CONCLUSION: The work is far from complete; strengthening communication between research and clinical realms, continued development and sharing of knowledge and tools, and improving access to richly characterized data sets are all required to diagnose the remaining molecularly undiagnosed patients.

A new phenotype of syndromic retinitis pigmentosa with myopathy is caused by mutations in retinol dehydrogenase 11.

Retinol dehydrogenase 11 (RDH11) is an 11-cis-retinol dehydrogenase that has a well-characterized, albeit auxiliary role in the retinoid cycle. Diseases caused by mutations in the RDH11 gene are very rare, and only one affected family with eye and intelligence involvement has been reported. In the present study, we describe the clinical and genetic findings in a Chinese patient with retinitis pigmentosa (RP), juvenile cataracts, intellectual disability, and myopathy. Trio-based whole-exome sequencing and whole genomic copy number variation detection were performed in this family, and compound heterozygous mutations were identified in RDH11 of the patient: c.938T>C (p.Leu313Pro) derived from the father and c.75-3C>A derived from the mother. Variant c.75-3C>A was confirmed to be a splice-site mutation by cDNA sequencing. It caused exon 2 skipping, resulting in a frameshift mutation and premature translation termination (p.Lys26Serfs*38). Moreover, we found mislocalization of RDH11 protein in muscle cells of the patient by using immunofluorescence staining. This is the first case reported in the Chinese population harboring mutations in RDH11 and revealing a new phenotype of syndromic RP with myopathy.

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.

High-throughput discovery of novel developmental phenotypes.

Approximately one-third of all mammalian genes are essential for life. Phenotypes resulting from knockouts of these genes in mice have provided tremendous insight into gene function and congenital disorders. As part of the International Mouse Phenotyping Consortium effort to generate and phenotypically characterize 5,000 knockout mouse lines, here we identify 410 lethal genes during the production of the first 1,751 unique gene knockouts. Using a standardized phenotyping platform that incorporates high-resolution 3D imaging, we identify phenotypes at multiple time points for previously uncharacterized genes and additional phenotypes for genes with previously reported mutant phenotypes. Unexpectedly, our analysis reveals that incomplete penetrance and variable expressivity are common even on a defined genetic background. In addition, we show that human disease genes are enriched for essential genes, thus providing a dataset that facilitates the prioritization and validation of mutations identified in clinical sequencing efforts.

Correction: Insights into the genetic epidemiology of Crohn's and rare diseases in the Ashkenazi Jewish population.

[This corrects the article DOI: 10.1371/journal.pgen.1007329.].

The Effect of ACTN3 Gene Doping on Skeletal Muscle Performance.

Loss of expression of ACTN3, due to homozygosity of the common null polymorphism (p.Arg577X), is underrepresented in elite sprint/power athletes and has been associated with reduced muscle mass and strength in humans and mice. To investigate ACTN3 gene dosage in performance and whether expression could enhance muscle force, we performed meta-analysis and expression studies. Our general meta-analysis using a Bayesian random effects model in elite sprint/power athlete cohorts demonstrated a consistent homozygous-group effect across studies (per allele OR = 1.4, 95% CI 1.3-1.6) but substantial heterogeneity in heterozygotes. In mouse muscle, rAAV-mediated gene transfer overexpressed and rescued α-actinin-3 expression. Contrary to expectation, in vivo "doping" of ACTN3 at low to moderate doses demonstrated an absence of any change in function. At high doses, ACTN3 is toxic and detrimental to force generation, to demonstrate gene doping with supposedly performance-enhancing isoforms of sarcomeric proteins can be detrimental for muscle function. Restoration of α-actinin-3 did not enhance muscle mass but highlighted the primary role of α-actinin-3 in modulating muscle metabolism with altered fatiguability. This is the first study to express a Z-disk protein in healthy skeletal muscle and measure the in vivo effect. The sensitive balance of the sarcomeric proteins and muscle function has relevant implications in areas of gene doping in performance and therapy for neuromuscular disease.

Variants in EXOSC9 Disrupt the RNA Exosome and Result in Cerebellar Atrophy with Spinal Motor Neuronopathy.

The exosome is a conserved multi-protein complex that is essential for correct RNA processing. Recessive variants in exosome components EXOSC3, EXOSC8, and RBM7 cause various constellations of pontocerebellar hypoplasia (PCH), spinal muscular atrophy (SMA), and central nervous system demyelination. Here, we report on four unrelated affected individuals with recessive variants in EXOSC9 and the effect of the variants on the function of the RNA exosome in vitro in affected individuals' fibroblasts and skeletal muscle and in vivo in zebrafish. The clinical presentation was severe, early-onset, progressive SMA-like motor neuronopathy, cerebellar atrophy, and in one affected individual, congenital fractures of the long bones. Three affected individuals of different ethnicity carried the homozygous c.41T>C (p.Leu14Pro) variant, whereas one affected individual was compound heterozygous for c.41T>C (p.Leu14Pro) and c.481C>T (p.Arg161∗). We detected reduced EXOSC9 in fibroblasts and skeletal muscle and observed a reduction of the whole multi-subunit exosome complex on blue-native polyacrylamide gel electrophoresis. RNA sequencing of fibroblasts and skeletal muscle detected significant >2-fold changes in genes involved in neuronal development and cerebellar and motor neuron degeneration, demonstrating the widespread effect of the variants. Morpholino oligonucleotide knockdown and CRISPR/Cas9-mediated mutagenesis of exosc9 in zebrafish recapitulated aspects of the human phenotype, as they have in other zebrafish models of exosomal disease. Specifically, portions of the cerebellum and hindbrain were absent, and motor neurons failed to develop and migrate properly. In summary, we show that variants in EXOSC9 result in a neurological syndrome combining cerebellar atrophy and spinal motoneuronopathy, thus expanding the list of human exosomopathies.

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.

Insights into the genetic epidemiology of Crohn's and rare diseases in the Ashkenazi Jewish population.

As part of a broader collaborative network of exome sequencing studies, we developed a jointly called data set of 5,685 Ashkenazi Jewish exomes. We make publicly available a resource of site and allele frequencies, which should serve as a reference for medical genetics in the Ashkenazim (hosted in part at https://ibd.broadinstitute.org, also available in gnomAD at http://gnomad.broadinstitute.org). We estimate that 34% of protein-coding alleles present in the Ashkenazi Jewish population at frequencies greater than 0.2% are significantly more frequent (mean 15-fold) than their maximum frequency observed in other reference populations. Arising via a well-described founder effect approximately 30 generations ago, this catalog of enriched alleles can contribute to differences in genetic risk and overall prevalence of diseases between populations. As validation we document 148 AJ enriched protein-altering alleles that overlap with "pathogenic" ClinVar alleles (table available at https://github.com/macarthur-lab/clinvar/blob/master/output/clinvar.tsv), including those that account for 10-100 fold differences in prevalence between AJ and non-AJ populations of some rare diseases, especially recessive conditions, including Gaucher disease (GBA, p.Asn409Ser, 8-fold enrichment); Canavan disease (ASPA, p.Glu285Ala, 12-fold enrichment); and Tay-Sachs disease (HEXA, c.1421+1G>C, 27-fold enrichment; p.Tyr427IlefsTer5, 12-fold enrichment). We next sought to use this catalog, of well-established relevance to Mendelian disease, to explore Crohn's disease, a common disease with an estimated two to four-fold excess prevalence in AJ. We specifically attempt to evaluate whether strong acting rare alleles, particularly protein-truncating or otherwise large effect-size alleles, enriched by the same founder-effect, contribute excess genetic risk to Crohn's disease in AJ, and find that ten rare genetic risk factors in NOD2 and LRRK2 are enriched in AJ (p < 0.005), including several novel contributing alleles, show evidence of association to CD. Independently, we find that genomewide common variant risk defined by GWAS shows a strong difference between AJ and non-AJ European control population samples (0.97 s.d. higher, p<10-16). Taken together, the results suggest coordinated selection in AJ population for higher CD risk alleles in general. The results and approach illustrate the value of exome sequencing data in case-control studies along with reference data sets like ExAC (sites VCF available via FTP at ftp.broadinstitute.org/pub/ExAC_release/release0.3/) to pinpoint genetic variation that contributes to variable disease predisposition across populations.