Piwi-interacting RNAs (piRNAs), a class of small non-coding RNAs (ncRNAs), play pivotal roles in maintaining the genomic stability and modulating biological processes such as growth and development via the regulation of gene expression. However, the piRNAs in the Asian honeybee (Apis cerana) are still largely unknown at present. In this current work, on the basis of previously gained high-quality small RNA-seq datasets, piRNAs in the larval gut of Apis cerana cerana, the nominated species of A. cerana, were identified for the first time, followed by an in-depth investigation of the regulatory roles of differentially expressed piRNAs (DEpiRNAs) in the developmental process of the A. c. cerana. Here, a total of 621 piRNAs were identified in A. c. cerana larval guts, among which 499 piRNAs were shared by 4-(Ac4 group), 5-(Ac5 group), and 6-day-old (Ac6 group) larval guts, while the numbers of unique ones equaled 79, 37, and 11, respectively. The piRNAs in each group ranged from 24 nucleotides (nt) to 33 nt in length, and the first base of the piRNAs had a cytosine (C) bias. Additionally, five up-regulated and five down-regulated piRNAs were identified in the Ac4 vs. Ac5 comparison group, nine of which could target 9011 mRNAs; these targets were involved in 41 GO terms and 137 pathways. Comparatively, 22 up-regulated piRNAs were detected in the Ac5 vs. Ac6 comparison group, 21 of which could target 28,969 mRNAs; these targets were engaged in 46 functional terms and 164 pathways. The results suggested an overall alteration of the expression pattern of piRNAs during the developmental process of A. c. cerana larvae. The regulatory network analysis showed that piR-ace-748815 and piR-ace-512574 in the Ac4 vs. Ac5 comparison group as well as piR-ace-716466 and piR-ace-828146 in the Ac5 vs. Ac6 comparison group were linked to the highest number of targets. Further investigation indicated that targets of DEpiRNAs in the abovementioned two comparison groups could be annotated to several growth and development-associated pathways, such as the Jak/STAT, TGF-β, and Wnt signaling pathways, indicating the involvement of DEpiRNAs in modulating larval gut development via these crucial pathways. Moreover, the expression trends of six randomly selected DEpiRNAs were verified using a combination of stem-loop RT-PCR and RT-qPCR. These results not only provide a novel insight into the development of the A. c. cerana larval gut, but also lay a foundation for uncovering the epigenetic mechanism underlying larval gut development.

Apis cerana is the original host for Nosema ceranae, a widespread fungal parasite resulting in honey bee nosemosis, which leads to severe losses to the apiculture industry throughout the world. However, knowledge of N. ceranae infecting eastern honey bees is extremely limited. Currently, the mechanism underlying N. ceranae infection is still largely unknown. Based on our previously gained high-quality transcriptome datasets derived from N. ceranae spores (NcCK group), N. ceranae infecting Apis cerana cerana workers at seven days post inoculation (dpi) and 10 dpi (NcT1 and NcT2 groups), comparative transcriptomic investigation was conducted in this work, with a focus on virulence factor-associated differentially expressed genes (DEGs). Microscopic observation showed that the midguts of A. c. cerana workers were effectively infected after inoculation with clean spores of N. ceranae. In total, 1411, 604, and 38 DEGs were identified from NcCK vs. NcT1, NcCK vs. NcT2, and NcT1 vs. NcT2 comparison groups. Venn analysis showed that 10 upregulated genes and nine downregulated ones were shared by the aforementioned comparison groups. The GO category indicated that these DEGs were involved in a series of functional terms relevant to biological process, cellular component, and molecular function such as metabolic process, cell part, and catalytic activity. Additionally, KEGG pathway analysis suggested that the DEGs were engaged in an array of pathways of great importance such as metabolic pathway, glycolysis, and the biosynthesis of secondary metabolites. Furthermore, expression clustering analysis demonstrated that the majority of genes encoding virulence factors such as ricin B lectins and polar tube proteins displayed apparent upregulation, whereas a few virulence factor-associated genes such as hexokinase gene and 6-phosphofructokinase gene presented downregulation during the fungal infection. Finally, the expression trend of 14 DEGs was confirmed by RT-qPCR, validating the reliability of our transcriptome datasets. These results together demonstrated that an overall alteration of the transcriptome of N. ceranae occurred during the infection of A. c. cerana workers, and most of the virulence factor-related genes were induced to activation to promote the fungal invasion. Our findings not only lay a foundation for clarifying the molecular mechanism underlying N. ceranae infection of eastern honey bee workers and microsporidian–host interaction.

