A hybrid approach of ensemble learning and grey wolf optimizer for DNA splice junction prediction.

DNA splice junction classification is a crucial job in computational biology. The challenge is to predict the junction type (IE, EI, or N) from a given DNA sequence. Predicting junction type is crucial for understanding gene expression patterns, disease causes, splicing regulation, and gene structure. The location of the regions where exons are joined, and introns are removed during RNA splicing is very difficult to determine because no universal rule guides this process. This study presents a two-layer hybrid approach inspired by ensemble learning to overcome this challenge. The first layer applies the grey wolf optimizer (GWO) for feature selection. GWO's exploration ability allows it to efficiently search a vast feature space, while its exploitation ability refines promising areas, thus leading to a more reliable feature selection. The selected features are then fed into the second layer, which employs a classification model trained on the retrieved features. Using cross-validation, the proposed method divides the DNA splice junction dataset into training and test sets, allowing for a thorough examination of the classifier's generalization ability. The ensemble model is trained on various partitions of the training set and tested on the remaining held-out fold. This process is performed for each fold, comprehensively evaluating the classifier's performance. We tested our method using the StatLog DNA dataset. Compared to various machine learning models for DNA splice junction prediction, the proposed GWO+SVM ensemble method achieved an accuracy of 96%. This finding suggests that the proposed ensemble hybrid approach is promising for DNA splice junction classification. The implementation code for the proposed approach is available at https://github.com/EFHamouda/DNA-splice-junction-prediction.

A deep intronic splice-altering AIRE variant causes APECED syndrome through antisense oligonucleotide-targetable pseudoexon inclusion.

Autoimmune polyendocrinopathy-candidiasis-ectodermal dystrophy (APECED) is a life-threatening monogenic autoimmune disorder primarily caused by biallelic deleterious variants in the autoimmune regulator (AIRE) gene. We prospectively evaluated 104 patients with clinically diagnosed APECED syndrome and identified 17 patients (16%) from 14 kindreds lacking biallelic AIRE variants in exons or flanking intronic regions; 15 had Puerto Rican ancestry. Through whole-genome sequencing, we identified a deep intronic AIRE variant (c.1504-818 G>A) cosegregating with the disease in all 17 patients. We developed a culture system of AIRE-expressing primary patient monocyte-derived dendritic cells and demonstrated that c.1504-818 G>A creates a cryptic splice site and activates inclusion of a 109-base pair frame-shifting pseudoexon. We also found low-level AIRE expression in patient-derived lymphoblastoid cell lines (LCLs) and confirmed pseudoexon inclusion in independent extrathymic AIRE-expressing cell lines. Through protein modeling and transcriptomic analyses of AIRE-transfected human embryonic kidney 293 and thymic epithelial cell 4D6 cells, we showed that this variant alters the carboxyl terminus of the AIRE protein, abrogating its function. Last, we developed an antisense oligonucleotide (ASO) that reversed pseudoexon inclusion and restored the normal AIRE transcript sequence in LCLs. Thus, our findings revealed c.1504-818 G>A as a founder APECED-causing AIRE variant in the Puerto Rican population and uncovered pseudoexon inclusion as an ASO-reversible genetic mechanism underlying APECED.

Host factor RBMX2 promotes epithelial cell apoptosis by downregulating APAF-1's Retention Intron after Mycobacterium bovis infection.

The Mycobacterium tuberculosis variant bovis (M. bovis) is a highly pathogenic environmental microorganism that causes bovine tuberculosis (bTB), a significant zoonotic disease. Currently, "test and culling" is the primary measure for controlling bTB, but it has been proven to be inadequate in animals due to their high susceptibility to the pathogen. Selective breeding for increased host resistance to bTB to reduce its prevalence is feasible. In this study, we found a vital host-dependent factor, RBMX2, that can potentially promote M. bovis infection. By knocking RBMX2 out, we investigated its function during M. bovis infection. Through transcriptome sequencing and alternative splicing transcriptome sequencing, we concluded that after M. bovis infection, embryo bovine lung (EBL) cells were significantly enriched in RNA splicing associated with apoptosis compared with wild-type EBL cells. Through protein/molecular docking, molecular dynamics simulations, and real-time quantitative PCR, we demonstrated that RBMX2 promotes the apoptosis of epithelial cells by upregulating and binding to apoptotic peptidase activating factor 1 (APAF-1), resulting in the alternative splicing of APAF-1 as a retention intron. To our knowledge, this is the first report of M. bovis affecting host epithelial cell apoptosis by hijacking RBMX2 to promote the intron splicing of downstream APAF-1. These findings may represent a significant contribution to the development of novel TB prevention and control strategies.

