Cross-kingdom noncoding RNA regulation facilitates Nosema bombycis proliferation

Abstract

Microsporidia are obligate intracellular fungal parasites that are widely distributed in natural and agricultural environments and represent an important biological hazard to economically important insects. Nosema bombycis (N. bombycis), the causative agent of pébrine disease, poses a serious threat to silkworm breeding and sericultural sustainability. However, the molecular mechanisms underlying microsporidian proliferation and host-pathogen interactions remain poorly understood. In this study, we systematically characterized the long noncoding RNAs (lncRNAs) and microRNA-like RNAs (milRNAs) of N. bombycis, an important pathogen of Bombyx mori (B. mori). Using the full-length transcriptome and small RNA sequencing data from B. mori midgut tissues infected with N. bombycis at key developmental stages, we identified 30 novel lncRNAs and 12 previously unreported milRNAs in N. bombycis. We further established a competing endogenous RNA regulatory network involving N. bombycis-derived lncRNAs, parasite-specific milRNAs, and host-derived microRNAs (miRNAs), revealing intricate molecular interactions that may be implicated in pathogenicity. Notably, NbLNC2914, a highly expressed lncRNA, was shown to function as a molecular sponge for the host miRNA bmo-miR-2808a-3p, modulating its expression and consequently regulating the downstream N. bombycis gene, NBO_58g0005. Functional validation assays demonstrated that this regulatory axis had an important impact on N. bombycis proliferation. Our findings reveal a previously unrecognized noncoding RNA (ncRNA)-mediated regulatory axis that facilitates N. bombycis proliferation within host cells. This study provides mechanistic insights into the molecular basis of microsporidia pathogenicity and highlights ncRNAs as potential molecular targets for mitigating microsporidia-associated biological hazards in agricultural ecosystems.

Keywords

Bombyx mori / Nosema bombycis / NbLNC2914 / Bmo-miR-2808a-3p / NBO_58g0005

References

[1] L. Hou, M. Wang, L. Zhu, M. Ning, J. Bi, J. Du, X. Kong, W. Gu, Q. Meng, Full—length transcriptome sequencing and comparative transcriptome analysis of eriocheir sinensis in response to infection by the microsporidian hepatospora eriocheir, Front. Cell Infect. Microbiol. 12 (2022) 997574, doi: 10.3389/fcimb.2022.997574.

[2] N. Corradi, K.L. Haag, J.F. Pombert, D. Ebert, P.J. Keeling, Draft genome sequence of the daphnia pathogen octosporea bayeri: insights into the gene content of a large microsporidian genome and a model for host—parasite interactions, Genome Biol. 10 (2009) R106, doi: 10.1186/gb-2009-10-10-r106.

[3] K.L. Haag, J.I. Larsson, D. Refardt, D. Ebert, Cytological and molecular description of hamiltosporidium tvaerminnensis gen. Et sp. nov., a microsporidian parasite of daphnia magna, and establishment of hamiltosporidium magnivora comb. nov, Parasitology 138 (2011) 447-462, doi: 10.1017/S0031182010001393.

[4] E. Peyretaillade, H. El Alaoui, M. Diogon, V. Polonais, N. Parisot, D.G. Biron, P. Peyret, F. Delbac, Extreme reduction and compaction of microsporidian genomes, Res. Microbiol. 162 (2011) 598-606, doi: 10.1016/j.resmic.2011.03.004.

[5] P.J. Keeling, N. Corradi, H.G. Morrison, K.L. Haag, D. Ebert, L.M. Weiss, D.E. Akiyoshi, S. Tzipori, The reduced genome of the parasitic microsporidian enterocytozoon bieneusi lacks genes for core carbon metabolism, Genome Biol. Evol. 2 (2010) 304-309, doi: 10.1093/gbe/evq022.

[6] S. Nakjang, T.A. Williams, E. Heinz, A.K. Watson, P.G. Foster, K.M. Sendra, S.E. Heaps, R.P. Hirt, T.M. Embley, Reduction and expansion in microsporidian genome evolution: new insights from comparative genomics, Genome Biol. Evol. 5 (2013) 2285-2303, doi: 10.1093/gbe/evt184.

[7] N. Corradi, J.F. Pombert, L. Farinelli, E.S. Didier, P.J. Keeling, The complete sequence of the smallest known nuclear genome from the microsporidian encephalitozoon intestinalis, Nat. Commun. 1 (2010) 77, doi: 10.1038/ncomms1082.

