MRSD: A quantitative approach for assessing suitability of RNA-seq in the investigation of mis-splicing in Mendelian disease.

Variable levels of gene expression between tissues complicates the use of RNA sequencing of patient biosamples to delineate the impact of genomic variants. Here, we describe a gene- and tissue-specific metric to inform the feasibility of RNA sequencing. This overcomes limitations of using expression values alone as a metric to predict RNA-sequencing utility. We have derived a metric, minimum required sequencing depth (MRSD), that estimates the depth of sequencing required from RNA sequencing to achieve user-specified sequencing coverage of a gene, transcript, or group of genes. We applied MRSD across four human biosamples: whole blood, lymphoblastoid cell lines (LCLs), skeletal muscle, and cultured fibroblasts. MRSD has high precision (90.1%-98.2%) and overcomes transcript region-specific sequencing biases. Applying MRSD scoring to established disease gene panels shows that fibroblasts, of these four biosamples, are the optimum source of RNA for 63.1% of gene panels. Using this approach, up to 67.8% of the variants of uncertain significance in ClinVar that are predicted to impact splicing could be assayed by RNA sequencing in at least one of the biosamples. We demonstrate the utility and benefits of MRSD as a metric to inform functional assessment of splicing aberrations, in particular in the context of Mendelian genetic disorders to improve diagnostic yield.

Global mapping of RNA homodimers in living cells.

RNA homodimerization is important for various physiological processes, including the assembly of membraneless organelles, RNA subcellular localization, and packaging of viral genomes. However, understanding RNA dimerization has been hampered by the lack of systematic in vivo detection methods. Here, we show that CLASH, PARIS, and other RNA proximity ligation methods detect RNA homodimers transcriptome-wide as "overlapping" chimeric reads that contain more than one copy of the same sequence. Analyzing published proximity ligation data sets, we show that RNA:RNA homodimers mediated by direct base-pairing are rare across the human transcriptome, but highly enriched in specific transcripts, including U8 snoRNA, U2 snRNA, and a subset of tRNAs. Mutations in the homodimerization domain of U8 snoRNA impede dimerization in vitro and disrupt zebrafish development in vivo, suggesting an evolutionarily conserved role of this domain. Analysis of virus-infected cells reveals homodimerization of SARS-CoV-2 and Zika genomes, mediated by specific palindromic sequences located within protein-coding regions of N gene in SARS-CoV-2 and NS2A gene in Zika. We speculate that regions of viral genomes involved in homodimerization may constitute effective targets for antiviral therapies.

Bi-allelic variants in the mitochondrial RNase P subunit PRORP cause mitochondrial tRNA processing defects and pleiotropic multisystem presentations.

Human mitochondrial RNase P (mt-RNase P) is responsible for 5' end processing of mitochondrial precursor tRNAs, a vital step in mitochondrial RNA maturation, and is comprised of three protein subunits: TRMT10C, SDR5C1 (HSD10), and PRORP. Pathogenic variants in TRMT10C and SDR5C1 are associated with distinct recessive or x-linked infantile onset disorders, resulting from defects in mitochondrial RNA processing. We report four unrelated families with multisystem disease associated with bi-allelic variants in PRORP, the metallonuclease subunit of mt-RNase P. Affected individuals presented with variable phenotypes comprising sensorineural hearing loss, primary ovarian insufficiency, developmental delay, and brain white matter changes. Fibroblasts from affected individuals in two families demonstrated decreased steady state levels of PRORP, an accumulation of unprocessed mitochondrial transcripts, and decreased steady state levels of mitochondrial-encoded proteins, which were rescued by introduction of the wild-type PRORP cDNA. In mt-tRNA processing assays performed with recombinant mt-RNase P proteins, the disease-associated variants resulted in diminished mitochondrial tRNA processing. Identification of disease-causing variants in PRORP indicates that pathogenic variants in all three subunits of mt-RNase P can cause mitochondrial dysfunction, each with distinct pleiotropic clinical presentations.