In Vitro Self-Circularization Methods Based on Self-Splicing Ribozyme.

In vitro circular RNA (circRNA) preparation methods have been gaining a lot of attention recently as several reports suggest that circRNAs are more stable, with better performances in cells and in vivo, than linear RNAs in various biomedical applications. Self-splicing ribozymes are considered a major in vitro circRNA generation method for biomedical applications due to their simplicity and efficiency in the circularization of the gene of interest. This review summarizes, updates, and discusses the recently developed self-circularization methods based on the self-splicing ribozyme, such as group I and II intron ribozymes, and the pros and cons of each method in preparing circRNA in vitro.

Deep Intronic ETFDH Variants Represent a Recurrent Pathogenic Event in Multiple Acyl-CoA Dehydrogenase Deficiency.

Multiple acyl-CoA dehydrogenase deficiency (MADD) is a rare inborn error of metabolism affecting fatty acid and amino acid oxidation with an incidence of 1 in 200,000 live births. MADD has three clinical phenotypes: severe neonatal-onset with or without congenital anomalies, and a milder late-onset form. Clinical diagnosis is supported by urinary organic acid and blood acylcarnitine analysis using tandem mass spectrometry in newborn screening programs. MADD is an autosomal recessive trait caused by biallelic mutations in the ETFA, ETFB, and ETFDH genes encoding the alpha and beta subunits of the electron transfer flavoprotein (ETF) and ETF-coenzyme Q oxidoreductase enzymes. Despite significant advancements in sequencing techniques, many patients remain undiagnosed, impacting their access to clinical care and genetic counseling. In this report, we achieved a definitive molecular diagnosis in a newborn by combining whole-genome sequencing (WGS) with RNA sequencing (RNA-seq). Whole-exome sequencing and next-generation gene panels fail to detect variants, possibly affecting splicing, in deep intronic regions. Here, we report a unique deep intronic mutation in intron 1 of the ETFDH gene, c.35-959A>G, in a patient with early-onset lethal MADD, resulting in pseudo-exon inclusion. The identified variant is the third mutation reported in this region, highlighting ETFDH intron 1 vulnerability. It cannot be excluded that these intronic sequence features may be more common in other genes than is currently believed. This study highlights the importance of incorporating RNA analysis into genome-wide testing to reveal the functional consequences of intronic mutations.

Global impact of unproductive splicing on human gene expression.

Alternative splicing (AS) in human genes is widely viewed as a mechanism for enhancing proteomic diversity. AS can also impact gene expression levels without increasing protein diversity by producing 'unproductive' transcripts that are targeted for rapid degradation by nonsense-mediated decay (NMD). However, the relative importance of this regulatory mechanism remains underexplored. To better understand the impact of AS-NMD relative to other regulatory mechanisms, we analyzed population-scale genomic data across eight molecular assays, covering various stages from transcription to cytoplasmic decay. We report threefold more unproductive splicing compared with prior estimates using steady-state RNA. This unproductive splicing compounds across multi-intronic genes, resulting in 15% of transcript molecules from protein-coding genes being unproductive. Leveraging genetic variation across cell lines, we find that GWAS trait-associated loci explained by AS are as often associated with NMD-induced expression level differences as with differences in protein isoform usage. Our findings suggest that much of the impact of AS is mediated by NMD-induced changes in gene expression rather than diversification of the proteome.

Sequestration of DBR1 to stress granules promotes lariat intronic RNAs accumulation for heat-stress tolerance.

