Shared geographic histories and dispersal contribute to congruent phylogenies between amphipods and their microsporidian parasites at regional and global scales

O short-branch Microsporidia, where art thou? Identifying diversity hotspots for future sampling

Short-branch Microsporidia were previously shown to form a basal grade within the expanded Microsporidia clade and to branch near the classical, long-branch Microsporidia. Although they share simpler versions of some morphological characteristics, they do not show accelerated evolutionary rates, making them ideal candidates to study the evolutionary trajectories that have led to long-branch microsporidian unique characteristics. However, most sequences assigned to the short-branch Microsporidia are undescribed, novel environmental lineages for which the identification requires knowledge of where they can be found. To direct future isolation, we used the EukBank database of the global UniEuk initiative that contains the majority of the publicly available environmental V4 SSU rRNA gene sequences of protists. The curated OTU table and corresponding metadata were used to evaluate the occurrence of short-branch Microsporidia across freshwater, hypersaline, marine benthic, marine pelagic, and terrestrial environments. Presence-absence analyses infer that short-branch Microsporidia are most abundant in freshwater and terrestrial environments, and alpha- and beta-diversity measures indicate that focusing our sampling effort on these two environments would cover a large part of their overall diversity. These results can be used to coordinate future isolation and sampling campaigns to better understand the enigmatic evolution of microsporidians' unique characteristics.

Testing the radiation cascade in postglacial radiations of whitefish and their parasites: founder events and host ecology drive parasite evolution

Reciprocal effects of adaptive radiations on the evolution of interspecific interactions, like parasitism, remain barely explored. We test whether the recent radiations of European whitefish (Coregonus spp.) across and within perialpine and subarctic lakes promote its parasite Proteocephalus fallax (Platyhelminthes: Cestoda) to undergo host repertoire expansion via opportunity and ecological fitting, or adaptive radiation by specialization. Using de novo genomic data, we examined P. fallax differentiation across lakes, within lakes across sympatric host species, and the contributions of host genetics versus host habitat use and trophic preferences. Whitefish intralake radiations prompted parasite host repertoire expansion in all lakes, whereas P. fallax differentiation remains incipient among sympatric fish hosts. Whitefish genetic differentiation per se did not explain the genetic differentiation among its parasite populations, ruling out codivergence with the host. Instead, incipient parasite differentiation was driven by whitefish phenotypic radiation in trophic preferences and habitat use in an arena of parasite opportunity and ecological fitting to utilize resources from emerging hosts. Whilst the whitefish radiation provides a substrate for the parasite to differentiate along the same water-depth ecological axis as Coregonus spp., the role of the intermediate hosts in parasite speciation may be overlooked. Parasite multiple-level ecological fitting to both fish and crustacean intermediate hosts resources may be responsible for parasite population substructure in Coregonus spp. We propose parasites' delayed arrival was key to the initial burst of postglacial intralake whitefish diversification, followed by opportunistic tapeworm host repertoire expansion and a delayed nonadaptive radiation cascade of incipient tapeworm differentiation. At the geographical scale, dispersal, founder events, and genetic drift following colonization of spatially heterogeneous landscapes drove strong parasite differentiation. We argue that these microevolutionary processes result in the mirroring of host-parasite phylogenies through phylogenetic tracking at macroevolutionary and geographical scales.

Assessing global drivers of parasite diversity: host diversity and body mass boost avian haemosporidian diversity

Biodiversity varies worldwide and is influenced by multiple factors, such as environmental stability and past historical events (e.g. Panama Isthmus). At the same time, organisms with unique life histories (e.g. parasites) are subject to unique selective pressures that structure their diversity patterns. Parasites represent one of the most successful life strategies, impacting, directly and indirectly, ecosystems by cascading effects on host fitness and survival. Here, I focused on a highly diverse, prevalent and cosmopolitan group of parasites (avian haemosporidians) to investigate the main drivers (e.g. host and environmental features) of regional parasite diversity on a global scale. To do so, I compiled data from 4 global datasets on (i) avian haemosporidian (malaria and malaria-like) parasites, (ii) bird species diversity, (iii) avian functional traits and (iv) climate data. Then, using generalized least square models, I evaluated the effect of host and environmental features on haemosporidian diversity. I found that haemosporidian diversity mirrors host regional diversity and that higher host body mass increases haemosporidian diversity. On the other hand, climatic conditions had no effect on haemosporidian diversity in any model. When evaluating Leucocytozoon parasites separately, I found parasite diversity was boosted by a higher proportion of migratory hosts. In conclusion, I demonstrated that haemosporidian parasite diversity is intrinsically associated with their hosts’ diversity and body mass.

