Hypervariable regions (V1–V9) of the dinoflagellate 18S rRNA using a large dataset for marker considerations

DNA sequencing methods have been used for the molecular taxonomic discrimination of dinoflagellate protists, particularly using partial 18S rRNA sequences. This study evaluated the taxonomic discrimination power of rRNA gene hypervariable regions (V1 to V9) in dinoflagellates from a large dataset. These included 77 dinoflagellate species (9 orders, 17 families, 40 genera). The complete 18S rRNA sequences of the dinoflagellates ranged from 1,787 to 1,813?bp in length, and consisted of eight V regions with a total combined length of 678 to 699?bp. Regions longer than 100?bp were recoded for V2, V4, and V8 regions; high nucleotide divergences were detected in V1, V2, and V4 regions. Statistic tests showed that the divergences of individual V regions were significantly different (t-test, P?<?0.05) compared with the complete 18S rRNA. The V2 region showed the highest score (83.5%) for PI sites. Moreover, intra-genus DNA similarities of the V2 were considerably low (<93%). Neighbor-joining analyses showed that phylogenetic resolution in the V2–V4 region was 1.32-fold higher than that of the complete 18S rRNA. These results demonstrate that V2 has the highest taxonomic resolving power within the 18S rRNA gene of dinoflagellates, suggesting the V2 and adjacent regions (e.g., V1 to V4) may be the best for marker considerations.     

Investigating Microbial Eukaryotic Diversity from a Global Census: Insights from a Comparison of Pyrotag and Full-Length Sequences of 18S rRNA Genes

Next-generation DNA sequencing (NGS) approaches are rapidly surpassing Sanger sequencing for characterizing the diversity of natural microbial communities. Despite this rapid transition, few comparisons exist between Sanger sequences and the generally much shorter reads of NGS. Operational taxonomic units (OTUs) derived from full-length (Sanger sequencing) and pyrotag (454 sequencing of the V9 hypervariable region) sequences of 18S rRNA genes from 10 global samples were analyzed in order to compare the resulting protistan community structures and species richness. Pyrotag OTUs called at 98% sequence similarity yielded numbers of OTUs that were similar overall to those for full-length sequences when the latter were called at 97% similarity. Singleton OTUs strongly influenced estimates of species richness but not the higher-level taxonomic composition of the community. The pyrotag and full-length sequence data sets had slightly different taxonomic compositions of rhizarians, stramenopiles, cryptophytes, and haptophytes, but the two data sets had similarly high compositions of alveolates. Pyrotag-based OTUs were often derived from sequences that mapped to multiple full-length OTUs at 100% similarity. Thus, pyrotags sequenced from a single hypervariable region might not be appropriate for establishing protistan species-level OTUs. However, nonmetric multidimensional scaling plots constructed with the two data sets yielded similar clusters, indicating that beta diversity analysis results were similar for the Sanger and NGS sequences. Short pyrotag sequences can provide holistic assessments of protistan communities, although care must be taken in interpreting the results. The longer reads (>500 bp) that are now becoming available through NGS should provide powerful tools for assessing the diversity of microbial eukaryotic assemblages.     

SSU‐rRNA Gene Sequencing Survey of Benthic Microbial Eukaryotes from Guaymas Basin Hydrothermal Vent

Microbial eukaryotes have important roles in marine food webs, but their diversity and activities in hydrothermal vent ecosystems are poorly characterized. In this study, we analyzed microbial eukaryotic communities associated with bacterial (Beggiatoa) mats in the 2,000 m deep‐sea Guaymas Basin hydrothermal vent system using 18S rRNA gene high‐throughput sequencing of the V4 region. We detected 6,954 distinct Operational Taxonomic Units (OTUs) across various mat systems. Of the sequences that aligned with known protistan phylotypes, most were affiliated with alveolates (especially dinoflagellates and ciliates) and cercozoans. OTU richness and community structure differed among sediment habitats (e.g. different mat types and cold sediments away from mats). Additionally, full‐length 18S rRNA genes amplified and cloned from single cells revealed the identities of some of the most commonly encountered, active ciliates in this hydrothermal vent ecosystem. Observations and experiments were also conducted to demonstrate that ciliates were trophically active and ingesting fluorescent bacteria or Beggiatoa trichomes. Our work suggests that the active and diverse protistan community at the Guaymas Basin hydrothermal vent ecosystem likely consumes substantial amounts of bacterial biomass, and that the different habitats, often defined by distances of just a few 10s of cm, select for particular assemblages and levels of diversity.     