Heat stress (HS) poses a significant challenge to plant survival, necessitating sophisticated molecular mechanisms to maintain cellular homeostasis. Here, we identify SICKLE (SIC) as a key modulator of HS responses in Arabidopsis (Arabidopsis thaliana). SIC is required for the sequestration of RNA DEBRANCHING ENZYME 1 (DBR1), a rate-limiting enzyme of lariat intronic RNA (lariRNA) decay, into stress granules (SGs). The sequestration of DBR1 by SIC enhances the accumulation of lariRNAs, branched circular RNAs derived from excised introns during pre-mRNA splicing, which in turn promote the transcription of their parental genes. Our findings further demonstrate that SIC-mediated DBR1 sequestration in SGs is crucial for plant HS tolerance, as deletion of the N-terminus of SIC (SIC1-244) impairs DBR1 sequestration and compromises plant response to HS. Overall, our study unveils a mechanism of transcriptional regulation in the HS response, where lariRNAs are enriched through DBR1 sequestration, ultimately promoting the transcription of heat stress tolerance genes.

An intronic copy number variation in Syntaxin 17 determines speed of greying and melanoma incidence in Grey horses.

The Greying with age phenotype in horses involves loss of hair pigmentation whereas skin pigmentation is not reduced, and a predisposition to melanoma. The causal mutation was initially reported as a duplication of a 4.6 kb intronic sequence in Syntaxin 17. The speed of greying varies considerably among Grey horses. Here we demonstrate the presence of two different Grey alleles, G2 carrying two tandem copies of the duplicated sequence and G3 carrying three. The latter is by far the most common allele, probably due to strong selection for the striking white phenotype. Our results reveal a remarkable dosage effect where the G3 allele is associated with fast greying and high incidence of melanoma whereas G2 is associated with slow greying and low incidence of melanoma. The copy number expansion transforms a weak enhancer to a strong melanocyte-specific enhancer that underlies hair greying (G2 and G3) and a drastically elevated risk of melanoma (G3 only). Our direct pedigree-based observation of the origin of a G2 allele from a G3 allele by copy number contraction demonstrates the dynamic evolution of this locus and provides the ultimate evidence for causality of the copy number variation of the 4.6 kb intronic sequence.

Investigating Splice Defects in USH2A Using Targeted Long-Read Sequencing.

Biallelic variants in USH2A are associated with retinitis pigmentosa (RP) and Type 2 Usher Syndrome (USH2), leading to impaired vision and, additionally, hearing loss in the latter. Although the introduction of next-generation sequencing into clinical diagnostics has led to a significant uplift in molecular diagnostic rates, many patients remain molecularly unsolved. It is thought that non-coding variants or variants of uncertain significance contribute significantly to this diagnostic gap. This study aims to demonstrate the clinical utility of the reverse transcription-polymerase chain reaction (RT-PCR)-Oxford Nanopore Technology (ONT) sequencing of USH2A mRNA transcripts from nasal epithelial cells to determine the splice-altering effect of candidate variants. Five affected individuals with USH2 or non-syndromic RP who had undergone whole genome sequencing were recruited for further investigation. All individuals had uncertain genotypes in USH2A, including deep intronic rare variants, c.8682-654C>G, c.9055+389G>A, and c.9959-2971C>T; a synonymous variant of uncertain significance, c.2139C>T; p.(Gly713=); and a predicted loss of function duplication spanning an intron/exon boundary, c.3812-3_3837dup p.(Met1280Ter). In silico assessment using SpliceAI provided splice-altering predictions for all candidate variants which were investigated using ONT sequencing. All predictions were found to be accurate; however, in the case of c.3812-3_3837dup, the outcome was a complex cryptic splicing pattern with predominant in-frame exon 18 skipping and a low level of exon 18 inclusion leading to the predicted stop gain. This study detected and functionally characterised simple and complex mis-splicing patterns in USH2A arising from previously unknown deep intronic variants and previously reported variants of uncertain significance, confirming the pathogenicity of the variants.

Engineering Oncogenic Hotspot Mutations on SF3B1 via CRISPR-Directed PRECIS Mutagenesis.

