Host and tissue variations overshadow the response of boreal moss‐associated fungal communities to increased nitrogen load

Human activity has more than doubled the amount of nitrogen entering the global nitrogen cycle, and the boreal forest biome is a nitrogen‐limited ecosystem sensitive to nitrogen load perturbation. Although bryophyte‐associated microbes contribute significantly to boreal forest ecosystem function, particularly in carbon and nitrogen cycling, little is known about their responses to anthropogenic global change. Amplicon pyrosequencing of the ITS2 region of rDNA was used to investigate how fungal communities associated with three bryophyte species responded to increased nitrogen loads in a long‐term fertilization experiment in a boreal Picea abies forest in southern Norway. Overall, OTU richness, community composition and the relative abundance of specific ecological guilds were primarily influenced by host species identity and tissue type. Although not the primary factor affecting fungal communities, nitrogen addition did impact the abundance of specific guilds of fungi and the resulting overall community composition. Increased nitrogen loads decreased ectomycorrhizal abundance, with Amphinema, Cortinarius, Russula and Tylospora OTUs responding negatively to fertilization. Pathogen abundance increased with fertilization, particularly in the moss pathogen Eocronartium. Saprophytic fungi were both positively and negatively impacted by the nitrogen addition, indicating a complex community level response. The overshadowing of the effects of increased nitrogen loads by variation related to host and tissue type highlights the complexity of bryophyte‐associated microbial communities and the intricate nature of their responses to anthropogenic global change.     

Coexisting cryptic species of the Litoditis marina complex (Nematoda) show differential resource use and have distinct microbiomes with high intraspecific variability

Differences in resource use or in tolerances to abiotic conditions are often invoked as potential mechanisms underlying the sympatric distribution of cryptic species. Additionally, the microbiome can provide physiological adaptations of the host to environmental conditions. We determined the intra‐ and interspecific variability of the microbiomes of three cryptic nematode species of the Litoditis marina species complex that co‐occur, but show differences in abiotic tolerances. Roche 454 pyrosequencing of the microbial 16S rRNA gene revealed distinct bacterial communities characterized by a substantial diversity (85–513 OTUs) and many rare OTUs. The core microbiome of each species contained only very few OTUs (2–6), and four OTUs were identified as potentially generating tolerance to abiotic conditions. A controlled experiment in which nematodes from two cryptic species (Pm1 and Pm3) were fed with either an E. coli suspension or a bacterial mix was performed, and the 16S rRNA gene was sequenced using the MiSeq technology. OTU richness was 10‐fold higher compared to the 454 data set and ranged between 1118 and 7864. This experiment confirmed the existence of species‐specific microbiomes, a core microbiome with few OTUs, and high interindividual variability. The offered food source affected the bacterial community and illustrated different feeding behaviour between the cryptic species, with Pm3 exhibiting a higher degree of selective feeding than Pm1. Morphologically similar species belonging to the same feeding guild (bacterivores) can thus have substantial differences in their associated microbiomes and feeding strategy, which in turn may have important ramifications for biodiversity–ecosystem functioning relationships.     

A common garden experiment with Porphyra umbilicalis (Rhodophyta) evaluates methods to study spatial differences in the macroalgal microbiome

While macroalgal microbiomes are the focus of many recent studies, there is little information about microbial spatial diversity across the thallus. Reliance on field material makes it difficult to discern whether recovered microbiomes belong to the host or its epiphytes, and technical comparisons of macroalgal samples for microbial studies are needed. Here, we use a common garden approach that avoids the problem of epiphytes, particularly at holdfasts, to examine the microbiome of Porphyra umbilicalis (strain Pum1). We used the V6 hypervariable region of the 16S rDNA with Illumina HiSeq sequencing and developed PNA clamps to block recovery of organelle V6 sequences. The common garden approach allowed us to determine differences in the microbiome at the holdfast versus blade margin. We found a notable increase in the relative abundance of Planctomycetes and Alphaproteobacteria at the holdfast, particularly of the possible symbiont Sulfitobacter sp. Nonadjacent 1.5 cm² samples of blade margin had microbiomes that were not statistically different. The most abundant phylum in the overall microbiome was Proteobacteria, followed by Bacteroidetes. Because phycologists often work in remote sites, we compared three stabilization and preparation techniques and found silica gel desiccation/bead‐beating and flash‐freezing/lyophilization/bead‐beating to be interchangeable. Core taxa (≥0.1% of sequences) across treatments were similar and accounted for ≥95% of all sequences. Finally, statistical conclusions for all comparisons were the same, regardless of which microbial community analysis tool was used: mothur or minimum entropy decomposition.     

