Tracking the assembly of nested parasite communities: Using β-diversity to understand variation in parasite richness and composition over time and scale

Hierarchical diversity partitioning of microscopic epibiont community on intertidal molluscan shells and inert surfaces over three geographic regions in Japan

Microscopic epibionts on molluscan shells are a component of the biodiversity of intertidal coastal areas. Because molluscan shells are discrete habitats for the epibiont community, and the molluscan basibionts belong to the local community, epibiont diversity can be evaluated hierarchically by basibiont categories including species. To evaluate the structure of epibiont diversity and effects of taxonomic resolution on the evaluation, epibionts on molluscan shells and inert surfaces were investigated at three geographically distant sites in Japan. In total, 94 species-level taxonomic units of epibionts were obtained from 31 basibiont molluscan species and inert surfaces (plastics and rock chips). The density and the species richness at the site of the lowest latitude were significantly lower than those at the other sites. The epibiont community differed between the three sites, although the major portion of the epibionts were diatoms. Between-site diversity contributed most of the total diversity of the species richness and Simpson diversity in the five levels of the hierarchical partitioning: sample (individual basibiont), basibiont species (molluscan species), surface group (bivalves, chitons + limpets, and globose gastropods), site, and the total. The taxonomic resolution did not markedly affect the variability of communities between the three sites, although the taxon richness was reduced to 51 in the genus-level analysis. The lower taxonomic resolution (genus level); however, increased the contribution of the within-sample and decreased the contribution of β diversities at the higher hierarchies, leading to a possible overestimation of biotic homogenization between the communities.

Diverging effects of host density and richness across biological scales drive diversity-disease outcomes

Understanding how biodiversity affects pathogen transmission remains an unresolved question due to the challenges in testing potential mechanisms in natural systems and how these mechanisms vary across biological scales. By quantifying transmission of an entire guild of parasites (larval trematodes) within 902 amphibian host communities, we show that the community-level drivers of infection depend critically on biological scale. At the individual host scale, increases in host richness led to fewer parasites per host for all parasite taxa, with no effect of host or predator densities. At the host community scale, however, the inhibitory effects of richness were counteracted by associated increases in total host density, leading to no overall change in parasite densities. Mechanistically, we find that while average host competence declined with increasing host richness, total community competence remained stable due to additive assembly patterns. These results help reconcile disease-diversity debates by empirically disentangling the roles of alternative ecological drivers of parasite transmission and how such effects depend on biological scale.

Nestedness of flea assemblages harboured by small mammalian hosts revisited: phylogenetic and functional nestedness do not follow compositional nestedness

We studied compositional, phylogenetic, and functional nestedness in the flea assemblages of 14 host species across regions. Our main questions were (a) are a host's flea assemblages compositionally, phylogenetically, or functionally nested? (b) Do similar processes drive these nestedness facets? (d) Are a host's biological traits associated with nestedness of its flea assemblages? Rows of host matrices were ordered by decreasing species richness/the sum of the branch lengths of a phylogenetic tree/functional dendrogram or by decreasing region area or by increasing distance from the centre of a host's geographic range. None of the matrices sorted by species richness/sum of branch lengths were nested from a compositional perspective, but they were significantly nested from phylogenetic and functional perspectives. Compositional, phylogenetic, and functional nestedness of matrices sorted by region area or by distance from the host's geographic range centre varied between hosts. In some hosts, flea assemblages were nested from all three perspectives independently of how matrix rows were sorted, whereas in other hosts, the occurrence of significant nestedness depended on the order of the matrix rows. The degree of phylogenetic and functional nestedness for matrices sorted by the sum of branch lengths was associated with a host species' morphoecological traits and the latitude of its geographic range. We conclude that consideration of nestedness based solely on species composition does not allow a comprehensive understanding of the patterns of parasite community structure. Nestedness should also be considered from phylogenetic and functional perspectives.

Host and parasite identity interact in scale-dependent fashion to determine parasite community structure