SF3B1 is the most recurrently mutated RNA splicing gene in cancer. However, research of its pathogenic role has been hindered by a lack of disease-relevant cell line models. Here, our study compared four genome engineering platforms to establish SF3B1 mutant cell lines: CRISPR-Cas9 editing, AAV homology-directed repair editing, base editing (ABEmax, ABE8e), and prime editing (PE2, PE3, PE5max). We showed that prime editing via PE5max achieved the most efficient SF3B1 K700E editing across a wide range of cell lines. Our approach was further refined by coupling prime editing with a fluorescent reporter that leverages a SF3B1 mutation-responsive synthetic intron to mark successfully edited cells. By applying this approach, called prime editing coupled intron-assisted selection (PRECIS), we introduced the K700E hotspot mutation into two chronic lymphocytic leukemia cell lines, HG-3 and MEC-1. We demonstrated that our PRECIS-engineered cells faithfully recapitulate known mutant SF3B1 phenotypes, including altered splicing, copy number variations, and cell-growth defect. Moreover, we discovered that the SF3B1 mutation can cause the loss of Y chromosome in chronic lymphocytic leukemia. Our results showcase that PRECIS is an efficient and generalizable method for engineering genetically faithful SF3B1 mutant models. Our approach provides new insights on the role of SF3B1 mutation in cancer and enables the generation of SF3B1 mutant cell lines in relevant cellular context. SIGNIFICANCE: This study developed an approach that can reliably and efficiently engineer SF3B1 mutation into different cellular contexts, thereby revealing novel roles of SF3B1 mutation in driving aberrant splicing, clonal evolution, and genome instability.

Splice modulation strategy applied to deep intronic variants in COL7A1 causing recessive dystrophic epidermolysis bullosa.

Recessive dystrophic epidermolysis bullosa (RDEB) is a rare and most often severe genetic disease characterized by recurrent blistering and erosions of the skin and mucous membranes after minor trauma, leading to major local and systemic complications. The disease is caused by loss-of-function variants in COL7A1 encoding type VII collagen (C7), the main component of anchoring fibrils, which form attachment structures stabilizing the cutaneous basement membrane zone. Alterations in C7 protein structure and/or expression lead to abnormal, rare or absent anchoring fibrils resulting in loss of dermal-epidermal adherence and skin blistering. To date, more than 1,200 distinct COL7A1 deleterious variants have been reported and 19% are splice variants. Here, we describe two RDEB patients for whom we identified two pathogenic deep intronic pathogenic variants in COL7A1. One of these variants (c.7795-97C > G) promotes the inclusion of a pseudoexon between exons 104 and 105 in the COL7A1 transcript, while the other causes partial or complete retention of intron 51. We used antisense oligonucleotide (ASO) mediated exon skipping to correct these aberrant splicing events in vitro. This led to increased normal mRNA splicing above 94% and restoration of C7 protein expression at a level (up to 56%) that should be sufficient to reverse the phenotype. This first report of exon skipping applied to counteract deep intronic variants in COL7A1 represents a promising therapeutic strategy for personalized medicine directed at patients with intronic variants at a distance of consensus splice sites.

De Novo Deep Intron ELANE Mutation Resulting in Severe Congenital Neutropenia.

Severe congenital neutropenia (SCN) comprises a diverse range of rare hematological disorders characterized by recurrent, often life-threatening infections that manifest within the first months of life. Mutations in the ELANE gene are the most prevalent cause of SCN. While over 230 mutations in ELANE have been documented, including substitutions, frameshifts, nonsense mutations, and splice site alterations, the occurrence of deep intronic mutations has not been previously reported. Herein, we present the case of a young girl who exhibited recurrent fever, respiratory infections, skin abscesses, and gingivitis shortly after birth. Laboratory analysis revealed markedly diminished neutrophil levels alongside elevated monocyte and eosinophil counts. Bone marrow examination disclosed a halt in myelopoiesis maturation. ELANE gene full-length sequencing identified a novel de novo deep intron mutation in ELANE (c.598 + 79G > T), subsequently confirmed by Sanger sequencing. cDNA sequencing of the patient demonstrated aberrant gene splicing. Utilizing a mini-gene splicing assay for ELANE intronic variants, we identified a mutant ELANE allele (c.597 + 1_597 + 83ins) leading to the creation of a premature termination codon (p.Gly200ValfsTer40). Confocal microscopy revealed heightened expression of myeloperoxidase and neutrophil elastase in the patient, suggesting a potential role for the unfolded protein response in the pathogenesis of the deep intron ELANE mutation. In summary, our findings illustrate the first reported instance of de novo deep intron ELANE mutations associated with SCN, underscoring the importance of exploring deep intronic regions in SCN patients lacking identifiable disease-causing gene mutations.