Host phylogeny and seasonality shapes avian haemosporidian prevalence in a Brazilian biodiverse and dry forest: the Caatinga

The relationships between host phylogenetics, functional traits and parasites in wildlife remain poorly understood in the Neotropics, especially in habitats with marked seasonal variation. Here, we examined the effect of seasonality and host functional traits on the prevalence of avian haemosporidians (Plasmodium and Haemoproteus) in the Brazilian Caatinga, a seasonally dry tropical forest. 933 birds were evaluated for haemosporidian infections. We found a high parasitism prevalence (51.2%), which was correlated with phylogenetic relatedness among avian species. Prevalence varied drastically among the 20 well-sampled species, ranging from 0 to 70%. Seasonality was the main factor associated with infections, but how this abiotic condition influenced parasite prevalence varied according to the host-parasite system. Plasmodium prevalence increased during the rainy season and, after excluding the large sample size of Columbiformes (n = 462/933), Plasmodium infection rate was maintained high in the wet season and showed a negative association with host body mass. No association was found between non-Columbiform bird prevalence and seasonality or body mass when evaluating both Plasmodium and Haemoproteus or only Haemoproteus infections. Parasite community was composed of 32 lineages including 7 new lineages. We evidenced that even dry domains can harbour a high prevalence and diversity of vector-borne parasites and pointed out seasonality as a ruling factor.

Effect of Host-Switching on the Ecological and Evolutionary Patterns of Parasites

Speciation via host-switching is a macroevolutionary process that emerges from a microevolutionary dynamic where individual parasites switch hosts, establish a new association, and reduce reproductive contact with the original parasite lineage. Phylogenetic distance and geographic distribution of the hosts have been shown to be determinants of the capacity and opportunity of the parasite to change hosts. Although speciation via host-switching has been reported in many host-parasite systems, its dynamic on the individual, population and community levels is poorly understood. Here we propose a theoretical model to simulate parasite evolution considering host-switching events on the microevolutionary scale, taking into account the macroevolutionary history of the hosts, to evaluate how host-switching can affect ecological and evolutionary patterns of parasites in empirical communities at regional and local scales. In the model, parasite individuals can switch hosts under variable intensity and have their evolution driven by mutation and genetic drift. Mating is sexual and only individuals that are sufficiently similar can produce offspring. We assumed that parasite evolution occurs at the same evolutionary time scale as their hosts, and that the intensity of host-switching decreases as the host species differentiate. Ecological and evolutionary patterns were characterized by the turnover of parasite species among host species, and parasite evolutionary tree imbalance respectively. We found a range of host-switching intensity that reproduces ecological and evolutionary patterns observed in empirical communities. Our results showed that turnover decreased as host-switching intensity increased, with low variation among the model replications. On the other hand, tree imbalance showed wide variation and non-monotonic tendency. We concluded that tree imbalance was sensitive to stochastic events, whereas turnover may be a good indicator of host-switching. We found that local communities corresponded to higher host-switching intensity when compared to regional communities, highlighting that spatial scale is a limitation for host-switching. [Dispersal of parasites, opportunity and capacity of interaction, phylogenetic conservatism, and community structure.].

Widespread infection, diversification and old host associations of Nosema Microsporidia in European freshwater gammarids (Amphipoda)

