Explore the vast range of titles available.

Predicting the impact of coding and noncoding variants on splicing is challenging, particularly in non-canonical splice sites, leading to missed diagnoses in patients. Existing splice prediction tools are complementary but knowing which to use for each splicing context remains difficult. Here, we describe Introme, which uses machine learning to integrate predictions from several splice detection tools, additional splicing rules, and gene architecture features to comprehensively evaluate the likelihood of a variant impacting splicing. Through extensive benchmarking across 21,000 splice-altering variants, Introme outperformed all tools (auPRC: 0.98) for the detection of clinically significant splice variants. Introme is available at https://github.com/CCICB/introme .

Background The presence of nuclear mitochondrial DNA (numtDNA) has been reported within several nuclear genomes. Next to mitochondrial protein-coding genes, numtDNA sequences also encode for mitochondrial tRNA genes. However, the biological roles of numtDNA remain elusive. Results Employing in silico analysis, we identify 281 mitochondrial tRNA homologs in the human genome, which we term nimtRNAs (nuclear intronic mitochondrial-derived tRNAs), being contained within introns of 76 nuclear host genes. Despite base changes in nimtRNAs when compared to their mtRNA homologs, a canonical tRNA cloverleaf structure is maintained. To address potential functions of intronic nimtRNAs, we insert them into introns of constitutive and alternative splicing reporters and demonstrate that nimtRNAs promote pre-mRNA splicing, dependent on the number and positioning of nimtRNA genes and splice site recognition efficiency. A mutational analysis reveals that the nimtRNA cloverleaf structure is required for the observed splicing increase. Utilizing a CRISPR/Cas9 approach, we show that a partial deletion of a single endogenous nimtRNA Lys within intron 28 of the PPFIBP1 gene decreases inclusion of the downstream-located exon 29 of the PPFIBP1 mRNA. By employing a pull-down approach followed by mass spectrometry, a 3′-splice site-associated protein network is identified, including KHDRBS1, which we show directly interacts with nimtRNA Tyr by an electrophoretic mobility shift assay. Conclusions We propose that nimtRNAs, along with associated protein factors, can act as a novel class of intronic splicing regulatory elements in the human genome by participating in the regulation of splicing.

The viral regulatory protein Tat is essential for establishing a productive transcription from the 5'-LTR promoter during the early phase of viral gene expression. Formation of the Tat-encoding mRNAs requires splicing at the viral 3'ss A3, which has previously been shown to be both negatively and positively regulated by the downstream splicing regulatory elements (SREs) ESS2p and ESE2/ESS2. However, using the novel RESCUE-type computational HEXplorer algorithm, we were recently able to identify another splicing enhancer (ESE(5807-5838), henceforth referred to as ESE tat ) located between ESS2p and ESE2/ESS2. Here we show that ESE tat has a great impact on viral tat-mRNA splicing and that it is fundamental for regulated 3'ss A3 usage. Mutational inactivation or locked nucleic acid (LNA)-directed masking of the ESE tat sequence in the context of a replication-competent virus was associated with a failure (i) to activate viral 3'ss A3 and (ii) to accumulate Tat-encoding mRNA species. Consequently, due to insufficient amounts of Tat protein efficient viral replication was drastically impaired. RNA in vitro binding assays revealed SRSF2 and SRSF6 as candidate splicing factors acting through ESE tat and ESE2 for 3'ss A3 activation. This notion was supported by coexpression experiments, in which wild-type, but not ESE tat -negative provirus responded to higher levels of SRSF2 and SRSF6 proteins with higher levels of tat-mRNA splicing. Remarkably, we could also find that SRSF6 overexpression established an antiviral state within provirus-transfected cells, efficiently blocking virus particle production. For the anti-HIV-1 activity the arginine-serine (RS)-rich domain of the splicing factor was dispensable. Based on our results, we propose that splicing at 3'ss A3 is dependent on binding of the enhancing SR proteins SRSF2 and SRSF6 to the ESE tat and ESE2 sequence. Mutational inactivation or interference specifically with ESE tat activity by LNA-directed masking seem to account for an early stage defect in viral gene expression, probably by cutting off the supply line of Tat that HIV needs to efficiently transcribe its genome.

