Dimethyl sulfoxide, an alternative for control of Nosema ceranae infection in honey bees (Apis mellifera).

Nosema ceranae is a microsporidian parasite that threatens current apiculture. N. ceranae-infected honey bees (Apis mellifera) exhibit morbid physiological impairments and reduced honey production, malnutrition, shorter life span, and higher mortality than healthy honey bees. In this study, we found that dimethyl sulfoxide (DMSO) could enhance the survival rate of N. ceranae-infected honey bees. Therefore, we investigated the effect of DMSO on N. ceranae-infected honey bees using comparative RNA sequencing analysis. Our results revealed that DMSO was able to affect several biochemical pathways, especially the metabolic-related pathways in N. ceranae-infected honey bees. Based on these findings, we conclude that DMSO may be a useful alternative for treating N. ceranae infection in apiculture.

Monitoring the Season-Prevalence Relationship of Vairimorpha ceranae in Honey Bees (Apis mellifera) over One Year and the Primary Assessment of Probiotic Treatment in Taichung, Taiwan.

Microsporidiosis, which is caused by the pathogen Vairimorpha ceranae, is a prevalent disease in the honey bee (Apis mellifera) and might lead to significant adult honey bee mortality. In this study, we conducted an annual survey of the mature spore load of V. ceranae in the guts of nurse bees and forager bees in the apiary of National Chung Hsing University (NCHU) in Taiwan. The results indicated that, on average, honey bees hosted approximately 2.13 × 106 mature spore counts (MSCs)/bee in their guts throughout the entire year. The highest number of MSCs was 6.28 × 106 MSCs/bee, which occurred in April 2020, and the lowest number of MSCs was 5.08 × 105 MSCs/bee, which occurred in November 2020. Furthermore, the guts of forager bees had significantly higher (>58%) MSCs than those of nurse bees. To evaluate the potential of the probiotic to treat microsporidiosis, the lactic acid bacterium Leuconostoc mesenteroides TBE-8 was applied to honey bee colonies. A significant reduction (>53%) in MSCs following probiotic treatment was observed, indicating the potential of probiotic treatment for managing microsporidiosis. This research provided information on V. ceranae MSCs in the honey bee gut at NCHU in Taiwan and the MSCs' correlation with the annual season. Furthermore, a potential probiotic treatment for microsporidiosis was assessed for future management.

Assessment of Aphidicidal Effect of Entomopathogenic Fungi against Parthenogenetic Insect, Mustard Aphid, Lipaphis erysimi (Kalt.).

The mustard aphid (L. erysimi) is a pest that infests various cruciferous crops and transmits plant viruses. To achieve eco-friendly pest management, entomopathogenic fungi (EPF) are potential microbial control agents for controlling this pest. Therefore, virulence screening of EPF isolates under Petri dish conditions is necessary before field application. However, the mustard aphid is a parthenogenetic insect, making it difficult to record data during Petri dish experiments. A modified system for detached-leaf bioassays was developed to address this issue, using a micro-sprayer to inoculate conidia onto aphids and prevent drowning by facilitating air-drying after spore suspension. The system maintained high relative humidity throughout the observation period, and the leaf disc remained fresh for over ten days, allowing parthenogenetic reproduction of the aphids. To prevent offspring buildup, a process of daily removal using a painting brush was implemented. This protocol demonstrates a stable system for evaluating the virulence of EPF isolates against mustard aphids or other aphids, enabling the selection of potential isolates for aphid control.

Investigation of the fall armyworm (Spodoptera frugiperda) gut microbiome and entomopathogenic fungus-induced pathobiome.

The gut microflora plays an important role in insect development and physiology. The gut bacterial microbiome of the fall armyworm (FAW), Spodoptera frugiperda, in both cornfield and laboratory-reared populations was investigated using a 16S metagenomic approach. The alpha- and beta-diversity of the cornfield FAW populations varied among sampling sites and were higher than those of the laboratory-reared FAW population, indicating that different diets and environments influence the gut bacterial composition. To better understand the interaction between the microbiome and entomopathogenic fungi (EPF), FAWs from organic and conventionally managed corn fields and from the laboratory-reared colony were inoculated with Beauveria bassiana NCHU-153 (Bb-NCHU-153). A longer median lethal time (LT50) was observed in the Bb-NCHU-153-infected cornfield FAW population than in the laboratory-reared FAWs. In terms of the microbiome, three Bb-NCHU-153-infected FAW groups showed different gut bacterial compositions compared to noninfected FAW. Further investigation of the cooccurrence network and linear discriminant analysis (LDA) of effect size (LEfSe) revealed that the enriched bacterial genera, such as Enterococcus, Serratia, Achromobacter, and Tsukamurella, in the gut might play the role of opportunistic pathogens after fungal infection; in contrast, some gut bacteria of Methylobacterium, Marinomonas, Paenochrobactrum, Pseudomonas, Acinetobacter, Delftia, Dietzia, Gordonia, Leucobacter, Paracoccus, and Stenotrophomonas might be probiotics against EPF infection. These results indicated that EPF infection can change the gut bacterial composition and lead to a pathobiome in the FAW and that some bacterial species might protect the FAW from EPF infection. These findings could be applied to the design of pathobiome-inducing biocontrol strategies.

