Intracellular diversity of the V4 and V9 regions of the 18S rRNA in marine protists (radiolarians) assessed by high-throughput sequencing

Comparative assessment of the intragenomic variations of dinoflagellate Tripos species through single-cell sequencing

Tripos is a large dinoflagellate genus widely distributed in the world's oceans. Morphology-based species identification is inconclusive due to high morphological intraspecific variability. Metabarcoding analysis has been demonstrated to be effective for species identification and tracking their spatiotemporal dynamics. However, accumulating evidence suggests high levels of intragenomic variations (IGVs) are common in many algae, leading to concerns about overinterpretation of molecular diversity in metabarcoding studies. In this project, we evaluated and compared IGVs in Tripos species by conducting the first high-throughput sequencing (HTS) of 18S rDNA V4 of Tripos single cells. High numbers of haplotypes (19-172) were identified in each of the 30 Tripos cells. Each cell contained one dominant haplotype with high relative abundance and many haplotypes with lower abundances. Thus, the presence of multiple minor haplotypes substantially overestimate the molecular diversity identified in metabarcoding analysis, which encompass not only interspecific and intraspecific diversities, but high levels of IGVs.

Diverse patterns of correspondence between protist metabarcodes and protist metagenome-assembled genomes

Two common approaches to study the composition of environmental protist communities are metabarcoding and metagenomics. Raw metabarcoding data are usually processed into Operational Taxonomic Units (OTUs) or amplicon sequence variants (ASVs) through clustering or denoising approaches, respectively. Analogous approaches are used to assemble metagenomic reads into metagenome-assembled genomes (MAGs). Understanding the correspondence between the data produced by these two approaches can help to integrate information between the datasets and to explain how metabarcoding OTUs and MAGs are related with the underlying biological entities they are hypothesised to represent. MAGs do not contain the commonly used barcoding loci, therefore sequence homology approaches cannot be used to match OTUs and MAGs. We made an attempt to match V9 metabarcoding OTUs from the 18S rRNA gene (V9 OTUs) and MAGs from the Tara Oceans expedition based on the correspondence of their relative abundances across the same set of samples. We evaluated several metrics for detecting correspondence between features in these two datasets and developed controls to filter artefacts of data structure and processing. After selecting the best-performing metrics, ranking the V9 OTU/MAG matches by their proportionality/correlation coefficients and applying a set of selection criteria, we identified candidate matches between V9 OTUs and MAGs. In some cases, V9 OTUs and MAGs could be matched with a one-to-one correspondence, implying that they likely represent the same underlying biological entity. More generally, matches we observed could be classified into 4 scenarios: one V9 OTU matches many MAGs; many V9 OTUs match many MAGs; many V9 OTUs match one MAG; one V9 OTU matches one MAG. Notably, we found some instances in which different OTU-MAG matches from the same taxonomic group were not classified in the same scenario, with all four scenarios possible even within the same taxonomic group, illustrating that factors beyond taxonomic lineage influence the relationship between OTUs and MAGs. Overall, each scenario produces a different interpretation of V9 OTUs, MAGs and how they compare in terms of the genomic and ecological diversity they represent.

Phylogenomic position of genetically diverse phagotrophic stramenopile flagellates in the sediment-associated MAST-6 lineage and a potentially halotolerant placididean

Unlike morphologically conspicuous ochrophytes, many flagellates belonging to basally branching stramenopiles are small and often overlooked. As a result, many of these lineages are known only through molecular surveys and identified as MArine STramenopiles (MAST), and remain largely uncharacterized at the cellular or genomic level. These likely phagotrophic flagellates are not only phylogenetically diverse, but also extremely abundant in some environments, making their characterization all the more important. MAST-6 is one example of a phylogenetically distinct group that has been known to be associated with sediments, but little else is known about it. Indeed, until the present study, only a single species from this group, Pseudophyllomitus vesiculosus (Pseudophyllomitidae), has been both formally described and associated with genomic information. Here, we describe four new species including two new genera of sediment-dwelling MAST-6, Vomastramonas tehuelche gen. et sp. nov., Mastreximonas tlaamin gen. et sp. nov., one undescribed Pseudophyllomitus sp., BSC2, and a new species belonging to Placididea, the potentially halotolerant Haloplacidia sinai sp. nov. We also provide two additional bikosian transcriptomes from a public culture collection, to allow for better phylogenetic reconstructions of deep-branching stramenopiles. With the SSU rRNA sequences of the new MAST-6 species, we investigate the phylogenetic diversity of the MAST-6 group and show a high relative abundance of MAST-6 related to M. tlaamin in samples across various depths and geographical locations. Using the new MAST-6 species, we also update the phylogenomic tree of stramenopiles, particularly focusing on the paraphyly of Bigyra.

Covariation between microeukaryotes and bacteria associated with Planorbidae snails

Background

Microbial communities associated with macroorganisms might affect host physiology and homeostasis. Bacteria are well studied in this context, but the diversity of microeukaryotes, as well as covariations with bacterial communities, remains almost unknown.

Methods

To study microeukaryotic communities associated with Planorbidae snails, we developed a blocking primer to reduce amplification of host DNA during metabarcoding analyses. Analyses of alpha and beta diversities were computed to describe microeukaryotes and bacteria using metabarcoding of 18S and 16S rRNA genes, respectively.

Results

Only three phyla (Amoebozoa, Opisthokonta and Alveolata) were dominant for microeukaryotes. Bacteria were more diverse with five dominant phyla (Proteobacteria, Bacteroidetes, Tenericutes, Planctomycetes and Actinobacteria). The composition of microeukaryotes and bacteria were correlated for the Biomphalaria glabrata species, but not for Planorbarius metidjensis. Network analysis highlighted clusters of covarying taxa. Among them, several links might reflect top-down control of bacterial populations by microeukaryotes, but also possible competition between microeukaryotes having opposite distributions (Lobosa and Ichthyosporea). The role of these taxa remains unknown, but we believe that the blocking primer developed herein offers new possibilities to study the hidden diversity of microeukaryotes within snail microbiota, and to shed light on their underestimated interactions with bacteria and hosts.

18S rDNA gene metabarcoding of microeukaryotes and epi-endophytes in the holobiome of seven species of large brown algae

