Molecular Detection, Quantification, and Diversity Evaluation of Microalgae

This study reviews the available molecular methods and new high-throughput technologies for their practical use in the molecular detection, quantification, and diversity assessment of microalgae. Molecular methods applied to other groups of organisms can be adopted for microalgal studies because they generally detect universal biomolecules, such as nucleic acids or proteins. These methods are primarily related to species detection and discrimination among various microalgae. Among current molecular methods, some molecular tools are highly valuable for small-scale detection [e.g., single-cell polymerase chain reaction (PCR), quantitative real-time PCR (qPCR), and biosensors], whereas others are more useful for large-scale, high-throughput detection [e.g., terminal restriction length polymorphism, isothermal nucleic acid sequence-based amplification, loop-mediated isothermal amplification, microarray, and next generation sequencing (NGS) techniques]. Each molecular technique has its own strengths in detecting microalgae, but they may sometimes have limitations in terms of detection of other organisms. Among current technologies, qPCR may be considered the best method for molecular quantification of microalgae. Metagenomic microalgal diversity can easily be achieved by 454 pyrosequencing rather than by the clone library method. Current NGS, third and fourth generation technologies pave the way for the high-throughput detection and quantification of microalgal diversity, and have significant potential for future use in field monitoring.     

Analysis of High-Throughput Sequencing and Annotation Strategies for Phage Genomes

Background

Bacterial viruses (phages) play a critical role in shaping microbial populations as they influence both host mortality and horizontal gene transfer. As such, they have a significant impact on local and global ecosystem function and human health. Despite their importance, little is known about the genomic diversity harbored in phages, as methods to capture complete phage genomes have been hampered by the lack of knowledge about the target genomes, and difficulties in generating sufficient quantities of genomic DNA for sequencing. Of the approximately 550 phage genomes currently available in the public domain, fewer than 5% are marine phage.

Methodology/Principal Findings

To advance the study of phage biology through comparative genomic approaches we used marine cyanophage as a model system. We compared DNA preparation methodologies (DNA extraction directly from either phage lysates or CsCl purified phage particles), and sequencing strategies that utilize either Sanger sequencing of a linker amplification shotgun library (LASL) or of a whole genome shotgun library (WGSL), or 454 pyrosequencing methods. We demonstrate that genomic DNA sample preparation directly from a phage lysate, combined with 454 pyrosequencing, is best suited for phage genome sequencing at scale, as this method is capable of capturing complete continuous genomes with high accuracy. In addition, we describe an automated annotation informatics pipeline that delivers high-quality annotation and yields few false positives and negatives in ORF calling.

Conclusions/Significance

These DNA preparation, sequencing and annotation strategies enable a high-throughput approach to the burgeoning field of phage genomics.     

Gene capture prediction and overlap estimation in EST sequencing from one or multiple libraries

Background  

In expressed sequence tag (EST) sequencing, we are often interested in how many genes we can capture in an EST sample of a targeted size. This information provides insights to sequencing efficiency in experimental design, as well as clues to the diversity of expressed genes in the tissue from which the library was constructed.     

Comparative analysis of eukaryotic marine microbial assemblages from 18S rRNA gene and gene transcript clone libraries by using different methods of extraction

Eukaryotic marine microbes play pivotal roles in biogeochemical nutrient cycling and ecosystem function, but studies that focus on the protistan biogeography and genetic diversity lag-behind studies of other microbes. 18S rRNA PCR amplification and clone library sequencing are commonly used to assess diversity that is culture independent. However, molecular methods are not without potential biases and artifacts. In this study, we compare the community composition of clone libraries generated from the same water sample collected at the San Pedro Ocean Time Series (SPOTs) station in the northwest Pacific Ocean. Community composition was assessed using different cell lysis methods (chemical and mechanical) and the extraction of different nucleic acids (DNA and RNA reverse transcribed to cDNA) to build Sanger ABI clone libraries. We describe specific biases for ecologically important phylogenetic groups resulting from differences in nucleic acid extraction methods that will inform future designs of eukaryotic diversity studies, regardless of the target sequencing platform planned.     

Environmental rRNA inventories miss over half of protistan diversity

Background  

