Using digital PCR to predict ciliate abundance from ribosomal RNA gene copy numbers

Ciliates play a key role in most ecosystems. Their abundance in natural samples is crucial for answering many ecological questions. Traditional methods of quantifying individual species, which rely on microscopy, are often labour-intensive, time-consuming and can be highly biassed. As a result, we investigated the potential of digital polymerase chain reaction (dPCR) for quantifying ciliates. A significant challenge in this process is the high variation in the copy number of the taxonomic marker gene (ribosomal RNA [rRNA]). We first quantified the rRNA gene copy numbers (GCN) of the model ciliate, Paramecium tetraurelia, during different stages of the cell cycle and growth phases. The per-cell rRNA GCN varied between approximately 11,000 and 130,000, averaging around 50,000 copies per cell. Despite these variations in per-cell rRNA GCN, we found a highly significant correlation between GCN and cell numbers. This is likely due to the coexistence of different cellular stages in an uncontrolled (environmental) ciliate population. Thanks to the high sensitivity of dPCR, we were able to detect the target gene in a sample that contained only a single cell. The dPCR approach presented here is a valuable addition to the molecular toolbox in protistan ecology. It may guide future studies in quantifying and monitoring the abundance of targeted (even rare) ciliates in natural samples.

Simultaneous analysis of seven 16S rRNA hypervariable gene regions increases efficiency in marine bacterial diversity detection

Environmental DNA sequencing is the gold standard to reveal microbial community structures. In most applications, a one-fragment PCR approach is applied to amplify a taxonomic marker gene, usually a hypervariable region of the 16S rRNA gene. We used a new reverse complement (RC)-PCR-based assay that amplifies seven out of the nine hypervariable regions of the 16S rRNA gene, to interrogate bacterial communities in sediment samples collected from different coastal marine sites with an impact gradient. In parallel, we employed a traditional one-fragment analysis of the hypervariable V3-V4 region to investigate whether the RC-PCR reveals more of the 'unseen' diversity obtained by the one-fragment approach. As a benchmark for the full deck of diversity, we subjected the samples to PCR-free metagenomic sequencing. None of the two PCR-based approaches recorded the full taxonomic repertoire obtained from the metagenomics datasets. However, the RC-PCR approach detected 2.8 times more bacterial genera compared to the near-saturation sequenced V3-V4 samples. RC-PCR is an ideal compromise between the standard one-fragment approach and metagenomics sequencing and may guide future environmental sequencing studies, in which bacterial diversity is a central subject.

Comparing quantile regression spline analyses and supervised machine learning for environmental quality assessment at coastal marine aquaculture installations

Organic enrichment associated with marine finfish aquaculture is a local stressor of marine coastal ecosystems. To maintain ecosystem services, the implementation of biomonitoring programs focusing on benthic diversity is required. Traditionally, impact-indices are determined by extracting and identifying benthic macroinvertebrates from samples. However, this is a time-consuming and expensive method with low upscaling potential. A more rapid, inexpensive, and robust method to infer the environmental quality of marine environments is eDNA metabarcoding of bacterial communities. To infer the environmental quality of coastal habitats from metabarcoding data, two taxonomy-free approaches have been successfully applied for different geographical regions and monitoring goals, namely quantile regression splines (QRS) and supervised machine learning (SML). However, their comparative performance remains untested for monitoring the impact of organic enrichment introduced by aquaculture on marine coastal environments. We compared the performance of QRS and SML using bacterial metabarcoding data to infer the environmental quality of 230 aquaculture samples collected from seven farms in Norway and seven farms in Scotland along an organic enrichment gradient. As a measure of environmental quality, we used the Infaunal Quality Index (IQI) calculated from benthic macrofauna data (reference index). The QRS analysis plotted the abundance of amplicon sequence variants (ASVs) as a function to the IQI from which the ASVs with a defined abundance peak were assigned to eco-groups and a molecular IQI was subsequently calculated. In contrast, the SML approach built a random forest model to directly predict the macrofauna-based IQI. Our results show that both QRS and SML perform well in inferring the environmental quality with 89% and 90% accuracy, respectively. For both geographic regions, there was high correspondence between the reference IQI and both the inferred molecular IQIs (p < 0.001), with the SML model showing a higher coefficient of determination compared to QRS. Among the 20 most important ASVs identified by the SML approach, 15 were congruent with the good quality spline ASV indicators identified via QRS for both Norwegian and Scottish salmon farms. More research on the response of the ASVs to organic enrichment and the co-influence of other environmental parameters is necessary to eventually select the most powerful stressor-specific indicators. Even though both approaches are promising to infer environmental quality based on metabarcoding data, SML showed to be more powerful in handling the natural variability. For the improvement of the SML model, addition of new samples is still required, as background noise introduced by high spatio-temporal variability can be reduced. Overall, we recommend the development of a powerful SML approach that will be onwards applied for monitoring the impact of aquaculture on marine ecosystems based on eDNA metabarcoding data.

Geographic characteristics and environmental variables determine the diversities and assembly of the algal communities in interconnected river-lake system

Algal blooms in lakes are a major hazard worldwide. Although various geographical and environmental patterns affect algal communities during river-lake transit, a thorough understanding of what patterns shape the algal communities is still rarely researched, particularly in complex interconnected river-lake systems. In this study, focusing on the most typical interconnected river-lake system in China, the Dongting Lake, we collected paired water and sediment samples in summer, when algal biomass and growth rate are at high levels. Based on 23S rRNA gene sequencing, we investigated the heterogeneity and the differences in assembly mechanisms of planktonic and benthic algae in Dongting Lake. Planktonic algae contained more Cyanobacteria and Cryptophyta, while sediment harbored higher proportions of Bacillariophyta and Chlorophyta. For planktonic algae, stochastic dispersal dominated the assembly of the communities. Upstream rivers and confluences were important sources of planktonic algae in lakes. Meanwhile, for benthic algae, deterministic environmental filtering shaped the communities, and the proportion of benthic algae exploded with increasing N:P ratio and Cu concentration until reaching thresholds of 1.5 and 0.013 g/kg respectively, and then started falling, showing non-linear responses. This study revealed the variability of different aspects of algal communities in different habitats, traced the main sources of planktonic algae, and identified the thresholds for benthic algal shifts in response to environmental filters. Hence, upstream and downstream monitoring as well as thresholds of environmental factors should be considered in further aquatic ecological monitoring or regulatory programs of harmful algal blooms in these complex systems.

Maximising signal-to-noise ratios in environmental DNA-based monitoring

