Search

IntroductionInsects share intimate relationships with microbes that play important roles in their biology. Yet our understanding of how host-bound microbial communities assemble and perpetuate over evolutionary time is limited. Ants host a wide range of microbes with diverse functions and are an emerging model for studying the evolution of insect microbiomes. Here, we ask whether phylogenetically related ant species have formed distinct and stable microbiomes.MethodsTo answer this question, we investigated the microbial communities associated with queens of 14 Formica species from five clades, using deep coverage 16S rRNA amplicon sequencing.ResultsWe reveal that Formica species and clades harbor highly defined microbial communities that are dominated by four bacteria genera: Wolbachia, Lactobacillus, Liliensternia, and Spiroplasma. Our analysis reveals that the composition of Formica microbiomes mirrors the phylogeny of the host, i.e., phylosymbiosis, in that related hosts harbor more similar microbial communities. In addition, we find there are significant correlations between microbe co-occurrences.DiscussionOur results demonstrate Formica ants carry microbial communities that recapitulate the phylogeny of their hosts. Our data suggests that the co-occurrence of different bacteria genera may at least in part be due to synergistic and antagonistic interactions between microbes. Additional factors potentially contributing to the phylosymbiotic signal are discussed, including host phylogenetic relatedness, host-microbe genetic compatibility, modes of transmission, and similarities in host ecologies (e.g., diets). Overall, our results support the growing body of evidence that microbial community composition closely depends on the phylogeny of their hosts, despite bacteria having diverse modes of transmission and localization within the host.

In the present study, an effort has been made to elucidate the gut bacterial diversity of twelve species of the family Araneidae under three subfamilies collected from 5 states of India along with their predicted metabolic role in functional metabolism. Further, we also compared the host species phylogeny based on partial cytochrome c oxidase subunit I (COI) sequences with the gut bacteria composition dendrogram to decipher the phylosymbiotic relationships. Analysis revealed the presence of 22 bacterial phyla, 145 families, and 364 genera in the gut, with Proteobacteria, Firmicutes, Actinobacteria, and Deinococcus-Thermus as the highest abundant phyla. Moreover, phylum Bacteriodetes was dominated only in Cyclosa mulmeinensis and Chlamydiae in Neoscona bengalensis. At the genus level, Bacillus, Acinetobacter, Cutibacterium, Pseudomonas, and Staphylococcus were the most dominant genera. Furthermore, the genus Prevotella was observed only in Cyclosa mulmeinensis, and endosymbiont Wolbachia only in Eriovixia laglaizei. The differential abundance analysis (DeSeq2) revealed the 19 significant ASVs represented by the genera like Acinetobacter, Vagoccoccus, Prevotella, Staphylococcus, Curvibacter, Corynebacterium, Paracoccus, Streptococcus, Microbacterium, and Pseudocitrobacter. The inter- and intra-subfamilies comparison based on diversity indices (alpha and beta diversity) revealed that the subfamily Araneinae have high richness and diversity than Argiopinae and Gasteracanthinae. The phylosymbiotic analysis revealed that there is no congruence between the gut bacteria composition dendrogram with their host phylogeny.

International audience; The degree of similarity between the microbiotas of host species often mirrors the phylogenetic proximity of the hosts. This pattern, referred to as phylosymbiosis, is widespread in animals and plants. While phylosymbiosis was initially interpreted as the signal of symbiotic transmission and coevolution between microbes and their hosts, it is now recognized that similar patterns can emerge even if the microbes are environmentally acquired. Distinguishing between these two scenarios, however, remains challenging. We recently developed HOME (host-microbiota evolution), a cophylogenetic model designed to detect vertically transmitted microbes and host switches from amplicon sequencing data. Here, we applied HOME to the microbiotas of Hawaiian spiders of the genus Ariamnes, which experienced a recent radiation on the archipelago. We demonstrate that although Hawaiian Ariamnes spiders display a significant phylosymbiosis, there is little evidence of microbial vertical transmission. Next, we performed simulations to validate the absence of transmitted microbes in Ariamnes spiders. We show that this is not due to a lack of detection power because of the low number of segregating sites or an effect of phylogenetically driven or geographically driven host switches. Ariamnes spiders and their associated microbes therefore provide an example of a pattern of phylosymbiosis likely emerging from processes other than vertical transmission. IMPORTANCE How host-associated microbiotas assemble and evolve is one of the outstanding questions of microbial ecology. Studies aiming at answering this question have repeatedly found a pattern of "phylosymbiosis," that is, a phylogenetic signal in the composition of host-associated microbiotas. While phylosymbiosis was often interpreted as evidence for vertical transmission and host-microbiota coevolution, simulations have now shown that it can emerge from other processes, including host filtering of environmentally acquired microbes. However, distinguishing the processes driving phylosymbiosis in nature remains challenging. We recently developed a cophylogenetic method that can detect vertical transmission. Here, we applied this method to the microbiotas of recently diverged spiders from the Hawaiian archipelago, which display a clear phylosymbiosis pattern. We found that none of the bacterial operational taxonomic units is vertically transmitted. We show with simulations that this result is not due to methodological artifacts. Thus, we provide a striking empirical example of phylosymbiosis emerging from processes other than vertical transmission. KEYWORDS microbiota, phylosymbiosis, vertical transmission, host filtering, Hawaiian arthropods, endosymbionts M ost multicellular organisms host complex microbial communities, referred to as microbiotas, which provide important functions to their hosts (1, 2). Although these microbial communities can fluctuate over short timescales and according to external variables such as animal diet or soil composition (3-5), microbiotas of host individuals from the

Most plants rely on specialized root-associated microbes to obtain essential nitrogen (N), yet not much is known about the evolutionary history of the rhizosphere–plant interaction. We conducted a common garden experiment to investigate the plant root–rhizosphere microbiome association using chloridoid grasses sampled from around the world and grown from seed in a greenhouse. We sought to test whether plants that are more closely related phylogenetically have more similar root bacterial microbiomes than plants that are more distantly related. Using metagenome sequencing, we found that there is a conserved core and a variable rhizosphere bacterial microbiome across the chloridoid grasses. Additionally, phylogenetic distance among the host plant species was correlated with bacterial community composition, suggesting the plant hosts prefer specific bacterial lineages. The functional potential for N utilization across microbiomes fluctuated extensively and mirrored variation in the microbial community composition across host plants. Variation in the bacterial potential for N fixation was strongly affected by the host plants’ phylogeny, whereas variation in N recycling, nitrification, and denitrification was unaffected. This study highlights the evolutionary linkage between the N fixation traits of the microbial community and the plant host and suggests that not all functional traits are equally important for plant–microbe associations.

