Modeling Taxa-Abundance Distributions in Microbial Communities using Environmental Sequence Data

We show that inferring the taxa-abundance distribution of a microbial community from small environmental samples alone is difficult. The difficulty stems from the disparity in scale between the number of genetic sequences that can be characterized and the number of individuals in communities that microbial ecologists aspire to describe. One solution is to calibrate and validate a mathematical model of microbial community assembly using the small samples and use the model to extrapolate to the taxa-abundance distribution for the population that is deemed to constitute a community. We demonstrate this approach by using a simple neutral community assembly model in which random immigrations, births, and deaths determine the relative abundance of taxa in a community. In doing so, we further develop a neutral theory to produce a taxa-abundance distribution for large communities that are typical of microbial communities. In addition, we highlight that the sampling uncertainties conspire to make the immigration rate calibrated on the basis of small samples very much higher than the true immigration rate. This scale dependence of model parameters is not unique to neutral theories; it is a generic problem in ecology that is particularly acute in microbial ecology. We argue that to overcome this, so that microbial ecologists can characterize large microbial communities from small samples, mathematical models that encapsulate sampling effects are required.     

Biogeography at the limits of life: Do extremophilic microbial communities show biogeographical regionalization?

Aim

Biogeographical regions are the fundamental geographical units for grouping Earth's biodiversity. Biogeographical regionalization has been demonstrated for many higher taxa, such as terrestrial plants and vertebrates, but not in microbial communities. Therefore, we sought to test empirically whether microbial communities, or taxa, show patterns consistent with biogeographical regionalization.

Location

Within halite (NaCl) crystals from coastal solar salterns of western Europe, the Mediterranean and east Africa.

Time period

Modern (2006–2013).

Major taxa studied

Archaea.

Methods

Using high‐throughput Illumina amplicon sequencing, we generated the most high‐resolution characterization of halite‐associated archaeal communities to date, using samples from 17 locations. We grouped communities into biogeographical clusters based on community turnover to test whether these communities show biogeographical regionalization. To examine whether individual taxa, rather than communities, show biogeographical patterns, we also tested whether the relative abundance of individual genera may be indicative of a community's biogeographical origins using machine learning methods, specifically random forest classification.

Results

We found that the rate of community turnover was greatest over subregional spatial scales (< 500 km), whereas at regional spatial scales the turnover was independent of geographical distance. Biogeographical clusters of communities were either not statistically robust or lacked spatial coherence, inconsistent with biogeographical regionalization. However, we identified several archaeal genera that were good indicators of biogeographical origin, providing classification error rates of < 10%.

Main conclusions

Overall, our results provide little support for the concept of biogeographical regions in these extremophilic microbial communities, despite the fact that some taxa do show biogeographical patterns. We suggest that variable dispersal ability among the halite‐associated Archaea may disrupt biogeographical patterns at the community level, preventing the formation of biogeographical regions. This means that the processes that lead to the formation of biogeographical regions operate differentially on individual microbial taxa rather than on entire communities.     

Identification of Specialists and Abundance-Occupancy Relationships among Intestinal Bacteria of Aves,Mammalia, and Actinopterygii

The coalescence of next-generation DNA sequencing methods, ecological perspectives, and bioinformatics analysis tools is rapidly advancing our understanding of the evolution and function of vertebrate-associated bacterial communities. Delineation of host-microbe associations has applied benefits ranging from clinical treatments to protecting our natural waters. Microbial communities follow some broad-scale patterns observed for macroorganisms, but it remains unclear how the specialization of intestinal vertebrate-associated communities to a particular host environment influences broad-scale patterns in microbial abundance and distribution. We analyzed the V6 region of 16S rRNA genes amplified from 106 fecal samples spanning Aves, Mammalia, and Actinopterygii (ray-finned fish). We investigated the interspecific abundance-occupancy relationship, where widespread taxa tend to be more abundant than narrowly distributed taxa, among operational taxonomic units (OTUs) within and among host species. In a separate analysis, we identified specialist OTUs that were highly abundant in a single host and rare in all other hosts by using a multinomial model without excluding undersampled OTUs a priori. We show that intestinal microbes in humans and other vertebrates display abundance-occupancy relationships, but because intestinal host-associated communities have undergone intense specialization, this trend is violated by a disproportionately large number of specialist taxa. Although it is difficult to distinguish the effects of dispersal limitations, host selection, historical contingency, and stochastic processes on community assembly, results suggest that intestinal bacteria can be shared among diverse hosts in ways that resemble the distribution of “free-living” bacteria in the extraintestinal environment.     