Man's impacts on global ecosystems are increasing and there is a growing demand that these activities be appropriately monitored. Monitoring requires measurement of a response metric ('signal') that changes maximally and consistently in response to the monitored activity irrespective of other factors ('noise'), thus maximising the signal-to-noise ratio. Indices derived from time-consuming morphology-based taxonomic identification of organisms are a core part of many monitoring programmes. Metabarcoding is an alternative to morphology-based identification and involves the sequencing of short fragments of DNA ('markers') from multiple taxa simultaneously. DNA suitable for metabarcoding includes that extracted from environmental samples (eDNA). Metabarcoding outputs DNA sequences that can be identified (annotated) by matching them against archived annotated sequences. However, sequences from most organisms are not archived - preventing annotation and potentially limiting metabarcoding in monitoring applications. Consequently, there is growing interest in using unannotated sequences as response metrics in monitoring programmes. We compared the sequences from three commonly used markers (16S (V3/V4 regions), 18S (V1/V2 regions) and COI) and, sampling along steep impact gradients, showed that the 16S and COI sequences were associated with the largest and smallest signal-to-noise ratio respectively. We trialled four separate, intuitive, noise-reduction approaches and demonstrated that removing less frequent sequences improved the signal-to-noise ratio, partitioning an additional 25 % from noise to explanatory factors in non-parametric ANOVA (NPA) and reducing dispersion in the data. For the 16S marker, retaining only the most frequently observed sequence, per sample, resulting in nine sequences across 150 samples, generated a near-maximal signal-to-noise ratio (95 % of the variance explained in NPA). We recommend that NPA, combined with rigorous elimination of less frequent sequences, be used to pre-filter sequences/taxa being used in monitoring applications. Our approach will simplify downstream analysis, for example the identification of key taxa and functional associations.

Environmental DNA-based profiling of benthic bacterial and eukaryote communities along a crude oil spill gradient in a coral reef in the Persian Gulf

Coral reef ecosystems in the Persian Gulf are frequently exposed to crude oil spills. We investigated benthic bacterial and eukaryote community structures at such coral reef sites subjected to different degrees of polycyclic aromatic hydrocarbon (PAH) pollution using environmental DNA (eDNA) metabarcoding. Both bacterial and eukaryote communities responded with pronounced shifts to crude oil pollution and distinguished control sites, moderately and heavily impacted sites with significant confidentiality. The observed community patterns were predominantly driven by Alphaproteobacteria and metazoans. Among these, we identified individual genera that were previously linked to oil spill stress, but also taxa, for which a link to hydrocarbon still remains to be established. Considering the lack of an early-warning system for the environmental status of coral reef ecosystems exposed to frequent crude-oil spills, our results encourage further research towards the development of an eDNA-based biomonitoring tool that exploits benthic bacterial and eukaryote communities as bioindicators.

Ultrastructure of Diophrys appendiculata and new systematic consideration of the euplotid family Uronychiidae (Protista, Ciliophora)

The ultrastructure of ciliates carries important cytological, taxonomical, and evolutionary signals for these single-celled eukaryotic organisms. However, little ultrastructural data have been accumulated for most ciliate groups with systematic problems. In the present work, a well-known marine uronychiid, Diophrys appendiculata, was investigated using electron microscopy and a comparison with, and a discussion considering, phylogenetic analyses were made. The new findings primarily show that: (i) this species lacks the typical alveolar plate, bears cortical ampule-like extrusomes, and has microtubular triads in the dorsal pellicle, and thus exhibits some ultrastructural features in common with most of its previously studied congeners; (ii) each adoral membranelle before the level of frontal cirrus II/2 contains three rows of kinetosomes and each membranelle after the level of frontal cirrus II/2 contains four rows, which might be related with morphogenesis and could be considered as a distinctive character of Diophrys; (iii) some structural details of the buccal field, such as the extra-pellicular fibrils, pellicle, pharyngeal disks and microtubular sheet, were documented. In addition, based on the ultrastructural comparison of representatives, we discuss the differentiation between the subfamilies Diophryinae and Uronychiinae. A hypothetical systematic relationship of members in the order Euplotida based on a wide range of data is also provided.

Environmental DNA metabarcoding for benthic monitoring: A review of sediment sampling and DNA extraction methods

Environmental DNA (eDNA) metabarcoding (parallel sequencing of DNA/RNA for identification of whole communities within a targeted group) is revolutionizing the field of aquatic biomonitoring. To date, most metabarcoding studies aiming to assess the ecological status of aquatic ecosystems have focused on water eDNA and macroinvertebrate bulk samples. However, the eDNA metabarcoding has also been applied to soft sediment samples, mainly for assessing microbial or meiofaunal biota. Compared to classical methodologies based on manual sorting and morphological identification of benthic taxa, eDNA metabarcoding offers potentially important advantages for assessing the environmental quality of sediments. The methods and protocols utilized for sediment eDNA metabarcoding can vary considerably among studies, and standardization efforts are needed to improve their robustness, comparability and use within regulatory frameworks. Here, we review the available information on eDNA metabarcoding applied to sediment samples, with a focus on sampling, preservation, and DNA extraction steps. We discuss challenges specific to sediment eDNA analysis, including the variety of different sources and states of eDNA and its persistence in the sediment. This paper aims to identify good-practice strategies and facilitate method harmonization for routine use of sediment eDNA in future benthic monitoring.

Molecular Data Reveal a Cryptic Diversity in the Genus Urotricha (Alveolata, Ciliophora, Prostomatida), a Key Player in Freshwater Lakes, With Remarks on Morphology, Food Preferences, and Distribution

Species of the ciliate genus Urotricha are key players in freshwater plankton communities. In the pelagial of lakes, about 20 urotrich species occur throughout an annual cycle, some of which play a pivotal role in aquatic food webs. For example, during the phytoplankton spring bloom, they consume a remarkable proportion of the algal production. In ecological studies, urotrich ciliates are usually merely identified to genus rank and grouped into size classes. This is unsatisfying considering the distinct autecological properties of individual species and their specific spatial and temporal distribution patterns. As a basis for future research, we characterized in detail four common urotrich morphotypes, i.e., specimens identified as U. furcata and tentatively as U. agilis, U. pseudofurcata, and U. castalia, using state-of-the-art methods. We used an integrative polyphasic approach, in which morphological studies (in vivo observation, silver staining methods, scanning electron microscopy) were linked with a molecular approach exploiting four different gene fragments as taxonomic DNA barcodes with different resolution potential (SSU rDNA, ITS-1, ITS-2, hypervariable V4 and V9 regions of the SSU rDNA). We shed light on the diversity of urotrich ciliates as well as on their global distribution patterns, and annual cycles. Additionally, we coupled individual species occurrences and environmental parameters, and subsequently modeled the distribution and occurrence, using logistic regressions. Furthermore, for one strain putatively identified as U. castalia, we ascertained the optimal cultivation media and food preferences. Thereby, our comprehensive view on these important freshwater ciliates that frequently occur in environmental high throughput sequencing datasets worldwide will allow future studies to better exploit protistan plankton data from lakes.