[8] A. Belkorchia, J.F. Pombert, V. Polonais, N. Parisot, F. Delbac, J.F. Brugere, P. Peyret, C. Gaspin, E. Peyretaillade, Comparative genomics of microsporidian genomes reveals a minimal non—coding RNA set and new insights for transcription in minimal eukaryotic genomes, DNA Res. 24 (2017) 251-260, doi: 10.1093/dnares/dsx002.

[9] M.L. Kent, R.W. Shaw, J.L.J.M.P.O.O. Sanders, Microsporidia in fish, (2014) 493-520.

[10] F.J. Slack, A.M. Chinnaiyan, The role of non—coding RNAs in oncology, Cell 179 (2019) 1033-1055, doi: 10.1016/j.cell.2019.10.017.

[11] M. Ghildiyal, P.D. Zamore, Small silencing RNAs: an expanding universe, Nat. Rev. Genet. 10 (2009) 94-108, doi: 10.1038/nrg2504.

[12] M.S. Hameed, Y. Ren, M. Tuda, A. Basit, N. Urooj, Role of argonaute proteins in RNAi pathway in plutella xylostella: a review, Gene 903 (2024) 148195, doi: 10.1016/j.gene.2024.148195.

[13] F. Di Gesualdo, S. Capaccioli, M. Lulli, A pathophysiological view of the long non—coding RNA world, Oncotarget 5 (2014) 10976-10996, doi: 10.18632/oncotarget.2770.

[14] H.C. Lee, L. Li, W. Gu, Z. Xue, S.K. Crosthwaite, A. Pertsemlidis, Z.A. Lewis, M. Freitag, E.U. Selker, C.C. Mello, Y. Liu, Diverse pathways generate microRNA—like RNAs and dicer—independent small interfering RNAs in fungi, Mol. Cell. 38 (2010) 803-814, doi: 10.1016/j.molcel.2010.04.005.

[15] S.S. Chang, Z. Zhang, Y. Liu, RNA interference pathways in fungi: mechanisms and functions, Annu. Rev. Microbiol. 66 (2012) 305-323, doi: 10.1146/annurev-micro-092611-150138.

[16] J. Tang, X. Chen, Y. Yan, J. Huang, C. Luo, H. Tom, L. Zheng, Comprehensive transcriptome profiling reveals abundant long non—coding RNAs associated with development of the rice false smut fungus, ustilaginoidea virens, Env. Microbiol. 23 (2021) 4998-5013, doi: 10.1111/1462—2920.15432.

[17] S. Dhingra, Role of non—coding RNAs in fungal pathogenesis and antifungal drug responses, Curr. Clin. Microbiol. Rep. 7 (2020) 133-141, doi: 10.1007/s40588-020-00151-7.

[18] P. Till, R.L. Mach, A.R. Mach—Aigner, A current view on long noncoding RNAs in yeast and filamentous fungi, Appl. Microbiol. Biotechnol. 102 (2018) 7319-7331, doi: 10.1007/s00253-018-9187-y.

[19] Q. Huang, J.D. Evans, Identification of microRNA—like small RNAs from fungal parasite nosema ceranae, J. Invertebr. Pathol. 133 (2016) 107-109, doi: 10.1016/j.jip.2015.12.005.

[20] S.S. Shao, W.Y. Yan, Q. Huang, Identification of novel miRNAs from the microsporidian parasite nosema ceranae, Infect. Genet. Evol. 93 (2021) 104930, doi: 10.1016/j.meegid.2021.104930.

[21] Z. Shen, Q. Yang, L. Luo, T. Li, Z. Ke, T. Li, J. Chen, X. Meng, H. Xiang, C. Li, Z. Zhou, P. Chen, G. Pan, Non—coding RNAs identification and regulatory networks in pathogen—host interaction in the microsporidia congenital infection, BMC Genom. 24 (2023) 420, doi: 10.1186/s12864-023-09490-3.

[22] Z. Dong, N. Zheng, C. Hu, B. Deng, W. Fang, Q. Wu, P. Chen, X. Huang, N. Gao, C. Lu, M. Pan, Nosema bombycis microRNA—like RNA 8 (Nb—milR8) increases fungal pathogenicity by modulating BmPEX16 gene expression in its host, bombyx mori, Microbiol. Spectr. 9 (2021) e0104821, doi: 10.1128/Spectrum.01048-21.

