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

Using multi-scale genomics to associate poorly annotated genes with rare diseases

BACKGROUND

Next-generation sequencing (NGS) has significantly transformed the landscape of identifying disease-causing genes associated with genetic disorders. However, a substantial portion of sequenced patients remains undiagnosed. This may be attributed not only to the challenges posed by harder-to-detect variants, such as non-coding and structural variations but also to the existence of variants in genes not previously associated with the patient's clinical phenotype. This study introduces EvORanker, an algorithm that integrates unbiased data from 1,028 eukaryotic genomes to link mutated genes to clinical phenotypes.

METHODS

EvORanker utilizes clinical data, multi-scale phylogenetic profiling, and other omics data to prioritize disease-associated genes. It was evaluated on solved exomes and simulated genomes, compared with existing methods, and applied to 6260 knockout genes with mouse phenotypes lacking human associations. Additionally, EvORanker was made accessible as a user-friendly web tool.

RESULTS

In the analyzed exomic cohort, EvORanker accurately identified the "true" disease gene as the top candidate in 69% of cases and within the top 5 candidates in 95% of cases, consistent with results from the simulated dataset. Notably, EvORanker outperformed existing methods, particularly for poorly annotated genes. In the case of the 6260 knockout genes with mouse phenotypes, EvORanker linked 41% of these genes to observed human disease phenotypes. Furthermore, in two unsolved cases, EvORanker successfully identified DLGAP2 and LPCAT3 as disease candidates for previously uncharacterized genetic syndromes.

CONCLUSIONS

We highlight clade-based phylogenetic profiling as a powerful systematic approach for prioritizing potential disease genes. Our study showcases the efficacy of EvORanker in associating poorly annotated genes to disease phenotypes observed in patients. The EvORanker server is freely available at https://ccanavati.shinyapps.io/EvORanker/ .

Multi-animal-model study reveals mutations in neural plasticity and nociception genes linked to excessive alcohol drinking

BACKGROUND

The basis for familial alcohol use disorder (AUD) remains an enigma due to various biological and societal confounds. The present study used three of the most adopted and documented rat models, combining the alcohol-preferring/non-alcohol-preferring (P/NP) lines and high alcohol-drinking/low alcohol-drinking (HAD/LAD) replicated lines, of AUD as examined through the lens of whole genomic analyses.

METHODS

We used complete genome sequencing of the P/NP lines and previously published sequences of the HAD/LAD replicates to enhance the discovery of variants associated with AUD and to remove confounding with genetic background and random genetic drift. Specifically, we used high-order statistical methods to search for genetic variants whose frequency changes in whole sets of gene ontologies corresponded with phenotypic changes in the direction of selection, that is, ethanol-drinking preference.

RESULTS

Our first finding was that in addition to variants causing translational changes, the principal genetic changes associated with drinking predisposition were silent mutations and mutations in the 3' untranslated regions (3'UTR) of genes. Neither of these types of mutations alters the amino acid sequence of the translated protein but they influence both the rate and conformation of gene transcription, including its stability and posttranslational events that alter gene efficacy. This finding argues for refocusing human genomic studies on changes in gene efficacy. Among the key ontologies identified were the central genes associated with the Na+ voltage-gated channels of neurons and glia (including the Scn1a, Scn2a, Scn2b, Scn3a, Scn7a, and Scn9a subtypes) and excitatory glutamatergic secretion (including Grm2 and Myo6), both of which are essential in neuroplasticity. In addition, we identified "Nociception or Sensory Perception of Pain," which contained variants in nociception (Arrb1, Ccl3, Ephb1) and enlist sodium (Scn1a, Scn2a, Scn2b, Scn3a, Scn7a), pain activation (Scn9a), and potassium channel (Kcna1) genes.

CONCLUSION

The multi-model analyses used herein reduced the confounding effects of random drift and the "founders" genetic background. The most differentiated bidirectionally selected genes across all three animal models were Scn9a, Scn1a, and Kcna, all of which are annotated in the nociception ontology. The complexity of neuroplasticity and nociception adds strength to the hypothesis that neuroplasticity and pain (physical or psychological) are prominent phenotypes genetically linked to the development of AUD.

