CircRNA-regulated immune responses of asian honey bee workers to microsporidian infection

Unraveling the modulatory manner and function of circRNAs in the Asian honey bee larval guts

Circular RNAs (circRNAs) are a class of non-coding RNAs (ncRNAs) that can participate in biological processes such as gene expression, growth, and development. However, little has been explored about the function of circRNAs in the development of Apis cerana larval guts. By using our previously gained deep sequencing data from the guts of A. cerana worker larvae at 4-, 5-, and 6-day-old (Ac4, Ac5, and Ac6 groups), the expression pattern and regulatory role of circular RNAs (circRNAs) during the development process was comprehensively investigated, with a focus on differentially expressed circRNAs (DEcircRNAs) relevant to immunity pathways and developmental signaling pathways, followed by validation of the binding relationships among a key competing endogenous RNA (ceRNA) axis. Here, 224 (158) DEcircRNAs were detected in the Ac4 vs. Ac5 (Ac5 vs. Ac6) comparison group. It's suggested that 172 (123) parental genes of DEcircRNAs were involved in 26 (20) GO terms such as developmental process and metabolic process and 138 (136) KEGG pathways like Hippo and Wnt signaling pathways. Additionally, ceRNA network analysis indicated that 21 (11) DEcircRNAs could target seven (three) DEmiRNAs, further targeting 324 (198) DEmRNAs. These DEmRNAs can be annotated to 33 (26) GO terms and 168 (200) KEGG pathways, including 12 (16) cellular and humoral immune pathways (endocytosis, lysosome, Jak-STAT, etc.) and 10 (nine) developmental signaling pathways (Hippo, mTOR, Hedgehog, etc.). Interestingly, DEcircRNAs in these two comparison groups could target the same ace-miR-6001-y, forming complex sub-networks. The results of PCR and Sanger sequencing confirmed the back-splicing sites within four randomly selected DEcircRNAs. RT-qPCR detection of these four DEcircRNAs verified the reliability of the used transcriptome data. The results of dual-luciferase reporter assay verified the binding relationships between novel_circ_001627 and ace-miR-6001-y and between ace-miR-6001-y and apterous-like. Our data demonstrated that DEcircRNAs were likely to modulate the developmental process of the A. cerana worker larval guts via regulation of parental gene transcription and ceRNA network, and novel_circ_001627/ace-miR-6001-y/apterous-like was a potential regulatory axis in the larval gut development. Findings from this work offer a basis and a candidate ceRNA axis for illustrating the circRNA-modulated mechanisms underlying the A. cerana larval guts.

A systemic study of cyenopyrafen in strawberry cultivation system: Efficacy, residue behavior, and impact on honeybees (Apis mellifera L.)

The pesticide application method is one of the important factors affecting its effectiveness and residues, and the risk of pesticides to non-target organisms. To elucidate the effect of application methods on the efficacy and residue of cyenopyrafen, and the toxic effects on pollinators honeybees in strawberry cultivation, the efficacy and residual behavior of cyenopyrafen were investigated using foliar spray and backward leaf spray in field trials. The results showed that the initial deposition of cyenopyrafen using backward leaf spray on target leaves reached 5.06-9.81 mg/kg at the dose of 67.5-101.25 g a.i./ha, which was higher than that using foliar spray (2.62-3.71 mg/kg). The half-lives of cyenopyrafen in leaves for foliar and backward leaf spray was 2.3-3.3 and 5.3-5.9 d, respectively. The residues (10 d) of cyenopyrafen in leaves after backward leaf spray was 1.41-3.02 mg/kg, which was higher than that after foliar spraying (0.25-0.37 mg/kg). It is the main reason for the better efficacy after backward leaf spray. However, the residues (10 d) in strawberry after backward leaf spray and foliar spray was 0.04-0.10 and < 0.01 mg/kg, which were well below the established maximum residue levels of cyenopyrafen in Japan and South Korea for food safety. To further investigate the effects of cyenopyrafen residues after backward leaf spray application on pollinator honeybees, sublethal effects of cyenopyrafen on honeybees were studied. The results indicated a significant inhibition in the detoxification metabolic enzymes of honeybees under continuous exposure of cyenopyrafen (0.54 and 5.4 mg/L) over 8 d. The cyenopyrafen exposure also alters the composition of honeybee gut microbiota, such as increasing the relative abundance of Rhizobiales and decreasing the relative abundance of Acetobacterales. The comprehensive data on cyenopyrafen provide basic theoretical for environmental and ecological risk assessment, while backward leaf spray proved to be effective and safe for strawberry cultivation.

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.

Novel Insights into the circRNA-Modulated Developmental Mechanism of Western Honey Bee Larval Guts