Brown algae (Phaeophyceae) are habitat-forming species in coastal ecosystems and include kelp forests and seaweed beds that support a wide diversity of marine life. Host-associated microbial communities are an integral part of phaeophyte biology, and whereas the bacterial microbial partners have received considerable attention, the microbial eukaryotes associated with brown algae have hardly been studied. Here, we used broadly targeted "pan-eukaryotic" primers (metabarcoding) to investigate brown algal-associated eukaryotes (the eukaryome). Using this approach, we aimed to investigate the eukaryome of seven large brown algae that are important and common species in coastal ecosystems. We also aimed to assess whether these macroalgae harbor novel eukaryotic diversity and to ascribe putative functional roles to the host-associated eukaryome based on taxonomic affiliation and phylogenetic placement. We detected a significant diversity of microeukaryotic and algal lineages associated with the brown algal species investigated. The operational taxonomic units (OTUs) were taxonomically assigned to 10 of the eukaryotic major supergroups, including taxonomic groups known to be associated with seaweeds as epibionts, endobionts, parasites, and commensals. Additionally, we revealed previously unrecorded sequence types, including novel phaeophyte OTUs, particularly in the Fucus spp. samples, that may represent fucoid genomic variants, sequencing artifacts, or undescribed epi-/endophytes. Our results provide baseline data and technical insights that will be useful for more comprehensive seaweed eukaryome studies investigating the evidently lineage-rich and functionally diverse symbionts of brown algae.

Intragenomic diversity of the V9 hypervariable domain in eukaryotes has little effect on metabarcoding

Metabarcoding revolutionized our understanding of diversity and ecology of microorganisms in different habitats. However, it is also associated with several inherent biases, one of which is associated with intragenomic diversity of a molecular barcode. Here, we compare intragenomic variability of the V9 region of the 18S rRNA gene in 19 eukaryotic phyla abundant in marine plankton. The level of intragenomic variability is comparable across all the phyla, and in most genomes and transcriptomes one V9 sequence and one OTU is predominant. However, most of the variability observed at the barcode level is probably caused by sequencing errors and is mitigated by using a denoising tool, DADA2. The SWARM algorithm commonly used in metabarcoding studies is not optimal for collapsing genuine and erroneous sequences into a single OTU, leading to an overestimation of diversity in metabarcoding data. For an unknown reason, SWARM inflates diversity of eupelagonemids more than that of other eukaryotes.

Distribution Patterns of Benthic Protist Communities Depending on Depth Revealed by Environmental Sequencing-From the Sublittoral to the Deep Sea

Protists are key components of the microbial food web in marine pelagic systems because they link algal and bacterial production to higher trophic levels. However, their functioning and bathymetric distribution in benthic deep-sea ecosystems are still only poorly understood. However, biogeographical patterns of communities can be coupled to the functioning of ecosystems and are therefore important to understand ecological and evolutionary processes. In this study, we investigated the diversity and distribution of benthic protist communities from the sublittoral down to the deep seafloor (50-2000 m) around three islands of the Azores in the North Atlantic Ocean. Using amplicon sequencing of the V9 region (18S rDNA) of 21 samples, we found that protist community compositions from different depths were significantly different. Three assemblages were separated along the following depths: 50 m, 150-500 m and 1000-2000 m, which indicate that deep-sea areas surrounding islands might act as isolating barriers for benthic protist species. A limited gene flow between the communities could favor speciation processes, leading to the unique protist communities found at the different investigated islands.

Single-Cell Genomics Reveals the Divergent Mitochondrial Genomes of Retaria (Foraminifera and Radiolaria)

Mitochondria originated from an ancient bacterial endosymbiont that underwent reductive evolution by gene loss and endosymbiont gene transfer to the nuclear genome. The diversity of mitochondrial genomes published to date has revealed that gene loss and transfer processes are ongoing in many lineages. Most well-studied eukaryotic lineages are represented in mitochondrial genome databases, except for the superphylum Retaria-the lineage comprising Foraminifera and Radiolaria. Using single-cell approaches, we determined two complete mitochondrial genomes of Foraminifera and two nearly complete mitochondrial genomes of radiolarians. We report the complete coding content of an additional 14 foram species. We show that foraminiferan and radiolarian mitochondrial genomes contain a nearly fully overlapping but reduced mitochondrial gene complement compared to other sequenced rhizarians. In contrast to animals and fungi, many protists encode a diverse set of proteins on their mitochondrial genomes, including several ribosomal genes; however, some aerobic eukaryotic lineages (euglenids, myzozoans, and chlamydomonas-like algae) have reduced mitochondrial gene content and lack all ribosomal genes. Similar to these reduced outliers, we show that retarian mitochondrial genomes lack ribosomal protein and tRNA genes, contain truncated and divergent small and large rRNA genes, and contain only 14 or 15 protein-coding genes, including nad1, -3, -4, -4L, -5, and -7, cob, cox1, -2, and -3, and atp1, -6, and -9, with forams and radiolarians additionally carrying nad2 and nad6, respectively. In radiolarian mitogenomes, a noncanonical genetic code was identified in which all three stop codons encode amino acids. Collectively, these results add to our understanding of mitochondrial genome evolution and fill in one of the last major gaps in mitochondrial sequence databases. IMPORTANCE We present the reduced mitochondrial genomes of Retaria, the rhizarian lineage comprising the phyla Foraminifera and Radiolaria. By applying single-cell genomic approaches, we found that foraminiferan and radiolarian mitochondrial genomes contain an overlapping but reduced mitochondrial gene complement compared to other sequenced rhizarians. An alternative genetic code was identified in radiolarian mitogenomes in which all three stop codons encode amino acids. Collectively, these results shed light on the divergent nature of the mitochondrial genomes from an ecologically important group, warranting further questions into the biological underpinnings of gene content variability and genetic code variation between mitochondrial genomes.

Intra-genomic rRNA gene variability of Nassellaria and Spumellaria (Rhizaria, Radiolaria) assessed by Sanger, MinION and Illumina sequencing

Ribosomal RNA (rRNA) genes are known to be valuable markers for the barcoding of eukaryotic life and its phylogenetic classification at various taxonomic levels. The large-scale exploration of environmental microbial diversity through metabarcoding approaches have been focused mainly on the V4 and V9 regions of the 18S rRNA gene. The accurate interpretation of such environmental surveys is hampered by technical (e.g., PCR and sequencing errors) and biological biases (e.g., intra-genomic variability). Here we explored the intra-genomic diversity of Nassellaria and Spumellaria specimens (Radiolaria) by comparing Sanger sequencing with Illumina and Oxford Nanopore Technologies (MinION). Our analysis determined that intra-genomic variability of Nassellaria and Spumellaria is generally low, yet some Spumellaria specimens showed two different copies of the V4 with <97% similarity. From the different sequencing methods, Illumina showed the highest number of contaminations (i.e., environmental DNA, cross-contamination, tag-jumping), revealed by its high sequencing depth; and MinION showed the highest sequencing rate error (~14%). Yet the long reads produced by MinION (~2900 bp) allowed accurate phylogenetic reconstruction studies. These results highlight the requirement for a careful interpretation of Illumina based metabarcoding studies, in particular regarding low abundant amplicons, and open future perspectives towards full-length rDNA environmental metabarcoding surveys.