Biallelic Loss-of-Function Variants in BICD1 Are Associated with Peripheral Neuropathy and Hearing Loss

Hearing loss and peripheral neuropathy are two clinical entities that are genetically and phenotypically heterogeneous and sometimes co-occurring. Using exome sequencing and targeted segregation analysis, we investigated the genetic etiology of peripheral neuropathy and hearing loss in a large Ashkenazi Jewish family. Moreover, we assessed the production of the candidate protein via western blotting of lysates from fibroblasts from an affected individual and an unaffected control. Pathogenic variants in known disease genes associated with hearing loss and peripheral neuropathy were excluded. A homozygous frameshift variant in the BICD1 gene, c.1683dup (p.(Arg562Thrfs*18)), was identified in the proband and segregated with hearing loss and peripheral neuropathy in the family. The BIDC1 RNA analysis from patient fibroblasts showed a modest reduction in gene transcripts compared to the controls. In contrast, protein could not be detected in fibroblasts from a homozygous c.1683dup individual, whereas BICD1 was detected in an unaffected individual. Our findings indicate that bi-allelic loss-of-function variants in BICD1 are associated with hearing loss and peripheral neuropathy. Definitive evidence that bi-allelic loss-of-function variants in BICD1 cause peripheral neuropathy and hearing loss will require the identification of other families and individuals with similar variants with the same phenotype.

Implementing computational methods in tandem with synonymous gene recoding for therapeutic development

Synonymous gene recoding, the substitution of synonymous variants into the genetic sequence, has been used to overcome many production limitations in therapeutic development. However, the safety and efficacy of recoded therapeutics can be difficult to evaluate because synonymous codon substitutions can result in subtle, yet impactful changes in protein features and require sensitive methods for detection. Given that computational approaches have made significant leaps in recent years, we propose that machine-learning (ML) tools may be leveraged to assess gene-recoded therapeutics and foresee an opportunity to adapt codon contexts to enhance some powerful existing tools. Here, we examine how synonymous gene recoding has been used to address challenges in therapeutic development, explain the biological mechanisms underlying its effects, and explore the application of computational platforms to improve the surveillance of functional variants in therapeutic design.

Novel compound heterozygous synonymous and missense variants in the MYO7A gene identified by next-generation sequencing in a Chinese family with nonsyndromic hearing loss

BACKGROUND

Variants in the MYO7A gene are increasingly identified among patients suffering from Usher syndrome type 1B (USH1B). However, such mutations are less commonly detected among patients suffering from nonsyndromic hearing loss (NSHL), including autosomal recessive deafness (DFNB2) and autosomal dominant deafness (DFNA11). This research attempts to clarify the genetic base of DFNB2 in a Chinese family and determine the pathogenicity of the identified mutations.

METHOD

Targeted next-generation sequencing (TGS) of 127 known deafness genes was performed for the 14-year-old proband. Then, Sanger sequencing was performed on the available family members. A minigene splicing assay was performed to verify the impact of the novel MYO7A synonymous variant. After performing targeted next-generation sequencing (TGS) of 127 existing hearing loss-related genes in a 14-year-old proband, Sanger sequencing was carried out on the available family members. Then, to confirm the influence of the novel MYO7A synonymous variants, a minigene splicing assay was performed.

RESULTS

Two heteroallelic mutants of MYO7A (NM_000260.3) were identified: a maternally inherited synonymous variant c.2904G > A (p.Glu968=) in exon 23 and a paternally inherited missense variant c.5994G > T (p.Trp1998Cys) in exon 44. The in vitro minigene expression indicated that c.2904G > A may result in skipping of exon 23 resulting in a truncated protein.

CONCLUSIONS

We reported a novel missense (c.5994G > T) and identified, for the first time, a novel pathogenic synonymous (c.2904G > A) variant within MYO7A in a patient with DFNB2. These findings enrich our understanding of the MYO7A variant spectrum of DFNB2 and can contribute to accurate genetic counseling and diagnosis of NSHL patients.

Genetic etiology of hearing loss in Iran