Circular RNAs (circRNAs) are a class of novel non-coding RNAs (ncRNAs) that play essential roles in the development and growth of vertebrates through multiple manners. However, the mechanism by which circRNAs modulate the honey bee gut development is currently poorly understood. Utilizing the transcriptome data we obtained earlier, the highly expressed circRNAs in the Apis mellifera worker 4-, 5-, and 6-day-old larval guts were analyzed, which was followed by an in-depth investigation of the expression pattern of circRNAs during the process of larval guts development and the potential regulatory roles of differentially expressed circRNAs (DEcircRNAs). In total, 1728 expressed circRNAs were detected in the A. mellifera larval guts. Among the most highly expressed 10 circRNAs, seven (novel_circ_000069, novel_circ_000027, novel_circ_000438, etc.) were shared by the 4-, 5-, and 6-day-old larval guts. In addition, 21 (46) up-regulated and 22 (27) down-regulated circRNAs were, respectively, screened in the Am4 vs. Am5 (Am5 vs. Am6) comparison groups. Additionally, nine DEcircRNAs, such as novel_circ_000340, novel_circ_000758 and novel_circ_001116, were shared by these two comparison groups. These DEcircRNAs were predicted to be transcribed from 14 and 29 parental genes; these were respectively annotated to 15 and 22 GO terms such as biological regulation and catalytic activity as well as 16 and 21 KEGG pathways such as dorsoventral axis formation and apoptosis. Moreover, a complicated competing endogenous RNA (ceRNA) network was observed; novel_circ_000838 in the Am4 vs. Am5 comparison group potentially targeted ame-miR-6000a-3p, further targeting 518 mRNAs engaged in several developmental signaling pathways (e.g., TGF-beta, hedgehog, and wnt signaling pathway) and immune pathways (e.g., phagosome, lysosome, and MAPK signaling pathway). The results demonstrated that the novel_circ_000838-ame-miR-6000a-3p axis may plays a critical regulatory part in the larval gut development and immunity. Furthermore, back-splicing sites of six randomly selected DEcircRNAs were amplified and verified by PCR; an RT-qPCR assay of these six DEcircRNAs confirmed the reliability of the used high-throughput sequencing data. Our findings provide a novel insight into the honey bee gut development and pave a way for illustration of the circRNA-modulated developmental mechanisms underlying the A. mellifera worker larval guts.

A comparative transcriptome analysis of how shrimp endure and adapt to long-term symbiosis with Enterocytozoon hepatopenaei infection

Enterocytozoon hepatopenaei (EHP) is a prevalent microsporidian pathogen responsible for hepatopancreatic microsporidiosis (HPM) in Litopenaeus vannamei. This infection not only leads to slowed growth in shrimp abut aslo inflicts substantial economic losses in the global aquaculture industry. However, the molecular mechanisms by which EHP influences the host during various infection stages remain unclear. This study employed comparative transcriptomics to examine the effects of EHP infection on Litopenaeus vannamei between early and late stage of infection groups. Utilizing transcriptomic approaches, we identified differentially expressed genes (DEGs) with notable biological significance through the COG, GO, KEGG, GSEA, and Mufzz time-series methodologies. The results reveal that EHP infection considerably influences host gene expression, with marked differences between early and late infection across distinct timeframes. Key processes such as detoxification, cell apoptosis, and lipid metabolism are pivotal during host-parasite interactions. Hexokinase and phosphatidic acid phosphatase emerge as key factors enabling invasion and sustained effects. Cytochrome P450 and glucose-6-phosphate dehydrogenase could facilitate infection progression. EHP significantly impacts growth, especially through ecdysteroids and 17β-estradiol dehydrogenase. By delineating stage-specific effects, we gain insights into interaction between EHP and Litopenaeus vannamei, showing how intracellular pathogens reprogram host defenses into mechanisms enabling long-term persistence. This study provides a deeper understanding of host-pathogen dynamics, emphasizing the interplay between detoxification, metabolism, immunity, apoptosis and growth regulation over the course of long-term symbiosis.

The Potential Biological Roles of Circular RNAs in the Immune Systems of Insects to Pathogen Invasion

Circular RNAs (circRNAs) are a newly discovered class of endogenously expressed non-coding RNAs (ncRNAs). They are highly stable, covalently closed molecules that frequently exhibit tissue-specific expression in eukaryotes. A small number of circRNAs are abundant and have been remarkably conserved throughout evolution. Numerous circRNAs are known to play important biological roles by acting as microRNAs (miRNAs) or protein inhibitors ('sponges'), by regulating the function of proteins, or by being translated themselves. CircRNAs have distinct cellular functions due to structural and production differences from mRNAs. Recent advances highlight the importance of characterizing circRNAs and their targets in a variety of insect species in order to fully understand how they contribute to the immune responses of these insects. Here, we focus on the recent advances in our understanding of the biogenesis of circRNAs, regulation of their abundance, and biological roles, such as serving as templates for translation and in the regulation of signaling pathways. We also discuss the emerging roles of circRNAs in regulating immune responses to various microbial pathogens. Furthermore, we describe the functions of circRNAs encoded by microbial pathogens that play in their hosts.