Diversity and ecology of Radiolaria in modern oceans

Among the many inhabitants of planktonic communities, several lineages have biomineralized intricate skeletons. These have existed for millions of years and include the Radiolaria, a group of marine protists, many of which bear delicate mineral skeletons of different natures. Radiolaria are well known for their paleontological signatures, but little is known about the ecology of modern assemblages. They are found from polar to tropical regions, in the sunlit layers of the ocean down to the deep and cold bathypelagic. They are closely involved in the biogeochemical cycles of silica, carbon and strontium sulfate, carrying important amounts of such elements to the deep ocean. However, relatively little is known on the actual extent of genetic diversity or biogeographic patterns. The rapid emergence and acceptance of molecular approaches have nevertheless led to major advances in our understanding of diversity within and evolutionary relationships between major radiolarian groups. Here, we review the state of knowledge relating to the classification, diversity and ecology of extant radiolarian orders, highlighting the substantial gaps in our understanding of the extent of their contribution to marine biodiversity and their role in marine food webs.

Coupling between Ribotypic and Phenotypic Traits of Protists across Life Cycle Stages and Temperatures

Relationships between ribotypic and phenotypic traits of protists across life cycle stages remain largely unknown. Herein, we used single cells of two soil and two marine ciliate species to examine phenotypic and ribotypic traits and their relationships across lag, log, plateau, cystic stages and temperatures. We found that Colpoda inflata and Colpoda steinii demonstrated allometric relationships between 18S ribosomal DNA (rDNA) copy number per cell (CNPC), cell volume (CV), and macronuclear volume across all life cycle stages. Integrating previously reported data of Euplotes vannus and Strombidium sulcatum indicated taxon-dependent rDNA CNPC-CV functions. Ciliate and prokaryote data analysis revealed that the rRNA CNPC followed a unified power-law function only if the rRNA-deficient resting cysts were not considered. Hence, a theoretical framework was proposed to estimate the relative quantity of resting cysts in the protistan populations with total cellular rDNA and rRNA copy numbers. Using rDNA CNPC was a better predictor of growth rate at a given temperature than rRNA CNPC and CV, suggesting replication of redundant rDNA operons as a key factor that slows cell division. Single-cell high-throughput sequencing and analysis after correcting sequencing errors revealed multiple rDNA and rRNA variants per cell. Both encystment and temperature affected the number of rDNA and rRNA variants in several cases. The divergence of rDNA and rRNA sequence in a single cell ranged from 1% to 10% depending on species. These findings have important implications for inferring cell-based biological traits (e.g., species richness, abundance and biomass, activity, and community structure) of protists using molecular approaches. IMPORTANCE Based on phenotypic traits, traditional surveys usually characterize organismal richness, abundance, biomass, and growth potential to describe diversity, organization, and function of protistan populations and communities. The rRNA gene (rDNA) and its transcripts have been widely used as molecular markers in ecological studies of protists. Nevertheless, the manner in which these molecules relate to cellular (organismal) and physiological traits remains poorly understood, which could lead to misinterpretations of protistan diversity and ecology. The current research highlights the dynamic nature of cellular rDNA and rRNA contents, which tightly couple with multiple phenotypic traits in ciliated protists. We demonstrate that quantity of resting cysts and maximum growth rate of a population can be theoretically estimated using ribotypic trait-based models. The intraindividual sequence polymorphisms of rDNA and rRNA can be influenced by encystment and temperature, which should be considered when interpreting species-level diversity and community structure of microbial eukaryotes.

Definition of a High-Resolution Molecular Marker for Tracking the Genetic Diversity of the Harmful Algal Species Eucampia zodiacus Through Comparative Analysis of Mitochondrial Genomes

The cosmopolitan phytoplankton species Eucampia zodiacus is a common harmful algal bloom (HAB) species that have been found to cause HABs in essentially all coastal regions except the Polar regions. However, molecular information for this HAB species is limited with only a few molecular markers. In this project, we constructed the mitochondrial genome (mtDNA) of E. zodiacus, which was also the first mtDNA constructed for any species in the order Hemiaulales that includes 145 reported species (including two additional HAB species Cerataulina bicornis and Cerataulina pelagica). Comparative analysis of eight E. zodiacus strains revealed that they could not be distinguished using common molecular markers, suggesting that common molecular markers do not have adequate resolution for distinguishing E. zodiacus strains. However, these E. zodiacus strains could be distinguished using whole mtDNAs, suggesting the presence of different genotypes due to evolutionary divergence. Through comparative analysis of the mtDNAs of multiple E. zodiacus strains, we identified a new molecular marker ezmt1 that could adequately distinguish different E. zodiacus strains isolated in various coastal regions in China. This molecular marker ezmt1, which was ∼400 bp in size, could be applied to identify causative genotypes during E. zodiacus HABs through tracking the dynamic changes of genetic diversity of E. zodiacus in HABs.

Core microbial communities of lacustrine microbialites sampled along an alkalinity gradient

Microbialites are usually carbonate-rich sedimentary rocks formed by the interplay of phylogenetically and metabolically complex microbial communities with their physicochemical environment. Yet, the biotic and abiotic determinants of microbialite formation remain poorly constrained. Here, we analysed the structure of prokaryotic and eukaryotic communities associated with microbialites occurring in several crater lakes of the Trans-Mexican volcanic belt along an alkalinity gradient. Microbialite size and community structure correlated with lake physicochemical parameters, notably alkalinity. Although microbial community composition varied across lake microbialites, major taxa-associated functions appeared quite stable with both, oxygenic and anoxygenic photosynthesis and, to less extent, sulphate reduction, as major putative carbonatogenic processes. Despite interlake microbialite community differences, we identified a microbial core of 247 operational taxonomic units conserved across lake microbialites, suggesting a prominent ecological role in microbialite formation. This core mostly encompassed Cyanobacteria and their typical associated taxa (Bacteroidetes, Planctomycetes) and diverse anoxygenic photosynthetic bacteria, notably Chloroflexi, Alphaproteobacteria (Rhodobacteriales, Rhodospirilalles), Gammaproteobacteria (Chromatiaceae) and minor proportions of Chlorobi. The conserved core represented up to 40% (relative abundance) of the total community in lakes Alchichica and Atexcac, displaying the highest alkalinities and the most conspicuous microbialites. Core microbialite communities associated with carbonatogenesis might be relevant for inorganic carbon sequestration purposes.

Planktonic foraminifera genomic variations reflect paleoceanographic changes in the Arctic: evidence from sedimentary ancient DNA