Amphibian skin fungal communities vary across host species and do not correlate with infection by a pathogenic fungus

Amphibian population declines caused by the fungus Batrachochytrium dendrobatidis (Bd) have prompted studies on the bacterial community that resides on amphibian skin. However, studies addressing the fungal portion of these symbiont communities have lagged behind. Using ITS1 amplicon sequencing, we examined the fungal portion of the skin microbiome of temperate and tropical amphibian species currently coexisting with Bd in nature. We assessed cooccurrence patterns between bacterial and fungal OTUs using a subset of samples for which bacterial 16S rRNA gene amplicon data were also available. We determined that fungal communities were dominated by members of the phyla Ascomycota and Basidiomycota, and also by Chytridiomycota in the most aquatic amphibian species. Alpha diversity of the fungal communities differed across host species, and fungal community structure differed across species and regions. However, we did not find a correlation between fungal diversity/community structure and Bd infection, though we did identify significant correlations between Bd and specific OTUs. Moreover, positive bacterial–fungal cooccurrences suggest that positive interactions between these organisms occur in the skin microbiome. Understanding the ecology of amphibian skin fungi, and their interactions with bacteria will complement our knowledge of the factors influencing community assembly and the overall function of these symbiont communities.     

Fungal infection alters the selection,dispersal and drift processes structuring the amphibian skin microbiome

Symbiotic microbial communities are important for host health, but the processes shaping these communities are poorly understood. Understanding how community assembly processes jointly affect microbial community composition is limited because inflexible community models rely on rejecting dispersal and drift before considering selection. We developed a flexible community assembly model based on neutral theory to ask: How do dispersal, drift and selection concurrently affect the microbiome across environmental gradients? We applied this approach to examine how a fungal pathogen affected the assembly processes structuring the amphibian skin microbiome. We found that the rejection of neutrality for the amphibian microbiome across a fungal gradient was not strictly due to selection processes, but was also a result of species‐specific changes in dispersal and drift. Our modelling framework brings the qualitative recognition that niche and neutral processes jointly structure microbiomes into quantitative focus, allowing for improved predictions of microbial community turnover across environmental gradients.     

Host population genetics and biogeography structure the microbiome of the sponge Cliona delitrix

Sponges occur across diverse marine biomes and host internal microbial communities that can provide critical ecological functions. While strong patterns of host specificity have been observed consistently in sponge microbiomes, the precise ecological relationships between hosts and their symbiotic microbial communities remain to be fully delineated. In the current study, we investigate the relative roles of host population genetics and biogeography in structuring the microbial communities hosted by the excavating sponge Cliona delitrix. A total of 53 samples, previously used to demarcate the population genetic structure of C. delitrix, were selected from two locations in the Caribbean Sea and from eight locations across the reefs of Florida and the Bahamas. Microbial community diversity and composition were measured using Illumina‐based high‐throughput sequencing of the 16S rRNA V4 region and related to host population structure and geographic distribution. Most operational taxonomic units (OTUs) specific to Cliona delitrix microbiomes were rare, while other OTUs were shared with congeneric hosts. Across a large regional scale (>1,000 km), geographic distance was associated with considerable variability of the sponge microbiome, suggesting a distance–decay relationship, but little impact over smaller spatial scales (<300 km) was observed. Host population structure had a moderate effect on the structure of these microbial communities, regardless of geographic distance. These results support the interplay between geographic, environmental, and host factors as forces determining the community structure of microbiomes associated with C. delitrix. Moreover, these data suggest that the mechanisms of host regulation can be observed at the population genetic scale, prior to the onset of speciation.     

Temporal Variation in Infection Prevalence by the Amphibian Chytrid Fungus in Three Species of Frogs at La Selva,Costa Rica

