为探究中华蜜蜂（Apis cerana cerana）与意大利蜜蜂（Apis mellifera ligustica）幼虫的球囊菌抗性差异产生的原因，利用Venn分析、GO分类分析、KEGG代谢通路分析以及Ka/ks分析对二者的转录组进行比较研究。Venn分析结果显示，中蜂与意蜂的同源基因有17 656 个，归属于8 111 个基因家族，二者的特有基因分别有1 158和241 个，分别归属于458和86 个基因家族。GO分类结果显示，中蜂和意蜂的特有基因分别富集在38和28 个GO条目。KEGG代谢通路富集分析结果显示，中蜂和意蜂的特有基因分别富集于96和21 个代谢通路；进一步分析发现中蜂的特有基因富集在内吞作用、溶酶体、泛素介导的蛋白水解和黑化作用等4 个细胞免疫通路，以及MAPK信号通路和Toll-like受体信号通路等2 个体液免疫通路，而意蜂仅有1 个特有基因富集在MAPK信号通路。Ka/Ks分析结果显示受到强烈正向选择、弱正向选择、中性选择和负向选择的单拷贝同源基因分别有37、281、2 630和3 269 个。对受到强烈正向选择的基因进行GO分类和KEGG代谢通路富集分析，GO结果显示中蜂和意蜂的单拷贝同源基因富集的GO 条目类型相同；KEGG结果显示中蜂的单拷贝基因仅富集在氧化磷酸化。上述结果表明免疫相关基因数量的差异是二者的球囊菌抗性差异产生的重要原因之一，中蜂在与球囊菌的协同进化过程中可能通过提高能量利用从而限制球囊菌的增殖。本研究结果可为阐明中蜂及意蜂幼虫的球囊菌抗性差异产生的分子机制提供参考。

Bombyx mori cypovirus (BmCPV) is one of the major viral pathogen for silkworm, and the genome of BmCPV is composed of 10 dsRNA segments. As construction system of recombinant BmCPV (rBmCPV) is scanty, researchers achieved little progress in studying gene function of BmCPV in recent decades. Here, 10 recombinant plasmids with a full-length cDNA of viral genome segments S1-S10 containing T7 promoter were constructed. After cotransfecting the BmN cells with the mixture of 10 in vitrotranscribed RNAs, pathological changes were observed. Real-time PCR and Western blot showed viral gene vp1 and structural proteins were expressed. It is found the genome of the rBmCPV is composed of 10 dsRNA segments similar to those of wild-type BmCPV. Moreover, viral particles and polyhedron with virions can be generated in the cotransfected cells and the injected silkworm midguts. These findings confirmed the formation of infective rBmCPV. Additionally, we found viable rBmCPV was generated by cotransfecting the mixture of in vitro-transcribed S1-S9 RNAs into the cultured cells, confirming polh was not essential for BmCPV replication. Moreover, an infectious rBmCPV expressing the DsRed protein was constructed based on this system. Further investigation showed S2 and S7 segments are indispensible for viral proliferation. Our findings demonstrated the construction system of rBmCPV can be utilized for exploring viral replication and pathogenesis, and investigated method for constructing rBmCPV will certainly facilitate developing novel biopesticides and expressing exogenous gene in the midgut of silkworm.

Nosema ceranae is a unicellular fungal parasite of honey bees and causes huge losses for apiculture. Until present, no study on N . ceranae long non-coding RNAs (lncRNAs) was documented. Here, we sequenced purified spores of N . ceranae using strand-specific library construction and high-throughput RNA sequencing technologies. In total, 83 novel lncRNAs were predicted from N . ceranae spore samples, including lncRNAs, long intergenic non-coding RNAs (lincRNAs), and sense lncRNAs. Moreover, these lncRNAs share similar characteristics with those identified in mammals and plants, such as shorter length and fewer exon number and transcript isoforms than protein-coding genes. Finally, the expression of 12 lncRNAs was confirmed with RT-PCR, confirming their true existence. To our knowledge, this is the first evidence of lncRNAs produced by a microsporidia species, offering novel insights into basic biology such as regulation of gene expression of this widespread taxonomic group.

