Massively parallel rRNA gene sequencing exacerbates the potential for biased community diversity comparisons due to variable library sizes

Technologies for massively parallel sequencing are revolutionizing microbial ecology and are vastly increasing the scale of ribosomal RNA (rRNA) gene studies. Although pyrosequencing has increased the breadth and depth of possible rRNA gene sampling, one drawback is that the number of reads obtained per sample is difficult to control. Pyrosequencing libraries typically vary widely in the number of sequences per sample, even within individual studies, and there is a need to revisit the behaviour of richness estimators and diversity indices with variable gene sequence library sizes. Multiple reports and review papers have demonstrated the bias in non-parametric richness estimators (e.g. Chao1 and ACE) and diversity indices when using clone libraries. However, we found that biased community comparisons are accumulating in the literature. Here we demonstrate the effects of sample size on Chao1, ACE, CatchAll, Shannon, Chao-Shen and Simpson's estimations specifically using pyrosequencing libraries. The need to equalize the number of reads being compared across libraries is reiterated, and investigators are directed towards available tools for making unbiased diversity comparisons.     

Potential of pmoA Amplicon Pyrosequencing for Methanotroph Diversity Studies

We analyzed the potential of pmoA amplicon pyrosequencing compared to that of Sanger sequencing with paddy soils as a model environment. We defined operational taxonomic unit (OTU) cutoff values of 7% and 18%, reflecting methanotrophic species and major phylogenetic pmoA lineages, respectively. Major lineages were already well covered by clone libraries; nevertheless, pyrosequencing provided a higher level of diversity at the species level.     

A Method for Studying Protistan Diversity Using Massively Parallel Sequencing of V9 Hypervariable Regions of Small-Subunit Ribosomal RNA Genes

Background

Massively parallel pyrosequencing of amplicons from the V6 hypervariable regions of small-subunit (SSU) ribosomal RNA (rRNA) genes is commonly used to assess diversity and richness in bacterial and archaeal populations. Recent advances in pyrosequencing technology provide read lengths of up to 240 nucleotides. Amplicon pyrosequencing can now be applied to longer variable regions of the SSU rRNA gene including the V9 region in eukaryotes.

Methodology/Principal Findings

We present a protocol for the amplicon pyrosequencing of V9 regions for eukaryotic environmental samples for biodiversity inventories and species richness estimation. The International Census of Marine Microbes (ICoMM) and the Microbial Inventory Research Across Diverse Aquatic Long Term Ecological Research Sites (MIRADA-LTERs) projects are already employing this protocol for tag sequencing of eukaryotic samples in a wide diversity of both marine and freshwater environments.

Conclusions/Significance

Massively parallel pyrosequencing of eukaryotic V9 hypervariable regions of SSU rRNA genes provides a means of estimating species richness from deeply-sampled populations and for discovering novel species from the environment.     

Analysis of the Rumen Bacterial Diversity under two Different Diet Conditions using Denaturing Gradient Gel Electrophoresis,Random Sequencing,and Statistical Ecology Approaches

Bacterial community structure and diversity in the rumen of steers in conditions of hay and corn diets was assessed by in vitro retrieval and analysis of the variable region (V3) of 16S rDNA. Two types of libraries were generated in this study: DGGE libraries, which further were analysed by excising, reamplification, and sequencing, and random shotgun sequence libraries. Phylogenetic and sequence similarity analyses of the resultant 68 clone sequences in DGGE libraries revealed the presence of 42 operational taxonomic units (OTUs) or phylotypes defined as having more than 97% of sequence similarity. One hundred and thirty four clone sequences in shotgun libraries were clustered into 72 phylotypes. The phylotype similarity, diversity, richness, and evenness in these libraries were estimated using a variety of diversity indices. In relation to diet, the corn-fed animals displayed more diverse and rich bacterial populations, which were mostly contributed by CFB-related phylotypes. Proteobacteria were also numerically prevalent on this diet (27%) but were represented by a few phylotypes thus diminishing the overall diversity and species richness values. On hay diet, the principal contributors to general diversity and species richness appeared to be low-G + C gram-positives. Although the ruminal Treponemaes were encountered only in hay-fed animals, their impact on species diversity on hay diet was low because of the limited number of phylotypes.     

