The Effects of Alignment Quality, Distance Calculation Method, Sequence Filtering, and Region on the Analysis of 16S rRNA Gene-Based Studies

Pyrosequencing of PCR-amplified fragments that target variable regions within the 16S rRNA gene has quickly become a powerful method for analyzing the membership and structure of microbial communities. This approach has revealed and introduced questions that were not fully appreciated by those carrying out traditional Sanger sequencing-based methods. These include the effects of alignment quality, the best method of calculating pairwise genetic distances for 16S rRNA genes, whether it is appropriate to filter variable regions, and how the choice of variable region relates to the genetic diversity observed in full-length sequences. I used a diverse collection of 13,501 high-quality full-length sequences to assess each of these questions. First, alignment quality had a significant impact on distance values and downstream analyses. Specifically, the greengenes alignment, which does a poor job of aligning variable regions, predicted higher genetic diversity, richness, and phylogenetic diversity than the SILVA and RDP-based alignments. Second, the effect of different gap treatments in determining pairwise genetic distances was strongly affected by the variation in sequence length for a region; however, the effect of different calculation methods was subtle when determining the sample's richness or phylogenetic diversity for a region. Third, applying a sequence mask to remove variable positions had a profound impact on genetic distances by muting the observed richness and phylogenetic diversity. Finally, the genetic distances calculated for each of the variable regions did a poor job of correlating with the full-length gene. Thus, while it is tempting to apply traditional cutoff levels derived for full-length sequences to these shorter sequences, it is not advisable. Analysis of β-diversity metrics showed that each of these factors can have a significant impact on the comparison of community membership and structure. Taken together, these results urge caution in the design and interpretation of analyses using pyrosequencing data

The Airborne Metagenome in an Indoor Urban Environment

The indoor atmosphere is an ecological unit that impacts on public health. To investigate the composition of organisms in this space, we applied culture-independent approaches to microbes harvested from the air of two densely populated urban buildings, from which we analyzed 80 megabases genomic DNA sequence and 6000 16S rDNA clones. The air microbiota is primarily bacteria, including potential opportunistic pathogens commonly isolated from human-inhabited environments such as hospitals, but none of the data contain matches to virulent pathogens or bioterror agents. Comparison of air samples with each other and nearby environments suggested that the indoor air microbes are not random transients from surrounding outdoor environments, but rather originate from indoor niches. Sequence annotation by gene function revealed specific adaptive capabilities enriched in the air environment, including genes potentially involved in resistance to desiccation and oxidative damage. This baseline index of air microbiota will be valuable for improving designs of surveillance for natural or man-made release of virulent pathogens

Culture-Independent Microbiological Analysis of Foley Urinary Catheter Biofilms

Background: Prevention of catheter-associated urinary tract infection (CAUTI), a leading cause of nosocomial disease, is complicated by the propensity of bacteria to form biofilms on indwelling medical devices [1,2,3,4,5]. Methodology/Principal Findings: To better understand the microbial diversity of these communities, we report the results of a culture-independent bacterial survey of Foley urinary catheters obtained from patients following total prostatectomy. Two patient subsets were analyzed, based on treatment or no treatment with systemic fluoroquinolone antibiotics during convalescence. Results indicate the presence of diverse polymicrobial assemblages that were most commonly observed in patients who did not receive systemic antibiotics. The communities typically contained both Gram-positive and Gramnegative microorganisms that included multiple potential pathogens. Conclusion/Significance: Prevention and treatment of CAUTI must take into consideration the possible polymicrobial nature of any particular infection

High prevalence of dna from non-H. pylori helicobacters in the gastric mucosa of venezuelan pet dogs and its histological alterations

Non-H. pylori helicobacters (NHPH) have been demonstrated as gastric spiral-shaped bacteria in specimens obtained from dogs; however, their roles in the pathogenesis of upper gastrointestinal disease have not yet been clearly established. The purpose of this study was to evaluate the prevalence of NHPH DNA in the gastric mucosa of dogs and its association with histopathology. Helicobacter was detected through histopathological techniques, PCR, and FISH analysis from fundic biopsies of twenty dogs with or without signs of gastrointestinal disease. PCR and FISH were based on partial 16S rRNA gene sequences. Nineteen dogs showed mild to marked gastritis in the fundus, and only one dog had a healthy gastric mucosa. NHPH DNA was detected in 18 dogs with gastritis and one with normal gastric mucosa. However, there was no significant correlation between the presence of NHPH DNA and the degree of gastritis. These results show a high prevalence of NHPH DNA in the gastric mucosa of dogs from Venezuela. Further studies are necessary to determine a possible association between a specific NHPH species and the degree of gastritis

