Effects of polymerase,template dilution and cycle number on PCR based 16 S rRNA diversity analysis using the deep sequencing method

Background  

The primer and amplicon length have been found to affect PCR based estimates of microbial diversity by pyrosequencing, while other PCR conditions have not been addressed using any deep sequencing method. The present study determined the effects of polymerase, template dilution and PCR cycle number using the Solexa platform.     

Effect of primer mismatch, annealing temperature and PCR cycle number on 16S rRNA gene-targetting bacterial community analysis

In the attempt to explore complex bacterial communities of environmental samples, primers hybridizing to phylogenetically highly conserved regions of 16S rRNA genes are widely used, but differential amplification is a recognized problem. The biases associated with preferential amplification of multitemplate PCR were investigated using 'universal' bacteria-specific primers, focusing on the effect of primer mismatch, annealing temperature and PCR cycle number. The distortion of the template-to-product ratio was measured using predefined template mixtures and environmental samples by terminal restriction fragment length polymorphism analysis. When a 1 : 1 genomic DNA template mixture of two strains was used, primer mismatches inherent in the 63F primer presented a serious bias, showing preferential amplification of the template containing the perfectly matching sequence. The extent of the preferential amplification showed an almost exponential relation with increasing annealing temperature from 47 to 61 degrees C. No negative effect of the various annealing temperatures was observed with the 27F primer, with no mismatches with the target sequences. The number of PCR cycles had little influence on the template-to-product ratios. As a result of additional tests on environmental samples, the use of a low annealing temperature is recommended in order to significantly reduce preferential amplification while maintaining the specificity of PCR.     

Composition of Acridid gut bacterial communities as revealed by 16S rRNA gene analysis

The gut bacterial community from four species of feral locusts and grasshoppers was determined by denaturing gradient gel electrophoresis (DGGE) analysis of bacterial 16S rRNA gene fragments. The study revealed an effect of phase polymorphism on gut bacterial diversity in brown locusts from South Africa. A single bacterial phylotype, consistent with Citrobacter sp. dominated the gut microbiota of two sympatric populations of Moroccan and Italian locusts in Spain. There was evidence for Wollbachia sp. in the meadow grasshopper caught locally in the UK. Sequence analysis of DGGE products did not reveal evidence for unculturable bacteria and homologies suggested that bacterial species were principally Gammaproteobacteria from the family Enterobacteriaceae similar to those recorded previously in laboratory reared locusts.     

Massive parallel 16S rRNA gene pyrosequencing reveals highly diverse fecal bacterial and fungal communities in healthy dogs and cats

This study evaluated the fecal microbiota of 12 healthy pet dogs and 12 pet cats using bacterial and fungal tag-encoded FLX-Titanium amplicon pyrosequencing. A total of 120,406 pyrosequencing reads for bacteria (mean 5017) and 5359 sequences (one pool each for dogs and cats) for fungi were analyzed. Additionally, group-specific 16S rRNA gene clone libraries for Bifidobacterium spp. and lactic acid-producing bacteria (LAB) were constructed. The most abundant bacterial phylum was Firmicutes, followed by Bacteroidetes in dogs and Actinobacteria in cats. The most prevalent bacterial class in dogs and cats was Clostridia, dominated by the genera Clostridium (clusters XIVa and XI) and Ruminococcus. At the genus level, 85 operational taxonomic units (OTUs) were identified in dogs and 113 OTUs in cats. Seventeen LAB and eight Bifidobacterium spp. were detected in canine feces. Ascomycota was the only fungal phylum detected in cats, while Ascomycota, Basidiomycota, Glomeromycota, and Zygomycota were identified in dogs. Nacaseomyces was the most abundant fungal genus in dogs; Saccharomyces and Aspergillus were predominant in cats. At the genus level, 33 different fungal OTUs were observed in dogs and 17 OTUs in cats. In conclusion, this study revealed a highly diverse bacterial and fungal microbiota in canine and feline feces.     

The structure of bacterial communities in the western Arctic Ocean as revealed by pyrosequencing of 16S rRNA genes