Deciphering the evolution of marine plankton is typically based on the study of microfossil groups. Cryptic speciation is common in these groups, and large intragenomic variations occur in ribosomal RNA genes of many morphospecies. In this study, we correlated the distribution of ribosomal amplicon sequence variants (ASVs) with paleoceanographic changes by analyzing the high-throughput sequence data assigned to Neogloboquadrina pachyderma in a 140,000-year-old sediment core from the Arctic Ocean. The sedimentary ancient DNA demonstrated the occurrence of various N. pachyderma ASVs whose occurrence and dominance varied through time. Most remarkable was the striking appearance of ASV18, which was nearly absent in older sediments but became dominant during the last glacial maximum and continues to persist today. Although the molecular ecology of planktonic foraminifera is still poorly known, the analysis of their intragenomic variations through time has the potential to provide new insight into the evolution of marine biodiversity and may lead to the development of new and important paleoceanographic proxies.

Unveiling protist diversity associated with the Pacific oyster Crassostrea gigas using blocking and excluding primers

BACKGROUND

Microbiome of macroorganisms might directly or indirectly influence host development and homeostasis. Many studies focused on the diversity and distribution of prokaryotes within these assemblages, but the eukaryotic microbial compartment remains underexplored so far.

RESULTS

To tackle this issue, we compared blocking and excluding primers to analyze microeukaryotic communities associated with Crassostrea gigas oysters. High-throughput sequencing of 18S rRNA genes variable loops revealed that excluding primers performed better by not amplifying oyster DNA, whereas the blocking primer did not totally prevent host contaminations. However, blocking and excluding primers showed similar pattern of alpha and beta diversities when protist communities were sequenced using metabarcoding. Alveolata, Stramenopiles and Archaeplastida were the main protist phyla associated with oysters. In particular, Codonellopsis, Cyclotella, Gymnodinium, Polarella, Trichodina, and Woloszynskia were the dominant genera. The potential pathogen Alexandrium was also found in high abundances within some samples.

CONCLUSIONS

Our study revealed the main protist taxa within oysters as well as the occurrence of potential oyster pathogens. These new primer sets are promising tools to better understand oyster homeostasis and disease development, such as the Pacific Oyster Mortality Syndrome (POMS) targeting juveniles.

Perspectives from Ten Years of Protist Studies by High-Throughput Metabarcoding

During the last decade, high-throughput metabarcoding became routine for analyzing protistan diversity and distributions in nature. Amid a multitude of exciting findings, scientists have also identified and addressed technical and biological limitations, although problems still exist for inference of meaningful taxonomic and ecological knowledge based on short DNA sequences. Given the extensive use of this approach, it is critical to settle our understanding on its strengths and weaknesses and to synthesize up-to-date methodological and conceptual trends. This article summarizes key scientific and technical findings, and identifies current and future directions in protist research that uses metabarcoding.

Comparative analyses of the V4 and V9 regions of 18S rDNA for the extant eukaryotic community using the Illumina platform

Illumina sequencing is a representative tool for understanding the massive diversity of microbial eukaryotes in natural ecosystems. Here, we investigated the eukaryotic community in a pond (salinity of 2–4) on Dokdo (island) in the East Sea, Korea, using Illumina sequencing with primer sets for the V4 and V9 regions of 18S rDNA from 2016 to 2018 for the first time. Totally, 1,413 operational taxonomic units (OTUs) and 915 OTUs were detected using the V9 and V4 primer sets, respectively. Taxonomic analyses of these OTUs revealed that although the V4 primer set failed to describe the extant diversity for some major sub-division groups, the V9 primer set represented their diversity. Moreover, the rare taxa with <1% of total reads were exclusively detected using V9 primer set. Hence, the diversity of the eukaryotic community can vary depending on the choice of primers. The Illumina sequencing data of the V9 region of 18S rDNA may be advantageous for estimating the richness of the eukaryotic community including a rare biosphere, whereas the simultaneous application of two biomarkers may be suitable for understanding the molecular phylogenetic relationships. We strongly recommend both biomarkers be used to assess the diversity and phylogenetic relationship within the eukaryotic community in natural samples.

Species detection and delineation in the marine planktonic diatoms Chaetoceros and Bacteriastrum through metabarcoding: making biological sense of haplotype diversity

High-throughput sequencing (HTS) metabarcoding is commonly applied to assess phytoplankton diversity. Usually, haplotypes are grouped into operational taxonomic units (OTUs) through clustering, whereby the resulting number of OTUs depends on chosen similarity thresholds. We applied, instead, a phylogenetic approach to infer taxa among 18S rDNA V4-metabarcode haplotypes gathered from 48 time-series samples using the marine planktonic diatoms Chaetoceros and Bacteriastrum as test case. The 73 recovered taxa comprised both solitary haplotypes and polytomies, the latter composed each of a highly abundant, dominant haplotype and one to several minor, peripheral haplotypes. The solitary and dominant haplotypes usually matched reference sequences, enabling species assignation of taxa. We hypothesise that the super-abundance of reads in dominant haplotypes results from the homogenization effect of concerted evolution. Reads of populous peripheral haplotypes and dominant haplotypes show comparable distribution patterns over the sample dates, suggesting that they are part of the same population. Many taxa revealed marked seasonality, with closely related ones generally showing distinct periodicity, whereas others occur year-round. Phylogenies inferred from metabarcode haplotypes enable delineation of biologically meaningful taxa, whereas OTUs resulting from clustering algorithms often deviate markedly from such taxa.

Widespread prevalence but contrasting patterns of intragenomic rRNA polymorphisms in nematodes: Implications for phylogeny, species delimitation and life history inference

Ribosomal RNA genes have long been a favoured locus in phylogenetic and metabarcoding studies. Within a genome, rRNA loci are organized as tandem repeated arrays and the copies are homogenized through the process of concerted evolution. However, some level of rRNA variation (intragenomic polymorphism) is known to persist and be maintained in the genomes of many species. In nematode worms, the extent of rRNA polymorphism (RP) across species and the evolutionary and life history factors that contribute to the maintenance of intragenomic RP is largely unknown. Here, we present an extensive analysis across 30 terrestrial nematode species representing a range of free-living and parasitic taxa isolated worldwide. Our results indicate that RP is common and widespread, ribosome function appears to be maintained despite mutational changes, and intragenomic variants are stable in the genome and neutrally evolving. However, levels of variation were varied widely across rRNA locus and species, with some taxa observed to lack RP entirely. Higher levels of RP were significantly correlated with shorter generation time and high reproductive rates, and population-level factors may play a role in the geographic and phylogenetic structuring of rRNA variants observed in genera such as Rotylenchulus and Pratylenchus. Although RP did not dramatically impact the clustering and recovery of taxa in mock metabarcoding analyses, the present study has significant implications for global biodiversity estimates of nematode species derived from environmental rRNA amplicon studies, as well as our understanding of the evolutionary and ecological factors shaping genetic diversity across the nematode Tree of Life.