基于分子标记的宏基因组16S rRNA基因高变区选择策略

随着新一代DNA测序技术出现,人们能够同时对多个DNA样本的宏基因组进行并行分析,尤其是以16S rRNA基因高变区为分子标记的测序已经成为微生物多样性研究最为简洁有效的方法. 目前二代高通量测序的读长不能覆盖16S rRNA基因的全长,需要选择一个有效的高变区进行测序.十多年来,对于16S rRNA基因高变区的选择策略没有统一的标准.本文分析了常用的高变区选择策略,指出不同环境条件是影响高变区选择的重要因素之一.在此基础上,提出了高变区选择的参考准则,同时建议应对选择的高变区进行有效评估.     

Culturing Heterotrophic Protists from the Baltic Sea: Mostly the “Usual Suspects” but a Few Novelties as Well

The study of cultured strains has a long tradition in protistological research and has greatly contributed to establishing the morphology, taxonomy, and ecology of many protist species. However, cultivation‐independent techniques, based on 18S rRNA gene sequences, have demonstrated that natural protistan assemblages mainly consist of hitherto uncultured protist lineages. This mismatch impedes the linkage of environmental diversity data with the biological features of cultured strains. Thus, novel taxa need to be obtained in culture to close this knowledge gap. In this study, traditional cultivation techniques were applied to samples from coastal surface waters and from deep oxygen‐depleted waters of the Baltic Sea. Based on 18S rRNA gene sequencing, 126 monoclonal cultures of heterotrophic protists were identified. The majority of the isolated strains were affiliated with already cultured and described taxa, mainly chrysophytes and bodonids. This was likely due to “culturing bias” but also to the eutrophic nature of the Baltic Sea. Nonetheless, ~ 12% of the isolates in our culture collection showed highly divergent 18S rRNA gene sequences compared to those of known organisms and thus may represent novel taxa, either at the species level or at the genus level. Moreover, we also obtained evidence that some of the isolated taxa are ecologically relevant, under certain conditions, in the Baltic Sea.     

Diversity estimates of microeukaryotes below the chemocline of the anoxic Mariager Fjord, Denmark

Microbial communities of extreme environments have often been assumed to have low species richness. We analysed 18S rRNA gene signatures in a sample collected below the chemocline of the anoxic Mariager Fjord in Denmark, and from these data we computed novel parametric and standard nonparametric estimates of protistan phylotype richness. Our results indicate unexpectedly high richness in this environment: at the 99.5% phylotype definition, our most conservative estimate was 568 phylotypes (+/-114, standard error). Phylogenetic analyses revealed that the sequences collected cover the majority of described lineages in the eukaryotic domain. Out of 384 sequences analysed, 307 were identified as protistan targets, none of which was identical to known sequences. However, based on what is known about species that are phylogenetically related to the Mariager sequences, most of the latter seem to belong to strictly or facultative anaerobe organisms. We also found signatures that together with other environmental 18S rRNA gene sequences represent environmental clades of possibly high taxonomic levels (class to kingdom level). One of these clades, consisting exclusively of sequences from anoxic sampling sites, branches at the base of the eukaryotic evolutionary tree among the earliest eukaryotic lineages. Assuming eukaryotic evolution under oxygen-depleted conditions, these sequences may represent immediate descendants of early eukaryotic ancestors.     

Phylogenetic analysis of the hypervariable region of the 18S rRNA gene of Cryptosporidium oocysts in feces of Canada geese (Branta canadensis): evidence for five novel genotypes