Effect of PCR amplicon size on assessments of clone library microbial diversity and community structure

PCR-based surveys of microbial communities commonly use regions of the small-subunit ribosomal RNA (SSU rRNA) gene to determine taxonomic membership and estimate total diversity. Here we show that the length of the target amplicon has a significant effect on assessments of microbial richness and community membership. Using operational taxonomic unit (OTU)- and taxonomy-based tools, we compared the V6 hypervariable region of the bacterial SSU rRNA gene of three amplicon libraries of c. 100, 400 and 1000 base pairs (bp) from each of two hydrothermal vent fluid samples. We found that the smallest amplicon libraries contained more unique sequences, higher diversity estimates and a different community structure than the other two libraries from each sample. We hypothesize that a combination of polymerase dissociation, cloning bias and mispriming due to secondary structure accounts for the differences. While this relationship is not linear, it is clear that the smallest amplicon libraries contained more different types of sequences, and accordingly, more diverse members of the community. Because divergent and lower abundant taxa can be more readily detected with smaller amplicons, they may provide better assessments of total community diversity and taxonomic membership than longer amplicons in molecular studies of microbial communities.     

Fatty acid methyl ester (FAME) profiles as measures of soil microbial community structure

Analysis of fatty acid methyl ester (FAME) profiles extracted from soils is a rapid and inexpensive procedure that holds great promise in describing soil microbial community structure without traditional reliance on selective culturing, which seems to severely underestimate community diversity. Interpretation of FAME profiles from environmental samples can be difficult because many fatty acids are common to different microorganisms and many fatty acids are extracted from each soil sample. We used principal components (PCA) and cluster analyses to identify similarities and differences among soil microbial communities described using FAME profiles. We also used PCA to identify particular FAMEs that characterized soil sample clusters. Fatty acids that are found only or primarily in particular microbial taxa-marker fatty acids-were used in conjunction with these analyses. We found that the majority of 162 soil samples taken from a conventionally-tilled corn field had similar FAME profiles but that about 20% of samples seemed to have relatively low, and that about 10% had relatively high, bacterial:fungal ratios. Using semivariance analysis we identified 21:0 iso as a new marker fatty acid. Concurrent use of geostatistical and FAME analyses may be a powerful means of revealing other potential marker FAMEs. When microbial communities from the same samples were cultured on R2A agar and their FAME profiles analyzed, there were many differences between FAME profiles of soil and plated communities, indicating that profiles of FAMEs extracted from soil reveal portions of the microbial community not culturable on R2A. When subjected to PCA, however, a small number of plated communities were found to be distinct due to some of the same profile characteristics (high in 12:0 iso, 15:0 and 17:1 ante A) that identified soil community FAME profiles as distinct. Semivariance analysis indicated that spatial distributions of soil microbial populations are maintained in a portion of the microbial community that is selected on laboratory media. These similarities between whole soil and plated community FAME profiles suggest that plated communities are not solely the result of selection by the growth medium, but reflect the distribution, in situ, of the dominant, culturable soil microbial populations.     

Local structuring factors of invertebrate communities in ephemeral freshwater rock pools and the influence of more permanent water bodies in the region

