The variability of the 16S rRNA gene in bacterial genomes and its consequences for bacterial community analyses

Microbial Population Differentials between Mucosal and Submucosal Intestinal Tissues in Advanced Crohn's Disease of the Ileum

Since Crohn's disease is a transmural disease, we hypothesized that examination of deep submucosal tissues directly involved in the inflammatory disease process may provide unique insights into bacterial populations transgressing intestinal barriers and bacterial populations more representative of the causes and agents of the disease. We performed deep 16s microbiota sequencing on isolated ilea mucosal and submucosal tissues on 20 patients with Crohn's disease and 15 non-inflammatory bowel disease controls with a depth of coverage averaging 81,500 sequences in each of the 70 DNA samples yielding an overall resolution down to 0.0001% of the bacterial population. Of the 4,802,328 total sequences generated, 98.9% or 4,749,183 sequences aligned with the Kingdom Bacteria that clustered into 8545 unique sequences with <3% divergence or operational taxonomic units enabling the identification of 401 genera and 698 tentative bacterial species. There were significant differences in all taxonomic levels between the submucosal microbiota in Crohn's disease compared to controls, including organisms of the Order Desulfovibrionales that were present within the submucosal tissues of most Crohn's disease patients but absent in the control group. A variety of organisms of the Phylum Firmicutes were increased in the subjacent submucosa as compared to the parallel mucosal tissue including Ruminococcus spp., Oscillospira spp., Pseudobutyrivibrio spp., and Tumebacillus spp. In addition, Propionibacterium spp. and Cloacibacterium spp. were increased as well as large increases in Proteobacteria including Parasutterella spp. and Methylobacterium spp. This is the first study to examine the microbial populations within submucosal tissues of patients with Crohn's disease and to compare microbial communities found deep within the submucosal tissues with those present on mucosal surfaces. Our data demonstrate the existence of a distinct submucosal microbiome and ecosystem that is not well reflected in the mucosa and/or downstream fecal material.

Bacterial diversity differences along an epigenic cave stream reveal evidence of community dynamics, succession, and stability

Unchanging physicochemical conditions and nutrient sources over long periods of time in cave and karst subsurface habitats, particularly aquifers, can support stable ecosystems, termed autochthonous microbial endokarst communities (AMEC). AMEC existence is unknown for other karst settings, such as epigenic cave streams. Conceptually, AMEC should not form in streams due to faster turnover rates and seasonal disturbances that have the capacity to transport large quantities of water and sediment and to change allochthonous nutrient and organic matter sources. Our goal was to investigate whether AMEC could form and persist in hydrologically active, epigenic cave streams. We analyzed bacterial diversity from cave water, sediments, and artificial substrates (Bio-Traps®) placed in the cave at upstream and downstream locations. Distinct communities existed for the water, sediments, and Bio-Trap® samplers. Throughout the study period, a subset of community members persisted in the water, regardless of hydrological disturbances. Stable habitat conditions based on flow regimes resulted in more than one contemporaneous, stable community throughout the epigenic cave stream. However, evidence for AMEC was insufficient for the cave water or sediments. Community succession, specifically as predictable exogenous heterotrophic microbial community succession, was evident from decreases in community richness from the Bio-Traps®, a peak in Bio-Trap® community biomass, and from changes in the composition of Bio-Trap® communities. The planktonic community was compositionally similar to Bio-Trap® initial colonizers, but the downstream Bio-Trap® community became more similar to the sediment community at the same location. These results can help in understanding the diversity of planktonic and attached microbial communities from karst, as well as microbial community dynamics, stability, and succession during disturbance or contamination responses over time.

Differential Decay of Wastewater Bacteria and Change of Microbial Communities in Beach Sand and Seawater Microcosms

Laboratory microcosm experiments were conducted to determine the decay kinetics of wastewater bacteria and the change of microbial communities in beach sand and seawater. Cultivation-based methods showed that common fecal indicator bacteria (FIBs; Escherichia coli, enterococci, and Clostridium perfringens) exhibited biphasic decay patterns in all microcosms. Enterococci and C. perfringens, but not E. coli, showed significantly smaller decay rates in beach sand than in seawater. Cultivation-independent qPCR quantification of 16S rRNA gene also showed significantly slower decrease of total bacterial densities in beach sand than in seawater. Microbial community analysis by next-generation sequencing (NGS) further illustrated that the decreasing relative abundance of wastewater bacteria was contrasted by the increase in indigenous beach sand and seawater microbiota, and the overall microbial community dynamics corresponded well with the decay of individual FIB populations. In summary, the differential decay of wastewater bacteria in beach sand and in seawater provides a kinetic explanation to the often-observed higher abundance of FIBs in beach sand, and the NGS-based microbial community analysis can provide valuable insights to understanding the fate of wastewater bacteria in the context of indigenous microbial communities in natural environments.

Analysis of the intestinal bacterial microbiota in maize- or sorghum-fed broiler chickens using real-time PCR

1. An experiment was conducted to study the effect of two different diets on zootechnical performance and the major bacterial groups in association with the host mucosa and dispersed in the lumen contents of the small intestine of broiler chickens. 2. The two experimental diets were maize or sorghum-based. In addition to the total bacteria, bacterial groups belonging to the Enterobacteriaceae (Enterococcus and Lactobacillus) were quantified by real-time PCR. 3. There were no differences in body weight gain and feed intake, but feed conversion ratio increased for sorghum-fed broilers at 21 and 42 d of age. 4. The Enterococcus group decreased in all gut segments from 7 to 42 d, while the Lactobacillus group increased in both ecosystems. In the ileal mucosa, the enterobacterial counts decreased from 7 to 42 d in the maize-based diet, but remained stable in the sorghum-based diet. 5. The results shed light on the spatial and temporal distribution of bacterial groups that play important physiological roles in the small intestine of chickens. Specifically, the increased Enterobacteria population in the ileum is consistent with the relatively poor feed conversion in sorghum-fed broilers.

Shifts in soil chemical properties and bacterial communities responding to biotransformed dry olive residue used as organic amendment

Dry olive residue (DOR) is a waste product derived from olive oil extraction and has been proposed as an organic amendment. However, it has been demonstrated that a pre-treatment, such as its transformation by saprophytic fungi, is required before DOR soil application. A greenhouse experiment was designed where 0 and 50 g kg(-1) of raw DOR (DOR), Coriolopsis floccosa-transformed DOR (CORDOR) and Fusarium oxysporum-transformed DOR (FUSDOR) were added to soil. Analyses of the soil chemical properties as well as the structure and relative abundance of bacterial and actinobacterial communities were conducted after 0, 30 and 60 days following amendment. The different amendments produced a slight decrease in soil pH and significant increases in carbon fractions, C/N ratios, phenols and K, with these increases being more significant after DOR application. Quantitative PCR assays of the 16S rRNA gene and PLFA analyses showed that all amendments favoured bacterial growth at 30 and 60 days, although actinobacterial proliferation was more evident after CORDOR and FUSDOR application at 60 days. Bacterial and actinobacterial DGGE multivariate analyses showed that the amendments produced structural changes in both communities, especially after 60 days of amendment. PLFA data analysis identified changes in soil microbial communities according to the amendment considered, with FUSDOR and CORDOR being less disruptive than DOR. Finally, integrated analysis of all data monitored in the present study enabled us to conclude that the greatest impact on soil properties was caused by DOR at 30 days and that soil showed some degree of resilience after this time.

Comparing K-mer based methods for improved classification of 16S sequences

Background

The need for precise and stable taxonomic classification is highly relevant in modern microbiology. Parallel to the explosion in the amount of sequence data accessible, there has also been a shift in focus for classification methods. Previously, alignment-based methods were the most applicable tools. Now, methods based on counting K-mers by sliding windows are the most interesting classification approach with respect to both speed and accuracy. Here, we present a systematic comparison on five different K-mer based classification methods for the 16S rRNA gene. The methods differ from each other both in data usage and modelling strategies. We have based our study on the commonly known and well-used naïve Bayes classifier from the RDP project, and four other methods were implemented and tested on two different data sets, on full-length sequences as well as fragments of typical read-length.