The main tool to discover novel microbial eukaryotes is the rRNA approach. This approach has important biases, including PCR discrimination against certain rRNA gene species, which makes molecular inventories skewed relative to the source communities. The degree of this bias has not been quantified, and it remains unclear whether species missed from clone libraries could be recovered by increasing sequencing efforts, or whether they cannot be detected in principle. Here we attempt to discriminate between these possibilities by statistically analysing four protistan inventories obtained using different general eukaryotic PCR primers.     

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.     

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.     

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.     

Removing Noise From Pyrosequenced Amplicons

Background  

In many environmental genomics applications a homologous region of DNA from a diverse sample is first amplified by PCR and then sequenced. The next generation sequencing technology, 454 pyrosequencing, has allowed much larger read numbers from PCR amplicons than ever before. This has revolutionised the study of microbial diversity as it is now possible to sequence a substantial fraction of the 16S rRNA genes in a community. However, there is a growing realisation that because of the large read numbers and the lack of consensus sequences it is vital to distinguish noise from true sequence diversity in this data. Otherwise this leads to inflated estimates of the number of types or operational taxonomic units (OTUs) present. Three sources of error are important: sequencing error, PCR single base substitutions and PCR chimeras. We present AmpliconNoise, a development of the PyroNoise algorithm that is capable of separately removing 454 sequencing errors and PCR single base errors. We also introduce a novel chimera removal program, Perseus, that exploits the sequence abundances associated with pyrosequencing data. We use data sets where samples of known diversity have been amplified and sequenced to quantify the effect of each of the sources of error on OTU inflation and to validate these algorithms.     

Characterization of protistan plankton diversity in ancient salt evaporation ponds located in a volcanic crater on the island Sal,Cape Verde

Salinity is an important factor when exploring the limits known for life. Therefore, hypersaline systems have attracted much attention in recent years. In this study, we investigated the protistan diversity and community composition in two natural salt evaporation ponds (27–30% salinity) located in an ancient volcanic crater on the Cape Verde island Sal using high-throughput DNA sequencing. Our study revealed a broad range of protistan taxa and a high taxonomic diversity within the Ciliophora, Dinophyceae, and Chlorophyta. We detected a total of 23 Dinophyceae families, although Dinophyceae were generally considered to be only this diverse in aquatic environments of less than 10% salinity. Moreover, we uncovered a high degree of genetic novelty in this habitat. The mean similarity of all detected OTUs to previously described sequences was only 93.6%. These findings strongly dispute the traditional view that extreme hypersaline environments generally maintain low protistan diversity. A meta-analysis covering our and previously published data from other inland and coastal salt ponds clearly showed that our samples clustered according to salinity and not biogeography. This result further supports the claim that salinity is a major transition boundary for protistan communities, regardless of their biogeographic origin.

     

Accurate multiplexing and filtering for high-throughput amplicon-sequencing

Tagging amplicons with tag sequences appended to PCR primers allow the multiplexing of numerous samples for high-throughput sequencing (HTS). This approach is routinely used in HTS-based diversity analyses, especially in microbial ecology and biomedical diagnostics. However, amplicon library preparation is subject to pervasive sample sequence cross-contaminations as a result of tag switching events referred to as mistagging. Here, we sequenced seven amplicon libraries prepared using various multiplexing designs in order to measure the magnitude of this phenomenon and its impact on diversity analyses. Up to 28.2% of the unique sequences correspond to undetectable (critical) mistags in single- or saturated double-tagging libraries. We show the advantage of multiplexing samples following Latin Square Designs in order to optimize the detection of mistags and maximize the information on their distribution across samples. We use this information in designs incorporating PCR replicates to filter the critical mistags and to recover the exact composition of mock community samples. Being parameter-free and data-driven, our approach can provide more accurate and reproducible HTS data sets, improving the reliability of their interpretations.     

Extension of Lander-Waterman theory for sequencing filtered DNA libraries

Background  

The degree to which conventional DNA sequencing techniques will be successful for highly repetitive genomes is unclear. Investigators are therefore considering various filtering methods to select against high-copy sequence in DNA clone libraries. The standard model for random sequencing, Lander-Waterman theory, does not account for two important issues in such libraries, discontinuities and position-based sampling biases (the so-called "edge effect"). We report an extension of the theory for analyzing such configurations.     

Unveiling trophic functions of uncultured protist taxa by incubation experiments in the brackish baltic sea

Background