Understanding the ecological assembly of parasite communities is critical to characterise how changing host and environmental landscapes will alter infection dynamics and outcomes. However, studies frequently assume that (a) closely related parasite species or those with identical life-history strategies are functionally equivalent, and (b) the same factors will drive infection dynamics for a single parasite across multiple host species, oversimplifying community assembly patterns. Here, we challenge these two assumptions using a naturally occurring host-parasite system, with the mussel Anodonta anatina infected by the digenean trematode Echinoparyphium recurvatum, and the snail Viviparus viviparus infected by both E. recurvatum and Echinostoma sp. By analysing the impact of temporal parasite dispersal, host species and size, and the impact of coinfection (moving from broader environmental factors to within-host dynamics), we show that neither assumption holds true, but at different ecological scales. The assumption that closely related parasites can be functionally grouped is challenged when considering dispersal to the host (i.e. larger scales), while the assumption that the same factors will drive infection dynamics for a single parasite across multiple host species is challenged when considering within-host interspecific competition (i.e. smaller scales). Our results demonstrate that host identity, parasite identity and ecological scale require simultaneous consideration in studies of parasite community composition and transmission.

Compositional and phylogenetic nestedness of host assemblages exploited by generalist ectoparasites across their geographic ranges: drivers and associations with ectoparasite traits

We investigated compositional and phylogenetic nestedness in the host assemblages of 26 host-generalist fleas across regions within the Palearctic. We asked the following questions: (i) are host assemblages exploited by a flea species compositionally or phylogenetically nested (=C-nested and P-nested, respectively) across regions?; (ii) if yes, what are the processes that generate nestedness, and does phylogenetic nestedness follow the same processes as compositional nestedness?; and (iii) are the biological traits of a flea species associated with its host assemblages' degree of nestedness? Nestedness was calculated for matrices with rows ordered either by decreasing region area (=a-matrices) or increasing distance from the centre of a flea's geographic range (d-matrices). Significant C-nestedness was found in either a- (three fleas) or d-matrices (three fleas) or both (10 fleas). Significant P-nestedness was detected in either a- (three fleas) or d-matrices (four fleas) or both (two fleas). In some but not other species, P-nestedness followed C-nestedness. The probability of C-nestedness to be significant, as well as its degree for d-matrices, was associated with a flea's morphoecological traits, whereas this was not the case for either a-matrices or the P-nestedness for either type of ordered matrices. We conclude that compositional, but not phylogenetic, nestedness is (i) generated by similar mechanisms in many flea species and (ii) may be simultaneously driven by different mechanisms in the same flea. In contrast, mechanisms promoting phylogenetic nestedness differ between flea species and seem to act separately.

Tapping into natural history collections to assess latitudinal gradients of parasite diversity

Parasites are key components of the biosphere not only due to their huge diversity, but also because they exert important influences on ecological processes. Nevertheless, we lack an understanding of the biogeographical patterns of parasite diversity. Here, we tap into the potential of biodiversity collections for understanding parasite biogeography. We assess species richness of supracommunities of helminth parasites infecting mammal assemblages in the Nearctic, and describe its relation to latitude, climate, host diversity, and land area. We compiled data from parasitology collections and assessed parasite diversity in Nearctic ecoregions for the entire parasite supracommunity of mammals in each ecoregion, as well as separately from carnivores and rodents to explore the effect of host taxonomic resolution on observed patterns. For carnivores, we found evidence of a negative latitudinal gradient, while parasites of rodents displayed no clear pattern. We found that parasite diversity was positively correlated with mean annual temperature and negatively correlated with seasonal precipitation. Parasite richness shows a diversity peak at intermediate host richness values and in carnivores correlates with temperature and seasonal precipitation. Rodent parasite diversity did not correlate with explored factors. Other researchers are encouraged to use parasitology collections to continue exploring patterns of parasite biogeography and macroecology.

Disentangling the determinants of symbiotic species richness in native and invasive gammarids (Crustacea, Amphipoda) of the Baltic region