The emerging infectious disease chytridiomycosis, caused by the amphibian chytrid fungus, Batrachochytrium dendrobatidis (Bd), is implicated in widespread population declines, extirpations, and extinctions of amphibians throughout the world. In the Neotropics, most amphibian declines have occurred in cool mid‐ to high‐elevation sites (> 400 m asl), and it is hypothesized that high temperatures limit the growth of Bd in lowland tropical sites, despite few data available on the distribution of Bd in lowland forests. Here, we report the results of a 12‐mo pathogen surveillance program for three common species of frogs at a warm lowland site in northeastern Costa Rica. We combine standard non‐invasive skin swabbing techniques with a quantitative polymerase chain reaction assay to analyze the infection prevalence and Bd load across a 1‐yr period. Our data indicate an overall Bd infection rate of 6.1 percent, but prevalence varies from < 5 percent in warmer months to a peak of 34.7 percent in the coolest months of the year. Despite very little seasonal variation in temperature (< 4°C), our data indicate strong seasonal variation in the prevalence of Bd, with highest prevalence of infection in months with coolest air temperatures. While it has been suggested that Bd is primarily a riparian fungus, we find no difference in prevalence of infection among our species despite considerable differences in affiliation of these species with water. Our study provides further evidence that infection by Bd is regulated by temperature and shows that warm temperatures in lowland forests may restrict, but not prevent, infection by Bd.     

Characterization of the cyanobacteria and associated bacterial community from an ephemeral wetland in New Zealand

New Zealand ephemeral wetlands are ecologically important, containing up to 12% of threatened native plant species and frequently exhibiting conspicuous cyanobacterial growth. In such environments, cyanobacteria and associated heterotrophs can influence primary production and nutrient cycling. Wetland communities, including bacteria, can be altered by increased nitrate and phosphate due to agricultural practices. We have characterized cyanobacteria from the Wairepo Kettleholes Conservation Area and their associated bacteria. Use of 16S rRNA amplicon sequencing identified several operational taxonomic units (OTUs) representing filamentous heterocystous and non‐heterocystous cyanobacterial taxa. One Nostoc OTU that formed macroscopic colonies dominated the cyanobacterial community. A diverse bacterial community was associated with the Nostoc colonies, including a core microbiome of 39 OTUs. Identity of the core microbiome associated with macroscopic Nostoc colonies was not changed by the addition of nutrients. One OTU was highly represented in all Nostoc colonies (27.6%–42.6% of reads) and phylogenetic analyses identified this OTU as belonging to the genus Sphingomonas. Scanning electron microscopy showed the absence of heterotrophic bacteria within the Nostoc colony but revealed a diverse community associated with the colonies on the external surface.     

Genetic variability,community structure,and horizontal transfer of endosymbionts among three Asia II‐Bemisia tabaci mitotypes in Pakistan

Endosymbionts associated with the whitefly Bemisia tabaci cryptic species are known to contribute to host fitness and environmental adaptation. The genetic diversity and population complexity were investigated for endosymbiont communities of B. tabaci occupying different micro‐environments in Pakistan. Mitotypes of B. tabaci were identified by comparative sequence analysis of the mitochondria cytochrome oxidase I (mtCOI) gene sequence. Whitefly mitotypes belonged to the Asia II‐1, ‐5, and ‐7 mitotypes of the Asia II major clade. The whitefly–endosymbiont communities were characterized based on 16S ribosomal RNA operational taxonomic unit (OTU) assignments, resulting in 43 OTUs. Most of the OTUs occurred in the Asia II‐1 and II‐7 mitotypes (r² = .9, p < .005), while the Asia II‐5 microbiome was less complex. The microbiome OTU groups were mitotype‐specific, clustering with a basis in phylogeographical distribution and the corresponding ecological niche of their whitefly host, suggesting mitotype‐microbiome co‐adaptation. The primary endosymbiont Portiera was represented by a single, highly homologous OTU (0%–0.67% divergence). Two of six Arsenophonus OTUs were uniquely associated with Asia II‐5 and ‐7, and one occurred exclusively in Asia II‐1, two only in Asia II‐5, and one in both Asia II‐1 and ‐7. Four other secondary endosymbionts, Cardinium, Hemipteriphilus, Rickettsia, and Wolbachia OTUs, were found at ≤29% frequencies. The most prevalent Arsenophonus OTU was found in all three Asia II mitotypes (55% frequency), whereas the same strain of Cardinium and Wolbachia was found in both Asia II‐1 and ‐5, and a single Hemipteriphilus OTU occurred in Asia II‐1 and ‐7. This pattern is indicative of horizontal transfer, suggestive of a proximity between mitotypes sufficient for gene flow at overlapping mitotype ecological niches.     

