High consistency between replicate 454 pyrosequencing analyses of ectomycorrhizal plant root samples

Host generalists dominate fungal communities associated with alpine knotweed roots: a study of Sebacinales

Bistorta vivipara is a widespread herbaceous perennial plant with a discontinuous pattern of distribution in arctic, alpine, subalpine and boreal habitats across the northern Hemisphere. Studies of the fungi associated with the roots of B. vivipara have mainly been conducted in arctic and alpine ecosystems. This study examined the fungal diversity and specificity from root tips of B. vivipara in two local mountain ecosystems as well as on a global scale. Sequences were generated by Sanger sequencing of the internal transcribed spacer (ITS) region followed by an analysis of accurately annotated nuclear segments including ITS1-5.8S-ITS2 sequences available from public databases. In total, 181 different UNITE species hypotheses (SHs) were detected to be fungi associated with B. vivipara, 73 of which occurred in the Bavarian Alps and nine in the Swabian Alps-with one SH shared among both mountains. In both sites as well as in additional public data, individuals of B. vivipara were found to contain phylogenetically diverse fungi, with the Basidiomycota, represented by the Thelephorales and Sebacinales, being the most dominant. A comparative analysis of the diversity of the Sebacinales associated with B. vivipara and other co-occurring plant genera showed that the highest number of sebacinoid SHs were associated with Quercus and Pinus, followed by Bistorta. A comparison of B. vivipara with plant families such as Ericaceae, Fagaceae, Orchidaceae, and Pinaceae showed a clear trend: Only a few species were specific to B. vivipara and a large number of SHs were shared with other co-occurring non-B. vivipara plant species. In Sebacinales, the majority of SHs associated with B. vivipara belonged to the ectomycorrhiza (ECM)-forming Sebacinaceae, with fewer SHs belonging to the Serendipitaceae encompassing diverse ericoid-orchid-ECM-endophytic associations. The large proportion of non-host-specific fungi able to form a symbiosis with other non-B. vivipara plants could suggest that the high fungal diversity in B. vivipara comes from an active recruitment of their associates from the co-occurring vegetation. The non-host-specificity suggests that this strategy may offer ecological advantages; specifically, linkages with generalist rather than specialist fungi. Proximity to co-occurring non-B. vivipara plants can maximise the fitness of B. vivipara, allowing more rapid and easy colonisation of the available habitats.

A Closer Examination of the 'Abundant-Center' for Ectomycorrhizal Fungal Community Associated With Picea crassifolia in China

A long-standing hypothesis in biogeography predicts that a species' abundance is highest at the center of its geographical range and decreases toward its edges. In this study, we test the abundant-center hypothesis of ectomycorrhizal (ECM) fungal communities associated with Picea crassifolia, an endemic species widely distributed in northwest China. We analyzed the taxonomic richness and the relative abundance of ECM fungi in four main distribution areas, from center to edges. In total, 234 species of ECM fungi were detected, and of these, 137 species were shared among all four sites. Inocybe, Sebacina, Tomentella, and Cortinarius were the dominant genera. ECM fungal richness and biodiversity were highest at the central and lower at peripheral sites. Our results indicated that ECM fungal species richness was consistent with the abundant-center hypothesis, while the relative abundances of individual fungal genera shifted inconsistently across the plant's range.

Seed quantity affects the fungal community composition detected using metabarcoding

Pest introductions via trade in tree seed may result from a lack of adequate survey and validation protocols. Developing better diagnostic protocols to identify potentially harmful pests and pathogens in forest tree seed is of critical importance. High-throughput sequencing-based barcoding and metabarcoding provide effective tools for screening potentially harmful organisms in various plant materials, including seeds. However, the sample size needed to detect the total microorganism diversity of a community is a major challenge in microbiome studies. In this work, we examined how increasing sample size (ranging between 100 and 1000 seeds) influences diversity of fungal communities detected by high throughput sequencing in Pinus sylvestris seeds. Our results showed that as sample size increased, fungal alpha diversity also increased. Beta-diversity estimators detected significant differences between the mycobiota from different samples. However, taxonomic and functional diversity were not correlated with sample size. In addition, we found that increasing the number of PCR replicates resulted in a higher abundance of plant pathogens. We concluded that for the purpose of screening for potentially harmful pathogens using HTS, greater efforts should be made to increase the sample size and replicates when testing tree seed.

Hidden Fungi: Combining Culture-Dependent and -Independent DNA Barcoding Reveals Inter-Plant Variation in Species Richness of Endophytic Root Fungi in Elymus repens

The root endophyte community of the grass species Elymus repens was investigated using both a culture-dependent approach and a direct amplicon sequencing method across five sites and from individual plants. There was much heterogeneity across the five sites and among individual plants. Focusing on one site, 349 OTUs were identified by direct amplicon sequencing but only 66 OTUs were cultured. The two approaches shared ten OTUs and the majority of cultured endophytes do not overlap with the amplicon dataset. Media influenced the cultured species richness and without the inclusion of 2% MEA and full-strength MEA, approximately half of the unique OTUs would not have been isolated using only PDA. Combining both culture-dependent and -independent methods for the most accurate determination of root fungal species richness is therefore recommended. High inter-plant variation in fungal species richness was demonstrated, which highlights the need to rethink the scale at which we describe endophyte communities.

Habitat specialisation controls ectomycorrhizal fungi above the treeline in the European Alps

Alpine habitats are one of the most vulnerable ecosystems to environmental change, however, little information is known about the drivers of plant-fungal interactions in these ecosystems and their resilience to climate change. We investigated the influence of the main drivers of ectomycorrhizal (EM) fungal communities along elevation and environmental gradients in the alpine zone of the European Alps and measured their degree of specialisation using network analysis. We sampled ectomycorrhizas of Dryas octopetala, Bistorta vivipara and Salix herbacea, and soil fungal communities at 28 locations across five countries, from the treeline to the nival zone. We found that: (1) EM fungal community composition, but not richness, changes along elevation, (2) there is no strong evidence of host specialisation, however, EM fungal networks in the alpine zone and within these, EM fungi associated with snowbed communities, are more specialised than in other alpine habitats, (3) plant host population structure does not influence EM fungal communities, and (4) most variability in EM fungal communities is explained by fine-scale changes in edaphic properties, like soil pH and total nitrogen. The higher specialisation and narrower ecological niches of these plant-fungal interactions in snowbed habitats make these habitats particularly vulnerable to environmental change in alpine ecosystems.