Microarray analysis and barcoded pyrosequencing provide consistent microbial profiles depending on the source of human intestinal samples

Large-scale and in-depth characterization of the intestinal microbiota necessitates application of high-throughput 16S rRNA gene-based technologies, such as barcoded pyrosequencing and phylogenetic microarray analysis. In this study, the two techniques were compared and contrasted for analysis of the bacterial composition in three fecal and three small intestinal samples from human individuals. As PCR remains a crucial step in sample preparation for both techniques, different forward primers were used for amplification to assess their impact on microbial profiling results. An average of 7,944 pyrosequences, spanning the V1 and V2 region of 16S rRNA genes, was obtained per sample. Although primer choice in barcoded pyrosequencing did not affect species richness and diversity estimates, detection of Actinobacteria strongly depended on the selected primer. Microbial profiles obtained by pyrosequencing and phylogenetic microarray analysis (HITChip) correlated strongly for fecal and ileal lumen samples but were less concordant for ileostomy effluent. Quantitative PCR was employed to investigate the deviations in profiling between pyrosequencing and HITChip analysis. Since cloning and sequencing of random 16S rRNA genes from ileostomy effluent confirmed the presence of novel intestinal phylotypes detected by pyrosequencing, especially those belonging to the Veillonella group, the divergence between pyrosequencing and the HITChip is likely due to the relatively low number of available 16S rRNA gene sequences of small intestinal origin in the DNA databases that were used for HITChip probe design. Overall, this study demonstrated that equivalent biological conclusions are obtained by high-throughput profiling of microbial communities, independent of technology or primer choice.     

乙肝治疗耐药基因突变的焦磷酸测序技术诊断

目的：建立焦磷酸测序技术检测拉米夫定和阿德福韦酯治疗乙肝所致乙肝病毒基因耐药突变的定量检测方法,为临床乙肝耐药诊断和治疗提供依据。方法：针对乙肝病毒DNA聚合酶基因序列上4个常见基因突变位点的6种突变形式,分别克隆构建野生型和突变型质粒作为标准品,应用生物信息学手段设计目标基因通用PCR引物和各突变点的焦磷酸测序引物,建立焦磷酸测序的突变检测方法。对接受拉米夫定、阿德福韦酯治疗的慢性乙型肝炎患者血清标本进行检测。结果：构建了乙肝病毒四种常见耐药性突变的标准株和变异株克隆,建立了分别或同时检测拉米夫定、阿德福韦酯耐药突变的焦磷酸测序方法,对68例临床耐药或疑似耐药的患者血清标本进行检测,双脱氧测序验证,检出拉米夫定耐药突变32例,阿德福韦酯耐药突变5例,其中焦磷酸测序检出20例为混合突变,而双脱氧测序显示为6例。结论：成功建立了焦磷酸测序定量检测拉米夫定、阿德福韦酯耐药基因突变的方法,构建了乙肝病毒耐药基因突变的标准质粒,为临床动态监测乙肝病毒变异病毒株、指导合理用药奠定了基础。     

Microbial communities in the surface mucopolysaccharide layer and the black band microbial mat of black band-diseased Siderastrea siderea

