Transcriptome of Nosema ceranae and Upregulated Microsporidia Genes during Its Infection of Western Honey Bee (Apis mellifera).

Nosema ceranae is one of the fungal parasites of Apis mellifera. It causes physical and behavioral effects in honey bees. However, only a few studies have reported on gene expression profiling during A. mellifera infection. In this study, the transcriptome profile of mature spores at each time point of infection (5, 10, and 20 days post-infection, d.p.i.) were investigated. Based on the transcriptome and expression profile analysis, a total of 878, 952, and 981 differentially expressed genes (DEGs) (fold change ≥ 2 or ≤ -2) were identified in N. ceranae spores (NcSp) at 5 d.p.i., 10 d.p.i., and 20 d.p.i., respectively. Moreover, 70 upregulated genes and 340 downregulated genes among common DEGs (so-called common DEGs) and 166 stage-specific genes at each stage of infection were identified. The Gene Ontology (GO) analysis indicated that the DEGs and corresponding common DEGs are involved in the functions of cytosol (GO:0005829), cytoplasm (GO:0005737), and ATP binding (GO:0005524). Furthermore, the pathway analysis found that the DEGs and common DEGs are involved in metabolism, environmental information processing, and organismal systems. Four upregulated common DEGs with higher fold-change values, highly associated with spore proteins and transcription factors, were selected for validation. In addition, the stage-specific genes are highly involved in the mechanism of pre-mRNA splicing according to GO enrichment analysis; thus, three of them showed high expression at each d.p.i. and were also subjected to validation. The relative gene expression levels showed a similar tendency as the transcriptome predictions at different d.p.i., revealing that the gene expression of N. ceranae during infection may be related to the mechanism of gene transcription, protein synthesis, and structural proteins. Our data suggest that the gene expression profiling of N. ceranae at the transcriptomic level could be a reference for the monitoring of nosemosis at the genetic level.

Sacbrood viruses cross-infection between Apis cerana and Apis mellifera: Rapid detection, viral dynamics, evolution and spillover risk assessment.

Recent outbreaks of sacbrood virus (SBV) have caused serious epizootic disease in Apis cerana populations across Asia including Taiwan. Earlier phylogenetic analyses showed that cross-infection of AcSBV and AmSBV in both A. cerana and A. mellifera seems common, raising a concern of cross-infection intensifying the risk of disease resurgence in A. cerana. In this study, we analyzed the dynamics of cross-infection in three different types of apiaries (A. mellifera-only, A. cerana-only and two species co-cultured apiaries) over one year in Taiwan. Using novel, genotype-specific primer sets, we showed that SBV infection status varies across apiaries: AmSBV-AM and AcSBV-AC were the major genotype in the A. mellifera-only and the A. cerana-only apiaries, respectively, while AmSBV-AC and AcSBV-AC were the dominant genotypes in the co-cultured apiaries. Interestingly, co-cultured apiaries were among the only apiary type that harbored all variants and dual infections (i.e., AC and AM genotype co-infection in a single sample), indicating the interactions between hosts may form a conduit for cross-infection. The cross-infection between the two honey bee species appears to occur in a regular cycle with temporal fluctuation of AmSBV-AC and AcSBV-AC prevalence synchronized to each other in the co-cultured apiaries. Artificial infection of AcSBV in A. mellifera workers showed the suppression of viral replication, suggesting the potential of A. mellifera serving as a AcSBV reservoir that may contribute to virus spillover. Furthermore, the survival rate of A. cerana larvae was significantly reduced after artificial infections of both SBVs, indicating fitness costs of cross-infection on A. cerana and thus a high risk of disease resurgence in co-cultured apiaries. Our field and laboratory data provide baseline information that facilitates understanding of the risk of SBV cross-infection, and highlights the urgent need of SBV monitoring in co-cultured apiaries.

Transcriptome-level assessment of the impact of deformed wing virus on honey bee larvae.

Deformed wing virus (DWV) prevalence is high in honey bee (Apis mellifera) populations. The virus infects honey bees through vertical and horizontal transmission, leading to behavioural changes, wing deformity, and early mortality. To better understand the impacts of viral infection in the larval stage of honey bees, artificially reared honey bee larvae were infected with DWV (1.55 × 1010 copies/per larva). No significant mortality occurred in infected honey bee larvae, while the survival rates decreased significantly at the pupal stage. Examination of DWV replication revealed that viral replication began at 2 days post inoculation (d.p.i.), increased dramatically to 4 d.p.i., and then continuously increased in the pupal stage. To better understand the impact of DWV on the larval stage, DWV-infected and control groups were subjected to transcriptomic analysis at 4 d.p.i. Two hundred fifty-five differentially expressed genes (DEGs) (fold change ≥ 2 or ≤ -2) were identified. Of these DEGs, 168 genes were downregulated, and 87 genes were upregulated. Gene Ontology (GO) analysis showed that 141 DEGs (55.3%) were categorized into molecular functions, cellular components and biological processes. One hundred eleven genes (38 upregulated and 73 downregulated) were annotated by KO (KEGG Orthology) pathway mapping and involved metabolic pathways, biosynthesis of secondary metabolites and glycine, serine and threonine metabolism pathways. Validation of DEGs was performed, and the related gene expression levels showed a similar tendency to the DEG predictions at 4 d.p.i.; cell wall integrity and stress response component 1 (wsc1), cuticular protein and myo-inositol 2-dehydrogenase (iolG) were significantly upregulated, and small conductance calcium-activated potassium channel protein (SK) was significantly downregulated at 4 d.p.i. Related gene expression levels at different d.p.i. revealed that these DEGs were significantly regulated from the larval stage to the pupal stage, indicating the potential impacts of gene expression levels from the larval to the pupal stages. Taken together, DWV infection in the honey bee larval stage potentially influences the gene expression levels from larvae to pupae and reduces the survival rate of the pupal stage. This information emphasizes the consequences of DWV prevalence in honey bee larvae for apiculture.