Functional and transcriptional profiling of non-coding RNAs in yeast reveal context-dependent phenotypes and in trans effects on the protein regulatory network.

Non-coding RNAs (ncRNAs), including the more recently identified Stable Unannotated Transcripts (SUTs) and Cryptic Unstable Transcripts (CUTs), are increasingly being shown to play pivotal roles in the transcriptional and post-transcriptional regulation of genes in eukaryotes. Here, we carried out a large-scale screening of ncRNAs in Saccharomyces cerevisiae, and provide evidence for SUT and CUT function. Phenotypic data on 372 ncRNA deletion strains in 23 different growth conditions were collected, identifying ncRNAs responsible for significant cellular fitness changes. Transcriptome profiles were assembled for 18 haploid ncRNA deletion mutants and 2 essential ncRNA heterozygous deletants. Guided by the resulting RNA-seq data we analysed the genome-wide dysregulation of protein coding genes and non-coding transcripts. Novel functional ncRNAs, SUT125, SUT126, SUT035 and SUT532 that act in trans by modulating transcription factors were identified. Furthermore, we described the impact of SUTs and CUTs in modulating coding gene expression in response to different environmental conditions, regulating important biological process such as respiration (SUT125, SUT126, SUT035, SUT432), steroid biosynthesis (CUT494, SUT053, SUT468) or rRNA processing (SUT075 and snR30). Overall, these data capture and integrate the regulatory and phenotypic network of ncRNAs and protein-coding genes, providing genome-wide evidence of the impact of ncRNAs on cellular homeostasis.

Analysis of U8 snoRNA Variants in Zebrafish Reveals How Bi-allelic Variants Cause Leukoencephalopathy with Calcifications and Cysts.

How mutations in the non-coding U8 snoRNA cause the neurological disorder leukoencephalopathy with calcifications and cysts (LCC) is poorly understood. Here, we report the generation of a mutant U8 animal model for interrogating LCC-associated pathology. Mutant U8 zebrafish exhibit defective central nervous system development, a disturbance of ribosomal RNA (rRNA) biogenesis and tp53 activation, which monitors ribosome biogenesis. Further, we demonstrate that fibroblasts from individuals with LCC are defective in rRNA processing. Human precursor-U8 (pre-U8) containing a 3' extension rescued mutant U8 zebrafish, and this result indicates conserved biological function. Analysis of LCC-associated U8 mutations in zebrafish revealed that one null and one functional allele contribute to LCC. We show that mutations in three nucleotides at the 5' end of pre-U8 alter the processing of the 3' extension, and we identify a previously unknown base-pairing interaction between the 5' end and the 3' extension of human pre-U8. Indeed, LCC-associated mutations in any one of seven nucleotides in the 5' end and 3' extension alter the processing of pre-U8, and these mutations are present on a single allele in almost all individuals with LCC identified to date. Given genetic data indicating that bi-allelic null U8 alleles are likely incompatible with human development, and that LCC is not caused by haploinsufficiency, the identification of hypomorphic misprocessing mutations that mediate viable embryogenesis furthers our understanding of LCC molecular pathology and cerebral vascular homeostasis.

New insights into Perrault syndrome, a clinically and genetically heterogeneous disorder.

Hearing loss and impaired fertility are common human disorders each with multiple genetic causes. Sometimes deafness and impaired fertility, which are the hallmarks of Perrault syndrome, co-occur in a person. Perrault syndrome is inherited as an autosomal recessive disorder characterized by bilateral mild to severe childhood sensorineural hearing loss with variable age of onset in both sexes and ovarian dysfunction in females who have a 46, XX karyotype. Since the initial clinical description of Perrault syndrome 70 years ago, the phenotype of some subjects may additionally involve developmental delay, intellectual deficit and other neurological disabilities, which can vary in severity in part dependent upon the genetic variants and the gene involved. Here, we review the molecular genetics and clinical phenotype of Perrault syndrome and focus on supporting evidence for the eight genes (CLPP, ERAL1, GGPS1, HARS2, HSD17B4, LARS2, RMND1, TWNK) associated with Perrault syndrome. Variants of these eight genes only account for approximately half of the individuals with clinical features of Perrault syndrome where the molecular genetic base remains under investigation. Additional environmental etiologies and novel Perrault disease-associated genes remain to be identified to account for unresolved cases. We also report a new genetic variant of CLPP, computational structural insight about CLPP and single cell RNAseq data for eight reported Perrault syndrome genes suggesting a common cellular pathophysiology for this disorder. Some unanswered questions are raised to kindle future research about Perrault syndrome.