Microbial community profiles and species composition associated with two black band-diseased colonies of the coral Siderastrea siderea were studied by 16S rRNA-targeted gene cloning, sequencing, and amplicon-length heterogeneity PCR (LH-PCR). Bacterial communities associated with the surface mucopolysaccharide layer (SML) of apparently healthy tissues of the infected colonies, together with samples of the black band disease (BBD) infections, were analyzed using the same techniques for comparison. Gene sequences, ranging from 424 to 1,537 bp, were retrieved from all positive clones (n = 43 to 48) in each of the four clone libraries generated and used for comparative sequence analysis. In addition to LH-PCR community profiling, all of the clone sequences were aligned with LH-PCR primer sequences, and the theoretical lengths of the amplicons were determined. Results revealed that the community profiles were significantly different between BBD and SML samples. The SML samples were dominated by gamma-proteobacteria (53 to 64%), followed by beta-proteobacteria (18 to 21%) and alpha-proteobacteria (5 to 11%). In contrast, both BBD clone libraries were dominated by alpha-proteobacteria (58 to 87%), followed by verrucomicrobia (2 to 10%) and 0 to 6% each of delta-proteobacteria, bacteroidetes, firmicutes, and cyanobacteria. Alphaproteobacterial sequence types related to the bacteria associated with toxin-producing dinoflagellates were observed in BBD clone libraries but were not found in the SML libraries. Similarly, sequences affiliated with the family Desulfobacteraceae and toxin-producing cyanobacteria, both believed to be involved in BBD pathogenesis, were found only in BBD libraries. These data provide evidence for an association of numerous toxin-producing heterotrophic microorganisms with BBD of corals.     

Comparison of bacterial composition between human saliva and dental unit water system

The bacterial compositions between the dental unit water system and human saliva were characterized and compared by direct sequence analysis of 16S rDNA clone libraries. Based on the species richness estimation, bacterial diversity in the dental unit water system (DUW) was more diverse than that of the human saliva (HS). The Chao1 estimates of species richness in HS and DUW samples were 12.0 and 72.4, respectively. The total numbers of OTUs observed in the combined libraries accounted for 83% (HS) and 59% (DUW) of the Chao1 diversity estimate as defined at the 80% similarity threshold. Based on the sequence analysis, the phylum Proteobacteria was the major group in both clone libraries at phylum level. DUW clone library contained 80.0% Proteobacteria, 8.0% Bacteroides, 4.0% Nitrospira, 4.0% Firmicutes, 2.0% Planctomycetes and 2.0% Acidobacteria. On the other hand, human saliva (HS) clone library contained 55.5% Proteobacteria, 36.1% Firmicutes and 8.4% Bacteroides. The majority of bacteria identified belonged to phylum Proteobacteria in both samples. In dental unit water system (DUW), Alphaproteobacteria was detected as the major group. There was no evidence of the bacterial contamination due to a dental treatment. Most sequences were related to microorganisms derived from biofilm in oligotrophic environments.     

PCR biases distort bacterial and archaeal community structure in pyrosequencing datasets

As 16S rRNA gene targeted massively parallel sequencing has become a common tool for microbial diversity investigations, numerous advances have been made to minimize the influence of sequencing and chimeric PCR artifacts through rigorous quality control measures. However, there has been little effort towards understanding the effect of multi-template PCR biases on microbial community structure. In this study, we used three bacterial and three archaeal mock communities consisting of, respectively, 33 bacterial and 24 archaeal 16S rRNA gene sequences combined in different proportions to compare the influences of (1) sequencing depth, (2) sequencing artifacts (sequencing errors and chimeric PCR artifacts), and (3) biases in multi-template PCR, towards the interpretation of community structure in pyrosequencing datasets. We also assessed the influence of each of these three variables on α- and β-diversity metrics that rely on the number of OTUs alone (richness) and those that include both membership and the relative abundance of detected OTUs (diversity). As part of this study, we redesigned bacterial and archaeal primer sets that target the V3-V5 region of the 16S rRNA gene, along with multiplexing barcodes, to permit simultaneous sequencing of PCR products from the two domains. We conclude that the benefits of deeper sequencing efforts extend beyond greater OTU detection and result in higher precision in β-diversity analyses by reducing the variability between replicate libraries, despite the presence of more sequencing artifacts. Additionally, spurious OTUs resulting from sequencing errors have a significant impact on richness or shared-richness based α- and β-diversity metrics, whereas metrics that utilize community structure (including both richness and relative abundance of OTUs) are minimally affected by spurious OTUs. However, the greatest obstacle towards accurately evaluating community structure are the errors in estimated mean relative abundance of each detected OTU due to biases associated with multi-template PCR reactions.     