Divergence thresholds and divergent biodiversity estimates: can metabarcoding reliably describe zooplankton communities?

DNA metabarcoding is a promising method for describing communities and estimating biodiversity. This approach uses high‐throughput sequencing of targeted markers to identify species in a complex sample. By convention, sequences are clustered at a predefined sequence divergence threshold (often 3%) into operational taxonomic units (OTUs) that serve as a proxy for species. However, variable levels of interspecific marker variation across taxonomic groups make clustering sequences from a phylogenetically diverse dataset into OTUs at a uniform threshold problematic. In this study, we use mock zooplankton communities to evaluate the accuracy of species richness estimates when following conventional protocols to cluster hypervariable sequences of the V4 region of the small subunit ribosomal RNA gene (18S) into OTUs. By including individually tagged single specimens and “populations” of various species in our communities, we examine the impact of intra‐ and interspecific diversity on OTU clustering. Communities consisting of single individuals per species generated a correspondence of 59–84% between OTU number and species richness at a 3% divergence threshold. However, when multiple individuals per species were included, the correspondence between OTU number and species richness dropped to 31–63%. Our results suggest that intraspecific variation in this marker can often exceed 3%, such that a single species does not always correspond to one OTU. We advocate the need to apply group‐specific divergence thresholds when analyzing complex and taxonomically diverse communities, but also encourage the development of additional filtering steps that allow identification of artifactual rRNA gene sequences or pseudogenes that may generate spurious OTUs.     

Microbiome responses during virulence adaptation by a phloem‐feeding insect to resistant near‐isogenic rice lines

The microbiomes of phloem‐feeding insects include functional bacteria and yeasts essential for herbivore survival and development. Changes in microbiome composition are implicated in virulence adaptation by herbivores to host plant species or host populations (including crop varieties). We examined patterns in adaptation by the green leafhopper, Nephotettix virescens, to near‐isogenic rice lines (NILs) with one or two resistance genes and the recurrent parent T65, without resistance genes. Only the line with two resistance genes was effective in reducing leafhopper fitness. After 20 generations on the resistant line, selected leafhoppers attained similar survival, weight gain, and egg laying to leafhoppers that were continually reared on the susceptible recurrent parent, indicating that they had adapted to the resistant host. By sequencing the 16s rRNA gene, we described the microbiome of leafhoppers from colonies associated with five collection sites, and continually reared or switched between NILs. The microbiomes included 69–119 OTUs of which 44 occurred in ≥90% of samples. Of these, 14 OTUs were assigned to the obligate symbiont Candidatus sulcia clade. After 20 generations of selection, collection site had a greater effect than host plant on microbiome composition. Six bacteria genera, including C. sulcia, were associated with leafhopper virulence. However, there was significant within‐treatment, site‐related variability in the prevalence of these taxa such that the mechanisms underlying their association with virulence remain to be determined. Our results imply that these taxa are associated with leafhopper nutrition. Ours is the first study to describe microbiome diversity and composition in rice leafhoppers. We discuss our results in light of the multiple functions of herbivore microbiomes during virulence adaptation in insect herbivores.     

Resistance,tolerance and environmental transmission dynamics determine host extinction risk in a load‐dependent amphibian disease

While disease‐induced extinction is generally considered rare, a number of recently emerging infectious diseases with load‐dependent pathology have led to extinction in wildlife populations. Transmission is a critical factor affecting disease‐induced extinction, but the relative importance of transmission compared to load‐dependent host resistance and tolerance is currently unknown. Using a combination of models and experiments on an amphibian species suffering extirpations from the fungal pathogen Batrachochytrium dendrobatidis (Bd), we show that while transmission from an environmental Bd reservoir increased the ability of Bd to invade an amphibian population and the extinction risk of that population, Bd‐induced extinction dynamics were far more sensitive to host resistance and tolerance than to Bd transmission. We demonstrate that this is a general result for load‐dependent pathogens, where non‐linear resistance and tolerance functions can interact such that small changes in these functions lead to drastic changes in extinction dynamics.     

