Identification of long non-coding RNAs in the chalkbrood disease pathogen Ascospheara apis

Dynamics and regulatory role of circRNAs in Asian honey bee larvae following fungal infection

Non-coding RNA (ncRNA) plays a vital part in the regulation of immune responses, growth, and development in plants and animals. Here, the identification, characteristic analysis, and molecular verification of circRNAs in Apis cerana cerana worker larval guts were conducted, followed by in-depth investigation of the expression pattern of larval circRNAs during Ascosphaera apis infection and exploration of the potential regulatory part of differentially expressed circRNAs (DEcircRNAs) in host immune responses. A total of 3178 circRNAs in the larval guts of A. c. cerana were identified, with a length distribution ranging from 15 to 96,007 nt. Additionally, 155, 95, and 86 DEcircRNAs were identified in the in the 4-, 5-, and 6-day-old larval guts following A. apis infection. These DEcircRNAs were predicted to target 29, 25, and 18 parental genes relevant to 12, 20, and 17 GO terms as well as 144, 114, and 61 KEGG pathways, including 5 cellular and 4 humoral immune pathways. Complex competing endogenous RNA (ceRNA) regulatory networks were detected as being formed among DEcircRNAs, DEmiRNAs, and DEmRNAs. The target DEmRNAs were engaged in 36, 47, and 47 GO terms as well as 331, 332, and 331 pathways, including 6 cellular and 6 humoral immune pathways. Further, 19 DEcircRNAs, 5 DEmiRNAs, and 3 mRNAs were included in the sub-networks relative to 3 antioxidant enzymes. Finally, back-splicing sites within 15 circRNAs and the difference in the 15 DEcircRNAs' expression between uninoculated and A. apis-inoculated larval guts were confirmed based on 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. c. cerana larvae against A. apis invasion. KEY POINTS: • The expression pattern of circRNAs was altered in the A. cerana worker larval guts following A. apis infection. • Back-splicing sites within 15 A. cerana circRNAs were verified using molecular approaches. DEcircRNAs potentially modulated immune responses and antioxidant enzymes in A. apis-challenged host guts.

First Characterization and Regulatory Function of piRNAs in the Apis mellifera Larval Response to Ascosphaera apis Invasion

piRNAs are a class of small non-coding RNAs that play essential roles in modulating gene expression and abundant biological processes. To decode the piRNA-regulated larval response of western honeybees (Apis mellifera) to Ascosphaera apis infection, the expression pattern of piRNAs in Apis mellifera ligustica larval guts after A. apis inoculation was analyzed based on previously obtained high-quality small RNA-seq datasets, followed by structural characterization, target prediction, regulatory network investigation, and functional dissection. Here, 504, 657, and 587 piRNAs were respectively identified in the 4-, 5-, and 6-day-old larval guts after inoculation with A. apis, with 411 ones shared. These piRNAs shared a similar length distribution and first base bias with mammal piRNAs. Additionally, 96, 103, and 143 DEpiRNAs were detected in the 4-, 5-, and 6-day-old comparison groups. Targets of the DEpiRNAs were engaged in diverse pathways such as the phosphatidylinositol signaling system, inositol phosphate metabolism, and Wnt signaling pathway. These targets were involved in three energy metabolism-related pathways, eight development-associated signaling pathways, and seven immune-relevant pathways such as the Jak-STAT signaling pathway. The expression trends of five randomly selected DEpiRNAs were verified using a combination of RT-PCR and RT-qPCR. The effective overexpression and knockdown of piR-ame-945760 in A. apis-infected larval guts were achieved by feeding a specific mimic and inhibitor. Furthermore, piR-ame-945760 negatively regulated the expression of two target immune mRNAs, SOCS5 and ARF1, in the larval gut during the A. apis infection. These findings indicated that the overall expression level of piRNAs was increased and the expression pattern of piRNAs in larval guts was altered due to the A. apis infection, DEpiRNAs were putative regulators in the A. apis-response of A. m. ligustica worker larvae. Our data provide not only a platform for the functional investigation of piRNAs in honeybees, especially in bee larvae, but also a foundation for illuminating the piRNA-involved mechanisms underlying the host response to the A. apis infection.

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.

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.

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.

Integrative Analysis of lncRNA-mRNA Co-expression Provides Novel Insights Into the Regulation of Developmental Transitions in Female Varroa destructor