Bacterial communities in the surface layer of the oceans consist of a few abundant phylotypes and many rare ones, most with unknown ecological functions and unclear roles in biogeochemical processes. To test hypotheses about relationships between abundant and rare phylotypes, we examined bacterial communities in the western Arctic Ocean using pyrosequence data of the V6 region of the 16S rRNA gene. Samples were collected from various locations in the Chukchi Sea, the Beaufort Sea and Franklin Bay in summer and winter. We found that bacterial communities differed between summer and winter at a few locations, but overall there was no significant difference between the two seasons in spite of large differences in biogeochemical properties. The sequence data suggested that abundant phylotypes remained abundant while rare phylotypes remained rare between the two seasons and among the Arctic regions examined here, arguing against the ‘seed bank’ hypothesis. Phylotype richness was calculated for various bacterial groups defined by sequence similarity or by phylogeny (phyla and proteobacterial classes). Abundant bacterial groups had higher within‐group diversity than rare groups, suggesting that the ecological success of a bacterial lineage depends on diversity rather than on the dominance of a few phylotypes. In these Arctic waters, in spite of dramatic variation in several biogeochemical properties, bacterial community structure was remarkably stable over time and among regions, and any variation was due to the abundant phylotypes rather than rare ones.     

Aberrant structures of fecal bacterial community in allergic infants profiled by 16S rRNA gene pyrosequencing

We investigated the correlation between fecal bacteria composition in early infancy and the prevalence of allergic diseases in late infancy. The fecal microbiota in the first 2?months was profiled using the 16S rRNA V6 short-tag sequences in the community and statistically compared between two groups of subjects who did and did not show allergic symptoms in the first 2?years (n?=?11 vs. 11). In the allergic group, genus Bacteroides at 1?month and genera Propionibacterium and Klebsiella at 2?months were more abundant, and genera Acinetobacter and Clostridium at 1?month were less abundant than in the nonallergic group. Allergic infants who showed high colonization of Bacteroides and/or Klebsiella showed less colonization of Clostridium perfringens/butyricum, suggesting antagonism between these bacterial groups in the gastrointestinal tract. It was also remarkable that the relative abundance of total Proteobacteria, excluding genus Klebsiella, was significantly lower in the allergic than in the nonallergic group at the age of 1?month. These results indicate that pyrosequence-based 16S rRNA gene profiling is valid to find the intestinal microbiotal disorder that correlates with allergy development in later life.     

Characterization of the fecal bacteria communities of forage-fed horses by pyrosequencing of 16S rRNA V4 gene amplicons

The diversity of the equine fecal bacterial community was evaluated using pyrosequencing of 16S rRNA gene amplicons. Fecal samples were obtained from horses fed cool-season grass hay. Fecal bacteria were characterized by amplifying the V4 region of bacterial 16S rRNA gene. Of 5898 mean unique sequences, a mean of 1510 operational taxonomic units were identified in the four fecal samples. Equine fecal bacterial richness was higher than that reported in humans, but lower than that reported in either cattle feces or soil. Bacterial classified sequences were assigned to 16 phyla, of which 10 were present in all samples. The largest number of reads belonged to Firmicutes (43.7% of total bacterial sequences), Verrucomicrobia (4.1%), Proteobacteria (3.8%), and Bacteroidetes (3.7%). The less abundant Actinobacteria, Cyanobacteria, and TM7 phyla presented here have not been previously described in the gut contents or feces of horses. Unclassified sequences represented 38.1% of total bacterial sequences; therefore, the equine fecal microbiome diversity is likely greater than that described. This is the first study to characterize the fecal bacterial community in horses by the use of 16S rRNA gene amplicon pyrosequencing, expanding our knowledge of the fecal microbiota of forage-fed horses.     

DNA extraction for 16S rRNA gene analysis to determine genetic diversity in deep sediment communities

A protocol was devised which permitted the extraction of DNA from deep marine sediments up to 503 m below the sea floor. These sediments have been laid down over the last 3 million years. 16S rRNA gene sequences were amplified from the DNA by the polymerase chain reaction. The details of the successful extraction and polymerase chain reaction methodology varied between samples from different depths. This emphasizes the attention to detail required to allow the diversity of bacteria in these deep sediments to be studied.     

Effect of genome size and rrn gene copy number on PCR amplification of 16S rRNA genes from a mixture of bacterial species.

