Prospecting environmental mycobacteria: combined molecular approaches reveal unprecedented diversity

Background: Environmental mycobacteria (EM) include species commonly found in various terrestrial and aquatic environments, encompassing animal and human pathogens in addition to saprophytes. Approximately 150 EM species can be separated into fast and slow growers based on sequence and copy number differences of their 16S rRNA genes. Cultivation methods are not appropriate for diversity studies; few studies have investigated EM diversity in soil despite their importance as potential reservoirs of pathogens and their hypothesized role in masking or blocking M. bovis BCG vaccine. Methods: We report here the development, optimization and validation of molecular assays targeting the 16S rRNA gene to assess diversity and prevalence of fast and slow growing EM in representative soils from semi tropical and temperate areas. New primer sets were designed also to target uniquely slow growing mycobacteria and used with PCR-DGGE, tag-encoded Titanium amplicon pyrosequencing and quantitative PCR. Results: PCR-DGGE and pyrosequencing provided a consensus of EM diversity; for example, a high abundance of pyrosequencing reads and DGGE bands corresponded to M. moriokaense, M. colombiense and M. riyadhense. As expected pyrosequencing provided more comprehensive information; additional prevalent species included M. chlorophenolicum, M. neglectum, M. gordonae, M. aemonae. Prevalence of the total Mycobacterium genus in the soil samples ranged from 2.3×107 to 2.7×108 gene targets g−1; slow growers prevalence from 2.9×105 to 1.2×107 cells g−1. Conclusions: This combined molecular approach enabled an unprecedented qualitative and quantitative assessment of EM across soil samples. Good concordance was found between methods and the bioinformatics analysis was validated by random resampling. Sequences from most pathogenic groups associated with slow growth were identified in extenso in all soils tested with a specific assay, allowing to unmask them from the Mycobacterium whole genus, in which, as minority members, they would have remained undetected

Comparative metagenomic, phylogenetic and physiological analyses of soil microbial communities across nitrogen gradients

Terrestrial ecosystems are receiving elevated inputs of nitrogen (N) from anthropogenic sources and understanding how these increases in N availability affect soil microbial communities is critical for predicting the associated effects on belowground ecosystems. We used a suite of approaches to analyze the structure and functional characteristics of soil microbial communities from replicated plots in two long-term N fertilization experiments located in contrasting systems. Pyrosequencing-based analyses of 16S rRNA genes revealed no significant effects of N fertilization on bacterial diversity, but significant effects on community composition at both sites; copiotrophic taxa (including members of the Proteobacteria and Bacteroidetes phyla) typically increased in relative abundance in the high N plots, with oligotrophic taxa (mainly Acidobacteria) exhibiting the opposite pattern. Consistent with the phylogenetic shifts under N fertilization, shotgun metagenomic sequencing revealed increases in the relative abundances of genes associated with DNA/RNA replication, electron transport and protein metabolism, increases that could be resolved even with the shallow shotgun metagenomic sequencing conducted here (average of 75 000 reads per sample). We also observed shifts in the catabolic capabilities of the communities across the N gradients that were significantly correlated with the phylogenetic and metagenomic responses, indicating possible linkages between the structure and functioning of soil microbial communities. Overall, our results suggest that N fertilization may, directly or indirectly, induce a shift in the predominant microbial life-history strategies, favoring a more active, copiotrophic microbial community, a pattern that parallels the often observed replacement of K-selected with r-selected plant species with elevated N

Characterization of the Oral Fungal Microbiome (Mycobiome) in Healthy Individuals

The oral microbiome–organisms residing in the oral cavity and their collective genome–are critical components of health and disease. The fungal component of the oral microbiota has not been characterized. In this study, we used a novel multitag pyrosequencing approach to characterize fungi present in the oral cavity of 20 healthy individuals, using the pan-fungal internal transcribed spacer (ITS) primers. Our results revealed the “basal” oral mycobiome profile of the enrolled individuals, and showed that across all the samples studied, the oral cavity contained 74 culturable and 11 non-culturable fungal genera. Among these genera, 39 were present in only one person, 16 genera were present in two participants, and 5 genera were present in three people, while 15 genera (including non-culturable organisms) were present in ≥4 (20%) participants. Candida species were the most frequent (isolated from 75% of participants), followed by Cladosporium (65%), Aureobasidium, Saccharomycetales (50% for both), Aspergillus (35%), Fusarium (30%), and Cryptococcus (20%). Four of these predominant genera are known to be pathogenic in humans. The low-abundance genera may represent environmental fungi present in the oral cavity and could simply be spores inhaled from the air or material ingested with food. Among the culturable genera, 61 were represented by one species each, while 13 genera comprised between 2 and 6 different species; the total number of species identified were 101. The number of species in the oral cavity of each individual ranged between 9 and 23. Principal component (PCO) analysis of the obtained data set followed by sample clustering and UniFrac analysis revealed that White males and Asian males clustered differently from each other, whereas both Asian and White females clustered together. This is the first study that identified the “basal mycobiome” of healthy individuals, and provides the basis for a detailed characterization of the oral mycobiome in health and disease

