To describe a microbe's physiology, including its metabolism, environmental roles, and growth characteristics, it must be grown in a laboratory culture. Unfortunately, many phylogenetically novel groups have never been cultured, so their physiologies have only been inferred from genomics and environmental characteristics. Although the diversity, or number of different taxonomic groups, of uncultured clades has been studied well, their global abundances, or numbers of cells in any given environment, have not been assessed. We quantified the degree of similarity of 16S rRNA gene sequences from diverse environments in publicly available metagenome and metatranscriptome databases, which we show have far less of the culture bias present in primer-amplified 16S rRNA gene surveys, to those of their nearest cultured relatives. Whether normalized to scaffold read depths or not, the highest abundances of metagenomic 16S rRNA gene sequences belong to phylogenetically novel uncultured groups in seawater, freshwater, terrestrial subsurface, soil, hypersaline environments, marine sediment, hot springs, hydrothermal vents, nonhuman hosts, snow, and bioreactors (22% to 87% uncultured genera to classes and 0% to 64% uncultured phyla). The exceptions were human and human-associated environments, which were dominated by cultured genera (45% to 97%). We estimate that uncultured genera and phyla could comprise 7.3 × 1029 (81%) and 2.2 × 1029 (25%) of microbial cells, respectively. Uncultured phyla were overrepresented in metatranscriptomes relative to metagenomes (46% to 84% of sequences in a given environment), suggesting that they are viable. Therefore, uncultured microbes, often from deeply phylogenetically divergent groups, dominate nonhuman environments on Earth, and their undiscovered physiologies may matter for Earth systems. IMPORTANCE In the past few decades, it has become apparent that most of the microbial diversity on Earth has never been characterized in laboratory cultures. We show that these unknown microbes, sometimes called "microbial dark matter," are numerically dominant in all major environments on Earth, with the exception of the human body, where most of the microbes have been cultured. We also estimate that about one-quarter of the population of microbial cells on Earth belong to phyla with no cultured relatives, suggesting that these never-before-studied organisms may be important for ecosystem functions. Author Video: An author video summary of this article is available.

The central aims of many host or environmental microbiome studies are to elucidate factors associated with microbial community compositions and to relate microbial features to outcomes. However, these aims are often complicated by difficulties stemming from high-dimensionality, non-normality, sparsity, and the compositional nature of microbiome data sets. A key tool in microbiome analysis is beta diversity, defined by the distances between microbial samples. Many different distance metrics have been proposed, all with varying discriminatory power on data with differing characteristics. Here, we propose a compositional beta diversity metric rooted in a centered log-ratio transformation and matrix completion called robust Aitchison PCA. We demonstrate the benefits of compositional transformations upstream of beta diversity calculations through simulations. Additionally, we demonstrate improved effect size, classification accuracy, and robustness to sequencing depth over the current methods on several decreased sample subsets of real microbiome data sets. Finally, we highlight the ability of this new beta diversity metric to retain the feature loadings linked to sample ordinations revealing salient intercommunity niche feature importance. IMPORTANCE By accounting for the sparse compositional nature of microbiome data sets, robust Aitchison PCA can yield high discriminatory power and salient feature ranking between microbial niches. The software to perform this analysis is available under an open-source license and can be obtained at https://github.com/biocore/DEICODE; additionally, a QIIME 2 plugin is provided to perform this analysis at https://library.qiime2.org/plugins/deicode/.

Background

The native microflora associated with mosquitoes have important roles in mosquito development and vector competence. Sequencing of bacterial V3 region from 16S rRNA genes across the developmental stages of Culex mosquitoes (early and late larval instars, pupae and adults) was used to test the hypothesis that bacteria found in the larval stage of Culex are transstadially transmitted to the adult stage, and to compare the microbiomes of field-collected versus laboratory-reared mosquitoes.

Results

Beta diversity analysis revealed that bacterial community structure differed among three life stages (larvae, pupae and adults) of Culex tarsalis. Although only ~2% of the total number of bacterial OTUs were found in all stages, sequences from these OTUs accounted for nearly 82% of the total bacterial sequences recovered from all stages. Thorsellia (Gammaproteobacteria) was the most abundant bacterial taxon found across all developmental stages of field-collected Culex mosquitoes, but was rare in mosquitoes from laboratory-reared colonies. The proportion of Thorsellia sequences in the microbiomes of mosquito life stages varied ontogenetically with the greatest proportions recovered from the pupae of C. tarsalis and the lowest from newly emerged adults. The microbiome of field-collected late instar larvae was not influenced significantly by differences in the microbiota of the habitat due to habitat age or biopesticide treatments. The microbiome diversity was the greatest in the early instar larvae and the lowest in laboratory-reared mosquitoes.

Conclusions

Bacterial communities in early instar C. tarsalis larvae were significantly more diverse when compared to late instar larvae, pupae and newly emerged adults. Some of the bacterial OTUs found in the early instar larvae were also found across developmental stages. Thorsellia dominated the bacterial communities in field-collected immature stages but occurred at much lower relative abundance in adults. Differences in microbiota observed in larval habitats did not influence bacterial community profiles of late instar larvae or adults. However, bacterial communities in laboratory-reared C. tarsalis larvae differed significantly from the field. Determining the role of Thorsellia in mosquitoes and its distribution across different species of mosquitoes warrants further investigation.