Hearing loss (HL) is an etiologically heterogeneous disorder that affects around 5% of the world's population. There has been an exponential increase in the identification of genes and variants responsible for hereditary HL over recent years. Iran, a country located in the Middle East, has a high prevalence of consanguineous marriages, so heterogeneous diseases such as HL are more common. Comprehensive studies using different strategies from linkage analysis to next-generation sequencing, especially exome-sequencing, have achieved significant success in identifying possible pathogens in deaf Iranian families. About 12% of non-syndromic autosomal recessive HL genes investigated to date, were first identified in families from Iran. Variations of 56 genes have been observed in families with NSHL in Iran. Variants in GJB2, SLC26A4, MYO15A, MYO7A, CDH23, and TMC1 account for 16.5%, 16.25%, 13.5%, 9.35%, 6.9% and 4.92%, cases of NSHL, respectively. In summary, there are also different diagnostic rates between studies conducted in Iran. In the comprehensive investigations conducted by the Genetic Research Center of the University of Social Welfare and Rehabilitation Sciences over the past 20 years, the overall diagnosis rate is about 80% while there are other studies with lower diagnostic rates which could reflect differences in project designs, sampling, and accuracy and validity of the methods used. Furthermore, there are several syndromic HHLs in Iran including, Waardenburg syndrome, BOR syndrome, Brown-Vialetto-Van Laere syndrome, Wolfram syndrome, among which Pendred and Usher syndromes are well-studied. These results are of importance for further investigation and elucidation of the molecular basis of HHL in Iran.

Comprehensive molecular-genetic analysis of mid-frequency sensorineural hearing loss

The genetic heterogeneity of sensorineural hearing loss (SNHL) is a major hurdle to the detection of disease-causing variants. We aimed to identify underlying causal genes associated with mid-frequency hearing loss (HL), which contributes to less than about 1% of SNHL cases, by whole exome sequencing (WES). Thirty families segregating mid-frequency SNHL, in whom biallelic GJB2 mutations had been previously excluded, were selected from among 851 families in our DNA repository of SNHL. DNA samples from the probands were subjected to WES analysis and searched for candidate variants associated with SNHL. We were able to identify the genetic aetiology in six probands (20%). In total, we found three pathogenic and three likely pathogenic variants in four genes (COL4A5, OTOGL, TECTA, TMPRSS3). One more proband was a compound heterozygote for a pathogenic variant and a variant of uncertain significance (VUS) in MYO15A gene. To date, MYO15A and TMPRSS3 have not yet been described in association with mid-frequency SNHL. In eight additional probands, eight candidate VUS variants were detected in five genes (DIAPH1, MYO7A, TECTA, TMC1, TSPEAR). Seven of these 16 variants have not yet been published or mentioned in the available databases. The most prevalent gene was TECTA, identified in 23% of all tested families. Furthermore, we confirmed the hypothesis that a substantive portion of cases with this conspicuous audiogram shape is a consequence of a genetic disorder.

Identification of Novel and Recurrent Variants in MYO15A in Ashkenazi Jewish Patients With Autosomal Recessive Nonsyndromic Hearing Loss

Hearing loss is a genetically and phenotypically heterogeneous disorder. The purpose of this study was to determine the genetic cause underlying hearing loss in four Ashkenazi Jewish families. We screened probands from each family using a combination of targeted mutation screening and exome sequencing to identifiy the genetic cause of hearing loss in each family. We identified four variants in MYO15A, two novel variants never previously linked to deafness (c.7212+5G>A and p.Leu2532ArgfsTer37) and two recurrent variants (p.Tyr2684His and p.Gly3287Gly). One family showed locus heterogeneity, segregrating two genetic forms of hearing loss. Mini-gene assays revealed the c.7212+5G>A variant results in abnormal splicing and is most likely a null allele. We show that families segregrating the p.Gly3287Gly variant show both inter and intra-familial phenotypic differences. These results add to the list of MYO15A deafness-causing variants, further confirm the pathogenicity of the p.Gly3287Gly variant and shed further light on the genetic etiology of hearing loss in the Ashkenazi Jewish population.