Background Cyanobacteria from the genus Microcystis can form large mucilaginous colonies with attached heterotrophic bacteria—their microbiome. However, the nature of the relationship between Microcystis and its microbiome remains unclear. Is it a long-term, evolutionarily stable association? Which partners benefit? Here we report the genomic diversity of 109 individual Microcystis colonies—including cyanobacteria and associated bacterial genomes—isolated in situ and without culture from Lake Champlain, Canada and Pampulha Reservoir, Brazil. Results We identified 14 distinct Microcystis genotypes from Canada, of which only two have been previously reported, and four genotypes specific to Brazil. Microcystis genetic diversity was much greater between than within colonies, consistent with colony growth by clonal expansion rather than aggregation of Microcystis cells. We also identified 72 bacterial species in the microbiome. Each Microcystis genotype had a distinct microbiome composition, and more closely related genotypes had more similar microbiomes. This pattern of phylosymbiosis could be explained by co-phylogeny in only two out of the nine most prevalent associated bacterial genera, Roseomonas and Rhodobacter. These phylogenetically associated genera could enrich the metabolic repertoire of Microcystis, for example by encoding the biosynthesis of complementary carotenoid molecules. In contrast, other colony-associated bacteria showed weaker signals of co-phylogeny, but stronger evidence of horizontal gene transfer with Microcystis. These observations suggest that acquired genes are more likely to be retained in both partners (Microcystis and members of its microbiome) when they are loosely associated, whereas one gene copy is sufficient when the association is physically tight and evolutionarily long-lasting. Conclusions We have introduced a method for culture-free isolation of single colonies from nature followed by metagenomic sequencing, which could be applied to other types of microbes. Together, our results expand the known genetic diversity of both Microcystis and its microbiome in natural settings, and support their long-term, specific, and potentially beneficial associations.

Ostreobium, the major algal symbiont of the coral skeleton, remains understudied despite extensive research on the coral holobiont. The enclosed nature of the coral skeleton might reduce the dispersal and exposure of residing bacteria to the outside environment, allowing stronger associations with the algae. Here, we describe the bacterial communities associated with cultured strains of 5 Ostreobium clades using 16S rRNA sequencing. We shed light on their likely physical associations by comparative analysis of three datasets generated to capture (1) all algae associated bacteria, (2) enriched tightly attached and potential intracellular bacteria, and (3) bacteria in spent media. Our data showed that while some bacteria may be loosely attached, some tend to be tightly attached or potentially intracellular. Although colonised with diverse bacteria, Ostreobium preferentially associated with 34 bacterial taxa revealing a core microbiome. These bacteria include known nitrogen cyclers, polysaccharide degraders, sulphate reducers, antimicrobial compound producers, methylotrophs, and vitamin B12 producers. By analysing co-occurrence networks of 16S rRNA datasets from Porites lutea and Paragoniastrea australensis skeleton samples, we show that the Ostreobium-bacterial associations present in the cultures are likely to also occur in their natural environment. Finally, our data show significant congruence between the Ostreobium phylogeny and the community composition of its tightly associated microbiome, largely due to the phylosymbiotic signal originating from the core bacterial taxa. This study offers insight into the Ostreobium microbiome and reveals preferential associations that warrant further testing from functional and evolutionary perspectives.

Insects share intimate relationships with microbes that play important roles in their biology. Yet our understanding of how host-bound microbial communities assemble and perpetuate over evolutionary time is limited. Ants host a wide range of microbes with diverse functions and are an emerging model for studying the evolution of insect microbiomes. Here, we ask whether phylogenetically related ant species have formed distinct and stable microbiomes. To answer this question, we investigated the microbial communities associated with queens of 14Formicaspecies from five clades, using deep coverage 16S rRNA amplicon sequencing. We reveal thatFormicaspecies and clades harbour highly defined microbial communities that are dominated by four bacteria genera:Wolbachia, Lactobacillus,Liliensternia,andSpiroplasma. Our analysis reveals that the composition ofFormicamicrobiomes mirrors the phylogeny of the host, i.e. phylosymbiosis, in that related hosts harbour more similar microbial communities. Our analysis also revealed significant correlations between microbe co-occurrences, which suggests that synergistic and antagonistic interactions may contribute to the phylosymbiotic signal. Additional factors potentially contributing to the phylosymbiotic signal are discussed, including host phylogenetic relatedness, host-microbe genetic compatibility, modes of transmission, and similarities in host ecologies (e.g., diets). Overall, our results support the growing body of evidence that microbial community composition closely depends on the phylogeny of their hosts, despite bacteria having diverse modes of transmission and localisation within the host.

Background Microbiomes can have profound impacts on host biology and evolution, but to date, remain vastly understudied in spiders despite their unique and diverse predatory adaptations. This study evaluates closely related species of spiders and their host-microbe relationships in the context of phylosymbiosis, an eco-evolutionary pattern where the microbial community profile parallels the phylogeny of closely related host species. Using 16S rRNA gene amplicon sequencing, we characterized the microbiomes of five species with known phylogenetic relationships from the family Theridiidae, including multiple closely related widow spiders (L. hesperus, L. mactans, L. geometricus, S. grossa, and P. tepidariorum). Results We compared whole animal and tissue-specific microbiomes (cephalothorax, fat bodies, venom glands, silk glands, and ovary) in the five species to better understand the relationship between spiders and their microbial symbionts. This showed a strong congruence of the microbiome beta-diversity of the whole spiders, cephalothorax, venom glands, and silk glands when compared to their host phylogeny. Our results support phylosymbiosis in these species and across their specialized tissues. The ovary tissue microbial dendrograms also parallel the widow phylogeny, suggesting vertical transfer of species-specific bacterial symbionts. By cross-validating with RNA sequencing data obtained from the venom glands, silk glands and ovaries of L. hesperus, L. geometricus, S. grossa, and P. tepidariorum we confirmed that several microbial symbionts of interest are viably active in the host. Conclusion Together these results provide evidence that supports the importance of host-microbe interactions and the significant role microbial communities may play in the evolution and adaptation of their hosts.

Background The vertebrate clade diverged into Chondrichthyes (sharks, rays, and chimeras) and Osteichthyes fishes (bony fishes) approximately 420 mya, with each group accumulating vast anatomical and physiological differences, including skin properties. The skin of Chondrichthyes fishes is covered in dermal denticles, whereas Osteichthyes fishes are covered in scales and are mucous rich. The divergence time among these two fish groups is hypothesized to result in predictable variation among symbionts. Here, using shotgun metagenomics, we test if patterns of diversity in the skin surface microbiome across the two fish clades match predictions made by phylosymbiosis theory. We hypothesize (1) the skin microbiome will be host and clade-specific, (2) evolutionary difference in elasmobranch and teleost will correspond with a concomitant increase in host-microbiome dissimilarity, and (3) the skin structure of the two groups will affect the taxonomic and functional composition of the microbiomes. Results We show that the taxonomic and functional composition of the microbiomes is host-specific. Teleost fish had lower average microbiome within clade similarity compared to among clade comparison, but their composition is not different among clade in a null based model. Elasmobranch’s average similarity within clade was not different than across clade and not different in a null based model of comparison. In the comparison of host distance with microbiome distance, we found that the taxonomic composition of the microbiome was related to host distance for the elasmobranchs, but not the teleost fishes. In comparison, the gene function composition was not related to the host-organism distance for elasmobranchs but was negatively correlated with host distance for teleost fishes. Conclusion Our results show the patterns of phylosymbiosis are not consistent across both fish clades, with the elasmobranchs showing phylosymbiosis, while the teleost fish are not. The discrepancy may be linked to alternative processes underpinning microbiome assemblage, including possible historical host-microbiome evolution of the elasmobranchs and convergent evolution in the teleost which filter specific microbial groups. Our comparison of the microbiomes among fishes represents an investigation into the microbial relationships of the oldest divergence of extant vertebrate hosts and reveals that microbial relationships are not consistent across evolutionary timescales.