Comparing sediment preservation methods for genomic biomonitoring of coastal marine ecosystems

To avoid loss of genetic information in environmental DNA (eDNA) field samples, the preservation of nucleic acids during field sampling is a critical step. In the development of standard operating procedures (SOPs) for eDNA-based compliance monitoring, the effect of different routinely used sediment preservations on biological community structures serving as bioindicators has gone untested. We compared eDNA metabarcoding results of marine bacterial communities from sample aliquots that were treated with a nucleic acid preservation solution (treated samples) and aliquots that were frozen without further treatment (non-treated samples). Sediment samples were obtained from coastal locations subjected to different stressors (aquaculture, urbanization, industry). DNA extraction efficiency, bacterial community profiles, and measures of alpha- and beta-diversity were highly congruent between treated and non-treated samples. As both preservation methods provide the same relevant information to environmental managers and regulators, we recommend the inclusion of both methods into SOPs for biomonitoring in marine coastal environments.

High salinity gradients and intermediate spatial scales shaped similar biogeographical and co-occurrence patterns of microeukaryotes in a tropical freshwater-saltwater ecosystem

Microeukaryotes play key ecological roles in the microbial web of aquatic ecosystems. However, large knowledge gaps urgently need to be filled regarding the biogeography with associated shaping mechanisms and co-occurrence patterns of microeukaryotes under freshwater-saltwater gradients, especially true in tropical regions. Here, we investigated microeukaryotes of six mixed freshwater-saltwater regions in the Pearl River Estuary and surrounding coasts in southern China, with salinity ranging 0.1-32.0% and distances spanned up to 500 km, using molecular ecological methods. Results indicate that the biogeography of abundant and rare microeukaryotic communities was similar, both their co-occurrence patterns and biogeographical patterns were driven by deterministic and stochastic processes. The environmental factors with higher selective pressure than dispersal limitation meant that the role of deterministic process in structuring communities was more significant than that of stochastic process, and salinity played important role in structuring both microeukaryotic communities and networks. The abundant communities had stronger influence on entire microeukaryotic communities and seemed to be more sensitive to environmental changes than their rare counterparts, while rare ones had stronger interspecific relationships. Finally, the geographic scale and environmental gradients of study regions should firstly be clarified in future research on the ecological processes of microeukaryotes before conclusions are drawn.

Ecosystems monitoring powered by environmental genomics: A review of current strategies with an implementation roadmap

A decade after environmental scientists integrated high-throughput sequencing technologies in their toolbox, the genomics-based monitoring of anthropogenic impacts on the biodiversity and functioning of ecosystems is yet to be implemented by regulatory frameworks. Despite the broadly acknowledged potential of environmental genomics to this end, technical limitations and conceptual issues still stand in the way of its broad application by end-users. In addition, the multiplicity of potential implementation strategies may contribute to a perception that the routine application of this methodology is premature or "in development", hence restraining regulators from binding these tools into legal frameworks. Here, we review recent implementations of environmental genomics-based methods, applied to the biomonitoring of ecosystems. By taking a general overview, without narrowing our perspective to particular habitats or groups of organisms, this paper aims to compare, review and discuss the strengths and limitations of four general implementation strategies of environmental genomics for monitoring: (a) Taxonomy-based analyses focused on identification of known bioindicators or described taxa; (b) De novo bioindicator analyses; (c) Structural community metrics including inferred ecological networks; and (d) Functional community metrics (metagenomics or metatranscriptomics). We emphasise the utility of the three latter strategies to integrate meiofauna and microorganisms that are not traditionally utilised in biomonitoring because of difficult taxonomic identification. Finally, we propose a roadmap for the implementation of environmental genomics into routine monitoring programmes that leverage recent analytical advancements, while pointing out current limitations and future research needs.

Global Trends of Benthic Bacterial Diversity and Community Composition Along Organic Enrichment Gradients of Salmon Farms

The analysis of benthic bacterial community structure has emerged as a powerful alternative to traditional microscopy-based taxonomic approaches to monitor aquaculture disturbance in coastal environments. However, local bacterial diversity and community composition vary with season, biogeographic region, hydrology, sediment texture, and aquafarm-specific parameters. Therefore, without an understanding of the inherent variation contained within community complexes, bacterial diversity surveys conducted at individual farms, countries, or specific seasons may not be able to infer global universal pictures of bacterial community diversity and composition at different degrees of aquaculture disturbance. We have analyzed environmental DNA (eDNA) metabarcodes (V3-V4 region of the hypervariable SSU rRNA gene) of 138 samples of different farms located in different major salmon-producing countries. For these samples, we identified universal bacterial core taxa that indicate high, moderate, and low aquaculture impact, regardless of sampling season, sampled country, seafloor substrate type, or local farming and environmental conditions. We also discuss bacterial taxon groups that are specific for individual local conditions. We then link the metabolic properties of the identified bacterial taxon groups to benthic processes, which provides a better understanding of universal benthic ecosystem function(ing) of coastal aquaculture sites. Our results may further guide the continuing development of a practical and generic bacterial eDNA-based environmental monitoring approach.

Identifying the minimum amplicon sequence depth to adequately predict classes in eDNA-based marine biomonitoring using supervised machine learning

Environmental DNA metabarcoding is a powerful approach for use in biomonitoring and impact assessments. Amplicon-based eDNA sequence data are characteristically highly divergent in sequencing depth (total reads per sample) as influenced inter alia by the number of samples simultaneously analyzed per sequencing run. The random forest (RF) machine learning algorithm has been successfully employed to accurately classify unknown samples into monitoring categories. To employ RF to eDNA data, and avoid sequencing-depth artifacts, sequence data across samples are normalized using rarefaction, a process that inherently loses information. The aim of this study was to inform future sampling designs in terms of the relationship between sampling depth and RF accuracy. We analyzed three published and one new bacterial amplicon datasets, using a RF, based initially on the maximal rarefied data available (minimum mean of > 30,000 reads across all datasets) to give our baseline performance. We then evaluated the RF classification success based on increasingly rarefied datasets. We found that extreme to moderate rarefaction (50-5000 sequences per sample) was sufficient to achieve prediction performance commensurate to the full data, depending on the classification task. We did not find that the number of classification classes, data balance across classes, or the total number of sequences or samples, were associated with predictive accuracy. We identified the ability of the training data to adequately characterize the classes being mapped as the most important criterion and discuss how this finding can inform future sampling design for eDNA based biomonitoring to reduce costs and computation time.