Results

The difference in classification error obtained by the methods seemed to be small, but they were stable and for both data sets tested. The Preprocessed nearest-neighbour (PLSNN) method performed best for full-length 16S rRNA sequences, significantly better than the naïve Bayes RDP method. On fragmented sequences the naïve Bayes Multinomial method performed best, significantly better than all other methods. For both data sets explored, and on both full-length and fragmented sequences, all the five methods reached an error-plateau.

Conclusions

We conclude that no K-mer based method is universally best for classifying both full-length sequences and fragments (reads). All methods approach an error plateau indicating improved training data is needed to improve classification from here. Classification errors occur most frequent for genera with few sequences present. For improving the taxonomy and testing new classification methods, the need for a better and more universal and robust training data set is crucial.

Inhibition of the endosymbiont "Candidatus Midichloria mitochondrii" during 16S rRNA gene profiling reveals potential pathogens in Ixodes ticks from Australia

Background

The Australian paralysis tick (Ixodes holocyclus) is of significant medical and veterinary importance as a cause of dermatological and neurological disease, yet there is currently limited information about the bacterial communities harboured by these ticks and the risk of infectious disease transmission to humans and domestic animals. Ongoing controversy about the presence of Borrelia burgdorferi sensu lato (the aetiological agent of Lyme disease) in Australia increases the need to accurately identify and characterise bacteria harboured by I. holocyclus ticks.

Methods

Universal PCR primers were used to amplify the V1-2 hyper-variable region of bacterial 16S rRNA genes present in DNA samples from I. holocyclus and I. ricinus ticks, collected in Australia and Germany respectively. The 16S amplicons were purified, sequenced on the Ion Torrent platform, and analysed in USEARCH, QIIME, and BLAST to assign genus and species-level taxonomy. Initial analysis of I. holocyclus and I. ricinus identified that > 95 % of the 16S sequences recovered belonged to the tick intracellular endosymbiont “Candidatus Midichloria mitochondrii” (CMM). A CMM-specific blocking primer was designed that decreased CMM sequences by approximately 96 % in both tick species and significantly increased the total detectable bacterial diversity, allowing identification of medically important bacterial pathogens that were previously masked by CMM.

Results

Borrelia burgdorferi sensu lato was identified in German I. ricinus, but not in Australian I. holocyclus ticks. However, bacteria of medical significance were detected in I. holocyclus ticks, including a Borrelia relapsing fever group sp., Bartonella henselae, novel “Candidatus Neoehrlichia” spp., Clostridium histolyticum, Rickettsia spp., and Leptospira inadai.

Conclusions

Abundant bacterial endosymbionts, such as CMM, limit the effectiveness of next-generation 16S bacterial community profiling in arthropods by masking less abundant bacteria, including pathogens. Specific blocking primers that inhibit endosymbiont 16S amplification during PCR are an effective way of reducing this limitation. Here, this strategy provided the first evidence of a relapsing fever Borrelia sp. and of novel “Candidatus Neoehrlichia” spp. in Australia. Our results raise new questions about tick-borne pathogens in I. holocyclus ticks.

Electronic supplementary material

The online version of this article (doi:10.1186/s13071-015-0958-3) contains supplementary material, which is available to authorized users.

Bandage: interactive visualization of de novo genome assemblies

Summary: Although de novo assembly graphs contain assembled contigs (nodes), the connections between those contigs (edges) are difficult for users to access. Bandage (a Bioinformatics Application for Navigating De novo Assembly Graphs Easily) is a tool for visualizing assembly graphs with connections. Users can zoom in to specific areas of the graph and interact with it by moving nodes, adding labels, changing colors and extracting sequences. BLAST searches can be performed within the Bandage graphical user interface and the hits are displayed as highlights in the graph. By displaying connections between contigs, Bandage presents new possibilities for analyzing de novo assemblies that are not possible through investigation of contigs alone.

Availability and implementation: Source code and binaries are freely available at https://github.com/rrwick/Bandage. Bandage is implemented in C++ and supported on Linux, OS X and Windows. A full feature list and screenshots are available at http://rrwick.github.io/Bandage.

Contact:rrwick@gmail.com

Supplementary information:Supplementary data are available at Bioinformatics online.

Phasing of single DNA molecules by massively parallel barcoding

High-throughput sequencing platforms mainly produce short-read data, resulting in a loss of phasing information for many of the genetic variants analysed. For certain applications, it is vital to know which variant alleles are connected to each individual DNA molecule. Here we demonstrate a method for massively parallel barcoding and phasing of single DNA molecules. First, a primer library with millions of uniquely barcoded beads is generated. When compartmentalized with single DNA molecules, the beads can be used to amplify and tag any target sequences of interest, enabling coupling of the biological information from multiple loci. We apply the assay to bacterial 16S sequencing and up to 94% of the hypothesized phasing events are shown to originate from single molecules. The method enables use of widely available short-read-sequencing platforms to study long single molecules within a complex sample, without losing phase information.

Integrated community profiling indicates long-term temporal stability of the predominant faecal microbiota in captive cheetahs

Understanding the symbiotic relationship between gut microbes and their animal host requires characterization of the core microbiota across populations and in time. Especially in captive populations of endangered wildlife species such as the cheetah (Acinonyx jubatus), this knowledge is a key element to enhance feeding strategies and reduce gastrointestinal disorders. In order to investigate the temporal stability of the intestinal microbiota in cheetahs under human care, we conducted a longitudinal study over a 3-year period with bimonthly faecal sampling of 5 cheetahs housed in two European zoos. For this purpose, an integrated 16S rRNA DGGE-clone library approach was used in combination with a series of real-time PCR assays. Our findings disclosed a stable faecal microbiota, beyond intestinal community variations that were detected between zoo sample sets or between animals. The core of this microbiota was dominated by members of Clostridium clusters I, XI and XIVa, with mean concentrations ranging from 7.5-9.2 log10 CFU/g faeces and with significant positive correlations between these clusters (P<0.05), and by Lactobacillaceae. Moving window analysis of DGGE profiles revealed 23.3-25.6% change between consecutive samples for four of the cheetahs. The fifth animal in the study suffered from intermediate episodes of vomiting and diarrhea during the monitoring period and exhibited remarkably more change (39.4%). This observation may reflect the temporary impact of perturbations such as the animal’s compromised health, antibiotic administration or a combination thereof, which temporarily altered the relative proportions of Clostridium clusters I and XIVa. In conclusion, this first long-term monitoring study of the faecal microbiota in feline strict carnivores not only reveals a remarkable compositional stability of this ecosystem, but also shows a qualitative and quantitative similarity in a defined set of faecal bacterial lineages across the five animals under study that may typify the core phylogenetic microbiome of cheetahs.

Profiling bacterial diversity in a limestone cave of the western Loess Plateau of China

Bacteria and archaea sustain subsurface cave ecosystems by dominating primary production and fueling biogeochemical cyclings, despite the permanent darkness and shortage of nutrients. However, the heterogeneity and underlying mechanism of microbial diversity in caves, in particular those well connect to surface environment are largely unexplored. In this study, we examined the bacterial abundance and composition in Jinjia Cave, a small and shallow limestone cave located on the western Loess Plateau of China, by enumerating and pyrosequencing small subunit rRNA genes. The results clearly reveal the contrasting bacterial community compositions in relation to cave habitat types, i.e., rock wall deposit, aquatic sediment, and sinkhole soil, which are differentially connected to the surface environment. The deposits on the cave walls were dominated by putative cave-specific bacterial lineages within the γ-Proteobacteria or Actinobacteria that are routinely found on cave rocks around the world. In addition, sequence identity with known functional groups suggests enrichments of chemolithotrophic bacteria potentially involved in autotrophic C fixation and inorganic N transformation on rock surfaces. By contrast, bacterial communities in aquatic sediments were more closely related to those in the overlying soils. This is consistent with the similarity in elemental composition between the cave sediment and the overlying soil, implicating the influence of mineral chemistry on cave microhabitat and bacterial composition. These findings provide compelling molecular evidence of the bacterial community heterogeneity in an East Asian cave, which might be controlled by both subsurface and surface environments.

