Regulatory Roles of Long Non-Coding RNAs Relevant to Antioxidant Enzymes and Immune Responses of Apis cerana Larvae Following Ascosphaera apis Invasion.

Long non-coding RNAs (lncRNAs) play an essential part in controlling gene expression and a variety of biological processes such as immune defense and stress-response. However, whether and how lncRNAs regulate responses of Apis cerana larvae to Ascosphaera apis invasion has remained unclear until now. Here, the identification and structural analysis of lncRNAs in the guts of A. cerana worker larvae were conducted, and the expression profile of larval lncRNAs during the A. apis infection process was then analyzed, followed by an investigation of the regulatory roles of differentially expressed lncRNAs (DElncRNAs) in the host response. In total, 76 sense lncRNAs, 836 antisense lncRNAs, 184 intron lncRNAs, 362 bidirectional lncRNAs, and 2181 intron lncRNAs were discovered in the larval guts. Additionally, 30 known and 9 novel lncRNAs were potential precursors for 36 and 11 miRNAs, respectively. In the three comparison groups, 386, 351, and 272 DElncRNAs were respectively identified, indicating the change in the overall expression pattern of host lncRNAs following the A. apis invasion. Analysis of cis-acting effect showed that DElncRNAs in the 4-, 5-, and 6-day-old comparison groups putatively regulated 55, 30, and 20 up- and down-stream genes, respectively, which were involved in a series of crucial functional terms and pathways, such as MAPK signaling pathway, and cell process. Analysis showed that 31, 8, and 11 DElncRNAs as potential antisense lncRNAs may interact with 26, 8, and 9 sense-strand mRNAs. Moreover, investigation of the competing endogenous RNA (ceRNA) network indicated that 148, 283, and 257 DElncRNAs were putatively regulated. The expression of target genes by targeting corresponding DEmiRNAs included those associated with antioxidant enzymes and immune responses. These results suggested that DElncRNAs played a potential part in the larval guts responding to the A. apis infection through a cis-acting manner and ceRNA mechanisms. Our findings deepen our understanding of interactions between A. cerana larvae and A. apis and offer a basis for clarifying the DElncRNA-mediated mechanisms underlying the host response to fungal invasion.

Diverse Regulatory Manners and Potential Roles of lncRNAs in the Developmental Process of Asian Honey Bee (Apis cerana) Larval Guts.

Long non-coding RNAs (lncRNAs) are crucial modulators in a variety of biological processes, such as gene expression, development, and immune defense. However, little is known about the function of lncRNAs in the development of Asian honey bee (Apis cerana) larval guts. Here, on the basis of our previously obtained deep-sequencing data from the 4-, 5-, and 6-day-old larval guts of A. cerana workers (Ac4, Ac5, and Ac6 groups), an in-depth transcriptome-wide investigation was conducted to decipher the expression pattern, regulatory manners, and potential roles of lncRNAs during the developmental process of A. cerana worker larval guts, followed by the verification of the relative expression of differentially expressed lncRNAs (DElncRNAs) and the targeting relationships within a competing endogenous RNA (ceRNA) axis. In the Ac4 vs. Ac5 and Ac5 vs. Ac6 comparison groups, 527 and 498 DElncRNAs were identified, respectively, which is suggestive of the dynamic expression of lncRNAs during the developmental process of larval guts. A cis-acting analysis showed that 330 and 393 neighboring genes of the aforementioned DElncRNAs were respectively involved in 29 and 32 functional terms, such as cellular processes and metabolic processes; these neighboring genes were also respectively engaged in 246 and 246 pathways such as the Hedgehog signaling pathway and the Wnt signaling pathway. Additionally, it was found that 79 and 76 DElncRNAs as potential antisense lncRNAs may, respectively, interact with 72 and 60 sense-strand mRNAs. An investigation of competing endogenous RNA (ceRNA) networks suggested that 75 (155) DElncRNAs in the Ac4 vs. Ac5 (Ac5 vs. Ac6) comparison group could target 7 (5) DEmiRNAs and further bind to 334 (248) DEmRNAs, which can be annotated to 33 (29) functional terms and 186 (210) pathways, including 12 (16) cellular- and humoral-immune pathways (lysosome pathway, necroptosis, MAPK signaling pathway, etc.) and 11 (10) development-associated signaling pathways (Wnt, Hippo, AMPK, etc.). The RT-qPCR detection of five randomly selected DElncRNAs confirmed the reliability of the used sequencing data. Moreover, the results of a dual-luciferase reporter assay were indicative of the binding relationship between MSTRG.11294.1 and miR-6001-y and between miR-6001-y and ncbi_107992440. These results demonstrate that DElncRNAs are likely to modulate the developmental process of larval guts via the regulation of the source genes' transcription, interaction with mRNAs, and ceRNA networks. Our findings not only yield new insights into the developmental mechanism underlying A. cerana larval guts, but also provide a candidate ceRNA axis for further functional dissection.