Prospects and challenges of implementing DNA metabarcoding for high-throughput insect surveillance

Trap-based surveillance strategies are widely used for monitoring of invasive insect species, aiming to detect newly arrived exotic taxa as well as track the population levels of established or endemic pests. Where these surveillance traps have low specificity and capture non-target endemic species in excess of the target pests, the need for extensive specimen sorting and identification creates a major diagnostic bottleneck. While the recent development of standardized molecular diagnostics has partly alleviated this requirement, the single specimen per reaction nature of these methods does not readily scale to the sheer number of insects trapped in surveillance programmes. Consequently, target lists are often restricted to a few high-priority pests, allowing unanticipated species to avoid detection and potentially establish populations. DNA metabarcoding has recently emerged as a method for conducting simultaneous, multi-species identification of complex mixed communities and may lend itself ideally to rapid diagnostics of bulk insect trap samples. Moreover, the high-throughput nature of recent sequencing platforms could enable the multiplexing of hundreds of diverse trap samples on a single flow cell, thereby providing the means to dramatically scale up insect surveillance in terms of both the quantity of traps that can be processed concurrently and number of pest species that can be targeted. In this review of the metabarcoding literature, we explore how DNA metabarcoding could be tailored to the detection of invasive insects in a surveillance context and highlight the unique technical and regulatory challenges that must be considered when implementing high-throughput sequencing technologies into sensitive diagnostic applications.

Coupled effects of land use pattern and hydrological regime on composition and diversity of riverine eukaryotic community in a coastal watershed of Southeast China

The coupled effects of land use pattern and hydrological regime on composition and diversity of riverine eukaryotic community are needed for understanding riverine ecosystem health and algal blooming mechanism. In-situ monitoring and 18S rRNA gene sequencing were used to investigate spatiotemporal variations of eukaryotic community in three types of watershed with different dominant land use type (i.e. urban, forest, and natural) during three seasons (i.e. dry, transition, and wet seasons) in a coastal watershed of Southeast China. Results showed that agricultural and urban watersheds had significantly higher diversity in dry and transition seasons, and higher richness in transition and wet seasons than those in forest watershed. The non-metric multidimensional scaling analysis further verified great spatiotemporal variations of eukaryotic community. Stramenopiles, Alveolata, Animalia, and Eukaryota, dominated the sequences reads for all sampling sites in three seasons. Agricultural watershed had the highest relative abundant of Animalia, whereas Eukaryota was the most abundant in urban watershed and forest watershed had the highest relative abundance of Stramenopiles and Alveolata. The RDA ordination showed that Builtup and streamflow were two most important factors for moderate taxa and abundant taxa, respectively. Variation partitioning revealed that land use pattern and hydrological regime together explained 54.4%, 61% and 67.2% variances of the composition of eukaryotic community. Among three sampling seasons, the relative contribution of land use pattern was higher than that of hydrological regime. The results of this investigation demonstrated how land use pattern and hydrological regime affected the composition and diversity of riverine eukaryotic community. The findings can provide a useful insight into the riverine eukaryotic communities and their underlying ecological mechanisms in coastal China watersheds.

High contribution of Rhizaria (Radiolaria) to vertical export in the California Current Ecosystem revealed by DNA metabarcoding

Passive sinking of particulate organic matter (POM) is the main mechanism through which the biological pump transports surface primary production to the ocean interior. However, the contribution and variability of different biological sources to vertical export is not fully understood. Here, we use DNA metabarcoding of the 18S rRNA gene and particle interceptor traps (PITs) to characterize the taxonomic composition of particles sinking out of the photic layer in the California Current Ecosystem (CCE), a productive system with high export potential. The PITs included formalin-fixed and 'live' traps to investigate eukaryotic communities involved in the export and remineralization of sinking particles. Sequences affiliated with Radiolaria dominated the eukaryotic assemblage in fixed traps (90%), with Dinophyta and Metazoa making minor contributions. The prominence of Radiolaria decreased drastically in live traps, possibly due to selective consumption by copepods, heterotrophic nanoflagellates, and phaeodarians that were heavily enriched in these traps. These patterns were consistent across the water masses surveyed extending from the coast to offshore, despite major differences in productivity and trophic structure of the epipelagic plankton community. Our findings identify Radiolaria as major actors in export fluxes in the CCE.

A census-based estimate of Earth's bacterial and archaeal diversity

The global diversity of Bacteria and Archaea, the most ancient and most widespread forms of life on Earth, is a subject of intense controversy. This controversy stems largely from the fact that existing estimates are entirely based on theoretical models or extrapolations from small and biased data sets. Here, in an attempt to census the bulk of Earth's bacterial and archaeal ("prokaryotic") clades and to estimate their overall global richness, we analyzed over 1.7 billion 16S ribosomal RNA amplicon sequences in the V4 hypervariable region obtained from 492 studies worldwide, covering a multitude of environments and using multiple alternative primers. From this data set, we recovered 739,880 prokaryotic operational taxonomic units (OTUs, 16S-V4 gene clusters at 97% similarity), a commonly used measure of microbial richness. Using several statistical approaches, we estimate that there exist globally about 0.8–1.6 million prokaryotic OTUs, of which we recovered somewhere between 47%–96%, representing >99.98% of prokaryotic cells. Consistent with this conclusion, our data set independently "recaptured" 91%–93% of 16S sequences from multiple previous global surveys, including PCR-independent metagenomic surveys. The distribution of relative OTU abundances is consistent with a log-normal model commonly observed in larger organisms; the total number of OTUs predicted by this model is also consistent with our global richness estimates. By combining our estimates with the ratio of full-length versus partial-length (V4) sequence diversity in the SILVA sequence database, we further estimate that there exist about 2.2–4.3 million full-length OTUs worldwide. When restricting our analysis to the Americas, while controlling for the number of studies, we obtain similar richness estimates as for the global data set, suggesting that most OTUs are globally distributed. Qualitatively similar results are also obtained for other 16S similarity thresholds (90%, 95%, and 99%). Our estimates constrain the extent of a poorly quantified rare microbial biosphere and refute recent predictions that there exist trillions of prokaryotic OTUs.

A massive survey of Earth's Bacteria and Archaea reveals that their diversity is orders of magnitude lower than previously thought. The study also indicates that extinctions played an important role in prokaryotic evolution.