The microsporidian genus Nosema is primarily known to infect insects of economic importance stimulating high research interest, while other hosts remain understudied. Nosema granulosis is one of the formally described Nosema species infecting amphipod crustaceans, being known to infect only two host species. Our first aim was to characterize Nosema spp. infections in different amphipod species from various European localities using the small subunit ribosomal DNA (SSU) marker. Second, we aimed to assess the phylogenetic diversity, host specificity and to explore the evolutionary history that may explain the diversity of gammarid-infecting Nosema lineages by performing a phylogenetic reconstruction based on RNA polymerase II subunit B1 (RPB1) gene sequences. For the host species Gammarus balcanicus, we also analyzed whether parasites were in excess in females to test for sex ratio distortion in relation with Nosema infection. We identified Nosema spp. in 316 individuals from nine amphipod species being widespread in Europe. The RPB1-based phylogenetic reconstruction using newly reported sequences and available data from other invertebrates identified 39 haplogroups being associated with amphipods. These haplogroups clustered into five clades (A-E) that did not form a single amphipod-infecting monophyletic group. Closely related sister clades C and D correspond to Nosema granulosis. Clades A, B and E might represent unknown Nosema species infecting amphipods. Host specificity seemed to be variable with some clades being restricted to single hosts, and some that could be found in several host species. We show that Nosema parasite richness in gammarid hosts is much higher than expected, illustrating the advantage of the use of RPB1 marker over SSU. Finally, we found no hint of sex ratio distortion in Nosema clade A infecting G. balcanicus. This study shows that Nosema spp. are abundant, widespread and diverse in European gammarids. Thus, Nosema is as diverse in aquatic as in terrestrial hosts.

Microsporidian diversity in the aquatic isopod Asellus aquaticus

We conducted a molecular survey on microsporidian diversity in different lineages (operational taxonomic units = OTUs) of Asellus aquaticus from 30 sites throughout Europe. Host body length was determined, and DNA was extracted from host tissue excluding the intestine and amplified by microsporidian-specific primers. In total, 247 A. aquaticus specimens were analysed from which 26.7% were PCR-positive for microsporidians, with significantly more infections in larger individuals. Prevalence ranged between 10 and 90%. At 9 sites, no microsporidians were detected. A significant relationship was found between the frequency of infected individuals and habitat type, as well as host OTU. The lowest proportion of infected individuals was detected in spring-habitats (8.7%, n = 46) and the highest in ponds (37.7%, n = 53). Proportion of infected individuals among host OTUs A, D and J was 31.7, 21.7 and 32.1%, respectively. No infections were detected in OTU F. Our results are, however, accompanied by a partially low sample size, as only a minimum of 5 individuals was available at a few locations. Overall, 17 different microsporidian molecular taxonomic units (MICMOTUs) were distinguished with 5 abundant isolates (found in 4–17 host individuals) while the remaining 12 MICMOTUs were “rare” and found only in 1–3 host individuals. No obvious spatio-genetic pattern could be observed. The MICMOTUs predominantly belonged to Nosematida and Enterocytozoonida. The present study shows that microsporidians in A. aquaticus are abundant and diverse but do not show obvious patterns related to host genetic lineages or geography.

Extremely divergent COI sequences within an amphipod species complex: A possible role for endosymbionts?

Some heritable endosymbionts can affect host mtDNA evolution in various ways. Amphipods host diverse endosymbionts, but whether their mtDNA has been influenced by these endosymbionts has yet to be considered. Here, we investigated the role of endosymbionts (microsporidians and Rickettsia) in explaining highly divergent COI sequences in Paracalliope fluviatilis species complex, the most common freshwater amphipods in New Zealand. We first contrasted phylogeographic patterns using COI, ITS, and 28S sequences. While molecular species delimitation methods based on 28S sequences supported 3-4 potential species (N, C, SA, and SB) among freshwater lineages, COI sequences supported 17-27 putative species reflecting high inter-population divergence. The deep divergence between NC and S lineages (~20%; 28S) and the substitution saturation on the 3rd codon position of COI detected even within one lineage (SA) indicate a very high level of morphological stasis. Interestingly, individuals infected and uninfected by Rickettsia comprised divergent COI lineages in one of four populations tested, suggesting a potential influence of endosymbionts in mtDNA patterns. We propose several plausible explanations for divergent COI lineages, although they would need further testing with multiple lines of evidence. Lastly, due to common morphological stasis and the presence of endosymbionts, phylogeographic patterns of amphipods based on mtDNA should be interpreted with caution.