Exploring Mitochondrial Heterogeneity and Evolutionary Dynamics in Thelephora ganbajun through Population Genomics.

Limited exploration in fungal mitochondrial genetics has uncovered diverse inheritance modes. The mitochondrial genomes are inherited uniparentally in the majority of sexual eukaryotes, our discovery of persistent mitochondrial heterogeneity within the natural population of the basidiomycete fungus Thelephora ganbajun represents a significant advance in understanding mitochondrial inheritance and evolution in eukaryotes. Here, we present a comprehensive analysis by sequencing and assembling the complete mitogenomes of 40 samples exhibiting diverse cox1 heterogeneity patterns from various geographical origins. Additionally, we identified heterogeneous variants in the nad5 gene, which, similar to cox1, displayed variability across multiple copies. Notably, our study reveals a distinct prevalence of introns and homing endonucleases in these heterogeneous genes. Furthermore, we detected potential instances of horizontal gene transfer involving homing endonucleases. Population genomic analyses underscore regional variations in mitochondrial genome composition among natural samples exhibiting heterogeneity. Thus, polymorphisms in heterogeneous genes, introns, and homing endonucleases significantly influence mitochondrial structure, structural variation, and evolutionary dynamics in this species. This study contributes valuable insights into mitochondrial genome architecture, population dynamics, and the evolutionary implications of mitochondrial heterogeneity in sexual eukaryotes.

MIPs: multi-locus intron polymorphisms in species identification and population genomics.

The study of species groups in which the presence of interspecific hybridization or introgression phenomena is known or suspected involves analysing shared bi-parentally inherited molecular markers. Current methods are based on different categories of markers among which the classical microsatellites or the more recent genome wide approaches for the analyses of thousands of SNPs or hundreds of microhaplotypes through high throughput sequencing. Our approach utilizes intron-targeted amplicon sequencing to characterise multi-locus intron polymorphisms (MIPs) and assess genetic diversity. These highly variable intron regions, combined with inter-specific transferable loci, serve as powerful multiple-SNP markers potentially suitable for various applications, from species and hybrid identification to population comparisons, without prior species knowledge. We developed the first panel of MIPs highly transferable across fish genomes, effectively distinguishing between species, even those closely related, and populations with different structures. MIPs offer versatile, hypervariable nuclear markers and promise to be especially useful when multiple nuclear loci must be genotyped across different species, such as for the monitoring of interspecific hybridization. Moreover, the relatively long sequences obtained ease the development of single-locus PCR-based diagnostic markers. This method, here demonstrated in teleost fishes, can be readily applied to other taxa, unlocking a new source of genetic variation.

A RORE-dependent Intronic Enhancer in the IL-7 Receptor-α Locus Controls Glucose Metabolism via Vγ4+ γδT17 Cells.

The IL-7R regulates the homeostasis, activation, and distribution of T cells in peripheral tissues. Although several transcriptional enhancers that regulate IL-7Rα expression in αß T cells have been identified, enhancers active in γδ T cells remain unknown. In this article, we discovered an evolutionarily conserved noncoding sequence (CNS) in intron 2 of the IL-7Rα-chain (IL-7Rα) locus and named this region CNS9. CNS9 contained a conserved retinoic acid receptor-related orphan receptor (ROR)-responsive element (RORE) and exerted RORγt-dependent enhancer activity in vitro. Mice harboring point mutations in the RORE in CNS9 (CNS9-RORmut) showed reduced IL-7Rα expression in IL-17-producing Vγ4+ γδ T cells. In addition, the cell number and IL-17A production of Vγ4+ γδ T cells were reduced in the adipose tissue of CNS9-RORmut mice. Consistent with the reduction in IL-17A, CNS9-RORmut mice exhibited decreased IL-33 expression in the adipose tissue, resulting in fewer regulatory T cells and glucose intolerance. The CNS9-ROR motif was partially responsible for IL-7Rα expression in RORγt+ regulatory T cells, whereas IL-7Rα expression was unaffected in RORγt-expressing Vγ2+ γδ T cells, Th17 cells, type 3 innate lymphoid cells, and invariant NKT cells. Our results indicate that CNS9 is a RORΕ-dependent, Vγ4+ γδ T cell-specific IL-7Rα enhancer that plays a critical role in adipose tissue homeostasis via regulatory T cells, suggesting that the evolutionarily conserved RORΕ in IL-7Rα intron 2 may influence the incidence of type 2 diabetes.