A Next Generation Semiconductor Based Sequencing Approach for the Identification of Meat Species in DNA Mixtures

The identification of the species of origin of meat and meat products is an important issue to prevent and detect frauds that might have economic, ethical and health implications. In this paper we evaluated the potential of the next generation semiconductor based sequencing technology (Ion Torrent Personal Genome Machine) for the identification of DNA from meat species (pig, horse, cattle, sheep, rabbit, chicken, turkey, pheasant, duck, goose and pigeon) as well as from human and rat in DNA mixtures through the sequencing of PCR products obtained from different couples of universal primers that amplify 12S and 16S rRNA mitochondrial DNA genes. Six libraries were produced including PCR products obtained separately from 13 species or from DNA mixtures containing DNA from all species or only avian or only mammalian species at equimolar concentration or at 1:10 or 1:50 ratios for pig and horse DNA. Sequencing obtained a total of 33,294,511 called nucleotides of which 29,109,688 with Q20 (87.43%) in a total of 215,944 reads. Different alignment algorithms were used to assign the species based on sequence data. Error rate calculated after confirmation of the obtained sequences by Sanger sequencing ranged from 0.0003 to 0.02 for the different species. Correlation about the number of reads per species between different libraries was high for mammalian species (0.97) and lower for avian species (0.70). PCR competition limited the efficiency of amplification and sequencing for avian species for some primer pairs. Detection of low level of pig and horse DNA was possible with reads obtained from different primer pairs. The sequencing of the products obtained from different universal PCR primers could be a useful strategy to overcome potential problems of amplification. Based on these results, the Ion Torrent technology can be applied for the identification of meat species in DNA mixtures.     

Characterization of monooxygenase gene diversity in benzene‐amended soils

Aim: To understand soil benzene monooxygenase gene diversity by clone library construction and microarray profiling. Methods and Results: A primer set was designed, and benzene monooxygenase gene diversity was characterized in two benzene‐amended soils. The dominant sequence types in the clone libraries were distinct between the two soils, and both sequences were assigned to novel clusters. Monooxygenase gene richness and diversity increased after benzene degradation. Oligonucleotide probes for microarray analysis were designed to detect a number of sequenced clones and reported monooxygenase genes. The microarray detected several genes that were not detected in the clone libraries of the same samples. Six probes were detected in more than one soil. Conclusions: The primer set designed in this study successfully detected diverse benzene monooxygenase genes. The level of diversity may have increased because the degradation of benzene differed from soil to soil. Microarrays have great potential in the comprehensive detection of gene richness as well as the elucidation of key genes for degradation. Significance and Impact of the Study: This study introduces a new primer set that may be used to identify diverse benzene monooxygenase genes in the environment; moreover, it demonstrates the potential of microarray technology in the profiling of environmental samples.     

Synthetic Statistical Approach Reveals a High Degree of Richness of Microbial Eukaryotes in an Anoxic Water Column

Molecular surveys suggest that communities of microbial eukaryotes are remarkably rich, because even large clone libraries seem to capture only a minority of species. This provides a qualitative picture of protistan richness but does not measure its real extent either locally or globally. Statistical analysis can estimate a community's richness, but the specific methods used to date are not always well grounded in statistical theory. Here we study a large protistan molecular survey from an anoxic water column in the Cariaco Basin (Caribbean Sea). We group individual 18S rRNA gene sequences into operational taxonomic units (OTUs) using different cutoff values for sequence similarity (99 to 50%) and systematically apply parametric models and nonparametric estimators to the OTU frequency data to estimate the total protistan diversity. The parametric models provided statistically sound estimates of protistan richness, with biologically meaningful standard errors, maximal data usage, and extensive model diagnostics and were preferable to the available nonparametric tools. Our clone library exceeded 700 clones but still covered only a minority of species and less than half of the larger protistan clades. Our estimates of total protistan richness portray the target community as very rich at all OTU levels, with hundreds of different populations apparently co-occurring in the small (3-liter) volume of our sample, as well as dozens of clades of the highest taxonomic order. These estimates are among the first for microbial eukaryotes that are obtained using state-of-the-art statistical methods and can serve as benchmark numbers for the local diversity of protists.     