Expanding the genotypic spectrum of TXNL4A variants in Burn-McKeown syndrome.

The developmental disorder Burn-McKeown Syndrome (BMKS) is characterised by choanal atresia and specific craniofacial features. BMKS is caused by biallelic variants in the pre-messenger RNA splicing factor TXNL4A. Most patients have a loss-of-function variant in trans with a 34-base pair (bp) deletion (type 1 Δ34) in the promoter region. Here, we identified two patients with BMKS. One individual has a TXNL4A c.93_94delCC, p.His32Argfs *21 variant combined with a type 1 Δ34 promoter deletion. The other has an intronic TXNL4A splice site variant (c.258-3C>G) and a type 1 Δ34 promoter deletion. We show the c.258-3C>G variant and a previously reported c.258-2A>G variant, cause skipping of the final exon of TXNL4A in a minigene splicing assay. Furthermore, we identify putative transcription factor binding sites within the 56 bp of the TXNL4A promoter affected by the type 1 and type 2 Δ34 and use dual luciferase assays to identify a 22 bp repeated motif essential for TXNL4A expression within this promoter region. We propose that additional variants affecting critical transcription factor binding nucleotides within the 22 bp repeated motif could be relevant to BMKS aetiology. Finally, our data emphasises the need to analyse the non-coding sequence in individuals where a single likely pathogenic coding variant is identified in an autosomal recessive disorder consistent with the clinical presentation.

Homozygous missense variants in BMPR15 can result in primary ovarian insufficiency.

RESEARCH QUESTION: Does a genetic condition underlie the diagnosis of primary ovarian insufficiency (POI) in a 13-year-old girl with primary amenorrhoea? DESIGN: A case report of a next-generation sequencing panel of 24 genes associated with syndromal and non-syndromal POI was conducted. RESULTS: A homozygous missense variant c.1076C>T, p.(Pro359Leu) in BMP15 was identified. CONCLUSIONS: The biallelic variant c.1076C >T, p.(Pro359Leu) in BMP15 is associated with primary ovarian failure.

Disease modeling of core pre-mRNA splicing factor haploinsufficiency.

The craniofacial disorder mandibulofacial dysostosis Guion-Almeida type is caused by haploinsufficiency of the U5 snRNP gene EFTUD2/SNU114. However, it is unclear how reduced expression of this core pre-mRNA splicing factor leads to craniofacial defects. Here we use a CRISPR-Cas9 nickase strategy to generate a human EFTUD2-knockdown cell line and show that reduced expression of EFTUD2 leads to diminished proliferative ability of these cells, increased sensitivity to endoplasmic reticulum (ER) stress and the mis-expression of several genes involved in the ER stress response. RNA-Seq analysis of the EFTUD2-knockdown cell line revealed transcriptome-wide changes in gene expression, with an enrichment for genes associated with processes involved in craniofacial development. Additionally, our RNA-Seq data identified widespread mis-splicing in EFTUD2-knockdown cells. Analysis of the functional and physical characteristics of mis-spliced pre-mRNAs highlighted conserved properties, including length and splice site strengths, of retained introns and skipped exons in our disease model. We also identified enriched processes associated with the affected genes, including cell death, cell and organ morphology and embryonic development. Together, these data support a model in which EFTUD2 haploinsufficiency leads to the mis-splicing of a distinct subset of pre-mRNAs with a widespread effect on gene expression, including altering the expression of ER stress response genes and genes involved in the development of the craniofacial region. The increased burden of unfolded proteins in the ER resulting from mis-splicing would exceed the capacity of the defective ER stress response, inducing apoptosis in cranial neural crest cells that would result in craniofacial abnormalities during development.