International audience; Planorbidae snails are the intermediate host for the trematode parasite of the Schistosoma genus which is responsible for schistosomiasis, a disease that affects both humans and cattle. The microbiota for Schistosoma has already been described as having an affect on host/parasite interactions, specifically through immunological interactions. Here we sought to characterize the microbiota composition of seven Planorbidae species and strains. Individual snail microbiota was determined using 16S rDNA amplicon sequencing. The bacterial composition was highly specific to host strain with limited inter-individual variation. In addition, it displayed complete congruence with host phylogeny, revealing a phylosymbiosis pattern. These results were confirmed in common garden, suggesting that the host highly constrains microbial composition. This study presents the first comparison of bacterial communities between several intermediate snail hosts of Schistosoma parasites, paving the way for further studies on the understanding of this tripartite interaction.

Ostreobium, the major algal symbiont of the coral skeleton, remains understudied despite extensive research on the coral holobiont. The enclosed nature of the coral skeleton might reduce the dispersal and exposure of residing bacteria to the outside environment, allowing stronger associations with the algae. Here, we describe the bacterial communities associated with cultured strains of 5Ostreobiumclades using 16S rRNA sequencing. We shed light on their likely physical associations by comparative analysis of three datasets generated to capture (1) all algae associated bacteria (2) enriched tightly attached and potential intracellular bacteria and (3) bacteria in spent media. Our data showed that while some bacteria may be loosely attached, some tend to be tightly attached or potentially intracellular. Although colonised with diverse bacteria,Ostreobiumpreferentially associated with 34 bacterial taxa revealing a core microbiome. These bacteria include taxa known as nitrogen cyclers, polysaccharide degraders, sulphate reducers, antimicrobial compound producers, methylotrophs and vitamin B12 producers. By analysing co-occurrence networks of 16S rRNA datasets fromPorites luteaandParagoniastrea australensisskeleton samples, we show that theOstreobium-bacterial associations present in the cultures are likely to also occur in their natural environment. Finally, our data show significant congruence between theOstreobiumphylogeny and the community composition of its tightly associated microbiome, largely due to the phylosymbiotic signal originating from the core bacterial taxa. This study offers insight into theOstreobiummicrobiome and reveals preferential associations that warrant further testing from functional and evolutionary perspectives.

Significance Skin forms a critical protective barrier between a mammal and its external environment. Baseline data on the mammalian skin microbiome elucidates which microorganisms are found on healthy skin and provides insight into mammalian evolutionary history. To our knowledge, this study represents the largest existing mammalian skin microbiome survey. Our findings demonstrate that human skin is distinct, not only from other Primates, but from all 10 mammalian orders sampled. Identifying significant similarities between branching of mammalian phylogenetic trees and relatedness trees for their corresponding microbial communities raises the possibility that mammals have experienced coevolution between skin microbiota and their corresponding host species., Skin is the largest organ of the body and represents the primary physical barrier between mammals and their external environment, yet the factors that govern skin microbial community composition among mammals are poorly understood. The objective of this research was to generate a skin microbiota baseline for members of the class Mammalia, testing the effects of host species, geographic location, body region, and biological sex. Skin from the back, torso, and inner thighs of 177 nonhuman mammals was sampled, representing individuals from 38 species and 10 mammalian orders. Animals were sampled from farms, zoos, households, and the wild. The DNA extracts from all skin swabs were amplified by PCR and sequenced, targeting the V3-V4 regions of bacterial and archaeal 16S rRNA genes. Previously published skin microbiome data from 20 human participants, sampled and sequenced using an identical protocol to the nonhuman mammals, were included to make this a comprehensive analysis. Human skin microbial communities were distinct and significantly less diverse than all other sampled mammalian orders. The factor most strongly associated with microbial community data for all samples was whether the host was a human. Within nonhuman samples, host taxonomic order was the most significant factor influencing skin microbiota, followed by the geographic location of the habitat. By comparing the congruence between host phylogeny and microbial community dendrograms, we observed that Artiodactyla (even-toed ungulates) and Perissodactyla (odd-toed ungulates) had significant congruence, providing evidence of phylosymbiosis between skin microbial communities and their hosts.

BackgroundWhile many studies have reported that the structure of the gut and skin microbiome is driven by both species-specific and habitat-specific factors, the relative importance of host-specific versus environmental factors in wild vertebrates remains poorly understood. The aim of this study was to determine the diversity and composition of fish skin, gut, and surrounding water bacterial communities (BCs) and assess the extent to which host habitat and phylogeny predict BC similarity. Skin swabs and gut samples from 334 fish belonging to 17 species were sampled in three Laurentian Great Lakes (LGLs) habitats (Detroit River, Lake Erie, Lake Ontario). We also collected and filtered water samples at the time of fish collection. We analyzed BC composition using 16S metabarcoding and tested for community variation. ResultsWe found that the water BC was distinct from the fish microbiome BCs, although the skin BC more closely resembled the water BC. We also found that environmental (sample location), habitat, fish diet and host species factors shape and promote divergence or convergence of the fish BCs. Since host species significantly affected both gut and skin BC (separately from host species effects), we tested for phylosymbiosis using pairwise host species phylogenetic distance versus BC dissimilarity. We found significant phylogenetic effects on BC dissimilarity, consistent with phylosymbiosis for both the fish skin and gut BCs, perhaps reflecting the longstanding co-evolutionary relationship between the host species and their microbiomes. ConclusionsAnalyzing the gut content and skin mucus BCs across diverse fish species in complex natural ecosystems such as the LGLs provides insights into the potential for habitat and species-specific effects on the microbiome, and ultimately the health, of the host.