Synthetic statistical approach reveals a high degree of richness of microbial eukaryotes in an anoxic water column

Molecular surveys suggest that communities of microbial eukaryotes are remarkably rich, because even large clone libraries seem to capture only a minority of species. This provides a qualitative picture of protistan richness but does not measure its real extent either locally or globally. Statistical analysis can estimate a community's richness, but the specific methods used to date are not always well grounded in statistical theory. Here we study a large protistan molecular survey from an anoxic water column in the Cariaco Basin (Caribbean Sea). We group individual 18S rRNA gene sequences into operational taxonomic units (OTUs) using different cutoff values for sequence similarity (99 to 50%) and systematically apply parametric models and nonparametric estimators to the OTU frequency data to estimate the total protistan diversity. The parametric models provided statistically sound estimates of protistan richness, with biologically meaningful standard errors, maximal data usage, and extensive model diagnostics and were preferable to the available nonparametric tools. Our clone library exceeded 700 clones but still covered only a minority of species and less than half of the larger protistan clades. Our estimates of total protistan richness portray the target community as very rich at all OTU levels, with hundreds of different populations apparently co-occurring in the small (3-liter) volume of our sample, as well as dozens of clades of the highest taxonomic order. These estimates are among the first for microbial eukaryotes that are obtained using state-of-the-art statistical methods and can serve as benchmark numbers for the local diversity of protists.     

Microbial Communities in the Surface Mucopolysaccharide Layer and the Black Band Microbial Mat of Black Band-Diseased Siderastrea siderea

Microbial community profiles and species composition associated with two black band-diseased colonies of the coral Siderastrea siderea were studied by 16S rRNA-targeted gene cloning, sequencing, and amplicon-length heterogeneity PCR (LH-PCR). Bacterial communities associated with the surface mucopolysaccharide layer (SML) of apparently healthy tissues of the infected colonies, together with samples of the black band disease (BBD) infections, were analyzed using the same techniques for comparison. Gene sequences, ranging from 424 to 1,537 bp, were retrieved from all positive clones (n = 43 to 48) in each of the four clone libraries generated and used for comparative sequence analysis. In addition to LH-PCR community profiling, all of the clone sequences were aligned with LH-PCR primer sequences, and the theoretical lengths of the amplicons were determined. Results revealed that the community profiles were significantly different between BBD and SML samples. The SML samples were dominated by γ-proteobacteria (53 to 64%), followed by β-proteobacteria (18 to 21%) and α-proteobacteria (5 to 11%). In contrast, both BBD clone libraries were dominated by α-proteobacteria (58 to 87%), followed by verrucomicrobia (2 to 10%) and 0 to 6% each of δ-proteobacteria, bacteroidetes, firmicutes, and cyanobacteria. Alphaproteobacterial sequence types related to the bacteria associated with toxin-producing dinoflagellates were observed in BBD clone libraries but were not found in the SML libraries. Similarly, sequences affiliated with the family Desulfobacteraceae and toxin-producing cyanobacteria, both believed to be involved in BBD pathogenesis, were found only in BBD libraries. These data provide evidence for an association of numerous toxin-producing heterotrophic microorganisms with BBD of corals.     

Generation of multimillion-sequence 16S rRNA gene libraries from complex microbial communities by assembling paired-end illumina reads