Currently, knowledge of circular RNAs (circRNAs) in insects including honeybee is extremely limited. Here, differential expression profiles and regulatory networks of circRNAs in the midguts of Apis cerana cerana workers were comprehensively investigated using transcriptome sequencing and bioinformatics. In total, 9589 circRNAs (201–800 nt in length) were identified from 8-day-old and 11-day-old workers’ midguts (Ac1 and Ac2); among them, 5916 (61.70%) A. cerana cerana circRNAs showed conservation with our previously indentified circRNAs in Apis mellifera ligucstica workers’ midguts (Xiong et al., Acta Entomologica Sinica 61:1363–1375, 2018). Five circRNAs were confirmed by RT-PCR and Sanger sequencing. Interestingly, novel_circ_003723, novel_circ_002714, novel_circ_002451, and novel_circ_001980 were highly expressed in both Ac1 and Ac2. In addition, the source genes of circRNAs were involved in 34 GO terms including organelle and cellular process and 141 pathways such as endocytosis and Wnt signaling pathway. Moreover, 55 DEcircRNAs including 34 upregulated and 21 downregulated circRNAs were identified in Ac2 compared with Ac1. circRNA-miRNA regulatory networks indicated that 1060 circRNAs can target 74 miRNAs; additionally, the DEcircRNA-miRNA-mRNA networks suggested that 13 downregulated circRNAs can bind to eight miRNAs and 29 miRNA-targeted mRNAs, while 16 upregulated circRNAs can link to 9 miRNAs and 29 miRNA-targeted mRNAs. These results indicated that DEcircRNAs as ceRNAs may play a comprehensive role in the growth, development, and metabolism of the worker’s midgut via regulating source genes and interacting with miRNAs. Notably, eight DEcircRNAs targeting miR-6001-y were likely to be key participants in the midgut development. Our findings not only offer a valuable resource for further studies on A. cerana cerana circRNA and novel insights into understanding the molecular mechanisms underlying the midgut development of eastern honeybee but also provide put

Ascosphaera apis is a widespread fungal pathogen of honeybee larvae, which causes heavy losses in apiculture. To date, knowledge about non-coding RNA (ncRNA) including circular RNA (circRNA) in A. apis is lacking. In this study, A. apis mycelia and spores were sequenced using RNA-seq technology. A total of 551 circRNAs were predicted on the basis of bioinformatic analyses, and most of the circRNAs were 200–600 bp in length, which were different from animal and plant circRNAs. In addition, the expression of six circRNAs in A. apis were confirmed using divergent and convergent PCR. Moreover, circRNA-microRNA regulation networks in A. apis were constructed, and further investigation showed that A. apis circRNAs could regulate gene expression by competitively binding miRNAs. GO and KEGG pathway enrichment analyses of the miRNAs target genes of circRNAs demonstrated that these A. apis circRNAs are likely to play key roles in metabolism, environmental response and gene expression.

Circular RNAs (circRNAs) are newly discovered endogenous non-coding RNAs (ncRNAs) that play key roles in microRNA function and transcriptional regulation. Though a large number of circRNAs had been identifed in animals and plants, however, little is known regarding circRNAs in Nosema ceranae, a widespread fungal parasite of honeybee. In this study, using deep sequencing technology and bioinformatic analysis, we predicted 204 circRNAs from N. ceranae spore samples, including 174 exonic circRNAs and 30 intergenic circRNAs. In addition, the expression of seven N. ceranae circRNAs was confrmed by RT-PCR assay. Furthermore, regulation networks of circRNAs were constructed, and 15 circRNAs were found to act as sponges of the corresponding three miRNAs. GO categorization and pathway enrichment analysis suggested that the circRNAs are likely to play signifcant roles in N. ceranae spore. This is the frst report of circRNAs generated by a microsporidia species. Our results provide novel insights into understanding the basic biology of N. ceranae.