Genomic Sequencing and Comparison of Sacbrood Viruses from Apis cerana and Apis mellifera in Taiwan.

Sacbrood virus (SBV) was the first identified bee virus and shown to cause serious epizootic infections in the population of Apis cerana in Taiwan in 2015. Herein, the whole genome sequences of SBVs in A. cerana and A. mellifera were decoded and designated AcSBV-TW and AmSBV-TW, respectively. The whole genomes of AcSBV-TW and AmSBV-TW were 8776 and 8885 bp, respectively, and shared 90% identity. Each viral genome encoded a polyprotein, which consisted of 2841 aa in AcSBV-TW and 2859 aa in AmSBV-TW, and these sequences shared 95% identity. Compared to 54 other SBVs, the structural protein and protease regions showed high variation, while the helicase was the most highly conserved region among SBVs. Moreover, a 17-amino-acid deletion was found in viral protein 1 (VP1) region of AcSBV-TW compared to AmSBV-TW. The phylogenetic analysis based on the polyprotein sequences and partial VP1 region indicated that AcSBV-TW was grouped into the SBV clade with the AC-genotype (17-aa deletion) and was closely related to AmSBV-SDLY and CSBV-FZ, while AmSBV-TW was grouped into the AM-genotype clade but branched independently from other AmSBVs, indicating that the divergent genomic characteristics of AmSBV-TW might be a consequence of geographic distance driving evolution, and AcSBV-TW was closely related to CSBV-FZ, which originated from China. This 17-amino-acid deletion could be found in either AcSBV or AmSBV in Taiwan, indicating cross-infection between the two viruses. Our data revealed geographic and host specificities between SBVs. The amino acid difference in the VP1 region might serve as a molecular marker for describing SBV cross-infection.

Screening of Differentially Expressed Microsporidia Genes from Nosema ceranae Infected Honey Bees by Suppression Subtractive Hybridization.

The microsporidium Nosema ceranae is a high prevalent parasite of the European honey bee (Apis mellifera). This parasite is spreading across the world into its novel host. The developmental process, and some mechanisms of N. ceranae-infected honey bees, has been studied thoroughly; however, few studies have been carried out in the mechanism of gene expression in N. ceranae during the infection process. We therefore performed the suppressive subtractive hybridization (SSH) approach to investigate the candidate genes of N. ceranae during its infection process. All 96 clones of infected (forward) and non-infected (reverse) library were dipped onto the membrane for hybridization. A total of 112 differentially expressed sequence tags (ESTs) had been sequenced. For the host responses, 20% of ESTs (13 ESTs, 10 genes, and 1 non-coding RNA) from the forward library and 93.6% of ESTs (44 ESTs, 28 genes) from the reverse library were identified as differentially expressed genes (DEGs) of the hosts. A high percentage of DEGs involved in catalytic activity and metabolic processes revealed that the host gene expression change after N. ceranae infection might lead to an unbalance of physiological mechanism. Among the ESTs from the forward library, 75.4% ESTs (49 ESTs belonged to 24 genes) were identified as N. ceranae genes. Out of 24 N. ceranae genes, nine DEGs were subject to real-time quantitative reverse transcription PCR (real-time qRT-PCR) for validation. The results indicated that these genes were highly expressed during N. ceranae infection. Among nine N. ceranae genes, one N. ceranae gene (AAJ76_1600052943) showed the highest expression level after infection. These identified differentially expressed genes from this SSH could provide information about the pathological effects of N. ceranae. Validation of nine up-regulated N. ceranae genes reveal high potential for the detection of early nosemosis in the field and provide insight for further applications.

Dynamics of Apis cerana Sacbrood Virus (AcSBV) Prevalence in Apis cerana (Hymenoptera: Apidae) in Northern Taiwan and Demonstration of its Infection in Apis mellifera (Hymenoptera: Apidae).