To assess genetic diversity in Cryptosporidium oocysts from Canada geese, 161 fecal samples from Canada geese in the United States were analyzed. Eleven (6.8%) were positive for Cryptosporidium spp. following nested PCR amplification of the hypervariable region of the 18S rRNA gene. Nine PCR products from geese were cloned and sequenced, and all nine diverged from previously reported Cryptosporidium 18S rRNA gene sequences. Five sequences were very similar or identical to each other but genetically distinct from that of Cryptosporidium baileyi; two were most closely related to, but genetically distinct from, the first five; and two were distinct from any other sequence analyzed. One additional sequence in the hypervariable region of the 18S rRNA gene isolated from a cormorant was identical to that of C. baileyi. Phylogenetic analysis provided evidence for new genotypes of Cryptosporidium species in Canada geese. Results of this study suggest that the taxonomy of Cryptosporidium species in geese is complex and that a more complete understanding of genetic diversity among these parasites will facilitate our understanding of oocyst sources and species in the environment.     

Variable copy number, intra-genomic heterogeneities and lateral transfers of the 16S rRNA gene in Pseudomonas

Even though the 16S rRNA gene is the most commonly used taxonomic marker in microbial ecology, its poor resolution is still not fully understood at the intra-genus level. In this work, the number of rRNA gene operons, intra-genomic heterogeneities and lateral transfers were investigated at a fine-scale resolution, throughout the Pseudomonas genus. In addition to nineteen sequenced Pseudomonas strains, we determined the 16S rRNA copy number in four other Pseudomonas strains by Southern hybridization and Pulsed-Field Gel Electrophoresis, and studied the intra-genomic heterogeneities by Denaturing Gradient Gel Electrophoresis and sequencing. Although the variable copy number (from four to seven) seems to be correlated with the evolutionary distance, some close strains in the P. fluorescens lineage showed a different number of 16S rRNA genes, whereas all the strains in the P. aeruginosa lineage displayed the same number of genes (four copies). Further study of the intra-genomic heterogeneities revealed that most of the Pseudomonas strains (15 out of 19 strains) had at least two different 16S rRNA alleles. A great difference (5 or 19 nucleotides, essentially grouped near the V1 hypervariable region) was observed only in two sequenced strains. In one of our strains studied (MFY30 strain), we found a difference of 12 nucleotides (grouped in the V3 hypervariable region) between copies of the 16S rRNA gene. Finally, occurrence of partial lateral transfers of the 16S rRNA gene was further investigated in 1803 full-length sequences of Pseudomonas available in the databases. Remarkably, we found that the two most variable regions (the V1 and V3 hypervariable regions) had probably been laterally transferred from another evolutionary distant Pseudomonas strain for at least 48.3 and 41.6% of the 16S rRNA sequences, respectively. In conclusion, we strongly recommend removing these regions of the 16S rRNA gene during the intra-genus diversity studies.     

Phylogenetic Analysis of the Hypervariable Region of the 18S rRNA Gene of Cryptosporidium Oocysts in Feces of Canada Geese (Branta canadensis): Evidence for Five Novel Genotypes

To assess genetic diversity in Cryptosporidium oocysts from Canada geese, 161 fecal samples from Canada geese in the United States were analyzed. Eleven (6.8%) were positive for Cryptosporidium spp. following nested PCR amplification of the hypervariable region of the 18S rRNA gene. Nine PCR products from geese were cloned and sequenced, and all nine diverged from previously reported Cryptosporidium 18S rRNA gene sequences. Five sequences were very similar or identical to each other but genetically distinct from that of Cryptosporidium baileyi; two were most closely related to, but genetically distinct from, the first five; and two were distinct from any other sequence analyzed. One additional sequence in the hypervariable region of the 18S rRNA gene isolated from a cormorant was identical to that of C. baileyi. Phylogenetic analysis provided evidence for new genotypes of Cryptosporidium species in Canada geese. Results of this study suggest that the taxonomy of Cryptosporidium species in geese is complex and that a more complete understanding of genetic diversity among these parasites will facilitate our understanding of oocyst sources and species in the environment.     

Protistan microbial observatory in the Cariaco Basin,Caribbean. I. Pyrosequencing vs Sanger insights into species richness