In order to assess the effect of genome size and number of 16S rRNA genes (rDNAs) on the quantities of PCR-generated partial 16S rDNA fragments, equimolar amounts of DNA from pairs of different species for which these parameters are known were subjected to gene amplification. The experimentally determined ratio of PCR products obtained, as determined by image analysis of SYBR-Green I-stained amplification products, corresponded well with the predicted ratio calculated from the number of rrn genes per equimolar amounts of DNA in mixtures of Escherichia coli and "Thermus thermophilus" and of Pseudomonas aeruginosa and "T. thermophilus." The values for the pair of Bacillus subtilis and "T. thermophilus" showed greater deviations from the predicted value. The dependence of the amount of 16S rDNA amplification product on these two parameters makes it impossible to quantify the number of species represented in 16S rDNA clone libraries of environmental samples as long as these two parameters are unknown for the species present.     

A comparison of rpoB and 16S rRNA as markers in pyrosequencing studies of bacterial diversity

Background

The 16S rRNA gene is the gold standard in molecular surveys of bacterial and archaeal diversity, but it has the disadvantages that it is often multiple-copy, has little resolution below the species level and cannot be readily interpreted in an evolutionary framework. We compared the 16S rRNA marker with the single-copy, protein-coding rpoB marker by amplifying and sequencing both from a single soil sample. Because the higher genetic resolution of the rpoB gene prohibits its use as a universal marker, we employed consensus-degenerate primers targeting the Proteobacteria.

Methodology/Principal Findings

Pyrosequencing can be problematic because of the poor resolution of homopolymer runs. As these erroneous runs disrupt the reading frame of protein-coding sequences, removal of sequences containing nonsense mutations was found to be a valuable filter in addition to flowgram-based denoising. Although both markers gave similar estimates of total diversity, the rpoB marker revealed more species, requiring an order of magnitude fewer reads to obtain 90% of the true diversity. The application of population genetic methods was demonstrated on a particularly abundant sequence cluster.

Conclusions/Significance

The rpoB marker can be a complement to the 16S rRNA marker for high throughput microbial diversity studies focusing on specific taxonomic groups. Additional error filtering is possible and tests for recombination or selection can be employed.     

Profiling of the bacteria responsible for pyogenic liver abscess by 16S rRNA gene pyrosequencing

Pyogenic liver abscess (PLA) is a severe disease with considerable mortality and is often polymicrobial. Understanding the pathogens that cause PLA is the basis for PLA treatment. Here, we profiled the bacterial composition in PLA fluid by pyrosequencing the 16S ribosomal RNA (rRNA) gene based on next-generation sequencing (NGS) technology to identify etiological agents of PLA and to provide information of their 16S rRNA sequences for application to DNA-based techniques in the hospital. Twenty patients with PLA who underwent percutaneous catheter drainage, abscess culture, and blood culture for isolates were included. Genomic DNAs from abscess fluids were subjected to polymerase chain reaction and pyrosequencing of the 16S rRNA gene with a 454 GS Junior System. The abscess and blood cultures were positive in nine (45%) and four (20%) patients, respectively. Pyrosequencing of 16S rRNA gene showed that 90% of the PLA fluid samples contained single or multiple genera of known bacteria such as Klebsiella, Fusobacterium, Streptococcus, Bacteroides, Prevotella, Peptostreptococcus, unassigned Enterobacteriaceae, and Dialister. Klebsiella was predominantly found in the PLA fluid samples. All samples that carried unassigned bacteria had 26.8% reads on average. We demonstrated that the occurrence of PLA was associated with eight known bacterial genera as well as unassigned bacteria and that 16S rRNA gene sequencing was more useful than conventional culture methods for accurate identification of bacterial pathogens from PLA.     

Phylogenetic analysis of the family Rhizobiaceae and related bacteria by sequencing of 16S rRNA gene using PCR and DNA sequencer

Abstract The 16S rRNA gene sequences of 19 strains covering 97% of the molecules were determined for the members of the family Rhizobiaceae and related bacteria by PCR and DNA sequencer. The three biovars of Agrobacterium were located separately, whereas Agrobacterium rubi clustered with A. tumefaciens . Phylogenetic locations for the species of the genera Rhizobium, Sinorhizobium, Agrobacterium, Phylobacterium, Mycoplana (M. dimorpha), Ochrobactrum, Brucella and Rochalimaea (a rickettsia) were intermingled with each other with the similarity values higher than 92%. The family Rhizobiaceae should be redefined including the above-mentioned genera despite the ability for plant association and nitrogen fixation. Bradyrhizobium japonicum and Mycoplana bullata were far remote from the other species and should be excluded from this family.     

