Circular RNA response of Apisceranacerana 6-day-old larvae to Ascosphaeraapis stress

• 1、College of Animal Sciences(College of Bee Science), Fujian Agriculture and Forestry University, Fuzhou 350002

Abstract： Ascosphaera apis exclusively infects honeybee larvae, leading to chalkbrood, a fungal disease damaging honeybee health and beekeeping industry. The objective of this study was to investigate the differential expression pattern of circular RNAs (circRNAs) and differentially expressed circRNAs (DEcircRNAs) involved in Apis cerana cerana 6-day-old larvae response to Ascosphaera apis stress, and the putative role of host A. apis-response. Normal and A. apis-infected 6-day-old larval guts of A. c. cerana (AcCK and AcT) were sequenced using linear RNA removed circRNA-seq technology. CircRNAs were identified using find_circ software, followed by summary of length and circulization type of circRNAs. DEcircRNAs were screened out following the standard of |log2(Fold change)|≥1 and P≤0.05. Source genes of DEcircRNAs were annotated to Gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) databases to gain function and pathway annotations. RT-qPCR was conducted to validate three randomly selected DEcircRNAs. Based on circRNA-seq, 76342570 and 68269362 raw reads were produced from AcCK and AcT, and 74524108 and 66974392 clean reads were obtained after strict quality control, with Q30 of 92.75% and 94%, and GC content of 54.31% and 54.90%, respectively. A total of 23648400 anchor reads were mapped to the reference genome of Apis cerana. In AcCK and AcT, 805 and 702 circRNAs were respectively identified; the lengths of these circRNAs were distributed among 201nt-1000 nt, and the most abundant circulization type was annotated exonic circRNA; however, the numbers of circRNAs distributed in various lengths and circulization types were different. There were 494 DEcircRNAs in AcCK vs AcT comparison group, including 257 up-regulated circRNAs and 237 down-regulated circRNAs, with novel_circ_000123 and novel_circ_000726 having the highest up-regulation and down-regulation levels. Source genes of the aforementioned DEcircRNAs were annotated to 11 biological process-associated terms, nine molecular function-associated terms and nine cellular component-associated terms. Additionally, these source genes were annotated to 73 pathways such as protein processing in endoplasmic reticulum, spliceosome, and oxidative phosphorylation. Further analysis demonstrated part of source genes were engaged in seven cellular immune pathways including endocytosis, phagocytosis, lysosome, ubiquitin-mediated proteolysis, insect hormone biosynthesis, metabolism of xenobiotics by cytochrome P450, and regulation of autophagy; and three humoral immune pathways including FoxO signaling pathway, MAPK signaling pathway and Jak-STAT signaling pathway. RT-qPCR result indicated the expression trend of three DEcircRNAs expression was in accordance with that in sequencing data. A. c. cerana 6-day-old larvae were likely to response to A. apis stress by altering the numbers of circRNAs distributed in various lengths and circulization types, specifically expressing some circRNAs, and regulating the expression of partial circRNAs; several DEcircRNAs such as novel_circ_000027, novel_circ_000127 and novel_circ_000312 may play a special role in host A. apis-response via regulation of oxidative phosphorylation, cellular and humoral immune pathways. Findings in this work offered novel insights into further understanding molecular mechanisms underlying A. c. cerana larval response to A. apis as well as host-pathogen interaction.

Keywords： Apis cerana cerana larva circular RNA non-coding RNA Ascosphaera apis immune response