Dispersal of alien species is a global problem threatening native biodiversity. Co-introduction of non-native parasites and pathogens adds to the severity of this threat, but this indirect impact has received less attention. To shed light on the key factors determining the richness of microorganisms in native and invasive host species, we compared symbiotic (parasitic and epibiotic) communities of gammarids across different habitats and localities along the Baltic coast of Poland. Seven gammarid species, two native and five invasive, were sampled from 16 freshwater and brackish localities. Sixty symbiotic species of microorganisms of nine phyla were identified. This taxonomically diverse species assemblage of symbionts allowed us to assess the effect of host translocation and regional ecological determinants driving assembly richness in the gammarid hosts. Our results revealed that (i) the current assemblages of symbionts of gammarid hosts in the Baltic region are formed by native and co-introduced species; (ii) species richness of the symbiotic community was higher in the native Gammarus pulex than in the invasive hosts, probably reflecting a process of species loss by invasive gammarids in the new area and the distinct habitat conditions occupied by G. pulex and invasive hosts; (iii) both host species and locality were key drivers shaping assembly composition of symbionts, whereas habitat condition (freshwater versus brackish) was a stronger determinant of communities than geographic distance; (iv) the dispersion patterns of the individual species richness of symbiotic communities were best described by Poisson distributions; in the case of an invasive host, the dispersion of the rich species diversity may switch to a right-skewed negative binomial distribution, suggesting a host-mediated regulation process. We believe this is the first analysis of the symbiotic species richness in native and invasive gammarid hosts in European waters based on original field data and a broad range of taxonomic groups including Microsporidia, Choanozoa, Ciliophora, Apicomplexa, Platyhelminthes, Nematoda, Nematomorha, Acanthocephala and Rotifera, to document the patterns of species composition and distribution.

Within-host and external environments differentially shape β-diversity across parasite life stages

Uncovering drivers of community assembly is a key aspect of learning how biological communities function. Drivers of community similarity can be especially useful in this task as they affect assemblage-level changes that lead to differences in species diversity between habitats. Concepts of β-diversity originally developed for use in free-living communities have been widely applied to parasite communities to gain insight into how infection risk changes with local conditions by comparing parasite communities across abiotic and biotic gradients. Factors shaping β-diversity in communities of immature parasites, such as larvae, are largely unknown. This is a key knowledge gap as larvae are frequently the infective life-stage and understanding variation in these larval communities is thus key for disease prevention. Our goal was to uncover links between β-diversity of parasite communities at different life stages; therefore, we used gastrointestinal nematodes infecting African buffalo in Kruger National Park, South Africa, to investigate within-host and extra-host drivers of adult and larval parasite community similarity. We employed a cross-sectional approach using PERMANOVA that examined each worm community at a single time point to assess independent drivers of β-diversity in larvae and adults as well as a longitudinal approach with path analysis where adult and larval communities from the same host were compared to better link drivers of β-diversity between these two life stages. Using the cross-sectional approach, we generally found that intrinsic, within-host traits had significant effects on β-diversity of adult nematode communities, while extrinsic, extra-host variables had significant effects on β-diversity of larval nematode communities. However, the longitudinal approach provided evidence that intrinsic, within-host factors affected the larval community indirectly via the adult community. Our results provide key data for the comparison of community-level processes where adult and immature stages inhabit vastly different habitats (i.e. within-host vs. abiotic environment). In the context of parasitism, this helps elucidate host infection risk via larval stages and the drivers that shape persistence of adult parasite assemblages, both of which are useful for predicting and preventing infectious disease.

Factors at multiple scales drive parasite community structure