Generalist parasites persist in degraded environments: a lesson learned from microsporidian diversity in amphipods

The present study provides new insight into suitable microsporidian–host associations. It relates regional and continental-wide host specialization in microsporidians infecting amphipods to degraded and recovering habitats across 2 German river catchments. It provides a unique opportunity to infer the persistence of parasites following anthropogenic disturbance and their establishment in restored rivers. Amphipods were collected in 31 sampling sites with differing degradation and restoration gradients. Specimens were morphologically (hosts) and molecularly identified (host and parasites). Amphipod diversity and abundance, microsporidian diversity, host phylogenetic specificity and continental-wide β-specificity were investigated and related to each other and/or environmental variables. Fourteen microsporidian molecular operational taxonomic units (MOTUs), mainly generalist parasites, infecting 6 amphipod MOTUs were detected, expanding the current knowledge on the host range by 17 interactions. There was no difference in microsporidian diversity and host specificity among restored and near-natural streams (Boye) or between those located in urban and rural areas (Kinzig). Similarly, microsporidian diversity was generally not influenced by water parameters. In the Boye catchment, host densities did not influence microsporidian MOTU richness across restored and near-natural sites. High host turnover across the geographical range suggests that neither environmental conditions nor host diversity plays a significant role in the establishment into restored areas. Host diversity and environmental parameters do not indicate the persistence and dispersal of phylogenetic host generalist microsporidians in environments that experienced anthropogenic disturbance. Instead, these might depend on more complex mechanisms such as the production of resistant spores, host switching and host dispersal acting individually or conjointly.

Migratory behaviour does not alter cophylogenetic congruence between avian hosts and their haemosporidian parasites

Parasites display various degrees of host specificity, reflecting different coevolutionary histories with their hosts. Avian hosts follow multiple migration patterns representing short but also long distances. As parasites infecting migratory birds are subjected to multiple environmental and biotic changes through their flyways, migration may disrupt or strengthen cophylogenetic congruence between hosts and parasites. On the one hand, parasites might adapt to a single migratory host, evolving to cope with the specific challenges associated with the multiple habitats occupied by the host. On the other, as migrants can introduce parasites into new habitats, higher rates of host switching could also disrupt cophylogenetic patterns. We analysed whether migratory behaviour shapes avian haemosporidian parasite–host cophylogenetic congruence by testing if contributions of host–parasite links to overall congruence differ among resident and short-, variable- and long-distance migrants globally and within South America only. On both scales, we found significant overall cophylogenetic congruence by testing whether overall congruence differed between haemosporidian lineages and bird species. However, we found no difference in contribution towards congruence among links involving resident vs migratory hosts in both models. Thus, migratory behaviour neither weakens nor strengthens bird–haemosporidian cophylogenetic congruence, suggesting that other avian host traits are more influential in generating phylogenetic congruence in this host–parasite system.

Two parasites in one host: spatiotemporal dynamics and co-occurrence of Microsporidia and Rickettsia in an amphipod host

Biological interactions can greatly influence the abundance of species. This is also true for parasitic species that share the same host. Microsporidia and Rickettsia are widespread intracellular parasites in populations of Paracalliope fluviatilis, the most common freshwater amphipods in New Zealand. Although both parasites coexist in many populations, it is unclear whether they interact with each other. Here, we investigated spatial−temporal dynamics and co-occurrence of the two parasites, Microsporidia and Rickettsia in P. fluviatilis hosts, across one annual cycle and in three different locations. Prevalence of both Microsporidia and Rickettsia changed over time. However, while the prevalence of Rickettsia varied significantly between sampling times, that of Microsporidia did not change significantly and remained relatively low. The two parasites therefore followed different temporal patterns. Also, the prevalence of both parasites differed among locations, though the two species reached their highest prevalence in different locations. Lastly, there was no evidence for positive or negative associations between the two parasite species; the presence of one parasite in an individual host does not appear to influence the probability of infection by the other parasite. Their respective prevalence may follow different patterns among populations on a larger spatial scale due to environmental heterogeneity across locations.