The contribution of deep intronic splice-altering variants to hereditary breast and ovarian cancer (HBOC) is unknown. Current computational in silico tools to predict spliceogenic variants leading to pseudoexons have limited efficiency. We assessed the performance of the SpliceAI tool combined with ESRseq scores to identify spliceogenic deep intronic variants by affecting cryptic sites or splicing regulatory elements (SREs) using literature and experimental datasets. Our results with 233 published deep intronic variants showed that SpliceAI, with a 0.05 threshold, predicts spliceogenic deep intronic variants affecting cryptic splice sites, but is less effective in detecting those affecting SREs. Next, we characterized the SRE profiles using ESRseq, showing that pseudoexons are significantly enriched in SRE-enhancers compared to adjacent intronic regions. Although the combination of SpliceAI with ESRseq scores (considering ∆ESRseq and SRE landscape) showed higher sensitivity, the global performance did not improve because of the higher number of false positives. The combination of both tools was tested in a tumor RNA dataset with 207 intronic variants disrupting splicing, showing a sensitivity of 86%. Following the pipeline, five spliceogenic deep intronic variants were experimentally identified from 33 variants in HBOC genes. Overall, our results provide a framework to detect deep intronic variants disrupting splicing.

For more than three decades, researchers have known that consensus splice sites alone are not sufficient regulatory elements to provide complex splicing regulation. Other regulators, so-called splicing regulatory elements (SREs) are needed. Most importantly, their sequence variants often underlie the development of various human disorders. However, due to their variable location and high degeneracy, these regulatory sequences are also very difficult to recognize and predict. Many different approaches aiming to identify SREs have been tried, often leading to the development of in silico prediction tools. While these tools were initially expected to be helpful to identify splicing-affecting mutations in genetic diagnostics, we are still quite far from meeting this goal. In fact, most of these tools are not able to accurately discern the SRE-affecting pathological variants from those not affecting splicing. Nonetheless, several recent evaluations have given appealing results (namely for EX-SKIP, ESRseq and Hexplorer predictors). In this review, we aim to summarize the history of the different approaches to SRE prediction, and provide additional validation of these tools based on patients’ clinical data. Finally, we evaluate their usefulness for diagnostic settings and discuss the challenges that have yet to be met.

Extensive alternative splicing (AS) of precursor mRNAs (pre-mRNAs) in multicellular eukaryotes increases the protein-coding capacity of a genome and allows novel ways to regulate gene expression. In flowering plants, up to 48% of intron-containing genes exhibit AS. However, the full extent of AS in plants is not yet known, as only a few high-throughput RNA-Seq studies have been performed. As the cost of obtaining RNA-Seq reads continues to fall, it is anticipated that huge amounts of plant sequence data will accumulate and help in obtaining a more complete picture of AS in plants. Although it is not an onerous task to obtain hundreds of millions of reads using high-throughput sequencing technologies, computational tools to accurately predict and visualize AS are still being developed and refined. This review will discuss the tools to predict and visualize transcriptome-wide AS in plants using short-reads and highlight their limitations. Comparative studies of AS events between plants and animals have revealed that there are major differences in the most prevalent types of AS events, suggesting that plants and animals differ in the way they recognize exons and introns. Extensive studies have been performed in animals to identify cis-elements involved in regulating AS, especially in exon skipping. However, few such studies have been carried out in plants. Here, we review the current state of research on splicing regulatory elements (SREs) and briefly discuss emerging experimental and computational tools to identify cis-elements involved in regulation of AS in plants. The availability of curated alternative splice forms in plants makes it possible to use computational tools to predict SREs involved in AS regulation, which can then be verified experimentally. Such studies will permit identification of plant-specific features involved in AS regulation and contribute to deciphering the splicing code in plants.

Lynch syndrome (LS) is one of the most common hereditary cancer syndrome in human populations, associated with germline variants in MLH1 , MSH2/EPCAM , MSH6 and PMS2 genes. The advent of next generation sequencing has proven a significant impact in germline variant detection in the causative genes; however, a large proportion of patients with clinical criteria still receive uncertain or negative results. PMS2 is the least frequent reported gene, associated with up to 15% of LS cases with late-onset disease and low penetrance phenotype; however, the proportion of PMS2 -LS cases is considered to be highly underestimated. In this context, our analysis aimed to improve the current diagnostic yield by focusing on missense and intronic PMS2 variants available in public clinical databases (ClinVar, LOVD). We performed an in silico assessment of the wild-type DNA sequence and the reported genetic variants, employing splicing bioinformatics tools known for their effectiveness in other genes. Splicing variants were predicted in silico and using GTEx short-read RNA expression data. Out of the 2384 missense variants discovered, 90% were classified with uncertain significance (VUS). 4.9% of missense variants were shown to have a potential splicing consequence (DS > 0.2) using SpliceAI. As described in the original publication, SpliceAI-visual was proven effective in annotation of short intronic variants (< 50 bp). Four short intronic variants were identified using SpliceAI-visual as potentially splicing disturbing, in spite of using a lower threshold (DS > 0.1). Exons 2, 3, 4, 5, 6, 7, 8, 11, 12 and 14 were consistently predicted in at least three out of eight software with weak canonical splice sites. Additionally, we noted that both Exonic Splicing Enhancers (ESEs) and Exonic Splicing Silencers (ESSs) contribute significantly to alternative splicing and exonic selection in PMS2 gene. Specifically, ESE motifs were consistently more abundant in highly expressed exons 5, 11 and 14, while ESS motifs played a fundamental role in exons 6, 7 and 10. Computational analysis performed in our study serves as a valuable filtering step for guiding further RNA experiments. Additional functional data is necessary to validate our findings.