Since 2016, Apis cerana sacbrood virus (AcSBV) has been recorded in Taiwan. It is epizootic in Apis cerana (Hymenoptera: Apidae) and causing serious loss of A. cerana. Herein, we performed a long-term survey of AcSBV prevalence in the populations of A. cerana in Northern Taiwan from January 2017 to July 2018. The surveillance of AcSBV prevalence in A. mellifera (Hymenoptera: Apidae) populations was starting and further confirmed by sequencing since April 2017; thus, these data were also included in this survey. In our survey, the average prevalence rates of AcSBV were 72 and 53% in A. cerana and A. mellifera, respectively, in 2017, which decreased to 45 and 27% in 2018. For the spatial analysis of AcSBV in two honey bee populations, Hsinchu showed the highest prevalence, followed by New Taipei, Yilan, Taipei, and Keelung, suggesting that AcSBV might have come from the southern part of Taiwan. Interestingly, the AcSBV prevalence rates from A. cerana and A. mellifera cocultured apiaries gradually synchronized. The result of phylogenetic analysis and comparison of the annual AcSBV prevalence in A. cerana-only, A. mellifera-only, and A. cerana/A. mellifera cocultured sample sites indicate cross-infection between A. cerana and A. mellifera; however, AcSBV may lose the advantage of virulence in A. mellifera. The evidence suggested that the transmission of AcSBV might occur among these two honey bee species in the field. Therefore, A. mellifera may serve as a guard species to monitor AcSBV in A. cerana, but the cross-infection still needs to be surveyed.

Genomic sequencing of Troides aeacus nucleopolyhedrovirus (TraeNPV) from golden birdwing larvae (Troides aeacus formosanus) to reveal defective Autographa californica NPV genomic features.

BACKGROUND: The golden birdwing butterfly (Troides aeacus formosanus) is a rarely observed species in Taiwan. Recently, a typical symptom of nuclear polyhedrosis was found in reared T. aeacus larvae. From the previous Kimura-2 parameter (K-2-P) analysis based on the nucleotide sequence of three genes in this isolate, polh, lef-8 and lef-9, the underlying virus did not belong to any known nucleopolyhedrovirus (NPV) species. Therefore, this NPV was provisionally named "TraeNPV". To understand this NPV, the nucleotide sequence of the whole TraeNPV genome was determined using next-generation sequencing (NGS) technology. RESULTS: The genome of TraeNPV is 125,477 bp in length with 144 putative open reading frames (ORFs) and its GC content is 40.45%. A phylogenetic analysis based on the 37 baculoviral core genes suggested that TraeNPV is a Group I NPV that is closely related to Autographa californica nucleopolyhedrovirus (AcMNPV). A genome-wide analysis showed that TraeNPV has some different features in its genome compared with other NPVs. Two novel ORFs (Ta75 and Ta139), three truncated ORFs (pcna, he65 and bro) and one duplicated ORF (38.7 K) were found in the TraeNPV genome; moreover, there are fewer homologous regions (hrs) than there are in AcMNPV, which shares eight hrs within the TraeNPV genome. TraeNPV shares similar genomic features with AcMNPV, including the gene content, gene arrangement and gene/genome identity, but TraeNPV lacks 15 homologous ORFs from AcMNPV in its genome, such as ctx, host cell-specific factor 1 (hcf-1), PNK/PNL, vp15, and apsup, which are involved in the auxiliary functions of alphabaculoviruses. CONCLUSIONS: Based on these data, TraeNPV would be clarified as a new NPV species with defective AcMNPV genomic features. The precise relationship between TraeNPV and other closely related NPV species were further investigated. This report could provide comprehensive information on TraeNPV for evolutionary insights into butterfly-infected NPV.

Characterization and functional assay of apsup (Lyxy105) from Lymantria xylina multiple nucleopolyhedrovirus (LyxyMNPV).

The baculoviral anti-apoptotic genes, p35 and iap (inhibitor of apoptosis), play important roles in the initiation of viral infection. Recently, a new anti-apoptotic gene (apoptosis suppressor, apsup) was identified in Lymantria dispar multiple nucleopolyhedrovirus (LdMNPV). An apsup homolog gene, Lyxy105 (ly-apsup), was also predicted in the Lymantria xylina multiple nucleopolyhedrovirus (LyxyMNPV) genome. In this study, we attempt to perform a gene expression analysis and a functional assay of ly-apsup to demonstrate its anti-apoptotic activity and identify the functional domain of this protein. The transcription of the ly-apsup gene region was detected from 12 h post-infection (hpi) and increased significantly after 24-72 hpi. Comparison of the putative amino acid sequences to those of 18 baculoviral homolog proteins showed high amino acid identity to the LdMNPV sequences. Moreover, five conserved protein domains (named as domains I-V) were found. Therefore, protein functional assays were conducted on full-length proteins and different truncation clones. The overexpression of each clone was confirmed by western blot analysis, and the data revealed that a cleavage of ~ 5 kDa at the N-terminal region of the full-length, domains I-IV (1-241) and I-III (1-178), proteins occurred. The results of the functional analysis showed that full-length Ly-apsup and Ly-apsup with domain I (1-70) could inhibit Drosophila-RPR protein (D-RPR)-induced and actinomycin D (ActD)-induced apoptoses. In addition, the domains I and I-II (1-126) regions showed higher anti-apoptotic activity than the other domains in both D-RPR-induced and ActD-induced cell apoptoses. In conclusion, domain I of Ly-apsup may play an important role in the anti-apoptotic activity of this protein; cleavage of the Ly-apsup N-terminus may lead to decreased anti-apoptotic activity.