Fine-scale diversity patterns in belowground microbial communities are consistent across kingdoms

The belowground environment is heterogeneous and complex at fine spatial scales. Physical structures, biotic components and abiotic conditions create a patchwork mosaic of potential niches for microbes. Questions remain about mechanisms and patterns of community assembly belowground, including: Do fungal and bacterial communities assemble differently? How do microbes reach the roots of host plants? Within a 4 m2 plot in alpine vegetation, high throughput sequencing of the 16S (bacteria) and ITS1 (fungal) ribosomal RNA genes was used to characterise microbial community composition in roots and adjacent soil of a viviparous host plant (Bistorta vivipara). At fine spatial scales, beta-diversity patterns in belowground bacterial and fungal communities were consistent, although compositional change was greater in bacteria than fungi. Spatial structure and distance-decay relationships were also similar for bacteria and fungi, with significant spatial structure detected at <50 cm among root- but not soil-associated microbes. Recruitment of root microbes from the soil community appeared limited at this sampling and sequencing depth. Possible explanations for this include recruitment from low-abundance populations of soil microbes, active recruitment from neighbouring plants and/or vertical transmission of symbionts to new clones, suggesting varied methods of microbial community assembly for viviparous plants. Our results suggest that even at relatively small spatial scales, deterministic processes play a significant role in belowground microbial community structure and assembly.

A single ectomycorrhizal plant root system includes a diverse and spatially structured fungal community

Although only a relatively small proportion of plant species form ectomycorrhizae with fungi, it is crucial for growth and survival for a number of widespread woody plant species. Few studies have attempted to investigate the fine scale spatial structure of entire root systems of adult ectomycorrhizal (EcM) plants. Here, we use the herbaceous perennial Bistorta vivipara to map the entire root system of an adult EcM plant and investigate the spatial structure of its root-associated fungi. All EcM root tips were sampled, mapped and identified using a direct PCR approach and Sanger sequencing of the internal transcribed spacer region. A total of 32.1% of all sampled root tips (739 of 2302) were successfully sequenced and clustered into 41 operational taxonomic units (OTUs). We observed a clear spatial structuring of the root-associated fungi within the root system. Clusters of individual OTUs were observed in the younger parts of the root system, consistent with observations of priority effects in previous studies, but were absent from the older parts of the root system. This may suggest a succession and fragmentation of the root-associated fungi even at a very fine scale, where competition likely comes into play at different successional stages within the root system.

Sun Exposure Shapes Functional Grouping of Fungi in Cryptoendolithic Antarctic Communities

Antarctic cryptoendolithic microbial communities dominate ice-free areas of continental Antarctica, among the harshest environments on Earth. The endolithic lifestyle is a remarkable adaptation to the exceptional environmental extremes of this area, which is considered the closest terrestrial example to conditions on Mars. Recent efforts have attempted to elucidate composition of these extremely adapted communities, but the functionality of these microbes have remained unexplored. We have tested for interactions between measured environmental characteristics, fungal community membership, and inferred functional classification of the fungi present and found altitude and sun exposure were primary factors. Sandstone rocks were collected in Victoria Land, Antarctica along an altitudinal gradient from 834 to 3100 m a.s.l.; differently sun-exposed rocks were selected to test the influence of this parameter on endolithic settlement. Metabarcoding targeting the fungal internal transcribed spacer region 1 (ITS1) was used to catalogue the species found in these communities. Functional profile of guilds found in the samples was associated to species using FUNGuild and variation in functional groups compared across sunlight exposure and altitude. Results revealed clear dominance of lichenized and stress-tolerant fungi in endolithic communities. The main variations in composition and abundance of functional groups among sites correlated to sun exposure, but not to altitude.

Arbuscular mycorrhizal fungal community composition is altered by long-term litter removal but not litter addition in a lowland tropical forest

Tropical forest productivity is sustained by the cycling of nutrients through decomposing organic matter. Arbuscular mycorrhizal (AM) fungi play a key role in the nutrition of tropical trees, yet there has been little experimental investigation into the role of AM fungi in nutrient cycling via decomposing organic material in tropical forests. We evaluated the responses of AM fungi in a long-term leaf litter addition and removal experiment in a tropical forest in Panama. We described AM fungal communities using 454-pyrosequencing, quantified the proportion of root length colonised by AM fungi using microscopy, and estimated AM fungal biomass using a lipid biomarker. AM fungal community composition was altered by litter removal but not litter addition. Root colonisation was substantially greater in the superficial organic layer compared with the mineral soil. Overall colonisation was lower in the litter removal treatment, which lacked an organic layer. There was no effect of litter manipulation on the concentration of the AM fungal lipid biomarker in the mineral soil. We hypothesise that reductions in organic matter brought about by litter removal may lead to AM fungi obtaining nutrients from recalcitrant organic or mineral sources in the soil, besides increasing fungal competition for progressively limited resources.

Diversity of Phototrophic Genes Suggests Multiple Bacteria May Be Able to Exploit Sunlight in Exposed Soils from the Sør Rondane Mountains, East Antarctica