[23] C. Franzen, A. Muller, Molecular techniques for detection, species differentiation, and phylogenetic analysis of microsporidia, Clin. Microbiol. Rev. 12 (1999) 243-285, doi: 10.1128/CMR.12.2.243.

[24] A. Cali, L.M. Weiss, P.M. Takvorian, Brachiola algerae spore membrane systems, their activity during extrusion, and a new structural entity, the multilayered interlaced network, associated with the polar tube and the sporoplasm, J. Eukaryot. Microbiol. 49 (2002) 164-174, doi: 10.1111/j.1550-7408.2002.tb00361.x.

[25] J. Luo, Q. He, J.Z. Xu, C. Xu, Y.Z. Han, H.L. Gao, X.Z. Meng, G.Q. Pan, T. Li, Z.Y. Zhou, Microsporidia infection upregulates host energy metabolism but maintains ATP homeostasis, J. Invertebr. Pathol. 186 (2021) 107596, doi: 10.1016/j.jip.2021.107596.

[26] M.H. Pan, X.J. Cai, M. Liu, J. Lv, H. Tang, J. Tan, C. Lu, Establishment and characterization of an ovarian cell line of the silkworm, bombyx mori, Tissue Cell. 42 (2010) 42-46, doi: 10.1016/j.tice.2009.07.002.

[27] Z. Dong, N. Gao, B. Deng, X. Huang, C. Hu, P. Chen, Q. Wu, C. Lu, M. Pan, Stable transformation of fluorescent proteins into nosema bombycis by electroporation, Parasit. Vectors 15 (2022) 141, doi: 10.1186/s13071-022-05236-4.

[28] D. Kim, G. Pertea, C. Trapnell, H. Pimentel, R. Kelley, S.L. Salzberg, TopHat2: accurate alignment of transcriptomes in the presence of insertions, deletions and gene fusions, Genome Biol. 14 (2013) R36, doi: 10.1186/gb-2013-14-4-r36.

[29] L. Kong, Y. Zhang, Z.Q. Ye, X.Q. Liu, S.Q. Zhao, L. Wei, G. Gao, CPC: assess the protein—coding potential of transcripts using sequence features and support vector machine, Nucleic. Acids. Res. 35 (2007) W345-W349, doi: 10.1093/nar/gkm391.

[30] L. Sun, H. Luo, D. Bu, G. Zhao, K. Yu, C. Zhang, Y. Liu, R. Chen, Y. Zhao, Utilizing sequence intrinsic composition to classify protein—coding and long non—coding transcripts, Nucleic. Acids. Res. 41 (2013) e166, doi: 10.1093/nar/gkt646.

[31] L. Wang, H.J. Park, S. Dasari, S. Wang, J.P. Kocher, W. Li, CPAT: coding—Potential assessment tool using an alignment—free logistic regression model, Nucleic. Acids. Res. 41 (2013) e74, doi: 10.1093/nar/gkt006.

[32] M. Punta, P.C. Coggill, R.Y. Eberhardt, J. Mistry, J. Tate, C. Boursnell, N. Pang, K. Forslund, G. Ceric, J. Clements, A. Heger, L. Holm, E.L. Sonnhammer, S.R. Eddy, A. Bateman, R.D. Finn, The pfam protein families database, Nucleic. Acids. Res. 40 (2012) D290-D301, doi: 10.1093/nar/gkr1065.

[33] B. Langmead, S.L. Salzberg, Fast gapped—read alignment with bowtie 2, Nat. Methods 9 (2012) 357-359, doi: 10.1038/nmeth.1923.

[34] R. Guo, D. Chen, C. Xiong, C. Hou, Y. Zheng, Z. Fu, Q. Liang, Q. Diao, L. Zhang, H. Wang, Z. Hou, D. Kumar, First identification of long non—coding RNAs in fungal parasite nosema ceranae, Apidologie (Celle) 49 (2018) 660-670, doi: 10.1007/s13592-018-0593-z.

[35] X. Xiong, C.J. Geden, Y. Tan, Y. Zhang, D. Zhang, J.H. Werren, X. Wang, Genome structure, evolution, and host shift of nosema, Biol. (Basel) 13 (2024) 952, doi: 10.3390/biology13110952.