Dirichlet Multinomial Mixtures: Generative Models for Microbial Metagenomics

We introduce Dirichlet multinomial mixtures (DMM) for the probabilistic modelling of microbial metagenomics data. This data can be represented as a frequency matrix giving the number of times each taxa is observed in each sample. The samples have different size, and the matrix is sparse, as communities are diverse and skewed to rare taxa. Most methods used previously to classify or cluster samples have ignored these features. We describe each community by a vector of taxa probabilities. These vectors are generated from one of a finite number of Dirichlet mixture components each with different hyperparameters. Observed samples are generated through multinomial sampling. The mixture components cluster communities into distinct ‘metacommunities’, and, hence, determine envirotypes or enterotypes, groups of communities with a similar composition. The model can also deduce the impact of a treatment and be used for classification. We wrote software for the fitting of DMM models using the ‘evidence framework’ (http://code.google.com/p/microbedmm/). This includes the Laplace approximation of the model evidence. We applied the DMM model to human gut microbe genera frequencies from Obese and Lean twins. From the model evidence four clusters fit this data best. Two clusters were dominated by Bacteroides and were homogenous; two had a more variable community composition. We could not find a significant impact of body mass on community structure. However, Obese twins were more likely to derive from the high variance clusters. We propose that obesity is not associated with a distinct microbiota but increases the chance that an individual derives from a disturbed enterotype. This is an example of the ‘Anna Karenina principle (AKP)’ applied to microbial communities: disturbed states having many more configurations than undisturbed. We verify this by showing that in a study of inflammatory bowel disease (IBD) phenotypes, ileal Crohn's disease (ICD) is associated with a more variable community

16S rRNA Gene-based Analysis of Fecal Microbiota from Preterm Infants with and without Necrotizing Enterocolitis

Neonatal necrotizing enterocolitis (NEC) is an inflammatory intestinal disorder affecting preterm infants. Intestinal bacteria play a key role; however no causative pathogen has been identified. The purpose of this study was to determine if there are differences in microbial patterns which may be critical to the development of this disease. Fecal samples from twenty preterm infants, ten with NEC and ten matched controls (including four twin pairs) were obtained from patients in a single site Level III neonatal intensive care unit. Bacterial DNA from individual fecal samples were PCR amplified and subjected to terminal restriction fragment length polymorphism analysis and library sequencing of the 16S rRNA gene to characterize diversity and structure of the enteric microbiota. The distribution of samples from NEC patients distinctly clustered separately from controls. Intestinal bacterial colonization in all preterm infants was notable for low diversity. Patients with NEC had even less diversity, an increase in abundance of Gammaproteobacteria, a decrease in other bacteria species, and had received a higher mean number of previous days of antibiotics. Our results suggest that NEC is associated with severe lack of microbiota diversity which may accentuate the impact of single dominant microorganisms favored by empiric and wide-spread use of antibiotics

Bacterial and Archaea Community Present in the Pine Barrens Forest of Long Island, NY: Unusually High Percentage of Ammonia Oxidizing Bacteria

Of the few preserved areas in the northeast of United States, the soil in the Pine Barrens Forests presents a harsh environment for the microorganisms to grow and survive. In the current study we report the use of clustering methods to scientifically select the sampling locations that would represent the entire forest and also report the microbial diversity present in various horizons of the soil. Sixty six sampling locations were selected across the forest and soils were collected from three horizons (sampling depths). The three horizons were 0–10 cm (Horizon O); 11–25 cm (Horizon A) and 26–40 cm (Horizon B). Based on the total microbial substrate utilization pattern and K-means clustering analysis, the soil in the Pine Barrens Forest can be classified into four distinct clusters at each of the three horizons. One soil sample from each of the four clusters were selected and archaeal and bacterial populations within the soil studied using pyrosequencing method. The results show the microbial communities present in each of these clusters are different. Within the microbial communities present, microorganisms involved in nitrogen cycle occupy a major fraction of microbial community in the soil. High level of diversity was observed for nitrogen fixing bacteria. In contrast, Nitrosovibrio and Nitrosocaldus spp are the single bacterial and archaeal population respectively carrying out ammonia oxidation in the soil

Environmental Shaping of Sponge Associated Archaeal Communities

Archaea are ubiquitous symbionts of marine sponges but their ecological roles and the influence of environmental factors on these associations are still poorly understood.We compared the diversity and composition of archaea associated with seawater and with the sponges Hymeniacidon heliophila, Paraleucilla magna and Petromica citrina in two distinct environments: Guanabara Bay, a highly impacted estuary in Rio de Janeiro, Brazil, and the nearby Cagarras Archipelago. For this we used metagenomic analyses of 16S rRNA and ammonia monooxygenase (amoA) gene libraries. Hymeniacidon heliophila was more abundant inside the bay, while P. magna was more abundant outside and P. citrina was only recorded at the Cagarras Archipelago. Principal Component Analysis plots (PCA) generated using pairwise unweighted UniFrac distances showed that the archaeal community structure of inner bay seawater and sponges was different from that of coastal Cagarras Archipelago. Rarefaction analyses showed that inner bay archaeaoplankton were more diverse than those from the Cagarras Archipelago. Only members of Crenarchaeota were found in sponge libraries, while in seawater both Crenarchaeota and Euryarchaeota were observed. Although most amoA archaeal genes detected in this study seem to be novel, some clones were affiliated to known ammonia oxidizers such as Nitrosopumilus maritimus and Cenarchaeum symbiosum.The composition and diversity of archaeal communities associated with pollution-tolerant sponge species can change in a range of few kilometers, probably influenced by eutrophication. The presence of archaeal amoA genes in Porifera suggests that Archaea are involved in the nitrogen cycle within the sponge holobiont, possibly increasing its resistance to anthropogenic impacts. The higher diversity of Crenarchaeota in the polluted area suggests that some marine sponges are able to change the composition of their associated archaeal communities, thereby improving their fitness in impacted environments

CORE: A Phylogenetically-Curated 16S rDNA Database of the Core Oral Microbiome

Comparing bacterial 16S rDNA sequences to GenBank and other large public databases via BLAST often provides results of little use for identification and taxonomic assignment of the organisms of interest. The human microbiome, and in particular the oral microbiome, includes many taxa, and accurate identification of sequence data is essential for studies of these communities. For this purpose, a phylogenetically curated 16S rDNA database of the core oral microbiome, CORE, was developed. The goal was to include a comprehensive and minimally redundant representation of the bacteria that regularly reside in the human oral cavity with computationally robust classification at the level of species and genus. Clades of cultivated and uncultivated taxa were formed based on sequence analyses using multiple criteria, including maximum-likelihood-based topology and bootstrap support, genetic distance, and previous naming. A number of classification inconsistencies for previously named species, especially at the level of genus, were resolved. The performance of the CORE database for identifying clinical sequences was compared to that of three publicly available databases, GenBank nr/nt, RDP and HOMD, using a set of sequencing reads that had not been used in creation of the database. CORE offered improved performance compared to other public databases for identification of human oral bacterial 16S sequences by a number of criteria. In addition, the CORE database and phylogenetic tree provide a framework for measures of community divergence, and the focused size of the database offers advantages of efficiency for BLAST searching of large datasets. The CORE database is available as a searchable interface and for download at http://microbiome.osu.edu

Spatial distribution and abundances of ammonia-oxidizing archaea (AOA) and ammonia-oxidizing bacteria (AOB) in mangrove sediments

We investigated the diversity, spatial distribution, and abundances of ammonia-oxidizing archaea (AOA) and ammonia-oxidizing bacteria (AOB) in sediment samples of different depths collected from a transect with different distances to mangrove forest in the territories of Hong Kong. Both the archaeal and bacterial amoA genes (encoding ammonia monooxygenase subunit A) from all samples supported distinct phylogenetic groups, indicating the presences of niche-specific AOA and AOB in mangrove sediments. The higher AOB abundances than AOA in mangrove sediments, especially in the vicinity of the mangrove trees, might indicate the more important role of AOB on nitrification. The spatial distribution showed that AOA had higher diversity and abundance in the surface layer sediments near the mangrove trees (0 and 10 m) but lower away from the mangrove trees (1,000 m), and communities of AOA could be clustered into surface and bottom sediment layer groups. In contrast, AOB showed a reverse distributed pattern, and its communities were grouped by the distances between sites and mangrove trees, indicating mangrove trees might have different influences on AOA and AOB community structures. Furthermore, the strong correlations among archaeal and bacterial amoA gene abundances and their ratio with NH4+, salinity, and pH of sediments indicated that these environmental factors have strong influences on AOA and AOB distributions in mangrove sediments. In addition, AOA diversity and abundances were significantly correlated with hzo gene abundances, which encodes the key enzyme for transformation of hydrazine into N2 in anaerobic ammonium-oxidizing (anammox) bacteria, indicating AOA and anammox bacteria may interact with each other or they are influenced by the same controlling factors, such as NH4+. The results provide a better understanding on using mangrove wetlands as biological treatment systems for removal of nutrients