We used three isolated clusters of small ephemeral rock pools on a sandstone flat in Utah to test the importance of local structuring processes on aquatic invertebrate communities. In the three clusters we characterized all ephemeral rock pools (total: 27) for their morphometry, and monitored their water quality, hydrology and community assemblage during a full hydrocycle. In each cluster we also sampled a set of more permanent interconnected freshwater systems positioned in a wash, draining the water from each cluster of rock pools. This design allowed additional testing for the potential role of more permanent water bodies in the region as source populations for the active dispersers and the effect on the community structure in the rock pools. Species richness and community composition in the rock pools correlated with level of permanence and the ammonia concentration. The length of the rock pool inundation cycle shaped community structure, most probably by inhibiting colonization by some taxa (e.g. tadpoles and insect larvae) through developmental constraints. The gradient in ammonia concentrations probably reflects differences in primary production. The more permanent water bodies in each wash differed both environmentally and in community composition from the connected set of rock pools. A limited set of active dispersers was observed in the rock pools. Our findings indicate that aquatic invertebrate communities in the ephemeral rock pools are mainly structured through habitat permanence, possibly linked with biotic interactions and primary production. Handling editor: S. I. Dodson     

Changes amid constancy: Flower and leaf microbiomes along land use gradients and between bioregions

Microbial communities inhabiting above-ground parts of plants affect their host's development, fitness and function. Although studies on plant-associated microbes are of growing interest, environmental drivers of flower microbiomes in particular are poorly characterized. In this study, we investigated flower and leaf epiphytic bacterial microbiomes of Ranunculus acris and Trifolium pratense using metabarcoding of 16S ribosomal DNA in three German bioregions and along land-use intensity gradients. Our data suggests that the structures of bacterial communities clearly differed between plant species and tissue types. Also, floral bacterial communities of R. acris showed higher variability in comparison to T. pratense. Bacteria usually associated with pollinators were found solely in flower samples, while bacteria usually associated with the rhizosphere were only present in high abundances on leaves. We identified Pseudomonadaceae, Enterobacteriaceae and Sphingomonadaceae as the most abundant taxa on flowers, while Sphingomonadaceae, Methylobacteriaceae and Cytophagaceae dominated bacterial communities on leaves. We found that bacterial communities did not differ between long-distant regions. However, there was a turnover within each bioregion between short-distant locations. High land use intensity caused phylogenetically less diverse and more homogenous bacterial communities with an exception of T. pratense flowers. This was associated with a loss of rare bacterial families. Intensification of mowing affected the bacterial communities associated with leaves of T. pratense and fertilization led to more homogenous flower and leaf communities of R. acris, while grazing had no effects on the bacterial community composition. However, dominant taxa were not affected by land use intensification. Despite that, we identified indicator taxa for regularly disturbed environments in flower microbiomes. In conclusion, our study contributes to the knowledge about microbial community structures of the phyllosphere and extends the understanding of their community dynamics with respect to biogeographical separation and anthropogenic changes of the environment.     

Use of multiplex terminal restriction fragment length polymorphism for rapid and simultaneous analysis of different components of the soil microbial community

A multiplex terminal restriction fragment length polymorphism (M-TRFLP) fingerprinting method was developed and validated for simultaneous analysis of the diversity and community structure of two or more microbial taxa (up to four taxa). The reproducibility and robustness of the method were examined using soil samples collected from different habitats. DNA was PCR amplified separately from soil samples using individual taxon-specific primers for bacteria, archaea, and fungi. The same samples were also subjected to a multiplex PCR with the primers for all three taxa. The terminal restriction fragment length polymorphism profiles generated for the two sets of PCR products were almost identical not only in terms of the presence of peaks but also in terms of the relative peak intensity. The M-TRFLP method was then used to investigate rhizosphere bacterial, fungal, and rhizobial/agrobacterial communities associated with the dwarf shrub Calluna vulgaris growing in either open moorland, a mature pine forest, or a transition zone between these two habitats containing naturally regenerating pine trees. Rhizosphere microbial communities associated with Vaccinium myrtillus collected from the native pine forest were also investigated. In this study, individual PCR products from the three taxa were also pooled before restriction digestion and fragment size analysis. The terminal restriction fragment length polymorphism profiles obtained with PCR products amplified individually and with multiplexed and pooled PCR products were found to be consistent with each other in terms of the number, position, and relative intensity of peaks. The results presented here confirm that M-TRFLP analysis is a highly reproducible and robust molecular tool for simultaneous investigation of multiple taxa, which allows more complete and higher resolution of microbial communities to be obtained more rapidly and economically.     

