Interactions Among Host-Parasite MicroRNAs During Nosema ceranae Proliferation in Apis mellifera

How to overcome constraints imposed by microsporidian genome features to ensure gene prediction?

Since the advent of sequencing techniques and due to their continuous evolution, it has become easier and less expensive to obtain the complete genome sequence of any organism. Nevertheless, to elucidate all biological processes governing organism development, quality annotation is essential. In genome annotation, predicting gene structure is one of the most important and captivating challenges for computational biology. This aspect of annotation requires continual optimization, particularly for genomes as unusual as those of microsporidia. Indeed, this group of fungal-related parasites exhibits specific features (highly reduced gene sizes, sequences with high rate of evolution) linked to their evolution as intracellular parasites, requiring the implementation of specific annotation approaches to consider all these features. This review aimed to outline these characteristics and to assess the increasingly efficient approaches and tools that have enhanced the accuracy of gene prediction for microsporidia, both in terms of sensitivity and specificity. Subsequently, a final part will be dedicated to postgenomic approaches aimed at reinforcing the annotation data generated by prediction software. These approaches include the characterization of other understudied genes, such as those encoding regulatory noncoding RNAs or very small proteins, which also play crucial roles in the life cycle of these microorganisms.

In-depth investigation of microRNA-mediated cross-kingdom regulation between Asian honey bee and microsporidian

Asian honey bee Apis cerana is the original host for Nosema ceranae, a unicellular fungal parasite that causes bee nosemosis throughout the world. Currently, interaction between A. cerana and N. ceranae is largely unknown. Our group previously prepared A. c. cerana workers’ midguts at 7 days post inoculation (dpi) and 10 dpi with N. ceranae spores as well as corresponding un-inoculated workers’ midguts, followed by cDNA library construction and a combination of RNAs-seq and small RNA-seq. Meanwhile, we previously prepared clean spores of N. ceranae, which were then subjected to cDNA library construction and deep sequencing. Here, based on the gained high-quality transcriptome datasets, N. ceranae differentially expressed mRNAs (DEmiRNAs) targeted by host DEmiRNAs, and A. c. cerana DEmRNAs targeted by microsporidian DEmiRNAs were deeply investigated, with a focus on targets involved in N. ceranae glycolysis/glyconeogenesis as well as virulence factors, and A. c. cerana energy metabolism and immune response. In A. c. cerana worker’s midguts at 7 (10) dpi (days post inoculation), eight (seven) up-regulated and six (two) down-regulated miRNAs were observed to target 97 (44) down-regulated and 60 (15) up-regulated N. ceranae mRNAs, respectively. Additionally, two up-regulated miRNAs (miR-60-y and miR-676-y) in host midgut at 7 dpi could target genes engaged in N. ceranae spore wall protein and glycolysis/gluconeogenesis, indicating potential host miRNA-mediated regulation of microsporidian virulence factor and energy metabolism. Meanwhile, in N. ceranae at 7 (10) dpi, 121 (110) up-regulated and 112 (104) down-regulated miRNAs were found to, respectively, target 343 (247) down-regulated and 138 (110) down-regulated mRNAs in A. c. cerana workers’ midguts. These targets in host were relevant to several crucial cellular and humoral immune pathways, such as phagasome, endocytosis, lysosomes, regulation of autophagy, and Jak–STAT signaling pathway, indicative of the involvement of N. ceranae DEmiRNAs in regulating these cellular and humoral immune pathways. In addition, N. ceranae miR-21-x was up-regulated at 7 dpi and had a target relative to oxidative phosphorylation, suggesting that miR-21-x may be used as a weapon to modulate this pivotal energy metabolism pathway. Furthermore, potential targeting relationships between two pairs of host DEmiRNAs-microsporidian DEmRNAs and two pairs of microsporidian DEmiRNAs-host DEmRNAs were validated using RT-qPCR. Our findings not only lay a foundation for exploring the molecular mechanism underlying cross-kingdom regulation between A. c. cerana workers and N. ceranae, but also offer valuable insights into Asian honey bee-microsporidian interaction.

Recent Advances in the Biocontrol of Nosemosis in Honey Bees (Apis mellifera L.)

Nosemosis is a disease triggered by the single-celled spore-forming fungi Nosema apis and Nosema ceranae, which can cause extensive colony losses in honey bees (Apis mellifera L.). Fumagillin is an effective antibiotic treatment to control nosemosis, but due to its toxicity, it is currently banned in many countries. Accordingly, in the beekeeping sector, there is a strong demand for alternative ecological methods that can be used for the prevention and therapeutic control of nosemosis in honey bee colonies. Numerous studies have shown that plant extracts, RNA interference (RNAi) and beneficial microbes could provide viable non-antibiotic alternatives. In this article, recent scientific advances in the biocontrol of nosemosis are summarized.