Whole-genome DNA methylome analysis of different developmental stages of the entomopathogenic fungus Beauveria bassiana NCHU-157 by nanopore sequencing.

The entomopathogenic fungus (EPF), Beauveria bassiana, is an important and commonly used EPF for microbial control. However, the role of DNA methylation has not been thoroughly studied. Therefore, the whole genomic DNA methylome of one promising EPF isolate, B. bassiana NCHU-157 (Bb-NCHU-157), was investigated by Oxford Nanopore Technologies (ONT). First, the whole genome of Bb-NCHU-157 was sequenced by next-generation sequencing (NGS) and ONT. The genome of Bb-NCHU-157 contains 16 contigs with 34.19 Mb and 50% GC content, which are composed of 10,848 putative protein-coding genes. Two putative DNA methyltransferases (DNMTs) were found, including Dim-2 and C-5 cytosine-specific DNA methylases. Both DNMTs showed higher expression levels in the mycelium stage than in the conidia stage, indicating that development of DNA methylation in Bb-NCHU-157 might occur in the mycelium stage. The global methylation level of the mycelium stage (5 mC = 4.56%, CG = 3.33%, CHG = 0.74%, CHH = 0.49%) was higher than that of the conidial stage (5 mC = 2.99%, CG = 1.99%, CHG = 0.63%, CHH = 0.37%) in both the gene and transposable element (TE) regions. Furthermore, the TE regions showed higher methylation frequencies than the gene regions, especially for CHH site methylation, suggesting regulation of genomic stabilization during mycelium development. In the gene regions, high methylation frequencies were found around the transcription start site (TSS) and transcription end site (TES). Moreover, CG and CHG methylation mainly occur in the promoter and intergenic regions, while CHH methylation occurs in the TE region. Among the methylated regions, 371, 661, and 756 differentially DNA methylated regions (DMRs) were hypermethylated in the mycelium in CG, CHG, and CHH, while only 13 and 7 DMRs were hypomethylated in the mycelium in CHG, and CHH, respectively. Genes located in the DMR shared the GO terms, DNA binding (GO: 0003677), and sequence-specific DNA binding (GO: 0043565) for hypermethylation in the mycelium, suggesting that methylation might regulate gene expression from the initial process. Evaluation of the DNA methylome in Bb-NCHU-157 by ONT provided new insight into this field. These data will be further validated, and epigenetic regulation during the development of B. bassiana will be explored.

Evaluation of the Potential Entomopathogenic Fungi Purpureocillium lilacinum and Fusarium verticillioides for Biological Control of Forcipomyia taiwana (Shiraki).

Forcipomyia taiwana (Diptera: Ceratopogonidae) is a nuisance blood-sucking pest to humans in Taiwan. An F. taiwana bite causes itching and redness and usually causes serious harassment to human outdoor activity. In terms of F. taiwana control, chemical pesticides are ineffective. Therefore, other efforts are needed. Fungal mycosis in the larvae, pupae, and emerging F. taiwana adults was found during the rearing of F. taiwana. In this study, six fungal isolates were isolated from infected cadavers and subjected to molecular identification. In addition, their biocontrol potential was evaluated against different life stages of F. taiwana. Based on the pathogenicity screening, two fungal isolates, NCHU-NPUST-175 and -178, which caused higher mortality on the fourth instar larvae of F. taiwana, were selected for virulence tests against different life stages of F. taiwana larvae. The results of the phylogenetic analysis indicated that the NCHU-NPUST-175 and -178 belonged to Purpureocillium lilacinum and Fusarium verticillioides, respectively. Bioassay against different life stages of F. taiwana with different spore concentrations (5 × 105 to 5 × 107 conidia/mL) revealed a dose-dependent effect on larvae for both fungal isolates, while only 38% and 50% mortality was found in highest concentration (5 × 107 conidia/mL) at fourth instar larvae by Pl-NCHU-NPUST-175 and Fv-NCHU-NPUST-178, respectively. Moreover, reductions in egg-hatching rate and adult emergence rate were found, when the last stage of F. taiwana was inoculated with both fungal isolates, indicating the ovicidal potential and the impact of entomopathogenic fungi on the development of F. taiwana. In conclusion, Pl-NCHU-NPUST-175 and Fv-NCHU-NPUST-178 showed larvicidal activity, ovicidal activity, and impact on adult emergence on F. taiwana.

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.