Ovine pedomics : the first study of the ovine foot 16S rRNA-based microbiome

We report the first study of the bacterial microbiome of ovine interdigital skin based on 16S rRNA by pyrosequencing and conventional cloning with Sanger-sequencing. Three flocks were selected, one a flock with no signs of footrot or interdigital dermatitis, a second flock with interdigital dermatitis alone and a third flock with both interdigital dermatitis and footrot. The sheep were classified as having either healthy interdigital skin (H), interdigital dermatitis (ID) or virulent footrot (VFR). The ovine interdigital skin bacterial community varied significantly by flock and clinical condition. The diversity and richness of operational taxonomic units was greater in tissue from sheep with ID than H or VFR affected sheep. Actinobacteria, Bacteriodetes, Firmicutes and Proteobacteria were the most abundant phyla comprising 25 genera. Peptostreptococcus, Corynebacterium and Staphylococcus were associated with H, ID and VFR respectively. Sequences of Dichelobacter nodosus, the causal agent of ovine footrot, were not amplified due to mismatches in the 16S rRNA universal forward primer (27F). A specific real time PCR assay was used to demonstrate the presence of D. nodosus which was detected in all samples including the flock with no signs of ID or VFR. Sheep with ID had significantly higher numbers of D. nodosus (104-109 cells/g tissue) than those with H or VFR feet

Application of Two-Part Statistics for Comparison of Sequence Variant Counts

Investigation of microbial communities, particularly human associated communities, is significantly enhanced by the vast amounts of sequence data produced by high throughput sequencing technologies. However, these data create high-dimensional complex data sets that consist of a large proportion of zeros, non-negative skewed counts, and frequently, limited number of samples. These features distinguish sequence data from other forms of high-dimensional data, and are not adequately addressed by statistical approaches in common use. Ultimately, medical studies may identify targeted interventions or treatments, but lack of analytic tools for feature selection and identification of taxa responsible for differences between groups, is hindering advancement. The objective of this paper is to examine the application of a two-part statistic to identify taxa that differ between two groups. The advantages of the two-part statistic over common statistical tests applied to sequence count datasets are discussed. Results from the t-test, the Wilcoxon test, and the two-part test are compared using sequence counts from microbial ecology studies in cystic fibrosis and from cenote samples. We show superior performance of the two-part statistic for analysis of sequence data. The improved performance in microbial ecology studies was independent of study type and sequence technology used

Randomized clinical trial to evaluate the effect of fecal microbiota transplant for initial Clostridium difficile infection in intestinal microbiome

Objective The aim of this study was to evaluate the impact of fecal donor-unrelated donor mix (FMT-FURM) transplantation as first-line therapy for C. difficile infection (CDI) in intestinal microbiome. Methods We designed an open, two-arm pilot study with oral vancomycin (250mg every 6 h for 10–14 days) or FMT-FURM as treatments for the first CDI episode in hospitalized adult patients in Hospital Universitario “Dr. Jose Eleuterio Gonzalez”. Patients were randomized by a closed envelope method in a 1: 1 ratio to either oral vancomycin or FMT-FURM. CDI resolution was considered when there was a reduction on the Bristol scale of at least 2 points, a reduction of at least 50% in the number of bowel movements, absence of fever, and resolution of abdominal pain (at least two criteria). From each patient, a fecal sample was obtained at days 0, 3, and 7 after treatment. Specimens were cultured to isolate C. difficile, and isolates were characterized by PCR. Susceptibility testing of isolates was performed using the agar dilution method. Fecal samples and FMT-FURM were analyzed by 16S rRNA sequencing. Results We included 19 patients; 10 in the vancomycin arm and 9 in the FMT-FURM arm. However, one of the patients in the vancomycin arm and two patients in the FMT-FURM arm were eliminated. Symptoms resolved in 8/9 patients (88.9%) in the vancomycin group, while symptoms resolved in 4/7 patients (57.1%) after the first FMT-FURM dose (P = 0.26) and in 5/7 patients (71.4%) after the second dose (P = 0.55). During the study, no adverse effects attributable to FMT-FURM were observed in patients. Twelve isolates were recovered, most isolates carried tcdB, tcdA, cdtA, and cdtB, with an 18-bp deletion in tcdC. All isolates were resistant to ciprofloxacin and moxifloxacin but susceptible to metronidazole, linezolid, fidaxomicin, and tetracycline. In the FMT-FURM group, the bacterial composition was dominated by Firmicutes, Bacteroidetes, and Proteobacteria at all-time points and the microbiota were remarkably stable over time. The vancomycin group showed a very different pattern of the microbial composition when comparing to the FMT-FURM group over time. Conclusion The results of this preliminary study showed that FMT-FURM for initial CDI is associated with specific bacterial communities that do not resemble the donors’ sample.Peer reviewedFinal Published versio