Gut microbiota have long attracted the interest of scientists due to their profound impact on the well-being of animals. A non-random pattern of microbial assembly that results in a parallelism between host phylogeny and microbial similarity is described as phylosymbiosis. Phylosymbiosis has been consistently observed in different clades of animal hosts, but there have been no studies on crustaceans. In this study, we investigated whether host phylogeny has an impact on the gut microbiota assemblages in decapod shrimps. We examined the gut microbial communities in 20 shrimp species from three families inhabiting distinct environments, using metabarcoding analyses of the V1-V3 hypervariable region of the 16S rRNA gene. Gut microbial communities varied within each shrimp group but were generally dominated by Proteobacteria. A prevalent phylosymbiotic pattern in shrimps was evidenced for the first time by the observations of (1) the distinguishability of microbial communities among species within each group, (2) a significantly lower intraspecific than interspecific gut microbial beta diversity across shrimp groups, (3) topological congruence between host phylogenetic trees and gut microbiota dendrograms, and (4) a correlation between host genetic distances and microbial dissimilarities. Consistent signals of phylosymbiosis were observed across all groups in dendrograms based on the unweighted UniFrac distances at 99% operational taxonomic units (OTUs) level and in Mantel tests based on the weighted UniFrac distances based on 97% OTUs and amplicon sequence variants. Penaeids exhibited phylosymbiosis in most tests, while phylosymbiotic signals in atyids and pandalids were only detected in fewer than half of the tests. A weak phylogenetic signal was detected in the predicted functions of the penaeid gut microbiota. However, the functional diversities of the two caridean groups were not significantly related to host phylogeny. Our observations of a parallelism in the taxonomy of the gut microbiota with host phylogeny for all shrimp groups examined and in the predicted functions for the penaeid shrimps indicate a tight host-microbial relationship during evolution.

Microbial communities can be critical for many metazoans, which can lead to the observation of phylosymbiosis with phylogenetically related species sharing similar microbial communities. Most of the previous studies on phylosymbiosis were conducted across the host families or genera. However, it is unclear whether the phylosymbiosis signal is still prevalent at lower taxonomic levels. In this study, 54 individuals from six species of the fig wasp genus Ceratosolen (Hymenoptera: Agaonidae) collected from nine natural populations and their associated microbiota were investigated. The fig wasp species were morphologically identified and further determined by mitochondrial CO1 gene fragments and nuclear ITS2 sequences, and the V4 region of 16S rRNA gene was sequenced to analyze the bacterial communities. The results suggest a significant positive correlation between host genetic characteristics and microbial diversity characteristics, indicating the phylosymbiosis signal between the phylogeny of insect hosts and the associated microbiota in the lower classification level within a genus. Moreover, we found that the endosymbiotic Wolbachia carried by fig wasps led to a decrease in bacterial diversity of host-associated microbial communities. This study contributes to our understanding of the role of host phylogeny, as well as the role of endosymbionts in shaping the host-associated microbial community.

Phylosymbiosis is an ecoevolutionary hypothesis and emerging pattern in animal-microbiota studies whereby the host phylogenetic relationships parallel the community relationships of the host-associated microbiota. A central prediction of phylosymbiosis is that closely related hosts exhibit a lower microbiota beta diversity than distantly related hosts. While phylosymbiosis has emerged as a widespread trend in a field often challenged to find trends across systems, two critical and understudied questions are whether or not phylosymbiosis is consequential to host biology and if adaptive evolutionary forces underpin the pattern. Here, using germfree rearing in the phylosymbiosis model Nasonia, we demonstrate that early life exposure to heat-inactivated microbiota from more distantly related species poses more severe developmental and survival costs than microbiota from closely related or the same species. This study advances a functional understanding of the consequences and potential selective pressures underpinning phylosymbiosis., Phylosymbiosis is defined as microbial community relationships that recapitulate the phylogeny of hosts. As evidence for phylosymbiosis rapidly accumulates in different vertebrate and invertebrate holobionts, a central question is what evolutionary forces cause this pattern. We use intra- and interspecific gut microbiota transplants to test for evidence of selective pressures that contribute to phylosymbiosis. We leverage three closely related species of the parasitoid wasp model Nasonia that recently diverged between 0.4 and 1 million years ago: N. vitripennis, N. giraulti, and N. longicornis. Upon exposure of germfree larvae to heat-inactivated microbiota from intra- or interspecific larvae, we measure larval growth, pupation rate, and adult reproductive capacity. We report three key findings: (i) larval growth significantly slows when hosts receive an interspecific versus intraspecific gut microbiota, (ii) marked decreases in pupation and resulting adult survival occur from interspecific gut microbiota exposure, and (iii) adult reproductive capacities including male fertility and longevity are unaffected by early life exposure to an interspecific microbiota. Overall, these findings reveal developmental and survival costs to Nasonia upon larval exposures to interspecific microbiota and provide evidence that selective pressures on phenotypes produced by host-microbiota interactions may underpin phylosymbiosis.

Summary Symbiotic association is universal in nature, and an array of symbionts play a crucial part in host life history. Aphids and their diverse symbionts have become a good model system to study insect‐symbiont interactions. Previous symbiotic diversity surveys have mainly focused on a few aphid clades, and the relative importance of different factors regulating microbial community structure is not well understood. In this study, we collected 65 colonies representing eight species of the aphid genus Mollitrichosiphum from different regions and plants in southern China and Nepal and characterized their microbial compositions using Illumina sequencing of the V3 − V4 hypervariable region of the 16S rRNA gene. We evaluated how microbiota varied across aphid species, geography and host plants and the correlation between microbial community structure and host aphid phylogeny. Heritable symbionts dominated the microbiota associated with Mollitrichosiphum, and multiple infections of secondary symbionts were prevalent. Ordination analyses and statistical tests highlighted the contribution of aphid species in shaping the structures of bacterial, symbiont and secondary symbiont communities. Moreover, we observed a significant correlation between Mollitrichosiphum aphid phylogeny and microbial community composition, providing evidence for a pattern of phylosymbiosis between natural aphid populations and their microbial associates.

Microbiome assemblages of plants and animals often show a degree of correlation with host phylogeny; an eco-evolutionary pattern known as phylosymbiosis. Using 16S rRNA gene sequencing to profile the microbiome, paired with COI, 18S rRNA and ITS1 host phylogenies, phylosymbiosis was investigated in four groups of coral reef invertebrates (scleractinian corals, octocorals, sponges and ascidians). We tested three commonly used metrics to evaluate the extent of phylosymbiosis: (a) intraspecific versus interspecific microbiome variation, (b) topological comparisons between host phylogeny and hierarchical clustering (dendrogram) of host-associated microbial communities, and (c) correlation of host phylogenetic distance with microbial community dissimilarity. In all instances, intraspecific variation in microbiome composition was significantly lower than interspecific variation. Similarly, topological congruency between host phylogeny and the associated microbial dendrogram was more significant than would be expected by chance across all groups, except when using unweighted UniFrac distance (compared with weighted UniFrac and Bray-Curtis dissimilarity). Interestingly, all but the ascidians showed a significant positive correlation between host phylogenetic distance and associated microbial dissimilarity. Our findings provide new perspectives on the diverse nature of marine phylosymbioses and the complex roles of the microbiome in the evolution of marine invertebrates.