Microbial life in exposed terrestrial surface layers in continental Antarctica is faced with extreme environmental conditions, including scarcity of organic matter. Bacteria in these exposed settings can therefore be expected to use alternative energy sources such as solar energy, abundant during the austral summer. Using Illumina MiSeq sequencing, we assessed the diversity and abundance of four conserved protein encoding genes involved in different key steps of light-harvesting pathways dependent on (bacterio)chlorophyll (pufM, bchL/chlL, and bchX genes) and rhodopsins (actinorhodopsin genes), in exposed soils from the Sør Rondane Mountains, East Antarctica. Analysis of pufM genes, encoding a subunit of the type 2 photochemical reaction center found in anoxygenic phototrophic bacteria, revealed a broad diversity, dominated by Roseobacter- and Loktanella-like sequences. The bchL and chlL, involved in (bacterio)chlorophyll synthesis, on the other hand, showed a high relative abundance of either cyanobacterial or green algal trebouxiophyceael chlL reads, depending on the sample, while most bchX sequences belonged mostly to previously unidentified phylotypes. Rhodopsin-containing phototrophic bacteria could not be detected in the samples. Our results, while suggesting that Cyanobacteria and green algae are the main phototrophic groups, show that light-harvesting bacteria are nevertheless very diverse in microbial communities in Antarctic soils.

Alpine bistort (Bistorta vivipara) in edge habitat associates with fewer but distinct ectomycorrhizal fungal species: a comparative study of three contrasting soil environments in Svalbard

Bistorta vivipara is a widespread arctic-alpine ectomycorrhizal (ECM) plant species. Recent findings suggest that fungal communities associated with B. vivipara roots appear random over short distances, but at larger scales, environmental filtering structure fungal communities. Habitats in highly stressful environments where specialist species with narrower niches may have an advantage represent unique opportunity to test the effect of environmental filtering. We utilised high-throughput amplicon sequencing to identify ECM communities associated with B. vivipara in Svalbard. We compared ECM communities in a core habitat where B. vivipara is frequent (Dryas-heath) with edge habitats representing extremes in terms of nutrient availability where B. vivipara is less frequent (bird-manured meadow and a nutrient-depleted mine tilling). Our analysis revealed that soil conditions in edge habitats favour less diverse but more distinct ECM fungal communities with functional traits adapted to local conditions. ECM richness was overall lower in both edge habitats, and the taxonomic compositions of ECM fungi were in line with our functional expectations. Stress-tolerant genera such as Laccaria and Hebeloma were abundant in nutrient-poor mine site whereas functional competitors genera such as Lactarius and Russula were dominant in the nutrient-rich bird-cliff site. Our results suggest that ECM communities in rare edge habitats are most likely not subsets of the larger pool of ECM fungi found in natural tundra, and they may represent a significant contribution to the overall diversity of ECM fungi in the Arctic.

Long-term increase in snow depth leads to compositional changes in arctic ectomycorrhizal fungal communities

Many arctic ecological processes are regulated by soil temperature that is tightly interconnected with snow cover distribution and persistence. Recently, various climate-induced changes have been observed in arctic tundra ecosystems, e.g. shrub expansion, resulting in reduction in albedo and greater C fixation in aboveground vegetation as well as increased rates of soil C mobilization by microbes. Importantly, the net effects of these shifts are unknown, in part because our understanding of belowground processes is limited. Here, we focus on the effects of increased snow depth, and as a consequence, increased winter soil temperature on ectomycorrhizal (ECM) fungal communities in dry and moist tundra. We analyzed deep DNA sequence data from soil samples taken at a long-term snow fence experiment in Northern Alaska. Our results indicate that, in contrast with previously observed responses of plants to increased snow depth at the same experimental site, the ECM fungal community of the dry tundra was more affected by deeper snow than the moist tundra community. In the dry tundra, both community richness and composition were significantly altered while in the moist tundra, only community composition changed significantly while richness did not. We observed a decrease in richness of Tomentella, Inocybe and other taxa adapted to scavenge the soil for labile N forms. On the other hand, richness of Cortinarius, and species with the ability to scavenge the soil for recalcitrant N forms, did not change. We further link ECM fungal traits with C soil pools. If future warmer atmospheric conditions lead to greater winter snow fall, changes in the ECM fungal community will likely influence C emissions and C fixation through altering N plant availability, fungal biomass and soil-plant C-N dynamics, ultimately determining important future interactions between the tundra biosphere and atmosphere.

The role of phosphorus, magnesium and potassium availability in soil fungal exploration of mineral nutrient sources in Norway spruce forests

We investigated fungal growth and community composition in buried meshbags, amended with apatite, biotite or hornblende, in Norway spruce (Picea abies) forests of varying nutrient status. Norway spruce needles and soil collected from forests overlying serpentinite had low levels of potassium and phosphorus, those from granite had low levels of magnesium, whereas those from amphibolite had comparably high levels of these nutrients. We assayed the fungal colonization of meshbags by measuring ergosterol content and fungal community with 454 sequencing of the internal transcribed spacer region. In addition, we measured fine root density. Fungal biomass was increased by apatite amendment across all plots and particularly on the K- and P-deficient serpentinite plots, whereas hornblende and biotite had no effect on fungal biomass on any plots. Fungal community (total fungal and ectomycorrhizal) composition was affected strongly by sampling location and soil depth, whereas mineral amendments had no effect on community composition. Fine root biomass was significantly correlated with fungal biomass. Ectomycorrhizal communities may respond to increased host-tree phosphorus demand by increased colonization of phosphorus-containing minerals, but this does not appear to translate to a shift in ectomycorrhizal community composition. This growth response to nutrient demand does not appear to exist for potassium or magnesium limitation.

Composition of fungal soil communities varies with plant abundance and geographic origin

Interactions of belowground fungal communities with exotic and native plant species may be important drivers of plant community structure in invaded grasslands. However, field surveys linking plant community structure with belowground fungal communities are missing. We investigated whether a selected number of abundant and relatively rare plants, either native or exotic, from an old-field site associate with different fungal communities. We also assessed whether these plants showed different symbiotic relationships with soil biota through their roots. We characterized the plant community and collected roots to investigate fungal communities using 454 pyrosequencing and assessed arbuscular mycorrhizal colonization and enemy-induced lesions. Differences in fungal communities were considered based on the assessment of α- and β diversity depending on plant 'abundance' and 'origin'. Plant abundance and origin determined the fungal community. Fungal richness was higher for native abundant as opposed to relatively rare native plant species. However, this was not observed for exotics of contrasting abundance. Regardless of their origin, β diversity was higher for rare than for abundant species. Abundant exotics in the community, which happen to be grasses, were the least mycorrhizal whereas rare natives were most susceptible to enemy attack. Our results suggest that compared with exotics, the relative abundance of remnant native plant species in our old-field site is still linked to the structure of belowground fungal communities. In contrast, exotic species may act as a disturbing agent contributing towards the homogenization of soil fungal communities, potentially changing feedback interactions.