Microbiome Profiles in Periodontitis in Relation to Host and Disease Characteristics

Periodontitis is an inflammatory condition that affects the supporting tissues surrounding teeth. The occurrence of periodontitis is associated with shifts in the structure of the communities that inhabit the gingival sulcus. Although great inter-subject variability in the subgingival microbiome has been observed in subjects with periodontitis, it is unclear whether distinct community types exist and if differences in microbial signatures correlate with host characteristics or with the variable clinical presentations of periodontitis. Therefore, in this study we explored the existence of different community types in periodontitis and their relationship with host demographic, medical and disease-related clinical characteristics. Clustering analyses of microbial abundance profiles suggested two types of communities (A and B) existed in the 34 subjects with periodontitis evaluated. Type B communities harbored greater proportions of certain periodontitis-associated taxa, including species historically associated with the disease, such as Porphyromonas gingivalis, Tannerella forsythia and Treponema denticola, and taxa recently linked to periodontitis. In contrast, subjects with type A communities had increased proportions of different periodontitis-associated species, and were also enriched for health-associated species and core taxa (those equally prevalent in health and periodontitis). Periodontitis subgingival clusters were not associated with demographic, medical or disease-specific clinical parameters other than periodontitis extent (proportion of sites affected), which positively correlated with the total proportion of cluster B signature taxa. In conclusion, two types of microbial communities were detected in subjects with periodontitis. Host demographics and underlying medical conditions did not correlate with these profiles, which instead appeared to be related to periodontitis extent, with type B communities present in more widespread disease cases. The two identified periodontitis profiles may represent distinct dysbiotic processes potentially requiring community-tailored therapeutic interventions.     

The effect of patch disturbance on stream invertebrate community structure: the influence of disturbance history

The effect of disturbance history on the recovery of benthic invertebrate communities following disturbance was investigated in four streams in the Southern Alps of New Zealand. Two of the streams had a history of fluctuating discharge and temperature while the others did not. Recovery from disturbance was tested experimentally using baskets of cobbles, a third of which were disturbed every week for 9 weeks, a further third every 3 weeks and the final third left undisturbed. Algal biiomass, number of invertebrate taxa and total number of invertebrates all declined in baskets disturbed more frequently. Although the relative abundance of some taxa declined with time since the last disturbance, no taxa showed a significant decline in absolute abundance. However, several taxa showed marked increases in relative abundance in the less disturbed treatments particularly at the more stable sites. In contrast to the predictions of ecological theory, numbers of taxa and total invertebrates appeared to recover more quickly in the more complex communities at the stable sites. However, if these communities are considered to represent only stable communities, they do support the view that more complex communities will be more resilient. Community structure at the stable sites was also more similar between baskets in the undisturbed treatment than at the unstable sites, suggesting communities had reached a constant state more quickly. The more rapid recovery of communities measured at the stable sites may have been a consequence of experimental scale; disturbed patches were only 0.045 m² in area and the higher densities of invertebrates at the stable sites meant a larger pool of colonists was available following each experimental disturbance. Nevertheless, ideas of stability in ecological theory and the scale of most spate events suggest this is the appropriate scale for examining community recovery. Furthermore, the larger pool of available colonists could not explain all the differences in community response, as patterns of change in community structure at the stable sites differed considerably more from those expected by purely random colonisation processes than at the unstable sites.     

A statistical procedure for the analysis of microbial communities based on phenotypic properties of isolates

