Assessing variants of uncertain significance implicated in hearing loss using a comprehensive deafness proteome.

Hearing loss is the leading sensory deficit, affecting ~ 5% of the population. It exhibits remarkable heterogeneity across 223 genes with 6328 pathogenic missense variants, making deafness-specific expertise a prerequisite for ascribing phenotypic consequences to genetic variants. Deafness-implicated variants are curated in the Deafness Variation Database (DVD) after classification by a genetic hearing loss expert panel and thorough informatics pipeline. However, seventy percent of the 128,167 missense variants in the DVD are "variants of uncertain significance" (VUS) due to insufficient evidence for classification. Here, we use the deep learning protein prediction algorithm, AlphaFold2, to curate structures for all DVD genes. We refine these structures with global optimization and the AMOEBA force field and use DDGun3D to predict folding free energy differences (∆∆GFold) for all DVD missense variants. We find that 5772 VUSs have a large, destabilizing ∆∆GFold that is consistent with pathogenic variants. When also filtered for CADD scores (> 25.7), we determine 3456 VUSs are likely pathogenic at a probability of 99.0%. Of the 224 genes in the DVD, 166 genes (74%) exhibit one or more missense variants predicted to cause a pathogenic change in protein folding stability. The VUSs prioritized here affect 119 patients (~ 3% of cases) sequenced by the OtoSCOPE targeted panel. Approximately half of these patients previously received an inconclusive report, and reclassification of these VUSs as pathogenic provides a new genetic diagnosis for six patients.

De novo variants are a common cause of genetic hearing loss.

PURPOSE: De novo variants (DNVs) are a well-recognized cause of genetic disorders. The contribution of DNVs to hearing loss (HL) is poorly characterized. We aimed to evaluate the rate of DNVs in HL-associated genes and assess their contribution to HL. METHODS: Targeted genomic enrichment and massively parallel sequencing were used for molecular testing of all exons and flanking intronic sequences of known HL-associated genes, with no exclusions on the basis of type of HL or clinical features. Segregation analysis was performed, and previous reports of DNVs in PubMed and ClinVar were reviewed to characterize the rate, distribution, and spectrum of DNVs in HL. RESULTS: DNVs were detected in 10% (24/238) of trios for whom segregation analysis was performed. Overall, DNVs were causative in at least ∼1% of probands for whom a genetic diagnosis was resolved, with marked variability based on inheritance mode and phenotype. DNVs of MITF were most common (21% of DNVs), followed by GATA3 (13%), STRC (13%), and ACTG1 (8%). Review of reported DNVs revealed gene-specific variability in contribution of DNV to the mutational spectrum of HL-associated genes. CONCLUSION: DNVs are a relatively common cause of genetic HL and must be considered in all cases of sporadic HL.

A comparative analysis of genetic hearing loss phenotypes in European/American and Japanese populations.

We present detailed comparative analyses to assess population-level differences in patterns of genetic deafness between European/American and Japanese cohorts with non-syndromic hearing loss. One thousand eighty-three audiometric test results (921 European/American and 162 Japanese) from members of 168 families (48 European/American and 120 Japanese) with non-syndromic hearing loss secondary to pathogenic variants in one of three genes (KCNQ4, TECTA, WFS1) were studied. Audioprofile characteristics, specific mutation types, and protein domains were considered in the comparative analyses. Our findings support differences in audioprofiles driven by both mutation type (non-truncating vs. truncating) and ethnic background. The former finding confirms data that ascribe a phenotypic consequence to different mutation types in KCNQ4; the latter finding suggests that there are ethnic-specific effects (genetic and/or environmental) that impact gene-specific audioprofiles for TECTA and WFS1. Identifying the drivers of ethnic differences will refine our understanding of phenotype-genotype relationships and the biology of hearing and deafness.

High-Throughput Genetic Testing for Thrombotic Microangiopathies and C3 Glomerulopathies.