The global diversity of Bacteria and Archaea ("prokaryotes"), the most ancient and most widespread forms of life on Earth, is subject to high uncertainty. Here, to estimate the global diversity of prokaryotes, we analyzed a large number of 16S ribosomal RNA gene sequences, found in all prokaryotes and commonly used to catalogue prokaryotic diversity. Sequences were obtained from a multitude of environments across thousands of geographic locations worldwide. From this data set, we recovered 739,880 prokaryotic operational taxonomic units (OTUs), i.e., 16S gene clusters sharing 97% similarity, roughly corresponding to prokaryotic species. Using several statistical approaches and through comparison with existing databases and previous independent surveys, we estimate that there exist globally between 0.8 and 1.6 million prokaryotic OTUs. When restricting our analysis to the Americas, while controlling for the number of studies, we obtain similar estimates as for the global data set, suggesting that most OTUs are not restricted to a single continent but are instead globally distributed. Our estimates constrain the extent of a commonly hypothesized but poorly quantified rare prokaryotic biosphere and refute recent predictions that there exists trillions of prokaryotic OTUs. Our findings also indicate that, contrary to common speculation, extinctions may strongly influence global prokaryotic diversity.

Efecto del compost en la diversidad de bacterias rizosféricas del cultivo de morera (Morus alba)

La producción orgánica de seda incluye la aplicación de compost como practica de cultivo en morera (Morus alba), sin embargo, el efecto de la fertilización orgánica en las poblaciones de bacterias rizosféricas no siempre es positivo. Para evaluar el efecto del compost en la diversidad de bacterias rizosféricas en cultivos de morera (Morus alba), se aplicaron 0, 0.25, 0.5 y 1 kg.m-2 de compost a parcelas con morera dispuestas en un diseño en bloques completos al azar. De cada parcela se extrajo ADN del suelo rizosférico a los 0, 5, 10 y 90 días de aplicado el compost y se amplificó la región V4 del gen ADNr 16S para su secuenciación y asignación taxonómica de los OTUS. Los índices de diversidad alpha mostraron la dominancia de algunos grupos taxonómicos, como los phyla Proteobacteria y Acidobacteria y los géneros Pseudomonas, Opitutus, Luteolibacter y Nitrospir. La diversidad beta indicó similitud entre las muestras influenciadas por la aplicación compost y el incremento de la diversidad en las parcelas muestreadas al final del experimento (90 días). Los grupos taxonómicos dominantes se caracterizan por su función en el ciclo del nitrógeno. Así, se concluyó que la aplicación de 1 kg.m-2 llevó al aumento de la humedad del suelo, el pH y la disponibilidad de nutrientes, lo que incremento la diversidad de bacterias rizosféricas con cambios positivos en composición, riqueza y abundancia en los niveles de orden, familia y género.

Microeukaryote community in a partial nitritation reactor prior to anammox and an insight into the potential of ciliates as performance bioindicators

An in-depth, long-term, multidisciplinary study was conducted in order to study the microeukaryote community in a partial nitritation (PN) reactor prior to anammox. The PN reactor operated with moving bed biofilm reactor (MBBR) technology, using plastic supports (carriers) for biofilm development. The microeukaryote community from the biofilm (BF) and the surrounding media (mixed liquor or ML) were analysed separately. Despite the physicochemical conditions under which the PN-MBBR operated (an average of 305.9±117mg TAN l^-1 and 328.4±131.9mg N-NO₂^- l^-1), up to 24 microeukaryotic taxa were observed by microscope. Microeukaryote species showed an uneven distribution in the PN-MBBR, thus suggesting the existence of two habitats: the BF, preferred by species with specific structures for adhering to a substrate, such as the stalked Peritrichia, and the ML, preferred by free-swimming or non-substrate dependent species. The results indicated that most ciliate population dynamics mainly responded to the nitrous acid and free ammonia concentrations and, to a lesser extent, to sCOD values. In the BF, variations in the population of Epistylis camprubii and Opercularia coarctata suggest the existence of competition between these species due to niche overlap. A V4 18S rDNA molecular survey (Illumina) was carried out for some samples with the aim of obtaining maximum coverage of the main eukaryote species that were microscopically detected throughout the study. The diversity and abundance data provided by both detection methods were compared. The study helped identify broader tolerance ranges of the microeukaryote taxa to the physicochemical parameters analysed.

Surface ocean metabarcoding confirms limited diversity in planktonic foraminifera but reveals unknown hyper-abundant lineages

Since the advent of DNA metabarcoding surveys, the planktonic realm is considered a treasure trove of diversity, inhabited by a small number of abundant taxa, and a hugely diverse and taxonomically uncharacterized consortium of rare species. Here we assess if the apparent underestimation of plankton diversity applies universally. We target planktonic foraminifera, a group of protists whose known morphological diversity is limited, taxonomically resolved and linked to ribosomal DNA barcodes. We generated a pyrosequencing dataset of ~100,000 partial 18S rRNA foraminiferal sequences from 32 size fractioned photic-zone plankton samples collected at 8 stations in the Indian and Atlantic Oceans during the Tara Oceans expedition (2009–2012). We identified 69 genetic types belonging to 41 morphotaxa in our metabarcoding dataset. The diversity saturated at local and regional scale as well as in the three size fractions and the two depths sampled indicating that the diversity of foraminifera is modest and finite. The large majority of the newly discovered lineages occur in the small size fraction, neglected by classical taxonomy. These unknown lineages dominate the bulk [>0.8 µm] size fraction, implying that a considerable part of the planktonic foraminifera community biomass has its origin in unknown lineages.

The diversity and biogeography of abundant and rare intertidal marine microeukaryotes explained by environment and dispersal limitation

Benthic microeukaryotes are key ecosystem drivers in marine sandy beaches, an important and dynamic environment; however, little is known about their diversity and biogeography on a large spatial scale. Here, we investigated the community composition and geographical distributions of benthic microeukaryotes using high-throughput sequencing of the 18S rRNA gene and quantified the contributions of environmental factors and spatial separation on the distribution patterns of both rare and abundant taxa. We collected 36 intertidal samples at 12 sandy beaches from four regions that spanned distances from 0.001 to 12,000 km. We found 12,890 operational taxonomic units (OTUs; 97% sequence identity level) including members of all eukaryotic super-groups and several phyla of uncertain position. Arthropoda and Diatomeae dominated the sequence reads in abundance, but Ciliophora and Discoba were the most diverse groups across all samples. About one-third of the OTUs could not be definitively classified at a similarity level of 80%, supporting the view that a large number of rare and minute marine species may have escaped previous characterization. We found generally similar geographical patterns for abundant and rare microeukaryotic sub-communities, and both showed a significant distance-decay similarity trend. Variation partitioning showed that both rare and abundant sub-communities exhibited a slightly stronger response to environmental factors than spatial (distance) factors. However, the abundant sub-community was strongly correlated with variations in spatial, environmental and sediment grain size factors (66% of variance explained), but the rare assemblage was not (16%). This suggests that different or more complex mechanisms generate and maintain diversity in the rare biosphere in this habitat.