Functional Characterization of Splice Variants in the Diagnosis of Albinism.

Albinism is a genetically heterogeneous disease in which 21 genes are known so far. Its inheritance mode is autosomal recessive except for one X-linked form. The molecular analysis of exonic sequences of these genes allows for about a 70% diagnostic rate. About half (15%) of the unsolved cases are heterozygous for one pathogenic or probably pathogenic variant. Assuming that the missing variant may be located in non-coding regions, we performed sequencing for 122 such heterozygous patients of either the whole genome (27 patients) or our NGS panel (95 patients) that includes, in addition to all exons of the 21 genes, the introns and flanking sequences of five genes, TYR, OCA2, SLC45A2, GPR143 and HPS1. Rare variants (MAF < 0.01) in trans to the first variant were tested by RT-PCR and/or minigene assay. Of the 14 variants tested, nine caused either exon skipping or the inclusion of a pseudoexon, allowing for the diagnosis of 11 patients. This represents 9.8% (12/122) supplementary diagnosis for formerly unsolved patients and 75% (12/16) of those in whom the candidate variant was in trans to the first variant. Of note, one missense variant was demonstrated to cause skipping of the exon in which it is located, thus shedding new light on its pathogenic mechanism. Searching for non-coding variants and testing them for an effect on RNA splicing is warranted in order to increase the diagnostic rate.

Intra-Phenotypic and -Genotypic Variations of Beauveria bassiana (Bals.) Vuill. Strains Infecting Aedes aegypti L. Adults.

Beauveria bassiana has potential for Aedes aegypti biological control. However, its efficacy depends on the strain's geographic location, host susceptibility, and virulence. The present study aimed to evaluate the effectiveness of B. bassiana strain BBPTG4 conidia in controlling Ae. aegypti adults and its detection via introns profile on exposed mosquito corpses. Morphologic characteristics among strains were highly similar. Comprehensive testing of these strains demonstrated that BBPT4 exhibited the ideal biological activity for Ae. aegypti control, with a median lethal time (TL50) of 7.5 d compared to ~3 d and ~10 d for BB01 and BB37 strains, respectively. Infected mosquitoes died after GHA and BBPTG4 exposure, and corpses were analyzed for infecting strains detection. Differences among the seven evaluated strains were determined, assessing five different insertion group I intron profiles in BBTG4, BB01, GHA, BB37, and BB02 strains. Mosquitoes infected by BBPTG4 and non-exposed (negative control) intron profiles were obtained. We detected the presence of introns in the BBPTG4 strain, which were not present in non-exposed mosquitoes. In conclusion, B. bassiana strains showed similarities in terms of their cultural and microscopic morphological characteristics and biologicals virulence level, but different intron profiles. BBPTG4 strain-infected Ae. aegypti adult corpses, showing specific amplicons, enabled us to identify B. bassiana at the strain level among infected mosquitoes. However, monitoring and detection of field-infected insects is essential for further verification.

Intron Retention of DDX39A Driven by SNRPD2 is a Crucial Splicing Axis for Oncogenic MYC/Spliceosome Program in Hepatocellular Carcinoma.

RNA splicing is a dynamic molecular process in response to environmental stimuli and is strictly regulated by the spliceosome. Sm proteins, constituents of the spliceosome, are key components that mediate splicing reactions; however, their potential role in hepatocellular carcinoma (HCC) is poorly understood. In the study, SNRPD2 (PD2) is found to be the most highly upregulated Sm protein in HCC and to act as an oncogene. PD2 modulates DDX39A intron retention together with HNRNPL to sustain the DDX39A short variant (39A_S) expression. Mechanistically, 39A_S can mediate MYC mRNA nuclear export to maintain high MYC protein expression, while MYC in turn potentiates PD2 transcription. Importantly, digitoxin can directly interact with PD2 and has a notable cancer-suppressive effect on HCC. The study reveals a novel mechanism by which DDX39A senses oncogenic MYC signaling and undergoes splicing via PD2 to form a positive feedback loop in HCC, which can be targeted by digitoxin.