Differential abundance testing is a critical task in microbiome studies that is complicated by the sparsity of data matrices. Here we adapt for microbiome studies a solution from the field of gene expression analysis to produce a new method, discrete false-discovery rate (DS-FDR), that greatly improves the power to detect differential taxa by exploiting the discreteness of the data. Additionally, DS-FDR is relatively robust to the number of noninformative features, and thus removes the problem of filtering taxonomy tables by an arbitrary abundance threshold. We show by using a combination of simulations and reanalysis of nine real-world microbiome data sets that this new method outperforms existing methods at the differential abundance testing task, producing a false-discovery rate that is up to threefold more accurate, and halves the number of samples required to find a given difference (thus increasing the efficiency of microbiome experiments considerably). We therefore expect DS-FDR to be widely applied in microbiome studies. IMPORTANCE DS-FDR can achieve higher statistical power to detect significant findings in sparse and noisy microbiome data compared to the commonly used Benjamini-Hochberg procedure and other FDR-controlling procedures.

Microbes are important for mosquito nutrition, growth, reproduction and control. In this study, we examined bacterial communities associated with larval mosquitoes and their habitats. Specifically, we characterized bacterial communities associated with late larval instars of the western encephalitis mosquito (Culextarsalis), the submerged portions of two emergent macrophytes (California bulrush, Schoenoplectuscalifornicus and alkali bulrush, Schoenoplectusmaritimus), and the associated water columns to investigate potential differential use of resources by mosquitoes in different wetland habitats. Using next-generation sequence data from 16S rRNA gene hypervariable regions, the alpha diversity of mosquito gut microbial communities did not differ between pond mesocosms containing distinct monotypic plants. Proteobacteria, dominated by the genus Thorsellia (Enterobacteriaceae), was the most abundant phylum recovered from C. tarsalis larvae. Approximately 49% of bacterial OTUs found in larval mosquitoes were identical to OTUs recovered from the water column and submerged portions of the two bulrushes. Plant and water samples were similar to one another, both being dominated by Actinobacteria, Bacteroidetes, Cyanobacteria, Proteobacteria and Verrucomicrobia phyla. Overall, the bacterial communities within C. tarsalis larvae were conserved and did not change across sampling dates and between two distinct plant habitats. Although Thorsellia spp. dominated mosquito gut communities, overlap of mosquito gut, plant and water-column OTUs likely reveal the effects of larval feeding. Future research will investigate the role of the key indicator groups of bacteria across the different developmental stages of this mosquito species.

Mock communities are an important tool for validating, optimizing, and comparing bioinformatics methods for microbial community analysis. We present mockrobiota, a public resource for sharing, validating, and documenting mock community data resources, available at http://caporaso-lab.github.io/mockrobiota/. The materials contained in mockrobiota include data set and sample metadata, expected composition data (taxonomy or gene annotations or reference sequences for mock community members), and links to raw data (e.g., raw sequence data) for each mock community data set. mockrobiota does not supply physical sample materials directly, but the data set metadata included for each mock community indicate whether physical sample materials are available. At the time of this writing, mockrobiota contains 11 mock community data sets with known species compositions, including bacterial, archaeal, and eukaryotic mock communities, analyzed by high-throughput marker gene sequencing. IMPORTANCE The availability of standard and public mock community data will facilitate ongoing method optimizations, comparisons across studies that share source data, and greater transparency and access and eliminate redundancy. These are also valuable resources for bioinformatics teaching and training. This dynamic resource is intended to expand and evolve to meet the changing needs of the omics community.

Although aquaria are common features of homes and other buildings, little is known about how environmental perturbations (i.e., tank cleaning, water changes, addition of habitat features) impact the diversity and succession of aquarium microbial communities. In this study, we sought to evaluate the hypotheses that newly established aquaria show clear microbial successional patterns over time and that common marine aquarium-conditioning practices, such as the addition of ocean-derived "live rocks" (defined as any "dead coral skeleton covered with crustose coralline algae" transferred into an aquarium from open ocean habitats) impact the diversity of microbial populations as well as nitrogen cycling in aquaria. We collected water chemistry data alongside water and sediment samples from two independent and newly established saltwater aquaria over a 3-month period. Microbial communities in samples were assessed by DNA extraction, amplification of the 16S rRNA gene, and Illumina MiSeq sequencing. Our results showed clear and replicable patterns of community succession in both aquaria, with the existence of multiple stable states for aquarium microbial assemblages. Notably, our results show that changes in aquarium microbial communities do not always correlate with water chemistry measurements and that operational taxonomic unit (OTU)-level patterns relevant to nitrogen cycling were not reported as statistically significant. Overall, our results demonstrate that aquarium perturbations have a substantial impact on microbial community profiles of aquarium water and sediment and that the addition of live rocks improves nutrient cycling by shifting aquarium communities toward a more typical saltwater assemblage of microbial taxa.IMPORTANCE Saltwater aquaria are living systems that support a complex biological community of fish, invertebrates, and microbes. The health and maintenance of saltwater tanks are pressing concerns for home hobbyists, zoos, and professionals in the aquarium trade; however, we do not yet understand the underlying microbial species interactions and community dynamics which contribute to tank setup and conditioning. This report provides a detailed view of ecological succession and changes in microbial community assemblages in two saltwater aquaria which were sampled over a 3-month period, from initial tank setup and conditioning with "live rocks" through subsequent tank cleanings and water replacement. Our results showed that microbial succession appeared to be consistent and replicable across both aquaria. However, changes in microbial communities did not always correlate with water chemistry measurements, and aquarium microbial communities appear to have shifted among multiple stable states without any obvious buildup of undesirable nitrogen compounds in the tank environment.