Systematic evaluation of genome sequencing for the diagnostic assessment of autism spectrum disorder and fetal structural anomalies.

Short-read genome sequencing (GS) holds the promise of becoming the primary diagnostic approach for the assessment of autism spectrum disorder (ASD) and fetal structural anomalies (FSAs). However, few studies have comprehensively evaluated its performance against current standard-of-care diagnostic tests: karyotype, chromosomal microarray (CMA), and exome sequencing (ES). To assess the clinical utility of GS, we compared its diagnostic yield against these three tests in 1,612 quartet families including an individual with ASD and in 295 prenatal families. Our GS analytic framework identified a diagnostic variant in 7.8% of ASD probands, almost 2-fold more than CMA (4.3%) and 3-fold more than ES (2.7%). However, when we systematically captured copy-number variants (CNVs) from the exome data, the diagnostic yield of ES (7.4%) was brought much closer to, but did not surpass, GS. Similarly, we estimated that GS could achieve an overall diagnostic yield of 46.1% in unselected FSAs, representing a 17.2% increased yield over karyotype, 14.1% over CMA, and 4.1% over ES with CNV calling or 36.1% increase without CNV discovery. Overall, GS provided an added diagnostic yield of 0.4% and 0.8% beyond the combination of all three standard-of-care tests in ASD and FSAs, respectively. This corresponded to nine GS unique diagnostic variants, including sequence variants in exons not captured by ES, structural variants (SVs) inaccessible to existing standard-of-care tests, and SVs where the resolution of GS changed variant classification. Overall, this large-scale evaluation demonstrated that GS significantly outperforms each individual standard-of-care test while also outperforming the combination of all three tests, thus warranting consideration as the first-tier diagnostic approach for the assessment of ASD and FSAs.

Harmonizing variant classification for return of results in the All of Us Research Program.

The All of Us Research Program (AoURP) is a historic effort to accelerate research and improve healthcare by generating and collating data from one million people in the United States. Participants will have the option to receive results from their genome analysis, including actionable findings in 59 gene-disorder pairs for which disorder-associated variants are recommended for return by the American College of Medical Genetics and Genomics. To ensure consistent reporting across the AoURP, in a prelaunch study the four participating clinical laboratories shared all variant classifications in the 59 genes of interest from their internal databases. Of the 11,813 unique variants classified by at least two of the four laboratories, classifications were concordant with regard to reportability for 99.1% (11,711), with only 0.9% (102) having reportability differences. Through variant reassessment, data sharing, and discussion of rationale, participating laboratories resolved all 102 reportable differences. These approaches will be maintained during routine AoU reporting to ensure continuous classification harmonization and consistent reporting within AoURP.

seqr: A web-based analysis and collaboration tool for rare disease genomics.

Exome and genome sequencing have become the tools of choice for rare disease diagnosis, leading to large amounts of data available for analyses. To identify causal variants in these datasets, powerful filtering and decision support tools that can be efficiently used by clinicians and researchers are required. To address this need, we developed seqr - an open-source, web-based tool for family-based monogenic disease analysis that allows researchers to work collaboratively to search and annotate genomic callsets. To date, seqr is being used in several research pipelines and one clinical diagnostic lab. In our own experience through the Broad Institute Center for Mendelian Genomics, seqr has enabled analyses of over 10,000 families, supporting the diagnosis of more than 3,800 individuals with rare disease and discovery of over 300 novel disease genes. Here, we describe a framework for genomic analysis in rare disease that leverages seqr's capabilities for variant filtration, annotation, and causal variant identification, as well as support for research collaboration and data sharing. The seqr platform is available as open source software, allowing low-cost participation in rare disease research, and a community effort to support diagnosis and gene discovery in rare disease.

Randomized prospective evaluation of genome sequencing versus standard-of-care as a first molecular diagnostic test.