Systematic Characterization and Regulatory Role of lncRNAs in Asian Honey Bees Responding to Microsporidian Infestation

Long noncoding RNAs (lncRNAs) are pivotal regulators in gene expression and diverse biological processes, such as immune defense and host-pathogen interactions. However, little is known about the roles of lncRNAs in the response of the Asian honey bee (Apis cerana) to microsporidian infestation. Based on our previously obtained high-quality transcriptome datasets from the midgut tissues of Apis cerana cerana workers at 7 days post inoculation (dpi) and 10 dpi with Nosema ceranae (AcT7 and AcT10 groups) and the corresponding un-inoculated midgut tissues (AcCK7 and AcCK10 groups), the transcriptome-wide identification and structural characterization of lncRNAs were conducted, and the differential expression pattern of lncRNAs was then analyzed, followed by investigation of the regulatory roles of differentially expressed lncRNAs (DElncRNAs) in host response. Here, 2365, 2322, 2487, and 1986 lncRNAs were, respectively, identified in the AcCK7, AcT7, AcCK7, and AcT10 groups. After removing redundant ones, a total of 3496 A. c. cerana lncRNAs were identified, which shared similar structural characteristics with those discovered in other animals and plants, such as shorter exons and introns than mRNAs. Additionally, 79 and 73 DElncRNAs were screened from the workers' midguts at 7 dpi and 10 dpi, respectively, indicating the alteration of the overall expression pattern of lncRNAs in host midguts after N. ceranae infestation. These DElncRNAs could, respectively, regulate 87 and 73 upstream and downstream genes, involving a suite of functional terms and pathways, such as metabolic process and Hippo signaling pathway. Additionally, 235 and 209 genes co-expressed with DElncRNAs were found to enrich in 29 and 27 terms, as well as 112 and 123 pathways, such as ABC transporters and the cAMP signaling pathway. Further, it was detected that 79 (73) DElncRNAs in the host midguts at 7 (10) dpi could target 321 (313) DEmiRNAs and further target 3631 (3130) DEmRNAs. TCONS_00024312 and XR_001765805.1 were potential precursors for ame-miR-315 and ame-miR-927, while TCONS_00006120 was the putative precursor for both ame-miR-87-1 and ame-miR-87-2. These results together suggested that DElncRNAs are likely to play regulatory roles in the host response to N. ceranae infestation through the regulation of neighboring genes via a cis-acting effect, modulation of co-expressed mRNAs via trans-acting effect, and control of downstream target genes' expression via competing endogenous RNA networks. Our findings provide a basis for disclosing the mechanism underlying DElncRNA-mediated host N. ceranae response and a new perspective into the interaction between A. c. cerana and N. ceranae.

Comparative transcriptome analysis of adult worker bees under short-term heat stress

High temperature affects behavior, physiology, survival, and the expression of related genes in adult honeybees. Apis mellifera is the common pollinator in greenhouse and is susceptible to high temperature stress. To further explore the molecular basis related to heat stress, we compared the transcriptome profiles of adult worker bees at 25 and 45°C, and detected the expression patterns of some differentially expressed genes (DEGs) in different tissues by q RT-PCR. Differential expression analysis showed that 277 DEGs were identified, including 167 genes upregulated and 110 genes downregulated after heat stress exposure in adult worker bees. In GO enrichment analysis, DEGs were mostly enriched for protein folding, unfold protein binding, and heme binding terms. Protein processing in endoplasmic reticulum and longevity regulating pathway-multiple species were significantly enriched in KEGG. The expression levels of 16 DEGs were consistent with the transcriptome results. The expression patterns of 9 DEGs in different tissues revealed high levels in the thorax, which was supposed that the thorax may be the most important part in the response to heat stress. This study provided valuable data for exploring the function of heat resistance-related genes.

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.

Current Knowledge on Bee Innate Immunity Based on Genomics and Transcriptomics

As important pollinators, bees play a critical role in maintaining the balance of the ecosystem and improving the yield and quality of crops. However, in recent years, the bee population has significantly declined due to various pathogens and environmental stressors including viruses, bacteria, parasites, and increased pesticide application. The above threats trigger or suppress the innate immunity of bees, their only immune defense system, which is essential to maintaining individual health and that of the colony. In addition, bees can be divided into solitary and eusocial bees based on their life traits, and eusocial bees possess special social immunities, such as grooming behavior, which cooperate with innate immunity to maintain the health of the colony. The omics approach gives us an opportunity to recognize the distinctive innate immunity of bees. In this regard, we summarize innate bee immunity from a genomic and transcriptomic perspective. The genetic characteristics of innate immunity were revealed by the multiple genomes of bees with different kinds of sociality, including honeybees, bumblebees, wasps, leaf-cutter bees, and so on. Further substantial transcriptomic data of different tissues from diverse bees directly present the activation or suppression of immune genes under the infestation of pathogens or toxicity of pesticides.