Microsporidia are a model clade for studying intracellular parasitism, being well-known for their streamlined genomes and their extreme life history. Although microsporidia are highly diverse and ecologically important to a broad range of hosts, previous research on genome architecture has focused primarily on the mammal-infecting genus Encephalitozoon. Here, we expand that work, testing the universality of the patterns observed in Encephalitozoon by investigating and comparing variation in genetic and epigenetic architectures in the high-quality genome assemblies of several major microsporidia clades. Our comparison of nine genomes, including the first genome assemblies of Binucleata daphniae, Gurleya vavrai, and Conglomerata obtusa, and revised, improved assemblies of Glugoides intestinalis, Mitosporidium daphniae, and Ordospora colligata, found limited conservation of genetic and epigenetic architecture across all microsporidia, although many genomic characteristics, such as nucleotide composition and repeat content, were shared between genomes of the same or related clades. For example, rRNA genes were hypermethylated in most species, but their position close to chromosome ends was only found in the Encephalitozoon and its sister clade. GC-content varied widely, linked to genome size, phylogenetic position, and activity of repeat elements. These findings enhance our insight into genome evolution and, consistent with findings from other systems, suggest epigenetic modification as a regulatory mechanism of gene expression and repeat element activity in microsporidia. Our comparative genome analysis reveals high variation in genetic and epigenetic architecture among microsporidia, despite all of them adapting to a parasitic lifestyle within host cells.

Crabs are a large subtaxon of the Arthropoda, the most diverse and species-rich metazoan group. Several outstanding questions remain regarding crab diversification, including about the genomic capacitors of physiological and morphological adaptation, that cannot be answered with available genomic resources. Physiologically and ecologically diverse Anomuran porcelain crabs offer a valuable model for investigating these questions and hence genomic resources of these crabs would be particularly useful. Here, we present the first two genome assemblies of congeneric and sympatric Anomuran porcelain crabs, Petrolisthes cinctipes and Petrolisthes manimaculis from different microhabitats. Pacific Biosciences high-fidelity sequencing led to genome assemblies of 1.5 and 0.9 Gb, with N50s of 706.7 and 218.9 Kb, respectively. Their assembly length difference can largely be attributed to the different levels of interspersed repeats in their assemblies: The larger genome of P. cinctipes has more repeats (1.12 Gb) than the smaller genome of P. manimaculis (0.54 Gb). For obtaining high-quality annotations of 44,543 and 40,315 protein-coding genes in P. cinctipes and P. manimaculis, respectively, we used RNA-seq as part of a larger annotation pipeline. Contrarily to the large-scale differences in repeat content, divergence levels between the two species as estimated from orthologous protein-coding genes are moderate. These two high-quality genome assemblies allow future studies to examine the role of environmental regulation of gene expression in the two focal species to better understand physiological response to climate change, and provide the foundation for studies in fine-scale genome evolution and diversification of crabs.