PURPOSE: To evaluate the diagnostic yield and clinical relevance of clinical genome sequencing (cGS) as a first genetic test for patients with suspected monogenic disorders. METHODS: We conducted a prospective randomized study with pediatric and adult patients recruited from genetics clinics at Massachusetts General Hospital who were undergoing planned genetic testing. Participants were randomized into two groups: standard-of-care genetic testing (SOC) only or SOC and cGS. RESULTS: Two hundred four participants were enrolled, 202 were randomized to one of the intervention arms, and 99 received cGS. In total, cGS returned 16 molecular diagnoses that fully or partially explained the indication for testing in 16 individuals (16.2% of the cohort, 95% confidence interval [CI] 8.9-23.4%), which was not significantly different from SOC (18.2%, 95% CI 10.6-25.8%, P = 0.71). An additional eight molecular diagnoses reported by cGS had uncertain relevance to the participant's phenotype. Nevertheless, referring providers considered 20/24 total cGS molecular diagnoses (83%) to be explanatory for clinical features or worthy of additional workup. CONCLUSION: cGS is technically suitable as a first genetic test. In our cohort, diagnostic yield was not significantly different from SOC. Further studies addressing other variant types and implementation challenges are needed to support feasibility and utility of broad-scale cGS adoption.

Mutation mapping and identification by whole-genome sequencing.

Genetic mapping of mutations in model systems has facilitated the identification of genes contributing to fundamental biological processes including human diseases. However, this approach has historically required the prior characterization of informative markers. Here we report a fast and cost-effective method for genetic mapping using next-generation sequencing that combines single nucleotide polymorphism discovery, mutation localization, and potential identification of causal sequence variants. In contrast to prior approaches, we have developed a hidden Markov model to narrowly define the mutation area by inferring recombination breakpoints of chromosomes in the mutant pool. In addition, we created an interactive online software resource to facilitate automated analysis of sequencing data and demonstrate its utility in the zebrafish and mouse models. Our novel methodology and online tools will make next-generation sequencing an easily applicable resource for mutation mapping in all model systems.

Rare ACTN2 Frameshift Variants Resulting in Protein Extension Cause Distal Myopathy and Hypertrophic Cardiomyopathy through Protein Aggregation.

Distal myopathies are a group of rare, inherited muscular disorders characterized by progressive loss of muscle fibers that begins in the distal parts of arms and legs. Recently, variants in a new disease gene, ACTN2 , have been shown to cause distal myopathy. ACTN2 , a gene previously only associated with cardiomyopathies, encodes alpha-actinin-2, a protein expressed in both cardiac and skeletal sarcomeres. The primary function of alpha-actinin-2 is to link actin and titin to the sarcomere Z-disk. New ACTN2 variants are continuously discovered, however, the clinical significance of many variants remains unknown. Thus, lack of clear genotype-phenotype correlations in ACTN2 -related diseases, actininopathies, persists. Objective: The objective of the study is to characterize the pathomechanisms underlying actininopathies. Methods: Functional characterization in C2C12 cell models of several ACTN2 variants is conducted, including frameshift and missense variants associated with dominant actininopathies. We assess the genotype-phenotype correlations of actininopathies using clinical data from several patients carrying these variants. Results: The results show that the missense variants associated with a recessive form of actininopathy do not cause detectable alpha-actinin-2 aggregates in the cell model. Conversely, dominant frameshift variants causing a protein extension do produce alpha-actinin-2 aggregates. Interpretation: The results suggest that alpha-actinin-2 aggregation is the disease mechanism underlying some dominant actininopathies, and thus we recommend that protein-extending frameshift variants in ACTN2 should be classified as pathogenic. However, this mechanism is likely elicited by only a limited number of variants. Alternative functional characterization methods should be explored to further investigate other molecular mechanisms underlying actininopathies.

Best practices for the interpretation and reporting of clinical whole genome sequencing.

Whole genome sequencing (WGS) shows promise as a first-tier diagnostic test for patients with rare genetic disorders. However, standards addressing the definition and deployment practice of a best-in-class test are lacking. To address these gaps, the Medical Genome Initiative, a consortium of leading health care and research organizations in the US and Canada, was formed to expand access to high quality clinical WGS by convening experts and publishing best practices. Here, we present best practice recommendations for the interpretation and reporting of clinical diagnostic WGS, including discussion of challenges and emerging approaches that will be critical to harness the full potential of this comprehensive test.

Endophenotype effect sizes support variant pathogenicity in monogenic disease susceptibility genes.