The role of splicing factors in retinitis pigmentosa: links to cilia.

Cilia are critical to numerous biological functions, both in development and everyday homeostatic processes. Diseases arising from genetic mutations that cause cilia dysfunction are termed ciliopathies. Several ubiquitously expressed splicing factors have been implicated in the condition Retinitis Pigmentosa (RP), a group of diseases characterised by the progressive degeneration of the retina. In many types of RP the disease affects the modified primary cilium of the photoreceptor cells and thus, these types of RP are considered ciliopathies. Here, we discuss sequence variants found within a number of these splicing factors, the resulting phenotypes, and the mechanisms underpinning disease pathology. Additionally, we discuss recent evidence investigating why RP patients with mutations in globally expressed splicing factors present with retina-specific phenotypes.

Biallelic loss of function variants in STAG3 result in primary ovarian insufficiency.

RESEARCH QUESTION: Does a genetic condition underlie the diagnosis of primary ovarian insufficiency (POI) in a 21-year-old woman with primary amenorrhoea? DESIGN: A karyotype and genetic testing for Fragile X syndrome was undertaken. A next-generation sequencing panel of 24 genes associated with syndromal and non-syndromal POI was conducted. RESULTS: A nonsense variant c.1336G>T, p.(Glu446Ter) and whole gene deletion in STAG3 were identified. CONCLUSIONS: Biallelic loss of function variants in STAG3 are associated with primary ovarian failure type 8 and are a rare cause of POI.

Large-scale profiling of noncoding RNA function in yeast.

Noncoding RNAs (ncRNAs) are emerging as key regulators of cellular function. We have exploited the recently developed barcoded ncRNA gene deletion strain collections in the yeast Saccharomyces cerevisiae to investigate the numerous ncRNAs in yeast with no known function. The ncRNA deletion collection contains deletions of tRNAs, snoRNAs, snRNAs, stable unannotated transcripts (SUTs), cryptic unstable transcripts (CUTs) and other annotated ncRNAs encompassing 532 different individual ncRNA deletions. We have profiled the fitness of the diploid heterozygous ncRNA deletion strain collection in six conditions using batch and continuous liquid culture, as well as the haploid ncRNA deletion strain collections arrayed individually onto solid rich media. These analyses revealed many novel environmental-specific haplo-insufficient and haplo-proficient phenotypes providing key information on the importance of each specific ncRNA in every condition. Co-fitness analysis using fitness data from the heterozygous ncRNA deletion strain collection identified two ncRNA groups required for growth during heat stress and nutrient deprivation. The extensive fitness data for each ncRNA deletion strain has been compiled into an easy to navigate database called Yeast ncRNA Analysis (YNCA). By expanding the original ncRNA deletion strain collection we identified four novel essential ncRNAs; SUT527, SUT075, SUT367 and SUT259/691. We defined the effects of each new essential ncRNA on adjacent gene expression in the heterozygote background identifying both repression and induction of nearby genes. Additionally, we discovered a function for SUT527 in the expression, 3' end formation and localization of SEC4, an essential protein coding mRNA. Finally, using plasmid complementation we rescued the SUT075 lethal phenotype revealing that this ncRNA acts in trans. Overall, our findings provide important new insights into the function of ncRNAs.

EFTUD2 missense variants disrupt protein function and splicing in mandibulofacial dysostosis Guion-Almeida type.

