Dynamic miRNA–host gene co-expression and functional regulation in response to salinity fluctuations during biological invasions

Abstract

The significantly negative impact of marine invasive species underscores the need to understand the dynamics of invasion success. MicroRNAs (miRNAs) play a crucial role in regulating gene expression in response to stresses during invasions. Using the invasive tunicate Ciona robusta as a model, here we aim to study intragenic miRNA–host gene co-expression and functional regulation in response to recurrent salinity challenges. Despite genomic nestedness, only 9% of miRNA–host gene pairs showed significant co-expression (p < 0.05, correlation coefficient > 0). Recurring stresses dynamically altered the co-expression, revealing distinct miRNA–host gene expression at different stress times and stages. These differentially expressed miRNAs (p_adj < 0.05, |log₂foldchange|> 1) regulated biological processes, including free amino acid metabolism, water channel function, and ion transport to maintain osmotic homeostasis. These functional regulations were specific to time and stage, targeting the same type of osmolytes through varied pathways. Our findings highlight the diverse regulatory roles of miRNAs in enabling rapid responses to environmental stresses during invasions, providing new insights into miRNA-driven phenotypic plasticity under changing conditions.

Keywords

Intragenic miRNA / Biological invasion / Co-expression / Environmental challenge

References

[1] Agranat-Tamir L, Shomron N, Sperling J, Sperling R. Interplay between pre-mRNA splicing and microrna biogenesis within the supraspliceosome. Nucleic Acids Res. 2014, 42: 4640-4651.

[2] Baskerville S, Bartel DP. Microarray profiling of micrornas reveals frequent coexpression with neighboring mirnas and host genes. RNA. 2005, 11: 241-247.

[3] Bell ML, Buvoli M, Leinwand LA. Uncoupling of expression of an intronic microrna and its myosin host gene by exon skipping. Mol Cell Biol. 2010, 30: 1937-1945.

[4] Bereza D, Shenkar N. Shipping voyage simulation reveals abiotic barriers to marine bioinvasions. Sci Total Environ. 2022, 837: 155741.

[5] Bolger AM, Lohse M, Usadel B. Trimmomatic: a flexible trimmer for Illumina sequence data. Bioinformatics. 2014, 302114-2120.

[6] Bossdorf O, Richards CL, Pigliucci M. Epigenetics for ecologists. Ecol Lett. 2008, 11: 106-115.

[7] Briski E, Bailey SA, Casas-Monroy O, DiBacco C, Kaczmarska I, Lawrence JE, Leichsenring J, Levings C, MacGillivary ML, McKindsey CW, Nasmith LE, Parenteau M, Piercey GE, Rivkin RB, Rochon A, Roy S, Simard N, Sun B, Way C, Weise AMet al. . Taxon- and vector-specific variation in species richness and abundance during the transport stage of biological invasions. Limnol Oceanogr. 2013, 58: 1361-1372.

[8] Brunetti R, Gissi C, Pennati R, Caicci F, Gasparini F, Manni L. Morphological evidence that the molecularly determined Ciona intestinalis type A and type B are different species: Ciona robusta and Ciona intestinalis. J Zool Syst Evol Res. 2015, 53: 186-193.

[9] Caldas P, Luz M, Baseggio S, Andrade R, Sobral D, Grosso AR. Transcription readthrough is prevalent in healthy human tissues and associated with inherent genomic features. Commun Biol. 2024, 7: 100.

[10] Cao J, Shi F. Comparative analysis of the aquaporin gene family in 12 fish species. Animals. 2019, 9233.

[11] Caputi L, Andreakis N, Mastrototaro F, Cirino P, Vassillo M, Sordino P. Cryptic speciation in a model invertebrate chordate. Proc Natl Acad Sci USA. 2007, 1049364-9369.

[12] Chen Y, Shenkar N, Ni P, Lin Y, Li S, Zhan A. Rapid microevolution during recent range expansion to harsh environments. BMC Evol Biol. 2018, 18: 1-13.

[13] Chen Z, Huang X, Fu R, Zhan A. Neighbours matter: effects of genomic organization on gene expression plasticity in response to environmental stresses during biological invasions. Comp Biochem Physiol D Genomics Proteomics. 2022, 42: 100992.

[14] Chen Y, Gao Y, Zhang Z, Zhan A. Multi-omics inform invasion risks under global climate change. Glob Chang Biol. 2024, 30e17588.