Accurate and efficient classification of variant pathogenicity is critical for research and clinical care. Using data from three large studies, we demonstrate that population-based associations between rare variants and quantitative endophenotypes for three monogenic diseases (low-density-lipoprotein cholesterol for familial hypercholesterolemia, electrocardiographic QTc interval for long QT syndrome, and glycosylated hemoglobin for maturity-onset diabetes of the young) provide evidence for variant pathogenicity. Effect sizes are associated with pathogenic ClinVar assertions (P < 0.001 for each trait) and discriminate pathogenic from non-pathogenic variants (area under the curve 0.82-0.84 across endophenotypes). An effect size threshold of ≥ 0.5 times the endophenotype standard deviation nominates up to 35% of rare variants of uncertain significance or not in ClinVar in disease susceptibility genes with pathogenic potential. We propose that variant associations with quantitative endophenotypes for monogenic diseases can provide evidence supporting pathogenicity.

Analysis of intragenic USH2A copy number variation unveils broad spectrum of unique and recurrent variants.

Given that all forms of Usher syndrome (USH) present with hearing loss in advance of retinal disease, the syndromic nature of the disorder is rarely appreciated when critical management decisions are being made. As a result, molecular diagnostics are crucial in guiding the management of USH patients. While 11 genes have been associated with USH, the USH2A gene is one of the largest contributors. Approximately 20% of suspected USH probands that undergo USH2A sequencing at our laboratory receive an inconclusive result due to the identification of a monoallelic disease-causing variant in USH2A. Many studies suggest that intragenic deletions and duplications represent an important USH2A variant type that can be missed by sequencing assays if supplemental algorithms or testing methods are not applied. To gain a comprehensive view of the contribution of USH2A CNVs to USH, we conducted prospective and retrospective screening in 700 hearing loss probands. Fourteen individuals with 11 unique USH2A CNVs are reported, including one pathogenic multi-exon duplication. Additionally, we mapped deletion breakpoints and performed a meta-analysis of USH2A CNVs to evaluate recurrence and underlying mechanisms. This analysis revealed breakpoint grouping within three introns, raising the possibility of CNV-susceptible regions within the gene. Overall, our data highlight the diversity of pathogenic CNVs in this gene, demonstrating that the comprehensive, high-resolution USH2A CNV analysis methods employed here are essential components of clinical genetic testing for USH.

Distribution and clinical impact of functional variants in 50,726 whole-exome sequences from the DiscovEHR study.

The DiscovEHR collaboration between the Regeneron Genetics Center and Geisinger Health System couples high-throughput sequencing to an integrated health care system using longitudinal electronic health records (EHRs). We sequenced the exomes of 50,726 adult participants in the DiscovEHR study to identify ~4.2 million rare single-nucleotide variants and insertion/deletion events, of which ~176,000 are predicted to result in a loss of gene function. Linking these data to EHR-derived clinical phenotypes, we find clinical associations supporting therapeutic targets, including genes encoding drug targets for lipid lowering, and identify previously unidentified rare alleles associated with lipid levels and other blood level traits. About 3.5% of individuals harbor deleterious variants in 76 clinically actionable genes. The DiscovEHR data set provides a blueprint for large-scale precision medicine initiatives and genomics-guided therapeutic discovery.

Cytoplasmic carboxypeptidase 5 regulates tubulin glutamylation and zebrafish cilia formation and function.

Glutamylation is a functionally important tubulin posttranslational modification enriched on stable microtubules of neuronal axons, mitotic spindles, centrioles, and cilia. In vertebrates, balanced activities of tubulin glutamyl ligase and cytoplasmic carboxypeptidase deglutamylase enzymes maintain organelle- and cell type-specific tubulin glutamylation patterns. Tubulin glutamylation in cilia is regulated via restricted subcellular localization or expression of tubulin glutamyl ligases (ttlls) and nonenzymatic proteins, including the zebrafish TPR repeat protein Fleer/Ift70. Here we analyze the expression patterns of ccp deglutamylase genes during zebrafish development and the effects of ccp gene knockdown on cilia formation, morphology, and tubulin glutamylation. The deglutamylases ccp2, ccp5, and ccp6 are expressed in ciliated cells, whereas ccp1 expression is restricted to the nervous system. Only ccp5 knockdown increases cilia tubulin glutamylation, induces ciliopathy phenotypes, including axis curvature, hydrocephalus, and pronephric cysts, and disrupts multicilia motility, suggesting that Ccp5 is the principal tubulin deglutamylase that maintains functional levels of cilia tubulin glutamylation. The ability of ccp5 knockdown to restore cilia tubulin glutamylation in fleer/ift70 mutants and rescue pronephric multicilia formation in both fleer- and ift88-deficient zebrafish indicates that tubulin glutamylation is a key driver of ciliogenesis.