The thrombotic microangiopathies (TMAs) and C3 glomerulopathies (C3Gs) include a spectrum of rare diseases such as atypical hemolytic uremic syndrome, thrombotic thrombocytopenic purpura, C3GN, and dense deposit disease, which share phenotypic similarities and underlying genetic commonalities. Variants in several genes contribute to the pathogenesis of these diseases, and identification of these variants may inform the diagnosis and treatment of affected patients. We have developed and validated a comprehensive genetic panel that screens all exons of all genes implicated in TMA and C3G. The closely integrated pipeline implemented includes targeted genomic enrichment, massively parallel sequencing, bioinformatic analysis, and a multidisciplinary conference to analyze identified variants in the context of each patient's specific phenotype. Herein, we present our 1-year experience with this panel, during which time we studied 193 patients. We identified 17 novel and 74 rare variants, which we classified as pathogenic (11), likely pathogenic (12), and of uncertain significance (68). Compared with controls, patients with C3G had a higher frequency of rare and novel variants in C3 convertase (C3 and CFB) and complement regulator (CFH, CFI, CFHR5, and CD46) genes (P<0.05). In contrast, patients with TMA had an increase in rare and novel variants only in complement regulator genes (P<0.01), a distinction consistent with differing sites of complement dysregulation in these two diseases. In summary, we were able to provide a positive genetic diagnosis in 43% and 41% of patients carrying the clinical diagnosis of C3G and TMA, respectively.

Comprehensive genetic testing in the clinical evaluation of 1119 patients with hearing loss.

Hearing loss is the most common sensory deficit in humans, affecting 1 in 500 newborns. Due to its genetic heterogeneity, comprehensive diagnostic testing has not previously been completed in a large multiethnic cohort. To determine the aggregate contribution inheritance makes to non-syndromic hearing loss, we performed comprehensive clinical genetic testing with targeted genomic enrichment and massively parallel sequencing on 1119 sequentially accrued patients. No patient was excluded based on phenotype, inheritance or previous testing. Testing resulted in identification of the underlying genetic cause for hearing loss in 440 patients (39%). Pathogenic variants were found in 49 genes and included missense variants (49%), large copy number changes (18%), small insertions and deletions (18%), nonsense variants (8%), splice-site alterations (6%), and promoter variants (<1%). The diagnostic rate varied considerably based on phenotype and was highest for patients with a positive family history of hearing loss or when the loss was congenital and symmetric. The spectrum of implicated genes showed wide ethnic variability. These findings support the more efficient utilization of medical resources through the development of evidence-based algorithms for the diagnosis of hearing loss.

Cordova: web-based management of genetic variation data.

UNLABELLED: Cordova is an out-of-the-box solution for building and maintaining an online database of genetic variations integrated with pathogenicity prediction results from popular algorithms. Our primary motivation for developing this system is to aid researchers and clinician-scientists in determining the clinical significance of genetic variations. To achieve this goal, Cordova provides an interface to review and manually or computationally curate genetic variation data as well as share it for clinical diagnostics and the advancement of research. AVAILABILITY AND IMPLEMENTATION: Cordova is open source under the MIT license and is freely available for download at https://github.com/clcg/cordova.

Hearing loss caused by a P2RX2 mutation identified in a MELAS family with a coexisting mitochondrial 3243AG mutation.

OBJECTIVES: We present a family with a mitochondrial DNA 3243A>G mutation resulting in mitochondrial myopathy, encephalopathy, lactic acidosis, and stroke-like episodes (MELAS), of which some members have hearing loss in which a novel mutation in the P2RX2 gene was identified. METHODS: One hundred ninety-four (194) Japanese subjects from unrelated families were enrolled in the study. Targeted genomic enrichment and massively parallel sequencing of all known nonsyndromic hearing loss genes were performed to identify the genetic causes of hearing loss. RESULTS: A novel mutation in the P2RX2 gene that corresponded to c.601G>A (p.Asp201Tyr) was identified. Two patients carried the mutation and had severe sensorineural hearing loss, while other members with MELAS (who did not carry the P2RX2 mutation) had normal hearing. CONCLUSION: This is the first case report of a diagnosis of hearing loss caused by P2RX2 mutation in patients with MELAS. A potential explanation is that a decrease in adenosine triphosphate (ATP) production due to MELAS with a mitochondrial 3243A>G mutation might suppress activation of P2X2 receptors. We also suggest that hearing loss caused by the P2RX2 mutation might be influenced by the decrease in ATP production due to MELAS.

Novel PTPRQ mutations identified in three congenital hearing loss patients with various types of hearing loss.