Pathogenic variants in the core spliceosome U5 small nuclear ribonucleoprotein gene EFTUD2/SNU114 cause the craniofacial disorder mandibulofacial dysostosis Guion-Almeida type (MFDGA). MFDGA-associated variants in EFTUD2 comprise large deletions encompassing EFTUD2, intragenic deletions and single nucleotide truncating or missense variants. These variants are predicted to result in haploinsufficiency by loss-of-function of the variant allele. While the contribution of deletions within EFTUD2 to allele loss-of-function are self-evident, the mechanisms by which missense variants are disease-causing have not been characterized functionally. Combining bioinformatics software prediction, yeast functional growth assays, and a minigene (MG) splicing assay, we have characterized how MFDGA missense variants result in EFTUD2 loss-of-function. Only four of 19 assessed missense variants cause EFTUD2 loss-of-function through altered protein function when modeled in yeast. Of the remaining 15 missense variants, five altered the normal splicing pattern of EFTUD2 pre-messenger RNA predominantly through exon skipping or cryptic splice site activation, leading to the introduction of a premature termination codon. Comparison of bioinformatic predictors for each missense variant revealed a disparity amongst different software packages and, in many cases, an inability to correctly predict changes in splicing subsequently determined by MG interrogation. This study highlights the need for laboratory-based validation of bioinformatic predictions for EFTUD2 missense variants.

A recurrent missense variant in HARS2 results in variable sensorineural hearing loss in three unrelated families.

HARS2 encodes mitochondrial histidyl-tRNA synthetase (HARS2), which links histidine to its cognate tRNA in the mitochondrial matrix. Biallelic variants in HARS2 are associated with Perrault syndrome, a rare recessive condition characterized by sensorineural hearing loss in both sexes and primary ovarian insufficiency in 46,XX females. Some individuals with Perrault syndrome have a broader phenotypic spectrum with neurological features, including ataxia and peripheral neuropathy. Here, we report a recurrent variant in HARS2 in association with sensorineural hearing loss. In affected individuals from three unrelated families, the variant HARS2 c.1439G>A p.(Arg480His) is present as a heterozygous variant in trans to a putative pathogenic variant. The low prevalence of the allele HARS2 c.1439G>A p.(Arg480His) in the general population and its presence in three families with hearing loss, confirm the pathogenicity of this variant and illustrate the presentation of Perrault syndrome as nonsyndromic hearing loss in males and prepubertal females.

A resource for functional profiling of noncoding RNA in the yeast Saccharomyces cerevisiae.

Eukaryotic genomes are extensively transcribed, generating many different RNAs with no known function. We have constructed 1502 molecular barcoded ncRNA gene deletion strains encompassing 443 ncRNAs in the yeast Saccharomyces cerevisiae as tools for ncRNA functional analysis. This resource includes deletions of small nuclear RNAs (snRNAs), transfer RNAs (tRNAs), small nucleolar RNAs (snoRNAs), and other annotated ncRNAs as well as the more recently identified stable unannotated transcripts (SUTs) and cryptic unstable transcripts (CUTs) whose functions are largely unknown. Specifically, deletions have been constructed for ncRNAs found in the intergenic regions, not overlapping genes or their promoters (i.e., at least 200 bp minimum distance from the closest gene start codon). The deletion strains carry molecular barcodes designed to be complementary with the protein gene deletion collection enabling parallel analysis experiments. These strains will be useful for the numerous genomic and molecular techniques that utilize deletion strains, including genome-wide phenotypic screens under different growth conditions, pooled chemogenomic screens with drugs or chemicals, synthetic genetic array analysis to uncover novel genetic interactions, and synthetic dosage lethality screens to analyze gene dosage. Overall, we created a valuable resource for the RNA community and for future ncRNA research.

A homozygous missense variant in CHRM3 associated with familial urinary bladder disease.

CHRM3 codes for the M3 muscarinic acetylcholine receptor that is located on the surface of smooth muscle cells of the detrusor, the muscle that effects urinary voiding. Previously, we reported brothers in a family affected by a congenital prune belly-like syndrome with mydriasis due to homozygous CHRM3 frameshift variants. In this study, we describe two sisters with bladders that failed to empty completely and pupils that failed to constrict fully in response to light, who are homozygous for the missense CHRM3 variant c.352G > A; p.(Gly118Arg). Samples were not available for genotyping from their brother, who had a history of multiple urinary tract infections and underwent surgical bladder draining in the first year of life. He died at the age of 6 years. This is the first independent report of biallelic variants in CHRM3 in a family with a rare serious bladder disorder associated with mydriasis and provides important evidence of this association.