Long-term warming alters richness and composition of taxonomic and functional groups of arctic fungi

Fungi, including symbionts, pathogens and decomposers, play crucial roles in community dynamics and nutrient cycling in terrestrial ecosystems. Despite their ecological importance, the response of most arctic fungi to climate warming is unknown, so are their potential roles in driving the observed and predicted changes in tundra communities. We carried out deep DNA sequencing of soil samples to study the long-term effects of experimental warming on fungal communities in dry heath and moist tussock tundra in Arctic Alaska. The data presented here indicate that fungal community composition responds strongly to warming in the moist tundra, but not in the dry tundra. While total fungal richness was not significantly affected by warming, there were clear correlations among operational taxonomic unit richness of various ecological and taxonomic groups and long-term warming. Richness of ectomycorrhizal, ericoid mycorrhizal and lichenized fungi generally decreased with warming, while richness of saprotrophic, plant and animal pathogenic, and root endophytic fungi tended to increase in the warmed plots. More importantly, various taxa within these functional guilds followed opposing trends that highlight the importance of species-specific responses to warming. We recommend that species-level ecological differences be taken into account in climate change and nutrient cycling studies that involve arctic fungi.

The Quantification of Representative Sequences pipeline for amplicon sequencing: case study on within-population ITS1 sequence variation in a microparasite infecting Daphnia

Next generation sequencing (NGS) platforms are replacing traditional molecular biology protocols like cloning and Sanger sequencing. However, accuracy of NGS platforms has rarely been measured when quantifying relative frequencies of genotypes or taxa within populations. Here we developed a new bioinformatic pipeline (QRS) that pools similar sequence variants and estimates their frequencies in NGS data sets from populations or communities. We tested whether the estimated frequency of representative sequences, generated by 454 amplicon sequencing, differs significantly from that obtained by Sanger sequencing of cloned PCR products. This was performed by analysing sequence variation of the highly variable first internal transcribed spacer (ITS1) of the ichthyosporean Caullerya mesnili, a microparasite of cladocerans of the genus Daphnia. This analysis also serves as a case example of the usage of this pipeline to study within-population variation. Additionally, a public Illumina data set was used to validate the pipeline on community-level data. Overall, there was a good correspondence in absolute frequencies of C. mesnili ITS1 sequences obtained from Sanger and 454 platforms. Furthermore, analyses of molecular variance (amova) revealed that population structure of C. mesnili differs across lakes and years independently of the sequencing platform. Our results support not only the usefulness of amplicon sequencing data for studies of within-population structure but also the successful application of the QRS pipeline on Illumina-generated data. The QRS pipeline is freely available together with its documentation under GNU Public Licence version 3 at http://code.google.com/p/quantification-representative-sequences.

Arctic fungal communities associated with roots of Bistorta vivipara do not respond to the same fine-scale edaphic gradients as the aboveground vegetation

Soil conditions and microclimate are important determinants of the fine-scale distribution of plant species in the Arctic, creating locally heterogeneous vegetation. We hypothesize that root-associated fungal (RAF) communities respond to the same fine-scale environmental gradients as the aboveground vegetation, creating a coherent pattern between aboveground vegetation and RAF. We explored how RAF communities of the ectomycorrhizal (ECM) plant Bistorta vivipara and aboveground vegetation structure of arctic plants were affected by biotic and abiotic variables at 0.3-3.0-m scales. RAF communities were determined using pyrosequencing. Composition and spatial structure of RAF and aboveground vegetation in relation to collected biotic and abiotic variables were analysed by ordination and semi-variance analyses. The vegetation was spatially structured along soil C and N gradients, whereas RAF lacked significant spatial structure. A weak relationship between RAF community composition and the cover of two ECM plants, B. vivipara and S. polaris, was found, and RAF richness increased with host root length and root weight. Results suggest that the fine-scale spatial structure of RAF communities of B. vivipara and the aboveground vegetation are driven by different factors. At fine spatial scales, neighbouring ECM plants may affect RAF community composition, whereas soil nutrients gradients structure the vegetation.

Host identity is a dominant driver of mycorrhizal fungal community composition during ecosystem development

Little is known about the response of arbuscular mycorrhizal fungal communities to ecosystem development. We use a long-term soil chronosequence that includes ecosystem progression and retrogression to quantify the importance of host plant identity as a factor driving fungal community composition during ecosystem development. We identified arbuscular mycorrhizal fungi and plant species from 50 individual roots from each of 10 sites spanning 5-120 000 yr of ecosystem age using terminal restriction fragment length polymorphism (T-RFLP), Sanger sequencing and pyrosequencing. Arbuscular mycorrhizal fungal communities were highly structured by ecosystem age. There was strong niche differentiation, with different groups of operational taxonomic units (OTUs) being characteristic of early succession, ecosystem progression and ecosystem retrogression. Fungal alpha diversity decreased with ecosystem age, whereas beta diversity was high at early stages and lower in subsequent stages. A total of 39% of the variance in fungal communities was explained by host plant and site age, 29% of which was attributed to host and the interaction between host and site (24% and 5%, respectively). The strong response of arbuscular mycorrhizal fungi to ecosystem development appears to be largely driven by plant host identity, supporting the concept that plant and fungal communities are tightly coupled rather than independently responding to habitat.

Summer temperature increase has distinct effects on the ectomycorrhizal fungal communities of moist tussock and dry tundra in Arctic Alaska