OBJECTIVES: We present 3 patients with congenital sensorineural hearing loss (SNHL) caused by novel PTPRQ mutations, including clinical manifestations and phenotypic features. METHODS: Two hundred twenty (220) Japanese subjects with SNHL from unrelated and nonconsanguineous families were enrolled in the study. Targeted genomic enrichment with massively parallel DNA sequencing of all known nonsyndromic hearing loss genes was performed to identify the genetic cause of hearing loss. RESULTS: Four novel causative PTPRQ mutations were identified in 3 cases. Case 1 had progressive profound SNHL with a homozygous nonsense mutation. Case 2 had nonprogressive profound SNHL with a compound heterozygous mutation (nonsense and missense mutation). Case 3 had nonprogressive moderate SNHL with a compound heterozygous mutation (missense and splice site mutation). Caloric test and vestibular evoked myogenic potential (VEMP) test showed vestibular dysfunction in Case 1. CONCLUSION: Hearing loss levels and progression among the present cases were varied, and there seem to be no obvious correlations between genotypes and the phenotypic features of their hearing loss. The PTPRQ mutations appeared to be responsible for vestibular dysfunction.

De novo mutation in X-linked hearing loss-associated POU3F4 in a sporadic case of congenital hearing loss.

OBJECTIVES: In this report, we present a male patient with no family history of hearing loss, in whom we identified a novel de novo mutation in the POU3F4 gene. METHODS: One hundred ninety-four (194) Japanese subjects from unrelated and nonconsanguineous families were enrolled in this study. We used targeted genomic enrichment and massively parallel sequencing of all known nonsyndromic hearing loss genes for identifying the genetic causes of hearing loss. RESULTS: A novel de novo frameshift mutation of POU3F4 to c.727_728insA (p.N244KfsX26) was identified. The patient was a 7-year-old male with congenital progressive hearing loss and inner ear deformity. Although the patient had received a cochlear implant, auditory skills were still limited. The patient also exhibited developmental delays similar to those previously associated with POU3F4 mutation. CONCLUSION: This is the first report of a mutation in POU3F4 causing hearing loss in a Japanese patient without a family history of hearing loss. This study underscores the importance of comprehensive genetic testing of patients with hearing loss for providing accurate prognostic information and guiding the optimal management of patient rehabilitation.

USH2 caused by GPR98 mutation diagnosed by massively parallel sequencing in advance of the occurrence of visual symptoms.

OBJECTIVE: We present 2 patients who were identified with mutations in the GPR98 gene that causes Usher syndrome type 2 (USH2). METHODS: One hundred ninety-four (194) Japanese subjects from unrelated families were enrolled in the study. Targeted genomic enrichment and massively parallel sequencing of all known nonsyndromic hearing loss genes were used to identify the genetic causes of hearing loss. RESULTS: We identified causative mutations in the GPR98 gene in 1 family (2 siblings). The patients had moderate sloping hearing loss, and no progression was observed over a period of 10 years. Fundus examinations were normal. However, electroretinograms revealed impaired responses in both patients. CONCLUSION: Early diagnosis of Usher syndrome has many advantages for patients and their families. This study supports the use of comprehensive genetic diagnosis for Usher syndrome, especially prior to the onset of visual symptoms, to provide the highest chance of diagnostic success in early life stages.

Rfam: Wikipedia, clans and the "decimal" release.

The Rfam database aims to catalogue non-coding RNAs through the use of sequence alignments and statistical profile models known as covariance models. In this contribution, we discuss the pros and cons of using the online encyclopedia, Wikipedia, as a source of community-derived annotation. We discuss the addition of groupings of related RNA families into clans and new developments to the website. Rfam is available on the Web at http://rfam.sanger.ac.uk.

Assessing Variants of Uncertain Significance Implicated in Hearing Loss Using a Comprehensive Deafness Proteome.

