Transcriptomic changes in the microsporidia proliferation and host responses in congenitally infected embryos and larvae.

Congenital infection caused by vertical transmission of microsporidia N. bombycis can result in severe economic losses in the silkworm-rearing industry. Whole-transcriptome analyses have revealed non-coding RNAs and their regulatory networks in N. bombycis infected embryos and larvae. However, transcriptomic changes in the microsporidia proliferation and host responses in congenitally infected embryos and larvae remains unclear. Here, we simultaneously compared the transcriptomes of N. bombycis and its host B. mori embryos of 5-day and larvae of 1-, 5- and 10-day during congenital infection. For the transcriptome of N. bombycis, a comparison of parasite expression patterns between congenital-infected embryos and larva showed most genes related to parasite central carbon metabolism were down-regulated in larvae during infection, whereas the majority of genes involved in parasite proliferation and growth were up-regulated. Interestingly, a large number of distinct or shared differentially expressed genes (DEGs) were revealed by the Venn diagram and heat map, many of them were connected to infection related factors such as Ricin B lectin, spore wall protein, polar tube protein, and polysaccharide deacetylase. For the transcriptome of B. mori infected with N. bombycis, beyond numerous DEGs related to DNA replication and repair, mRNA surveillance pathway, RNA transport, protein biosynthesis, and proteolysis, with the progression of infection, a large number of DEGs related to immune and infection pathways, including phagocytosis, apoptosis, TNF, Toll-like receptor, NF-kappa B, Fc epsilon RI, and some diseases, were successively identified. In contrast, most genes associated with the insulin signaling pathway, 2-oxacarboxylic acid metabolism, amino acid biosynthesis, and lipid metabolisms were up-regulated in larvae compared to those in embryos. Furthermore, dozens of distinct and three shared DEGs that were involved in the epigenetic regulations, such as polycomb, histone-lysine-specific demethylases, and histone-lysine-N-methyltransferases, were identified via the Venn diagram and heat maps. Notably, many DEGs of host and parasite associated with lipid-related metabolisms were verified by RT-qPCR. Taken together, simultaneous transcriptomic analyses of both host and parasite genes lead to a better understanding of changes in the microsporidia proliferation and host responses in embryos and larvae in N. bombycis congenital infection.

Non-coding RNAs identification and regulatory networks in pathogen-host interaction in the microsporidia congenital infection.

BACKGROUND: The interaction networks between coding and non-coding RNAs (ncRNAs) including long non-coding RNA (lncRNA), covalently closed circular RNA (circRNA) and miRNA are significant to elucidate molecular processes of biological activities and interactions between host and pathogen. Congenital infection caused by vertical transmission of microsporidia N. bombycis can result in severe economic losses in the silkworm-feeding industry. However, little is known about ncRNAs that take place in the microsporidia congenital infection. Here we conducted whole-transcriptome RNA-Seq analyses to identify ncRNAs and regulatory networks for both N. bombycis and host including silkworm embryos and larvae during the microsporidia congenital infection. RESULTS: A total of 4,171 mRNAs, 403 lncRNA, 62 circRNAs, and 284 miRNAs encoded by N. bombycis were identified, among which some differentially expressed genes formed cross-talk and are involved in N. bombycis proliferation and infection. For instance, a lncRNA/circRNA competing endogenous RNA (ceRNA) network including 18 lncRNAs, one circRNA, and 20 miRNAs was constructed to describe 14 key parasites genes regulation, such as polar tube protein 3 (PTP3), ricin-B-lectin, spore wall protein 4 (SWP4), and heat shock protein 90 (HSP90). Regarding host silkworm upon N. bombycis congenital infection, a total of 14,889 mRNAs, 3,038 lncRNAs, 19,039 circRNAs, and 3,413 miRNAs were predicted based on silkworm genome with many differentially expressed coding and non-coding genes during distinct developmental stages. Different species of RNAs form interacting network to modulate silkworm biological processes, such as growth, metamorphosis and immune responses. Furthermore, a lncRNA/circRNA ceRNA network consisting of 140 lncRNAs, five circRNA, and seven miRNAs are constructed hypothetically to describe eight key host genes regulation, such as Toll-6, Serpin-6, inducible nitric oxide synthase (iNOS) and Caspase-8. Notably, cross-species analyses indicate that parasite and host miRNAs play a vital role in pathogen-host interaction in the microsporidia congenital infection. CONCLUSION: This is the first comprehensive pan-transcriptome study inclusive of both N. bombycis and its host silkworm with a specific focus on the microsporidia congenital infection, and show that ncRNA-mediated regulation plays a vital role in the microsporidia congenital infection, which provides a new insight into understanding the basic biology of microsporidia and pathogen-host interaction.