Arctic regions are experiencing the greatest rates of climate warming on the planet and marked changes have already been observed in terrestrial arctic ecosystems. While most studies have focused on the effects of warming on arctic vegetation and nutrient cycling, little is known about how belowground communities, such as fungi root-associated, respond to warming. Here, we investigate how long-term summer warming affects ectomycorrhizal (ECM) fungal communities. We used Ion Torrent sequencing of the rDNA internal transcribed spacer 2 (ITS2) region to compare ECM fungal communities in plots with and without long-term experimental warming in both dry and moist tussock tundra. Cortinarius was the most OTU-rich genus in the moist tundra, while the most diverse genus in the dry tundra was Tomentella. On the diversity level, in the moist tundra we found significant differences in community composition, and a sharp decrease in the richness of ECM fungi due to warming. On the functional level, our results indicate that warming induces shifts in the extramatrical properties of the communities, where the species with medium-distance exploration type seem to be favored with potential implications for the mobilization of different nutrient pools in the soil. In the dry tundra, neither community richness nor community composition was significantly altered by warming, similar to what had been observed in ECM host plants. There was, however, a marginally significant increase in OTUs identified as ECM fungi with the medium-distance exploration type in the warmed plots. Linking our findings of decreasing richness with previous results of increasing ECM fungal biomass suggests that certain ECM species are favored by warming and may become more abundant, while many other species may go locally extinct due to direct or indirect effects of warming. Such compositional shifts in the community might affect nutrient cycling and soil organic C storage.

Characterization of the bacterial communities of life stages of free living lone star ticks (Amblyomma americanum)

The lone star tick (Amblyomma americanum) is an abundant and aggressive biter of humans, domestic animals, and wildlife in the southeastern-central USA and an important vector of several known and suspected zoonotic bacterial pathogens. However, the biological drivers of bacterial community variation in this tick are still poorly defined. Knowing the community context in which tick-borne bacterial pathogens exist and evolve is required to fully understand the ecology and immunobiology of the ticks and to design effective public health and veterinary interventions. We performed a metagenomic survey of the bacterial communities of questing A. americanum and tested 131 individuals (66 nymphs, 24 males, and 41 females) from five sites in three states. Pyrosequencing was performed with barcoded eubacterial primers targeting variable 16S rRNA gene regions 5–3. The bacterial communities were dominated by Rickettsia (likely R. amblyommii) and an obligate Coxiella symbiont, together accounting for 6.7–100% of sequences per tick. DNAs from Midichloria, Borrelia, Wolbachia, Ehrlichia, Pseudomonas, or unidentified Bacillales, Enterobacteriaceae, or Rhizobiales groups were also detected frequently. Wolbachia and Midichloria significantly co-occurred in Georgia (p<0.00001), but not in other states. The significance of the Midichloria-Wolbachia co-occurrence is unknown. Among ticks collected in Georgia, nymphs differed from adults in both the composition (p = 0.002) and structure (p = 0.002) of their bacterial communities. Adults differed only in their community structure (p = 0.002) with males containing more Rickettsia and females containing more Coxiella. Comparisons among adult ticks collected in New York and North Carolina supported the findings from the Georgia collection despite differences in geography, collection date, and sample handling, implying that the differences detected are consistent attributes. The data also suggest that some members of the bacterial community change during the tick life cycle and that some sex-specific attributes may be detectable in nymphs.

The contribution of DNA metabarcoding to fungal conservation: diversity assessment, habitat partitioning and mapping red-listed fungi in protected coastal Salix repens communities in the Netherlands

Western European coastal sand dunes are highly important for nature conservation. Communities of the creeping willow (Salix repens) represent one of the most characteristic and diverse vegetation types in the dunes. We report here the results of the first kingdom-wide fungal diversity assessment in S. repens coastal dune vegetation. We carried out massively parallel pyrosequencing of ITS rDNA from soil samples taken at ten sites in an extended area of joined nature reserves located along the North Sea coast of the Netherlands, representing habitats with varying soil pH and moisture levels. Fungal communities in Salix repens beds are highly diverse and we detected 1211 non-singleton fungal 97% sequence similarity OTUs after analyzing 688,434 ITS2 rDNA sequences. Our comparison along a north-south transect indicated strong correlation between soil pH and fungal community composition. The total fungal richness and the number OTUs of most fungal taxonomic groups negatively correlated with higher soil pH, with some exceptions. With regard to ecological groups, dark-septate endophytic fungi were more diverse in acidic soils, ectomycorrhizal fungi were represented by more OTUs in calcareous sites, while detected arbuscular mycorrhizal genera fungi showed opposing trends regarding pH. Furthermore, we detected numerous red listed species in our samples often from previously unknown locations, indicating that some of the fungal species currently considered rare may be more abundant in Dutch S. repens communities than previously thought.

Influence of artifact removal on rare species recovery in natural complex communities using high-throughput sequencing

Large-scale high-throughput sequencing techniques are rapidly becoming popular methods to profile complex communities and have generated deep insights into community biodiversity. However, several technical problems, especially sequencing artifacts such as nucleotide calling errors, could artificially inflate biodiversity estimates. Sequence filtering for artifact removal is a conventional method for deleting error-prone sequences from high-throughput sequencing data. As rare species represented by low-abundance sequences in datasets may be sensitive to artifact removal process, the influence of artifact removal on rare species recovery has not been well evaluated in natural complex communities. Here we employed both internal (reliable operational taxonomic units selected from communities themselves) and external (indicator species spiked into communities) references to evaluate the influence of artifact removal on rare species recovery using 454 pyrosequencing of complex plankton communities collected from both freshwater and marine habitats. Multiple analyses revealed three clear patterns: 1) rare species were eliminated during sequence filtering process at all tested filtering stringencies, 2) more rare taxa were eliminated as filtering stringencies increased, and 3) elimination of rare species intensified as biomass of a species in a community was reduced. Our results suggest that cautions be applied when processing high-throughput sequencing data, especially for rare taxa detection for conservation of species at risk and for rapid response programs targeting non-indigenous species. Establishment of both internal and external references proposed here provides a practical strategy to evaluate artifact removal process.