[15] Crystal-Ornelas R, Hudgins EJ, Cuthbert RN, Haubrock PJ, Fantle-Lepczyk J, Angulo E, Kramer AM, Ballesteros-Mejia L, Leroy B, Leung B. Economic costs of biological invasions within North America. NeoBiota. 2021, 67485-510.

[16] Cuthbert RN, Pattison Z, Taylor NG, Verbrugge L, Diagne C, Ahmed DA, Leroy B, Angulo E, Briski E, Capinha C. Global economic costs of aquatic invasive alien species. Sci Total Environ. 2021, 775: 145238.

[17] Diagne C, Leroy B, Vaissière AC, Gozlan RE, Roiz D, Jarić I, Salles J, Bradshaw C, Courchamp F. High and rising economic costs of biological invasions worldwide. Nature. 2021, 592571-576.

[18] Friedländer MR, Mackowiak SD, Li N, Chen W, Rajewsky N. miRDeep2 accurately identifies known and hundreds of novel microRNA genes in seven animal clades. Nucleic Acids Res. 2012, 40: 37-52.

[19] Friedman RC, Farh KKH, Burge CB, Bartel DP. Most mammalian mRNAs are conserved targets of microRNAs. Genome Res. 2009, 19: 92-105.

[20] Fu R, Huang X, Chen Y, Chen Z, Zhan A. Interactive regulations of dynamic methylation and transcriptional responses to recurring environmental stresses during biological invasions. Front Mar Sci. 2021, 8: 800745.

[21] Fu R, Huang X, Zhan A. Identification of DNA (de) methylation-related genes and their transcriptional response to environmental challenges in an invasive model ascidian. Gene. 2021, 768145331.

[22] Gao X, Qiao Y, Han D, Zhang Y, Ma N. Enemy or partner: relationship between intronic micrornas and their host genes. IUBMB Life. 2012, 64835-840.

[23] Georgakilas G, Vlachos IS, Paraskevopoulou MD, Yang P, Zhang Y, Economides AN, Hatzigeorgiou AG. MicroTSS: accurate microRNA transcription start site identification reveals a significant number of divergent pri-miRNAs. Nat Commun. 2014, 5: 5700.

[24] He C, Li Z, Chen P, Huang H, Hurst LD, Chen J. Young intragenic miRNAs are less coexpressed with host genes than old ones: implications of miRNA–host gene coevolution. Nucleic Acids Res. 2012, 404002-4012.

[25] Hill MS, Vande ZP, Wittkopp PJ. Molecular and evolutionary processes generating variation in gene expression. Nat Rev Genet. 2021, 22: 203-215.

[26] Huang X, Zhan A. Highly dynamic transcriptional reprogramming and shorter isoform shifts under acute stresses during biological invasions. RNA Biol. 2021, 18: 340-353.

[27] Huang X, Li S, Ni P, Gao Y, Jiang B, Zhou Z, Zhan A. Rapid response to changing environments during biological invasions: DNA methylation perspectives. Mol Ecol. 2017, 26: 6621-6633.

[28] Huang X, Li S, Gao Y, Zhan A. Genome-wide identification, characterization and expression analyses of heat shock protein-related genes in a highly invasive ascidian Ciona savignyi. Front Physiol. 2018, 9: 1043.

[29] Huang X, Li H, Shenkar N, Zhan A. Multidimensional plasticity jointly contributes to rapid acclimation to environmental challenges during biological invasions. RNA. 2023, 29675-690.

[30] Huang X, Li H, Zhan A. Interplays between cis-and trans-acting factors for alternative splicing in response to environmental changes during biological invasions of ascidians. Int J Mol Sci. 2023, 24: 14921.

[31] Jong CJ, Sandal P, Schaffer SW. The role of taurine in mitochondria health: more than just an antioxidant. Molecules. 2021, 264913.

[32] Jonsson B, Jonsson N. Phenotypic plasticity and epigenetics of fish: embryo temperature affects later-developing life-history traits. Aquat Biol. 2019, 2821-32.

[33] Kim YK, Kim VN. Processing of intronic microRNAs. EMBO J. 2007, 26: 775-783.

[34] Kristensen H, Haldrup C, Strand S, Mundbjerg K, Mortensen MM, Thorsen K, Ostenfeld MS, Wild PJ, Arsov C, Goering W, Visakorpi T. Hypermethylation of the GABRE ~miR-452 ~miR-224 promoter in prostate cancer predicts biochemical recurrence after radical prostatectomy. Clin Cancer Res. 2014, 20: 2169-2181.

