Sex-specific asymmetry within the cloacal microbiota of the striped plateau lizard, <Emphasis Type="Italic">Sceloporus virgatus</Emphasis>

The structure and diversity of microbial communities in wild vertebrate populations remain poorly understood, but are expected to have important consequences for individual survival and reproductive success. For instance, recent work has demonstrated that cloacal microbe assemblages of wild birds are related to the phenotypic quality of the host. To contribute to this field of study, we examined the composition and diversity of the cloacal microbiota of free-ranging striped plateau lizards, Sceloporus virgatus, using 16s rRNA-based culture independent techniques. Our dataset, generated from cloacal swabs of six males and six females, and based on twenty five 16s rRNA clones from each sample, revealed (i) low overall microbial diversity, (ii) a striking sex asymmetry in microbial community composition with males displaying cloacal microbiota more typical of gastrointestinal residents found in other organisms, while females display only gammaproteobacterial phylotypes, (iii) a significant sex difference in microbial community structure, with females having significantly lower microbial diversity and richness than do males, and (iv) that the diversity of the female microbial community is negatively correlated to her ectoparasitic mite load. It is not yet clear if the female-specific paucity of cloacal microbial diversity is due to host function or microbe-microbe interactions, or whether the relationship to female mite load is causal, however these findings are expected to have relevance to the species’ life history and ecology. Although the diversity of microbiota from humans, mice, birds, zebrafish, and invertebrates is widely investigated, this is one of only a few reports in the literature describing the cloacal microbiota of a wild vertebrate, and is perhaps the first report for wild reptiles that utilizes culture-independent techniques.     

Microbial community diversity and composition varies with habitat characteristics and biofilm function in macrophyte‐rich streams