Compound heterozygosity of low-frequency promoter deletions and rare loss-of-function mutations in TXNL4A causes Burn-McKeown syndrome.

Mutations in components of the major spliceosome have been described in disorders with craniofacial anomalies, e.g., Nager syndrome and mandibulofacial dysostosis type Guion-Almeida. The U5 spliceosomal complex of eight highly conserved proteins is critical for pre-mRNA splicing. We identified biallelic mutations in TXNL4A, a member of this complex, in individuals with Burn-McKeown syndrome (BMKS). This rare condition is characterized by bilateral choanal atresia, hearing loss, cleft lip and/or palate, and other craniofacial dysmorphisms. Mutations were found in 9 of 11 affected families. In 8 families, affected individuals carried a rare loss-of-function mutation (nonsense, frameshift, or microdeletion) on one allele and a low-frequency 34 bp deletion (allele frequency 0.76%) in the core promoter region on the other allele. In a single highly consanguineous family, formerly diagnosed as oculo-oto-facial dysplasia, the four affected individuals were homozygous for a 34 bp promoter deletion, which differed from the promoter deletion in the other families. Reporter gene and in vivo assays showed that the promoter deletions led to reduced expression of TXNL4A. Depletion of TXNL4A (Dib1) in yeast demonstrated reduced assembly of the tri-snRNP complex. Our results indicate that BMKS is an autosomal-recessive condition, which is frequently caused by compound heterozygosity of low-frequency promoter deletions in combination with very rare loss-of-function mutations.

Non-coding RNAs and disease: the classical ncRNAs make a comeback.

Many human diseases have been attributed to mutation in the protein coding regions of the human genome. The protein coding portion of the human genome, however, is very small compared with the non-coding portion of the genome. As such, there are a disproportionate number of diseases attributed to the coding compared with the non-coding portion of the genome. It is now clear that the non-coding portion of the genome produces many functional non-coding RNAs and these RNAs are slowly being linked to human diseases. Here we discuss examples where mutation in classical non-coding RNAs have been attributed to human disease and identify the future potential for the non-coding portion of the genome in disease biology.

Remodeling of U2-U6 snRNA helix I during pre-mRNA splicing by Prp16 and the NineTeen Complex protein Cwc2.

Removal of intron regions from pre-messenger RNA (pre-mRNA) requires spliceosome assembly with pre-mRNA, then subsequent spliceosome remodeling to allow activation for the two steps of intron removal. Spliceosome remodeling is carried out through the action of DExD/H-box ATPases that modulate RNA-RNA and protein-RNA interactions. The ATPase Prp16 remodels the spliceosome between the first and second steps of splicing by catalyzing release of first step factors Yju2 and Cwc25 as well as destabilizing U2-U6 snRNA helix I. How Prp16 destabilizes U2-U6 helix I is not clear. We show that the NineTeen Complex (NTC) protein Cwc2 displays genetic interactions with the U6 ACAGAGA, the U6 internal stem loop (ISL) and the U2-U6 helix I, all RNA elements that form the spliceosome active site. We find that one function of Cwc2 is to stabilize U2-U6 snRNA helix I during splicing. Cwc2 also functionally cooperates with the NTC protein Isy1/NTC30. Mutation in Cwc2 can suppress the cold sensitive phenotype of the prp16-302 mutation indicating a functional link between Cwc2 and Prp16. Specifically the prp16-302 mutation in Prp16 stabilizes Cwc2 interactions with U6 snRNA and destabilizes Cwc2 interactions with pre-mRNA, indicating antagonistic functions of Cwc2 and Prp16. We propose that Cwc2 is a target for Prp16-mediated spliceosome remodeling during pre-mRNA splicing.