# Splicing regulation and intron evolution in the short-intron ciliate model of endosymbiosis Paramecium bursaria

**Authors:** Thi Ngan Giang Nguyen, Md Mostafa Kamal, Chien-Ling Lin, Jun-Yi Leu

PMC · DOI: 10.1093/nar/gkag063 · Nucleic Acids Research · 2026-02-02

## TL;DR

This study explores how splicing regulation and intron evolution help the ciliate Paramecium bursaria adapt during endosymbiosis with algae.

## Contribution

The research identifies splicing as a key mechanism in host adaptation during endosymbiosis and reveals evolutionary traits of short introns.

## Key findings

- Splicing, especially of 5′ proximal introns, enhances gene expression in symbiotic cells.
- 883 genes with differentially spliced introns are enriched in transmembrane transporters for nutrient exchange.
- Conserved introns show efficient splicing with lower GC content and uniform length.

## Abstract

The integration of symbionts into host cells during endosymbiosis significantly alters gene expression and cell physiology. Though alternative splicing facilitates cellular adaptation through rapid modulation of gene expression and protein isoform diversity, its regulatory role during endosymbiosis remains poorly understood. Paramecium bursaria, which harbors hundreds of Chlorella variabilis algae within its cytoplasm, offers a powerful model to study splicing during endosymbiosis, especially given its exceptionally short introns (median ∼24 nt). Using time-course RNA sequencing of symbiotic and aposymbiotic cells, we found that splicing, especially of 5′ proximal introns, enhances gene expression. Moreover, we identified 883 genes with differentially spliced introns, particularly enriched in transmembrane transporters essential for establishing nutrient exchange between a host cell and algal symbionts. Splicing regulation correlated with expression changes in conserved spliceosome components, implicating that these factors act as splicing enhancers or repressors during symbiosis. By exploring intron orthology across ciliates, we found that conserved introns exhibited more efficient splicing, characterized by lower GC content and uniform length, suggesting that intron evolution favors features that optimize expression. Our study reveals how splicing contributes to host adaptation during endosymbiosis and highlights the evolutionary dynamics of short introns in eukaryotes.

Graphical Abstract

## Linked entities

- **Species:** Paramecium bursaria (taxon 74790), Chlorella variabilis (taxon 554065)

## Full-text entities