ame-miR-34 Modulates the Larval Body Weight and Immune Response of Apis mellifera Workers to Ascosphara apis Invasion.

MiRNAs are critical regulators of numerous physiological and pathological processes. Ascosphaera apis exclusively infects bee larvae and causes chalkbrood disease. However, the function and mechanism of miRNAs in the bee larval response to A. apis infection is poorly understood. Here, ame-miR-34, a previously predicted miRNA involved in the response of Apis mellifera larvae to A. apis invasion, was subjected to molecular validation, and overexpression and knockdown were then conducted to explore the regulatory functions of ame-miR-34 in larval body weight and immune response. Stem-loop RT-PCR and Sanger sequencing confirmed the authenticity of ame-miR-34 in the larval gut of A. mellifera. RT-qPCR results demonstrated that compared with that in the uninfected larval guts, the expression level of ame-miR-34 was significantly downregulated (p < 0.001) in the guts of A. apis-infected 4-, 5-, and 6-day-old larvae, indicative of the remarkable suppression of host ame-miR-34 due to A. apis infection. In comparison with the corresponding negative control (NC) groups, the expression level of ame-miR-34 in the larval guts in the mimic-miR-34 group was significantly upregulated (p < 0.001), while that in the inhibitor-miR-34 group was significantly downregulated (p < 0.01). Similarly, effective overexpression and knockdown of ame-miR-34 were achieved. In addition, the body weights of 5- and 6-day-old larvae were significantly increased compared with those in the mimic-NC group; the weights of 5-day-old larvae in the inhibitor-miR-34 group were significantly decreased in comparison with those in the inhibitor-NC group, while the weights of 4- and 6-day-old larvae in the inhibitor-miR-34 group were significantly increased, indicating the involvement of ame-miR-34 in modulating larval body weight. Furthermore, the expression levels of both hsp and abct in the guts of A. apis-infected 4-, 5-, and 6-day-old larvae were significantly upregulated after ame-miR-34 overexpression. In contrast, after ame-miR-34 knockdown, the expression levels of the aforementioned two key genes in the A. apis-infected 4-, 5-, and 6-day-old larval guts were significantly downregulated. Together, the results demonstrated that effective overexpression and knockdown of ame-miR-34 in both noninfected and A. apis-infected A. mellifera larval guts could be achieved by the feeding method, and ame-miR-34 exerted a regulatory function in the host immune response to A. apis invasion through positive regulation of the expression of hsp and abct. Our findings not only provide a valuable reference for the functional investigation of bee larval miRNAs but also reveal the regulatory role of ame-miR-34 in A. mellifera larval weight and immune response. Additionally, the results of this study may provide a promising molecular target for the treatment of chalkbrood disease.

Unveiling the circRNA-Mediated Immune Responses of Western Honey Bee Larvae to Ascosphaera apis Invasion.

Western honey bee (Apis mellifera), a eusocial insect with a superior economic and ecological value, is widely used in the beekeeping industry throughout the world. As a new class of non-coding RNAs (ncRNAs), circular RNAs (circRNAs) participate in the modulation of considerable biological processes, such as the immune response via diverse manners. Here, the identification, characteristic investigation, and molecular verification of circRNAs in the Apis mellifera ligustica larval guts were conducted, and the expression pattern of larval circRNAs during the Ascosphaera apis infection was analyzed, followed by the exploration of the potential regulatory part of differentially expressed circRNAs (DEcircRNAs) in host immune responses. A total of 2083 circRNAs in the larval guts of A. m. ligustcia were identified, with a length distribution ranging from 106 nt to 92,798 nt. Among these, exonic circRNAs were the most abundant type and LG1 was the most distributed chromosome. Additionally, 25, 14, and 30 up-regulated circRNAs as well as 26, 25, and 62 down-regulated ones were identified in the A. apis-inoculated 4-, 5-, and 6-day-old larval guts in comparison with the corresponding un-inoculated larval guts. These DEcircRNAs were predicted to target 35, 70, and 129 source genes, which were relative to 12, 23, and 20 GO terms as well as 11, 10, and 27 KEGG pathways, including 5 cellular and humoral immune pathways containing apoptosis, autophagy, endocytosis, MAPK, Toll, and Imd signaling pathways. Furthermore, complex competing endogenous RNA (ceRNA) regulatory networks were detected to be formed among DEcircRNAs, DEmiRNAs, and DEmRNAs. The Target DEmRNAs were engaged in 24, 20, and 25 functional terms as well as 62, 80, and 159 pathways, including several vital immune defense-associated pathways, namely the lysosome, endocytosis, phagosome, autophagy, apoptosis, MAPK, Jak-STAT, Toll, and Imd signaling pathways. Finally, back-splicing sites within 15 circRNAs and the difference in the 9 DEcircRNAs' expression between un-inoculated and A. apis-inoculated larval guts were confirmed utilizing molecular methods. These findings not only enrich our understanding of bee host-fungal pathogen interactions, but also lay a foundation for illuminating the mechanism underlying the DEcircRNA-mediated immune defense of A. m. ligustica larvae against A. apis invasion.