Variations in 16S rRNA-based microbiome profiling between pyrosequencing runs and between pyrosequencing facilities

Pyrosequencing of 16S rRNA gene amplicons on the 454 FLX Titanium platform has been widely used to analyze microbiomes in various environments. However, different results may stem from variations among sequencing runs or among sequencing facilities. This study aimed to evaluate these variations between different pyrosequencing runs by sequencing 16S rRNA gene amplicon libraries generated from three sets of rumen samples twice each on the 454 FLX Titanium system at two independent sequencing facilities. Similar relative abundances were found for predominant taxa represented by large numbers of sequence reads but not for minor taxa represented by small numbers of sequence reads. The two sequencing facilities revealed different bacterial profiles with respect to both predominant taxa and minor taxa, including the most predominant genus Prevotella, the family Lachnospiraceae, and the phylum Proteobacteria. Differences in primers used to generate amplicon libraries may be a major source of variations in microbiome profiling. Because different primers and regions of 16S rRNA genes are often used by different researchers, significant variations likely exist among studies. Quantitative interpretation for relative abundance of taxa, especially minor taxa, from prevalence of sequence reads and comparisons of results from different studies should be done with caution.

Different bacterial communities in ectomycorrhizae and surrounding soil

Several eukaryotic symbioses have shown to host a rich diversity of prokaryotes that interact with their hosts. Here, we study bacterial communities associated with ectomycorrhizal root systems of Bistorta vivipara compared to bacterial communities in bulk soil using pyrosequencing of 16S rRNA amplicons. A high richness of Operational Taxonomic Units (OTUs) was found in plant roots (3,571 OTUs) and surrounding soil (3,476 OTUs). The community composition differed markedly between these two environments. Actinobacteria, Armatimonadetes, Chloroflexi and OTUs unclassified at phylum level were significantly more abundant in plant roots than in soil. A large proportion of the OTUs, especially those in plant roots, presented low similarity to Sanger 16S rRNA reference sequences, suggesting novel bacterial diversity in ectomycorrhizae. Furthermore, the bacterial communities of the plant roots were spatially structured up to a distance of 60 cm, which may be explained by bacteria using fungal hyphae as a transport vector. The analyzed ectomycorrhizae presents a distinct microbiome, which likely influence the functioning of the plant-fungus symbiosis.

Improving accuracy of DNA diet estimates using food tissue control materials and an evaluation of proxies for digestion bias

Ecologists are increasingly interested in quantifying consumer diets based on food DNA in dietary samples and high-throughput sequencing of marker genes. It is tempting to assume that food DNA sequence proportions recovered from diet samples are representative of consumer's diet proportions, despite the fact that captive feeding studies do not support that assumption. Here, we examine the idea of sequencing control materials of known composition along with dietary samples in order to correct for technical biases introduced during amplicon sequencing and biological biases such as variable gene copy number. Using the Ion Torrent PGM(©) , we sequenced prey DNA amplified from scats of captive harbour seals (Phoca vitulina) fed a constant diet including three fish species in known proportions. Alongside, we sequenced a prey tissue mix matching the seals' diet to generate tissue correction factors (TCFs). TCFs improved the diet estimates (based on sequence proportions) for all species and reduced the average estimate error from 28 ± 15% (uncorrected) to 14 ± 9% (TCF-corrected). The experimental design also allowed us to infer the magnitude of prey-specific digestion biases and calculate digestion correction factors (DCFs). The DCFs were compared with possible proxies for differential digestion (e.g. fish protein%, fish lipid%) revealing a strong relationship between the DCFs and percent lipid of the fish prey, suggesting prey-specific corrections based on lipid content would produce accurate diet estimates in this study system. These findings demonstrate the value of parallel sequencing of food tissue mixtures in diet studies and offer new directions for future research in quantitative DNA diet analysis.

454 pyrosequencing to describe microbial eukaryotic community composition, diversity and relative abundance: a test for marine haptophytes

Next generation sequencing of ribosomal DNA is increasingly used to assess the diversity and structure of microbial communities. Here we test the ability of 454 pyrosequencing to detect the number of species present, and assess the relative abundance in terms of cell numbers and biomass of protists in the phylum Haptophyta. We used a mock community consisting of equal number of cells of 11 haptophyte species and compared targeting DNA and RNA/cDNA, and two different V4 SSU rDNA haptophyte-biased primer pairs. Further, we tested four different bioinformatic filtering methods to reduce errors in the resulting sequence dataset. With sequencing depth of 11000–20000 reads and targeting cDNA with Haptophyta specific primers Hap454 we detected all 11 species. A rarefaction analysis of expected number of species recovered as a function of sampling depth suggested that minimum 1400 reads were required here to recover all species in the mock community. Relative read abundance did not correlate to relative cell numbers. Although the species represented with the largest biomass was also proportionally most abundant among the reads, there was generally a weak correlation between proportional read abundance and proportional biomass of the different species, both with DNA and cDNA as template. The 454 sequencing generated considerable spurious diversity, and more with cDNA than DNA as template. With initial filtering based only on match with barcode and primer we observed 100-fold more operational taxonomic units (OTUs) at 99% similarity than the number of species present in the mock community. Filtering based on quality scores, or denoising with PyroNoise resulted in ten times more OTU_99% than the number of species. Denoising with AmpliconNoise reduced the number of OTU_99% to match the number of species present in the mock community. Based on our analyses, we propose a strategy to more accurately depict haptophyte diversity using 454 pyrosequencing.

Substantial compositional turnover of fungal communities in an alpine ridge-to-snowbed gradient