Understanding how ecological communities are assembled remains a key goal of ecosystem ecology. Because communities are hierarchical, factors acting at multiple scales can contribute to patterns of community structure. Parasites provide a natural system to explore this idea, as they exist as discrete communities within host individuals, which are themselves part of a community and metacommunity. We aimed to understand the relative contribution of multi-scale drivers in parasite community assembly and assess how patterns at one level may mask those occurring at another. Specifically, we wanted to disentangle patterns caused by passive sampling from those determined by ecological drivers, and how these vary with scale. We applied a Markov Random Fields model and assessed measures of β-diversity and nestedness for 420 replicate parasite infracommunities (parasite assemblages in host individuals) across two freshwater mussel host species, three sites and two time periods, comparing our results to simulations from four different ecologically relevant null models. We showed that β-diversity between sites (explaining 25% of variation in parasite distribution) and host species (41%) is greater than expected, and β-diversity between individual hosts is smaller than expected, even after accounting for parasite prevalence and characteristics of host individuals. Furthermore, parasite communities were significantly less nested than expected once parasite prevalence and host characteristics were both accounted for, but more nested than expected otherwise, suggesting a degree of modularity at the within-host level that is masked if underlying host and parasite characteristics are not taken into account. The Markov Random Fields model provided evidence for possible competitive within-host parasite interactions, providing a mechanism for the observed infracommunity modularity. An integrative approach that examines factors at multiple scales is necessary to understand the composition of ecological communities. Furthermore, patterns at one level can alter the interpretation of ecologically important drivers at another if variation at higher scales is not accounted for.

Spatial Factors Outperform Local Environmental and Geo-Climatic Variables in Structuring Multiple Facets of Stream Macroinvertebrates' β-Diversity

Simple Summary

One of the key targets of community ecology and biogeography concerns revealing the variability and underlying drivers of biodiversity. Most current studies understand biodiversity based on taxonomic information alone. Our study was based on macroinvertebrates from 179 stream sampling sites in the Hun-Tai River Basin in Northeastern China. The correlation of different facets of β-diversity was compared while revealing the relative contribution of multiple abiotic factors (i.e., local environmental, geo-climatic, and spatial factors) to shaping β-diversity based on taxonomic, functional, and phylogenetic information. The results showed that functional β-diversity provides important complementary information to taxonomic and phylogenetic β-diversity. Moreover, spatial factors outperform local environmental and geo-climatic variables in structuring multiple facets of stream macroinvertebrates’ β-diversity. Our study provides guidance for future conservation studies of watershed biodiversity, as well as implications for future studies of β-diversity.

Abstract

One of the key targets of community ecology and biogeography concerns revealing the variability and underlying drivers of biodiversity. Most current studies understand biodiversity based on taxonomic information alone, but few studies have shown the relative contributions of multiple abiotic factors in shaping biodiversity based on taxonomic, functional, and phylogenetic information. We collected 179 samples of macroinvertebrates in the Hun-Tai River Basin. We validated the complementarity between the three facets and components of β-diversity using the Mantel test. Distance-based redundancy analysis and variance partitioning were applied to explore the comparative importance of local environmental, geo-climatic, and spatial factors on each facet and component of β-diversity. Our study found that taxonomic and phylogenetic total β-diversity was mainly forced by turnover, while functional total β-diversity was largely contributed by nestedness. There is a strong correlation between taxonomic and phylogenetic β-diversity. However, the correlations of functional with both taxonomic and phylogenetic β-diversity were relatively weak. The findings of variation partitioning suggested that distinct facets and components of macroinvertebrates’ β-diversity were impacted by abiotic factors to varying degrees. The contribution of spatial factors was greater than that of the local environment and geo-climatic factors for taxonomic, functional, and phylogenetic β-diversity. Thus, studying different facets and components of β-diversity allows a clearer comprehension of the influence of abiotic factors on diversity patterns. Therefore, future research should investigate patterns and mechanisms of β-diversity from taxonomic, functional, and phylogenetic perspectives.

Component Endoparasite Communities Mirror Life-History Specialization in Syntopic Reed Frogs (Hyperolius spp.)