Censusing marine eukaryotic diversity in the twenty-first century

The ocean constitutes one of the vastest and richest biomes on our planet. Most recent estimations, all based on indirect approaches, suggest that there are millions of marine eukaryotic species. Moreover, a large majority of these are small (less than 1 mm), cryptic and still unknown to science. However, this knowledge gap, caused by the lack of diagnostic morphological features in small organisms and the limited sampling of the global ocean, is currently being filled, thanks to new DNA-based approaches. The molecular technique of PCR amplification of homologous gene regions combined with high-throughput sequencing, routinely used to census unculturable prokaryotes, is now also being used to characterize whole communities of marine eukaryotes. Here, we review how this methodological advancement has helped to better quantify the magnitude and patterns of marine eukaryotic diversity, with an emphasis on taxonomic groups previously largely overlooked. We then discuss obstacles remaining to achieve a global understanding of marine eukaryotic diversity. In particular, we argue that 18S variable regions do not provide sufficient taxonomic resolution to census marine life, and suggest combining broad eukaryotic surveys targeting the 18S rRNA region with more taxon-focused analyses of hypervariable regions to improve our understanding of the diversity of species, the functional units of marine ecosystems.

This article is part of the themed issue ‘From DNA barcodes to biomes’.

Evaluation of the reproducibility of amplicon sequencing with Illumina MiSeq platform

Illumina’s MiSeq has become the dominant platform for gene amplicon sequencing in microbial ecology studies; however, various technical concerns, such as reproducibility, still exist. To assess reproducibility, 16S rRNA gene amplicons from 18 soil samples of a reciprocal transplantation experiment were sequenced on an Illumina MiSeq. The V4 region of 16S rRNA gene from each sample was sequenced in triplicate with each replicate having a unique barcode. The average OTU overlap, without considering sequence abundance, at a rarefaction level of 10,323 sequences was 33.4±2.1% and 20.2±1.7% between two and among three technical replicates, respectively. When OTU sequence abundance was considered, the average sequence abundance weighted OTU overlap was 85.6±1.6% and 81.2±2.1% for two and three replicates, respectively. Removing singletons significantly increased the overlap for both (~1–3%, p<0.001). Increasing the sequencing depth to 160,000 reads by deep sequencing increased OTU overlap both when sequence abundance was considered (95%) and when not (44%). However, if singletons were not removed the overlap between two technical replicates (not considering sequence abundance) plateaus at 39% with 30,000 sequences. Diversity measures were not affected by the low overlap as α-diversities were similar among technical replicates while β-diversities (Bray-Curtis) were much smaller among technical replicates than among treatment replicates (e.g., 0.269 vs. 0.374). Higher diversity coverage, but lower OTU overlap, was observed when replicates were sequenced in separate runs. Detrended correspondence analysis indicated that while there was considerable variation among technical replicates, the reproducibility was sufficient for detecting treatment effects for the samples examined. These results suggest that although there is variation among technical replicates, amplicon sequencing on MiSeq is useful for analyzing microbial community structure if used appropriately and with caution. For example, including technical replicates, removing spurious sequences and unrepresentative OTUs, using a clustering method with a high stringency for OTU generation, estimating treatment effects at higher taxonomic levels, and adapting the unique molecular identifier (UMI) and other newly developed methods to lower PCR and sequencing error and to identify true low abundance rare species all can increase reproducibility.

Biogeography and diversity of Collodaria (Radiolaria) in the global ocean

Collodaria are heterotrophic marine protists that exist either as large colonies composed of hundreds of cells or as large solitary cells. All described species so far harbour intracellular microalgae as photosymbionts. Although recent environmental diversity surveys based on molecular methods demonstrated their consistently high contribution to planktonic communities and their worldwide occurrence, our understanding of their diversity and biogeography is still very limited. Here we estimated the 18S ribosomal DNA (rDNA) gene copies per collodarian cell for solitary (5770±1960 small subunit (SSU) rDNA copies) and colonial specimens (37 474±17 799 SSU rDNA copies, for each individual cell within a colony) using single-specimen quantitative PCR. We then investigated the environmental diversity of Collodaria within the photic zone through the metabarcoding survey from the Tara Oceans expedition and found that the two collodarian families Collosphaeridae and Sphaerozoidae contributed the most to the collodarian diversity and encompassed mostly cosmopolitan taxa. Although the biogeographical patterns were homogeneous within each biogeochemical biome considered, we observed that coastal biomes were consistently less diverse than oceanic biomes and were dominated by the Sphaerozoidae while the Collosphaeridae were dominant in the open oceans. The significant relationships with six environmental variables suggest that collodarian diversity is influenced by the trophic status of oceanic provinces and increased towards more oligotrophic regions.

Environmental Sequencing Provides Reasonable Estimates of the Relative Abundance of Specific Picoeukaryotes

High-throughput sequencing (HTS) is revolutionizing environmental surveys of microbial diversity in the three domains of life by providing detailed information on which taxa are present in microbial assemblages. However, it is still unclear how the relative abundance of specific taxa gathered by HTS correlates with cell abundances. Here, we quantified the relative cell abundance of 6 picoeukaryotic taxa in 13 planktonic samples from 6 European coastal sites using epifluorescence microscopy on tyramide signal amplification-fluorescence in situ hybridization preparations. These relative abundance values were then compared with HTS data obtained in three separate molecular surveys: 454 sequencing of the V4 region of the 18S ribosomal DNA (rDNA) using DNA and RNA extracts (DNA-V4 and cDNA-V4) and Illumina sequencing of the V9 region (cDNA-V9). The microscopic and molecular signals were generally correlated, indicating that a relative increase in specific 18S rDNA was the result of a large proportion of cells in the given taxa. Despite these positive correlations, the slopes often deviated from 1, precluding a direct translation of sequences to cells. Our data highlighted clear differences depending on the nucleic acid template or the 18S rDNA region targeted. Thus, the molecular signal obtained using cDNA templates was always closer to relative cell abundances, while the V4 and V9 regions gave better results depending on the taxa. Our data support the quantitative use of HTS data but warn about considering it as a direct proxy of cell abundances.

IMPORTANCE

Direct studies on marine picoeukaryotes by epifluorescence microscopy are problematic due to the lack of morphological features and due to the limited number and poor resolution of specific phylogenetic probes used in fluorescence in situ hybridization (FISH) routines. As a consequence, there is an increasing use of molecular methods, including high-throughput sequencing (HTS), to study marine microbial diversity. HTS can provide a detailed picture of the taxa present in a community and can reveal diversity not evident using other methods, but it is still unclear what the meaning of the sequence abundance in a given taxon is. Our aim is to investigate the correspondence between the relative HTS signal and relative cell abundances in selected picoeukaryotic taxa. Environmental sequencing provides reasonable estimates of the relative abundance of specific taxa. Better results are obtained when using RNA extracts as the templates, while the region of 18S ribosomal DNA had different influences depending on the taxa assayed.