Community assembly of a euryhaline fish microbiome during salinity acclimation

Microbiomes play a critical role in promoting a range of host functions. Microbiome function, in turn, is dependent on its community composition. Yet, how microbiome taxa are assembled from their regional species pool remains unclear. Many possible drivers have been hypothesized, including deterministic processes of competition, stochastic processes of colonization and migration, and physiological ‘host‐effect’ habitat filters. The contribution of each to assembly in nascent or perturbed microbiomes is important for understanding host–microbe interactions and host health. In this study, we characterized the bacterial communities in a euryhaline fish and the surrounding tank water during salinity acclimation. To assess the relative influence of stochastic versus deterministic processes in fish microbiome assembly, we manipulated the bacterial species pool around each fish by changing the salinity of aquarium water. Our results show a complete and repeatable turnover of dominant bacterial taxa in the microbiomes from individuals of the same species after acclimation to the same salinity. We show that changes in fish microbiomes are not correlated with corresponding changes to abundant taxa in tank water communities and that the dominant taxa in fish microbiomes are rare in the aquatic surroundings, and vice versa. Our results suggest that bacterial taxa best able to compete within the unique host environment at a given salinity appropriate the most niche space, independent of their relative abundance in tank water communities. In this experiment, deterministic processes appear to drive fish microbiome assembly, with little evidence for stochastic colonization.     

Diversity and temporal dynamics of primate milk microbiomes

Milk is inhabited by a community of bacteria and is one of the first postnatal sources of microbial exposure for mammalian young. Bacteria in breast milk may enhance immune development, improve intestinal health, and stimulate the gut‐brain axis for infants. Variation in milk microbiome structure (e.g., operational taxonomic unit [OTU] diversity, community composition) may lead to different infant developmental outcomes. Milk microbiome structure may depend on evolutionary processes acting at the host species level and ecological processes occurring over lactation time, among others. We quantified milk microbiomes using 16S rRNA high‐throughput sequencing for nine primate species and for six primate mothers sampled over lactation. Our data set included humans (Homo sapiens, Philippines and USA) and eight nonhuman primate species living in captivity (bonobo [Pan paniscus], chimpanzee [Pan troglodytes], western lowland gorilla [Gorilla gorilla gorilla], Bornean orangutan [Pongo pygmaeus], Sumatran orangutan [Pongo abelii], rhesus macaque [Macaca mulatta], owl monkey [Aotus nancymaae]) and in the wild (mantled howler monkey [Alouatta palliata]). For a subset of the data, we paired microbiome data with nutrient and hormone assay results to quantify the effect of milk chemistry on milk microbiomes. We detected a core primate milk microbiome of seven bacterial OTUs indicating a robust relationship between these bacteria and primate species. Milk microbiomes differed among primate species with rhesus macaques, humans and mantled howler monkeys having notably distinct milk microbiomes. Gross energy in milk from protein and fat explained some of the variations in microbiome composition among species. Microbiome composition changed in a predictable manner for three primate mothers over lactation time, suggesting that different bacterial communities may be selected for as the infant ages. Our results contribute to understanding ecological and evolutionary relationships between bacteria and primate hosts, which can have applied benefits for humans and endangered primates in our care.     

Patterns of Symbiodinium (Dinophyceae) diversity and assemblages among diverse hosts and the coral reef environment of Lizard Island,Australia