Metagenomic analysis of the impact of nitrofurantoin treatment on the human faecal microbiota

OBJECTIVES

The objective was to study changes in the faecal microbiota of patients with uncomplicated urinary tract infections (UTIs) treated with nitrofurantoin and of non-treated healthy controls using 16S rRNA analysis.

METHODS

Serial stool samples were collected from patients receiving nitrofurantoin treatment at different timepoints [before treatment (day 1; T1), within 48 h of end of treatment (days 5-15; T2) and 28 days after treatment (days 31-43; T3)], as well as from healthy controls. Direct DNA extraction (PowerSoil DNA Isolation Kit, MoBio Laboratories, Carlsbad, CA, USA) from stool samples was followed by pyrosequencing (454 GS FLX Titanium) of the V3-V5 region of the bacterial 16S rRNA gene.

RESULTS

Among UTI patients, mean proportions of the Actinobacteria phylum increased by 19.6% in the first follow-up sample (T2) in comparison with the pretreatment baseline stool sample (T1) (P = 0.026). However, proportions of Actinobacteria reversed to 'normal' pre-antibiotic levels, with a mean difference of 1.0% compared with baseline proportions, in the second follow-up sample (T3). The increase in Actinobacteria was specifically due to an increase in the Bifidobacteriaceae family (Bifidobacterium genus), which constituted 81.0% (95% CI ±7.4%) of this phylum.

CONCLUSIONS

No significant impact was observed other than a temporary increase in the beneficial Bifidobacterium genus following nitrofurantoin treatment, which supports its reintroduction into clinical use.

Quantitative profiling of colorectal cancer-associated bacteria reveals associations between fusobacterium spp., enterotoxigenic Bacteroides fragilis (ETBF) and clinicopathological features of colorectal cancer

Various studies have presented clinical or in vitro evidence linking bacteria to colorectal cancer, but these bacteria have not previously been concurrently quantified by qPCR in a single cohort. We quantify these bacteria (Fusobacterium spp., Streptococcus gallolyticus, Enterococcus faecalis, Enterotoxigenic Bacteroides fragilis (ETBF), Enteropathogenic Escherichia coli (EPEC), and afaC- or pks-positive E. coli) in paired tumour and normal tissue samples from 55 colorectal cancer patients. We further investigate the relationship between a) the presence and b) the level of colonisation of each bacterial species with site and stage of disease, age, gender, ethnicity and MSI-status. With the exception of S. gallolyticus, we detected all bacteria profiled here in both tumour and normal samples at varying frequencies. ETBF (FDR = 0.001 and 0.002 for normal and tumour samples) and afaC-positive E. coli (FDR = 0.03, normal samples) were significantly enriched in the colon compared to the rectum. ETBF (FDR = 0.04 and 0.002 for normal and tumour samples, respectively) and Fusobacterium spp. (FDR = 0.03 tumour samples) levels were significantly higher in late stage (III/IV) colorectal cancers. Fusobacterium was by far the most common bacteria detected, occurring in 82% and 81% of paired tumour and normal samples. Fusobacterium was also the only bacterium that was significantly higher in tumour compared to normal samples (p = 6e-5). We also identified significant associations between high-level colonisation by Fusobacterium and MSI-H (FDR = 0.05), age (FDR = 0.03) or pks-positive E. coli (FDR = 0.01). Furthermore, we exclusively identified atypical EPEC in our cohort, which has not been previously reported in association with colorectal cancer. By quantifying colorectal cancer-associated bacteria across a single cohort, we uncovered inter- and intra-individual patterns of colonization not previously recognized, as well as important associations with clinicopathological features, especially in the case of Fusobacterium and ETBF.

International interlaboratory study comparing single organism 16S rRNA gene sequencing data: Beyond consensus sequence comparisons

This study presents the results from an interlaboratory sequencing study for which we developed a novel high-resolution method for comparing data from different sequencing platforms for a multi-copy, paralogous gene. The combination of PCR amplification and 16S ribosomal RNA gene (16S rRNA) sequencing has revolutionized bacteriology by enabling rapid identification, frequently without the need for culture. To assess variability between laboratories in sequencing 16S rRNA, six laboratories sequenced the gene encoding the 16S rRNA from Escherichia coli O157:H7 strain EDL933 and Listeria monocytogenes serovar 4b strain NCTC11994. Participants performed sequencing methods and protocols available in their laboratories: Sanger sequencing, Roche 454 pyrosequencing(®), or Ion Torrent PGM(®). The sequencing data were evaluated on three levels: (1) identity of biologically conserved position, (2) ratio of 16S rRNA gene copies featuring identified variants, and (3) the collection of variant combinations in a set of 16S rRNA gene copies. The same set of biologically conserved positions was identified for each sequencing method. Analytical methods using Bayesian and maximum likelihood statistics were developed to estimate variant copy ratios, which describe the ratio of nucleotides at each identified biologically variable position, as well as the likely set of variant combinations present in 16S rRNA gene copies. Our results indicate that estimated variant copy ratios at biologically variable positions were only reproducible for high throughput sequencing methods. Furthermore, the likely variant combination set was only reproducible with increased sequencing depth and longer read lengths. We also demonstrate novel methods for evaluating variable positions when comparing multi-copy gene sequence data from multiple laboratories generated using multiple sequencing technologies.

Oligotyping reveals stronger relationship of organic soil bacterial community structure with N-amendments and soil chemistry in comparison to that of mineral soil at Harvard Forest, MA, USA

The impact of chronic nitrogen amendments on bacterial communities was evaluated at Harvard Forest, Petersham, MA, USA. Thirty soil samples (3 treatments × 2 soil horizons × 5 subplots) were collected in 2009 from untreated (control), low nitrogen-amended (LN; 50 kg NH4NO3 ha(-1) yr(-1)) and high nitrogen-amended (HN; 150 kg NH4NO3 ha(-1) yr(-1)) plots. PCR-amplified partial 16S rRNA gene sequences made from soil DNA were subjected to pyrosequencing (Turlapati et al., 2013) and analyses using oligotyping. The parameters M (the minimum count of the most abundant unique sequence in an oligotype) and s (the minimum number of samples in which an oligotype is expected to be present) had to be optimized for forest soils because of high diversity and the presence of rare organisms. Comparative analyses of the pyrosequencing data by oligotyping and operational taxonomic unit clustering tools indicated that the former yields more refined units of taxonomy with sequence similarity of ≥99.5%. Sequences affiliated with four new phyla and 73 genera were identified in the present study as compared to 27 genera reported earlier from the same data (Turlapati et al., 2013). Significant rearrangements in the bacterial community structure were observed with N-amendments revealing the presence of additional genera in N-amended plots with the absence of some that were present in the control plots. Permutational MANOVA analyses indicated significant variation associated with soil horizon and N treatment for a majority of the phyla. In most cases soil horizon partitioned more variation relative to treatment and treatment effects were more evident for the organic (Org) horizon. Mantel test results for Org soil showed significant positive correlations between bacterial communities and most soil parameters including NH4 and NO3. In mineral soil, correlations were seen only with pH, NH4, and NO3. Regardless of the pipeline used, a major hindrance for such a study remains to be the lack of reference databases for forest soils.

The hidden 'mycobacteriome' of the human healthy oral cavity and upper respiratory tract