- **Genes:** DDX23 (DEAD-box helicase 23) [NCBI Gene 9416] {aka PRPF28, SNRNP100, U5-100K, U5-100KD, prp28}, LINC02605 (long intergenic non-protein coding RNA 2605) [NCBI Gene 112935892] {aka AS, IL-7, IL-7-AS}, TXNL4A (thioredoxin like 4A) [NCBI Gene 10907] {aka BMKS, DIB1, DIM1, SNRNP15, TXNL4, U5-15kD}, SNRNP70 (small nuclear ribonucleoprotein U1 subunit 70) [NCBI Gene 6625] {aka RNPU1Z, RPU1, SNRP70, Snp1, U1-70K, U170K}, MYB (MYB proto-oncogene, transcription factor) [NCBI Gene 4602] {aka Cmyb, c-myb, c-myb_CDS, efg}, RPS11 (ribosomal protein S11) [NCBI Gene 6205] {aka S11, uS17}, DDX46 (DEAD-box helicase 46) [NCBI Gene 9879] {aka PRPF5, Prp5}, GTF2B (general transcription factor IIB) [NCBI Gene 2959] {aka TF2B, TFIIB}, LUC7L3 (LUC7 like 3 pre-mRNA splicing factor) [NCBI Gene 51747] {aka CRA, CREAP-1, CROP, LUC7A, OA48-18, hLuc7A}, CD2BP2 (CD2 cytoplasmic tail binding protein 2) [NCBI Gene 10421] {aka FWP010, LIN1, PPP1R59, Snu40, U5-52K}, SMARCA5 (SNF2 related chromatin remodeling ATPase 5) [NCBI Gene 8467] {aka ISWI, SNF2H, WCRF135, hISWI, hSNF2H}, RPL22 (ribosomal protein L22) [NCBI Gene 6146] {aka EAP, HBP15, HBP15/L22, L22, eL22}, PTBP1 (polypyrimidine tract binding protein 1) [NCBI Gene 5725] {aka HNRNP-I, HNRNPI, HNRPI, PTB, PTB-1, PTB-T}, HSPA4 (heat shock protein family A (Hsp70) member 4) [NCBI Gene 3308] {aka APG-2, HEL-S-5a, HS24/P52, HSPH2, RY, hsp70}, SMARCA1 (SNF2 related chromatin remodeling ATPase 1) [NCBI Gene 6594] {aka ISWI, NURF140, SNF2L, SNF2L1, SNF2LB, SNF2LT}, SF1 (splicing factor 1) [NCBI Gene 7536] {aka BBP, D11S636, MBBP, ZCCHC25, ZFM1, ZNF162}, LSM1 (LSM1 homolog, mRNA degradation associated) [NCBI Gene 27257] {aka CASM, FICUS, YJL124C}, HNRNPK (heterogeneous nuclear ribonucleoprotein K) [NCBI Gene 3190] {aka AUKS, CSBP, HNRPK, TUNP}, RPS20 (ribosomal protein S20) [NCBI Gene 6224] {aka S20, uS10}, GTF2H5 (general transcription factor IIH subunit 5) [NCBI Gene 404672] {aka C6orf175, TFB5, TFIIH, TGF2H5, TTD, TTD-A}, TBP (TATA-box binding protein) [NCBI Gene 6908] {aka GTF2D, GTF2D1, HDL4, SCA17, TBP1, TFIID}, SNRPC (small nuclear ribonucleoprotein polypeptide C) [NCBI Gene 6631] {aka U1C, Yhc1}, U2AF1 (U2 small nuclear RNA auxiliary factor 1) [NCBI Gene 7307] {aka FP793, RN, RNU2AF1, U2AF35, U2AFBP}, IGKV3-7 (immunoglobulin kappa variable 3-7 (non-functional)) [NCBI Gene 28915] {aka IGKV37, L10, L10a, Vh}, SNRPD1 (small nuclear ribonucleoprotein D1 polypeptide) [NCBI Gene 6632] {aka HsT2456, SMD1, SNRPD, Sm-D1}, LUC7L (LUC7 like) [NCBI Gene 55692] {aka LUC7B1, Luc7, SR+89, hLuc7B1}
- **Diseases:** A3SS (MESH:D009371), AS (MESH:C536589), DSIs (MESH:D012734)
- **Chemicals:** GC (MESH:C057580), cycloheximide (MESH:D003513), glutamine (MESH:D005973), Boston lettuce (-)
- **Species:** Paramecium caudatum (species) [taxon 5885], Stylonychia lemnae (species) [taxon 5949], Klebsiella pneumoniae (species) [taxon 573], Homo sapiens (human, species) [taxon 9606], Chlorella variabilis (species) [taxon 554065], Caenorhabditis elegans (species) [taxon 6239], Mus musculus (house mouse, species) [taxon 10090], Saccharomyces cerevisiae (baker's yeast, species) [taxon 4932], PX clade (clade) [taxon 569578], Drosophila melanogaster (fruit fly, species) [taxon 7227], Paramecium tetraurelia (species) [taxon 5888], Tetrahymena thermophila (species) [taxon 5911], Tetrahymena malaccensis (species) [taxon 5901], Ichthyophthirius multifiliis (species) [taxon 5932], Paramecium bursaria (species) [taxon 74790], C. elegans [taxon 328850], Danio rerio (leopard danio, species) [taxon 7955], Arabidopsis thaliana (mouse-ear cress, species) [taxon 3702], Hydra viridissima (green hydra, species) [taxon 6082], Stentor coeruleus (species) [taxon 5963]
- **Mutations:** GT-AG 5
- **Cell lines:** S2 — Drosophila melanogaster (Fruit fly), Spontaneously immortalized cell line (CVCL_Z232)

## Full text

_Full body text omitted from this summary view._ Fetch the complete paper as Markdown: https://tomesphere.com/paper/PMC12862383/full.md

## Figures

6 figures with captions in the complete paper: https://tomesphere.com/paper/PMC12862383/full.md

## References

119 references — full list in the complete paper: https://tomesphere.com/paper/PMC12862383/full.md

---
Source: https://tomesphere.com/paper/PMC12862383