Microbial diversity and distribution are topics of intensive research. In two companion papers in this issue, we describe the results of the Cariaco Microbial Observatory (Caribbean Sea, Venezuela). The Basin contains the largest body of marine anoxic water, and presents an opportunity to study protistan communities across biogeochemical gradients. In the first paper, we survey 18S ribosomal RNA (rRNA) gene sequence diversity using both Sanger- and pyrosequencing-based approaches, employing multiple PCR primers, and state-of-the-art statistical analyses to estimate microbial richness missed by the survey. Sampling the Basin at three stations, in two seasons, and at four depths with distinct biogeochemical regimes, we obtained the largest, and arguably the least biased collection of over 6000 nearly full-length protistan rRNA gene sequences from a given oceanographic regime to date, and over 80 000 pyrosequencing tags. These represent all major and many minor protistan taxa, at frequencies globally similar between the two sequence collections. This large data set provided, via the recently developed parametric modeling, the first statistically sound prediction of the total size of protistan richness in a large and varied environment, such as the Cariaco Basin: over 36 000 species, defined as almost full-length 18S rRNA gene sequence clusters sharing over 99% sequence homology. This richness is a small fraction of the grand total of known protists (over 100 000–500 000 species), suggesting a degree of protistan endemism.     

Seasonal diversity of planktonic protists in Southwestern Alberta rivers over a 1-year period as revealed by terminal restriction fragment length polymorphism and 18S rRNA gene library analyses

The temporal dynamics of planktonic protists in river water have received limited attention despite their ecological significance and recent studies linking phagotrophic protists to the persistence of human-pathogenic bacteria. Using molecular-based techniques targeting the 18S rRNA gene, we studied the seasonal diversity of planktonic protists in Southwestern Alberta rivers (Oldman River Basin) over a 1-year period. Nonmetric multidimensional scaling analysis of terminal restriction fragment length polymorphism (T-RFLP) data revealed distinct shifts in protistan community profiles that corresponded to season rather than geographical location. Community structures were examined by using clone library analysis; HaeIII restriction profiles of 18S rRNA gene amplicons were used to remove prevalent solanaceous plant clones prior to sequencing. Sanger sequencing of the V1-to-V3 region of the 18S rRNA gene libraries from spring, summer, fall, and winter supported the T-RFLP results and showed marked seasonal differences in the protistan community structure. The spring library was dominated by Chloroplastidae (29.8%), Centrohelida (28.1%), and Alveolata (25.5%), while the summer and fall libraries contained primarily fungal clones (83.0% and 88.0%, respectively). Alveolata (35.6%), Euglenozoa (24.4%), Chloroplastida (15.6%), and Fungi (15.6%) dominated the winter library. These data demonstrate that planktonic protists, including protozoa, are abundant in river water in Southwestern Alberta and that conspicuous seasonal shifts occur in the community structure.     

Methodology of protistan discovery: from rRNA detection to quality scanning electron microscope images

Each year, thousands of new protistan 18S rRNA sequences are detected in environmental samples. Many of these sequences are molecular signatures of new protistan species, classes, and/or kingdoms that have never been seen before. The main goal of this study was to enable visualization of these novel organisms and to conduct quality ultrastructural examination. We achieved this goal by modifying standard procedures for cell fixation, fluorescence in situ hybridization, and scanning electron microscopy (SEM) and by making these methodologies work in concert. As a result, the same individual cell can now be detected by 18S rRNA-targeted fluorochrome-labeled probes and then viewed by SEM to reveal its diagnostic morphological characteristics. The method was successfully tested on a wide range of protists (alveolates, stramenopiles, kinetoplastids, and cryptomonads). The new methodology thus opens a way for fine microscopy studies of many organisms previously known exclusively by their 18S rRNA sequences.     

Microeukaryotic diversity in marine environments,an analysis of surface layer sediments from the East Sea

Molecular techniques, based on clone library of 18S rRNA gene, were employed to ascertain the diversity of microeukaryotic organisms in sediments from the East Sea. A total of 261 clones were recovered from surface sediments. Most of the clone sequences (90%) were affiliated with protists, dominated by Ciliates (18%) and Dinoflagellates (19%) of Alveolates, phototrophic Stramenopiles (11%), and Cercozoa (20%). Many of the clones were related to uncultivated eukaryotes clones retrieved from anoxic environments with several highly divergent 18S rRNA gene sequences. However, no clones were related to cultivated obligate anaerobic protists. Protistan communities between subsurface layers of 1 and 9 cm shared 23% of total phylotypes which comprised 64% of total clones retrieved. Analysis of diversity indices and rarefaction curve showed that the protistan community within the 1 cm layer exhibited higher diversity than the 9 cm layer. Our results imply that diverse protists remain to be uncovered within marine benthic environments.     