Comparison of three clustering approaches for detecting novel environmental microbial diversity

Discovery of novel diversity in high-throughput sequencing studies is an important aspect in environmental microbial ecology. To evaluate the effects that amplicon clustering methods have on the discovery of novel diversity, we clustered an environmental marine high-throughput sequencing dataset of protist amplicons together with reference sequences from the taxonomically curated Protist Ribosomal Reference (PR²) database using three de novo approaches: sequence similarity networks, USEARCH, and Swarm. The potentially novel diversity uncovered by each clustering approach differed drastically in the number of operational taxonomic units (OTUs) and in the number of environmental amplicons in these novel diversity OTUs. Global pairwise alignment comparisons revealed that numerous amplicons classified as potentially novel by USEARCH and Swarm were more than 97% similar to references of PR². Using shortest path analyses on sequence similarity network OTUs and Swarm OTUs we found additional novel diversity within OTUs that would have gone unnoticed without further exploiting their underlying network topologies. These results demonstrate that graph theory provides powerful tools for microbial ecology and the analysis of environmental high-throughput sequencing datasets. Furthermore, sequence similarity networks were most accurate in delineating novel diversity from previously discovered diversity.

Marine protist diversity in European coastal waters and sediments as revealed by high-throughput sequencing

Although protists are critical components of marine ecosystems, they are still poorly characterized. Here we analysed the taxonomic diversity of planktonic and benthic protist communities collected in six distant European coastal sites. Environmental deoxyribonucleic acid (DNA) and ribonucleic acid (RNA) from three size fractions (pico-, nano- and micro/mesoplankton), as well as from dissolved DNA and surface sediments were used as templates for tag pyrosequencing of the V4 region of the 18S ribosomal DNA. Beta-diversity analyses split the protist community structure into three main clusters: picoplankton-nanoplankton-dissolved DNA, micro/mesoplankton and sediments. Within each cluster, protist communities from the same site and time clustered together, while communities from the same site but different seasons were unrelated. Both DNA and RNA-based surveys provided similar relative abundances for most class-level taxonomic groups. Yet, particular groups were overrepresented in one of the two templates, such as marine alveolates (MALV)-I and MALV-II that were much more abundant in DNA surveys. Overall, the groups displaying the highest relative contribution were Dinophyceae, Diatomea, Ciliophora and Acantharia. Also, well represented were Mamiellophyceae, Cryptomonadales, marine alveolates and marine stramenopiles in the picoplankton, and Monadofilosa and basal Fungi in sediments. Our extensive and systematic sequencing of geographically separated sites provides the most comprehensive molecular description of coastal marine protist diversity to date.

Divergence thresholds and divergent biodiversity estimates: can metabarcoding reliably describe zooplankton communities?

DNA metabarcoding is a promising method for describing communities and estimating biodiversity. This approach uses high-throughput sequencing of targeted markers to identify species in a complex sample. By convention, sequences are clustered at a predefined sequence divergence threshold (often 3%) into operational taxonomic units (OTUs) that serve as a proxy for species. However, variable levels of interspecific marker variation across taxonomic groups make clustering sequences from a phylogenetically diverse dataset into OTUs at a uniform threshold problematic. In this study, we use mock zooplankton communities to evaluate the accuracy of species richness estimates when following conventional protocols to cluster hypervariable sequences of the V4 region of the small subunit ribosomal RNA gene (18S) into OTUs. By including individually tagged single specimens and "populations" of various species in our communities, we examine the impact of intra- and interspecific diversity on OTU clustering. Communities consisting of single individuals per species generated a correspondence of 59-84% between OTU number and species richness at a 3% divergence threshold. However, when multiple individuals per species were included, the correspondence between OTU number and species richness dropped to 31-63%. Our results suggest that intraspecific variation in this marker can often exceed 3%, such that a single species does not always correspond to one OTU. We advocate the need to apply group-specific divergence thresholds when analyzing complex and taxonomically diverse communities, but also encourage the development of additional filtering steps that allow identification of artifactual rRNA gene sequences or pseudogenes that may generate spurious OTUs.

Estimating Protistan Diversity Using High-Throughput Sequencing

Sequencing hypervariable regions from the 18S rRNA gene is commonly employed to characterize protistan biodiversity, yet there are concerns that short reads do not provide the same taxonomic resolution as full-length sequences. A total of 7,432 full-length sequences were used to perform an in silico analysis of how sequences of various lengths and target regions impact downstream ecological interpretations. Sequences that were longer than 400 nucleotides and included the V4 hypervariable region generated results similar to those derived from full-length 18S rRNA gene sequences. Present high-throughput sequencing capabilities are approaching protistan diversity estimation comparable to whole gene sequences.

High-throughput metagenomic technologies for complex microbial community analysis: open and closed formats

Understanding the structure, functions, activities and dynamics of microbial communities in natural environments is one of the grand challenges of 21st century science. To address this challenge, over the past decade, numerous technologies have been developed for interrogating microbial communities, of which some are amenable to exploratory work (e.g., high-throughput sequencing and phenotypic screening) and others depend on reference genes or genomes (e.g., phylogenetic and functional gene arrays). Here, we provide a critical review and synthesis of the most commonly applied “open-format” and “closed-format” detection technologies. We discuss their characteristics, advantages, and disadvantages within the context of environmental applications and focus on analysis of complex microbial systems, such as those in soils, in which diversity is high and reference genomes are few. In addition, we discuss crucial issues and considerations associated with applying complementary high-throughput molecular technologies to address important ecological questions.

Early detection of eukaryotic communities from marine biofilm using high-throughput sequencing: an assessment of different sampling devices

Marine biofilms are precursors for colonization by larger fouling organisms, including non-indigenous species (NIS). In this study, high-throughput sequencing (HTS) of 18S rRNA metabarcodes was used to investigate four sampling methods (modified syringe, sterilized sponge, underwater tape and sterilized swab) for characterizing eukaryotic communities in marine biofilms. Perspex™ plates were sampled in and out of water. DNA collected with tape did not amplify. Otherwise, there were no statistical differences in communities among the remaining three sampling devices or between the two environments. Sterilized sponges are recommended for ease of use underwater. In-depth HTS analysis identified diverse eukaryotic communities, dominated by Metazoa and Chromoalveolata. Among the latter, diatoms (Bacillariophyceae) were particularly abundant (33% of reads assigned to Chromalveolata). The NIS Ciona savignyi was detected in all samples. The application of HTS in marine biofilm surveillance could facilitate early detection of NIS, improving the probability of successful eradication.