Characterization of Trapped Lignin-Degrading Microbes in Tropical Forest Soil

Lignin is often the most difficult portion of plant biomass to degrade, with fungi generally thought to dominate during late stage decomposition. Lignin in feedstock plant material represents a barrier to more efficient plant biomass conversion and can also hinder enzymatic access to cellulose, which is critical for biofuels production. Tropical rain forest soils in Puerto Rico are characterized by frequent anoxic conditions and fluctuating redox, suggesting the presence of lignin-degrading organisms and mechanisms that are different from known fungal decomposers and oxygen-dependent enzyme activities. We explored microbial lignin-degraders by burying bio-traps containing lignin-amended and unamended biosep beads in the soil for 1, 4, 13 and 30 weeks. At each time point, phenol oxidase and peroxidase enzyme activity was found to be elevated in the lignin-amended versus the unamended beads, while cellulolytic enzyme activities were significantly depressed in lignin-amended beads. Quantitative PCR of bacterial communities showed more bacterial colonization in the lignin-amended compared to the unamended beads after one and four weeks, suggesting that the lignin supported increased bacterial abundance. The microbial community was analyzed by small subunit 16S ribosomal RNA genes using microarray (PhyloChip) and by high-throughput amplicon pyrosequencing based on universal primers targeting bacterial, archaeal, and eukaryotic communities. Community trends were significantly affected by time and the presence of lignin on the beads. Lignin-amended beads have higher relative abundances of representatives from the phyla Actinobacteria, Firmicutes, Acidobacteria and Proteobacteria compared to unamended beads. This study suggests that in low and fluctuating redox soils, bacteria could play a role in anaerobic lignin decomposition

Phylometrics: a pipeline for inferring phylogenetic trees from a sequence relationship network perspective

<p>Abstract</p> <p>Background</p> <p>Comparative sequence analysis of the 16S rRNA gene is frequently used to characterize the microbial diversity of environmental samples. However, sequence similarities do not always imply functional or evolutionary relatedness due to many factors, including unequal rates of change and convergence. Thus, relying on top BLASTN hits for phylogenetic studies may misrepresent the diversity of these constituents. Furthermore, attempts to circumvent this issue by including a large number of BLASTN hits per sequence in one tree to explore their relatedness presents other problems. For instance, the multiple sequence alignment will be poor and computationally costly if not relying on manual alignment, and it may be difficult to derive meaningful relationships from the resulting tree. Analyzing sequence relationship networks within collective BLASTN results, however, reveal sequences that are closely related despite low rank.</p> <p>Results</p> <p>We have developed a web application, Phylometrics, that relies on networks of collective BLASTN results (rather than single BLASTN hits) to facilitate the process of building phylogenetic trees in an automated, high-throughput fashion while offering novel tools to find sequences that are of significant phylogenetic interest with minimal human involvement. The application, which can be installed locally in a laboratory or hosted remotely, utilizes a simple wizard-style format to guide the user through the pipeline without necessitating a background in programming. Furthermore, Phylometrics implements an independent job queuing system that enables users to continue to use the system while jobs are run with little or no degradation in performance. </p> <p>Conclusions</p> <p>Phylometrics provides a novel data mining method to screen supplied DNA sequences and to identify sequences that are of significant phylogenetic interest using powerful analytical tools. Sequences that are identified as being similar to a number of supplied sequences may provide key insights into their functional or evolutionary relatedness. Users require the same basic computer skills as for navigating most internet applications.</p