Next-generation cophylogeny: unravelling eco-evolutionary processes

A fundamental question in evolutionary biology is how microevolutionary processes translate into species diversification. Cophylogeny provides an appropriate framework to address this for symbiotic associations, but historically has been primarily limited to unveiling patterns. We argue that it is essential to integrate advances from ecology and evolutionary biology into cophylogeny, to gain greater mechanistic insights and transform cophylogeny into a platform to advance understanding of interspecific interactions and diversification more widely. We discuss key directions, such as incorporating trait reconstruction and considering multiple scales of network organization, and highlight recent developments for implementation. A new quantitative framework is proposed to allow integration of relevant information, such as quantitative traits and assessment of the contribution of individual mechanisms to cophylogenetic patterns.

Revisiting the phylogeny of microsporidia

Canonical microsporidians are a group of obligate intracellular parasites of a wide range of hosts comprising ~1,300 species of >220 genera. Microsporidians are related to fungi, and many characterised and uncharacterized groups closely related to them have been discovered recently, filling the knowledge gaps between them. These groups assigned to the superphylum Opisthosporidia have provided several important insights into the evolution of diverse intracellular parasitic lineages within the tree of eukaryotes. The most studied among opisthosporidians, canonical microsporidians, were known to science more than 160 years ago, however, the classification of canonical Microsporidia has been challenging due to common morphological homoplasy, and accelerated evolutionary rates. Instead of morphological characters, ssrRNA sequences have been used as the primary data for the classification of canonical microsporidians. Previous studies have produced a useful backbone of the microsporidian phylogeny, but provided only some nodal support, causing some confusion. Here, we reconstructed phylogenetic trees of canonical microsporidians using Bayesian and Maximum Likelihood inferences. We included rRNA sequences of 126 described/named genera, by far the broadest taxon coverage to date. Overall, our trees show similar topology and recovered four of the five main clades demonstrated in previous studies (Clades 1, 3, 4 and 5). Family level clades were well resolved within each major clade, but many were discordant with the recently revised classification. Therefore, revision and some reshuffling, especially within and between Clades 1 and 3 are required. We also reconstructed phylogenetic trees of Opisthosporidia to better integrate the evolutionary history of canonical microsporidians in a broader context. We discuss several traits shared only by canonical microsporidians that may have contributed to their striking ecological success in diverse metazoans. More targeted studies on the neglected host groups will be of value for a better understanding of the evolutionary history of these interesting intracellular parasites.

Wide geographic distribution of overlooked parasites: Rare Microsporidia in Gammarus balcanicus, a species complex with a high rate of endemism

Parasites and other symbionts deeply influence host organisms, and no living organism can be considered to have evolved independent of its symbionts. The first step towards understanding symbiotic influences upon host organisms is a strong supporting knowledge of parasite/symbiont diversity. Parasites of freshwater amphipods are diverse, with Microsporidia being a major group. These intracellular parasites impact gammarid fitness in different ways, ranging from reduced fitness to increased fecundity. Many Microsporidia have been recorded using molecular data, with multiple taxa pending formal taxonomic description. While some parasites are common, others are known only through sporadic records of single infections. In this study, we focus on rare/sporadic microsporidian infections within Gammarus balcanicus, a host species complex with a high level of endemism. In addition to enriching our knowledge on Microsporidia parasite diversity in amphipod hosts, we test whether these symbionts are specific to G. balcanicus or if they are the same taxa infecting other gammarid species. Of 2231 hosts from 87 sites, we catalogued 29 sequences of "rare" Microsporidia clustering into 19 haplogroups. These haplogroups cluster into 11 lineages: four pre-described taxa (Cucumispora roeselum, C. ornata, C. dikerogammari and Enterocytospora artemiae) and seven 'Molecular Operational Taxonomic Units', which are known from previously published studies to infect other European amphipod species. Our study significantly widens the geographic range of these Microsporidia and expands the known spectrum of hosts infected. Our results suggest that these parasites are ancient infections of European gammarids. For some host-parasite systems, we hypothesize that the common parasite ancestors that infected the hosts' common ancestors, diversified alongside host diversification. For others, we observe Microsporidia taxa with wide host ranges that do not follow host phylogeny.