A Novel Application of 3D Printing Technology Facilitating Shell Wound Healing of Freshwater Turtle.

Numerous cases and a shortage of resources usually limit wild animal rescue. New technology might save these severely injured wild animals from euthanasia by easing the requirement of intensive medication. Three-dimensional (3D) technologies provide precise and accurate results that improve the quality of medical applications. These 3D tools have become relatively low-cost and accessible in recent years. In the medical field of exotic animals, turtle shell defects are highly challenging because of inevitable water immersion. This report is the first attempt to apply the combination of 3D scanning, computer-aided design (CAD), and 3D printing to make a device that protects the wound from exposure to water or infection sources. The presented techniques successfully extricate a wild freshwater turtle from an extensive shell defect within a short period. Integration of multiple sciences to 3D technology can provide a facile model for veterinary medical applications.

Low-Level Fluvalinate Treatment in the Larval Stage Induces Impaired Olfactory Associative Behavior of Honey Bee Workers in the Field.

Fluvalinate is a widely used insecticide for varroa mite control in apiculture. While most beekeepers have ignored the effects of low levels of fluvalinate on bees, this study aims to demonstrate its effects at very low concentrations. We first used fluvalinate doses ranging from 0.4 to 400 ng/larva to monitor the capping, pupation, and emergence rates of larval bees. Second, we used the honey bees' proboscis extension reflex reaction to test the learning ability of adult bees that were exposed to fluvalinate doses from 0.004 to 4 ng/larva in the larval stage. The brood-capped rate of larvae decreased dramatically when the dose was increased to 40 ng/larva. Although no significant effect was observed on brood-capping, pupation, and eclosion rates with a dose of 4 ng/larva, we found that the olfactory associative behavior of adult bees was impaired when they were treated with sublethal doses from 0.004 to 4 ng/larva in the larval stage. These findings suggest that a sublethal dose of fluvalinate given to larvae affects the subsequent associative ability of adult honey bee workers. Thus, a very low dose may affect the survival conditions of the entire colony.

The complete mitochondrial genome of Attacus atlas formosanus Villiard, 1969 (Lepidoptera: Saturniidae).

The complete mitochondrial genome (mitogenome) of Attacus atlas formosanus (Villiard, 1969) is 15,280 bp in length, with the typical gene content and arrangement usually observed in Insecta. It contains 13 protein-coding genes, 22 transfer RNA genes, two ribosomal RNA genes, and one AT-rich region. The overall nucleotide composition of the mitogenome was 39.8% A, 12.9% C, 7.7% G, and 39.6% T, with an A + T bias of 79.4%. Phylogenetic analyses of 23 species in Saturniidae and 3 species in Bombycidae by Bayesian inference showed that A. atlas formosanus belonged to the Tribe Attacini, closely related to Tribe Saturniini. Besides, A. atlas formosanus is closely related to A. atlas with 99% sequence identity. This result well supported the taxonomic position of Saturniidae and their close relationship with the family Bombycidae.

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.

Isolation and Selection of Entomopathogenic Fungi from Soil Samples and Evaluation of Fungal Virulence against Insect Pests.

Entomopathogenic fungi (EPF) are one of the microbial control agents for integrated pest management. To control local or invasive pests, it is important to isolate and select indigenous EPF. Therefore, the soil bait method combined with the insect bait (mealworm, Tenebrio molitor) system was used in this study with some modifications. The isolated EPF were then subjected to the virulence test against the agricultural pest Spodoptera litura. Furthermore, the potential EPF strains were subjected to morphological and molecular identifications. In addition, the conidia production and thermotolerance assay were performed for the promising EPF strains and compared; these data were further substituted into the formula of effective conidia number (ECN) for laboratory ranking. The soil bait-mealworm system and the ECN formula can be improved by replacing insect species and integrating more stress factors for the evaluation of commercialization and field application. This protocol provides a quick and efficient approach for EPF selection and will improve the research on biological control agents.

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.

Deformed Wing Virus in Two Widespread Invasive Ants: Geographical Distribution, Prevalence, and Phylogeny.