Lake Ecosystem Robustness and Resilience Inferred from a Climate-Stressed Protistan Plankton Network

Network analyses of biological communities allow for identifying potential consequences of climate change on the resilience of ecosystems and their robustness to resist stressors. Using DNA metabarcoding datasets from a three-year-sampling (73 samples), we constructed the protistan plankton co-occurrence network of Lake Zurich, a model lake ecosystem subjected to climate change. Despite several documentations of dramatic lake warming in Lake Zurich, our study provides an unprecedented perspective by linking changes in biotic association patterns to climate stress. Water temperature belonged to the strongest environmental parameters splitting the data into two distinct seasonal networks (October–April; May–September). The expected ecological niche of phytoplankton, weakened through nutrient depletion because of permanent thermal stratification and through parasitic fungi, was occupied by the cyanobacterium Planktothrix rubescens and mixotrophic nanoflagellates. Instead of phytoplankton, bacteria and nanoflagellates were the main prey organisms associated with key predators (ciliates), which contrasts traditional views of biological associations in lake plankton. In a species extinction scenario, the warm season network emerged as more vulnerable than the cold season network, indicating a time-lagged effect of warmer winter temperatures on the communities. We conclude that climate stressors compromise lake ecosystem robustness and resilience through species replacement, richness differences, and succession as indicated by key network properties.

Expanding ecological assessment by integrating microorganisms into routine freshwater biomonitoring

Bioindication has become an indispensable part of water quality monitoring in most countries of the world, with the presence and abundance of bioindicator taxa, mostly multicellular eukaryotes, used for biotic indices. In contrast, microbes (bacteria, archaea and protists) are seldom used as bioindicators in routine assessments, although they have been recognized for their importance in environmental processes. Recently, the use of molecular methods has revealed unexpected diversity within known functional groups and novel metabolic pathways that are particularly important in energy and nutrient cycling. In various habitats, microbial communities respond to eutrophication, metals, and natural or anthropogenic organic pollutants through changes in diversity and function. In this review, we evaluated the common trends in these changes, documenting that they have value as bioindicators and can be used not only for monitoring but also for improving our understanding of the major processes in lotic and lentic environments. Current knowledge provides a solid foundation for exploiting microbial taxa, community structures and diversity, as well as functional genes, in novel monitoring programs. These microbial community measures can also be combined into biotic indices, improving the resolution of individual bioindicators. Here, we assess particular molecular approaches complemented by advanced bioinformatic analysis, as these are the most promising with respect to detailed bioindication value. We conclude that microbial community dynamics are a missing link important for our understanding of rapid changes in the structure and function of aquatic ecosystems, and should be addressed in the future environmental monitoring of freshwater ecosystems.

Aquatic food webs in deep temperate lakes: Key species establish through their autecological versatility

Microbial planktonic communities are the basis of food webs in aquatic ecosystems since they contribute substantially to primary production and nutrient recycling. Network analyses of DNA metabarcoding data sets emerged as a powerful tool to untangle the complex ecological relationships among the key players in food webs. In this study, we evaluated co-occurrence networks constructed from time-series metabarcoding data sets (12 months, biweekly sampling) of protistan plankton communities in surface layers (epilimnion) and bottom waters (hypolimnion) of two temperate deep lakes, Lake Mondsee (Austria) and Lake Zurich (Switzerland). Lake Zurich plankton communities were less tightly connected, more fragmented and had a higher susceptibility to a species extinction scenario compared to Lake Mondsee communities. We interpret these results as a lower robustness of Lake Zurich protistan plankton to environmental stressors, especially stressors resulting from climate change. In all networks, the phylum Ciliophora contributed the highest number of nodes, among them several in key positions of the networks. Associations in ciliate-specific subnetworks resembled autecological species-specific traits that indicate adaptions to specific environmental conditions. We demonstrate the strength of co-occurrence network analyses to deepen our understanding of plankton community dynamics in lakes and indicate biotic relationships, which resulted in new hypotheses that may guide future research in climate-stressed ecosystems.

Morphological, Phylogenetic and Ecophysiological Characterization of a New Ciliate, Platynematum rossellomorai n. sp. (Oligohymenophorea, Scuticociliatia), Detected in a Hypersaline Pond on Mallorca, Spain

With highly specialized morphology and unexplored functional capacities, ciliates from extreme habitats are drawing increasing attention. During a microbial investigation of a solar saltern pond (salinity 240‰) on Mallorca, Spain, a previously unknown scuticociliate, Platynematum rossellomorai n. sp. was isolated, cultured and studied using a tripartite approach consisting of a morphological description, a molecular analysis and an ecophysiological characterization. The ciliate has distinct morphological characteristics and its main diagnostic features include a large anteriorly positioned oral area (occupying almost half of the body length), two caudal cilia and a small number of somatic kineties. However, due to the most important generic feature of Cinetochilidae, the consistency of the arrangement of the adoral membranes, the ciliate is classified as a new member of the genus Platynematum. Its 18S rRNA gene sequence shows a sequence similarity of 91.0% to the closest deposited relative, Platynematum salinarum, and a phylogenetic analysis reveals a close relationship to other members of the family Cinetochilidae Perty, 1852. Growth experiments identify the ciliate as a borderline halophile, with a tolerance range between 180 and 280‰ salinity. The ciliate apparently accumulates the compatible solutes glycine betaine and ectoine to counterbalance osmotic stress, however, other osmoregulatory mechanisms are not excluded.

Supervised machine learning is superior to indicator value inference in monitoring the environmental impacts of salmon aquaculture using eDNA metabarcodes