Aphids and their symbionts represent an outstanding model for studies of insect–symbiont interactions. The aphid microbiota can be shaped by aphid species, geography and host plants. However, the relative importance of phylogenetic and ecological factors in shaping microbial community structures is not well understood. Using Illumina sequencing of the V3–V4 hypervariable region of the 16S rRNA gene, we characterized the microbial compositions of 215 aphid colonies representing 53 species of the aphid subfamily Greenideinae from different regions and plants in China, Nepal, and Vietnam. The primary endosymbiont Buchnera aphidicola and secondary symbiont Serratia symbiotica dominated the microbiota of Greenideinae. We simultaneously explored the relative contribution of host identity (i.e., aphid genus and aphid species), geography and host plant to the structures of bacterial, symbiont and secondary symbiont communities. Ordination analyses and statistical tests highlighted the strongest impact of aphid species on the microbial flora in Greenideinae. Furthermore, we found a phylosymbiosis pattern in natural Greenideinae populations, in which the aphid phylogeny was positively correlated with microbial community dissimilarities. These findings will advance our knowledge of host-associated microbiota assembly across both host phylogenetic and ecological contexts.

Viruses are the most abundant biological entity on the planet and interact with microbial communities with which they associate. The virome of animals is often dominated by bacterial viruses, known as bacteriophages or phages, which can (re)structure bacterial communities potentially vital to the animal host. Beta diversity relationships of animal-associated bacterial communities in laboratory and wild populations frequently parallel animal phylogenetic relationships, a pattern termed phylosymbiosis. However, little is known about whether viral communities also exhibit this eco-evolutionary pattern. Metagenomics of purified viruses from recently diverged species of Nasonia parasitoid wasps reared in the lab indicates for the first time that the community relationships of the virome can also exhibit complete phylosymbiosis. Therefore, viruses, particularly bacteriophages here, may also be influenced by animal evolutionary changes either directly or indirectly through the tripartite interactions among hosts, bacteria, and phage communities. Moreover, we report several new bacteriophage genomes from the common gut bacteria in Nasonia., Phylosymbiosis was recently proposed to describe the eco-evolutionary pattern whereby the ecological relatedness (e.g., beta diversity relationships) of host-associated microbial communities parallels the phylogeny of the host species. Representing the most abundant biological entities on the planet and common members of the animal-associated microbiome, viruses can be influential members of host-associated microbial communities that may recapitulate, reinforce, or ablate phylosymbiosis. Here we sequence the metagenomes of purified viral communities from three different parasitic wasp Nasonia species, one cytonuclear introgression line of Nasonia, and the flour moth outgroup Ephestia kuehniella. Results demonstrate complete phylosymbiosis between the viral metagenome and insect phylogeny. Across all Nasonia contigs, 69% of the genes in the viral metagenomes are either new to the databases or uncharacterized, yet over 99% of the contigs have at least one gene with similarity to a known sequence. The core Nasonia virome spans 21% of the total contigs, and the majority of that core is likely derived from induced prophages residing in the genomes of common Nasonia-associated bacterial genera: Proteus, Providencia, and Morganella. We also assemble the first complete viral particle genomes from Nasonia-associated gut bacteria. Taken together, results reveal the first complete evidence for phylosymbiosis in viral metagenomes, new genome sequences of viral particles from Nasonia-associated gut bacteria, and a large set of novel or uncharacterized genes in the Nasonia virome. This work suggests that phylosymbiosis at the host-microbiome level will likely extend to the host-virome level in other systems as well. IMPORTANCE Viruses are the most abundant biological entity on the planet and interact with microbial communities with which they associate. The virome of animals is often dominated by bacterial viruses, known as bacteriophages or phages, which can (re)structure bacterial communities potentially vital to the animal host. Beta diversity relationships of animal-associated bacterial communities in laboratory and wild populations frequently parallel animal phylogenetic relationships, a pattern termed phylosymbiosis. However, little is known about whether viral communities also exhibit this eco-evolutionary pattern. Metagenomics of purified viruses from recently diverged species of Nasonia parasitoid wasps reared in the lab indicates for the first time that the community relationships of the virome can also exhibit complete phylosymbiosis. Therefore, viruses, particularly bacteriophages here, may also be influenced by animal evolutionary changes either directly or indirectly through the tripartite interactions among hosts, bacteria, and phage communities. Moreover, we report several new bacteriophage genomes from the common gut bacteria in Nasonia.

Animals and microorganisms often establish close ecological relationships. However, much of our knowledge about animal microbiomes comes from two deeply studied groups: vertebrates and arthropods. To understand interactions on a broader scale of diversity, we characterized the bacterial microbiomes of close to 1,000 microscopic marine invertebrates from 21 phyla, spanning most of the remaining tree of metazoans. Samples were collected from five temperate and tropical locations covering three marine habitats (sediment, water column and intertidal macroalgae) and bacterial microbiomes were characterized using 16S ribosomal RNA gene sequencing. Our data show that, despite their size, these animals harbour bacterial communities that differ from those in the surrounding environment. Distantly related but coexisting invertebrates tend to share many of the same bacteria, suggesting that guilds of microorganisms preferentially associated with animals, but not tied to any specific host lineage, are the main drivers of the ecological relationship. Host identity is a minor factor shaping these microbiomes, which do not show the same correlation with host phylogeny, or 'phylosymbiosis', observed in many large animals. Hence, the current debate on the varying strength of phylosymbiosis within selected lineages should be reframed to account for the possibility that such a pattern might be the exception rather than the rule.

Scleractinian corals’ microbial symbionts influence host health, yet how coral microbiomes assembled over evolution is not well understood. We survey bacterial and archaeal communities in phylogenetically diverse Australian corals representing more than 425 million years of diversification. We show that coral microbiomes are anatomically compartmentalized in both modern microbial ecology and evolutionary assembly. Coral mucus, tissue, and skeleton microbiomes differ in microbial community composition, richness, and response to host vs. environmental drivers. We also find evidence of coral-microbe phylosymbiosis, in which coral microbiome composition and richness reflect coral phylogeny. Surprisingly, the coral skeleton represents the most biodiverse coral microbiome, and also shows the strongest evidence of phylosymbiosis. Interactions between bacterial and coral phylogeny significantly influence the abundance of four groups of bacteria–including Endozoicomonas-like bacteria, which divide into host-generalist and host-specific subclades. Together these results trace microbial symbiosis across anatomy during the evolution of a basal animal lineage., Associations between corals and symbiotic microorganisms could be driven by the environment or shared evolutionary history. Here, the authors examine relationships between coral phylogenies and associated microbiomes, finding evidence of phylosymbiosis in microbes from coral skeleton and tissue, but not mucus.