This protocol was designed for DNA extraction of about 50 adult female Daphnia magna, it should also work for 1-150 animals but with adjusted reagent volumes. For achieving HMW DNA or maximizing yield, some modifications are indicated as substeps. Before DNA extraction, animals can be freed from microbes using antibiotics (https://www.evolution.unibas.ch/ebert/lab/daphnia_dna.htm) and should be dehydrated and snap-frozen with liquid nitrogen for archiving HMW DNA.

Two important characteristics of metapopulations are extinction-(re)colonization dynamics and gene flow between subpopulations. These processes can cause strong shifts in genome-wide allele frequencies that are generally not observed in "classical" (large, stable, and panmictic) populations. Subpopulations founded by one or a few individuals, the so-called propagule model, are initially expected to show intermediate allele frequencies at polymorphic sites until natural selection and genetic drift drive allele frequencies toward a mutation-selection-drift equilibrium characterized by a negative exponential-like distribution of the site frequency spectrum. We followed changes in site frequency spectrum distribution in a natural metapopulation of the cyclically parthenogenetic pond-dwelling microcrustacean Daphnia magna using biannual pool-seq samples collected over a 5-yr period from 118 ponds occupied by subpopulations of known age. As expected under the propagule model, site frequency spectra in newly founded subpopulations trended toward intermediate allele frequencies and shifted toward right-skewed distributions as the populations aged. Immigration and subsequent hybrid vigor altered this dynamic. We show that the analysis of site frequency spectrum dynamics is a powerful approach to understand evolution in metapopulations. It allowed us to disentangle evolutionary processes occurring in a natural metapopulation, where many subpopulations evolve in parallel. Thereby, stochastic processes like founder and immigration events lead to a pattern of subpopulation divergence, while genetic drift leads to converging site frequency spectrum distributions in the persisting subpopulations. The observed processes are well explained by the propagule model and highlight that metapopulations evolve differently from classical populations.

Genomics is a powerful toolkit for unravelling how evolutionary processes drive organisms’ small- and large-scale genetic variation. Several outstanding questions remain concerning the evolution of genome size and architecture, especially in intracellular parasites. Microsporidia became a model for this field of study as they exhibit genome size variation of more than an order of magnitude. Here, we discuss evolution in the large-genome microsporidium Hamiltosporidium tvaerminnensis, a parasite of a water flea.

Microsporidia are intracellular parasitic fungi whose genomes rank among the smallest of all known eukaryotes. A number of outstanding questions remain concerning the evolution of their large-scale variation in genome architecture, responsible for genome size variation of more than an order of magnitude. This genome report presents the first near-chromosomal assembly of a large-genome microsporidium, Hamiltosporidium tvaerminnensis. Combined Oxford Nanopore, Pacific Biosciences (PacBio), and Illumina sequencing led to a genome assembly of 17 contigs, 11 of which represent complete chromosomes. Our assembly is 21.64 Mb in length, has an N50 of 1.44 Mb, and consists of 39.56% interspersed repeats. We introduce a novel approach in microsporidia, PacBio Iso-Seq, as part of a larger annotation pipeline for obtaining high-quality annotations of 3,573 protein-coding genes. Based on direct evidence from the full-length Iso-Seq transcripts, we present evidence for alternative polyadenylation and variation in splicing efficiency, which are potential regulation mechanisms for gene expression in microsporidia. The generated high-quality genome assembly is a necessary resource for comparative genomics that will help elucidate the evolution of genome architecture in response to intracellular parasitism.

The determinants of variation in a species' genome-wide nucleotide diversity include historical, environmental, and stochastic aspects. This diversity can inform us about the species' past and present evolutionary dynamics. In parasites, the mode of transmission and the interactions with the host might supersede the effects of these aspects in shaping parasite genomic diversity. We used genomic samples from 10 populations of the microsporidian parasite Ordospora colligata to investigate present genomic diversity and how it was shaped by evolutionary processes, specifically, the role of phylogeography, co-phylogeography (with the host), natural selection, and transmission mode. Although very closely related microsporidia cause diseases in humans, O. colligata is specific to the freshwater crustacean Daphnia magna and has one of the smallest known eukaryotic genomes. We found an overlapping phylogeography between O. colligata and its host highlighting the long-term, intimate relationship between them. The observed geographic distribution reflects previous findings that O. colligata exhibits adaptations to colder habitats, which differentiates it from other microsporidian gut parasites of D. magna predominantly found in warmer areas. The co-phylogeography allowed us to calibrate the O. colligata phylogeny and thus estimate its mutation rate. We identified several genetic regions under potential selection. Our whole-genome study provides insights into the evolution of one of the most reduced eukaryotic genomes and shows how different processes shape genomic diversity of an obligate parasite.

The dynamics of extinction and (re)colonization in habitat patches are characterizing features of dynamic metapopulations, causing them to evolve differently than large, stable populations. The propagule model, which assumes genetic bottlenecks during colonization, posits that newly founded subpopulations have low genetic diversity and are genetically highly differentiated from each other. Immigration may then increase diversity and decrease differentiation between subpopulations. Thus, older and/or less isolated subpopulations are expected to have higher genetic diversity and less genetic differentiation. We tested this theory using whole-genome pool-sequencing to characterize nucleotide diversity and differentiation in 60 subpopulations of a natural metapopulation of the cyclical parthenogen Daphnia magna. For comparison, we characterized diversity in a single, large, and stable D. magna population. We found reduced (synonymous) genomic diversity, a proxy for effective population size, weak purifying selection, and low rates of adaptive evolution in the metapopulation compared with the large, stable population. These differences suggest that genetic bottlenecks during colonization reduce effective population sizes, which leads to strong genetic drift and reduced selection efficacy in the metapopulation. Consistent with the propagule model, we found lower diversity and increased differentiation in younger and also in more isolated subpopulations. Our study sheds light on the genomic consequences of extinction-(re)colonization dynamics to an unprecedented degree, giving strong support for the propagule model. We demonstrate that the metapopulation evolves differently from a large, stable population and that evolution is largely driven by genetic drift.

Analysing variation in a species' genomic diversity can provide insights into its historical demography, biogeography and population structure, and thus its ecology and evolution. Although such studies are rarely undertaken for parasites, they can be highly revealing because of the parasite's co-evolutionary relationships with hosts. Modes of reproduction and transmission are thought to be strong determinants of genomic diversity for parasites and vary widely among microsporidia (fungal-related intracellular parasites), which are known to have high intraspecific genetic diversity and interspecific variation in genome architecture. Here we explore genomic variation in the microsporidium Hamiltosporidium, a parasite of the freshwater crustacean Daphnia magna, looking especially at which factors contribute to nucleotide variation. Genomic samples from 18 Eurasian populations and a new, long-read-based reference genome were used to determine the roles that reproduction mode, transmission mode and geography play in determining population structure and demographic history. We demonstrate two main Hamiltosporidium tvaerminnensis lineages and a pattern of isolation-by-distance, but note an absence of congruence between these two parasite lineages and the two Eurasian host lineages. We suggest a comparatively recent parasite spread through Northern Eurasian host populations after a change from vertical to mixed-mode transmission and the loss of sexual reproduction. While gaining knowledge about the ecology and evolution of this focal parasite, we also identify common features that shape variation in genomic diversity for many parasites, such as distinct modes of reproduction and the intertwining of host-parasite demographies.

The determinants of variation in a species’ genome-wide nucleotide diversity include historical, environmental, and stochastic aspects. This diversity can inform us about the species’ past and present evolutionary dynamics. In parasites, the mode of transmission and the interactions with the host might supersede the effects of these aspects in shaping parasite genomic diversity. We used genomic samples from ten populations of the microsporidian parasite Ordospora colligata to investigate present genomic diversity and how it was shaped by evolutionary processes, specifically, the role of phylogeography, co-phylogeography (with the host), natural selection, and transmission mode. Although very closely related microsporidia cause diseases in humans, O. colligata is specific to the freshwater crustacean Daphnia magna and has one of the smallest known eukaryotic genomes. We found an overlapping phylogeography between O. colligata and its host highlighting the long-term, intimate relationship between them. The observed geographic distribution reflects previous findings that O. colligata exhibits adaptations to colder habitats, which differentiates it from other microsporidian gut parasites of D. magna predominantly found in warmer areas. The co-phylogeography allowed us to calibrate the O. colligata phylogeny and thus estimate its mutation rate. We found patterns of more efficient purifying selection in O. colligata relative to other microsporidia sharing the same host, which likely allowed this parasite to maintain its very compact genome. We also identified regions under potential selection related to coevolution including the ribosomal protein L24, a leucyl-tRNA synthetase, and a putative ABC-like lipid transport protein. Our whole-genome study provides insights into the evolution of one of the most reduced eukaryotic genomes and shows how different processes shape genomic diversity of an obligate parasite. Author summary Microsporidia are intracellular parasites that infect vertebrates, invertebrates, and even unicellular organisms. Due to their high variation in many aspects of life history and genomics, microsporidia have become a model clade for understanding evolutionary processes related to intracellular parasitism. However, the evolution of extreme genomic architectures in microsporidia and the coevolution with their hosts is still under-surveyed, especially given their role in human disease. Here, we study past and present evolutionary dynamics in a microsporidian species with one of the smallest known eukaryotic genomes, O. colligata . Close relatives of O. colligata cause death and disease in humans and agriculturally important animals. We show that purifying selection helped maintaining its reduced, compact genome and corroborate hypotheses about the evolution of different genome sizes in microsporidia. Importantly, we utilize the highly resolved phylogeny of its host to estimate the parasite’s mutation rate. This methodology allowed us to establish the first mutation rate estimate for a microsporidium, an estimate which is within the range of mutation rates estimated for phylogenetically related, non-parasitic fungi. Our study exemplifies how the combined knowledge about a species’ biology, ecology, and genomic diversity helps to resolve its evolutionary dynamics, in particular when phylogenomic information can be brought to bear for both host and parasite.