Varroa destructor is a major pathogenic driver of the Western honeybee colony losses globally. Understanding the developmental regulation of V. destructor is critical to develop effective control measures. Development is a complex biological process regulated by numerous genes and long non-coding RNAs (lncRNAs); however, the underlying regulation of lncRNAs in the development of V. destructor remains unknown. In this study, we analyzed the RNA sequencing (RNA-Seq) data derived from the four stages of female V. destructor in the reproductive phase (i.e., egg, protonymph, deutonymph, and adult). The identified differentially expressed mRNAs and lncRNAs exhibited a stage-specific pattern during developmental transitions. Further functional enrichment established that fat digestion and absorption, ATP-binding cassette (ABC) transporters, mitogen-activated protein kinase (MAPK) signaling pathway, and ubiquitin-proteasome pathway play key roles in the maturation of female V. destructor. Moreover, the lncRNAs and mRNAs of some pivotal genes were significantly upregulated at the deutonymph stage, such as cuticle protein 65/6.4/63/38 and mucin 5AC, suggesting that deutonymph is the key stage of metamorphosis development and pathogen resistance acquisition for female V. destructor. Our study provides novel insights into a foundational understanding of V. destructor biology.

Identification of long noncoding RNAs reveals the effects of dinotefuran on the brain in Apis mellifera (Hymenopptera: Apidae)

Background

Dinotefuran (CAS No. 165252–70-0), a neonicotinoid insecticide, has been used to protect various crops against invertebrate pests and has been associated with numerous negative sublethal effects on honey bees. Long noncoding RNAs (lncRNAs) play important roles in mediating various biological and pathological processes, involving transcriptional and gene regulation. The effects of dinotefuran on lncRNA expression and lncRNA function in the honey bee brain are still obscure.

Results

Through RNA sequencing, a comprehensive analysis of lncRNAs and mRNAs was performed following exposure to 0.01 mg/L dinotefuran for 1, 5, and 10 d. In total, 312 lncRNAs and 1341 mRNAs, 347 lncRNAs and 1458 mRNAs, and 345 lncRNAs and 1155 mRNAs were found to be differentially expressed (DE) on days 1, 5 and 10, respectively. Gene set enrichment analysis (GSEA) indicated that the dinotefuran-treated group showed enrichment in carbohydrate and protein metabolism and immune-inflammatory responses such as glycine, serine and threonine metabolism, pentose and glucuronate interconversion, and Hippo and transforming growth factor-β (TGF-β) signaling pathways. Moreover, the DE lncRNA TCONS_00086519 was shown by fluorescence in situ hybridization (FISH) to be distributed mainly in the cytoplasm, suggesting that it may serve as a competing endogenous RNA and a regulatory factor in the immune response to dinotefuran.

Conclusion

This study characterized the expression profile of lncRNAs upon exposure to neonicotinoid insecticides in young adult honey bees and provided a framework for further study of the role of lncRNAs in honey bee growth and the immune response.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07811-y.

LncRNA LINC00963 promotes osteogenic differentiation of hBMSCs and alleviates osteoporosis progression by targeting miRNA-760/ETS1 axis

Although long non-coding RNA LINC00963 has been reported to play a crucial regulatory role in osteoporosis (OP), its specific mechanism has not been well studied. Cell viability of human bone marrow mesenchymal stem cells (hBMSCs) transfected with short hairpin RNA targeting LINC00963 (sh-LINC00963) and negative control (sh-NC) was analysed by cell counting kit-8 (CCK-8) assay. Alkaline phosphatase (ALP) activity in hBMSCs transfected with sh-LINC00963 and sh-NC after induction by osteogenic medium (OM) on day 7 was detected. The protein expression levels of osteocalcin (OCN) and osteopontin (OPN) in hBMSCs transfected with sh-LINC00963 and sh-NC during OM induction on day 3 were detected by western blot. The relationship among LINC00963, miR-760, and E26 transformation specific-1 (ETS1) was determined by bioinformatics analysis, luciferase reporter assay, and RNA-binding protein immunoprecipitation (RIP) assay. A rat model with OP was established to confirm the role of LINC00963 in vivo. The expression level of LINC00963 was much lower in hBMSCs isolated from the discarded femoral head tissues of OP patients compared with that in health patients. Meanwhile, the expression level of LINC00963 was significantly increased and the expression level of miR-760 was decreased in hBMSCs during osteogenic induction. LINC00963 could bind to the 3'-untranslated region (3'-UTR) of miR-760 and negatively regulate the expression of miR-760, then promote the osteogenic differentiation in hBMSCs. ETS1 was identified as a target of miR-760. Moreover, overexpression of LINC00963 obviously reduced bone mineral density (BMD) of the left femur in OP rats and alleviated OP progression in vivo. Our results demonstrated that LINC00963 positively regulated the expression of ETS1 by directly targeting miR-760, and then promoted osteogenic differentiation of hBMSCs in vitro, and also attenuated OP progression in vivo, suggesting that LINC00963 might be a potential therapeutic target for OP.