Large‐scale environmental disturbances may impact both partners in coral host–Symbiodinium systems. Elucidation of the assembly patterns in such complex and interdependent communities may enable better prediction of environmental impacts across coral reef ecosystems. In this study, we investigated how the community composition and diversity of dinoflagellate symbionts in the genus Symbiodinium were distributed among 12 host species from six taxonomic orders (Actinaria, Alcyonacea, Miliolida, Porifera, Rhizostoma, Scleractinia) and in the reef water and sediments at Lizard Island, Great Barrier Reef before the 3rd Global Coral Bleaching Event. 454 pyrosequencing of the ITS2 region of Symbiodinium yielded 83 operational taxonomic units (OTUs) at a 97% similarity cut‐off. Approximately half of the Symbiodinium OTUs from reef water or sediments were also present in symbio. OTUs belonged to six clades (A‐D, F‐G), but community structure was uneven. The two most abundant OTUs (100% matches to types C1 and A3) comprised 91% of reads and OTU C1 was shared by all species. However, sequence‐based analysis of these dominant OTUs revealed host species specificity, suggesting that genetic similarity cut‐offs of Symbiodinium ITS2 data sets need careful evaluation. Of the less abundant OTUs, roughly half occurred at only one site or in one species and the background Symbiodinium communities were distinct between individual samples. We conclude that sampling multiple host taxa with differing life history traits will be critical to fully understand the symbiont diversity of a given system and to predict coral ecosystem responses to environmental change and disturbance considering the differential stress response of the taxa within.     

Manipulation of the seagrass-associated microbiome reduces disease severity

Host-associated microbes influence host health and function and can be a first line of defence against infections. While research increasingly shows that terrestrial plant microbiomes contribute to bacterial, fungal, and oomycete disease resistance, no comparable experimental work has investigated marine plant microbiomes or more diverse disease agents. We test the hypothesis that the eelgrass (Zostera marina) leaf microbiome increases resistance to seagrass wasting disease. From field eelgrass with paired diseased and asymptomatic tissue, 16S rRNA gene amplicon sequencing revealed that bacterial composition and richness varied markedly between diseased and asymptomatic tissue in one of the two years. This suggests that the influence of disease on eelgrass microbial communities may vary with environmental conditions. We next experimentally reduced the eelgrass microbiome with antibiotics and bleach, then inoculated plants with Labyrinthula zosterae, the causative agent of wasting disease. We detected significantly higher disease severity in eelgrass with a native microbiome than an experimentally reduced microbiome. Our results over multiple experiments do not support a protective role of the eelgrass microbiome against L. zosterae. Further studies of these marine host–microbe–pathogen relationships may continue to show new relationships between plant microbiomes and diseases.     

Gut microbiomes of free‐ranging and captive Namibian cheetahs: Diversity,putative functions and occurrence of potential pathogens

Although the significance of the gut microbiome for host health is well acknowledged, the impact of host traits and environmental factors on the interindividual variation of gut microbiomes of wildlife species is not well understood. Such information is essential; however, as changes in the composition of these microbial communities beyond the natural range might cause dysbiosis leading to increased susceptibility to infections. We examined the potential influence of sex, age, genetic relatedness, spatial tactics and the environment on the natural range of the gut microbiome diversity in free‐ranging Namibian cheetahs (Acinonyx jubatus). We further explored the impact of an altered diet and frequent contact with roaming dogs and cats on the occurrence of potential bacterial pathogens by comparing free‐ranging and captive individuals living under the same climatic conditions. Abundance patterns of particular bacterial genera differed between the sexes, and bacterial diversity and richness were higher in older (>3.5 years) than in younger individuals. In contrast, male spatial tactics, which probably influence host exposure to environmental bacteria, had no discernible effect on the gut microbiome. The profound resemblance of the gut microbiome of kin in contrast to nonkin suggests a predominant role of genetics in shaping bacterial community characteristics and functional similarities. We also detected various Operational Taxonomic Units (OTUs) assigned to potential pathogenic bacteria known to cause diseases in humans and wildlife species, such as Helicobacter spp., and Clostridium perfringens. Captive individuals did not differ in their microbial alpha diversity but exhibited higher abundances of OTUs related to potential pathogenic bacteria and shifts in disease‐associated functional pathways. Our study emphasizes the need to integrate ecological, genetic and pathogenic aspects to improve our comprehension of the main drivers of natural variation and shifts in gut microbial communities possibly affecting host health. This knowledge is essential for in situ and ex situ conservation management.     

Host‐derived population genomics data provides insights into bacterial and diatom composition of the killer whale skin