[35] Kushawaha AK, Khan A, Sopory SK, Sanan-Mishra N. Priming by high temperature stress induces microRNA regulated heat shock modules indicating their involvement in thermopriming response in rice. Life. 2021, 11: 291.

[36] Langmead B, Salzberg SL. Fast gapped-read alignment with Bowtie 2. Nat Methods. 2012, 9357-359.

[37] Li H, Huang X, Zhan A. Stress memory of recurrent environmental challenges in marine invasive species: Ciona robusta as a case study. Front Physiol. 2020, 11: 94.

[38] Li H, Huang X, Zhan A. Context-dependent antioxidant defense system (ADS)-based stress memory in response to recurrent environmental challenges in congeneric invasive species. Mar Life Sci Technol. 2024, 6315-330.

[39] Liu B, Shyr Y, Cai J, Liu Q. Interplay between miRNAs and host genes and their role in cancer. Briefings Funct Genom. 2019, 18: 255-266.

[40] Liu Y, Wen H, Qi X, Zhang X, Zhang K, Fan H, Tian Y, Hu Y, Li Y. Genome-wide identification of the Na⁺/H⁺ exchanger gene family in Lateolabrax maculatus and its involvement in salinity regulation. Comp Biochem Physiol D Genomics Proteomics. 2019, 29: 286-298.

[41] Liu B, Shyr Y, Liu Q. Pan-cancer analysis reveals common and specific relationships between intragenic miRNAs and their host genes. Biomedicines. 2021, 9: 1263.

[42] Liu X, Huang W, Liu Y, Zhan A. Perspectives of invasive alien species management in China. Ecol Appl. 2024, 34e2926.

[43] Lou F, Gao T, Han Z. Effect of salinity fluctuation on the transcriptome of the Japanese mantis shrimp Oratosquilla oratoria. Int J Biol Macromol. 2019, 1401202-1213.

[44] Love MI, Huber W, Anders S. Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol. 2014, 15: 1-21.

[45] Martin M. Cutadapt removes adapter sequences from high-throughput sequencing reads. Embnet J. 2011, 17: 10-12.

[46] Melamed Z, Levy A, Ashwal-Fluss R, Lev-Maor G, Mekahel K, Atias N, Gilad S, Sharan R, Levy C, Kadener S. Alternative splicing regulates biogenesis of miRNAs located across exon-intron junctions. Mol Cell. 2013, 50: 869-881.

[47] Monteys AM, Spengler RM, Wan J, Tecedor L, Lennox KA, Xing Y, Davidson BL. Structure and activity of putative intronic miRNA promoters. RNA. 2010, 16: 495-505.

[48] Pourmozaffar S, Tamadoni Jahromi S, Rameshi H, Sadeghi A, Bagheri T, Behzadi S, Gozari M, Zahedi M, Abrari L, Abrari Lazarjani S. The role of salinity in physiological responses of bivalves. Rev Aquacult. 2020, 12: 1548-1566.

[49] Pu C, Zhan A. Epigenetic divergence of key genes associated with water temperature and salinity in a highly invasive model ascidian. Biol Invasions. 2017, 19: 2015-2028.

[50] Pu M, Chen J, Tao Z, Miao L, Qi X, Wang Y, Ren J. Regulatory network of miRNA on its target: coordination between transcriptional and post-transcriptional regulation of gene expression. Cell Mol Life Sci. 2019, 76: 441-451.

[51] Ramalingam P, Palanichamy JK, Singh A, Das P, Bhagat M, Kassab MA, Sinha S, Chattopadhyay P. Biogenesis of intronic miRNAs located in clusters by independent transcription and alternative splicing. RNA. 2014, 20: 76-87.

[52] Rani V, Sengar RS. Biogenesis and mechanisms of microRNA-mediated gene regulation. Biotechnol Bioeng. 2022, 119: 685-692.

[53] Revelle W, Revelle MW. Package ‘psych’. Comprehensive R Archive Netw. 2015, 337: 161-165

[54] Richards CL, Bossdorf O, Muth NZ, Gurevitch J, Pigliucci M. Jack of all trades, master of some? On the role of phenotypic plasticity in plant invasions. Ecol Lett. 2006, 9: 981-993.