Identification and Expression Profiling of Circulating MicroRNAs in Serum of Cysticercus pisiformis-Infected Rabbits

Cysticercus pisiformis (C. pisiformis), the larval form of Taenia pisiformis, parasitize mainly the liver, omentum and mesentery of rabbits and cause huge economic losses in the rabbit breeding industry. MicroRNA (miRNA), a short non-coding RNA, is widely and stably distributed in the plasma and serum. Numerous data demonstrates that, after parasitic infection, miRNAs become the key regulatory factor for controlling host biological processes. However, the roles of serum miRNAs in C. pisiformis-infected rabbits have not been elucidated. In this study, we compared miRNA expression profiles between the C. pisiformis-infected and healthy rabbit serum using RNA-seq. A total of 192 miRNAs were differentially expressed (fold change ≥ 2 and p < 0.05), including 79 up- and 113 downregulated miRNAs. These data were verified by qRT-PCR (real time quantitative polymerase chain reaction) analysis. Additionally, GO analysis showed that the target genes of these dysregulated miRNAs were most enriched in cellular, single-organism and metabolic processes. KEGG pathway analysis showed that these miRNAs target genes were involved in PI3K-Akt, viral carcinogenesis and B cell receptor signaling pathways. Interestingly, after aligning clean reads to the T. pisiformis genome, four (miR-124-3p_3, miR-124-3p_4, miR-124a and novel-miR1) T. pisiformis-derived miRNAs were found. Of these, novel-miR1was upregulated in different periods after C. pisiformis infection, which was verified qRT-PCR, and pre- novel-miR-1 was amplified from the cysticerci by RT-PCR, implying novel-miR-1 was derived from C. pisiformis and has great potential for the diagnosis of Cysticercosis pisiformis infection. This is the first investigation of miRNA expression profile and function in the serum of rabbits infected by C. pisiformis, providing fundamental data for developing diagnostic targets for Cysticercosis pisiformis.

Immunity and physiological changes in adult honey bees (Apis mellifera) infected with Nosema ceranae: The natural colony environment

Nosema ceranae is a microsporidium that infects Apis mellifera, causing diverse physiological and behavioral alterations. Given the existence of individual and social mechanisms to reduce infection and fungal spread in the colony, bees may respond differently to infection depending on their rearing conditions. In this study, we investigated the effect of N. ceranae in honey bee foragers naturally infected with different fungal loads in a tropical region. In addition, we explored the effects of N. ceranae artificially infected young bees placed in a healthy colony under field conditions. Honey bees naturally infected with higher loads of N. ceranae showed downregulation of genes from Toll and IMD immune pathways and antimicrobial peptide (AMP) genes, but hemolymph total protein amount and Vitellogenin (Vg) titers were not affected. Artificially infected bees spread N. ceranae to the controls in the colony, but fungal loads were generally lower than those observed in cages, probably because of social immunity. Although no significant changes in mRNA levels of AMP-encoding were observed, N. ceranae artificially infected bees showed downregulation of miR-989 (an immune-related microRNA), lower vitellogenin gene expression, and decreased hemolymph Vg titers. Our results demonstrate for the first time that natural infection by N. ceranae suppresses the immune system of honey bee foragers in the field. This parasite is detrimental to the immune system of young and old bees, and disease spread, mitigation and containment will depend on the colony environment.

Screening of Differentially Expressed Microsporidia Genes from Nosema ceranae Infected Honey Bees by Suppression Subtractive Hybridization

The microsporidium Nosema ceranae is a high prevalent parasite of the European honey bee (Apis mellifera). This parasite is spreading across the world into its novel host. The developmental process, and some mechanisms of N. ceranae-infected honey bees, has been studied thoroughly; however, few studies have been carried out in the mechanism of gene expression in N. ceranae during the infection process. We therefore performed the suppressive subtractive hybridization (SSH) approach to investigate the candidate genes of N. ceranae during its infection process. All 96 clones of infected (forward) and non-infected (reverse) library were dipped onto the membrane for hybridization. A total of 112 differentially expressed sequence tags (ESTs) had been sequenced. For the host responses, 20% of ESTs (13 ESTs, 10 genes, and 1 non-coding RNA) from the forward library and 93.6% of ESTs (44 ESTs, 28 genes) from the reverse library were identified as differentially expressed genes (DEGs) of the hosts. A high percentage of DEGs involved in catalytic activity and metabolic processes revealed that the host gene expression change after N. ceranae infection might lead to an unbalance of physiological mechanism. Among the ESTs from the forward library, 75.4% ESTs (49 ESTs belonged to 24 genes) were identified as N. ceranae genes. Out of 24 N. ceranae genes, nine DEGs were subject to real-time quantitative reverse transcription PCR (real-time qRT-PCR) for validation. The results indicated that these genes were highly expressed during N. ceranae infection. Among nine N. ceranae genes, one N. ceranae gene (AAJ76_1600052943) showed the highest expression level after infection. These identified differentially expressed genes from this SSH could provide information about the pathological effects of N. ceranae. Validation of nine up-regulated N. ceranae genes reveal high potential for the detection of early nosemosis in the field and provide insight for further applications.