Increasing anthropogenic impact and global change effects on natural ecosystems has prompted the development of less expensive and more efficient bioassessments methodologies. One promising approach is the integration of DNA metabarcoding in environmental monitoring. A critical step in this process is the inference of ecological quality (EQ) status from identified molecular bioindicator signatures that mirror environmental classification based on standard macroinvertebrate surveys. The most promising approaches to infer EQ from biotic indices (BI) are supervised machine learning (SML) and the calculation of indicator values (IndVal). In this study we compared the performance of both approaches using DNA metabarcodes of bacteria and ciliates as bioindicators obtained from 152 samples collected from seven Norwegian salmon farms. Results from standard macroinvertebrate-monitoring of the same samples were used as reference to compare the accuracy of both approaches. First, SML outperformed the IndVal approach to infer EQ from eDNA metabarcodes. The Random Forest (RF) algorithm appeared to be less sensitive to noisy data (a typical feature of massive environmental sequence data sets) and uneven data coverage across EQ classes (a typical feature of environmental compliance monitoring scheme) compared to a widely used method to infer IndVals for the calculation of a BI. Second, bacteria allowed for a more accurate EQ assessment than ciliate eDNA metabarcodes. For the implementation of DNA metabarcoding into routine monitoring programmes to assess EQ around salmon aquaculture cages, we therefore recommend bacterial DNA metabarcodes in combination with SML to classify EQ categories based on molecular signatures.

Improving eDNA-based protist diversity assessments using networks of amplicon sequence variants

Effective and precise grouping of highly similar sequences remains a major bottleneck in the evaluation of high-throughput sequencing datasets. Amplicon sequence variants (ASVs) offer a promising alternative that may supersede the widely used operational taxonomic units (OTUs) in environmental sequencing studies. We compared the performance of a recently developed pipeline based on the algorithm DADA2 for obtaining ASVs against a pipeline based on the algorithm SWARM for obtaining OTUs. Illumina-sequencing of 29 individual ciliate species resulted in up to 11 ASVs per species, while SWARM produced up to 19 OTUs per species. To improve the congruency between species diversity and molecular diversity, we applied sequence similarity networks (SSNs) for second-level sequence grouping into network sequence clusters (NSCs). At 100% sequence similarity in SWARM-SSNs, NSC numbers decreased from 7.9-fold overestimation without abundance filter, to 4.5-fold overestimation when an abundance filter was applied. For the DADA2-SSN approach, NSC numbers decreased from 3.5-fold to 3-fold overestimation. Rand index cluster analyses predicted best binning results between 97% and 94% sequence similarity for both DADA2-SSNs and SWARM-SSNs. Depending on the ecological questions addressed in an environmental sequencing study with protists we recommend ASVs as replacement for OTUs, best in combination with SSNs.

Redescription of Dexiotricha colpidiopsis (Kahl, 1926) Jankowski, 1964 (Ciliophora, Oligohymenophorea) from a Hot Spring in Iceland with Identification Key for Dexiotricha species

We isolated an encysted ciliate from a geothermal field in Iceland. The morphological features of this isolate fit the descriptions of Dexiotricha colpidiopsis Kahl, 1926) Jankowski, 1964 very well. These comprise body shape and size in vivo, the number of somatic kineties, and the positions of macronucleus and contractile vacuole. Using state-of-the-art taxonomic methods, the species is redescribed, including phylogenetic analyses of the small subunit ribosomal RNA (SSU rRNA) gene as molecular marker. In the phylogenetic analyses, D. colpidiopsis clusters with the three available SSU rRNA gene sequences of congeners, suggesting a monophyly of the genus Dexiotricha. Its closest relative in phylogenetic analyses is D. elliptica, which also shows a high morphological similarity. This is the first record of a Dexiotricha species from a hot spring, indicating a wide temperature tolerance of this species at least in the encysted state. The new findings on D. colpidiopsis are included in a briefly revision of the scuticociliate genus Dexiotricha and an identification key to the species.

Induction and genetic identification of a callus-like growth developed in the brown alga Fucus vesiculosus

Induction of an axenic filamentous-like callus growth from the brown algae Fucus vesiculosus is described. Different treatments were investigated in various combinations to develop axenic cultures based on identification of surface symbionts via 18S ribosomal RNA. Moreover, viability was confirmed after such processes by 2,3,5-triphenyl tetrazolium chloride assay that demonstrated an average viability of 29%, relative to nonsterilized explants. After six weeks of a phototrophic cultivation on artificial sea water-12-nitrilotriacetic acid (0.5% w/v agar), a filamentous-like callus growth was observed, which was identified genetically through its mitochondrial DNA after subculturing. Achievement of confirmed marine callus cultures might enrich old previously established blue biotechnology techniques and open new chances for cultivation of brown algae for production of good manufacturing practice-compliant bioproducts.

Seasonality of Planktonic Freshwater Ciliates: Are Analyses Based on V9 Regions of the 18S rRNA Gene Correlated With Morphospecies Counts?

Ciliates represent central nodes in freshwater planktonic food webs, and many species show pronounced seasonality, with short-lived maxima of a few dominant taxa while many being rare or ephemeral. These observations are primarily based on morphospecies counting methods, which, however, have limitations concerning the amount and volume of samples that can be processed. For high sampling frequencies at large scales, high throughput sequencing (HTS) of freshwater ciliates seems to be a promising tool. However, several studies reported large discrepancy between species abundance determinations by molecular compared to morphological means. Therefore, we compared ciliate DNA metabarcodes (V9 regions of the 18S rRNA gene) with morphospecies counts for a 3-year study (Lake Zurich, Switzerland; biweekly sampling, n = 74). In addition, we isolated, cultivated and sequenced the 18S rRNA gene of twelve selected ciliate species that served as seeds for HTS analyses. This workflow allowed for a detailed comparison of V9 data with microscopic analyses by quantitative protargol staining (QPS). The dynamics of V9 read abundances over the seasonal cycle corroborated well with morphospecies population patterns. Annual successions of rare and ephemeral species were more adequately characterized by V9 reads than by QPS. However, numbers of species specific sequence reads only partly reflected rank orders seen by counts. In contrast, biomass-based assemblage compositions showed higher similarity to V9 read numbers, probably indicating a relation between cell sizes and numbers / sizes of macronuclei (or 18S rRNA operons). Full-length 18S rRNA sequences of ciliates assigned to certain morphospecies are urgently needed for barcoding approaches as planktonic taxa are still poorly represented in public databases and the interpretation of HTS data depends on profound reference sequences. Through linking operational taxonomic units (OTUs) with known morphospecies, we can use the deep knowledge about the autecology of these species.

Morphology, Morphogenesis and Molecular Phylogeny of a New Obligate Halophile Ciliate, Schmidtiella ultrahalophila gen. nov., spec. nov. (Ciliophora, Hypotrichia) Isolated from a Volcanic Crater on Sal (Cape Verde Islands)