Hearing loss is the leading sensory deficit, affecting ~ 5% of the population. It exhibits remarkable heterogeneity across 223 genes with 6,328 pathogenic missense variants, making deafness-specific expertise a prerequisite for ascribing phenotypic consequences to genetic variants. Deafness-implicated variants are curated in the Deafness Variation Database (DVD) after classification by a genetic hearing loss expert panel and thorough informatics pipeline. However, seventy percent of the 128,167 missense variants in the DVD are "variants of uncertain significance" (VUS) due to insufficient evidence for classification. Here, we use the deep learning protein prediction algorithm, AlphaFold2, to curate structures for all DVD genes. We refine these structures with global optimization and the AMOEBA force field and use DDGun3D to predict folding free energy differences (∆∆G Fold ) for all DVD missense variants. We find that 5,772 VUSs have a large, destabilizing ∆∆G Fold that is consistent with pathogenic variants. When also filtered for CADD scores (> 25.7), we determine 3,456 VUSs are likely pathogenic at a probability of 99.0%. These VUSs affect 119 patients (~ 3% of cases) sequenced by the OtoSCOPE targeted panel. Approximately half of these patients previously received an inconclusive report, and reclassification of these VUSs as pathogenic provides a new genetic diagnosis for six patients.

Fast filtering for RNA homology search.

