{"status":"ok","message-type":"work","message-version":"1.0.0","message":{"indexed":{"date-parts":[[2026,9,17]],"date-time":"2026-09-17T09:04:48Z","timestamp":1789635888570,"version":"4.0.0"},"reference-count":47,"publisher":"Frontiers Media SA","license":[{"start":{"date-parts":[[2026,9,17]],"date-time":"2026-09-17T00:00:00Z","timestamp":1789603200000},"content-version":"vor","delay-in-days":0,"URL":"https:\/\/creativecommons.org\/licenses\/by\/4.0\/"}],"content-domain":{"domain":["frontiersin.org"],"crossmark-restriction":true},"short-container-title":["Front. Bioinform."],"abstract":"<jats:sec>\n                    <jats:title>Introduction<\/jats:title>\n                    <jats:p>Alternative splicing allows a single gene to produce multiple messenger RNA (mRNA) variants by differential RNA processing, resulting in the translation of distinct protein isoforms. Intron retention (IR) is a specific type of alternative splicing in which introns remain unspliced in the mature mRNA. The exact regulatory code behind intron splicing remains not fully deciphered, and unraveling these mechanisms is critical to uncovering the foundations of genetic disorders: a substantial fraction of disease-causing mutations disrupt intron splicing.<\/jats:p>\n                  <\/jats:sec>\n                  <jats:sec>\n                    <jats:title>Methods<\/jats:title>\n                    <jats:p>\n                      We applied explainable machine learning (xML) models to intronic sequences from four species of the ciliate genus Tetrahymena (\n                      <jats:italic>T. thermophila<\/jats:italic>\n                      ,\n                      <jats:italic>T. malaccensis<\/jats:italic>\n                      ,\n                      <jats:italic>T. borealis<\/jats:italic>\n                      and\n                      <jats:italic>T. elliotti<\/jats:italic>\n                      ), which allow the intron retention process to be analysed without tissue-specific confounders. Retained introns (RIs) and constitutively spliced introns (CSIs) were extracted through an automated Galaxy pipeline and described by several features of the intronic sequences, including the absence of repetitive nucleotide motifs \u2013 quantified as entropy \u2013 the GC content, the complexity of the secondary structures as estimated by the Lempel\u2013Ziv (LZ) measure, and the effective distance between the branch point and the 3\u2032 splice site. Four classifiers were compared by 10-fold cross-validation within the training set and evaluated on an independent test set.\n                    <\/jats:p>\n                  <\/jats:sec>\n                  <jats:sec>\n                    <jats:title>Results<\/jats:title>\n                    <jats:p>Across 289.6 k analysed introns, the C5.0 decision tree gave the best and most consistent performance, ranking first in three of the four species (accuracy 0.948\u20130.984). The key distinguishing features of retained introns were a higher entropy (Shannon entropy above 1.67, versus below 1.53 for constitutively spliced introns, against a theoretical maximum of 2) and more complex secondary structures near the 3\u2032 splice sites, together with a greater effective branch-point-to-3\u2032-splice-site distance and a relative enrichment of the GC dinucleotide.<\/jats:p>\n                  <\/jats:sec>\n                  <jats:sec>\n                    <jats:title>Discussion<\/jats:title>\n                    <jats:p>These features converge on weaker, less accessible splicing signals, which may impair splice-site recognition and bias the outcome towards retention. Our work identifies key sequence-level features predictive of intron retention in Tetrahymena and could provide a generalizable, explainable ML framework for investigating splicing regulation in other eukaryotes.<\/jats:p>\n                  <\/jats:sec>","DOI":"10.3389\/fbinf.2026.1889469","type":"journal-article","created":{"date-parts":[[2026,9,17]],"date-time":"2026-09-17T08:38:11Z","timestamp":1789634291000},"update-policy":"https:\/\/doi.org\/10.3389\/crossmark-policy","source":"Crossref","is-referenced-by-count":0,"title":["Entropy-driven machine learning for deciphering mRNA rearrangement by splicing in complex microorganisms"],"prefix":"10.3389","volume":"6","author":[{"given":"Alessio","family":"Mancini","sequence":"first","affiliation":[{"name":"School of Biosciences and Veterinary Medicine, Biosciences and Biotechnology Division, University of Camerino","place":["Camerino, Italy"]},{"name":"Clinical Laboratory Unit, Senigallia Hospital, Azienda Sanitaria Territoriale (AST)","place":["Ancona, Italy"]}],"role":[{"vocabulary":"crossref","role":"author"}]},{"given":"Emanuela","family":"Merelli","sequence":"additional","affiliation":[{"name":"School of Science and Technology, Computer Science Division, University of Camerino","place":["Camerino, Italy"]}],"role":[{"vocabulary":"crossref","role":"author"}]},{"given":"Marco","family":"Piangerelli","sequence":"additional","affiliation":[{"name":"School of Science and Technology, Computer Science Division, University of Camerino","place":["Camerino, Italy"]},{"name":"R&D Division, Vici & C. S.p.A., Santarcangelo di Romagna","place":["(RN), Italy"]}],"role":[{"vocabulary":"crossref","role":"author"}]},{"given":"Sandra","family":"Pucciarelli","sequence":"additional","affiliation":[{"name":"School of Biosciences and Veterinary Medicine, Biosciences and Biotechnology Division, University of Camerino","place":["Camerino, Italy"]}],"role":[{"vocabulary":"crossref","role":"author"}]},{"given":"Leonardo","family":"Vito","sequence":"additional","affiliation":[{"name":"School of Science and Technology, Computer Science Division, University of Camerino","place":["Camerino, Italy"]},{"name":"R&D Division, Logical System Srl","place":["Jesi (AN), Italy"]}],"role":[{"vocabulary":"crossref","role":"author"}]}],"member":"1965","published-online":{"date-parts":[[2026,9,17]]},"reference":[{"key":"B1","doi-asserted-by":"publisher","first-page":"W537","DOI":"10.1093\/nar\/gky379","article-title":"The galaxy platform for accessible, reproducible and collaborative biomedical analyses: 2018 update","volume":"46","author":"Afgan","year":"2018","journal-title":"Nucleic Acids Res."},{"key":"B2","doi-asserted-by":"publisher","first-page":"53","DOI":"10.1038\/nature09000","article-title":"Deciphering the splicing code","volume":"465","author":"Barash","year":"2010","journal-title":"Nature"},{"key":"B3","doi-asserted-by":"publisher","first-page":"2114","DOI":"10.1093\/bioinformatics\/btu170","article-title":"Trimmomatic: a flexible trimmer for illumina sequence data","volume":"30","author":"Bolger","year":"2014","journal-title":"Bioinformatics"},{"key":"B4","doi-asserted-by":"publisher","first-page":"1774","DOI":"10.1101\/gr.177790.114","article-title":"Widespread intron retention in mammals functionally tunes transcriptomes","volume":"24","author":"Braunschweig","year":"2014","journal-title":"Genome Res."},{"key":"B5","doi-asserted-by":"publisher","first-page":"321","DOI":"10.1613\/jair.953","article-title":"SMOTE: synthetic minority over-sampling technique","volume":"16","author":"Chawla","year":"2002","journal-title":"J. Artif. Intell. Res."},{"key":"B6","doi-asserted-by":"publisher","first-page":"e1001016","DOI":"10.1371\/journal.pcbi.1001016","article-title":"Genome-wide association between branch point properties and alternative splicing","volume":"6","author":"Corvelo","year":"2010","journal-title":"PLoS Comput. Biol."},{"key":"B7","doi-asserted-by":"publisher","first-page":"e0161476","DOI":"10.1371\/journal.pone.0161476","article-title":"Bimodal length distribution of ciliate introns reveals two distinct splice-site recognition modes","volume":"11","author":"Cs\u0171r\u00f6s","year":"2016","journal-title":"PLoS ONE"},{"key":"B8","doi-asserted-by":"publisher","first-page":"15","DOI":"10.1093\/bioinformatics\/bts635","article-title":"STAR: ultrafast universal RNA-seq aligner","volume":"29","author":"Dobin","year":"2013","journal-title":"Bioinformatics"},{"key":"B9","doi-asserted-by":"publisher","first-page":"5866","DOI":"10.1093\/hmg\/ddu309","article-title":"Multiple evidence strands suggest that there may be as few as 19,000 human protein-coding genes","volume":"23","author":"Ezkurdia","year":"2014","journal-title":"Hum. Mol. Genet."},{"key":"B10","doi-asserted-by":"publisher","first-page":"W297","DOI":"10.1093\/nar\/gkm311","article-title":"ASTALAVISTA: dynamic and flexible analysis of alternative splicing events in custom gene datasets","volume":"35","author":"Foissac","year":"2007","journal-title":"Nucleic Acids Res."},{"key":"B11","doi-asserted-by":"publisher","first-page":"771","DOI":"10.11517\/jjsai.14.5_771","article-title":"A short introduction to boosting","volume":"14","author":"Freund","year":"1999","journal-title":"J. Jpn. Soc. Artif. Intell."},{"key":"B12","doi-asserted-by":"publisher","first-page":"451","DOI":"10.1146\/annurev-cellbio-100617-062826","article-title":"Regulation of neuronal differentiation, function, and plasticity by alternative splicing","volume":"34","author":"Furlanis","year":"2018","journal-title":"Annu. Rev. Cell Dev. Biol."},{"key":"B13","doi-asserted-by":"publisher","first-page":"601","DOI":"10.1126\/science.aat5993","article-title":"How does temperature determine sex?","volume":"360","author":"Georges","year":"2018","journal-title":"Science"},{"key":"B14","doi-asserted-by":"publisher","first-page":"R1","DOI":"10.1186\/gb-2006-7-1-r1","article-title":"A class of human exons with predicted distant branch points revealed by analysis of AG dinucleotide exclusion zones","volume":"7","author":"Gooding","year":"2006","journal-title":"Genome Biol."},{"key":"B15","doi-asserted-by":"publisher","first-page":"994","DOI":"10.1093\/bioinformatics\/15.12.994","article-title":"On the complexity measures of genetic sequences","volume":"15","author":"Gusev","year":"1999","journal-title":"Bioinformatics"},{"key":"B16","doi-asserted-by":"publisher","first-page":"1988","DOI":"10.1038\/sj.emboj.7600677","article-title":"Intronic CA-repeat and CA-rich elements: a new class of regulators of mammalian alternative splicing","volume":"24","author":"Hui","year":"2005","journal-title":"EMBO J."},{"key":"B17","doi-asserted-by":"publisher","first-page":"2716","DOI":"10.1093\/nar\/gkz121","article-title":"Does co-transcriptional regulation of alternative splicing mediate plant stress responses?","volume":"47","author":"Jabre","year":"2019","journal-title":"Nucleic Acids Res."},{"key":"B18","doi-asserted-by":"publisher","first-page":"283","DOI":"10.1016\/S0168-9525(00)89076-9","article-title":"Dinucleotide relative abundance extremes: a genomic signature","volume":"11","author":"Karlin","year":"1995","journal-title":"Trends Genet."},{"key":"B19","doi-asserted-by":"publisher","first-page":"291","DOI":"10.1146\/annurev-biochem-060614-034316","article-title":"Mechanisms and regulation of alternative pre-mRNA splicing","volume":"84","author":"Lee","year":"2015","journal-title":"Annu. Rev. Biochem."},{"key":"B20","doi-asserted-by":"publisher","first-page":"128","DOI":"10.1186\/s12864-020-6541-0","article-title":"iREAD: a tool for intron retention detection from RNA-seq data","volume":"21","author":"Li","year":"2020","journal-title":"BMC Genomics"},{"key":"B21","doi-asserted-by":"publisher","first-page":"26","DOI":"10.1186\/1748-7188-6-26","article-title":"ViennaRNA package 2.0","volume":"6","author":"Lorenz","year":"2011","journal-title":"Algorithms Mol. Biol."},{"key":"B22","doi-asserted-by":"publisher","first-page":"113","DOI":"10.1038\/s41586-022-04633-0","article-title":"Intron-mediated induction of phenotypic heterogeneity","volume":"605","author":"Luka\u010di\u0161in","year":"2022","journal-title":"Nature"},{"key":"B23","doi-asserted-by":"publisher","first-page":"442","DOI":"10.1186\/s12859-018-2422-9","article-title":"CiliateGEM: an open-project and a tool for predictions of ciliate metabolic variations and experimental condition design","volume":"19","author":"Mancini","year":"2018","journal-title":"BMC Bioinforma."},{"key":"B24","doi-asserted-by":"publisher","first-page":"347","DOI":"10.1186\/s12859-020-03566-7","article-title":"Machine learning models predicting multidrug resistant urinary tract infections using \u201cDSaaS\u201d","volume":"21","author":"Mancini","year":"2020","journal-title":"BMC Bioinforma."},{"key":"B25","doi-asserted-by":"publisher","first-page":"e104049","DOI":"10.1371\/journal.pone.0104049","article-title":"Comparative analyses between retained introns and constitutively spliced introns in Arabidopsis thaliana using random forest and support vector machine","volume":"9","author":"Mao","year":"2014","journal-title":"PLoS One"},{"key":"B26","doi-asserted-by":"publisher","first-page":"4562","DOI":"10.1128\/MCB.17.8.4562","article-title":"G triplets located throughout a class of small vertebrate introns enforce intron borders and regulate splice site selection","volume":"17","author":"McCullough","year":"1997","journal-title":"Mol. Cell. Biol."},{"key":"B27","doi-asserted-by":"publisher","first-page":"9225","DOI":"10.1128\/MCB.20.24.9225-9235.2000","article-title":"An intronic splicing enhancer binds U1 snRNPs to enhance splicing and select 5\u2032 splice sites","volume":"20","author":"McCullough","year":"2000","journal-title":"Mol. Cell. Biol."},{"key":"B28","doi-asserted-by":"publisher","first-page":"92","DOI":"10.1007\/s10618-012-0295-5","article-title":"Training and assessing classification rules with imbalanced data","volume":"28","author":"Menardi","year":"2014","journal-title":"Data Min. Knowl. Discov."},{"key":"B29","doi-asserted-by":"publisher","first-page":"11497","DOI":"10.1093\/nar\/gkz1068","article-title":"The changing paradigm of intron retention: regulation, ramifications and recipes","volume":"47","author":"Monteuuis","year":"2019","journal-title":"Nucleic Acids Res."},{"key":"B30","doi-asserted-by":"publisher","first-page":"457","DOI":"10.1038\/nature08909","article-title":"Expansion of the eukaryotic proteome by alternative splicing","volume":"463","author":"Nilsen","year":"2010","journal-title":"Nature"},{"key":"B31","doi-asserted-by":"publisher","first-page":"W628","DOI":"10.1093\/nar\/gkh466","article-title":"Complexity: an internet resource for analysis of DNA sequence complexity","volume":"32","author":"Orlov","year":"2004","journal-title":"Nucleic Acids Res."},{"key":"B32","doi-asserted-by":"publisher","first-page":"18","DOI":"10.5120\/20639-3318","article-title":"C5.0 algorithm to improved decision tree with feature selection and reduced error pruning","volume":"117","author":"Pandya","year":"2015","journal-title":"Int. J. Comput. Appl."},{"key":"B33","doi-asserted-by":"publisher","first-page":"115965","DOI":"10.1016\/j.envpol.2020.115965","article-title":"Transcriptomic responses to silver nanoparticles in the freshwater unicellular eukaryote Tetrahymena thermophila","volume":"269","author":"Piersanti","year":"2021","journal-title":"Environ. Pollut."},{"key":"B34","doi-asserted-by":"publisher","first-page":"341","DOI":"10.1007\/978-1-62703-980-2_25","article-title":"Approaches to link RNA secondary structures with splicing regulation","volume":"1126","author":"Plass","year":"2014","journal-title":"Plant Epigenetics Epigenomics Methods Mol. Biol."},{"key":"B35","doi-asserted-by":"publisher","first-page":"e01287","DOI":"10.7554\/eLife.01287","article-title":"Regulated aggregative multicellularity in a close unicellular relative of metazoa","volume":"2","author":"Seb\u00e9-Pedr\u00f3s","year":"2013","journal-title":"eLife"},{"key":"B36","doi-asserted-by":"publisher","first-page":"920","DOI":"10.1093\/bioinformatics\/btx688","article-title":"Machine learning annotation of human branchpoints","volume":"34","author":"Signal","year":"2018","journal-title":"Bioinformatics"},{"key":"B37","doi-asserted-by":"publisher","first-page":"4939","DOI":"10.1128\/mcb.13.8.4939-4952.1993","article-title":"Scanning and competition between AGs are involved in 3\u2032 splice site selection in mammalian introns","volume":"13","author":"Smith","year":"1993","journal-title":"Mol. Cell. Biol."},{"key":"B38","doi-asserted-by":"publisher","first-page":"199","DOI":"10.1023\/B:STCO.0000035301.49549.88","article-title":"A tutorial on support vector regression","volume":"14","author":"Smola","year":"2004","journal-title":"Stat. Comput."},{"key":"B39","doi-asserted-by":"publisher","first-page":"562","DOI":"10.1038\/nprot.2012.016","article-title":"Differential gene and transcript expression analysis of RNA-seq experiments with TopHat and cufflinks","volume":"7","author":"Trapnell","year":"2012","journal-title":"Nat. Protoc."},{"key":"B40","doi-asserted-by":"publisher","first-page":"2100","DOI":"10.1038\/s41467-017-02007-z","article-title":"Molecular basis of differential 3\u2032 splice site sensitivity to anti-tumor drugs targeting U2 snRNP","volume":"8","author":"Vigevani","year":"2017","journal-title":"Nat. Commun."},{"key":"B41","doi-asserted-by":"publisher","first-page":"10","DOI":"10.1186\/s40708-024-00222-1","article-title":"Interpreting artificial intelligence models: a systematic review on the application of LIME and SHAP in Alzheimer\u2019s disease detection","volume":"11","author":"Vimbi","year":"2024","journal-title":"Brain Inf."},{"key":"B42","doi-asserted-by":"publisher","first-page":"802","DOI":"10.1261\/rna.876308","article-title":"Splicing regulation: from a parts list of regulatory elements to an integrated splicing code","volume":"14","author":"Wang","year":"2008","journal-title":"RNA"},{"key":"B43","doi-asserted-by":"publisher","first-page":"789","DOI":"10.1016\/j.tig.2022.03.017","article-title":"Intron retention: importance, challenges, and opportunities","volume":"38","author":"Wong","year":"2022","journal-title":"Trends Genet."},{"key":"B44","doi-asserted-by":"publisher","first-page":"1012","DOI":"10.1016\/j.molcel.2018.07.033","article-title":"Quantitative activity profile and context dependence of all human 5\u2032 splice sites","volume":"71","author":"Wong","year":"2018","journal-title":"Mol. Cell"},{"key":"B45","doi-asserted-by":"publisher","first-page":"e30630","DOI":"10.1371\/journal.pone.0030630","article-title":"Transcriptome analysis of the model protozoan, Tetrahymena thermophila, using deep RNA sequencing","volume":"7","author":"Xiong","year":"2012","journal-title":"PLoS ONE"},{"key":"B46","doi-asserted-by":"publisher","first-page":"377","DOI":"10.1089\/1066527041410418","article-title":"Maximum entropy modeling of short sequence motifs with applications to RNA splicing signals","volume":"11","author":"Yeo","year":"2004","journal-title":"J. Comput. Biol."},{"key":"B47","doi-asserted-by":"publisher","first-page":"439","DOI":"10.1126\/science.1191150","article-title":"Integrative modeling defines the nova splicing-regulatory network and its combinatorial controls","volume":"329","author":"Zhang","year":"2010","journal-title":"Science"}],"container-title":["Frontiers in Bioinformatics"],"original-title":[],"link":[{"URL":"https:\/\/www.frontiersin.org\/articles\/10.3389\/fbinf.2026.1889469\/full","content-type":"unspecified","content-version":"vor","intended-application":"similarity-checking"}],"deposited":{"date-parts":[[2026,9,17]],"date-time":"2026-09-17T08:38:13Z","timestamp":1789634293000},"score":1,"resource":{"primary":{"URL":"https:\/\/www.frontiersin.org\/articles\/10.3389\/fbinf.2026.1889469\/full"}},"subtitle":[],"short-title":[],"issued":{"date-parts":[[2026,9,17]]},"references-count":47,"alternative-id":["10.3389\/fbinf.2026.1889469"],"URL":"https:\/\/doi.org\/10.3389\/fbinf.2026.1889469","relation":{},"ISSN":["2673-7647"],"issn-type":[{"value":"2673-7647","type":"electronic"}],"subject":[],"published":{"date-parts":[[2026,9,17]]},"article-number":"1889469"}}