In numerous animal clades, the evolutionary history of host species drives patterns of gut microbial community structure, resulting in more divergent microbiota with increasing phylogenetic distance between hosts. This phenomenon, termed phylosymbiosis, has been observed in diverse evolutionary lineages, but has been difficult to detect in birds. Previous tests of phylosymbiosis among birds have been conducted using wild individuals, and thus interspecific differences in diet and environment may have masked a phylogenetic signal. Therefore, we tested for phylosymbiosis among all 15 species of cranes (family Gruidae) housed in the same captive environment and maintained on identical diets. 16S rRNA sequencing revealed that crane species harbour distinct gut microbiota. Overall, we detected marginally significant patterns of phylosymbiosis, the strength of which was increased when including the estimates of absolute microbial abundance (rather than relative abundance) derived from microbial densities determined by flow cytometry. Using this approach, we detected the statistically significant signatures of phylosymbiosis only after removing male cranes from our analysis, suggesting that using mixed-sex animal cohorts may prevent the detection of phylosymbiosis. Though weak compared with mammals (and especially insects), these results provide evidence of phylosymbiosis in birds. We discuss the potential differences between birds and mammals, such as transmission routes and host filtering, that may underlie the differences in the strength of phylosymbiosis.

In numerous animal clades, the evolutionary history of host species drives patterns of gut microbial community structure, resulting in more divergent microbiota with increased phylogenetic distance between hosts. This phenomenon, termed phylosymbiosis, has been observed in diverse evolutionary lineages, but has been difficult to detect in birds. Previous tests of phylosymbiosis among birds have been conducted using wild individuals, and thus interspecific differences in diet and environment may have masked a phylogenetic signal. Therefore, we tested for phylosymbiosis among all 15 species of cranes (family Gruidae) housed in the same captive environment and maintained on identical diets. 16S rRNA sequencing revealed that crane species harbour distinct gut microbiota. Overall, we detected the marginally significant patterns of phylosymbiosis, the strength of which was increased when including the estimates of absolute microbial abundance (rather than relative abundance) derived from microbial densities determined by flow cytometry. Using this approach, we detected the statistically significant signatures of phylosymbiosis only after removing male cranes from our analysis, suggesting that using mixed-sex animal cohorts may prevent the detection of phylosymbiosis. Though weak compared to mammals (and especially insects), these results provide the evidence of phylosymbiosis in birds. We discuss the potential differences between birds and mammals, such as transmission routes and host filtering, that may underlie the differences in the strength of phylosymbiosis.

Phylosymbiosis is a cross-system trend whereby microbial community relationships recapitulate the host phylogeny. In Nasonia parasitoid wasps, phylosymbiosis occurs throughout development, is distinguishable between sexes, and benefits host development and survival. Moreover, the microbiome shifts in hybrids as a rare Proteus bacteria in the microbiome becomes dominant. The larval hybrids then catastrophically succumb to bacterial-assisted lethality and reproductive isolation between the species. Two important questions for understanding phylosymbiosis and bacterial-assisted lethality in hybrids are: (i) Do the Nasonia bacterial genomes differ from other animal isolates and (ii) Are the hybrid bacterial genomes the same as those in the parental species? Here we report the cultivation, whole genome sequencing, and comparative analyses of the most abundant gut bacteria in Nasonia larvae, Providencia rettgeri and Proteus mirabilis. Characterization of new isolates shows Proteus mirabilis forms a more robust biofilm than Providencia rettgeri and when grown in co-culture, Proteus mirabilis significantly outcompetes Providencia rettgeri. Providencia rettgeri genomes from Nasonia are similar to each other and more divergent to pathogenic, human-associates strains. Proteus mirabilis from N. vitripennis, N. giraulti, and their hybrid offspring are nearly identical and relatively distinct from human isolates. These results indicate that members of the larval gut microbiome within Nasonia are most similar to each other, and the strain of the dominant Proteus mirabilis in hybrids is resident in parental species. Holobiont interactions between shared, resident members of the wasp microbiome and the host underpin phylosymbiosis and hybrid breakdown.IMPORTANCEAnimal and plant hosts often establish intimate relationships with their microbiomes. In varied environments, closely-related host species share more similar microbiomes, a pattern termed phylosymbiosis. When phylosymbiosis is functionally significant and beneficial, microbial transplants between host species or host hybridization can have detrimental consequences on host biology. In the Nasonia parasitoid wasp genus that contains a phylosymbiotic gut community, both effects occur and provide evidence for selective pressures on the holobiont. Here, we show that bacterial genomes in Nasonia differ from other environments and harbor genes with unique functions that may regulate phylosymbiotic relationships. Furthermore, the bacteria in hybrids are identical to parental species, thus supporting a hologenomic tenet that the same members of the microbiome and the host genome impact phylosymbiosis, hybrid breakdown, and speciation.

Microbiomes are vital in promoting nutrient cycling and plant growth in soil ecosystems. However, microbiomes face adverse effects from multiple persistent pollutants, including per- and poly-fluoroalkyl substances (PFASs). PFASs threaten the fertility and health of soil ecosystems, yet the response of microbial community stability and trophic transfer efficiencies to PFASs is still poorly understood. This study explored the spatial patterns of PFASs in topsoil environments from the West Taihu Lake Basin of China and links their presence to soil microbial community stability at compositional and functional levels. Our results revealed that PFBA (13.87%), PFTrDA (11.63%), PFDoA (11.02%), PFOA (10.99%), and PFOS (10.39%) contributed the most to the spatial occurrence of PFASs. Soil properties, including salinity (14.47%), uniformity (9.68%), dissolved inorganic carbon (8.62%), and clay content (8.18%), affected PFASs distribution the most. In soil microbiomes, eukaryotic taxa had wider niche breadths and stronger community stability than prokaryotes when exposed to PFASs (p 0.05). The presence of PFBA and PFHpA inhibited the functional stability of archaeal and bacterial communities (p 0.05). PFBA and PFPeA reduced the structural stability of heterotrophic bacteria and Myxobacteria, respectively (p 0.05). Based on null modeling, PFPeA significantly regulated the assembly processes of most microbial sub-communities (p 0.01). The trophic transferring efficiencies of autotrophic bacteria to metazoan organisms were directly stimulated by PFASs (p 0.05), and the potential trophic transferring efficiencies of methanogenic archaea to protozoa were inhibited by PFASs (p 0.05). This study highlighted the potential contributions of PFASs to soil microbial community stability and food webs during ecological soil management.