Our knowledge of the phylogeny and diversity of aquatic protists is rapidly increasing due to molecular surveys and next-generation sequencing approaches. This has led to a considerable discrepancy between the taxa known from cultures and those known from environmental 18S rRNA gene sequences. Hence, it is generally difficult to assign ecological functions to new taxa detected by culture-independent molecular approaches.

Methodology/Principal Findings

A combination of unamended dark incubations and 18S rRNA sequencing was chosen to link molecular diversity data of uncultured protists with heterotrophic, presumably bacterivorous, growth. The incubations, conducted with Baltic Sea brackish water, resulted in a consistent shift from a protistan community dominated by phototrophs to one in which heterotrophs predominated. This was determined on the basis of cell abundance and 18S rRNA sequences derived from fingerprint analysis and clone libraries. The bulk of enriched phylotypes after incubation were related to hitherto uncultured marine taxa within chrysophytes, ochrophytes, choanoflagellates, cercozoans, and picobiliphytes, mostly represented in recently established or here defined environmental clades. Their growth in the dark, together with coinciding results from studies with a similar objective, provides evidence that these uncultured taxa represent heterotrophic or mixotrophic species.

Conclusions/Significance

These findings shed some light into the trophic role of diverse uncultured protists especially within functionally heterogeneous groups (e.g., chrysophytes, ochrophytes) and groups that appear to be puzzling with regard to their nutrition (picobiliphytes). Additionally, our results indicate that the heterotrophic flagellate community in the southwestern Baltic Sea is dominated by species of marine origin. The combination of unamended incubations with molecular diversity analysis is thus confirmed as a promising approach to explore the trophic mode of environmentally relevant protist taxa for which only sequence data are currently available.     

Evaluation of genomic high-throughput sequencing data generated on Illumina HiSeq and Genome Analyzer systems

Background  

The generation and analysis of high-throughput sequencing data are becoming a major component of many studies in molecular biology and medical research. Illumina's Genome Analyzer (GA) and HiSeq instruments are currently the most widely used sequencing devices. Here, we comprehensively evaluate properties of genomic HiSeq and GAIIx data derived from two plant genomes and one virus, with read lengths of 95 to 150 bases.     

Systematic cross-validation of 454 sequencing and pyrosequencing for the exact quantification of DNA methylation patterns with single CpG resolution

Background  

New high-throughput sequencing technologies promise a very sensitive and high-resolution analysis of DNA methylation patterns in quantitative terms. However, a detailed and comprehensive comparison with existing validated DNA methylation analysis methods is not yet available. Therefore, a systematic cross-validation of 454 sequencing and conventional pyrosequencing, both of which offer exact quantification of methylation levels with a single CpG dinucleotide resolution, was performed.     

Multivariate Cutoff Level Analysis (MultiCoLA) of large community data sets

High-throughput sequencing techniques are becoming attractive to molecular biologists and ecologists as they provide a time- and cost-effective way to explore diversity patterns in environmental samples at an unprecedented resolution. An issue common to many studies is the definition of what fractions of a data set should be considered as rare or dominant. Yet this question has neither been satisfactorily addressed, nor is the impact of such definition on data set structure and interpretation been fully evaluated. Here we propose a strategy, MultiCoLA (Multivariate Cutoff Level Analysis), to systematically assess the impact of various abundance or rarity cutoff levels on the resulting data set structure and on the consistency of the further ecological interpretation. We applied MultiCoLA to a 454 massively parallel tag sequencing data set of V6 ribosomal sequences from marine microbes in temperate coastal sands. Consistent ecological patterns were maintained after removing up to 35–40% rare sequences and similar patterns of beta diversity were observed after denoising the data set by using a preclustering algorithm of 454 flowgrams. This example validates the importance of exploring the impact of the definition of rarity in large community data sets. Future applications can be foreseen for data sets from different types of habitats, e.g. other marine environments, soil and human microbiota.     

Distinct protistan assemblages characterize the euphotic zone and deep sea (2500 m) of the western North Atlantic (Sargasso Sea and Gulf Stream)