A new hypotrichous ciliate, Schmidtiella ultrahalophila gen. nov., spec. nov., was isolated from a solar saltern on the island of Sal, Cape Verde. The possession of only one short dorsal kinety clearly distinguishes S. ultrahalophila from other known hypotrichous genera and species. Further diagnostic characters include: a flexible and slender body, an average size of 85 × 15 μm in vivo; a bipartite adoral zone with two hypertrophied frontal adoral membranelles and nine to twelve ventral adoral membranelles; three frontal, one parabuccal, two frontoventral, two or three postoral ventral, and two or three frontoterminal cirri; and marginal cirral rows variable in number, usually one on each side. Ontogenetic data indicate the following: the frontal-ventral cirri originate from six or five anlagen; the proter inherits the parental adoral zone; the frontal and ventral cirri originate from five or six anlagen; and the marginal cirral rows and the dorsal kinety tend to originate intrakinetally. Additional marginal rows are rarely derived from de novo anlagen. Based on its morphology, morphogenesis and its SSU rRNA phylogenetic placement, the new species should be assigned to the order Sporadotrichida Fauré-Fremiet, 1961. Due to low taxon sampling, however, its exact position in this order remains enigmatic.

A fundamental difference between macrobiota and microbial eukaryotes: protistan plankton has a species maximum in the freshwater-marine transition zone of the Baltic Sea

Remane's Artenminimum at the horohalinicum is a fundamental concept in ecology to describe and explain the distribution of organisms along salinity gradients. However, a recent metadata analysis challenged this concept for protists, proposing a species maximum in brackish waters. Due to data bias, this literature-based investigation was highly discussed. Reliable data verifying or rejecting the species minimum for protists in brackish waters were critically lacking. Here, we sampled a pronounced salinity gradient along a west-east transect in the Baltic Sea and analysed protistan plankton communities using high-throughput eDNA metabarcoding. A strong salinity barrier at the upper limit of the horohalinicum and 10 psu appeared to select for significant shifts in protistan community structures, with dinoflagellates being dominant at lower salinities, and dictyochophytes and diatoms being keyplayers at higher salinities. Also in vertical water column gradients in deeper basins (Kiel Bight, Arkona and Bornholm Basin) appeared salinity as significant environmental determinant influencing alpha- and beta-diversity patterns. Importantly, alpha-diversity indices revealed species maxima in brackish waters, that is, indeed contrasting Remane's Artenminimum concept. Statistical analyses confirmed salinity as the major driving force for protistan community structuring with high significance. This suggests that macrobiota and microbial eukaryotes follow fundamentally different rules regarding diversity patterns in the transition zone from freshwater to marine waters.

Towards an eDNA metabarcode-based performance indicator for full-scale municipal wastewater treatment plants

Knowledge-driven management for wastewater treatment plant (WWTP) performance becomes increasingly important considering the globally growing production of wastewater and the rising demand of clean water supply. Even though the potential of microbial organisms (bacteria and protists) as bioindicators for WWTP performance is well known, it is far from being fully exploited for routine monitoring programs. Therefore, we here used massive sequencing of environmental (e)DNA metabarcodes from bacterial (V3-V4 region of the SSU rRNA gene) and eukaryote (V9 region of the SSU rRNA gene) communities in 21 activated sludge samples obtained from full-scale municipal WWTPs in Germany. Microbial community patterns were compared to standard WWTP operating parameters and two traditionally used WWTP performance indicators (Sludge Biotic Index and Sludge Index). Both indices showed low concordance and hardly correlated with chemical WWTP performance parameters nor did they correlate with microbial community structures. In contrast, microbial community profiles significantly correlated with WWTP performance parameters and operating conditions of the plants under study. Therefore, eDNA metabarcode profiles of whole microbial communities indicate the performance of WWTP and can provide useful information for management strategies. We here suggest a strategy for the development of an eDNA metabarcode based bioindicator system, which can be implemented in future standard monitoring programs for WWTP performance and effluent quality.

A Proposed Timescale for the Evolution of Armophorean Ciliates: Clevelandellids Diversify More Rapidly Than Metopids

Members of the class Armophorea occur in microaerophilic and anaerobic habitats, including the digestive tract of invertebrates and vertebrates. Phylogenetic kinships of metopid and clevelandellid armophoreans conflict with traditional morphology-based classifications. To reconcile their relationships and understand their morphological evolution and diversification, we utilized the molecular clock theory as well as information contained in the estimated time trees and morphology of extant taxa. The radiation of the last common ancestor of metopids and clevelandellids very likely occurred during the Paleozoic and crown diversification of the endosymbiotic clevelandellids dates back to the Mesozoic. According to diversification analyses, endosymbiotic clevelandellids have higher net diversification rates than predominantly free-living metopids. Their cladogenic success was very likely associated with sharply isolated ecological niches constituted by their hosts. Conflicts between traditional classifications and molecular phylogenies of metopids and clevelandellids very likely come from processes, leading to further diversification without extinction of ancestral lineages as well as from morphological plesiomorphies incorrectly classified as apomorphies. Our study thus suggests that diversification processes and reconstruction of ancestral morphologies improve the understanding of paraphyly which occurs in groups of organisms with an apparently long evolutionary history and when speciation prevails over extinction.

Consistent patterns of high alpha and low beta diversity in tropical parasitic and free-living protists

Tropical animals and plants are known to have high alpha diversity within forests, but low beta diversity between forests. By contrast, it is unknown whether microbes inhabiting the same ecosystems exhibit similar biogeographic patterns. To evaluate the biogeographies of tropical protists, we used metabarcoding data of species sampled in the soils of three lowland Neotropical rainforests. Taxa-area and distance-decay relationships for three of the dominant protist taxa and their subtaxa were estimated at both the OTU and phylogenetic levels, with presence-absence and abundance-based measures. These estimates were compared to null models. High local alpha and low regional beta diversity patterns were consistently found for both the parasitic Apicomplexa and the largely free-living Cercozoa and Ciliophora. Similar to animals and plants, the protists showed spatial structures between forests at the OTU and phylogenetic levels, and only at the phylogenetic level within forests. These results suggest that the biogeographies of macro- and micro-organismal eukaryotes in lowland Neotropical rainforests are partially structured by the same general processes. However, and unlike the animals and plants, the protist OTUs did not exhibit spatial structures within forests, which hinders our ability to estimate the local and regional diversity of protists in tropical forests.

Environmental status assessment using biological traits analyses and functional diversity indices of benthic ciliate communities