Construction and heterologous overexpression of two chimeric proteins carrying outer hydrophilic loops of Vairimorpha ceranae and Nosema bombycis ATP/ADP carriers

Microsporidia Nosema bombycis and Vairimorpha ceranae cause destructive epizootics of honey bees and silkworms. Insufficient efficiency of the antibiotic fumagillin against V. ceranae, its toxicity and the absence of effective methods of N. bombycis treatment demand the discovery of novel strategies to suppress infections of domesticated insects. RNA interference is one such novel treatment strategy. Another one implies that the intracellular development of microsporidia may be suppressed by single-chain antibodies (scFv fragments) against functionally important parasite proteins. Important components of microsporidian metabolism are non-mitochondrial, plastidic-bacterial ATP/ADP carriers. These membrane transporters import host-derived ATP and provide the capacity to pathogens for energy parasitism. Here, we analyzed membrane topology of four V. ceranae and three N. bombycis ATP/ADP transporters to construct two fusion proteins carrying their outer hydrophilic loops contacting with infected host cell cytoplasm. Interestingly, full-size genes of N. bombycis transporters may be derived from the Asian swallowtail Papilio xuthus genome sequencing project. Synthesis of the artificial genes was followed by overexpression of recombinant proteins in E. coli as insoluble inclusion bodies. The gene fragments encoding the loops of individual transporters were also effectively expressed in bacteria. The chimeric antigens may be used to construct immune libraries or select microsporidia-suppressing scFv fragments from synthetic, semisynthetic, naïve and immune antibody libraries. A further expression of such antibodies in insect cells may increase their resistance to microsporidial infections.

Proliferation characteristics of the intracellular microsporidian pathogen Nosema bombycis in congenitally infected embryos

Nosema bombycis is an obligate intracellular pathogen that can be transmitted vertically from infected females to eggs, resulting in congenital infections in embryos. Here we investigated the proliferation characteristics of N. bombycis in silkworm embryos using a histopathological approach and deep RNA sequencing. We found that N. bombycis proliferated mainly around yolk granules at the early stage of the embryonic development, 1-2 days post oviposition (dpo). At 4-6 dpo, a portion of N. bombycis in different stages adjacent to the embryo were packaged into the newly formed intestinal lumen, while the remaining parasites continued to proliferate around yolk granules. In the newly hatched larvae (9 dpo), the newly formed spores accumulated in the gut lumen and immediately were released into the environment via the faeces. Transcriptional profiling of N. bombycis further confirmed multiplication of N. bombycis throughout every stage of embryonic development. Additionally, the increased transcriptional level of spore wall proteins and polar tube proteins from 4 dpo indicated an active formation of mature spores. Taken together, our results have provided a characterization of the proliferation of this intracellular microsporidian pathogen in congenitally infected embryos leading to vertical transmission.

Comparative Identification of MicroRNAs in Apis cerana cerana Workers' Midguts in Responseto Nosema ceranae Invasion

Here, the expression profiles and differentially expressed miRNAs (DEmiRNAs) in the midguts of Apis cerana cerana workers at 7 d and 10 d post-inoculation (dpi) with N. ceranae were investigated via small RNA sequencing and bioinformatics. Five hundred and twenty nine (529) known miRNAs and 25 novel miRNAs were identified in this study, and the expression of 16 predicted miRNAs was confirmed by Stem-loop RT-PCR. A total of 14 DEmiRNAs were detected in the midgut at 7 dpi, including eight up-regulated and six down-regulated miRNAs, while 12 DEmiRNAs were observed in the midgut at 10 dpi, including nine up-regulated and three down-regulated ones. Additionally, five DEmiRNAs were shared, while nine and seven DEmiRNAs were specifically expressed in midguts at 7 dpi and 10 dpi. Gene ontology analysis suggested some DEmiRNAs and corresponding target mRNAs were involved in various functions including immune system processes and response to stimulus. KEGG pathway analysis shed light on the potential functions of some DEmiRNAs in regulating target mRNAs engaged in material and energy metabolisms, cellular immunity and the humoral immune system. Further investigation demonstrated a complex regulation network between DEmiRNAs and their target mRNAs, with miR-598-y, miR-252-y, miR-92-x and miR-3654-y at the center. Our results can facilitate future exploration of the regulatory roles of miRNAs in host responses to N. ceranae, and provide potential candidates for further investigation of the molecular mechanisms underlying eastern honeybee-microsporidian interactions.