[55] Satou Y, Mineta K, Ogasawara M, Sasakura Y, Shoguchi E, Ueno K, Yamada L, Matsumoto J, Wasserscheid J, Dewar K. Improved genome assembly and evidence-based global gene model set for the chordate Ciona intestinalis: new insight into intron and operon populations. Genome Biol. 2008, 9: 1-11.

[56] Schanen BC, Li X. Transcriptional regulation of mammalian miRNA genes. Genomics. 2011, 97: 1-6.

[57] Shenkar N, Shmuel Y, Huchon D. The invasive ascidian Ciona robusta recorded from a Red Sea marina. Mar Biodivers. 2018, 48: 2211-2214.

[58] Simberloff D, Martin JL, Genovesi P, Maris V, Wardle DA, Aronson J, Courchamp F, Galil B, García-Berthou E, Pascal M. Impacts of biological invasions: what's what and the way forward. Trends Ecol Evol. 2013, 28: 58-66.

[59] Sokolov EP, Sokolova IM. Compatible osmolytes modulate mitochondrial function in a marine osmoconformer Crassostrea gigas (Thunberg, 1793). Mitochondrion. 2019, 45: 29-37.

[60] Solís FJT, Fernandes JM, Sarropoulou E, Monzón IF. Noncoding RNAs in fish physiology and development: miRNAs as a cornerstone in gene networks. 2022, Cellular and Molecular Approaches. In Fish Biology, Elsevier105159

[61] Sundaram GM, Common JE, Gopal FE, Srikanta S, Lakshman K, Lunny DP, Lim TC, Tanavde V, Lane EB, Sampath P. See-saw’expression of microRNA-198 and FSTL1 from a single transcript in wound healing. Nature. 2013, 495: 103-106.

[62] Towle DW, Henry RP, Terwilliger NB. Microarray-detected changes in gene expression in gills of green crabs (Carcinus maenas) upon dilution of environmental salinity. Comp Biochem Physiol D Genomics Proteomics. 2011, 6: 115-125.

[63] Turbelin AJ, Cuthbert RN, Essl F, Haubrock PJ, Ricciardi A, Courchamp F. Biological invasions are as costly as natural hazards. Perspect Ecol Conserv. 2023, 21: 143-150

[64] Veronese A, Visone R, Consiglio J, Acunzo M, Lupini L, Kim T, Ferracin M, Lovat F, Miotto E, Balatti V. Mutated β-catenin evades a microRNA-dependent regulatory loop. Proc Natl Acad Sci USA. 2011, 108: 4840-4845.

[65] Verri T, Terova G, Romano A, Barca A, Pisani P, Storelli C, Saroglia M. The SoLute carrier (SLC) family series in teleost fish. Funct Genom Aquacult. 2012, 24: 219-320.

[66] Wang J, Xiang G, Zhang K, Zhou Y. Expression signatures of intragenic miRNAs and their corresponding host genes in myeloid leukemia cells. Biotechnol Lett. 2012, 34: 2007-2015.

[67] Wickham H. ggplot2. Wiley Interdiscip Rev Comput Stat. 2011, 3: 180-185.

[68] Wu CT, Chiou CY, Chiu HC, Yang UC. Fine-tuning of microRNA-mediated repression of mRNA by splicing-regulated and highly repressive microRNA recognition element. BMC Genomics. 2013, 14: 1-12.

[69] Yan W, Fu R, Huang X, Zhan A. Complementary roles of DNA methylation and miRNA in regulating gene expression under environmental stress in biological invasions. Evol Appl. 2025, 18. e70178

[70] Yan W, Fu R, Huang X, Zhan A. Dynamic and functional miRNA-mediated gene regulations in response to recurrent environmental challenges during biological invasions. Mol Ecol. 2025, 34. e17749

[71] Yancey PH. Organic osmolytes as compatible, metabolic and counteracting cytoprotectants in high osmolarity and other stresses. J Exp Biol. 2005, 208: 2819-2830.

[72] Zeidler M, Hüttenhofer A, Kress M, Kummer KK. Intragenic microRNAs autoregulate their host genes in both direct and indirect ways - a cross-species analysis. Cells. 2020, 9: 232.

[73] Zhan A. Multi-omics-driven adaptive management of biological invasions: toward a proactive, predictive, and integrative framework. Biol Divers. 2025, 2163-192.

[74] Zhan A, Briski E, Bock DG, Ghabooli S, MacIsaac HJ. Ascidians as models for studying invasion success. Mar Biol. 2015, 162: 2449-2470.

RIGHTS & PERMISSIONS

The Author(s)