Effect of Ascosphaera apis Infestation on the Activities of Four Antioxidant Enzymes in Asian Honey Bee Larval Guts.

Ascosphaera apis infects exclusively bee larvae and causes chalkbrood, a lethal fungal disease that results in a sharp reduction in adult bees and colony productivity. However, little is known about the effect of A. apis infestation on the activities of antioxidant enzymes in bee larvae. Here, A. apis spores were purified and used to inoculate Asian honey bee (Apis cerana) larvae, followed by the detection of the host survival rate and an evaluation of the activities of four major antioxidant enzymes. At 6 days after inoculation (dpi) with A. apis spores, obvious symptoms of chalkbrood disease similar to what occurs in Apis mellifera larvae were observed. PCR identification verified the A. apis infection of A. cerana larvae. Additionally, the survival rate of larvae inoculated with A. apis was high at 1−2 dpi, which sharply decreased to 4.16% at 4 dpi and which reached 0% at 5 dpi, whereas that of uninoculated larvae was always high at 1~8 dpi, with an average survival rate of 95.37%, indicating the negative impact of A. apis infection on larval survival. As compared with those in the corresponding uninoculated groups, the superoxide dismutase (SOD) and catalase (CAT) activities in the 5- and 6-day-old larval guts in the A. apis−inoculated groups were significantly decreased (p < 0.05) and the glutathione S-transferase (GST) activity in the 4- and 5-day-old larval guts was significantly increased (p < 0.05), which suggests that the inhibition of SOD and CAT activities and the activation of GST activity in the larval guts was caused by A. apis infestation. In comparison with that in the corresponding uninoculated groups, the polyphenol oxidase (PPO) activity was significantly increased (p < 0.05) in the 5-day-old larval gut but significantly reduced (p < 0.01) in the 6-day-old larval gut, indicating that the PPO activity in the larval guts was first enhanced and then suppressed. Our findings not only unravel the response of A. cerana larvae to A. apis infestation from a biochemical perspective but also offer a valuable insight into the interaction between Asian honey bee larvae and A. apis.

Comprehensive investigation and regulatory function of lncRNAs engaged in western honey bee larval immune response to Ascosphaera apis invasion.

Ascosphaera apis is a fungal pathogen that exclusively infects bee larvae, causing chalkbrood disease, which results in severe damage for beekeeping industry. Long non-coding RNAs (lncRNAs) are versatile regulators in various biological processes such as immune defense and host-pathogen interaction. However, expression pattern and regulatory role of lncRNAs involved in immune response of bee host to A. apis invasion is still very limited. Here, the gut tissues of Apis mellifera ligustica 4-, 5-, and 6-day-old larvae inoculated by A. apis spores (AmT1, AmT2, and AmT3 groups) and corresponding un-inoculated larval guts (AmCK1, AmCK2, and AmCK3 groups) were prepared and subjected to deep sequencing, followed by identification of lncRNAs, analysis of differentially expressed lncRNAs (DElncRNAs), and investigation of competing endogenous RNA (ceRNA) network. In total, 3,746 A. m. ligustica lncRNAs were identified, including 78 sense lncRNAs, 891 antisense lncRNAs, 1,893 intergenic lncRNAs, 346 bidirectional lncRNAs, and 210 intronic lncRNAs. In the 4-, 5-, and 6- comparison groups, 357, 236, and 505 DElncRNAs were discovered. Additionally, 217, 129, and 272 DElncRNAs were respectively predicted to regulate neighboring genes via cis-acting manner, and these targets were associated with a series of GO terms and KEGG pathways of great importance, such as response to stimulus and Jak-STAT signaling pathway. Moreover, 197, 95, and 356 DElncRNAs were observed to target 10, eight, and 21 DEmiRNAs and further target 147, 79, and 315 DEmRNAs, forming complex regulatory networks. Further investigation suggested that these targets were engaged in several key cellular and humoral immune pathways, such as phagosome and MAPK signaling pathway. Ultimately, the expression trends of nine randomly selected DElncRNAs were verified by RT-qPCR, confirming the authenticity and reliability of our transcriptome data. Findings in this current work not only provide candidate DElncRNAs for functional study, but also lay a foundation for unclosing the mechanism underlying DElncRNA-regulated larval immune responses to A. apis invasion.