The main gradient in vascular plant, bryophyte and lichen species composition in alpine areas, structured by the topographic gradient from wind-exposed ridges to snowbeds, has been extensively studied. Tolerance to environmental stress, resulting from wind abrasion and desiccation towards windswept ridges or reduced growing season due to prolonged snow cover towards snowbeds, is an important ecological mechanism in this gradient. The extent to which belowground fungal communities are structured by the same topographic gradient and the eventual mechanisms involved are less well known. In this study, we analysed variation in fungal diversity and community composition associated with roots of the ectomycorrhizal plant Bistorta vivipara along the ridge-to-snowbed gradient. We collected root samples from fifty B. vivipara plants in ten plots in an alpine area in central Norway. The fungal communities were analysed using 454 pyrosequencing analyses of tag-encoded ITS1 amplicons. A distinct gradient in the fungal community composition was found that coincided with variation from ridge to snowbeds. This gradient was paralleled by change in soil content of carbon, nitrogen and phosphorus. A large proportion (66%) of the detected 801 nonsingleton operational taxonomic units (OTUs) were ascomycetes, while basidiomycetes dominated quantitatively (i.e. with respect to number of reads). Numerous fungal OTUs, many with taxonomic affinity to Sebacinales, Cortinarius and Meliniomyces, showed distinct affinities either to ridge or to snowbed plots, indicating habitat specialization. The compositional turnover of fungal communities along the gradient was not paralleled by a gradient in species richness.

Fungal community analysis by high-throughput sequencing of amplified markers--a user's guide

Novel high-throughput sequencing methods outperform earlier approaches in terms of resolution and magnitude. They enable identification and relative quantification of community members and offer new insights into fungal community ecology. These methods are currently taking over as the primary tool to assess fungal communities of plant-associated endophytes, pathogens, and mycorrhizal symbionts, as well as free-living saprotrophs. Taking advantage of the collective experience of six research groups, we here review the different stages involved in fungal community analysis, from field sampling via laboratory procedures to bioinformatics and data interpretation. We discuss potential pitfalls, alternatives, and solutions. Highlighted topics are challenges involved in: obtaining representative DNA/RNA samples and replicates that encompass the targeted variation in community composition, selection of marker regions and primers, options for amplification and multiplexing, handling of sequencing errors, and taxonomic identification. Without awareness of methodological biases, limitations of markers, and bioinformatics challenges, large-scale sequencing projects risk yielding artificial results and misleading conclusions.

Random sampling process leads to overestimation of β-diversity of microbial communities

The site-to-site variability in species composition, known as β-diversity, is crucial to understanding spatiotemporal patterns of species diversity and the mechanisms controlling community composition and structure. However, quantifying β-diversity in microbial ecology using sequencing-based technologies is a great challenge because of a high number of sequencing errors, bias, and poor reproducibility and quantification. Herein, based on general sampling theory, a mathematical framework is first developed for simulating the effects of random sampling processes on quantifying β-diversity when the community size is known or unknown. Also, using an analogous ball example under Poisson sampling with limited sampling efforts, the developed mathematical framework can exactly predict the low reproducibility among technically replicate samples from the same community of a certain species abundance distribution, which provides explicit evidences of random sampling processes as the main factor causing high percentages of technical variations. In addition, the predicted values under Poisson random sampling were highly consistent with the observed low percentages of operational taxonomic unit (OTU) overlap (<30% and <20% for two and three tags, respectively, based on both Jaccard and Bray-Curtis dissimilarity indexes), further supporting the hypothesis that the poor reproducibility among technical replicates is due to the artifacts associated with random sampling processes. Finally, a mathematical framework was developed for predicting sampling efforts to achieve a desired overlap among replicate samples. Our modeling simulations predict that several orders of magnitude more sequencing efforts are needed to achieve desired high technical reproducibility. These results suggest that great caution needs to be taken in quantifying and interpreting β-diversity for microbial community analysis using next-generation sequencing technologies.

Due to the vast diversity and uncultivated status of the majority of microorganisms, microbial detection, characterization, and quantitation are of great challenge. Although large-scale metagenome sequencing technology such as PCR-based amplicon sequencing has revolutionized the studies of microbial communities, it suffers from several inherent drawbacks, such as a high number of sequencing errors, biases, poor quantitation, and very high percentages of technical variations, which could greatly overestimate microbial biodiversity. Based on general sampling theory, this study provided the first explicit evidence to demonstrate the importance of random sampling processes in estimating microbial β-diversity, which has not been adequately recognized and addressed in microbial ecology. Since most ecological studies are involved in random sampling, the conclusions learned from this study should also be applicable to other ecological studies in general. In summary, the results presented in this study should have important implications for examining microbial biodiversity to address both basic theoretical and applied management questions.

Bacterial population in intestines of the black tiger shrimp (Penaeus monodon) under different growth stages

Intestinal bacterial communities in aquaculture have been drawn to attention due to potential benefit to their hosts. To identify core intestinal bacteria in the black tiger shrimp (Penaeus monodon), bacterial populations of disease-free shrimp were characterized from intestines of four developmental stages (15-day-old post larvae (PL15), 1- (J1), 2- (J2), and 3-month-old (J3) juveniles) using pyrosequencing, real-time PCR and denaturing gradient gel electrophoresis (DGGE) approaches. A total of 25,121 pyrosequencing reads (reading length = 442±24 bases) were obtained, which were categorized by barcode for PL15 (7,045 sequences), J1 (3,055 sequences), J2 (13,130 sequences) and J3 (1,890 sequences). Bacteria in the phyla Bacteroides, Firmicutes and Proteobacteria were found in intestines at all four growth stages. There were 88, 14, 27, and 20 bacterial genera associated with the intestinal tract of PL15, J1, J2 and J3, respectively. Pyrosequencing analysis revealed that Proteobacteria (class Gammaproteobacteria) was a dominant bacteria group with a relative abundance of 89% for PL15 and 99% for J1, J2 and J3. Real-time PCR assay also confirmed that Gammaproteobacteria had the highest relative abundance in intestines from all growth stages. Intestinal bacterial communities from the three juvenile stages were more similar to each other than that of the PL shrimp based on PCA analyses of pyrosequencing results and their DGGE profiles. This study provides descriptive bacterial communities associated to the black tiger shrimp intestines during these growth development stages in rearing facilities.

Short-read assembly of full-length 16S amplicons reveals bacterial diversity in subsurface sediments