Honeybees are susceptible to a variety of diseases, including chalkbrood, which is capable of causing huge losses of both the number of bees and colony productivity. This research is designed to characterize the transcriptome profiles of Ascosphaera apis-treated and un-treated larval guts of Apis mellifera ligustica in an attempt to unravel the molecular mechanism underlying the immune responses of western honeybee larval guts to mycosis. In this study, 24, 296 and 2157 genes were observed to be differentially expressed in A. apis-treated Apis mellifera (4-, 5- and 6-day-old) compared with un-treated larval guts. Moreover, the expression patterns of differentially expressed genes (DEGs) were examined via trend analysis, and subsequently, gene ontology analysis and KEGG pathway enrichment analysis were conducted for DEGs involved in up- and down-regulated profiles. Immunityrelated pathways were selected for further analysis, and our results demonstrated that a total of 13 and 50 DEGs were annotated in the humoral immune-related and cellular immune-related pathways, respectively. Additionally, we observed that many DEGs up-regulated in treated guts were part of cellular immune pathways, such as the lysosome, ubiquitin mediated proteolysis, and insect hormone biosynthesis pathways and were induced by A. apis invasion. However, more down-regulated DEGs were restrained. Surprisingly, a majority of DEGs within the Toll-like receptor signaling pathway, and the MAPK signaling pathway were up-regulated in treated guts, while all but two genes involved in the NF-κB signaling pathway were down-regulated, which suggested that most genes involved in humoral immune-related pathways were activated in response to the invasive fungal pathogen. This study's findings provide valuable information regarding the investigation of the molecular mechanism of immunity defenses of A. m. ligustica larval guts to infection with A. apis.

Long non-coding RNAs (lncRNAs) have been proven to play pivotal roles in a wide range of biological processes in animals and plants. Here, we performed whole transcriptome strand-specific RNA sequencing of normal midguts of Apis mellifera ligustica workers (Am7CK, Am10CK) and Nosema ceranae-inoculated midguts (Am7T, Am10T), and 4749 conserved lncRNAs and 1604 novel lncRNAs were identified. 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 offer a rich genetic resource for further investigation of the functional roles of lncRNAs involved in the A. m. ligustica response to N. ceranae infection.

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 downregulated 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 honeybeemicrosporidian interactions.

Chalkbrood is the most common fungal disease in honeybees. The objective of this study was to reveal immune responses in the Apis cerana cerana larval gut following Ascosphaera apis invasion. Combining a previously assembled transcriptome of A. c. cerana larval gut and the high-throughput sequencing data obtained in this study, 6152 differentially expressed genes (DEGs) were clustered into eight profiles. Trend analysis showed three significant up-regulated profiles (p ≤ 0.05) and three down-regulated profiles. Gene Ontology (GO) term analysis suggested that DEGs within significant up-regulated and down-regulated clusters were enriched in 46 and 38 functional groups, respectively. Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis indicated a majority of DEGs were involved in ribosome structure or function, carbon metabolism, biosynthesis of amino acids, and oxidative phosphorylation. In addition, 142 and 14 DEGs were annotated in the cellular immune- and humoral immune-related pathways, respectively. Further investigation indicated that DEGs up-regulated in cellular immune and humoral immune pathways outnumbered those that were downregulated. Moreover, immune responses of A. c. cerana and Apis mellifera ligustica larvae were compared and studied to decipher resistance of eastern honeybee larvae to A. apis. These results demonstrated that a large number of genes involved in immunity-related pathways were activated by A. apis. Our findings provided valuable information for elucidating the molecular mechanisms underlying immune responses of A. c. cerana larvae to A. apis infection and pathogen-host interactions during chalkbrood infection.