MOTIVATION: Homology search for RNAs can use secondary structure information to increase power by modeling base pairs, as in covariance models, but the resulting computational costs are high. Typical acceleration strategies rely on at least one filtering stage using sequence-only search. RESULTS: Here we present the multi-segment CYK (MSCYK) filter, which implements a heuristic of ungapped structural alignment for RNA homology search. Compared to gapped alignment, this approximation has lower computation time requirements (O(N4) reduced to O(N³), and space requirements (O(N³) reduced to O(N²). A vector-parallel implementation of this method gives up to 100-fold speed-up; vector-parallel implementations of standard gapped alignment at two levels of precision give 3- and 6-fold speed-ups. These approaches are combined to create a filtering pipeline that scores RNA secondary structure at all stages, with results that are synergistic with existing methods.

Rfam: updates to the RNA families database.

Rfam is a collection of RNA sequence families, represented by multiple sequence alignments and covariance models (CMs). The primary aim of Rfam is to annotate new members of known RNA families on nucleotide sequences, particularly complete genomes, using sensitive BLAST filters in combination with CMs. A minority of families with a very broad taxonomic range (e.g. tRNA and rRNA) provide the majority of the sequence annotations, whilst the majority of Rfam families (e.g. snoRNAs and miRNAs) have a limited taxonomic range and provide a limited number of annotations. Recent improvements to the website, methodologies and data used by Rfam are discussed. Rfam is freely available on the Web at http://rfam.sanger.ac.uk/and http://rfam.janelia.org/.

A synonymous variant in MYO15A enriched in the Ashkenazi Jewish population causes autosomal recessive hearing loss due to abnormal splicing.

Nonsyndromic hearing loss is genetically heterogeneous. Despite comprehensive genetic testing, many cases remain unsolved because the clinical significance of identified variants is uncertain or because biallelic pathogenic variants are not identified for presumed autosomal recessive cases. Common synonymous variants are often disregarded. Determining the pathogenicity of synonymous variants may improve genetic diagnosis. We report a synonymous variant c.9861 C > T/p.(Gly3287=) in MYO15A in homozygosity or compound heterozygosity with another pathogenic or likely pathogenic MYO15A variant in 10 unrelated families with nonsyndromic sensorineural hearing loss. Biallelic variants in MYO15A were identified in 21 affected and were absent in 22 unaffected siblings. A mini-gene assay confirms that the synonymous variant leads to abnormal splicing. The variant is enriched in the Ashkenazi Jewish population. Individuals carrying biallelic variants involving c.9861 C > T often exhibit progressive post-lingual hearing loss distinct from the congenital profound deafness typically associated with biallelic loss-of-function MYO15A variants. This study establishes the pathogenicity of the c.9861 C > T variant in MYO15A and expands the phenotypic spectrum of MYO15A-related hearing loss. Our work also highlights the importance of multicenter collaboration and data sharing to establish the pathogenicity of a relatively common synonymous variant for improved diagnosis and management of hearing loss.

Local RNA structure alignment with incomplete sequence.

MOTIVATION: Accuracy of automated structural RNA alignment is improved by using models that consider not only primary sequence but also secondary structure information. However, current RNA structural alignment approaches tend to perform poorly on incomplete sequence fragments, such as single reads from metagenomic environmental surveys, because nucleotides that are expected to be base paired are missing. RESULTS: We present a local RNA structural alignment algorithm, trCYK, for aligning and scoring incomplete sequences under a model using primary sequence conservation and secondary structure information when possible. The trCYK algorithm improves alignment accuracy and coverage of sequence fragments of structural RNAs in simulated metagenomic shotgun datasets. AVAILABILITY: The source code for Infernal 1.0, which includes trCYK, is available at http://infernal.janelia.org.

Infernal 1.0: inference of RNA alignments.

SUMMARY: INFERNAL builds consensus RNA secondary structure profiles called covariance models (CMs), and uses them to search nucleic acid sequence databases for homologous RNAs, or to create new sequence- and structure-based multiple sequence alignments. AVAILABILITY: Source code, documentation and benchmark downloadable from http://infernal.janelia.org. INFERNAL is freely licensed under the GNU GPLv3 and should be portable to any POSIX-compliant operating system, including Linux and Mac OS/X.

Is it Usher syndrome? Collaborative diagnosis and molecular genetics of patients with visual impairment and hearing loss.

Background: Usher syndrome is the most common hereditary syndrome combining deafness and blindness. In the 2017 National Child Count of Children and Youth who are Deaf-Blind, Usher syndrome represented 329 of 10,000 children, but there were also at least 70 other etiologies of deaf-blindness documented. The purpose of this study was to analyze the work-up and ultimate diagnoses of 21 consecutive families who presented to the Genetic Eye-Ear Clinic (GEEC) at the University of Iowa. Our hypothesis was that most families referred to the GEEC would have initial and final diagnoses of Usher syndrome.Materials and Methods: Patients were identified through an IRB approved retrospective chart review of referrals to the GEEC between 2012 and 2019. Details about each patient's history, exam, and clinical and genetic work-up were recorded.Results: From 2012 to 2019, 21 families (25 patients) were referred to the collaborative GEEC. Overall molecular diagnostic rate in this cohort was 14/21 (67%). Evaluation resulted in a change of diagnosis in 11/21 (52%) families. Ultimately, there were eleven unique diagnoses including hereditary, non-hereditary, and independent causes of combined visual impairment and hearing loss. The most common diagnosis was Usher syndrome, which represented 6/21 (29%) families.Conclusions: Providing a correct diagnosis for patients with visual impairment and hearing loss can be challenging for clinicians and their patients, but it can greatly improve clinical care and outcomes. We recommend an algorithm that includes multidisciplinary collaboration, careful clinical evaluation, strategic molecular testing, and consideration of a broad differential diagnosis.

Insights into the Biology of Hearing and Deafness Revealed by Single-Cell RNA Sequencing.

Single-cell RNA sequencing is a powerful tool by which to characterize the transcriptional profile of low-abundance cell types, but its application to the inner ear has been hampered by the bony labyrinth, tissue sparsity, and difficulty dissociating the ultra-rare cells of the membranous cochlea. Herein, we present a method to isolate individual inner hair cells (IHCs), outer hair cells (OHCs), and Deiters' cells (DCs) from the murine cochlea at any post-natal time point. We harvested more than 200 murine IHCs, OHCs, and DCs from post-natal days 15 (p15) to 228 (p228) and leveraged both short- and long-read single-cell RNA sequencing to profile transcript abundance and structure. Our results provide insights into the expression profiles of these cells and document an unappreciated complexity in isoform variety in deafness-associated genes. This refined view of transcription in the organ of Corti improves our understanding of the biology of hearing and deafness.

Comprehensive Genetic Testing for Deafness from Fresh and Archived Dried Blood Spots.

Comprehensive genetic testing has become integral in the evaluation of children with deafness, but the amount of blood required to obtain DNA can be prohibitive in newborns. Dried blood spots (DBSs) are routinely collected and would provide an alternative source of DNA. Our objective was to evaluate the use of DBSs for comprehensive genetic testing for deafness. DNA derived from fresh and archived DBS samples was compared with DNA from whole blood. We performed next-generation sequencing of all known deafness genes in 4 DBS samples: 2 positive controls, an unknown sample, and a negative control. The DBS-derived DNA was of sufficient quantity and quality for clinical testing. In the 2 positive control samples, pathogenic variants were identified; in the negative control, no pathogenic variants were found; and in the unknown sample, homozygous deletion of the OTOA gene was identified as the cause of deafness. This pilot study shows that comprehensive genetic testing for deafness is feasible with fresh and/or archived DBSs.