Here, a comparative transcriptome investigation was conducted based on high-quality deep sequencing data from the midguts of Apis cerana cerana workers at 7 d post-inoculation (dpi) and 10 dpi with Nosema ceranae and corresponding un-inoculated midguts. PCR identification and microscopic observation of paraffin sections confirmed the effective infection of A. c. cerana worker by N. ceranae. In total, 1127 and 957 N. ceranae-responsive genes were identified in the infected midguts at 7 dpi and 10 dpi, respectively. RT-qPCR results validated the reliability of our transcriptome data. GO categorization indicated the differentially expressed genes (DEGs) were respectively engaged in 34 and 33 functional terms associated with biological processes, cellular components, and molecular functions. Additionally, KEGG pathway enrichment analysis showed that DEGs at 7 dpi and 10 dpi could be enriched in 231 and 226 pathways, respectively. Moreover, DEGs in workers’ midguts at both 7 dpi and 10 dpi were involved in six cellular immune pathways such as autophagy and phagosome and three humoral immune pathways such as the Toll/Imd signaling pathway and Jak-STAT signaling pathway. In addition, one up-regulated gene (XM_017055397.1) was enriched in the NF-κB signaling pathway in the workers’ midgut at 10 dpi. Further investigation suggested the majority of these DEGs were engaged in only one immune pathway, while a small number of DEGs were simultaneously involved in two immune pathways. These results together demonstrated that the overall gene expression profile in host midgut was altered by N. ceranae infection and some of the host immune pathways were induced to activation during fungal infection, whereas some others were suppressed via host–pathogen interaction. Our findings offer a basis for clarification of the mechanism underlying the immune response of A. c. cerana workers to N. ceranae infection, but also provide novel insights into eastern honeybee-microsporodian interaction.

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.

Aspergillus flavus and its secondary metabolites aflatoxins pose a significant threat to the health of humans, animals, and plants. Therefore, there is an urgent need to control A. flavus contamination. AfverB plays a key role in the aflatoxin gene cluster; however, its function and mechanism in fungal development and virulence remain poorly understood. In this study, we constructed afVerB gene deletion mutants (∆afVerB−1 and ∆afVerB−2) and two CYP domain mutants (afVerB^∆D1 and afVerB^∆D2) through homologous recombination. Phenotype analysis revealed that, via its two CYP domains, AfVerB is deeply involved in fungal morphogenesis and aflatoxin synthesis. Insect and crop colonization models revealed that AfVerB plays a key role in the fungus’s ability to infect hosts, and stress experiments discovered that AfVerB plays a significant role in the response to various environmental stresses, which explains why AfVerB is a key factor in fungal infection to some extent. RT-qPCR analysis demonstrated that AfVerB performs its bio-function through corresponding regulatory factors. We ultimately discovered that AfVerB is deeply involved in cell membrane stress stability, thereby participating in the regulation of fungal drug resistance (sensitive to AMB and resistant to VOR in this study). The CYP domain of AfVerB, particularly its second CYP domain, is crucial for the execution of its biological functions. This study elucidated the regulatory mechanisms by which AfVerB regulates fungal pathogenicity and aflatoxin biosynthesis, providing potential strategies for controlling A. flavus and its aflatoxin contamination.

Aspergillus flavus Link, 1809 is a pathogenic fungus widely present in the environment. It can infect plants and also acts as an opportunistic pathogen affecting humans and other animals. The aflatoxins (AFs), it produces, can cause cancers such as liver cancer. Therefore, in-depth research into the pathogenic mechanisms of A. flavus is crucial. Arp8 (Actin-like protein Arp8) is a unique subunit within the chromatin remodeling complex INO80, regulating processes including chromatin remodeling. However, the biological function of Arp8 in A. flavus remains unclear. This study constructed A. flavus arp8 knockout (Δarp8) and complementation (Com-arp8) strains via homologous recombination. Subsequent research revealed that following the deletion of arp8, A. flavus exhibits a reduction of approximately 51% in conidia production, complete abrogation of sclerotia formation, and significantly impairment of aflatoxin B1 (AFB1) biosynthesis. Crop grain colonization and Bombyx mori Linnaeus, 1758 infection models demonstrated that Arp8 plays a crucial role in A. flavus ability to infect hosts. Environmental stress experiments identified Arp8 as a vital factor for A. flavus in response to various environmental stresses. Quantitative RT-PCR (qRT-PCR) analysis indicated Arp8 achieves its biological functions through corresponding regulatory factors. This study elucidates the biological functions of Arp8 in A. flavus growth and development, pathogenicity, and aflatoxin synthesis, laying a foundation to illuminate the mechanisms of A. flavus pathogenicity and AFs production.