A novel statistical procedure for the analysis of microbial communities based on phenotypic properties of randomly collected isolates is presented and discussed. The procedure allows the representation of the microbial communities as a set of ellipses in a bidimensional graph. This representation is obtained by the following steps: (a) measurement of a set of binary phenotypic properties for n isolates belonging to k samples, each representing a different community; (b) repeated sampling by bootstrapping of the m samples, thus obtaining, for each community, i subsamples of j isolates; (c) calculation of the frequency of positive results for each test for each subsample; (d) calculation of the matrix of Euclidean distances between the k x i frequency vectors; (e) use of multidimensional scaling (MDS) to obtain a representation in two dimensions of the distance relationships between the frequency vectors; (f) plotting of the 95% confidence ellipses for the i frequency vectors of each of the k communities. By using both simple, synthetic microbial communities, and samples of lactic acid bacteria isolated from natural microbial communities (sourdoughs, compressed yeast, fermented sausages), it was demonstrated that the position and shape of the ellipses are clearly related to the composition of the community, while the relationship between the size of the ellipses and the phenotypical diversity of the community is less straightforward: while communities with very different diversity (measured with the Functional Evenness index and the mean taxonomic distance) had ellipses that were very different in size, there was no strict proportionality between the size of the ellipse and the diversity of the community. Nevertheless, the representation of microbial communities obtained by bootstrapping and multidimensional scaling appears to be superior to the more usual representation based on tabulation of the frequencies of isolates belonging to different clusters.     

Development of the human infant intestinal microbiota

Almost immediately after a human being is born, so too is a new microbial ecosystem, one that resides in that person's gastrointestinal tract. Although it is a universal and integral part of human biology, the temporal progression of this process, the sources of the microbes that make up the ecosystem, how and why it varies from one infant to another, and how the composition of this ecosystem influences human physiology, development, and disease are still poorly understood. As a step toward systematically investigating these questions, we designed a microarray to detect and quantitate the small subunit ribosomal RNA (SSU rRNA) gene sequences of most currently recognized species and taxonomic groups of bacteria. We used this microarray, along with sequencing of cloned libraries of PCR-amplified SSU rDNA, to profile the microbial communities in an average of 26 stool samples each from 14 healthy, full-term human infants, including a pair of dizygotic twins, beginning with the first stool after birth and continuing at defined intervals throughout the first year of life. To investigate possible origins of the infant microbiota, we also profiled vaginal and milk samples from most of the mothers, and stool samples from all of the mothers, most of the fathers, and two siblings. The composition and temporal patterns of the microbial communities varied widely from baby to baby. Despite considerable temporal variation, the distinct features of each baby's microbial community were recognizable for intervals of weeks to months. The strikingly parallel temporal patterns of the twins suggested that incidental environmental exposures play a major role in determining the distinctive characteristics of the microbial community in each baby. By the end of the first year of life, the idiosyncratic microbial ecosystems in each baby, although still distinct, had converged toward a profile characteristic of the adult gastrointestinal tract.     

Antibiotic disturbance affects aquatic microbial community composition and food web interactions but not community resilience

Notwithstanding the fundamental role that environmental microbes play for ecosystem functioning, data on how microbes react to disturbances are still scarce, and most factors that confer stability to microbial communities are unknown. In this context, antibiotic discharge into the environment is considered a worldwide threat for ecosystems with potential risks to human health. We therefore tested resilience of microbial communities challenged by the presence of an antibiotic. In a continuous culture experiment, we compared the abundance, composition and diversity of microbial communities undisturbed or disturbed by the constant addiction of tetracycline in low (10 µg/L) or intermediate (100 µg/L) concentration (press disturbance). Further, the bacterial communities in the three treatments had to face the sudden pulse disturbance of adding an allochthonous bacterium (Escherichia coli). Tetracycline, even at low concentrations, affected microbial communities by changing their phylogenetic composition and causing cell aggregation. This, however, did not coincide with a reduced microbial diversity, but was mainly caused by a shift in dominance of specific bacterial families. Moreover, the less disturbed community (10 µg/L tetracycline) was sometimes more similar to the control and sometimes more similar to heavily disturbed community (100 µg/L tetracycline). All in all, we could not see a pattern where the communities disturbed with antibiotics were less resilient to a second disturbance introducing E. coli, but they seemed to be able to buffer the input of the allochthonous strain in a similar manner as the control.     