Recent exploration into the interactions and relationship between hosts and their microbiota has revealed a connection between many aspects of the host's biology, health and associated micro‐organisms. Whereas amplicon sequencing has traditionally been used to characterize the microbiome, the increasing number of published population genomics data sets offers an underexploited opportunity to study microbial profiles from the host shotgun sequencing data. Here, we use sequence data originally generated from killer whale Orcinus orca skin biopsies for population genomics, to characterize the skin microbiome and investigate how host social and geographical factors influence the microbial community composition. Having identified 845 microbial taxa from 2.4 million reads that did not map to the killer whale reference genome, we found that both ecotypic and geographical factors influence community composition of killer whale skin microbiomes. Furthermore, we uncovered key taxa that drive the microbiome community composition and showed that they are embedded in unique networks, one of which is tentatively linked to diatom presence and poor skin condition. Community composition differed between Antarctic killer whales with and without diatom coverage, suggesting that the previously reported episodic migrations of Antarctic killer whales to warmer waters associated with skin turnover may control the effects of potentially pathogenic bacteria such as Tenacibaculum dicentrarchi. Our work demonstrates the feasibility of microbiome studies from host shotgun sequencing data and highlights the importance of metagenomics in understanding the relationship between host and microbial ecology.     

Wood‐rotting basidiomycetes associated with declining native trees in timber‐harvesting compartments of the Garden Route National Park of South Africa

Trees in the Garden Route National Park (GRNP) indigenous forests in South Africa are selectively harvested for timber based on criteria that include signs and symptoms induced by wood‐rotting fungi. However, virtually nothing is known regarding the identity and host associations of these macro‐fungi in this natural ecosystem. Surveys were conducted in three harvesting compartments in the GRNP to investigate the taxonomic affiliation and species richness of these fungi on standing and recently harvested trees. Samples were collected from basidiomes on infected trees and tree stumps, and from diseased tissues on symptomatic trees. Phylogenetic analyses using ITS sequences characterized the isolates obtained into 26 Operational Taxonomic Units (OTUs) belonging to 17 genera after clustering the sequences at a 97% identity threshold. Ganoderma (Ganodermataceae) and Inonotus (Hymenochaetaceae) were the most species‐rich genera and the Bloukrans compartment, with 22 OTUs, showed the highest species richness. A fungus (OTU1) affiliated with Ganoderma pfeifferi was the most abundant in the surveyed areas. Its predominance was also evidenced on host trees since it occurred on 15 of the 20 tree species sampled, with Olea capensis subsp. macrocarpa (Oleaceae) being the most colonized host. Given the wide variety of wood‐rotting basidiomycetes revealed by this study and particularly the preponderance of species with pathogenic potential, more attention should be given to better understand their ecological role in this natural ecosystem as well as the effects of logging that may enhance their dissemination or negatively affect their diversity and the health of trees in the region.     

The amphibian microbiome exhibits poor resilience following pathogen-induced disturbance

Infectious pathogens can disrupt the microbiome in addition to directly affecting the host. Impacts of disease may be dependent on the ability of the microbiome to recover from such disturbance, yet remarkably little is known about microbiome recovery after disease, particularly in nonhuman animals. We assessed the resilience of the amphibian skin microbial community after disturbance by the pathogen, Batrachochytrium dendrobatidis (Bd). Skin microbial communities of laboratory-reared mountain yellow-legged frogs were tracked through three experimental phases: prior to Bd infection, after Bd infection (disturbance), and after clearing Bd infection (recovery period). Bd infection disturbed microbiome composition and altered the relative abundances of several dominant bacterial taxa. After Bd infection, frogs were treated with an antifungal drug that cleared Bd infection, but this did not lead to recovery of microbiome composition (measured as Unifrac distance) or relative abundances of dominant bacterial groups. These results indicate that Bd infection can lead to an alternate stable state in the microbiome of sensitive amphibians, or that microbiome recovery is extremely slow—in either case resilience is low. Furthermore, antifungal treatment and clearance of Bd infection had the additional effect of reducing microbial community variability, which we hypothesize results from similarity across frogs in the taxa that colonize community vacancies resulting from the removal of Bd. Our results indicate that the skin microbiota of mountain yellow-legged frogs has low resilience following Bd-induced disturbance and is further altered by the process of clearing Bd infection, which may have implications for the conservation of this endangered amphibian.Subject terms: Microbial ecology, Community ecology