Biofilms in streams play an integral role in ecosystem processes and function yet few studies have investigated the broad diversity of these complex prokaryotic and eukaryotic microbial communities. Physical habitat characteristics can affect the composition and abundance of microorganisms in these biofilms by creating microhabitats. Here we describe the prokaryotic and eukaryotic microbial diversity of biofilms in sand and macrophyte habitats (i.e. epipsammon and epiphyton, respectively) in five macrophyte‐rich streams in Jutland, Denmark. The macrophyte species varied in growth morphology, C:N stoichiometry, and preferred stream habitat, providing a range in environmental conditions for the epiphyton. Among all habitats and streams, the prokaryotic communities were dominated by common phyla, including Alphaproteobacteria, Bacteriodetes, and Gammaproteobacteria, while the eukaryotic communities were dominated by Stramenopiles (i.e. diatoms). For both the prokaryotes and eukaryotes, the epipsammon were consistently the most diverse communities and the epiphytic communities were generally similar among the four macrophyte species. However, the communities on the least complex macrophyte, Sparganium emersum, had the lowest richness and evenness and fewest unique OTUs, whereas the macrophyte with the most morphological complexity, Callitriche spp., had the highest number of unique OTUs. In general, the microbial taxa were ubiquitously distributed across the relatively homogeneous Danish landscape as determined by measuring the similarity among communities (i.e. Sørensen similarity index). Furthermore, we found significant correlations between microbial diversity (i.e. Chao1 rarefied richness and Pielou's evenness) and biofilm structure and function (i.e. C:N ratio and ammonium uptake efficiency, respectively); communities with higher richness and evenness had higher C:N ratios and lower uptake efficiency. In addition to describing the prokaryotic and eukaryotic community composition in stream biofilms, our study indicates that 1) physical habitat characteristics influence microbial diversity and 2) the variation in microbial diversity may dictate the structural and functional characteristics of stream biofilm communities.     

Variation in gut microbial communities and its association with pathogen infection in wild bumble bees (Bombus)

Bacterial gut symbiont communities are critical for the health of many insect species. However, little is known about how microbial communities vary among host species or how they respond to anthropogenic disturbances. Bacterial communities that differ in richness or composition may vary in their ability to provide nutrients or defenses. We used deep sequencing to investigate gut microbiota of three species in the genus Bombus (bumble bees). Bombus are among the most economically and ecologically important non-managed pollinators. Some species have experienced dramatic declines, probably due to pathogens and land-use change. We examined variation within and across bee species and between semi-natural and conventional agricultural habitats. We categorized as ‘core bacteria'' any operational taxonomic units (OTUs) with closest hits to sequences previously found exclusively or primarily in the guts of honey bees and bumble bees (genera Apis and Bombus). Microbial community composition differed among bee species. Richness, defined as number of bacterial OTUs, was highest for B. bimaculatus and B. impatiens. For B. bimaculatus, this was due to high richness of non-core bacteria. We found little effect of habitat on microbial communities. Richness of non-core bacteria was negatively associated with bacterial abundance in individual bees, possibly due to deeper sampling of non-core bacteria in bees with low populations of core bacteria. Infection by the gut parasite Crithidia was negatively associated with abundance of the core bacterium Gilliamella and positively associated with richness of non-core bacteria. Our results indicate that Bombus species have distinctive gut communities, and community-level variation is associated with pathogen infection.     

Midgut fungal and bacterial microbiota of Aedes triseriatus and Aedes japonicus shift in response to La Crosse virus infection

Understanding how midgut microbial communities of field‐collected mosquitoes interact with pathogens is critical for controlling vector infection and disease. We used 16S rRNA and internal transcribed spacer sequencing to characterize the midgut bacterial and fungal communities of adult females of Aedes triseriatus and Aedes japonicus collected as pupae in tree holes, plastic bins and waste tires and their response to La Crosse virus (LACV) infection. For both mosquito species and across all habitat and virus treatments, a total of 62 bacterial operational taxonomic units (OTUs) from six phyla and 21 fungal OTUs from two phyla were identified. The majority of bacterial (92%) and fungal (71%) OTUs were shared between the mosquito species; however, several OTUs were unique to each species. Bacterial and fungal communities of individuals that took either infectious or noninfectious bloodmeals were less diverse and more homogeneous compared to those of newly emerged adults. Interestingly, LACV‐infected A. triseriatus and A. japonicus had higher bacterial richness and lower fungal richness compared to individuals that took a noninfectious bloodmeal, suggesting that viral infection was associated with an increase in bacterial OTUs and a decrease in fungal OTUs. For both mosquito species, several OTUs were identified that had both high fidelity and specificity to mosquito midguts that were infected with LACV. Overall, these findings demonstrate that bacterial and fungal communities that reside in mosquito midguts respond to host diet and viral infection and could play a role in modulating vector susceptibility to LACV.     

Metagenomic and Small-Subunit rRNA Analyses Reveal the Genetic Diversity of Bacteria, Archaea, Fungi, and Viruses in Soil

Recent studies have highlighted the surprising richness of soil bacterial communities; however, bacteria are not the only microorganisms found in soil. To our knowledge, no study has compared the diversities of the four major microbial taxa, i.e., bacteria, archaea, fungi, and viruses, from an individual soil sample. We used metagenomic and small-subunit RNA-based sequence analysis techniques to compare the estimated richness and evenness of these groups in prairie, desert, and rainforest soils. By grouping sequences at the 97% sequence similarity level (an operational taxonomic unit [OTU]), we found that the archaeal and fungal communities were consistently less even than the bacterial communities. Although total richness levels are difficult to estimate with a high degree of certainty, the estimated number of unique archaeal or fungal OTUs appears to rival or exceed the number of unique bacterial OTUs in each of the collected soils. In this first study to comprehensively survey viral communities using a metagenomic approach, we found that soil viruses are taxonomically diverse and distinct from the communities of viruses found in other environments that have been surveyed using a similar approach. Within each of the four microbial groups, we observed minimal taxonomic overlap between sites, suggesting that soil archaea, bacteria, fungi, and viruses are globally as well as locally diverse.     

Metagenomic and small-subunit rRNA analyses reveal the genetic diversity of bacteria, archaea, fungi, and viruses in soil

Recent studies have highlighted the surprising richness of soil bacterial communities; however, bacteria are not the only microorganisms found in soil. To our knowledge, no study has compared the diversities of the four major microbial taxa, i.e., bacteria, archaea, fungi, and viruses, from an individual soil sample. We used metagenomic and small-subunit RNA-based sequence analysis techniques to compare the estimated richness and evenness of these groups in prairie, desert, and rainforest soils. By grouping sequences at the 97% sequence similarity level (an operational taxonomic unit [OTU]), we found that the archaeal and fungal communities were consistently less even than the bacterial communities. Although total richness levels are difficult to estimate with a high degree of certainty, the estimated number of unique archaeal or fungal OTUs appears to rival or exceed the number of unique bacterial OTUs in each of the collected soils. In this first study to comprehensively survey viral communities using a metagenomic approach, we found that soil viruses are taxonomically diverse and distinct from the communities of viruses found in other environments that have been surveyed using a similar approach. Within each of the four microbial groups, we observed minimal taxonomic overlap between sites, suggesting that soil archaea, bacteria, fungi, and viruses are globally as well as locally diverse.     

Toward a census of bacteria in soil

For more than a century, microbiologists have sought to determine the species richness of bacteria in soil, but the extreme complexity and unknown structure of soil microbial communities have obscured the answer. We developed a statistical model that makes the problem of estimating richness statistically accessible by evaluating the characteristics of samples drawn from simulated communities with parametric community distributions. We identified simulated communities with rank-abundance distributions that followed a truncated lognormal distribution whose samples resembled the structure of 16S rRNA gene sequence collections made using Alaskan and Minnesotan soils. The simulated communities constructed based on the distribution of 16S rRNA gene sequences sampled from the Alaskan and Minnesotan soils had a richness of 5,000 and 2,000 operational taxonomic units (OTUs), respectively, where an OTU represents a collection of sequences not more than 3% distant from each other. To sample each of these OTUs in the Alaskan 16S rRNA gene library at least twice, 480,000 sequences would be required; however, to estimate the richness of the simulated communities using nonparametric richness estimators would require only 18,000 sequences. Quantifying the richness of complex environments such as soil is an important step in building an ecological framework. We have shown that generating sufficient sequence data to do so requires less sequencing effort than completely sequencing a bacterial genome.     

Pyrosequencing reveals highly diverse and species-specific microbial communities in sponges from the Red Sea

Marine sponges are associated with a remarkable array of microorganisms. Using a tag pyrosequencing technology, this study was the first to investigate in depth the microbial communities associated with three Red Sea sponges, Hyrtios erectus, Stylissa carteri and Xestospongia testudinaria. We revealed highly diverse sponge-associated bacterial communities with up to 1000 microbial operational taxonomic units (OTUs) and richness estimates of up to 2000 species. Altogether, 26 bacterial phyla were detected from the Red Sea sponges, 11 of which were absent from the surrounding sea water and 4 were recorded in sponges for the first time. Up to 100 OTUs with richness estimates of up to 300 archaeal species were revealed from a single sponge species. This is by far the highest archaeal diversity ever recorded for sponges. A non-negligible proportion of unclassified reads was observed in sponges. Our results demonstrated that the sponge-associated microbial communities remained highly consistent in the same sponge species from different locations, although they varied at different degrees among different sponge species. A significant proportion of the tag sequences from the sponges could be assigned to one of the sponge-specific clusters previously defined. In addition, the sponge-associated microbial communities were consistently divergent from those present in the surrounding sea water. Our results suggest that the Red Sea sponges possess highly sponge-specific or even sponge-species-specific microbial communities that are resistant to environmental disturbance, and much of their microbial diversity remains to be explored.     

Ironing out the wrinkles in the rare biosphere through improved OTU clustering

Deep sequencing of PCR amplicon libraries facilitates the detection of low‐abundance populations in environmental DNA surveys of complex microbial communities. At the same time, deep sequencing can lead to overestimates of microbial diversity through the generation of low‐frequency, error‐prone reads. Even with sequencing error rates below 0.005 per nucleotide position, the common method of generating operational taxonomic units (OTUs) by multiple sequence alignment and complete‐linkage clustering significantly increases the number of predicted OTUs and inflates richness estimates. We show that a 2% single‐linkage preclustering methodology followed by an average‐linkage clustering based on pairwise alignments more accurately predicts expected OTUs in both single and pooled template preparations of known taxonomic composition. This new clustering method can reduce the OTU richness in environmental samples by as much as 30–60% but does not reduce the fraction of OTUs in long‐tailed rank abundance curves that defines the rare biosphere.     

Impact of Virioplankton on Archaeal and Bacterial Community Richness as Assessed in Seawater Batch Cultures

During cruises in the tropical Atlantic Ocean (January to February 2000) and the southern North Sea (December 2000), experiments were conducted to monitor the impact of virioplankton on archaeal and bacterial community richness. Prokaryotic cells equivalent to 10 to 100% of the in situ abundance were inoculated into virus-free seawater, and viruses equivalent to 35 to 360% of the in situ abundance were added. Batch cultures with microwave-inactivated viruses and without viruses served as controls. The apparent richness of archaeal and bacterial communities was determined by terminal restriction fragment length polymorphism (T-RFLP) analysis of PCR-amplified 16S rRNA gene fragments. Although the estimated richness of the prokaryotic communities generally was greatly reduced within the first 24 h of incubation due to confinement, the effects of virus amendment were detected at the level of individual operational taxonomic units (OTUs) in the T-RFLP patterns of both groups, Archaea and Bacteria. One group of OTUs was detected in the control samples but was absent from the virus-treated samples. This negative response of OTUs to virus amendment probably was caused by viral lysis. Additionally, we found OTUs not responding to the amendments, and several OTUs exhibited variable responses to the addition of inactive or active viruses. Therefore, we conclude that individual members of pelagic archaeal and bacterial communities can be differently affected by the presence of virioplankton.     

Investigating Microbial Eukaryotic Diversity from a Global Census: Insights from a Comparison of Pyrotag and Full-Length Sequences of 18S rRNA Genes

Next-generation DNA sequencing (NGS) approaches are rapidly surpassing Sanger sequencing for characterizing the diversity of natural microbial communities. Despite this rapid transition, few comparisons exist between Sanger sequences and the generally much shorter reads of NGS. Operational taxonomic units (OTUs) derived from full-length (Sanger sequencing) and pyrotag (454 sequencing of the V9 hypervariable region) sequences of 18S rRNA genes from 10 global samples were analyzed in order to compare the resulting protistan community structures and species richness. Pyrotag OTUs called at 98% sequence similarity yielded numbers of OTUs that were similar overall to those for full-length sequences when the latter were called at 97% similarity. Singleton OTUs strongly influenced estimates of species richness but not the higher-level taxonomic composition of the community. The pyrotag and full-length sequence data sets had slightly different taxonomic compositions of rhizarians, stramenopiles, cryptophytes, and haptophytes, but the two data sets had similarly high compositions of alveolates. Pyrotag-based OTUs were often derived from sequences that mapped to multiple full-length OTUs at 100% similarity. Thus, pyrotags sequenced from a single hypervariable region might not be appropriate for establishing protistan species-level OTUs. However, nonmetric multidimensional scaling plots constructed with the two data sets yielded similar clusters, indicating that beta diversity analysis results were similar for the Sanger and NGS sequences. Short pyrotag sequences can provide holistic assessments of protistan communities, although care must be taken in interpreting the results. The longer reads (>500 bp) that are now becoming available through NGS should provide powerful tools for assessing the diversity of microbial eukaryotic assemblages.     

Global diversity and biogeography of deep-sea pelagic prokaryotes

The deep-sea is the largest biome of the biosphere, and contains more than half of the whole ocean''s microbes. Uncovering their general patterns of diversity and community structure at a global scale remains a great challenge, as only fragmentary information of deep-sea microbial diversity exists based on regional-scale studies. Here we report the first globally comprehensive survey of the prokaryotic communities inhabiting the bathypelagic ocean using high-throughput sequencing of the 16S rRNA gene. This work identifies the dominant prokaryotes in the pelagic deep ocean and reveals that 50% of the operational taxonomic units (OTUs) belong to previously unknown prokaryotic taxa, most of which are rare and appear in just a few samples. We show that whereas the local richness of communities is comparable to that observed in previous regional studies, the global pool of prokaryotic taxa detected is modest (~3600 OTUs), as a high proportion of OTUs are shared among samples. The water masses appear to act as clear drivers of the geographical distribution of both particle-attached and free-living prokaryotes. In addition, we show that the deep-oceanic basins in which the bathypelagic realm is divided contain different particle-attached (but not free-living) microbial communities. The combination of the aging of the water masses and a lack of complete dispersal are identified as the main drivers for this biogeographical pattern. All together, we identify the potential of the deep ocean as a reservoir of still unknown biological diversity with a higher degree of spatial complexity than hitherto considered.     

Evaluation of non-lethal gut microbiome sampling methods in a passerine bird

Gut microbial communities play critical roles in the biological functions of their host, such as mediating nutrient absorption, digesting food components the host cannot, and offering protection against enteric pathogens. Extensive research on gut microbial communities has been conducted on mammals, including humans and rodents, but much less work has been done in birds. Furthermore, much of the research on host–microbe interactions make use of faecal samples and rectal/cloacal swabs as a proxy for intestinal samples, which can be difficult to obtain directly. However, little is known about the overlap between the microbial communities of the gut, faeces and swabs, which limits interpretability of results based on faecal samples and swabs. To address this gap in knowledge, we compared the microbiome from five sample types – proventriculus, small intestine, large intestine, cloacal swabs and faeces – across individual Zebra Finches Taeniopygia guttata housed in constant conditions with a standardized diet. We compared diversity and community composition through 16S rRNA gene sequencing. Our results show that microbial communities from both cloacal swabs and faeces were distinct from proventriculus and small intestinal samples, but generally indistinguishable from large intestinal samples, indicating that these non-lethal samples may be useful proxies for large intestinal bacterial communities. Gaining insight into non-invasive sampling techniques for passerines has implications for studies of gut microbial diversity and abundance in wild bird populations. Furthermore, reliable non-lethal sampling is necessary for experiments where repeated sampling is required.     

Introducing DOTUR, a Computer Program for Defining Operational Taxonomic Units and Estimating Species Richness

Although copious qualitative information describes the members of the diverse microbial communities on Earth, statistical approaches for quantifying and comparing the numbers and compositions of lineages in communities are lacking. We present a method that addresses the challenge of assigning sequences to operational taxonomic units (OTUs) based on the genetic distances between sequences. We developed a computer program, DOTUR, which assigns sequences to OTUs by using either the furthest, average, or nearest neighbor algorithm for each distance level. DOTUR uses the frequency at which each OTU is observed to construct rarefaction and collector's curves for various measures of richness and diversity. We analyzed 16S rRNA gene libraries derived from Scottish and Amazonian soils and the Sargasso Sea with DOTUR, which assigned sequences to OTUs rapidly and reliably based on the genetic distances between sequences and identified previous inconsistencies and errors in assigning sequences to OTUs. An analysis of the two 16S rRNA gene libraries from soil demonstrated that they do not contain enough sequences to support a claim that they contain different numbers of bacterial lineages with statistical confidence (P > 0.05), nor do they contain enough sequences to provide a robust estimate of species richness when an OTU is defined as containing sequences that are no more than 3% different from each other. In contrast, the richness of OTUs at the 3% level in the Sargasso Sea collection began to plateau after the sampling of 690 sequences. We anticipate that an equivalent extent of sampling for soil would require sampling more than 10,000 sequences, almost 100 times the size of typical sequence collections obtained from soil.     

Metabarcoding of fungal communities associated with bark beetles

Many species of fungi are closely allied with bark beetles, including many tree pathogens, but their species richness and patterns of distribution remain largely unknown. We established a protocol for metabarcoding of fungal communities directly from total genomic DNA extracted from individual beetles, showing that the ITS3/4 primer pair selectively amplifies the fungal ITS. Using three specimens of bark beetle from different species, we assess the fungal diversity associated with these specimens and the repeatability of these estimates in PCRs conducted with different primer tags. The combined replicates produced 727 fungal Operational Taxonomic Units (OTUs) for the specimen of Hylastes ater, 435 OTUs for Tomicus piniperda, and 294 OTUs for Trypodendron lineatum, while individual PCR reactions produced on average only 229, 54, and 31 OTUs for the three specimens, respectively. Yet, communities from PCR replicates were very similar in pairwise comparisons, in particular when considering species abundance, but differed greatly among the three beetle specimens. Different primer tags or the inclusion of amplicons in separate libraries did not impact the species composition. The ITS2 sequences were identified with the Lowest Common Ancestor approach and correspond to diverse lineages of fungi, including Ophiostomaceae and Leotiomycetes widely found to be tree pathogens. We conclude that Illumina MiSeq metabarcoding reliably captures fungal diversity associated with bark beetles, although numerous PCR replicates are recommended for an exhaustive sample. Direct PCR from beetle DNA extractions provides a rapid method for future surveys of fungal species diversity and their associations with bark beetles and environmental variables.     

Patterns of fungal communities among and within decaying logs, revealed by 454 sequencing

Owing to previous methodological limitations, knowledge about the fine-scale distribution of fungal mycelia in decaying logs is limited. We investigated fungal communities in decaying Norway spruce logs at various spatial scales at two environmentally different locations in Sweden. On the basis of 454 pyrosequencing of the ITS2 region of rDNA, 1914 operational taxonomic units (OTUs) were detected in 353 samples. The communities differed significantly among logs, but the physical distance between logs was not found to have a significant effect on whether fungal communities had any resemblance to each other. Within a log, samples that were closer together generally had communities that showed more resemblance to each other than those that were further apart. OTUs characteristic for particular positions on the logs could be identified. In general, these OTUs did not overlap with the most abundant OTUs, and their ecological role was often unknown. Only a few OTUs were detected in the majority of logs, whereas numerous OTUs were rare and present in only one or a few logs. Wood-decaying Basidiomycetes were often represented by higher sequence reads in individual logs than Ascomycete OTUs, suggesting that Basidiomycete mycelia spread out more rapidly when established. OTU richness tended to increase with the decay stage of the sample; however, the known wood decayers were most abundant in less-decomposed samples. The fungi identified in the logs represented different ecological strategies. Our findings differ from previously published sporocarp studies, indicating that the highly abundant fruiting species may respond to environment in different ways than the rest of the fungal community.     

Next Generation Sequencing to Define Prokaryotic and Fungal Diversity in the Bovine Rumen

A combination of Sanger and 454 sequences of small subunit rRNA loci were used to interrogate microbial diversity in the bovine rumen of 12 cows consuming a forage diet. Observed bacterial species richness, based on the V1–V3 region of the 16S rRNA gene, was between 1,903 to 2,432 species-level operational taxonomic units (OTUs) when 5,520 reads were sampled per animal. Eighty percent of species-level OTUs were dominated by members of the order Clostridiales, Bacteroidales, Erysipelotrichales and unclassified TM7. Abundance of Prevotella species varied widely among the 12 animals. Archaeal species richness, also based on 16S rRNA, was between 8 and 13 OTUs, representing 5 genera. The majority of archaeal OTUs (84%) found in this study were previously observed in public databases with only two new OTUs discovered. Observed rumen fungal species richness, based on the 18S rRNA gene, was between 21 and 40 OTUs with 98.4–99.9% of OTUs represented by more than one read, using Good’s coverage. Examination of the fungal community identified numerous novel groups. Prevotella and Tannerella were overrepresented in the liquid fraction of the rumen while Butyrivibrio and Blautia were significantly overrepresented in the solid fraction of the rumen. No statistical difference was observed between the liquid and solid fractions in biodiversity of archaea and fungi. The survey of microbial communities and analysis of cross-domain correlations suggested there is a far greater extent of microbial diversity in the bovine rumen than previously appreciated, and that next generation sequencing technologies promise to reveal novel species, interactions and pathways that can be studied further in order to better understand how rumen microbial community structure and function affects ruminant feed efficiency, biofuel production, and environmental impact.     

Inter‐annual and body topographic consistency in the plumage bacterial load of Great Tits

ABSTRACT Plumage bacteria may play an important role in shaping the life histories of birds. However, to design suitable experiments to examine causal relationships between plumage bacteria and the fitness of host birds, natural variation in plumage bacterial communities needs to be better understood. We examined within‐individual consistency of plumage bacterial contamination in Great Tits (Parus major), comparing different body regions (ventral vs. dorsal) and comparing contamination between years. Numbers of free‐living and attached bacteria and the species richness of feather‐degrading bacterial assemblages were studied using flow cytometry and ribosomal intergenic spacer analysis (RISA). Numbers of both types of bacteria were higher on dorsal than on ventral feathers. Numbers of free‐living, but not attached, bacteria on the two body regions were highly positively correlated. There was also a strong within‐individual correlation between numbers of attached bacteria during the same breeding stages in different years. These results suggest that, despite variation in absolute levels of feather bacterial loads between years and different body regions, sampling individual birds can provide reliable estimates of relative levels of bacterial contamination, as long as sampling time and body region are carefully standardized.     

Introducing DOTUR, a computer program for defining operational taxonomic units and estimating species richness

Although copious qualitative information describes the members of the diverse microbial communities on Earth, statistical approaches for quantifying and comparing the numbers and compositions of lineages in communities are lacking. We present a method that addresses the challenge of assigning sequences to operational taxonomic units (OTUs) based on the genetic distances between sequences. We developed a computer program, DOTUR, which assigns sequences to OTUs by using either the furthest, average, or nearest neighbor algorithm for each distance level. DOTUR uses the frequency at which each OTU is observed to construct rarefaction and collector's curves for various measures of richness and diversity. We analyzed 16S rRNA gene libraries derived from Scottish and Amazonian soils and the Sargasso Sea with DOTUR, which assigned sequences to OTUs rapidly and reliably based on the genetic distances between sequences and identified previous inconsistencies and errors in assigning sequences to OTUs. An analysis of the two 16S rRNA gene libraries from soil demonstrated that they do not contain enough sequences to support a claim that they contain different numbers of bacterial lineages with statistical confidence (P > 0.05), nor do they contain enough sequences to provide a robust estimate of species richness when an OTU is defined as containing sequences that are no more than 3% different from each other. In contrast, the richness of OTUs at the 3% level in the Sargasso Sea collection began to plateau after the sampling of 690 sequences. We anticipate that an equivalent extent of sampling for soil would require sampling more than 10,000 sequences, almost 100 times the size of typical sequence collections obtained from soil.     

Microbial community composition but not diversity changes along succession in arctic sand dunes

The generality of increasing diversity of fungi and bacteria across arctic sand dune succession was tested. Microbial communities were examined by high‐throughput sequencing of 16S rRNA genes (bacteria) and internal transcribed spacer (ITS) regions (fungi). We studied four microbial compartments (inside leaf, inside root, rhizosphere and bulk soil) and characterized microbes associated with a single plant species (Deschampsia flexuosa) across two sand dune successional stages (early and late). Bacterial richness increased across succession in bulk soil and leaf endosphere. In contrast, soil fungal richness remained constant while root endosphere fungal richness increased across succession. There was, however, no significant difference in Shannon diversity indices between early and late successional stage in any compartment. There was a significant difference in the composition of microbial communities between early and late successional stage in all compartments, although the major microbial OTUs were shared between early and late successional stage. Co‐occurrence network analysis revealed successional stage‐specific microbial groups. There were more co‐occurring modules in early successional stage than in late stage. Altogether, these results emphasize that succession strongly affects distribution of microbial species, but not microbial diversity in arctic sand dune ecosystem and that fungi and bacteria may not follow the same successional trajectories.