Microbial communities host unparalleled taxonomic diversity. Adequate characterization of environmental and host-associated samples remains a challenge for microbiologists, despite the advent of 16S rRNA gene sequencing. In order to increase the depth of sampling for diverse bacterial communities, we developed a method for sequencing and assembling millions of paired-end reads from the 16S rRNA gene (spanning the V3 region; ～200 nucleotides) by using an Illumina genome analyzer. To confirm reproducibility and to identify a suitable computational pipeline for data analysis, sequence libraries were prepared in duplicate for both a defined mixture of DNAs from known cultured bacterial isolates (>1 million postassembly sequences) and an Arctic tundra soil sample (>6 million postassembly sequences). The Illumina 16S rRNA gene libraries represent a substantial increase in number of sequences over all extant next-generation sequencing approaches (e.g., 454 pyrosequencing), while the assembly of paired-end 125-base reads offers a methodological advantage by incorporating an initial quality control step for each 16S rRNA gene sequence. This method incorporates indexed primers to enable the characterization of multiple microbial communities in a single flow cell lane, may be modified readily to target other variable regions or genes, and demonstrates unprecedented and economical access to DNAs from organisms that exist at low relative abundances.     

High-Density Universal 16S rRNA Microarray Analysis Reveals Broader Diversity than Typical Clone Library When Sampling the Environment

Molecular approaches aimed at detection of a broad-range of prokaryotes in the environment routinely rely on classifying heterogeneous 16S rRNA genes amplified by polymerase chain reaction (PCR) using primers with broad specificity. The general method of sampling and categorizing DNA has been to clone then sequence the PCR products. However, the number of clones required to adequately catalog the majority of taxa in a sample is unwieldy. Alternatively, hybridizing target sequences to a universal 16S rRNA gene microarray may provide a more rapid and comprehensive view of prokaryotic community composition. This study investigated the breadth and accuracy of a microarray in detecting diverse 16S rRNA gene sequence types compared to clone-and-sequencing using three environmental samples: urban aerosol, subsurface soil, and subsurface water. PCR products generated from universal 16S rRNA gene-targeted primers were classified by using either the clone-and-sequence method or by hybridization to a novel high-density microarray of 297,851 probes complementary to 842 prokaryotic subfamilies. The three clone libraries comprised 1391 high-quality sequences. Approximately 8% of the clones could not be placed into a known subfamily and were considered novel. The microarray results confirmed the majority of clone-detected subfamilies and additionally demonstrated greater amplicon diversity extending into phyla not observed by the cloning method. Sequences matching operational taxonomic units within the phyla Nitrospira, Planctomycetes, and TM7, which were uniquely detected by the array, were verified with specific primers and subsequent amplicon sequencing. Subfamily richness detected by the array corresponded well with nonparametric richness predictions extrapolated from clone libraries except in the water community where clone-based richness predictions were greatly exceeded. It was concluded that although the microarray is unreliable in identifying novel prokaryotic taxa, it reveals greater diversity in environmental samples than sequencing a typically sized clone library. Furthermore, the microarray allowed samples to be rapidly evaluated with replication, a significant advantage in studies of microbial ecology.     

Probe-based negative selection for underrepresented phylotypes in large environmental clone libraries

Studies based on cloning and sequencing to investigate microbial diversity in a vast range of samples has become widespread in recent years. Results have revealed immense microbial diversity in many different environments, but also dominance of a few sequence types in the constructed clone libraries. Here we describe a method to enrich the clone libraries by avoiding sequencing of known, abundant sequence types, instead focusing on novel, rare ones. The protocol is based on gridding the PCR products from clone libraries on membranes and hybridisation of species-specific probes. Clones that do not give positive hybridisation results are sequenced. This method was used for fungal clone libraries from compost samples. Altogether 1536 clones were gridded and six probes used. From these clones, 59% hybridised with a probe, and therefore, only 41% of the clones were sequenced. In addition, 384 samples were sequenced to verify the hybridisation results. The numbers of false-negative (5.2%) and false-positive (3.9%) hybridisations were low. This method provides a mean of lowering the costs of sequencing projects and speeding up the process of characterising microbial diversity in environmental samples. The method is especially suitable for samples with a few dominating sequence types.     