The incidence of opportunistic non-tuberculous mycobacteria (NTM) infections has increased considerably in the past decades causing an array of infections, including respiratory and soft-tissue infections. NTM are ubiquitous and can be found in numerous environments, including households and water plants. However, NTM have not been reported to be associated with the healthy human oral microbiome. Since the oral cavity and upper respiratory track are the main ports of entry of microorganisms into the human body, elucidating NTM diversity and prevalence will assist in the assessment of the potential risks of infection elicited by these opportunistic pathogens. Here, we report the identification of a ‘non-tuberculous mycobacteriome’ in healthy individuals. We employed a modified DNA extraction procedure in conjunction with mycobacterial-specific primers to screen niches in the oral cavity (buccal mucosa and dental plaque) and upper respiratory tract (nostrils and oropharynx) of 10 healthy subjects. A total of 50 prevalent operational taxonomic units sequenced on MiSeq (Illumina) using 16S rRNA V3–V4 region were detected across all screened niches, showing the presence of diverse NTM communities. NTM DNA was detected in the nostrils of all 10 subjects, in buccal mucosa of 8 subjects, in the oropharynx of 7 subjects, and in the dental plaques of 5 subjects. Results from quantitative PCR showed each individual harbored 10³–10⁴ predicted NTM per each screened niche. The modification of standard DNA isolation methods to increase sensitivity toward mycobacterial species represents an important step to advance the knowledge of the oral as well as the overall human microbiome. These findings clearly reveal for the first time that healthy individuals harbor a ‘non-tuberculous mycobacteriome’ in their oral cavity and upper respiratory tract and may have important implications in our understanding of infections caused by NTM.

Long-term forest soil warming alters microbial communities in temperate forest soils

Soil microbes are major drivers of soil carbon cycling, yet we lack an understanding of how climate warming will affect microbial communities. Three ongoing field studies at the Harvard Forest Long-term Ecological Research (LTER) site (Petersham, MA) have warmed soils 5°C above ambient temperatures for 5, 8, and 20 years. We used this chronosequence to test the hypothesis that soil microbial communities have changed in response to chronic warming. Bacterial community composition was studied using Illumina sequencing of the 16S ribosomal RNA gene, and bacterial and fungal abundance were assessed using quantitative PCR. Only the 20-year warmed site exhibited significant change in bacterial community structure in the organic soil horizon, with no significant changes in the mineral soil. The dominant taxa, abundant at 0.1% or greater, represented 0.3% of the richness but nearly 50% of the observations (sequences). Individual members of the Actinobacteria, Alphaproteobacteria and Acidobacteria showed strong warming responses, with one Actinomycete decreasing from 4.5 to 1% relative abundance with warming. Ribosomal RNA copy number can obfuscate community profiles, but is also correlated with maximum growth rate or trophic strategy among bacteria. Ribosomal RNA copy number correction did not affect community profiles, but rRNA copy number was significantly decreased in warming plots compared to controls. Increased bacterial evenness, shifting beta diversity, decreased fungal abundance and increased abundance of bacteria with low rRNA operon copy number, including Alphaproteobacteria and Acidobacteria, together suggest that more or alternative niche space is being created over the course of long-term warming.

Emergence shapes the structure of the seed microbiota

Seeds carry complex microbial communities, which may exert beneficial or deleterious effects on plant growth and plant health. To date, the composition of microbial communities associated with seeds has been explored mainly through culture-based diversity studies and therefore remains largely unknown. In this work, we analyzed the structures of the seed microbiotas of different plants from the family Brassicaceae and their dynamics during germination and emergence through sequencing of three molecular markers: the ITS1 region of the fungal internal transcribed spacer, the V4 region of 16S rRNA gene, and a species-specific bacterial marker based on a fragment of gyrB. Sequence analyses revealed important variations in microbial community composition between seed samples. Moreover, we found that emergence strongly influences the structure of the microbiota, with a marked reduction of bacterial and fungal diversity. This shift in the microbial community composition is mostly due to an increase in the relative abundance of some bacterial and fungal taxa possessing fast-growing abilities. Altogether, our results provide an estimation of the role of the seed as a source of inoculum for the seedling, which is crucial for practical applications in developing new strategies of inoculation for disease prevention.

High-resolution melt analysis for rapid comparison of bacterial community compositions

In the study of bacterial community composition, 16S rRNA gene amplicon sequencing is today among the preferred methods of analysis. The cost of nucleotide sequence analysis, including requisite computational and bioinformatic steps, however, takes up a large part of many research budgets. High-resolution melt (HRM) analysis is the study of the melt behavior of specific PCR products. Here we describe a novel high-throughput approach in which we used HRM analysis targeting the 16S rRNA gene to rapidly screen multiple complex samples for differences in bacterial community composition. We hypothesized that HRM analysis of amplified 16S rRNA genes from a soil ecosystem could be used as a screening tool to identify changes in bacterial community structure. This hypothesis was tested using a soil microcosm setup exposed to a total of six treatments representing different combinations of pesticide and fertilization treatments. The HRM analysis identified a shift in the bacterial community composition in two of the treatments, both including the soil fumigant Basamid GR. These results were confirmed with both denaturing gradient gel electrophoresis (DGGE) analysis and 454-based 16S rRNA gene amplicon sequencing. HRM analysis was shown to be a fast, high-throughput technique that can serve as an effective alternative to gel-based screening methods to monitor microbial community composition.

The role of the gut microbiome in the healthy adult status

The gut microbiome, which hosts up to 1000 bacterial species that encode about 5 million genes, perform many of the functions required for host physiology and survival. Consequently, it is also known as "our forgotten organ". The recent development of next-generation sequencing technologies has greatly improved metagenomic research. In particular, it has increased our knowledge about the microbiome and its mutually beneficial relationships with the human host. Microbial colonization begins immediately at birth. Although influenced by a variety of stimuli, namely, diet, physical activity, travel, illness, hormonal cycles and therapies, the microbiome is practically stable in healthy adults. This suggests that the microbiome plays a role in the maintenance of a healthy state in adulthood. Quantitative and qualitative alterations in the composition of the gut microbiome could lead to pathological dysbiosis, and have been related to an increasing number of intestinal and extra-intestinal diseases. With the increase in knowledge about gut microbiome functions, it is becoming increasingly more possible to develop novel diagnostic, prognostic and, most important, therapeutic strategies based on microbiome manipulation.

Effects of photoperiod on nutrient removal, biomass production, and algal-bacterial population dynamics in lab-scale photobioreactors treating municipal wastewater

Effects of photoperiod were investigated in lab-scale photobioreactors containing algal-bacterial consortia to reduce organic nutrients from municipal wastewater. Under three photoperiod conditions (12 h:12 h, 36 h:12 h, and 60 h:12 h dark–light cycles), nutrient removals and biomass productions were measured along with monitoring microbial population dynamics. After a batch operation for 12 days, 59–80% carbon, 35–88% nitrogen, and 43–89% phosphorus were removed from influents, respectively. In this study, carbon removal was related positively to the length of dark cycles, while nitrogen and phosphorus removals inversely. On the contrast, the highest microbial biomass in terms of chlorophyll a, dry cell weight, and algal/bacterial rRNA gene markers was produced under the 12 h:12 h dark–light cycle among the three photoperiods. The results showed 1) simultaneous growths between algae and bacteria in the microbial consortia and 2) efficient nitrogen and phosphorus removals along with high microbial biomass production under prolonged light conditions. Statistical analyses indicated that carbon removal was significantly related to the ratio of bacteria to algae in the microbial consortia along with prolonged dark conditions (p < 0.05). In addition, the ratio of nitrogen removal to phosphorus removal decreased significantly under prolonged dark conditions (p < 0.001). These results indicated that the photoperiod condition has remarkable impacts on adjusting nutrient removal, producing microbial biomass, and altering algal-bacterial population dynamics. Therefore, the control of photoperiod was suggested as an important operating parameter in the algal wastewater treatment.

Global biogeography of Alnus-associated Frankia actinobacteria