Dictyocoela microsporidia diversity and co-diversification with their host, a gammarid species complex (Crustacea, Amphipoda) with an old history of divergence and high endemic diversity

BACKGROUND

Although the processes of co-evolution between parasites and their hosts are well known, evidence of co-speciation remains scarce. Microsporidian intracellular parasites, due to intimate relationships with their hosts and mixed mode of transmission (horizontal but also vertical, from mother to offspring), may represent an interesting biological model for investigating co-speciation. Amphipod crustaceans, especially gammarids, are regular hosts of microsporidian parasites, in particular the Dictyocoela spp., which have so far been found limited to these amphipods and are known to use a vertical mode of transmission. The amphipod genus Gammarus has a diversification history spanning the last 50-60 Mya and an extensive cryptic diversity in most of the nominal species. Here, we investigated the degree of co-diversification between Dictyocoela and Gammarus balcanicus, an amphipod with high degrees of ancient cryptic diversification and lineage endemism, by examining the genetic diversity of these parasites over the entire geographic range of the host. We hypothesised that the strong host diversification and vertical transmission of Dictyocoela would promote co-diversification.

RESULTS

Using the parasite SSU rDNA as a molecular marker, analyzing 2225 host specimens from 88 sites covering whole host range, we found 31 haplogroups of Dictyocoela, 30 of which were novel, belonging to four Dictyocoela species already known to infect other Gammarus spp. The relationships between Dictyocoela and gammarids is therefore ancient, with the speciation in parasites preceding those of the hosts. Each novel haplogroup was nevertheless specific to G. balcanicus, leaving the possibility for subsequent co-diversification process during host diversification. A Procrustean Approach to Co-phylogeny (PACo) analysis revealed that diversification of Dictyocoela was not random with respect to that of the host. We found high degrees of congruence between the diversification of G. balcanicus and that of Dictyocoela roeselum and D. muelleri. However, we also found some incongruences between host and Dictyocoela phylogenies, e.g. in D. duebenum, probably due to host shifts between different G. balcanicus cryptic lineages.

CONCLUSION

The evolutionary history of Dictyocoela and Gammarus balcanicus represents an example of an overall host-parasite co-diversification, including cases of host shifts.

Widespread Torix Rickettsia in New Zealand amphipods and the use of blocking primers to rescue host COI sequences

Endosymbionts and intracellular parasites are common in arthropod hosts. As a consequence, (co)amplification of untargeted bacterial sequences has been occasionally reported as a common problem in DNA barcoding. While identifying amphipod species with universal COI primers, we unexpectedly detected rickettsial endosymbionts belonging to the Torix group. To map the distribution and diversity of Rickettsia species among amphipod hosts, we conducted a nationwide molecular screening of seven families of New Zealand freshwater amphipods. In addition to uncovering a diversity of Torix Rickettsia species across multiple amphipod populations from three different families, our research indicates that: (1) detecting Torix Rickettsia with universal primers is not uncommon, (2) obtaining 'Rickettsia COI sequences' from many host individuals is highly likely when a population is infected, and (3) obtaining 'host COI' may not be possible with a conventional PCR if an individual is infected. Because Rickettsia COI is highly conserved across diverse host taxa, we were able to design blocking primers that can be used in a wide range of host species infected with Torix Rickettsia. We propose the use of blocking primers to circumvent problems caused by unwanted amplification of Rickettsia and to obtain targeted host COI sequences for DNA barcoding, population genetics, and phylogeographic studies.