Comparative study of the validity of three regions of the 18S‐rRNA gene for massively parallel sequencing‐based monitoring of the planktonic eukaryote community

The nuclear 18S‐rRNA gene has been used as a metabarcoding marker in massively parallel sequencing (MPS)‐based environmental surveys for plankton biodiversity research. However, different hypervariable regions have been used in different studies, and their utility has been debated among researchers. In this study, detailed investigations into 18S‐rRNA were carried out; we investigated the effective number of sequences deposited in international nucleotide sequence databases (INSDs), the amplification bias, and the amplicon sequence variability among the three variable regions, V1–3, V4–5 and V7–9, using in silico polymerase chain reaction (PCR) amplification based on INSDs. We also examined the primer universality and the taxonomic identification power, using MPS‐based environmental surveys in the Sea of Okhotsk, to determine which region is more useful for MPS‐based monitoring. The primer universality was not significantly different among the three regions, but the number of sequences deposited in INSDs was markedly larger for the V4–5 region than for the other two regions. The sequence variability was significantly different, with the highest variability in the V1–3 region, followed by the V7–9 region, and the lowest variability in the V4–5 region. The results of the MPS‐based environmental surveys showed significantly higher identification power in the V1–3 and V7–9 regions than in the V4–5 region, but no significant difference was detected between the V1–3 and V7–9 regions. We therefore conclude that the V1–3 region will be the most suitable for future MPS‐based monitoring of natural eukaryote communities, as the number of sequences deposited in INSDs increases.     

Application of the hypervariable region of the 16S rDNA sequence as an index for the rapid identification of species in the genus Alicyclobacillus

A comparison of the 16S rRNA gene (rDNA) sequences of seven type strains belonging to different Alicyclobacillus species (i.e., five validated species, one proposed species and one genomic species) suggested that the 5' end hypervariable region (259-273 bases in length) of 16S rDNA was specific for the respective type strains. Further phylogenetic analysis based on DNA sequences of the hypervariable region using 24 Alicyclobacillus strains revealed that the strains could be categorized into five species and the A. acidocaldarius-Alicyclobacillus genomic species 1 group. The hypervariable region was highly conserved among the five species: A. acidiphilus, A. acidoterrestris, A. cycloheptanicus, A. herbarius, and A. hesperidum. The strains in the A. acidocaldarius-Alicyclobacillus genomic species 1 group were subdivided into two clusters (Clusters I and II) based on DNA sequences in the hypervariable region. On the basis of phenotypic characteristics, chemotaxonomic and phylogenetic analyses, and DNA-DNA hybridization data, strains in Cluster I were grouped as Alicyclobacillus genomic species 1 and strains in Cluster II were re-identified as A. acidocaldarius, thereby demonstrating that the hypervariable regions were also highly conserved within these two species. These results suggest that as is the case with Bacillus, the hypervariable region is significantly species-specific in the genus Alicyclobacillus to distinguish Alicyclobacillus species by DNA sequence comparison of the hypervariable region.     

Soil Bacterial Diversity Screening Using Single 16S rRNA Gene V Regions Coupled with Multi-Million Read Generating Sequencing Technologies

The novel multi-million read generating sequencing technologies are very promising for resolving the immense soil 16S rRNA gene bacterial diversity. Yet they have a limited maximum sequence length screening ability, restricting studies in screening DNA stretches of single 16S rRNA gene hypervariable (V) regions. The aim of the present study was to assess the effects of properties of four consecutive V regions (V3-6) on commonly applied analytical methodologies in bacterial ecology studies. Using an in silico approach, the performance of each V region was compared with the complete 16S rRNA gene stretch. We assessed related properties of the soil derived bacterial sequence collection of the Ribosomal Database Project (RDP) database and concomitantly performed simulations based on published datasets. Results indicate that overall the most prominent V region for soil bacterial diversity studies was V3, even though it was outperformed in some of the tests. Despite its high performance during most tests, V4 was less conserved along flanking sites, thus reducing its ability for bacterial diversity coverage. V5 performed well in the non-redundant RDP database based analysis. However V5 did not resemble the full-length 16S rRNA gene sequence results as well as V3 and V4 did when the natural sequence frequency and occurrence approximation was considered in the virtual experiment. Although, the highly conserved flanking sequence regions of V6 provide the ability to amplify partial 16S rRNA gene sequences from very diverse owners, it was demonstrated that V6 was the least informative compared to the rest examined V regions. Our results indicate that environment specific database exploration and theoretical assessment of the experimental approach are strongly suggested in 16S rRNA gene based bacterial diversity studies.     