Genome-Wide Identification of Long Non-Coding RNAs and Their Regulatory Networks Involved in Apis mellifera ligustica Response to Nosema ceranae Infection

Long non-coding RNAs (lncRNAs) are a diverse class of transcripts that structurally resemble mRNAs but do not encode proteins, and lncRNAs have been proven to play pivotal roles in a wide range of biological processes in animals and plants. However, knowledge of expression patterns and potential roles of honeybee lncRNA response to Nosema ceranae infection is completely unknown. Here, we performed whole transcriptome strand-specific RNA sequencing of normal midguts of Apis mellifera ligustica workers (Am7CK, Am10CK) and N. ceranae-inoculated midguts (Am7T, Am10T), followed by comprehensive analyses using bioinformatic and molecular approaches. A total of 6353 A. m. ligustica lncRNAs were identified, including 4749 conserved lncRNAs and 1604 novel lncRNAs. These lncRNAs had minimal sequence similarities with other known lncRNAs in other species; however, their structural features were similar to counterparts in mammals and plants, including shorter exon and intron length, lower exon number, and lower expression level, compared with protein-coding transcripts. Further, 111 and 146 N. ceranae-responsive lncRNAs were identified from midguts at 7-days post-inoculation (dpi) and 10 dpi compared with control midguts. Twelve differentially expressed lncRNAs (DElncRNAs) were shared by Am7CK vs. Am7T and Am10CK vs. Am10T comparison groups, while the numbers of unique DElncRNAs were 99 and 134, respectively. Functional annotation and pathway analysis showed that the DElncRNAs may regulate the expression of neighboring genes by acting in cis and trans fashion. Moreover, we discovered 27 lncRNAs harboring eight known miRNA precursors and 513 lncRNAs harboring 2257 novel miRNA precursors. Additionally, hundreds of DElncRNAs and their target miRNAs were found to form complex competitive endogenous RNA (ceRNA) networks, suggesting that these DElncRNAs may act as miRNA sponges. Furthermore, DElncRNA-miRNA-mRNA networks were constructed and investigated, the results demonstrated that a portion of the DElncRNAs were likely to participate in regulating the host material and energy metabolism as well as cellular and humoral immune host responses to N. ceranae invasion. Our findings revealed here offer not only a rich genetic resource for further investigation of the functional roles of lncRNAs involved in the A. m. ligustica response to N. ceranae infection, but also a novel insight into understanding the host-pathogen interaction during honeybee microsporidiosis.

Scientific Advances in Controlling Nosema ceranae (Microsporidia) Infections in Honey Bees (Apis mellifera)

Honey bees (Apis mellifera) are agriculturally important pollinators that have been recently at risk to severe colony losses. A variety of parasites and pathogens have been linked to colony decline, including the microsporidian parasite Nosema ceranae. While fumagillin has been used to control nosemosis in managed honey bee colonies for decades, research shows that this antibiotic poses a toxic threat and that its efficacy against N. ceranae is uncertain. There is certainly a demand for a new veterinary medication to treat honey bee colonies infected with N. ceranae. In this review, recent scientific advances in controlling N. ceranae infections in honey bees are summarized.

Evolution of Dicer and Argonaute orthologs in microsporidian parasites

Microsporidia are a group of intracellular parasites which infect animal hosts. The infection can broadly influence the hosts' metabolism, growth as well as immune responses. Recently, a functional RNAi pathway was suggested from the microsporidia parasite Nosema ceranae, whereby the gene Dicer showed strong impact on spore proliferation. Based on sequenced microsporidian species, the RNAi gene orthologs have only been annotated for a few species. In order to study the selection of RNAi gene Dicer and Argonaute orthologs from microsporidian genomes, a phylogenetic analysis was performed based on single copy orthologs of 21 microsporidian parasite species. Of the 21 studied parasite species, 11 parasite species maintained Dicer and Argonaute orthologs, which were further used to build the gene trees. The gene Dicer and Argonaute orthologs were either both maintained or both lost. The topology structures between the phylogenetic Dicer, Argonaute and species trees were consistent. The results suggest that the gene Dicer and Argonaute were selected as a unit, which were selectively maintained/lost during the lineage divergence. The study provides general insights on the selection of RNAi pathway in microsporidian parasites and the evolution of parasitism.