The Effect of Artificial Media and Temperature on the Growth and Development of the Honey Bee Brood Pathogen Ascosphaera apis

Ascosphaera apis is a causative agent of chalkbrood, which is one of the most widespread honey bee diseases. In our experiments, the influence of several artificial media and cultivation under different temperatures was evaluated. Concretely, the radial growth of separated mating types was measured, reproductive structures in a Neubauer hemocytometer chamber were counted simultaneously, and the morphometry of spore cysts and spore balls was assessed. The complex set of experiments determined suitable cultivation conditions. A specific pattern between reproductive structure size and temperature was found. The optimal temperature for both mating types was 30 °C. SDA and YGPSA media are suitable for fast mycelial growth. Moreover, the effect of bee brood on fungus growth and development in vitro was investigated by modification of culture medium. The newly modified medium PDA-BB4 was most effective for the production of the reproductive structures. The result suggests that honey bee brood provides necessary nutrients for proper fungus development during in vitro cultivation. As there is no registered therapeutic agent against chalkbrood in most countries, including the European Union, the assessment of A. apis growth and development in different conditions could help to understand fungus pathogenesis and thus control chalkbrood disease.

Genome-Wide Identification of Long Non-coding RNAs in the Gravid Ectoparasite Varroa destructor

Long non-coding RNAs (lncRNAs) emerge as critical regulators across a wide variety of biological functions in living organisms. However, to date, no systematic characterization of lncRNAs has been investigated in the ectoparasitic mite Varroa destructor, the most severe biotic threat to honey bees worldwide. Here, we performed an initial genome-wide identification of lncRNAs in V. destructor via high-throughput sequencing technology and reported, for the first time, the transcriptomic landscape of lncRNAs in the devastating parasite. By means of a lncRNA identification pipeline, 6,645 novel lncRNA transcripts, encoded by 3,897 gene loci, were identified, including 2,066 sense lncRNAs, 2,772 lincRNAs, and 1,807 lncNATs. Compared with protein-coding mRNAs, V. destructor lncRNAs are shorter in terms of full length, as well as of the ORF length, contain less exons, and express at lower level. GO term and KEGG pathway enrichment analyses of the lncRNA target genes demonstrated that these predicted lncRNAs may be potentially responsible for the regulatory functions of cellular and biological progresses in the reproductive phase of V. destructor. To our knowledge, this is the first catalog of lncRNA profile in the parasitiformes species, providing a valuable resource for genetic and genomic studies. Understanding the characteristics and features of lncRNAs in V. destructor would promote sustainable parasite control.

Reconstruction and functional annotation of Ascosphaera apis full-length transcriptome utilizing PacBio long reads combined with Illumina short reads

Ascosphaera apis is a widespread fungal pathogen of honeybee larvae that results in chalkbrood disease, leading to heavy losses for the beekeeping industry in China and many other countries. This work was aimed at generating a full-length transcriptome of A. apis using PacBio single-molecule real-time (SMRT) sequencing. Here, more than 23.97 Gb of clean reads was generated from long-read sequencing of A. apis mycelia, including 464,043 circular consensus sequences (CCS) and 394,142 full-length non-chimeric (FLNC) reads. In total, we identified 174,095 high-confidence transcripts covering 5141 known genes with an average length of 2728 bp. We also discovered 2405 genic loci and 11,623 isoforms that have not been annotated yet within the current reference genome. Additionally, 16,049, 10,682, 4520 and 7253 of the discovered transcripts have annotations in the Non-redundant protein (Nr), Clusters of Eukaryotic Orthologous Groups (KOG), Gene Ontology (GO), and Kyoto Encyclopedia of Genes and Genomes (KEGG) databases. Moreover, 1205 long non-coding RNAs (lncRNAs) were identified, which have less exons, shorter exon and intron lengths, shorter transcript lengths, lower GC percent, lower expression levels, and fewer alternative splicing (AS) evens, compared with protein-coding transcripts. A total of 253 members from 17 transcription factor (TF) families were identified from our transcript datasets. Finally, the expression of A. apis isoforms was validated using a molecular approach. Overall, this is the first report of a full-length transcriptome of entomogenous fungi including A. apis. Our data offer a comprehensive set of reference transcripts and hence contributes to improving the genome annotation and transcriptomic study of A. apis.

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.