Target region selection is a critical determinant of community fingerprints generated by 16S pyrosequencing

Pyrosequencing of 16S rRNA genes allows for in-depth characterization of complex microbial communities. Although it is known that primer selection can influence the profile of a community generated by sequencing, the extent and severity of this bias on deep-sequencing methodologies is not well elucidated. We tested the hypothesis that the hypervariable region targeted for sequencing and primer degeneracy play important roles in influencing the composition of 16S pyrotag communities. Subgingival plaque from deep sites of current smokers with chronic periodontitis was analyzed using Sanger sequencing and pyrosequencing using 4 primer pairs. Greater numbers of species were detected by pyrosequencing than by Sanger sequencing. Rare taxa constituted nearly 6% of each pyrotag community and less than 1% of the Sanger sequencing community. However, the different target regions selected for pyrosequencing did not demonstrate a significant difference in the number of rare and abundant taxa detected. The genera Prevotella, Fusobacterium, Streptococcus, Granulicatella, Bacteroides, Porphyromonas and Treponema were abundant when the V1-V3 region was targeted, while Streptococcus, Treponema, Prevotella, Eubacterium, Porphyromonas, Campylobacter and Enterococcus predominated in the community generated by V4-V6 primers, and the most numerous genera in the V7-V9 community were Veillonella, Streptococcus, Eubacterium, Enterococcus, Treponema, Catonella and Selenomonas. Targeting the V4-V6 region failed to detect the genus Fusobacterium, while the taxa Selenomonas, TM7 and Mycoplasma were not detected by the V7-V9 primer pairs. The communities generated by degenerate and non-degenerate primers did not demonstrate significant differences. Averaging the community fingerprints generated by V1-V3 and V7-V9 primers provided results similar to Sanger sequencing, while allowing a significantly greater depth of coverage than is possible with Sanger sequencing. It is therefore important to use primers targeted to these two regions of the 16S rRNA gene in all deep-sequencing efforts to obtain representational characterization of complex microbial communities.     

Statistical approaches for estimating actinobacterial diversity in marine sediments

Bacterial diversity in a deep-sea sediment was investigated by constructing actinobacterium-specific 16S ribosomal DNA (rDNA) clone libraries from sediment sections taken 5 to 12, 15 to 18, and 43 to 46 cm below the sea floor at a depth of 3,814 m. Clones were placed into operational taxonomic unit (OTU) groups with >/= 99% 16S rDNA sequence similarity; the cutoff value for an OTU was derived by comparing 16S rRNA homology with DNA-DNA reassociation values for members of the class Actinobacteria. Diversity statistics were used to determine how the level of dominance, species richness, and genetic diversity varied with sediment depth. The reciprocal of Simpson's index (1/D) indicated that the pattern of diversity shifted toward dominance from uniformity with increasing sediment depth. Nonparametric estimation of the species richness in the 5- to 12-, 15- to 18-, and 43- to 46-cm sediment sections revealed a trend of decreasing species number with depth, 1,406, 308, and 212 OTUs, respectively. Application of the LIBSHUFF program indicated that the 5- to 12-cm clone library was composed of OTUs significantly (P = 0.001) different from those of the 15- to 18- and 43- to 46-cm libraries. F(ST) and phylogenetic grouping of taxa (P tests) were both significant (P < 0.00001 and P < 0.001, respectively), indicating that genetic diversity decreased with sediment depth and that each sediment community harbored unique phylogenetic lineages. It was also shown that even nonconservative OTU definitions result in severe underestimation of species richness; unique phylogenetic clades detected in one OTU group suggest that OTUs do not correspond to real ecological groups sensu Palys (T. Palys, L. K. Nakamura, and F. M. Cohan, Int. J. Syst. Bacteriol. 47:1145-1156, 1997). Mechanisms responsible for diversity and their implications are discussed.