[36] S. Pu, Y. Fang, Y. Yang, Q. Qu, M. Liu, J. Lian, X. Tang, Z. Shen, P. Qian, Identification of long non—coding RNAs in response to microsporidia infection in silkworm, Bombyx mori, J. Econ. Entomol. 117 (2024) 772-781, doi: 10.1093/jee/toae072.

[37] R.T. Chen, Y. Xiao, Z. Liu, L.L. Li, Y. Lu, P. Jiao, Y.G. Miao, Three vital RNA functions and interactions in the process of silk gland apoptosis in silkworm bombyx mori, Arch. Insect Biochem. Physiol. 100 (2019) e21511, doi: 10.1002/arch.21511.

[38] G. Pan, J. Xu, T. Li, Q. Xia, S.L. Liu, G. Zhang, S. Li, C. Li, H. Liu, L. Yang, T. Liu, X. Zhang, Z. Wu, W. Fan, X. Dang, H. Xiang, M. Tao, Y. Li, J. Hu, Z. Li, L. Lin, J. Luo, L. Geng, L. Wang, M. Long, Y. Wan, N. He, Z. Zhang, C. Lu, P.J. Keeling, J. Wang, Z. Xiang, Z. Zhou, Comparative genomics of parasitic silkworm microsporidia reveal an association between genome expansion and host adaptation, BMC Genom. 14 (2013) 186, doi: 10.1186/1471-2164-14-186.

[39] J. Xu, G. Pan, L. Fang, J. Li, X. Tian, T. Li, Z. Zhou, Z. Xiang, The varying microsporidian genome: existence of long—terminal repeat retrotransposon in domesticated silkworm parasite nosema bombycis, Int. J. Parasitol. 36 (2006) 1049-1056, doi: 10.1016/j.ijpara.2006.04.010.

[40] Y. Chen, J.S. Pettis, Y. Zhao, X. Liu, L.J. Tallon, L.D. Sadzewicz, R. Li, H. Zheng, S. Huang, X. Zhang, M.C. Hamilton, S.F. Pernal, A.P. Melathopoulos, X. Yan, J.D. Evans, Genome sequencing and comparative genomics of honey bee microsporidia, nosema apis reveal novel insights into host—parasite interactions, BMC Genom. 14 (2013) 451, doi: 10.1186/1471-2164-14-451.

[41] M. Rodriguez—Lopez, S. Anver, C. Cotobal, S. Kamrad, M. Malecki, C. Correia—Melo, M. Hoti, S. Townsend, S. Marguerat, S.K. Pong, M.Y. Wu, L. Montemayor, M. Howell, M. Ralser, J. Bahler, Functional profiling of long intergenic non—coding RNAs in fission yeast, Elife 11 (2022), doi: 10.7554/eLife.76000.

[42] J. Li, X. Liu, Z. Yin, Z. Hu, K.Q. Zhang, An overview on identification and regulatory mechanisms of long non—coding RNAs in fungi, Front. Microbiol. 12 (2021) 638617, doi: 10.3389/fmicb.2021.638617.

[43] H. Liu, M. Li, S. Cai, X. He, Y. Shao, X. Lu, Ricin—B—lectin enhances microsporidia nosema bombycis infection in BmN cells from silkworm bombyx mori, Acta Biochim. Biophys. Sin. (Shanghai) 48 (2016) 1050-1057, doi: 10.1093/abbs/gmw093.

[44] B. Zhang, H. Liang, H. Zou, J. Lu, M. Zhang, B. Liang, Comprehensive analysis of the lncRNAs, mRNAs, and miRNAs implicated in the immune response of pinctada fucata martensii to vibrio parahaemolyticus, Fish. Shellfish. Immunol. 130 (2022) 132-140, doi: 10.1016/j.fsi.2022.09.006.

[45] H.M. Glad, S.M. Tralamazza, D. Croll, The expression landscape and pangenome of long non—coding RNA in the fungal wheat pathogen zymoseptoria tritici, Microb. Genom. 9 (2023), doi: 10.1099/mgen.0.001136.

[46] B. Han, P.M. Takvorian, L.M. Weiss, Invasion of host cells by microsporidia, Front. Microbiol. 11 (2020) 172, doi: 10.3389/fmicb.2020.00172.