The gut microbiome of animals, which serves important functions but can also contain potential pathogens, is to varying degrees under host genetic control. This can generate signals of phylosymbiosis, whereby gut microbiome composition matches host phylogenetic structure. However, the genetic mechanisms that generate phylosymbiosis and the scale at which they act remain unclear. Two non‐mutually exclusive hypotheses are that phylosymbiosis is driven by immunogenetic regions such as the major histocompatibility complex (MHC) controlling microbial composition, or by spatial structuring of neutral host genetic diversity via founder effects, genetic drift, or isolation by distance. Alternatively, associations between microbes and host phylogeny may be generated by their spatial autocorrelation across landscapes, rather than the direct effects of host genetics. In this study, we collected MHC, microsatellite, and gut microbiome data from separate individuals belonging to the Galápagos mockingbird species complex, which consists of four allopatrically distributed species. We applied multiple regression with distance matrices and Bayesian inference to test for correlations between average genetic and microbiome similarity across nine islands for which all three levels of data were available. Clustering of individuals by species was strongest when measured with microsatellite markers and weakest for gut microbiome distributions, with intermediate clustering of MHC allele frequencies. We found that while correlations between island‐averaged gut microbiome composition and both microsatellite and MHC dissimilarity existed across species, these relationships were greatly weakened when accounting for geographic distance. Overall, our study finds little support for large‐scale control of gut microbiome composition by neutral or adaptive genetic regions across closely related bird phylogenies, although this does not preclude the possibility that host genetics shapes gut microbiome at the individual level., This study uses the Galápagos mockingbird species complex to test the relationship between gut microbiome composition and host genetic structure at the population level based on both MHC and microsatellite genotypes. We found no relationship between microbiome composition and genetic distance based on either MHC or microsatellite markers. This suggests that, at this geographic scale, host population genetics exerts little influence on the gut microbiome.

Gut microbiota have long attracted the interest of scientists due to their profound impact on the well-being of animals. A non-random pattern of microbial assembly that results in a parallelism between host phylogeny and microbial similarity is described as phylosymbiosis. Phylosymbiosis has been consistently observed in different clades of animal hosts, but there have been no studies on crustaceans. In this study, we investigated whether host phylogeny has an impact on the gut microbiota assemblages in decapod shrimps. We examined the gut microbial communities in 20 shrimp species from three families inhabiting distinct environments, using metabarcoding analyses of the V1-V3 hypervariable region of the 16S rRNA gene. Gut microbial communities varied within each shrimp group but were generally dominated by Proteobacteria. A prevalent phylosymbiotic pattern in shrimps was evidenced for the first time by the observations of (1) the distinguishability of microbial communities among species within each group, (2) a significantly lower intraspecific than interspecific gut microbial beta diversity across shrimp groups, (3) topological congruence between host phylogenetic trees and gut microbiota dendrograms, and (4) a correlation between host genetic distances and microbial dissimilarities. Consistent signals of phylosymbiosis were observed across all groups in dendrograms based on the unweighted UniFrac distances at 99% operational taxonomic units (OTUs) level and in Mantel tests based on the weighted UniFrac distances based on 97% OTUs and amplicon sequence variants. Penaeids exhibited phylosymbiosis in most tests, while phylosymbiotic signals in atyids and pandalids were only detected in fewer than half of the tests. A weak phylogenetic signal was detected in the predicted functions of the penaeid gut microbiota. However, the functional diversities of the two caridean groups were not significantly related to host phylogeny. Our observations of a parallelism in the taxonomy of the gut microbiota with host phylogeny for all shrimp groups examined and in the predicted functions for the penaeid shrimps indicate a tight host-microbial relationship during evolution.

Phylosymbiosis, where similarities in host-associated microbial communities recapitulate the phylogeny of their hosts, is a newly recognized yet pervasive pattern in the field of host–microbe interactions. While phylosymbiosis has been documented across many systems, we still have a poor understanding of the mechanisms that underlie this emergent pattern. Host selection of the microbiome is a widely cited mechanism, yet other basic ecological and evolutionary processes (dispersal, drift and diversification) may also be at play. This paper discusses the roles that each of these processes and their interactions may play in yielding phylosymbiotic signals across hosts. Finally, this paper will identify open questions and methods that are required to better understand the relative contributions of these basic processes to phylosymbiosis. Given that phylosymbiosis has been shown to relate to functional components of host fitness, understanding the processes that contribute to these patterns will be important for our understanding of the ecology and evolution of host–microbe interactions. This article is part of the theme issue ‘Conceptual challenges in microbial community ecology’.

The microbiome of plants and animals are often correlated with host phylogeny; an eco-evolutionary pattern known as phylosymbiosis. While the underlying mechanisms of this pattern are still debated, phylosymbiosis is consistently being observed as more complex hosts are studied. Using 16S rRNA sequencing to profile the microbiome, paired with COI and 18S sequences for host phylogeny, we show that diverse groups of coral reef invertebrates encompassing scleractinian coral, octocoral, sponges and ascidians exhibit phylosymbiosis.

Cyanobacteria from the genus Microcystis can form large mucilaginous colonies with attached heterotrophic bacteria – their microbiome. However, the nature of the relationship between Microcystis and its microbiome remains unclear. Is it a long-term, evolutionarily stable association? Which partners benefit? Here we report the genomic diversity of 109 individual Microcystis colonies – including cyanobacteria and associated bacterial genomes – isolated in situ and without culture from Lake Champlain, Canada and Pampulha Reservoir, Brazil. We found 14 distinct Microcystis genotypes from Canada, of which only two have been previously reported, and four genotypes specific to Brazil. Microcystis genetic diversity was much greater between than within colonies, consistent with colony growth by clonal expansion rather than aggregation of Microcystis cells. We also identified 72 bacterial species in the microbiome. Each Microcystis genotype had a distinct microbiome composition, and more closely-related genotypes had more similar microbiomes. This pattern of phylosymbiosis could be explained by co-phylogeny in two out of the nine most prevalent associated bacterial genera, Roseomonas and Rhodobacter, suggesting long-term evolutionary associations. Roseomonas and Rhodobacter genomes encode functions which could complement the metabolic repertoire of Microcystis, such as cobalamin and carotenoid biosynthesis, and nitrogen fixation. In contrast, other colony-associated bacteria showed weaker signals of co-phylogeny, but stronger evidence of horizontal gene transfer with Microcystis. These observations suggest that acquired genes are more likely to be retained in both partners (Microcystis and members of its microbiome) when they are loosely associated, whereas one gene copy is sufficient when the association is physically tight and evolutionarily long-lasting.

Phylosymbiosis was recently formulated to support a hypothesis-driven framework for the characterization of a new, cross-system trend in host-associated microbiomes. Defining phylosymbiosis as “microbial community relationships that recapitulate the phylogeny of their host”, we review the relevant literature and data in the last decade, emphasizing frequently used methods and regular patterns observed in analyses. Quantitative support for phylosymbiosis is provided by statistical methods evaluating higher microbiome variation between host species than within host species, topological similarities between the host phylogeny and microbiome dendrogram, and a positive association between host genetic relationships and microbiome beta diversity. Significant degrees of phylosymbiosis are prevalent, but not universal, in microbiomes of plants and animals from terrestrial and aquatic habitats. Consistent with natural selection shaping phylosymbiosis, microbiome transplant experiments demonstrate reduced host performance and/or fitness upon host-microbiome mismatches. Hybridization can also disrupt phylosymbiotic microbiomes and cause hybrid pathologies. The pervasiveness of phylosymbiosis carries several important implications for advancing knowledge of eco-evolutionary processes that impact host-microbiome interactions and future applications of precision microbiology. Important future steps will be to examine phylosymbiosis beyond bacterial communities, apply evolutionary modeling for an increasingly sophisticated understanding of phylosymbiosis, and unravel the host and microbial mechanisms that contribute to the pattern. This review serves as a gateway to experimental, conceptual, and quantitative themes of phylosymbiosis and outlines opportunities ripe for investigations from a diversity of disciplines.