Macroecological patterns of microbes have received relatively little attention until recently. This study aimed to disentangle the determinants of the global biogeographic community of Alnus-associated actinobacteria belonging to the Frankia alni complex. By determining a global sequence similarity threshold for the nitrogenase reductase (nifH) gene, we separated Frankia into operational taxonomic units (OTUs) and tested the relative effects of Alnus phylogeny, geographic relatedness, and climatic and edaphic variables on community composition at the global scale. Based on the optimal nifH gene sequence similarity threshold of 99.3%, we distinguished 43 Frankia OTUs from root systems of 22 Alnus species on four continents. Host phylogeny was the main determinant of Frankia OTU-based community composition, but there was no effect on the phylogenetic structure of Frankia. Biogeographic analyses revealed the strongest cross-continental links over the Beringian land bridge. Despite the facultative symbiotic nature of Frankia, phylogenetic relations among Alnus species play a prominent role in structuring root-associated Frankia communities and their biogeographic patterns. Our results suggest that Alnus species exert strong phylogenetically determined selection pressure on compatible Actinobacteria.

rrnDB: improved tools for interpreting rRNA gene abundance in bacteria and archaea and a new foundation for future development

Microbiologists utilize ribosomal RNA genes as molecular markers of taxonomy in surveys of microbial communities. rRNA genes are often co-located as part of an rrn operon, and multiple copies of this operon are present in genomes across the microbial tree of life. rrn copy number variability provides valuable insight into microbial life history, but introduces systematic bias when measuring community composition in molecular surveys. Here we present an update to the ribosomal RNA operon copy number database (rrnDB), a publicly available, curated resource for copy number information for bacteria and archaea. The redesigned rrnDB (http://rrndb.umms.med.umich.edu/) brings a substantial increase in the number of genomes described, improved curation, mapping of genomes to both NCBI and RDP taxonomies, and refined tools for querying and analyzing these data. With these changes, the rrnDB is better positioned to remain a comprehensive resource under the torrent of microbial genome sequencing. The enhanced rrnDB will contribute to the analysis of molecular surveys and to research linking genomic characteristics to life history.

Marine aerosol as a possible source for endotoxins in coastal areas

Marine aerosols, that are very common in the highly populated coastal cities and communities, may contain biological constituents. Some of this biological fraction of marine aerosols, such as cyanobacteria and plankton debris, may influence human health by inflammation and allergic reactions when inhaled. In this study we identify and compare sources for endotoxins sampled on filters in an on-shore and more-inland site. Filter analysis included endotoxin content, total bacteria, gram-negative bacteria and cyanobacteria genome concentrations as well as ion content in order to identify possible sources for the endotoxins. Satellite images of chlorophyll-a levels and back trajectory analysis were used to further study the cyanobacteria blooms in the sea, close to the trajectory of the sampled air. The highest endotoxin concentrations found in the shoreline site were during winter (3.23±0.17 EU/m(3)), together with the highest cyanobacteria genome (1065.5 genome/m(3)). The elevated endotoxin concentrations were significantly correlated with cyanobacterial levels scaled to the presence of marine aerosol (r=0.90), as well as to chlorophyll-a (r=0.96). Filters sampled further inland showed lower and non-significant correlation between endotoxin and cyanobacteria (r=0.70, P value=0.19), suggesting decrease in marine-originated endotoxin, with possible contributions from other sources of gram-negative non-cyanobacteria. We conclude that marine cyanobacteria may be a dominant contributor to elevated endotoxin levels in coastal areas.

Multilocus sequence typing of Mycoplasma bovis reveals host-specific genotypes in cattle versus bison

Mycoplasma bovis is a primary agent of mastitis, pneumonia and arthritis in cattle and the bacterium most frequently isolated from the polymicrobial syndrome known as bovine respiratory disease complex. Recently, M. bovis has emerged as a significant health problem in bison, causing necrotic pharyngitis, pneumonia, dystocia and abortion. Whether isolates from cattle and bison comprise genetically distinct populations is unknown. This study describes the development of a highly discriminatory multilocus sequencing typing (MLST) method for M. bovis and its use to investigate the population structure of the bacterium. Genome sequences from six M. bovis isolates were used for selection of gene targets. Seven of 44 housekeeping genes initially evaluated were selected as targets on the basis of sequence variability and distribution within the genome. For each gene target sequence, four to seven alleles could be distinguished that collectively define 32 sequence types (STs) from a collection of 94 cattle isolates and 42 bison isolates. A phylogeny based on concatenated target gene sequences of each isolate revealed that bison isolates are genetically distinct from strains that infect cattle, suggesting recent disease outbreaks in bison may be due to the emergence of unique genetic variants. No correlation was found between ST and disease presentation or geographic origin. MLST data reported here were used to populate a newly created and publicly available, curated database to which researchers can contribute. The MLST scheme and database provide novel tools for exploring the population structure of M. bovis and tracking the evolution and spread of strains.

Bacterial diversity dynamics associated with different diets and different primer pairs in the rumen of Kankrej cattle

The ruminal microbiome in herbivores plays a dominant role in the digestion of lignocellulose and has potential to improve animal productivity. Kankrej cattle, a popular native breed of the Indian subcontinent, were used to investigate the effect of different dietary treatments on the bacterial diversity in ruminal fractions using different primer pairs. Two groups of four cows were assigned to two primary diets of either dry or green forages. Each group was fed one of three dietary treatments for six weeks each. Dietary treatments were; K1 (50% dry/green roughage: 50% concentrate), K2 (75% dry/green roughage: 25% concentrate) and K3 (100% dry/green roughage). Rumen samples were collected using stomach tube at the end of each dietary period and separated into solid and liquid fractions. The DNA was extracted and amplified for V1–V3, V4–V5 and V6–V8 hypervariable regions using P1, P2 and P3 primer pairs, sequenced on a 454 Roche platform and analyzed using QIIME. Community compositions and the abundance of most bacterial lineages were driven by interactions between primer pair, dietary treatment and fraction. The most abundant bacterial phyla identified were Bacteroidetes and Firmicutes however, the abundance of these phyla varied between different primer pairs; in each primer pair the abundance was dependent on the dietary treatment and fraction. The abundance of Bacteroidetes in cattle receiving K1 treatment indicate their diverse functional capabilities in the digestion of both carbohydrate and protein while the predominance of Firmicutes in the K2 and K3 treatments signifies their metabolic role in fibre digestion. It is apparent that both liquid and solid fractions had distinct bacterial community patterns (P<0.001) congruent to changes in the dietary treatments. It can be concluded that the P1 primer pair flanking the V1–V3 hyper-variable region provided greater species richness and diversity of bacterial populations in the rumen of Kankrej cattle.

Dual investigation of methanogenic processes by quantitative PCR and quantitative microscopic fingerprinting

Monitoring of methanogenic communities in anaerobic digesters using molecular-based methods is very attractive but can be cost-intensive. A new and fast quantification method by microscopic image analysis was developed to accompany molecular-based methods. This digitalized method, called quantitative microscopic fingerprinting (QMF), enables quantification of active methanogenic cells (N mL(-1)) by their characteristic auto-fluorescence based on coenzyme F420 . QMF was applied to analyze the methanogenic communities in three biogas plant samples, and the results were compared with the relative proportion of gene copy numbers obtained with the quantitative PCR (qPCR). Analysis of QMF demonstrated dominance of Methanomicrobiales and Methanobacteriales in relation to the total methanogenic community in digesters operating at high ammonia concentrations, which corresponded to the results established by qPCR. Absolute microbial counts by QMF and the numbers obtained by qPCR were not always comparable. On the other hand, the restricted morphological analysis by QMF was enhanced by the capability of qPCR to identify microbes. Consequently, dual investigations of both methods are proposed to improve monitoring of anaerobic digesters. For a rough estimation of the methanogenic composition in anaerobic digesters, the QMF method seems to be a promising approach for the rapid detection of microbial changes.

Seasonality and depth distribution of the abundance and activity of ammonia oxidizing microorganisms in marine coastal sediments (North Sea)

Microbial processes such as nitrification and anaerobic ammonium oxidation (anammox) are important for nitrogen cycling in marine sediments. Seasonal variations of archaeal and bacterial ammonia oxidizers (AOA and AOB) and anammox bacteria, as well as the environmental factors affecting these groups, are not well studied. We have examined the seasonal and depth distribution of the abundance and potential activity of these microbial groups in coastal marine sediments of the southern North Sea. This was achieved by quantifying specific intact polar lipids as well as the abundance and gene expression of their 16S rRNA gene, the ammonia monooxygenase subunit A (amoA) gene of AOA and AOB, and the hydrazine synthase (hzsA) gene of anammox bacteria. AOA, AOB, and anammox bacteria were detected and transcriptionally active down to 12 cm sediment depth. In all seasons, the abundance of AOA was higher compared to the AOB abundance suggesting that AOA play a more dominant role in aerobic ammonia oxidation in these sediments. Anammox bacteria were abundant and active even in oxygenated and bioturbated parts of the sediment. The abundance of AOA and AOB was relatively stable with depth and over the seasonal cycle, while anammox bacteria abundance and transcriptional activity were highest in August. North Sea sediments thus seem to provide a common, stable, ecological niche for AOA, AOB, and anammox bacteria.

Ether- and ester-bound iso-diabolic acid and other lipids in members of acidobacteria subdivision 4

Recently, iso-diabolic acid (13,16-dimethyl octacosanedioic acid) has been identified as a major membrane-spanning lipid of subdivisions 1 and 3 of the Acidobacteria, a highly diverse phylum within the Bacteria. This finding pointed to the Acidobacteria as a potential source for the bacterial glycerol dialkyl glycerol tetraethers that occur ubiquitously in peat, soil, lakes, and hot springs. Here, we examined the lipid composition of seven phylogenetically divergent strains of subdivision 4 of the Acidobacteria, a bacterial group that is commonly encountered in soil. Acid hydrolysis of total cell material released iso-diabolic acid derivatives in substantial quantities (11 to 48% of all fatty acids). In contrast to subdivisions 1 and 3 of the Acidobacteria, 6 out of the 7 species of subdivision 4 (excepting "Candidatus Chloracidobacterium thermophilum") contained iso-diabolic acid ether bound to a glycerol in larger fractional abundance than iso-diabolic acid itself. This is in agreement with the analysis of intact polar lipids (IPLs) by high-performance liquid chromatography-mass spectrometry (HPLC-MS), which showed the dominance of mixed ether-ester glycerides. iso-Diabolic acid-containing IPLs were not identified, because these IPLs are not released with a Bligh-Dyer extraction, as observed before when studying lipid compositions of subdivisions 1 and 3 of the Acidobacteria. The presence of ether bonds in the membrane lipids does not seem to be an adaptation to temperature, because the five mesophilic isolates contained a larger amount of ether lipids than the thermophile "Ca. Chloracidobacterium thermophilum." Furthermore, experiments with Pyrinomonas methylaliphatogenes did not reveal a major influence of growth temperature over the 50 to 69°C range.

Impact of next generation sequencing techniques in food microbiology

Understanding the Maxam-Gilbert and Sanger sequencing as the first generation, in recent years there has been an explosion of newly-developed sequencing strategies, which are usually referred to as next generation sequencing (NGS) techniques. NGS techniques have high-throughputs and produce thousands or even millions of sequences at the same time. These sequences allow for the accurate identification of microbial taxa, including uncultivable organisms and those present in small numbers. In specific applications, NGS provides a complete inventory of all microbial operons and genes present or being expressed under different study conditions. NGS techniques are revolutionizing the field of microbial ecology and have recently been used to examine several food ecosystems. After a short introduction to the most common NGS systems and platforms, this review addresses how NGS techniques have been employed in the study of food microbiota and food fermentations, and discusses their limits and perspectives. The most important findings are reviewed, including those made in the study of the microbiota of milk, fermented dairy products, and plant-, meat- and fish-derived fermented foods. The knowledge that can be gained on microbial diversity, population structure and population dynamics via the use of these technologies could be vital in improving the monitoring and manipulation of foods and fermented food products. They should also improve their safety.

Chemolithotrophic primary production in a subglacial ecosystem

Glacial comminution of bedrock generates fresh mineral surfaces capable of sustaining chemotrophic microbial communities under the dark conditions that pervade subglacial habitats. Geochemical and isotopic evidence suggests that pyrite oxidation is a dominant weathering process generating protons that drive mineral dissolution in many subglacial systems. Here, we provide evidence correlating pyrite oxidation with chemosynthetic primary productivity and carbonate dissolution in subglacial sediments sampled from Robertson Glacier (RG), Alberta, Canada. Quantification and sequencing of ribulose-1,5-bisphosphate carboxylase/oxygenase (RuBisCO) transcripts suggest that populations closely affiliated with Sideroxydans lithotrophicus, an iron sulfide-oxidizing autotrophic bacterium, are abundant constituents of microbial communities at RG. Microcosm experiments indicate sulfate production during biological assimilation of radiolabeled bicarbonate. Geochemical analyses of subglacial meltwater indicate that increases in sulfate levels are associated with increased calcite and dolomite dissolution. Collectively, these data suggest a role for biological pyrite oxidation in driving primary productivity and mineral dissolution in a subglacial environment and provide the first rate estimate for bicarbonate assimilation in these ecosystems. Evidence for lithotrophic primary production in this contemporary subglacial environment provides a plausible mechanism to explain how subglacial communities could be sustained in near-isolation from the atmosphere during glacial-interglacial cycles.

Microbiome analysis - from technical advances to biological relevance

The development of culture-independent techniques and next-generation sequencing has led to a staggering rise in the number of microbiome studies over the last decade. Although it remains important to identify the taxa of microbes present in a variety of environmental samples, including the gut microbiomes of healthy and diseased individuals, the next stage of microbiome research will need to focus on uncovering the role of the microbiome rather than its mere composition. Here, we introduce techniques that go beyond identifying the taxa present within a sample and examine the biological function of the microbiome or the host-microbiome interaction.

Intertidal epilithic bacteria diversity changes along a naturally occurring carbon dioxide and pH gradient

Intertidal epilithic bacteria communities are important components of coastal ecosystems, yet few studies have assessed their diversity and how it may be affected by changing environmental parameters. Submarine CO2 seeps produce localised areas of CO2-enriched seawater with reduced pH levels. We utilised the seawater pH/CO2 gradient at Levante Bay (Italy) to test the hypothesis that epilithic bacteria communities are modified by exposure to seawater with the varying chemical parameters. Biofilms were sampled from three sites exposed to seawater with different pH/CO2 levels and diversity determined using high-throughput sequencing of 16S rRNA genes. Seawater pCO2 concentrations were increased from ambient at site 1 to 621 μatm at site 2 and 1654 μatm site 3, similar to the predicated future oceans beyond 2050 and 2150, respectively. Alpha diversity of total bacteria communities and Cyanobacteria communities was significantly different between sites (anova P < 0.05). Comparison between sites showed that bacteria communities and Cyanobacteria communities were significantly different (anosim P < 0.01; permanova P < 0.01). Proteobacteria, Bacteroidetes and Cyanobacteria dominated all communities; however, there were differences between sites in the relative abundance of specific orders. This study provides the most detailed assessment of intertidal epilithic bacteria diversity and shows that diversity is significantly different along a seawater pH/CO2 gradient. This information supports the evaluation of the impacts of future ocean acidification on coastal marine ecosystems.

Mycobacterial phylogenomics: an enhanced method for gene turnover analysis reveals uneven levels of gene gain and loss among species and gene families

Species of the genus Mycobacterium differ in several features, from geographic ranges, and degree of pathogenicity, to ecological and host preferences. The recent availability of several fully sequenced genomes for a number of these species enabled the comparative study of the genetic determinants of this wide lifestyle diversity. Here, we applied two complementary phylogenetic-based approaches using information from 19 Mycobacterium genomes to obtain a more comprehensive view of the evolution of this genus. First, we inferred the phylogenetic relationships using two new approaches, one based on a Mycobacterium-specific amino acid substitution matrix and the other on a gene content dissimilarity matrix. Then, we utilized our recently developed gain-and-death stochastic models to study gene turnover dynamics in this genus in a maximum-likelihood framework. We uncovered a scenario that differs markedly from traditional 16S rRNA data and improves upon recent phylogenomic approaches. We also found that the rates of gene gain and death are high and unevenly distributed both across species and across gene families, further supporting the utility of the new models of rate heterogeneity applied in a phylogenetic context. Finally, the functional annotation of the most expanded or contracted gene families revealed that the transposable elements and the fatty acid metabolism-related gene families are the most important drivers of gene content evolution in Mycobacterium.

Effect of dried distillers’ grains and solubles when replacing corn or soybean meal on rumen microbial growth in vitro as measured using DNA as a microbial marker

Castillo-Lopez, E., Klopfenstein, T. J., Fernando, S. C. and Kononoff, P. J. 2014. Effect of dried distillers’ grains and solubles when replacing corn or soybean meal on rumen microbial growth in vitro as measured using DNA as a microbial marker. Can. J. Anim. Sci. 94: 349–356. The objectives were to evaluate the use of rDNA markers to measure the effects of dried distillers’ grains with solubles (DDGS) and the potential treatment×time interaction on microbial crude protein (MCP) synthesis in vitro and secondly to measure the contribution of yeast based protein originating from DDGS. Treatments were: (1) CONT, control with no DDGS, but with alfalfa hay, corn silage, ground corn (GC) and soybean meal (SBM) included at 25% (DM basis); (2) LOWCORN, 20% DDGS (DM basis) replacing GC; (3) LOWSBM, 20% DDGS (DM basis) replacing SBM; and (4) LOWCORNSBM, 20% DDGS (DM basis) replacing 10% GC and 10% SBM. Treatments (0.5 g) were incubated in 50 mL of inoculum in duplicate. At 0, 4, 16, 32, 48 and 96 h of fermentation total DNA was extracted from each treatment and MCP was measured using rDNA markers. The sum of bacterial crude protein (BCP) and protozoal crude protein (PCP) was considered as MCP. Data were analyzed as a completely randomized design. The treatment×time interaction was tested and the SLICE option was included to evaluate the effect of treatment at each fermentation time point. There was a tendency to a treatment×time interaction (P=0.07) for MCP. Specifically, at 16 h, LOWCORNSBM yielded greater (P<0.05) MCP compared to either CONT or LOWCORN with estimates of 68.5, 33.8 and 23.3±8.9 mg g–1DM, for LOWCORNSBM, CONT and LOWCORN, respectively. At 48 h, however, LOWCORN yielded greater MCP (P<0.05) compared with LOWSBM with estimates of 72.2 and 32.5±8.9 mg g–1DM, for LOWCORN and LOWSBM, respectively. Yeast crude protein (YCP) was not affected (P=0.21) and averaged 0.04±0.02 mg g–1of substrate (DM basis). Overall, rDNA markers were effective for quantifying MCP, but further research on the methodology is needed. With DDGS inclusion, MCP was maintained; however, yeast cells were extensively degraded during fermentation.

Ecotype diversity and conversion in Photobacterium profundum strains

Photobacterium profundum is a cosmopolitan marine bacterium capable of growth at low temperature and high hydrostatic pressure. Multiple strains of P. profundum have been isolated from different depths of the ocean and display remarkable differences in their physiological responses to pressure. The genome sequence of the deep-sea piezopsychrophilic strain Photobacterium profundum SS9 has provided some clues regarding the genetic features required for growth in the deep sea. The sequenced genome of Photobacterium profundum strain 3TCK, a non-piezophilic strain isolated from a shallow-water environment, is now available and its analysis expands the identification of unique genomic features that correlate to environmental differences and define the Hutchinsonian niche of each strain. These differences range from variations in gene content to specific gene sequences under positive selection. Genome plasticity between Photobacterium bathytypes was investigated when strain 3TCK-specific genes involved in photorepair were introduced to SS9, demonstrating that horizontal gene transfer can provide a mechanism for rapid colonisation of new environments.

Meta'omic analytic techniques for studying the intestinal microbiome

Nucleotide sequencing has become increasingly common and affordable, and is now a vital tool for studies of the human microbiome. Comprehensive microbial community surveys such as MetaHit and the Human Microbiome Project have described the composition and molecular functional profile of the healthy (normal) intestinal microbiome. This knowledge will increase our ability to analyze host and microbial DNA (genome) and RNA (transcriptome) sequences. Bioinformatic and statistical tools then can be used to identify dysbioses that might cause disease, and potential treatments. Analyses that identify perturbations in specific molecules can leverage thousands of culture-based isolate genomes to contextualize culture-independent sequences, or may integrate sequence data with whole-community functional assays such as metaproteomic or metabolomic analyses. We review the state of available systems-level models for studies of the intestinal microbiome, along with analytic techniques and tools that can be used to determine its functional capabilities in healthy and unhealthy individuals.

Decline in presumptively protective gut bacterial species and metabolites are paradoxically associated with disease improvement in pediatric Crohn's disease during enteral nutrition

BACKGROUND

The gut microbiota is implicated in the pathogenesis of Crohn's disease (CD). Exclusive enteral nutrition (EEN) is a successful treatment, but its mode of action remains unknown. This study assessed serial changes in the fecal microbiota milieu during EEN.

METHODS

Five fecal samples were collected from CD children: 4 during EEN (start, 15, 30, end EEN approximately 60 days) and the fifth on habitual diet. Two samples were collected from healthy control subjects. Fecal pH, bacterial metabolites, global microbial diversity abundance, composition stability, and quantitative changes of total and 7 major bacterial groups previously implicated in CD were measured.

RESULTS

Overall, 68 samples were from 15 CD children and 40 from 21 control subjects. Fecal pH and total sulfide increased and butyric acid decreased during EEN (all P < 0.05). Global bacterial diversity abundance decreased (P < 0.05); a higher degree of microbiota composition stability was seen in control subjects than in CD children during EEN (at P ≤ 0.008). Faecalibacterium prausnitzii spp concentration significantly decreased after 30 days on EEN (P < 0.05). In patients who responded to EEN, the magnitude of the observed changes was greater and the concentration of Bacteroides/Prevotella group decreased (P < 0.05). All these changes reverted to pretreatment levels on free diet, and EEN microbiota diversity increased when the children returned to their free diet.

CONCLUSIONS

EEN impacts on gut microbiota composition and changes fecal metabolic activity. It is difficult to infer a causative association between such changes and disease improvement, but the results do challenge the current perception of a protective role for F. prausnitzii in CD.

Effect of hillslope position and manure application rates on the persistence of fecal source tracking indicators in an agricultural soil

The influence of liquid dairy manure (LDM) application rates (12.5 and 25 kL ha) and soil type on the decay rates of library-independent fecal source tracking markers (host-associated and mitochondrial DNA) and persistent (>58 d) population structure was examined in a field study. The soils compared were an Aquic Haplorthod and a Typic Haplorthod in Nova Scotia, Canada, that differed according to landscape position and soil moisture regime. Soil type and LDM application rate did not influence decay rates (0.045-0.057 d). population structure, in terms of the occurrence of abundance of strain types, varied according to soil type ( = 0.012) but did not vary by LDM application rate ( = 0.121). Decay of ruminant-specific (BacR), bovine-specific (CowM2), and mitochondrial DNA (AcytB) markers was analyzed for 13 d after LDM application. The decay rates of BacR were greater under high-LDM application rates (0.281-0.358 d) versus low-LDM application rates (0.212-0.236 d) but were unaffected by soil type. No decay rates could be calculated for the CowM2 marker because it was undetectable within 6 d after manure application. Decay rates for AcytB were lower for the Aquic Haplorthod (0.088-0.100 d), with higher moisture status compared with the Typic Haplorthod (0.135 d). Further investigation into the decay of fecal source tracking indicators in agricultural field soils is warranted to assess the influence of soil type and agronomic practice on the differential decay of relevant markers and the likelihood of transport in runoff.

Eisenbergiella tayi gen. nov., sp. nov., isolated from human blood

A catalase-positive, rod-shaped, non-proteolytic, non-motile, anaerobic bacterial strain, designated B086562T, was isolated from a blood culture of an 84-year-old male patient in Israel. According to 16S rRNA gene sequence phylogeny, this strain has no known close relatives among recognized bacteria but should be placed within the family Lachnospiraceae . The most closely related recognized bacteria were from the ‘ Clostridium clostridioforme group’: C. clostridioforme (92.4 %) and Clostridium bolteae (92.3 %). The isolate produced butyrate, lactate, acetate and succinate as major metabolic end products. The major fatty acids were C16 : 0 and C18 : 1 cis 9 DMA and the DNA G+C content was 46.0 mol%. On the basis of the phenotypic properties and phylogenetic distinctiveness, the blood isolate represents a novel species of a new genus in the family Lachnospiraceae , for which the name Eisenbergiella tayi gen. nov., sp. nov. is proposed. The type strain of Eisenbergiella tayi is B086562T ( = LMG 27400T = DSM 26961T = ATCC BAA-2558T).

A novel molecular typing method of Mycobacteria based on DNA barcoding visualization

Different subtypes of Mycobacterium tuberculosis (MTB) may induce diverse severe human infections, and some of their symptoms are similar to other pathogenes, e.g. Nontuberculosis mycobacteria (NTM). So determination of mycobacterium subtypes facilitates the effective control of MTB infection and proliferation. This study exploits a novel DNA barcoding visualization method for molecular typing of 17 mycobacteria genomes published in the NCBI prokaryotic genome database. Three mycobacterium genes (Rv0279c, Rv3508 and Rv3514) from the PE/PPE family of MT Band were detected to best represent the inter-strain pathogenetic variations. An accurate and fast MTB substrain typing method was proposed based on the combination of the aforementioned three biomarker genes and the 16S rRNA gene. The protocol of establishing a bacterial substrain typing system used in this study may also be applied to the other pathogenes.

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.

Microbial genomics, transcriptomics and proteomics: new discoveries in decomposition research using complementary methods

Molecular methods for the analysis of biomolecules have undergone rapid technological development in the last decade. The advent of next-generation sequencing methods and improvements in instrumental resolution enabled the analysis of complex transcriptome, proteome and metabolome data, as well as a detailed annotation of microbial genomes. The mechanisms of decomposition by model fungi have been described in unprecedented detail by the combination of genome sequencing, transcriptomics and proteomics. The increasing number of available genomes for fungi and bacteria shows that the genetic potential for decomposition of organic matter is widespread among taxonomically diverse microbial taxa, while expression studies document the importance of the regulation of expression in decomposition efficiency. Importantly, high-throughput methods of nucleic acid analysis used for the analysis of metagenomes and metatranscriptomes indicate the high diversity of decomposer communities in natural habitats and their taxonomic composition. Today, the metaproteomics of natural habitats is of interest. In combination with advanced analytical techniques to explore the products of decomposition and the accumulation of information on the genomes of environmentally relevant microorganisms, advanced methods in microbial ecophysiology should increase our understanding of the complex processes of organic matter transformation.

Genome degeneration and adaptation in a nascent stage of symbiosis

Symbiotic associations between animals and microbes are ubiquitous in nature, with an estimated 15% of all insect species harboring intracellular bacterial symbionts. Most bacterial symbionts share many genomic features including small genomes, nucleotide composition bias, high coding density, and a paucity of mobile DNA, consistent with long-term host association. In this study, we focus on the early stages of genome degeneration in a recently derived insect-bacterial mutualistic intracellular association. We present the complete genome sequence and annotation of Sitophilus oryzae primary endosymbiont (SOPE). We also present the finished genome sequence and annotation of strain HS, a close free-living relative of SOPE and other insect symbionts of the Sodalis-allied clade, whose gene inventory is expected to closely resemble the putative ancestor of this group. Structural, functional, and evolutionary analyses indicate that SOPE has undergone extensive adaptation toward an insect-associated lifestyle in a very short time period. The genome of SOPE is large in size when compared with many ancient bacterial symbionts; however, almost half of the protein-coding genes in SOPE are pseudogenes. There is also evidence for relaxed selection on the remaining intact protein-coding genes. Comparative analyses of the whole-genome sequence of strain HS and SOPE highlight numerous genomic rearrangements, duplications, and deletions facilitated by a recent expansion of insertions sequence elements, some of which appear to have catalyzed adaptive changes. Functional metabolic predictions suggest that SOPE has lost the ability to synthesize several essential amino acids and vitamins. Analyses of the bacterial cell envelope and genes encoding secretion systems suggest that these structures and elements have become simplified in the transition to a mutualistic association.

CopyRighter: a rapid tool for improving the accuracy of microbial community profiles through lineage-specific gene copy number correction

BACKGROUND

Culture-independent molecular surveys targeting conserved marker genes, most notably 16S rRNA, to assess microbial diversity remain semi-quantitative due to variations in the number of gene copies between species.

RESULTS

Based on 2,900 sequenced reference genomes, we show that 16S rRNA gene copy number (GCN) is strongly linked to microbial phylogenetic taxonomy, potentially under-representing Archaea in amplicon microbial profiles. Using this relationship, we inferred the GCN of all bacterial and archaeal lineages in the Greengenes database within a phylogenetic framework. We created CopyRighter, new software which uses these estimates to correct 16S rRNA amplicon microbial profiles and associated quantitative (q)PCR total abundance. CopyRighter parses microbial profiles and, because GCN estimates are pre-computed for all taxa in the reference taxonomy, rapidly corrects GCN bias. Software validation with in silico and in vitro mock communities indicated that GCN correction results in more accurate estimates of microbial relative abundance and improves the agreement between metagenomic and amplicon profiles. Analyses of human-associated and anaerobic digester microbiomes illustrate that correction makes tangible changes to estimates of qPCR total abundance, α and β diversity, and can significantly change biological interpretation. For example, human gut microbiomes from twins were reclassified into three rather than two enterotypes after GCN correction.

CONCLUSIONS

The CopyRighter bioinformatic tools permits rapid correction of GCN in microbial surveys, resulting in improved estimates of microbial abundance, α and β diversity.

POGO-DB--a database of pairwise-comparisons of genomes and conserved orthologous genes

POGO-DB (http://pogo.ece.drexel.edu/) provides an easy platform for comparative microbial genomics. POGO-DB allows users to compare genomes using pre-computed metrics that were derived from extensive computationally intensive BLAST comparisons of >2000 microbes. These metrics include (i) average protein sequence identity across all orthologs shared by two genomes, (ii) genomic fluidity (a measure of gene content dissimilarity), (iii) number of ‘orthologs’ shared between two genomes, (iv) pairwise identity of the 16S ribosomal RNA genes and (v) pairwise identity of an additional 73 marker genes present in >90% prokaryotes. Users can visualize these metrics against each other in a 2D plot for exploratory analysis of genome similarity and of how different aspects of genome similarity relate to each other. The results of these comparisons are fully downloadable. In addition, users can download raw BLAST results for all or user-selected comparisons. Therefore, we provide users with full flexibility to carry out their own downstream analyses, by creating easy access to data that would normally require heavy computational resources to generate. POGO-DB should prove highly useful for researchers interested in comparative microbiology and benefit the microbiome/metagenomic communities by providing the information needed to select suitable phylogenetic marker genes within particular lineages.