VITCOMIC2: visualization tool for the phylogenetic composition of microbial communities based on 16S rRNA gene amplicons and metagenomic shotgun sequencing

Background

The 16S rRNA gene-based amplicon sequencing analysis is widely used to determine the taxonomic composition of microbial communities. Once the taxonomic composition of each community is obtained, evolutionary relationships among taxa are inferred by a phylogenetic tree. Thus, the combined representation of taxonomic composition and phylogenetic relationships among taxa is a powerful method for understanding microbial community structure; however, applying phylogenetic tree-based representation with information on the abundance of thousands or more taxa in each community is a difficult task. For this purpose, we previously developed the tool VITCOMIC (VIsualization tool for Taxonomic COmpositions of MIcrobial Community), which is based on the genome-sequenced microbes’ phylogenetic information. Here, we introduce VITCOMIC2, which incorporates substantive improvements over VITCOMIC that were necessary to address several issues associated with 16S rRNA gene-based analysis of microbial communities.

Results

We developed VITCOMIC2 to provide (i) sequence identity searches against broad reference taxa including uncultured taxa; (ii) normalization of 16S rRNA gene copy number differences among taxa; (iii) rapid sequence identity searches by applying the graphics processing unit-based sequence identity search tool CLAST; (iv) accurate taxonomic composition inference and nearly full-length 16S rRNA gene sequence reconstructions for metagenomic shotgun sequencing; and (v) an interactive user interface for simultaneous representation of the taxonomic composition of microbial communities and phylogenetic relationships among taxa. We validated the accuracy of processes (ii) and (iv) by using metagenomic shotgun sequencing data from a mock microbial community.

Conclusions

The improvements incorporated into VITCOMIC2 enable users to acquire an intuitive understanding of microbial community composition based on the 16S rRNA gene sequence data obtained from both metagenomic shotgun and amplicon sequencing.

     

Use of Multiplex Terminal Restriction Fragment Length Polymorphism for Rapid and Simultaneous Analysis of Different Components of the Soil Microbial Community▿

A multiplex terminal restriction fragment length polymorphism (M-TRFLP) fingerprinting method was developed and validated for simultaneous analysis of the diversity and community structure of two or more microbial taxa (up to four taxa). The reproducibility and robustness of the method were examined using soil samples collected from different habitats. DNA was PCR amplified separately from soil samples using individual taxon-specific primers for bacteria, archaea, and fungi. The same samples were also subjected to a multiplex PCR with the primers for all three taxa. The terminal restriction fragment length polymorphism profiles generated for the two sets of PCR products were almost identical not only in terms of the presence of peaks but also in terms of the relative peak intensity. The M-TRFLP method was then used to investigate rhizosphere bacterial, fungal, and rhizobial/agrobacterial communities associated with the dwarf shrub Calluna vulgaris growing in either open moorland, a mature pine forest, or a transition zone between these two habitats containing naturally regenerating pine trees. Rhizosphere microbial communities associated with Vaccinium myrtillus collected from the native pine forest were also investigated. In this study, individual PCR products from the three taxa were also pooled before restriction digestion and fragment size analysis. The terminal restriction fragment length polymorphism profiles obtained with PCR products amplified individually and with multiplexed and pooled PCR products were found to be consistent with each other in terms of the number, position, and relative intensity of peaks. The results presented here confirm that M-TRFLP analysis is a highly reproducible and robust molecular tool for simultaneous investigation of multiple taxa, which allows more complete and higher resolution of microbial communities to be obtained more rapidly and economically.     

Rare but active taxa contribute to community dynamics of benthic biofilms in glacier‐fed streams