Spillover of honey bee viruses have posed a significant threat to pollination services, triggering substantial effort in determining the host range of the viruses as an attempt to understand the transmission dynamics. Previous studies have reported infection of honey bee viruses in ants, raising the concern of ants serving as a reservoir host. Most of these studies, however, are restricted to a single, local ant population. We assessed the status (geographical distribution/prevalence/viral replication) and phylogenetic relationships of honey bee viruses in ants across the Asia-Pacific region, using deformed wing virus (DWV) and two widespread invasive ants, Paratrechina longicornis and Anoplolepis gracilipes, as the study system. DWV was detected in both ant species, with differential geographical distribution patterns and prevenance levels between them. These metrics, however, are consistent across the geographical range of the same ant species. Active replication was only evident in P. longicornis. We also showed that ant-associated DWV is genetically similar to that isolated from Asian populations of honey bees, suggesting that local acquisition of DWV by the invasive ants may have been common at least in some of our sampled regions. Transmission efficiency of DWV to local arthropods mediated by ant, however, may vary across ant species.

Comparison of gut microbiota of healthy and diseased walking sticks, Phasmotaenia lanyuhensis.

Research on gut microbiota of phytophagous insects has shown to be important for the physiological functions of insect hosts; however, little is known about the changes in gut microbiota when they are suffering from environmental stress or pathogen infections. During rearing of Phasmotaenia lanyuhensis (Phasmatodea: Phasmatidae), sluggish locomotion was usually followed by the death of the insect with a symptom of melanization in the front part of the abdomen. Therefore, the abnormal individuals were initially classified into moribund, light- and serious-symptom based on the level of abnormal physiological circumstances and melanization. The gut microbiota of these samples were further investigated by 16S metagenomic sequencing and the differences in bacterial abundance and structure of bacterial community were analyzed. A decrease in microbiota diversity was observed in the diseased P. lanyuhensis, with the abundance of phyla Proteobacteria and Firmicute relatively higher compared to those without symptom. Interestingly, principal component analysis based on the bacterial richness was correlated to the level of melanization symptom in the diseased P. lanyuhensis, suggested the change in bacterial microbiota involved in this abnormal circumstance. However, the factor that caused the initial alternation of microbiota remains to be identified. Additionally, the lack of bacterial diversity (i.e., absence of Meiothermus and Nubsella spp.) in P. lanyuhensis might reduce the fitness for surviving. This report provided the comprehensive microbiota analysis for P. lanyuhensis and concluded that either the relative abundance or the bacterial diversity of microbiota in the insect digestive system may influence the physiological functions of phytophagous insects.

Complete Genome Sequence of a Novel Putative RNA Virus, RiPV-2, from the Bean Bug Riptortus pedestris.

A novel putative single-stranded RNA virus was discovered from the transcriptome of a bean bug, Riptortus pedestris, infected with the entomopathogenic fungus Beaveria bassiana JEF-007. The complete genome sequence was 9,915 nucleotides long and encoded a 2,916-amino-acid polyprotein. This virus belonged to Iflaviridae based on phylogenetic analysis and was named RiPV-2.

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.

Diversity of Fungal DNA Methyltransferases and Their Association With DNA Methylation Patterns.

DNA methyltransferases (DNMTs) are a group of proteins that catalyze DNA methylation by transferring a methyl group to DNA. The genetic variation in DNMTs results in differential DNA methylation patterns associated with various biological processes. In fungal species, DNMTs and their DNA methylation profiles were found to be very diverse and have gained many research interests. We reviewed fungal DNMTs in terms of their biological functions, protein domain structures, and their associated epigenetic regulations compared to those known in plant and animal systems. In addition, we summarized recent reports on potential RNA-directed DNA methylation (RdDM) related to DNMT5 in fungi. We surveyed up to 40 fungal species with published genome-wide DNA methylation profiles (methylomes) and presented the associations between the specific patterns of fungal DNA methylation and their DNMTs based on a phylogenetic tree of protein domain structures. For example, the main DNMTs in Basidiomycota, DNMT1 with RFD domain + DNMT5, contributing to CG methylation preference, were distinct from RID + Dim-2 in Ascomycota, resulting in a non-CG methylation preference. Lastly, we revealed that the dynamic methylation involved in fungal life stage changes was particularly low in mycelium and DNA methylation was preferentially located in transposable elements (TEs). This review comprehensively discussed fungal DNMTs and methylomes and their connection with fungal development and taxonomy to present the diverse usages of DNA methylation in fungal genomes.