In this study, we tested the hypothesis that the functional diversity of benthic ciliates has high potential to monitor marine ecological status. Therefore, we investigated the spatial and temporal variation of functional diversity of benthic ciliates in the Yangtze Estuary during one year using biological traits analyses and functional diversity indices. Traits and community compositions showed clear spatial and temporal variations. Among a variety of biological traits, feeding type and body size emerged as strongest predictable variables. Functional divergence (FDiv) had an advantage over two other functional diversity indices, as well as over classical diversity measures (i.e. richness, evenness, Shannon-Wiener) to infer environmental status. Significant correlations between biological traits, FDiv and environmental variables (i.e. nutrients, temperature, salinity) suggested that functional diversity of benthic ciliates might be used as a bio-indicator in environmental status assessments. Further mandatory researches need to implement functional diversity of ciliates in routine monitoring programs were discussed.

Environmental DNA metabarcoding of benthic bacterial communities indicates the benthic footprint of salmon aquaculture

We evaluated benthic bacterial communities as bioindicators in environmental impact assessments of salmon aquaculture, a rapidly growing sector of seafood industry. Sediment samples (n=72) were collected from below salmon cages towards distant reference sites. Bacterial community profiles inferred from DNA metabarcodes were compared to reference data from standard macrofauna biomonitoring surveys of the same samples. Deltaproteobacteria were predominant in immediate vicinity of the salmon cages. Along the transect, significant shifts in bacterial community structures were observed with Gammaproteobacteria dominating the less-impacted sites. Alpha- and beta-diversity measures of bacterial communities correlated significantly with macrofauna diversity metrics and with five ecological status indices. Benthic bacterial communities mirror the reaction of macrofauna bioindicators to environmental disturbances caused by salmon farming. The implementation of bacterial eDNA metabarcoding in future Strategic Framework Directives is an alternative cost-effective high-throughput biomonitoring solution, providing a basis for management strategies in a matter of days rather than months.

Ciliate diversity and distribution patterns in the sediments of a seamount and adjacent abyssal plains in the tropical Western Pacific Ocean

Background

Benthic ciliates and the environmental factors shaping their distribution are far from being completely understood. Likewise, deep-sea systems are amongst the least understood ecosystems on Earth. In this study, using high-throughput DNA sequencing, we investigated the diversity and community composition of benthic ciliates in different sediment layers of a seamount and an adjacent abyssal plain in the tropical Western Pacific Ocean with water depths ranging between 813 m and 4566 m. Statistical analyses were used to assess shifts in ciliate communities across vertical sediment gradients and water depth.

Results

Nine out of 12 ciliate classes were detected in the different sediment samples, with Litostomatea accounting for the most diverse group, followed by Plagiopylea and Oligohymenophorea. The novelty of ciliate genetic diversity was extremely high, with a mean similarity of 93.25% to previously described sequences. On a sediment depth gradient, ciliate community structure was more similar within the upper sediment layers (0-1 and 9-10 cm) compared to the lower sediment layers (19-20 and 29-30 cm) at each site. Some unknown ciliate taxa which were absent from the surface sediments were found in deeper sediments layers. On a water depth gradient, the proportion of unique OTUs was between 42.2% and 54.3%, and that of OTUs shared by all sites around 14%. However, alpha diversity of the different ciliate communities was relatively stable in the surface layers along the water depth gradient, and about 78% of the ciliate OTUs retrieved from the surface layer of the shallowest site were shared with the surface layers of sites deeper than 3800 m. Correlation analyses did not reveal any significant effects of measured environmental factors on ciliate community composition and structure.

Conclusions

We revealed an obvious variation in ciliate community along a sediment depth gradient in the seamount and the adjacent abyssal plain and showed that water depth is a less important factor shaping ciliate distribution in deep-sea sediments unlike observed for benthic ciliates in shallow seafloors. Additionally, an extremely high genetic novelty of ciliate diversity was found in these habitats, which points to a hot spot for the discovery of new ciliate species.

Electronic supplementary material

The online version of this article (10.1186/s12866-017-1103-6) contains supplementary material, which is available to authorized users.

A Contribution to the Morphology and Phylogeny of Chlamydodon, with Three New Species from China (Ciliophora, Cyrtophoria)

Three new cyrtophorian ciliates isolated from coastal areas of China were described based on morphological and genetic data. The Chlamydodon mnemosyne-like species Chlamydodon similis sp. n. differs from its congeners mainly by its number of somatic kineties. Chlamydodon oligochaetus sp. n. is distinguished from its congeners mainly by having fewer somatic kineties, and/or an elongated body shape. Chlamydodon crassidens sp. n. is characterized mainly by an inverted triangular body shape, a posteriorly interrupted cross-striated band (5-6 μm wide), and a large cytostome. Moreover, we provided small-subunit (SSU) rDNA sequences of C. similis sp. n. and C. oligochaetus sp. n. Maximum likelihood (ML) and Bayesian inference (BI) consistently placed C. similis sp. n. as a sister to C. paramnemosyne, but showed different branching position of C. oligochaetus sp. n., which may be due to a low taxon sampling in the Chlamydodontidae and/or an insufficient resolution of the marker gene at species level.

Transition boundaries for protistan species turnover in hypersaline waters of different biogeographic regions

The identification of environmental barriers which govern species distribution is a fundamental concern in ecology. Even though salt was previously identified as a major transition boundary for micro- and macroorganisms alike, the salinities causing species turnover in protistan communities are unknown. We investigated 4.5 million high-quality protistan metabarcodes (V4 region of the SSU rDNA) obtained from 24 shallow salt ponds (salinities 4%-44%) from South America and Europe. Statistical analyses of protistan community profiles identified four salinity classes, which strongly selected for different protistan communities: 4-9%, 14-24%, 27-36% and 38-44%. The proportion of organisms unknown to science is highest in the 14-24% salinity class, showing that environments within this salinity range are an unappreciated reservoir of as yet undiscovered organisms. Distinct higher-rank taxon groups dominated in the four salinity classes in terms of diversity. As increasing salinities require different cellular responses to cope with salt, our results suggest that different evolutionary lineages of protists have evolved distinct haloadaptation strategies. Salinity appears to be a stronger selection factor for the structuring of protistan communities than geography. Yet, we find a higher degree of endemism in shallow salt ponds compared with less isolated ecosystems such as the open ocean. Thus, rules for biogeographic structuring of protistan communities are not universal, but depend on the ecosystem under consideration.

Beyond the "Code": A Guide to the Description and Documentation of Biodiversity in Ciliated Protists (Alveolata, Ciliophora)