Organization of ribosomal RNA gene repeats of the mouse.

The organization of the ribosomal RNA (rRNA) genes of the mouse was determined by Southern blot hybridization using cloned rDNA fragments as probes, which could encompass the entire spacer region between two rRNA gene regions. The rRNA genes are organized into tandem repeats of nearly uniform length of about 44 kb. The heterogeneity detected in the nontranscribed spacer appears to be caused by its sequence rather than its length difference. At least three kinds of repetitive sequences are present in the non-transcribed spacer region; two of them are located 13 kb upstream from the 5'-end of 18S RNA gene and the other located 1 to 4 kb downstream from the 3'-end of 28S RNA gene.     

Microeukaryote community patterns along an O2/H2S gradient in a supersulfidic anoxic fjord (Framvaren, Norway)

To resolve the fine-scale architecture of anoxic protistan communities, we conducted a cultivation-independent 18S rRNA survey in the superanoxic Framvaren Fjord in Norway. We generated three clone libraries along the steep O(2)/H(2)S gradient, using the multiple-primer approach. Of 1,100 clones analyzed, 753 proved to be high-quality protistan target sequences. These sequences were grouped into 92 phylotypes, which displayed high protistan diversity in the fjord (17 major eukaryotic phyla). Only a few were closely related to known taxa. Several sequences were dissimilar to all previously described sequences and occupied a basal position in the inferred phylogenies, suggesting that the sequences recovered were derived from novel, deeply divergent eukaryotes. We detected sequence clades with evolutionary importance (for example, clades in the euglenozoa) and clades that seem to be specifically adapted to anoxic environments, challenging the hypothesis that the global dispersal of protists is uniform. Moreover, with the detection of clones affiliated with jakobid flagellates, we present evidence that primitive descendants of early eukaryotes are present in this anoxic environment. To estimate sample coverage and phylotype richness, we used parametric and nonparametric statistical methods. The results show that although our data set is one of the largest published inventories, our sample missed a substantial proportion of the protistan diversity. Nevertheless, statistical and phylogenetic analyses of the three libraries revealed the fine-scale architecture of anoxic protistan communities, which may exhibit adaptation to different environmental conditions along the O(2)/H(2)S gradient.     

Microeukaryote Community Patterns along an O2/H2S Gradient in a Supersulfidic Anoxic Fjord (Framvaren, Norway)

To resolve the fine-scale architecture of anoxic protistan communities, we conducted a cultivation-independent 18S rRNA survey in the superanoxic Framvaren Fjord in Norway. We generated three clone libraries along the steep O₂/H₂S gradient, using the multiple-primer approach. Of 1,100 clones analyzed, 753 proved to be high-quality protistan target sequences. These sequences were grouped into 92 phylotypes, which displayed high protistan diversity in the fjord (17 major eukaryotic phyla). Only a few were closely related to known taxa. Several sequences were dissimilar to all previously described sequences and occupied a basal position in the inferred phylogenies, suggesting that the sequences recovered were derived from novel, deeply divergent eukaryotes. We detected sequence clades with evolutionary importance (for example, clades in the euglenozoa) and clades that seem to be specifically adapted to anoxic environments, challenging the hypothesis that the global dispersal of protists is uniform. Moreover, with the detection of clones affiliated with jakobid flagellates, we present evidence that primitive descendants of early eukaryotes are present in this anoxic environment. To estimate sample coverage and phylotype richness, we used parametric and nonparametric statistical methods. The results show that although our data set is one of the largest published inventories, our sample missed a substantial proportion of the protistan diversity. Nevertheless, statistical and phylogenetic analyses of the three libraries revealed the fine-scale architecture of anoxic protistan communities, which may exhibit adaptation to different environmental conditions along the O₂/H₂S gradient.