Abstract
Standing genetic variation is the predominant source acted on by selection. Organisms with high genetic diversity generally show faster responses toward environmental change. Nosema ceranae is a microsporidian parasite of honey bees, infecting midgut epithelial cells. High genetic diversity has been found in this parasite, but the mechanism for the parasite to maintain this diversity remains unclear. This study involved continuous inoculation of N. ceranae to honey bees. We found that the parasites slowly increased genetic diversity over three continuous inoculations. The number of lost single nucleotide variants (SNVs) was balanced with novel SNVs, which were mainly embedded in coding regions. Classic allele frequency oscillation was found at the regional level along the genome, and the associated genes were enriched in apoptosis regulation and ATP binding. The ratio of synonymous and non-synonymous substitution suggests a purifying selection, and our results provide novel insights into the evolutionary dynamics in microsporidian parasites.
Introduction
Microsporidia comprise a large group of obligate intracellular parasites, which forms an early branch in the fungal kingdom (Capella-Gutiérrez et al., 2012; Choi and Kim, 2017). Microsporidian genomes were generally compact, ranging between 2 and 5 Mbp (Ndikumana et al., 2017), arguably because they rely on host cells for much of their metabolic processes. Interestingly, one exceptionally large genome of 51 Mbp has been reported, and this species maintains transportable elements, which are exceedingly rare in microsporidia generally (Desjardins et al., 2015). Conventional mitochondria are lost during evolution, leaving tiny mitochondrially derived organelles called mitosomes (Burri et al., 2006). Microsporidian infections have been reported to inhibit host apoptosis, suppress host immune responses, and enhance host metabolism (Scanlon et al., 1999; del Aguila et al., 2006; Antunez et al., 2009; Cuomo et al., 2012; Higes et al., 2013; Martín-Hernández et al., 2017). More recently, a parasitic miRNA was found to be a virulence factor for one microsporidian to establish infection in silkworms (Dong et al., 2021).
Nosema ceranae is a microsporidian parasite that infects epithelial cells of honey bee midgut tissue (Higes et al., 2007). N. ceranae was first described in the Asian honey bee Apis cerana, which has established infection in the European honey bee Apis mellifera globally (Fries et al., 1996; Higes et al., 2006; Huang et al., 2007). The genome size of N. ceranae is 8.8 Mbps, with 2,280 protein-coding genes (Huang et al., 2021). The infection starts with the ingestion of N. ceranae spore-contaminated nectar. A parasite proliferation cycle lasts approximately 4 days, defined by life stages from meronts and sporonts to mature spores (Higes et al., 2007). Infected cells eventually burst to release a large number of spores. Parasites can then germinate to infect neighboring healthy epithelial cells or can be expelled into the environment through feces. Honey bees infected by N. ceranae showed reduced life span, enhanced energetic burden, and impaired navigation and memory (Antunez et al., 2009; Mayack and Naug, 2009; Eiri et al., 2015; Gage et al., 2018). In extreme conditions, the infection may lead to colony collapse (Higes et al., 2008; Botías et al., 2013). Using a set of genetic markers, high genetic diversity of N. ceranae has been found. Surprisingly, this genetic diversity can be even higher within a colony than among colonies (Gómez-Moracho et al., 2014, 2015). Additionally, the genetic variation of N. ceranae was minor among geographically distant locations (Pelin et al., 2015). The host shapes the parasite diversity, and the parasite drives the host diversification (Kamiya et al., 2014; Betts et al., 2018; McNew et al., 2021). It remains unclear how N. ceranae could maintain high levels of polymorphism. We hypothesized that mutation and frequency-dependent selection together shape the adaptation of the parasite. To validate the above hypothesis, the single nucleotide variants (SNVs) of N. ceranae were quantified over three successive generations to assess the frequency of novel and lost SNVs, and allele oscillation from parental to offspring spores.
Materials and methods
Parasite isolation, inoculation, and sample collection
To harvest sufficient spores, 300 honey bee foragers were randomly collected from five heavily infected honey bee (A. mellifera) colonies in the experimental apiary of Jiangxi Agricultural University following a regular survey. Midguts of these honey bees were dissected and homogenized to isolate N. ceranae spores via centrifugation (Fries et al., 2013). Then, the spores were purified using Percoll gradient centrifugation and confirmed as pure N. ceranae by species-specific fragment size of a conventional PCR (Chen et al., 2013; Fries et al., 2013). The purified spores were defined as F0_spores. A subset of these F0_spores was used to extract genomic DNA using cetyltrimethylammonium bromide (CTAB) (Chen et al., 2013). The remaining F0_spores were used to inoculate newly emerged honey bees. Sealed honey bee (A. mellifera) brood frames were kept in an incubator (34 ± 1°C, 60% relative humidity). In total, 200 newly emerged honey bee workers were individually fed with 2 μl sucrose solution containing 105N. ceranae F0 spores.
Additionally, 200 newly emerged honey bee workers were fed with 2 μl sucrose solution as the F1 uninfected group. During the experiment, cohorts were divided into 4 cups containing 50 bees each and maintained on 50% sucrose solution ad libitum in the incubator. To harvest sufficient spores to infect a new batch of honey bees, the midgut was dissected from individual honey bees at 14 days postinfection (dpi), and F1 spores were collected. A subset of F1 spores was again used to extract genomic DNA, and the remaining F1 spores were used to inoculate 100 newly emerged honey bees due to a limited number of spores. In total, 100 newly emerged honey bee workers were fed with 2 μl sucrose solution as the F2 uninfected group. At 14 dpi, the midgut was dissected from individual honey bees and then pooled for spore purification, defined as F2 spores. The experiment was duplicated. Overall, two libraries for F0_spores, two libraries for F1 spores, and two libraries for F2 spores were prepared and sequenced on Illumina Hiseq 2000 platform in Novogene Co.
Bioinformatic and statistics
The quality of DNA sequencing reads was processed with Fastp and validated with FastQC (Chen et al., 2018). DNA sequencing reads were aligned to the N. ceranae genome (Ncer 3.0, GCA_004919615.1) using BWA with default parameters (Li and Durbin, 2009; Huang et al., 2021). The SNVs were identified and annotated using the Picard-GATK-SNPEFF pipeline (Van der Auwera et al., 2013). The SNVs with a quality score of less than 50 and locus alignment coverage of less than 10 were removed. Enrichment GO terms were analyzed using TopGO, R (Alexa and Rahnenfuhrer, 2021). The SNVs among F0, F1, and F2 spores were analyzed using chi-squared tests implemented using R (R Core Team, 2013). The lost and novel SNVs among the two generations were analyzed using paired t-test, implemented using R (R Core Team, 2013). The ratio of πnon–synonymous/πsynonymous (πa/πs) was calculated using SNPGenie (Nelson et al., 2015). The genome diversity π, Watterson’s θ, and corrected Tajima’s D were calculated using Popoolation (Kofler et al., 2011a). The fixation index FST was calculated using Popoolation2 (Kofler et al., 2011b). The recombination rate was calculated based on the allele frequency using ReLERNN (Adrion et al., 2020). DNA sequencing reads have been deposited in NCBI under BioProject PRJNA784016.
Results
The number of novel and lost single nucleotide variants were balanced between parental and offspring spores
Nosema ceranae spores were not found in the uninfected groups. For the spores purified from infected honey bees, 104 ± 7 million reads (mean ± SD) were aligned to the N. ceranae genome (150 bp paired reads, over 1,000 genome coverage), which accounted for 93.2% of all sequenced reads. In F0 spores, 97,279 and 98,314 SNVs were identified in two replicates with 96.8% overlap, reflecting approximately 11 SNVs per Kbp (Figure 1). Additionally, over 95% SNVs were heterozygous. Within coding regions, 7,352 non-synonymous and 8,695 synonymous SNVs were identified in F0 spores. From F0 to F1 spores, 103 synonymous and 200 non-synonymous SNVs were lost. Comparatively, 314 synonymous and 350 non-synonymous novel SNVs were gained in F1 spores. From F1 to F2 spores, 302 synonymous and 314 non-synonymous SNVs were lost. Meanwhile, 982 synonymous and 1,228 non-synonymous novel SNVs were gained. The dynamic of non-synonymous SNVs was significantly stronger than synonymous SNVs from F0 to F2 spores (paired t-test, P < 0.01). Among the novel non-synonymous SNVs in F1 spores, 198 SNVs were maintained in F2 spores located within 54 genes. Using the genome as the background, the 54 genes were significantly enriched in the metabolism pathway (Fisher’s exact test, P < 0.001). After normalizing the contig length, the contig SMUP01000049.1 (13.1 novel SNVs per 10 kbp), SMUP01000091.1 (8.6 novel SNVs per 10 kbp), SMUP01000014.1 (8.0 novel SNVs per 10 kbp), SMUP01000071.1 (6.7 novel SNVs per 10 kbp), and SMUP01000003.1 (6.2 novel SNVs per 10 kbp) showed a significantly higher tendency to generate novel SNVs compared with the genome (one sample t-test, adjusted P < 0.05) (Figure 2A). There were 269 genes embedded in the above contigs, which were enriched in the cell cycle (GO: 0010971, adjusted P < 0.05) and ATP binding (GO: 0005524, adjusted P < 0.001).
FIGURE 1
FIGURE 2
Genome recombination and allele frequency oscillation from parental to offspring spores
A total of 87,129 SNVs were inherited from F0 to F2 spores in both replicates. Over 95% of loci were heterozygous with an average alternative allele frequency of 0.38. Contigs with lower recombination rates were rarely found. Most contigs were homogeneous with 5.14e-9 cross-over per bp (Figure 2B). To quantify the allele oscillation, the allele frequency variation was calculated between parents and offspring, denoted as coefficient S. When offspring inherits the same allele frequency from their parents, the coefficient S equals 0. In our data, the coefficient S significantly deviated from F0 to F1, as well as from F1 to F2 in both replicates (one sample t-test, adjusted P < 0.001) (Figure 3). By comparing with the genome, the genes associated with the selected loci were significantly enriched in negative regulation of apoptosis (GO: 0043066, adjusted P < 0.05) and ATP binding (GO: 0005524, adjusted P < 0.001).
FIGURE 3
Genetic differentiation accumulated over generations
Genome diversity (π) slowly increased from F0 (0.0084 ± 0.0001) and F1 (0.0087 ± 0.0001) to F2 (0.0088 ± 0.0001). Tajima’s D was positive, and the ratio of πa/πs was less than 1 for the parasites in all three generations. The fixation index FST was 0.0015 ± 0.0001 between F0 and F1 spores, which increased by twofold (0.0042 ± 0.0001) between F1 and F2 spores. After two parasite generations, the index FST increased by fourfold from F0 to F2 spores (Table 1).
TABLE 1
| Parasites | π a/π s | π | Waterson’s θ | Tajima’s D | F0 (FST) | F1 (FST) |
| F0 | 0.2027 ± 0.0122 | 0.0084 ± 0.0001 | 0.0080 ± 0.0001 | 0.7980 ± 0.0159 | ||
| F1 | 0.2009 ± 0.0107 | 0.0087 ± 0.0001 | 0.0088 ± 0.0001 | 0.3367 ± 0.0152 | 0.0015 ± 0.0001 | |
| F2 | 0.2034 ± 0.0117 | 0.0088 ± 0.0001 | 0.0085 ± 0.0001 | 0.7194 ± 0.0157 | 0.0071 ± 0.0001 | 0.0042 ± 0.0001 |
Population genetic statistics for the three generations of Nosema ceranae spores (mean ± SE).
πa indicates the diversity at non-synonymous sites, and πs indicates the diversity at synonymous sites.
The fixation index FST of the parasites in three generations was pairwise analyzed, which increased substantially with each generation.
Discussion
Host-parasite coevolution is a process of reciprocal allelic oscillation between antagonist species and is often characterized by the Red Queen Hypothesis, which is a central theorem in coevolutionary biology (Liow et al., 2011; Brockhurst et al., 2014; Rabajante et al., 2016; Martínez-López et al., 2021). High-standing genetic variation allows fast responses to environmental change (Lai et al., 2019). N. ceranae is a native parasite of Eastern honey bees (A. cerana), which has successfully established widespread infections in Western honey bees (A. mellifera). The parasite showed high virulence in the novel host, and infected honey bees showed a shortened life span (Eiri et al., 2015), suppressed apoptosis (Antunez et al., 2009; Kurze et al., 2015; Martín-Hernández et al., 2017), and arguably are prone to colony collapses (Higes et al., 2008; Botías et al., 2013). Using a set of genetic markers, the diversity of this parasite was found higher within a honey bee colony than among colonies (Gómez-Moracho et al., 2014, 2015). In our data, the number of identified SNVs was congruent with a previous study (Ke et al., 2021). The numbers of inherited, lost, and novel SNVs were highly congruent between the two replicates, suggesting that selection is not random at the colony level. Additionally, we found minor variance between the replicates, suggesting that genetic diversity is higher within a colony than among colonies, as indicated previously (Gómez-Moracho et al., 2014; Pelin et al., 2015). Previously, the recombination event was reported in N. ceranae (Pelin et al., 2015). In our study, the calculated recombination rate was orders of magnitude smaller than those of other fungal pathogens (Heinzelmann et al., 2020; Wyka et al., 2022). Inhibiting apoptosis and ATP acquisition were essential to the success of N. ceranae infection. When apoptosis was facilitated, infected honey bees were tolerant of the infection (Kurze et al., 2015). Suppressing the expression of ATP transporters also reduced N. ceranae proliferation (Paldi et al., 2010). In our data, the alleles associated with apoptosis and ATP binding were under constant oscillation, which supports the antagonistic allele oscillation under the Red Queen hypothesis (Papkou et al., 2019). Anthropocene could also shape the genome selection of the parasite (Ptaszyńska et al., 2021), as well as the foraging behavior of honey bees. For example, parasites can quickly disperse through commonly visited flowers, which may decrease colony level variance (Graystock et al., 2015). The transporting pollination activities may further enhance the parasite dispersal at a larger scale (Proesmans et al., 2021). Furthermore, plant diversity can strongly impact the honey bee diet quality, which may shape the susceptibility to pathogens (Ptaszyńska et al., 2021). Genetic drift and selection decrease genetic variation, while mutation, gene flow, and horizontal gene transfer increase genetic variation (Schulte et al., 2013; Ellegren and Galtier, 2016). In our data, a strong fraction of SNVs was lost during the bottleneck caused by a single inoculation. Meanwhile, mutations led to ample novel SNVs. As the parasite inoculation was employed in lab conditions, the impacts of gene flow were minimal. This suggests that a balance between lost and novel SNVs maintains the high genetic diversity in N. ceranae. In our data, the accumulating π suggests a slow accumulation of genetic diversity over generations. Genome diversity π was slightly higher than Watterson’s θ, leading to a small positive Tajima’s D, which suggests a small number of low-frequency alleles (Fu and Li, 1993; Stajich and Hahn, 2005). The ratio of πa/πs was less than 1 in all three parasite generations, which suggests that the parasite genome is going through purifying selection as found in historical samples (Huang et al., 2016). The fixation index increased fourfold from F0 to F2 spores, suggesting a genome diversification (Jackson et al., 2015). The parasite can transmit among bee species, which might have posed additional selective pressure on the parasite to adapt to multiple hosts. In yeast, static environments favor small numbers of fit genotypes, dramatically decreasing genetic diversity. In contrast, fluctuating environments enrich genotypes and partially contribute to the maintenance of genetic diversity through balancing selection (Abdul-Rahman et al., 2021). Additionally, parasites show elevated diversity after infecting diverse host populations compared with monocultures (Ekroth et al., 2021), suggesting that host genotype shapes parasite diversity. To adapt to a highly polyandrous honey bees, Nosema parasites may be under the pressure to maintain high polymorphism. In Denmark, honey bees tolerant toward N. ceranae infection showed positive selection (Huang et al., 2014).
N. ceranae is a notorious parasite that weakens honey bee health. In this study, we partially explained the mechanism underlying this parasite’s observed high genetic diversity. We need to point out that multiple proliferation cycles may have occurred for each sampling generation because the infection period was three times more than a typical parasite life cycle. We found that a balance between the novel and lost SNVs drives the maintenance of genetic diversity in this species. Additionally, the allele frequencies of genes essential for parasite infection showed oscillation, supporting the Red Queen hypothesis. Our study enhances the understanding of the genetic structure, facilitated by the global spread of this parasite while providing a blueprint for parasite control.
Statements
Data availability statement
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found below: https://www.ncbi.nlm.nih.gov/, PRJNA784016.
Author contributions
XW conducted the experiment. QH designed the experiment. XW, JZ, JE, and QH organized the manuscript. All authors contributed to the article and approved the submitted version.
Funding
This project was supported by the initiation package of Jiangxi Agricultural University (050014/923230722) and the National Natural Science Foundation of China #32260862.
Acknowledgments
We appreciate the technical support from Zhihao Wu, Bing Liu, and Lizheng Zhang.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
References
1
Abdul-RahmanF.TranchinaD.GreshamD. (2021). Fluctuating environments maintain genetic diversity through neutral fitness effects and balancing selection.Mol. Biol. Evol.384362–4375. 10.1093/molbev/msab173
2
AdrionJ. R.GallowayJ. G.KernA. D. (2020). Predicting the landscape of recombination using deep learning.Mol. Biol. Evol.371790–1808. 10.1093/molbev/msaa038
3
AlexaA.RahnenfuhrerJ. (2021). TopGO: Enrichment Analysis for Gene Ontology.Brighton, MA: Galaxy Integrated.
4
AntunezK.Martin-HernandezR.PrietoL.MeanaA.ZuninoP.HigesM. (2009). Immune suppression in the honey bee (Apis mellifera) following infection by Nosema ceranae (Microsporidia).Environ. Microbiol.112284–2290. 10.1111/j.1462-2920.2009.01953.x
5
BettsA.GrayC.ZelekM.MacLeanR. C.KingK. C. (2018). High parasite diversity accelerates host adaptation and diversification.Science360907–911. 10.1126/science.aam9974
6
BotíasC.Martín-HernándezR.BarriosL.MeanaA.HigesM. (2013). Nosema spp. infection and its negative effects on honey bees (Apis mellifera iberiensis) at the colony level.Vet. Res.44:25. 10.1186/1297-9716-44-25
7
BrockhurstM. A.ChapmanT.KingK. C.MankJ. E.PatersonS.HurstG. D. D. (2014). Running with the Red Queen: the role of biotic conflicts in evolution.Proc. Biol. Sci.281:20141382. 10.1098/rspb.2014.1382
8
BurriL.WilliamsB. A. P.BursacD.LithgowT.KeelingP. J. (2006). Microsporidian mitosomes retain elements of the general mitochondrial targeting system.Proc. Natl. Acad. Sci. U.S.A.10315916–15920. 10.1073/pnas.0604109103
9
Capella-GutiérrezS.Marcet-HoubenM.GabaldónT. (2012). Phylogenomics supports microsporidia as the earliest diverging clade of sequenced fungi.BMC Biol.10:47. 10.1186/1741-7007-10-47
10
ChenS.ZhouY.ChenY.GuJ. (2018). fastp: an ultra-fast all-in-one FASTQ preprocessor.Bioinformatics34i884–i890. 10.1093/bioinformatics/bty560
11
ChenY. P.PettisJ. S.ZhaoY.LiuX.TallonL. J.SadzewiczL. D.et al (2013). Genome sequencing and comparative genomics of honey bee microsporidia, Nosema apis reveal novel insights into host-parasite interactions.BMC Genomics14:451. 10.1186/1471-2164-14-451
12
ChoiJ.KimS.-H. (2017). A genome tree of life for the fungi kingdom.Proc. Natl. Acad. Sci. U.S.A.1149391–9396. 10.1073/pnas.1711939114
13
CuomoC. A.DesjardinsC. A.BakowskiM. A.GoldbergJ.MaA. T.BecnelJ. J.et al (2012). Microsporidian genome analysis reveals evolutionary strategies for obligate intracellular growth.Genome Res.222478–2488. 10.1101/gr.142802.112
14
del AguilaC.IzquierdoF.GranjaA. G.HurtadoC.FenoyS.FresnoM.et al (2006). Encephalitozoon microsporidia modulates p53-mediated apoptosis in infected cells.Int. J. Parasitol.36869–876. 10.1016/j.ijpara.2006.04.002
15
DesjardinsC. A.SanscrainteN. D.GoldbergJ. M.HeimanD.YoungS.ZengQ.et al (2015). Contrasting host-pathogen interactions and genome evolution in two generalist and specialist microsporidian pathogens of mosquitoes.Nat. Commun.6:7121. 10.1038/ncomms8121
16
DongZ.ZhengN.HuC.DengB.FangW.WuQ.et al (2021). Nosema bombycis microRNA-like RNA 8 (Nb-milR8) increases fungal pathogenicity by modulating BmPEX16 gene expression in its host, Bombyx mori.Microbiol. Spectr.9:e0104821. 10.1128/Spectrum.01048-21
17
EiriD. M.SuwannapongG.EndlerM.NiehJ. C. (2015). Nosema ceranae can infect honey bee larvae and reduces subsequent adult longevity.PLoS One10:e0126330. 10.1371/journal.pone.0126330
18
EkrothA. K. E.GerthM.StevensE. J.FordS. A.KingK. C. (2021). Host genotype and genetic diversity shape the evolution of a novel bacterial infection.ISME J.152146–2157. 10.1038/s41396-021-00911-3
19
EllegrenH.GaltierN. (2016). Determinants of genetic diversity.Nat. Rev. Genet.17422–433. 10.1038/nrg.2016.58
20
FriesI.ChauzatM.-P.ChenY.-P.DoubletV.GenerschE.GisderS.et al (2013). Standard methods for Nosema research.J. Apic. Res.521–28. 10.3896/IBRA.1.52.1.14
21
FriesI.FengF.da SilvaA.SlemendaS. B.PieniazekN. J. (1996). Nosema ceranae n. sp. (Microspora, Nosematidae), morphological and molecular characterization of a microsporidian parasite of the Asian honey bee Apis cerana (Hymenoptera, Apidae).Eur. J. Protistol.32356–365. 10.1016/S0932-4739(96)80059-9
22
FuY. X.LiW. H. (1993). Statistical tests of neutrality of mutations.Genetics133693–709. 10.1093/genetics/133.3.693
23
GageS. L.KramerC.CalleS.CarrollM.HeienM.DeGrandi-HoffmanG. (2018). Nosema ceranae parasitism impacts olfactory learning and memory and neurochemistry in honey bees (Apis mellifera).J. Exp. Biol.221jeb161489. 10.1242/jeb.161489
24
Gómez-MorachoT.BartoloméC.BelloX.Martín-HernándezR.HigesM.MasideX. (2015). Recent worldwide expansion of Nosema ceranae (Microsporidia) in Apis mellifera populations inferred from multilocus patterns of genetic variation.Infect. Genet. Evol.3187–94. 10.1016/j.meegid.2015.01.002
25
Gómez-MorachoT.MasideX.Martin-HernandezR.HigesM.BartolomeC. (2014). High levels of genetic diversity in Nosema ceranae within Apis mellifera colonies.Parasitology141475–481. 10.1017/S0031182013001790
26
GraystockP.GoulsonD.HughesW. O. H. (2015). Parasites in bloom: flowers aid dispersal and transmission of pollinator parasites within and between bee species.Proc. Biol. Sci.282:20151371. 10.1098/rspb.2015.1371
27
HeinzelmannR.RiglingD.SiposG.MünsterkötterM.CrollD. (2020). Chromosomal assembly and analyses of genome-wide recombination rates in the forest pathogenic fungus Armillaria ostoyae.Heredity124699–713. 10.1038/s41437-020-0306-z
28
HigesM.Garcia-PalenciaP.Martin-HernandezR.MeanaA. (2007). Experimental infection of Apis mellifera honeybees with Nosema ceranae (Microsporidia).J. Invertebr. Pathol.94211–217. 10.1016/j.jip.2006.11.001
29
HigesM.JuarranzA.Dias-AlmeidaJ.LucenaS.BotiasC.MeanaA.et al (2013). Apoptosis in the pathogenesis of Nosema ceranae (Microsporidia: Nosematidae) in honey bees (Apis mellifera).Environ. Microbiol. Rep.5530–536. 10.1111/1758-2229.12059
30
HigesM.MartinR.MeanaA. (2006). Nosema ceranae, a new microsporidian parasite in honeybees in Europe.J. Invertebr. Pathol.9293–95. 10.1016/j.jip.2006.02.005
31
HigesM.Martín-HernándezR.BotíasC.BailónE. G.González-PortoA. V.BarriosL.et al (2008). How natural infection by Nosema ceranae causes honeybee colony collapse.Environ. Microbiol.102659–2669. 10.1111/j.1462-2920.2008.01687.x
32
HuangQ.ChenY. P.WangR. W.ChengS.EvansJ. D. (2016). Host-parasite interactions and purifying selection in a microsporidian parasite of honey bees.PLoS One11:e0147549. 10.1371/journal.pone.0147549
33
HuangQ.LattorffH. M. G.KrygerP.Le ConteY.MoritzR. F. A. (2014). A selective sweep in a microsporidian parasite Nosema-tolerant honeybee population, Apis mellifera.Anim. Genet.45267–273. 10.1111/age.12114
34
HuangQ.WuZ. H.LiW. F.GuoR.XuJ. S.DangX. Q.et al (2021). Genome and evolutionary analysis of Nosema ceranae: a microsporidian parasite of honey bees.Front. Microbiol.12:645353. 10.3389/fmicb.2021.645353
35
HuangW.-F.JiangJ.-H.ChenY.-W.WangC.-H. (2007). A Nosema ceranae isolate from the honeybee Apis mellifera.Apidologie3830–37. 10.1051/apido:2006054
36
JacksonB. C.CamposJ. L.ZengK. (2015). The effects of purifying selection on patterns of genetic differentiation between Drosophila melanogaster populations.Heredity114163–174. 10.1038/hdy.2014.80
37
KamiyaT.O’DwyerK.NakagawaS.PoulinR. (2014). Host diversity drives parasite diversity: meta-analytical insights into patterns and causal mechanisms.Ecography37689–697. 10.1111/j.1600-0587.2013.00571.x
38
KeL.YanW. Y.ZhangL. Z.ZengZ. J.EvansJ. D.HuangQ. (2021). Honey bee habitat sharing enhances gene flow of the parasite Nosema ceranae.Microb. Ecol.831105–1111. 10.1007/s00248-021-01827-3
39
KoflerR.Orozco-terWengelP.De MaioN.PandeyR. V.NolteV.FutschikA.et al (2011a). PoPoolation: a toolbox for population genetic analysis of next generation sequencing data from pooled individuals.PLoS One6:e15925. 10.1371/journal.pone.0015925
40
KoflerR.PandeyR. V.SchlöttererC. (2011b). PoPoolation2: identifying differentiation between populations using sequencing of pooled DNA samples (Pool-Seq).Bioinformatics273435–3436. 10.1093/bioinformatics/btr589
41
KurzeC.Le ConteY.DussaubatC.ErlerS.KrygerP.LewkowskiO.et al (2015). Nosema tolerant honeybees (Apis mellifera) escape parasitic manipulation of apoptosis.PLoS One10:e0140174. 10.1371/journal.pone.0140174
42
LaiY.-T.YeungC. K. L.OmlandK. E.PangE.-L.HaoY.LiaoB.-Y.et al (2019). Standing genetic variation as the predominant source for adaptation of a songbird.Proc. Natl. Acad. Sci. U.S.A.1162152–2157. 10.1073/pnas.1813597116
43
LiH.DurbinR. (2009). Fast and accurate short read alignment with Burrows-Wheeler transform.Bioinformatics251754–1760. 10.1093/bioinformatics/btp324
44
LiowL. H.Van ValenL.StensethN. C. (2011). Red Queen: from populations to taxa and communities.Trends Ecol. Evol.26349–358. 10.1016/j.tree.2011.03.016
45
Martínez-LópezV.RuizC.MuñozI.OrnosaC.HigesM.Martín-HernándezR.et al (2021). Detection of microsporidia in pollinator communities of a mediterranean biodiversity hotspot for wild bees.Microb. Ecol.84638–642. 10.1007/s00248-021-01854-0
46
Martín-HernándezR.HigesM.SagastumeS.JuarranzÁDias-AlmeidaJ.BudgeG. E.et al (2017). Microsporidia infection impacts the host cell’s cycle and reduces host cell apoptosis.PLoS One12:e0170183. 10.1371/journal.pone.0170183
47
MayackC.NaugD. (2009). Energetic stress in the honeybee Apis mellifera from Nosema ceranae infection.J. Invertebr. Pathol.100185–188. 10.1016/j.jip.2008.12.001
48
McNewS. M.BarrowL. N.WilliamsonJ. L.GalenS. C.SkeenH. R.DuBayS. G.et al (2021). Contrasting drivers of diversity in hosts and parasites across the tropical Andes.Proc. Natl. Acad. Sci. U.S.A.118:e2010714118. 10.1073/pnas.2010714118
49
NdikumanaS.PelinA.WilliotA.SandersJ. L.KentM.CorradiN. (2017). Genome analysis of Pseudoloma neurophilia: a microsporidian Parasite of Zebrafish (Danio rerio).J. Eukaryot. Microbiol.6418–30. 10.1111/jeu.12331
50
NelsonC. W.MonclaL. H.HughesA. L. (2015). SNPGenie: estimating evolutionary parameters to detect natural selection using pooled next-generation sequencing data.Bioinformatics313709–3711. 10.1093/bioinformatics/btv449
51
PaldiN.GlickE.OlivaM.ZilberbergY.AubinL.PettisJ.et al (2010). Effective gene silencing in a microsporidian parasite associated with honeybee (Apis mellifera) colony declines.Appl. Environ. Microbiol.765960–5964. 10.1128/AEM.01067-10
52
PapkouA.GuzellaT.YangW.KoepperS.PeesB.SchalkowskiR.et al (2019). The genomic basis of Red Queen dynamics during rapid reciprocal host–pathogen coevolution.Proc. Natl. Acad. Sci. U.S.A.116923–928. 10.1073/pnas.1810402116
53
PelinA.SelmanM.Aris-BrosouS.FarinelliL.CorradiN. (2015). Genome analyses suggest the presence of polyploidy and recent human-driven expansions in eight global populations of the honeybee pathogen Nosema ceranae.Environ. Microbiol.174443–4458. 10.1111/1462-2920.12883
54
ProesmansW.AlbrechtM.GajdaA.NeumannP.PaxtonR. J.PiozM.et al (2021). Pathways for novel epidemiology: plant-pollinator-pathogen networks and global change.Trends Ecol. Evol.36623–636. 10.1016/j.tree.2021.03.006
55
PtaszyńskaA. A.LatochP.HurdP. J.PolaszekA.Michalska-MadejJ.GrochowalskiŁet al (2021). Amplicon sequencing of variable 16S rRNA from Bacteria and ITS2 regions from fungi and plants, reveals honeybee susceptibility to diseases results from their forage availability under anthropogenic landscapes.Pathogens10:381. 10.3390/pathogens10030381
56
R Core Team (2013). R: A Language and Environment for Statistical Computing.Vienna: R Foundation for Statistical Computing.
57
RabajanteJ. F.TubayJ. M.ItoH.UeharaT.KakishimaS.MoritaS.et al (2016). Host-parasite Red Queen dynamics with phase-locked rare genotypes.Sci. Adv.2:e1501548. 10.1126/sciadv.1501548
58
ScanlonM.LeitchG. J.ShawA. P.MouraH.VisvesvaraG. S. (1999). Susceptibility to apoptosis is reduced in the Microsporidia-infected host cell.J. Eukaryot. Microbiol.4634–35.
59
SchulteR. D.MakusC.SchulenburgH. (2013). Host-parasite coevolution favours parasite genetic diversity and horizontal gene transfer.J. Evol. Biol.261836–1840. 10.1111/jeb.12174
60
StajichJ. E.HahnM. W. (2005). Disentangling the effects of demography and selection in human history.Mol. Biol. Evol.2263–73. 10.1093/molbev/msh252
61
Van der AuweraG. A.CarneiroM. O.HartlC.PoplinR.Del AngelG.Levy-MoonshineA.et al (2013). From FastQ data to high confidence variant calls: the genome analysis toolkit best practices pipeline.Curr. Protoc. Bioinforma.431–33. 10.1002/0471250953.bi1110s43
62
WykaS.MondoS.LiuM.NalamV.BrodersK. (2022). A large accessory genome and high recombination rates may influence global distribution and broad host range of the fungal plant pathogen Claviceps purpurea.PLoS One17:e0263496. 10.1371/journal.pone.0263496
Summary
Keywords
microsporidian, SNV, selection, honey bee, mutation
Citation
Wei X, Zheng J, Evans JD and Huang Q (2022) Transgenerational genomic analyses reveal allelic oscillation and purifying selection in a gut parasite Nosema ceranae. Front. Microbiol. 13:927892. doi: 10.3389/fmicb.2022.927892
Received
25 April 2022
Accepted
11 October 2022
Published
01 November 2022
Volume
13 - 2022
Edited by
Andrew Paul Jackson, University of Liverpool, United Kingdom
Reviewed by
Christopher Mayack, Sabancı University, Turkey; Rui Guo, Fujian Agriculture and Forestry University, China
Updates
Copyright
© 2022 Wei, Zheng, Evans and Huang.
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Qiang Huang, qiang-huang@live.com
This article was submitted to Evolutionary and Genomic Microbiology, a section of the journal Frontiers in Microbiology
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.