In microbial ecology, a fundamental question relates to how community diversity and composition change in response to perturbation. Most studies have had limited ability to deeply sample community structure (e.g. Sanger-sequenced 16S rRNA libraries), or have had limited taxonomic resolution (e.g. studies based on 16S rRNA hypervariable region sequencing). Here, we combine the higher taxonomic resolution of near-full-length 16S rRNA gene amplicons with the economics and sensitivity of short-read sequencing to assay the abundance and identity of organisms that represent as little as 0.01% of sediment bacterial communities. We used a new version of EMIRGE optimized for large data size to reconstruct near-full-length 16S rRNA genes from amplicons sheared and sequenced with Illumina technology. The approach allowed us to differentiate the community composition among samples acquired before perturbation, after acetate amendment shifted the predominant metabolism to iron reduction, and once sulfate reduction began. Results were highly reproducible across technical replicates, and identified specific taxa that responded to the perturbation. All samples contain very high alpha diversity and abundant organisms from phyla without cultivated representatives. Surprisingly, at the time points measured, there was no strong loss of evenness, despite the selective pressure of acetate amendment and change in the terminal electron accepting process. However, community membership was altered significantly. The method allows for sensitive, accurate profiling of the “long tail” of low abundance organisms that exist in many microbial communities, and can resolve population dynamics in response to environmental change.

ITS1 versus ITS2 as DNA metabarcodes for fungi

The nuclear ribosomal Internal Transcribed Spacer ITS region is widely used as a DNA metabarcoding marker to characterize the diversity and composition of fungal communities. In amplicon pyrosequencing studies of fungal diversity, one of the spacers ITS1 or ITS2 of the ITS region is normally used. In this methodological study we evaluate the usability of ITS1 vs. ITS2 as a DNA metabarcoding marker for fungi. We analyse three data sets: two comprising ITS1 and ITS2 sequences of known taxonomic affiliations and a third comprising ITS1 and ITS2 environmental amplicon pyrosequencing data. Clustering analyses of sequences with known taxonomy using the bioinformatics pipeline ClustEx revealed that a 97% similarity cut-off represent a reasonable threshold for estimating the number of known species in the data sets for both ITS1 and ITS2. However, no single threshold value worked well for all fungi at the same time within the curated UNITE database, and we found that the Operational Taxonomic Unit (OTU) concept is not easily translated into the level of species because many species are distributed over several clusters. Clustering analyses of the 134 692 ITS1 and ITS2 pyrosequences using a 97% similarity cut-off revealed a high similarity between the two data sets when it comes to taxonomic coverage. Although some groups are under- or unrepresented in the two data sets due to, e.g. primer mismatches, our results indicate that ITS1 and ITS2 to a large extent yield similar results when used as DNA metabarcodes for fungi.

Seasonal trends in the biomass and structure of bryophyte-associated fungal communities explored by 454 pyrosequencing

Bryophytes are a dominant vegetation component of the boreal forest, but little is known about their associated fungal communities, including seasonal variation within them. Seasonal variation in the fungal biomass and composition of fungal communities associated with three widespread boreal bryophytes was investigated using HPLC assays of ergosterol and amplicon pyrosequencing of the internal transcribed spacer 2 (ITS2) region of rDNA. The bryophyte phyllosphere community was dominated by Ascomycota. Fungal biomass did not decline appreciably in winter (P=0.272). Significant host-specific patterns in seasonal variation of biomass were detected (P=0.003). Although seasonal effects were not the primary factors structuring community composition, collection date significantly explained (P=0.001) variation not attributed to locality, host, and tissue. Community homogenization and a reduction in turnover occurred with the onset of frost events and subzero air and soil temperatures. Fluctuations in the relative abundance of particular fungal groups seem to reflect the nature of their association with mosses, although conclusions are drawn with caution because of potential methodological bias. The moss-associated fungal community is dynamic, exhibiting seasonal turnover in composition and relative abundance of different fungal groups, and significant fungal biomass is present year-round, suggesting a winter-active fungal community.

Bacterial associates of two Caribbean coral species reveal species-specific distribution and geographic variability

Scleractinian corals harbor microorganisms that form dynamic associations with the coral host and exhibit substantial genetic and ecological diversity. Microbial associates may provide defense against pathogens and serve as bioindicators of changing environmental conditions. Here we describe the bacterial assemblages associated with two of the most common and phylogenetically divergent reef-building corals in the Caribbean, Montastraea faveolata and Porites astreoides. Contrasting life history strategies and disease susceptibilities indicate potential differences in their microbiota and immune function that may in part drive changes in the composition of coral reef communities. The ribotype structure and diversity of coral-associated bacteria within the surface mucosal layer (SML) of healthy corals were assessed using denaturing gradient gel electrophoresis (DGGE) fingerprinting and 454 bar-coded pyrosequencing. Corals were sampled at disparate Caribbean locations representing various levels of anthropogenic impact. We demonstrate here that M. faveolata and P. astreoides harbor distinct, host-specific bacteria but that specificity varies by species and site. P. astreoides generally hosts a bacterial assemblage of low diversity that is largely dominated by one bacterial genus, Endozoicomonas, within the order Oceanospirillales. The bacterial assemblages associated with M. faveolata are significantly more diverse and exhibit higher specificity at the family level than P. astreoides assemblages. Both corals have more bacterial diversity and higher abundances of disease-related bacteria at sites closer to the mainland than at those furthest away. The most diverse bacterial taxa and highest relative abundance of disease-associated bacteria were seen for corals near St. Thomas, U.S. Virgin Islands (USVI) (2.5 km from shore), and the least diverse taxa and lowest relative abundance were seen for corals near our most pristine site in Belize (20 km from shore). We conclude that the two coral species studied harbor distinct bacterial assemblages within the SML, but the degree to which each species maintains specific microbial associations varies both within each site and across large spatial scales. The taxonomic scale (i.e., phylum versus genus) at which scientists examine coral-microbe associations, in addition to host-elicited factors and environmental fluctuations, must be considered carefully in future studies of the coral holobiont.