Protistan diversity was characterized at three locations in the western North Atlantic (Sargasso Sea and Gulf Stream) by sequencing 18S rRNA genes in samples from euphotic (< or = 125 m) and bathypelagic depths (2500 m). A total of 923 partial-length protistan sequences were analysed, revealing 324 distinct operational taxonomic units (OTUs) determined by an automated OTU-calling program set to 95% sequence similarity. Most OTUs were comprised of only one or two sequences suggesting a large but rare pool of protistan diversity. Many OTUs from both depth strata were associated with recently described novel alveolate and stramenopile lineages while many OTUs from the bathypelagic were affiliated with Acantharea, Polycystinea and Euglenozoa and were not observed in euphotic zone libraries. Protistan assemblages from the euphotic zone and the deep sea were largely composed of distinct OTUs; only 28 of the 324 protistan OTUs were detected in both shallow and deep sea clone libraries. The diversity of protistan assemblages in the deep sea was distinctly lower than the diversity of euphotic zone assemblages. Protistan assemblages from the Gulf Stream were the most diverse for either depth strata. Overall, protistan assemblages from different stations but comparable depths were more similar than the assemblages from different depths at the same station. These data suggest that particular groups of protistan OTUs formed distinct 'shallow' and 'deep-sea' assemblages across widely spaced oceanic locales.     

Graph-based clustering and characterization of repetitive sequences in next-generation sequencing data

Background  

The investigation of plant genome structure and evolution requires comprehensive characterization of repetitive sequences that make up the majority of higher plant nuclear DNA. Since genome-wide characterization of repetitive elements is complicated by their high abundance and diversity, novel approaches based on massively-parallel sequencing are being adapted to facilitate the analysis. It has recently been demonstrated that the low-pass genome sequencing provided by a single 454 sequencing reaction is sufficient to capture information about all major repeat families, thus providing the opportunity for efficient repeat investigation in a wide range of species. However, the development of appropriate data mining tools is required in order to fully utilize this sequencing data for repeat characterization.     

Rapid shifts in the structure and composition of a protistan assemblage during bottle incubations affect estimates of total protistan species richness

Changes in the structure and composition of a protistan community were characterized through the analysis of small-subunit ribosomal RNA gene (18S) sequences for a 3-day bottle incubation using a single sample collected in the western North Atlantic. Cloning and sequencing was used to investigate changes in perceived species richness and diversity as a consequence of environmental perturbation. The treatments included a control (unamended seawater), inorganic nutrient enrichment, and enrichment with a complex organic mixture. Five clone libraries were constructed and analyzed at the time of collection (t-0?h) and after 24 (t-24?h) and 72 (t-72?h) h for the control, and at t-72?h for the inorganic and organic enrichments, resulting in an analysis of 1,626 partial 18S rDNA sequences that clustered into 238 operational taxonomic units (OTUs). Analysis of the clone libraries revealed that protistan assemblages were highly dynamic and changed substantially at both the OTU level and higher taxonomic classifications during time frames consistent with many oceanographic methods used for measuring biological rates. Changes were most dramatic in enrichments, which yielded community compositions that were strongly dominated by one or a few taxa. Changes in community structure during incubation dramatically influenced estimates of species richness, which were substantially lower with longer incubation and especially with amendment, even though all incubated samples originated from the same aliquot of seawater. Containment and enrichment of the seawater sample led to the detection of otherwise undetected protistan taxa, suggesting that characterization of protistan diversity in a sample only at the time of collection could lead to an underrepresentation of unique taxa. Additionally, the rapid increase in the relative abundance of some members of the ??rare biosphere?? in our results implies an ecological importance of at least some of the taxa comprising the ??rare biosphere.??     

Sequencing of 6.7 Mb of the melon genome using a BAC pooling strategy

Background  

Cucumis melo (melon) belongs to the Cucurbitaceae family, whose economic importance among horticulture crops is second only to Solanaceae. Melon has a high intra-specific genetic variation, morphologic diversity and a small genome size (454 Mb), which make it suitable for a great variety of molecular and genetic studies. A number of genetic and genomic resources have already been developed, such as several genetic maps, BAC genomic libraries, a BAC-based physical map and EST collections. Sequence information would be invaluable to complete the picture of the melon genomic landscape, furthering our understanding of this species' evolution from its relatives and providing an important genetic tool. However, to this day there is little sequence data available, only a few melon genes and genomic regions are deposited in public databases. The development of massively parallel sequencing methods allows envisaging new strategies to obtain long fragments of genomic sequence at higher speed and lower cost than previous Sanger-based methods.