Phylosymbiosis refers to a congruent pattern between the similarity of microbiomes of different species and the branching pattern of the host phylogeny. Phylosymbiosis has been detected in a variety of vertebrate and invertebrate hosts, but has only been assessed in geographically isolated populations. We tested for phylosymbiosis in eight (sub)species of western chipmunks with overlapping ranges and ecological niches; we used a nuclear (Acrosin) and a mitochondrial (CYTB) phylogenetic marker because there are many instances of mitochondrial introgression in chipmunks. We predicted that similarity among microbiomes increases with: (1) increasing host mitochondrial relatedness, (2) increasing host nuclear genome relatedness and (3) decreasing geographic distance among hosts. We did not find statistical evidence supporting phylosymbiosis in western chipmunks. Furthermore, in contrast to studies of other mammalian microbiomes, similarity of chipmunk microbiomes is not predominantly determined by host species. Sampling site explained most variation in microbiome composition, indicating an important role of local environment in shaping microbiomes. Fecal microbiomes of chipmunks were dominated by Bacteroidetes (72.2%), followed by Firmicutes (24.5%), which is one of the highest abundances of Bacteroidetes detected in wild mammals. Future work will need to elucidate the effects of habitat, ecology and host genomics on chipmunk microbiomes.

Host-associated microbial communities consist of stable and transient members that can assemble through purely stochastic processes associated with the environment or by interactions with the host. Phylosymbiosis predicts that if host-microbiota interactions impact assembly patterns, then one conceivable outcome is concordance between host evolutionary histories (phylogeny) and the ecological similarities in microbial community structures (microbiota dendrogram). This assembly pattern has been demonstrated in several clades of animal hosts in laboratory and natural populations, but in vertebrates, it has only been investigated using samples from faeces or the distal colon. Here, we collected the contents of five gut regions from seven rodent species and inventoried the bacterial communities by sequencing the 16S rRNA gene. We investigated how community structures varied across gut regions and whether the pattern of phylosymbiosis was present along the length of the gut. Gut communities varied by host species and gut region, with Oscillospira and Ruminococcus being more abundant in the stomach and hindgut regions. Gut microbial communities were highly distinguishable by host species across all gut regions, with the strength of the discrimination increasing along the length of the gut. Last, the pattern of phylosymbiosis was found in all five gut regions, as well as faeces. Aspects of the gut environment, such as oxygen levels, production of antimicrobials or other factors, may shift microbial communities across gut regions. However, regardless of these differences, host species maintain distinguishable, phylosymbiotic assemblages of microbes that may have functional impacts for the host.

Additional file 4. Supplementary figures.

Additional file 1. Supplementary methods. DNA extraction colonies.

Additional file 19: Table S15. Colonies sequenced per each sampling time point.

Additional file 10: Figure S4. Phylosymbiosis between Widow Phylogeny and Microbiota Dendrogram Relationships. Congruency between widow related spider phylogeny and beta-diversity dendrograms was measured via Robinson-Foulds Cluster (R-F Cluster) metric, where a value of 0 = complete congruence and 2 = incomplete incongruence. Matching Cluster (MC) Metric also accounts for incongruency between closely related branches of rooted trees and a value of 0 also indicates complete congruence between compared trees. Pink asterisks are utilized to denote phylosymbiotic relationships between the widow phylogeny and microbial dendrogram.

Perfluorinated compounds (PFCs) pose serious threats to aquatic ecosystems, especially their microbial communities. However, little is known about the phylosymbiosis of aquatic fungal and viridiplantae communities in response to PFC accumulation. We quantified the distribution of 14 PFCs in rivers and found that PFBA was dominant in the transition from water to sediment. High through-put sequencing revealed that phyla Ascomycota, Basidiomycota, Anthophyta, and Chlorophyta were the predominant in eukaryotic community. The effects of PFCs on spatial community coalescence at taxonomic and phylogenetic levels (p 0.05) were revealed. Fungal community coalescence triggered the spatial assembly of fungal and viridiplantae communities in riverine environments (p 0.05). Null modeling indicated that PFBA, PFTrDA and PFOS, etc, mediated phylogenetic assembly (p 0.05) and stochastic processes (86.67-100%) maintain phylogenetic turnover in the fungal community. Meanwhile, variable selection (27.78-54.44%) explained the viridiplantae community assemblage. Finally, we identified fungal genera Hannaella, Naganishia, Purpureocillium and Stachybotrys as indicators for PFC pollution (p 0.001). These results help explain the effects of PFCs on riverine ecological remediation.

Additional file 4: Table S3. Comparative mock community composition.

Additional file 11: Table S6. Microbial transcripts detected in RNA sequencing datasets from silk, venom and/ or ovary glands.

Additional file 1: Figure S1. Study Design Flowchart. Overview of methods utilized to isolate, evaluate, and compare the microbiomes of our host spider species and their tissue samples.

Additional file 3: Table S2. Sequencing Primers.

Additional file 2: Table S1. PCR 16S rRNA gene V1-V2 Primer Sequences.

Additional file 5: Figure S2. QIIME 2 Microbial Alpha-Diversity Plots. Alpha-Rarefaction results for Faith-PD (a), Observed OTUs (b), and Shannon Index (c) with Depth of Coverage at 8000 rarified sequences.

Additional file 7: Figure S4. QIIME 2 Alpha Diversity Group Significance Testing. Box plots of alphadiversity tests grouped by base sample type. Overall Kruskal-Wallis p-value = 0.00028 (for all groups). (a) Observed OTUs, (b) Faith Phylogenetic-Diversity Index, (c) Shannon Index, (d) Pielou Evenness Index.

Additional file 1 Supplemental Table 1. All pairwise comparisons among species across all community dimensions. P (perm) is the calculated p-value based on permuted values. BC similarity is Bray-Curtis similarity. Phylogenetic distance is the KR distance method.

Additional file 9: Table S5. Listings of ovary gland (a) and fat tissue (b) specific microbiota.

Additional file 6: Figure S3. Taxa Plots per Tissue Type – Level 3, Class. Major microbial constituents by class (> 2%) within each whole spider and tissue sample for each species (replicates grouped via mean-ceiling). Sample types as follows: Whole Spiders, Cephalothorax, Venom Glands, Ovary, Silk Glands, and Fat Tissue.

Additional file 8: Tables S4. Listings of venom gland (a) and silk gland (b) specific microbiota.