Recent advances in molecular technology have revolutionized research on all aspects of the biology of organisms, including ciliates, and created unprecedented opportunities for pursuing a more integrative approach to investigations of biodiversity. However, this goal is complicated by large gaps and inconsistencies that still exist in the foundation of basic information about biodiversity of ciliates. The present paper reviews issues relating to the taxonomy of ciliates and presents specific recommendations for best practice in the observation and documentation of their biodiversity. This effort stems from a workshop that explored ways to implement six Grand Challenges proposed by the International Research Coordination Network for Biodiversity of Ciliates (IRCN-BC). As part of its commitment to strengthening the knowledge base that supports research on biodiversity of ciliates, the IRCN-BC proposes to populate The Ciliate Guide, an online database, with biodiversity-related data and metadata to create a resource that will facilitate accurate taxonomic identifications and promote sharing of data.

Benthic protists: the under-charted majority

Marine protist diversity inventories have largely focused on planktonic environments, while benthic protists have received relatively little attention. We therefore hypothesize that current diversity surveys have only skimmed the surface of protist diversity in marine sediments, which may harbor greater diversity than planktonic environments. We tested this by analyzing sequences of the hypervariable V4 18S rRNA from benthic and planktonic protist communities sampled in European coastal regions. Despite a similar number of OTUs in both realms, richness estimations indicated that we recovered at least 70% of the diversity in planktonic protist communities, but only 33% in benthic communities. There was also little overlap of OTUs between planktonic and benthic communities, as well as between separate benthic communities. We argue that these patterns reflect the heterogeneity and diversity of benthic habitats. A comparison of all OTUs against the Protist Ribosomal Reference database showed that a higher proportion of benthic than planktonic protist diversity is missing from public databases; similar results were obtained by comparing all OTUs against environmental references from NCBI's Short Read Archive. We suggest that the benthic realm may therefore be the world's largest reservoir of marine protist diversity, with most taxa at present undescribed.

Microbial eukaryote plankton communities of high-mountain lakes from three continents exhibit strong biogeographic patterns

Microbial eukaryotes hold a key role in aquatic ecosystem functioning. Yet, their diversity in freshwater lakes, particularly in high-mountain lakes, is relatively unknown compared with the marine environment. Low nutrient availability, low water temperature and high ultraviolet radiation make most high-mountain lakes extremely challenging habitats for life and require specific molecular and physiological adaptations. We therefore expected that these ecosystems support a plankton diversity that differs notably from other freshwater lakes. In addition, we hypothesized that the communities under study exhibit geographic structuring. Our rationale was that geographic dispersal of small-sized eukaryotes in high-mountain lakes over continental distances seems difficult. We analysed hypervariable V4 fragments of the SSU rRNA gene to compare the genetic microbial eukaryote diversity in high-mountain lakes located in the European Alps, the Chilean Altiplano and the Ethiopian Bale Mountains. Microbial eukaryotes were not globally distributed corroborating patterns found for bacteria, multicellular animals and plants. Instead, the plankton community composition emerged as a highly specific fingerprint of a geographic region even on higher taxonomic levels. The intraregional heterogeneity of the investigated lakes was mirrored in shifts in microbial eukaryote community structure, which, however, was much less pronounced compared with interregional beta-diversity. Statistical analyses revealed that on a regional scale, environmental factors are strong predictors for plankton community structures in high-mountain lakes. While on long-distance scales (>10 000 km), isolation by distance is the most plausible scenario, on intermediate scales (up to 6000 km), both contemporary environmental factors and historical contingencies interact to shift plankton community structures.

Assessing Low-Intensity Relationships in Complex Networks

Many large network data sets are noisy and contain links representing low-intensity relationships that are difficult to differentiate from random interactions. This is especially relevant for high-throughput data from systems biology, large-scale ecological data, but also for Web 2.0 data on human interactions. In these networks with missing and spurious links, it is possible to refine the data based on the principle of structural similarity, which assesses the shared neighborhood of two nodes. By using similarity measures to globally rank all possible links and choosing the top-ranked pairs, true links can be validated, missing links inferred, and spurious observations removed. While many similarity measures have been proposed to this end, there is no general consensus on which one to use. In this article, we first contribute a set of benchmarks for complex networks from three different settings (e-commerce, systems biology, and social networks) and thus enable a quantitative performance analysis of classic node similarity measures. Based on this, we then propose a new methodology for link assessment called z* that assesses the statistical significance of the number of their common neighbors by comparison with the expected value in a suitably chosen random graph model and which is a consistently top-performing algorithm for all benchmarks. In addition to a global ranking of links, we also use this method to identify the most similar neighbors of each single node in a local ranking, thereby showing the versatility of the method in two distinct scenarios and augmenting its applicability. Finally, we perform an exploratory analysis on an oceanographic plankton data set and find that the distribution of microbes follows similar biogeographic rules as those of macroorganisms, a result that rejects the global dispersal hypothesis for microbes.

Description of the Halophile Euplotes qatarensis nov. spec. (Ciliophora, Spirotrichea, Euplotida) Isolated from the Hypersaline Khor Al-Adaid Lagoon in Qatar

The morphology, ontogenesis, and phylogenetic relationships of a halophile euplotid ciliates, Euplotes qatarensis nov. spec., isolated from the Khor Al-Adaid Lagoon in Qatar were investigated based on live observation as well as protargol- and silver nitrate-impregnated methods. The new species is characterised by a combination of features: the halophile habitat, a cell size of 50-65 × 33-40 μm, seven dorsal ridges, 10 commonly sized frontoventral cirri, two widely spaced marginal cirri, 10 dorsolateral kineties, and a double silverline pattern. The morphogenesis is similar to that of its congeners: (i) the oral primordium develops hypoapokinetally and the parental oral apparatus is retained; (ii) the frontoventral-transverse field of five streaks gives rise to the frontal, ventral, and transverse cirri, but not to the cirri I/1 and the marginal cirri; (iii) the dorsal somatic ciliature develops by intrakinetal proliferation of basal bodies in two anlagen per kinety that are just anterior and posterior to the future division furrow; (iv) the caudal cirri are formed by the two rightmost dorsolateral kineties. The SSU rDNA sequence of E. qatarensis branches with full support in the Euplotopsis elegans-Euplotes nobilii-Euplotopsis raikovi clade. The closest related publicly available SSU rDNA sequence is the one of E. nobilii, with which E. qatarensis has 93.4% sequence similarity. Euplotes parawoodruffi Song & Bradbury, 1997 is transferred to the genus Euplotoides based on the absence of frontoventral cirrus VI/3.