Accurate pre-mRNA splicing requires primary splicing signals, including the splice sites, a polypyrimidine tract, and a branch site, other splicing-regulatory elements (SREs). The SREs include exonic splicing enhancers (ESEs), exonic splicing silencers (ESSs), intronic splicing enhancers (ISEs), and intronic splicing silencers (ISSs), which are typically located near the splice sites. However, it is unclear to what extent splicing-driven selective pressure constrains exonic and intronic sequences, especially those distant from the splice sites. Here, we studied the distribution of SREs in human genes in terms of DNA strand-asymmetry patterns. Under a neutral evolution model, each mononucleotide or oligonucleotide should have a symmetric (Chargaff's second parity rule), or weakly asymmetric yet uniform, distribution throughout a pre-mRNA transcript. However, we found that large sets of unbiased, experimentally determined SREs show a distinct strand-asymmetry pattern that is inconsistent with the neutral evolution model, and reflects their functional roles in splicing. ESEs are selected in exons and depleted in introns and vice versa for ESSs. Surprisingly, this trend extends into deep intronic sequences, accounting for one third of the genome. Selection is detectable even at the mononucleotide level, so that the asymmetric base compositions of exons and introns are predictive of ESEs and ESSs. We developed a method that effectively predicts SREs based on strand asymmetry, expanding the current catalog of SREs. Our results suggest that human genes have been optimized for exon and intron discrimination through an RNA landscape shaped during evolution.

Transcription of the HIV-1 provirus generates a viral pre-mRNA, which is alternatively spliced into more than 50 HIV-1 mRNAs encoding all viral proteins. Regulation of viral alternative splice site usage includes the presence of splicing regulatory elements (SREs) which can dramatically impact RNA expression and HIV-1 replication when mutated. Recently, we were able to show that two viral SREs, GI3-2 and ESEtat, are important players in the generation of viral vif, vpr and tat mRNAs. Furthermore, we demonstrated that masking these SREs by transfected locked nucleic acid (LNA) mixmers affect the viral splicing pattern and viral particle production. With regard to the development of future therapeutic LNA mixmer-based antiretroviral approaches, we delivered the GI3-2 and the ESEtat LNA mixmers “nakedly”, without the use of transfection reagents (gymnosis) into HIV-1 infected cells. Surprisingly, we observed that gymnotically-delivered LNA mixmers accumulated in the cytoplasm, and seemed to co-localize with GW bodies and induced degradation of mRNAs containing their LNA target sequence. The GI3-2 and the ESEtat LNA-mediated RNA degradation resulted in abrogation of viral replication in HIV-1 infected Jurkat and PM1 cells as well as in PBMCs.

Background Cystinuria is an autosomal recessive disorder of dibasic amino acid transport in the kidney and the intestine leading to increased urinary cystine excretion and nephrolithiasis. Two genes, SLC3A1 and SLC7A9, coding respectively for rBAT and b0,+AT, account for the genetic basis of cystinuria. Methods This study reports the clinical and molecular characterization of a French cohort including 112 cystinuria patients and 25 relatives from 99 families. Molecular screening was performed using sequencing and Quantitative Multiplex PCR of Short Fluorescent Fragments analyses. Functional minigene‐based assays have been used to characterize splicing variants. Results Eighty‐eight pathogenic nucleotide changes were identified in SLC3A1 (63) and SLC7A9 (25) genes, of which 42 were novel. Interestingly, 17% (15/88) and 11% (10/88) of the total number of variants correspond, respectively, to large‐scale rearrangements and splicing mutations. Functional minigene‐based assays were performed for six variants located outside the most conserved sequences of the splice sites; three variants affect splice sites, while three others modify exonic splicing regulatory elements (ESR), in good agreement with a new in silico prediction based on ΔtESRseq values. Conclusion This report expands the spectrum of SLC3A1 and SLC7A9 variants and supports that digenic inheritance is unlikely. Furthermore, it highlights the relevance of assessing large‐scale rearrangements and splicing mutations to fully characterize cystinuria patients at the molecular level. This study addressing cystinuria adds inputs to the clinical knowledge, expands the spectrum of SLC3A1 and SLC7A9 variants and supports that digenic inheritance is unlikely. Furthermore, it highlights the relevance of assessing large‐scale rearrangements and splicing mutations to fully characterize cystinuria patients at the molecular level. Moreover, this work illustrates that special consideration should be given to the potential splicing effects of exonic variants, whatever their predicted consequences at the protein level (missense, synonymous, or even nonsense).