Distinct ectomycorrhizospheres share similar bacterial communities as revealed by pyrosequencing-based analysis of 16S rRNA genes

Analysis of the 16S rRNA gene sequences generated from Xerocomus pruinatus and Scleroderma citrinum ectomycorrhizospheres revealed that similar bacterial communities inhabited the two ectomycorrhizospheres in terms of phyla and genera, with an enrichment of the Burkholderia genus. Compared to the bulk soil habitat, ectomycorrhizospheres hosted significantly more Alpha-, Beta-, and Gammaproteobacteria.     

Copy number of the 16S rRNA gene in Rickettsia prowazekii.

The obligate intracellular parasite, Rickettsia prowazekii, is a slowly growing bacterium with a doubling time of 8 to 12 h. The copy number of the 16S rRNA gene in the rickettsial chromosome was determined to be one. Genomic DNA from R. prowazekii was digested either by a variety of restriction enzymes known not to cut at any site in the rickettsial 16S rRNA gene or by a combination of these noncutting enzymes and SmaI, which cuts the gene only once. Only one DNA fragment in these digests hybridized to a biotinylated probe containing a portion of the rickettsial 16S rRNA gene. Moreover, the density of the rickettsial 16S rRNA gene fragment after hybridization was equal to the density of each of the seven 16S rRNA gene fragments in Escherichia coli.     

16S rRNA gene and 16S-23S rRNA gene internal transcribed spacer sequences analysis of the genus Myxococcus.

Phylogenetic relationships of the species belonging to the genus Myxococcus were elucidated based on the sequences of 16S rRNA genes and 16S-23S rRNA gene internal transcribed spacer (ITS) regions. The Myxococcus species were consequently classified into four distinct groups. The type strain of Myxococcus coralloides occupied an independent position (Group 1); it has been recently reclassified as Corallococcus coralloides. Group 2 comprised the type strains of both Myxococcus virescens and Myxococcus xanthus, and some strains assigned to Myxococcus flavescens. The type strain of M. flavescens was contained in Group 3 along with the strains of Myxococcus fulvus. Group 4 included the strains belonging to C. coralloides, M. fulvus, and M. stipitatus. The type strain of M. fulvus that was allocated outside Group 4 in the 16S rRNA gene tree belonged to Group 3 in the ITS tree. These results strongly suggest that the morphological characteristics of Myxococcus species are not consistent with the phylogenetic relationships. The Myxococcus species must therefore be redefined according to the phylogenetic relationships revealed in this study.     

Illumina sequencing of bacterial 16S rDNA and 16S rRNA reveals seasonal and species-specific variation in bacterial communities in four moss species

In order to better understand the factors that influence bacterial diversity and community composition in moss-associated bacteria, a study of bacterial communities in four moss species collected in three seasons was carried out via high-throughput sequencing of 16S rDNA and 16S rRNA. Moss species included Cratoneuron filicinum, Pylaisiella polyantha, Campyliadelphus polygamum, and Grimmia pilifera, with samples collected in May, July, and October 2015 from rocks at Beijing Songshan National Nature Reserve. In total, the bacterial richness and diversity were high regardless of moss species, sampling season, or data source (DNA vs. RNA). Bacterial sequences were assigned to a total of 558 OTUs and 279 genera in 16 phyla. Proteobacteria and Actinobacteria were the two most abundant phyla, and Cellvibrio, Lapillicoccus, Jatrophihabitans, Friedmanniella, Oligoflexus, and Bosea the most common genera in the samples. A clustering algorithm and principal coordinate analysis revealed that C. filicinum and C. polygamum had similar bacterial communities, as did P. polyantha and G. pilifera. Metabolically active bacteria showed the same pattern in addition to seasonal variation: bacterial communities were most similar in summer and autumn, looking at each moss species separately. In contrast, DNA profiles lacked obvious seasonal dynamics. A partial least squares discriminant analysis identified three groups of samples that correlated with differences in moss species resources. Although bacterial community composition did vary with the sampling season and data source, these were not the most important factors influencing bacterial communities. Previous reports exhibited that mosses have been widely used in biomonitoring of air pollution by enriching some substances or elements in the moss-tag technique and the abundant moss associated bacteria might also be important components involved in the related biological processes. Thus, this survey not only enhanced our understanding of the factors which influence microbial communities in mosses but also would be helpful for better use and development of the moss-tag technique in the environmental biomonitoring.