Most of our knowledge on the processes structuring parasite communities in amphibians originate from temperate-zone taxa, whereas Afrotropical communities have been neglected so far. We found evidence that ecological fitting of the hosts and, probably, differential immune response may influence the variation in parasite species richness, prevalence, and infestation intensity of East African frogs Hyperolius kivuensis and H. viridiflavus. The most closely related host species share the same macrohabitat (that implies the same pool of potential parasites), but differ in microhabitat preference, so that a comparative analyses of syntopic and allopatric populations is expedient to reveal ecological fitting. We detected 11 parasite species (one annelid, four nematodes, five trematodes, one cestode) and two endocommensal species (protozoans). The component parasite communities included 4–5 helminth species in H. kivuensis and 6–8 in the more aquatic H. viridiflavus, supporting the hypothesis that trematode diversity increases with the amount of time spent in water. Five parasite species (Orneoascaris chrysanthemoides, Clinostomum chabaudi, an undetermined echinostomatid) and two protozoans (Nyctotheroides sp., and Protoopalina sp.) are shared among the syntopic amphibian populations. This finding indicates a similar susceptibility of these amphibians to infestation from the local parasite pool. Yet, the low prevalence of single- and multi-species infestations in H. kivuensis indicates that parasite clearing by its immune response is probably more effective and prominent than in H. viridiflavus. Therefore, H. viridiflavus suffered from significantly reduced short-term survival due to the infection. Thus, we conclude that the processes structuring component parasite communities in amphibians do not differ generally between temperate-zone and Afrotropical host species, but they do in the magnitude of ecological fitting.

Abundance data applied to a novel model invertebrate host shed new light on parasite community assembly in nature

Understanding how environmental drivers influence the assembly of parasite communities, in addition to how parasites may interact at an infracommunity level, are fundamental requirements for the study of parasite ecology. Knowledge of how parasite communities are assembled will help to predict the risk of parasitism for hosts, and model how parasite communities may change under variable conditions. However, studies frequently rely on presence-absence data and examine multiple host species or sites, metrics which may be too coarse to characterise nuanced within-host patterns. We utilised a novel host system, the freshwater mussel Anodonta anatina, to investigate the drivers of community structure and explore parasite interactions. In addition, we aimed to highlight consistencies and inconsistencies between PA and abundance data. Our analysis incorporated 14 parasite taxa and 720 replicate infracommunities. Using Redundancy Analysis, a joint species distribution model and a Markov random field approach, we modelled the impact of both host-level and environment-level characteristics on parasite structure, as well as parasite-parasite correlations after accounting for all other factors. This approach was repeated for both the presence and abundance of all parasites. We demonstrated that the regional species pool, individual host characteristics (mussel length and gravidity) and predicted parasite-parasite interactions are all important but to varying degrees across parasite species, suggesting that applying generalities to parasite community construction is too simplistic. Furthermore, we showed that PA data fail to capture important density-dependent effects of parasite load for parasites with high abundance, and in general performs poorly for high-intensity parasites. Host and parasite traits, as well as broader environmental factors, all contribute to parasite community structure, emphasising that an integrated approach is required to study community assembly. However, care must be taken with the data used to infer patterns, as presence-absence data may lead to incorrect ecological inference.

Scale-dependent effects of host patch traits on species composition in a stickleback parasite metacommunity

A core goal of ecology is to understand the abiotic and biotic variables that regulate species distributions and community composition. A major obstacle is that the rules governing species distributions can change with spatial scale. Here, we illustrate this point using data from a spatially nested metacommunity of parasites infecting a metapopulation of threespine stickleback fish from 34 lakes on Vancouver Island, British Columbia. Like most parasite metacommunities, the composition of stickleback parasites differs among host individuals within each host population, and differs between host populations. The distribution of each parasite taxon depends, to varying degrees, on individual host traits (e.g., mass, diet) and on host-population characteristics (e.g., lake size, mean host mass, mean diet). However, in most cases in this data set, a given parasite was regulated by different factors at the host-individual and host-population scales, leading to scale-dependent patterns of parasite-species co-occurrence.