Glaciers harbour diverse microorganisms, which upon ice melt can be released downstream. In glacier‐fed streams microorganisms can attach to stones or sediments to form benthic biofilms. We used 454‐pyrosequencing to explore the bulk (16S rDNA) and putatively active (16S rRNA) microbial communities of stone and sediment biofilms across 26 glacier‐fed streams. We found differences in community composition between bulk and active communities among streams and a stronger congruence between biofilm types. Relative abundances of rRNA and rDNA were positively correlated across different taxa and taxonomic levels, but at lower taxonomic levels, the higher abundance in either the active or the bulk communities became more apparent. Here, environmental variables played a minor role in structuring active communities. However, we found a large number of rare taxa with higher relative abundances in rRNA compared with rDNA. This suggests that rare taxa contribute disproportionately to microbial community dynamics in glacier‐fed streams. Our findings propose that high community turnover, where taxa repeatedly enter and leave the ‘seed bank’, contributes to the maintenance of microbial biodiversity in harsh ecosystems with continuous environmental perturbations, such as glacier‐fed streams.     

Accuracy, Reproducibility, and Interpretation of Fatty Acid Methyl Ester Profiles of Model Bacterial Communities

We determined the accuracy and reproducibility of whole-community fatty acid methyl ester (FAME) analysis with two model bacterial communities differing in composition by using the Microbial ID, Inc. (MIDI), system. The biomass, taxonomic structure, and expected MIDI-FAME profiles under a variety of environmental conditions were known for these model communities a priori. Not all members of each community could be detected in the composite profile because of lack of fatty acid “signatures” in some isolates or because of variations (approximately fivefold) in fatty acid yield across taxa. MIDI-FAME profiles of replicate subsamples of a given community were similar in terms of fatty acid yield per unit of community dry weight and relative proportions of specific fatty acids. Principal-components analysis (PCA) of MIDI-FAME profiles resulted in a clear separation of the two different communities and a clustering of replicates of each community from two separate experiments on the first PCA axis. The first PCA axis accounted for 57.1% of the variance in the data and was correlated with fatty acids that varied significantly between communities and reflected the underlying community taxonomic structure. On the basis of our data, community fatty acid profiles can be used to assess the relative similarities and differences of microbial communities that differ in taxonomic composition. However, detailed interpretation of community fatty acid profiles in terms of biomass or community taxonomic composition must be viewed with caution until our knowledge of the quantitative and qualitative distribution of fatty acids over a wide variety of taxa and the effects of growth conditions on fatty acid profiles is more extensive.     

Unexpected biodiversity of ciliates in marine samples from below the photic zone

Marine microbial eukaryotes play critical roles in planktonic food webs and have been described as most diverse in the photic zone where productivity is high. We used high‐throughput sequencing (HTS) to analyse the spatial distribution of planktonic ciliate diversity from shallow waters (<30 m depth) to beyond the continental shelf (>800 m depth) along a 163 km transect off the coast of New England, USA. We focus on ciliates in the subclasses Oligotrichia and Choreotrichia (class Spirotrichea), as these taxa are major components of marine food webs. We did not observe the decrease of diversity below the photic zone expected based on productivity and previous analyses. Instead, we saw an increase of diversity with depth. We also observed that the ciliate communities assessed by HTS cluster by depth layer and degree of water column stratification, suggesting that community assembly is driven by environmental factors. Across our samples, abundant OTUs tend to match previously characterized morphospecies while rare OTUs are more often undescribed, consistent with the idea that species in the rare biosphere remain to be characterized by microscopy. Finally, samples taken below the photic zone also reveal the prevalence of two uncharacterized (i.e. lacking sequenced morphospecies) clades – clusters X₁ and X₂ – that are enriched within the nano‐sized fraction (2–10 μm) and are defined by deletions within the region of the SSU‐rDNA analysed here. Together, these data reinforce that we still have much to learn about microbial diversity in marine ecosystems, especially in deep‐waters that may be a reservoir for rare species and uncharacterized taxa.