A molecular approach to identify active microbes in environmental eukaryote clone libraries

Assessment of responses of cultured benthic foraminiferal communities to copper pollution through environmental RNA metabarcoding analysis

Copper (Cu) is an essential trace element and can accumulate in sedimentary environments. However, beyond a certain concentration threshold, Cu can pose toxicity risks to biota. Benthic foraminifera are widely distributed and exhibit characteristics such as high biodiversity, short life cycle, good preservation potential, and high pollutant sensitivity, making them good bioindicators. In this study, we investigated the responses of benthic foraminiferal communities to various Cu concentrations through environmental RNA metabarcoding under controlled laboratory conditions. Some foraminiferal genera-Glandulina, Pyrgo, Quinqueloculina, Sorites, and Vellaria-exhibited sensitivity to high Cu concentrations. Other genera-Ammonia, Operculina, Reophax, Bolivina, Glabratella, Miliammina, and Parasorites-displayed tolerance to high Cu concentrations. Moreover, high Cu levels led to reduced living benthic foraminifera abundance and interspecific interactions, albeit without inducing test deformities. Our findings provide a novel and efficient approach for exploring the responses of benthic foraminiferal communities to heavy metal pollution.

Vulnerability of Arctic Ocean microbial eukaryotes to sea ice loss

The Arctic Ocean (AO) is changing at an unprecedented rate, with ongoing sea ice loss, warming and freshening impacting the extent and duration of primary productivity over summer months. Surface microbial eukaryotes are vulnerable to such changes, but basic knowledge of the spatial variability of surface communities is limited. Here, we sampled microbial eukaryotes in surface waters of the Beaufort Sea from four contrasting environments: the Canada Basin (open ocean), the Mackenzie Trough (river-influenced), the Nuvuk region (coastal) and the under-ice system of the Canada Basin. Microbial community structure and composition varied significantly among the systems, with the most phylogenetically diverse communities being found in the more coastal systems. Further analysis of environmental factors showed potential vulnerability to change in the most specialised community, which was found in the samples taken in water immediately beneath the sea ice, and where the community was distinguished by rare species. In the context of ongoing sea ice loss, specialised ice-associated microbial assemblages may transition towards more generalist assemblages, with implications for the eventual loss of biodiversity and associated ecosystem function in the Arctic Ocean.

Microbial eukaryotic predation pressure and biomass at deep-sea hydrothermal vents

Deep-sea hydrothermal vent geochemistry shapes the foundation of the microbial food web by fueling chemolithoautotrophic microbial activity. Microbial eukaryotes (or protists) play a critical role in hydrothermal vent food webs as consumers and hosts of symbiotic bacteria, and as a nutritional source to higher trophic levels. We measured microbial eukaryotic cell abundance and predation pressure in low-temperature diffuse hydrothermal fluids at the Von Damm and Piccard vent fields along the Mid-Cayman Rise in the Western Caribbean Sea. We present findings from experiments performed under in situ pressure that show cell abundances and grazing rates higher than those done at 1 atmosphere (shipboard ambient pressure); this trend was attributed to the impact of depressurization on cell integrity. A relationship between the protistan grazing rate, prey cell abundance, and temperature of end-member hydrothermal vent fluid was observed at both vent fields, regardless of experimental approach. Our results show substantial protistan biomass at hydrothermally fueled microbial food webs, and when coupled with improved grazing estimates, suggest an important contribution of grazers to the local carbon export and supply of nutrient resources to the deep ocean.

Globally-distributed microbial eukaryotes exhibit endemism at deep-sea hydrothermal vents

Single-celled microbial eukaryotes inhabit deep-sea hydrothermal vent environments and play critical ecological roles in the vent-associated microbial food web. 18S rRNA amplicon sequencing of diffuse venting fluids from four geographically- and geochemically-distinct hydrothermal vent fields was applied to investigate community diversity patterns among protistan assemblages. The four vent fields include Axial Seamount at the Juan de Fuca Ridge, Sea Cliff and Apollo at the Gorda Ridge, all in the NE Pacific Ocean, and Piccard and Von Damm at the Mid-Cayman Rise in the Caribbean Sea. We describe species diversity patterns with respect to hydrothermal vent field and sample type, identify putative vent endemic microbial eukaryotes, and test how vent fluid geochemistry may influence microbial community diversity. At a semi-global scale, microbial eukaryotic communities at deep-sea vents were composed of similar proportions of dinoflagellates, ciliates, Rhizaria, and stramenopiles. Individual vent fields supported distinct and highly diverse assemblages of protists that included potentially endemic or novel vent-associated strains. These findings represent a census of deep-sea hydrothermal vent protistan communities. Protistan diversity, which is shaped by the hydrothermal vent environment at a local scale, ultimately influences the vent-associated microbial food web and the broader deep-sea carbon cycle.

Improve Niche Colonization and Microbial Interactions for Organohalide-Respiring-Bacteria-Mediated Remediation of Chloroethene-Contaminated Sites

Organohalide-respiring bacteria (OHRB)-mediated reductive dehalogenation is promising in in situ bioremediation of chloroethene-contaminated sites. The bioremediation efficiency of this approach is largely determined by the successful colonization of fastidious OHRB, which is highly dependent on the presence of proper growth niches and microbial interactions. In this study, based on two ecological principles (i.e., Priority Effects and Coexistence Theory), three strategies were developed to enhance niche colonization of OHRB, which were tested both in laboratory experiments and field applications: (i) preinoculation of a niche-preparing culture (NPC, being mainly constituted of fermenting bacteria and methanogens); (ii) staggered fermentation; and (iii) increased inoculation of CE40 (a Dehalococcoides-containing tetrachloroethene-to-ethene dechlorinating enrichment culture). Batch experimental results show significantly higher dechlorination efficiencies, as well as lower concentrations of volatile fatty acids (VFAs) and methane, in experimental sets with staggered fermentation and niche-preconditioning with NPC for 4 days (CE40_NPC-4) relative to control sets. Accordingly, a comparatively higher abundance of Dehalococcoides as major OHRB, together with a lower abundance of fermenting bacteria and methanogens, was observed in CE40_NPC-4 with staggered fermentation, which indicated the balanced syntrophic and competitive interactions between OHRB and other populations for the efficient dechlorination. Further experiments with microbial source tracking analyses suggested enhanced colonization of OHRB by increasing the inoculation ratio of CE40. The optimized conditions for enhanced colonization of OHRB were successfully employed for field bioremediation of trichloroethene (TCE, 0.3-1.4 mM)- and vinyl chloride (VC, ∼0.04 mM)-contaminated sites, resulting in 96.6% TCE and 99.7% VC dechlorination to ethene within 5 and 3 months, respectively. This study provides ecological principles-guided strategies for efficient bioremediation of chloroethene-contaminated sites, which may be also employed for removal of other emerging organohalide pollutants.

Spatial dynamics of active microeukaryotes along a latitudinal gradient: Diversity, assembly process, and co-occurrence relationships

Recent global warming is profoundly and increasingly influencing the Arctic ecosystem. Understanding how microeukaryote communities respond to changes in the Arctic Ocean is crucial for understanding their roles in the biogeochemical cycles of nutrients and elements. Between July 22 and August 19, 2016, during cruise ARA07, seawater samples were collected along a latitudinal transect extending from the East Sea of Korea to the central Arctic Ocean. Environmental RNA was extracted and the V4 hypervariable regions of the reverse transcribed SSU rRNA were amplified. The sequences generated by high throughput sequencing were clustered into zero-radius OTUs (ZOTUs), and the taxonomic identities of each ZOTU were assigned using SINTAX against the PR2 database. Thus, the diversity, community composition, and co-occurrence networks of size fractionated microeukaryotes were revealed. The present study found: 1) the alpha diversity of pico- and nano-sized microeukaryotes showed a latitudinal diversity gradient; 2) three distinct communities were identified, i.e., the Leg-A, Leg-B surface, and Leg-B subsurface chlorophyll a maximum (SCM) groups; 3) distinct network structure and composition were found in the three groups; and 4) water temperature was identified as the primary factor driving both the alpha and beta diversities of microeukaryotes. This study conducted a comprehensive and systematic survey of active microeukaryotes along a latitudinal gradient, elucidated the diversity, community composition, co-occurrence relationships, and community assembly processes among major microeukaryote assemblages, and will help shed more light on our understanding of the responses of microeukaryote communities to the changing Arctic Ocean.

Environmental RNA outperforms eDNA metabarcoding in assessing impact of marine pollution: A chromium-spiked mesocosm test

Environmental (e)DNA metabarcoding holds great promise for biomonitoring and ecotoxicological applications. However, few studies have compared the performance of eDNA versus eRNA metabarcoding in assessing organismal response to marine pollution, in experimental conditions. Here, we performed a chromium (Cr)-spiked mesocosm experimental test on benthic foraminiferal community to investigate the effects on species diversity by analysing both eDNA and eRNA metabarcoding data across different Cr concentrations in the sediment. Foraminiferal diversity in the eRNA data showed a significant negative correlation with the Cr concentration in the sediment, while a positive response was observed in the eDNA data. The foraminiferal OTUs exhibited a higher turnover rate in eRNA than in the eDNA-derived community. Furthermore, in the eRNA samples, OTUs abundance was significantly affected by the Cr gradient in the sediment (Pseudo-R² = 0.28, p = 0.05), while no significant trend was observed in the eDNA samples. The correlation between Cr concentration and foraminiferal diversity in eRNA datasets was stronger when the less abundant OTUs (<100 reads) were removed and the analyses were conducted exclusively on OTUs shared between eRNA and eDNA datasets. This indicates the importance of metabarcoding data filtering to capture ecological impacts, in addition to using the putatively active organisms in the eRNA dataset. The comparative analyses on foraminiferal diversity revealed that eRNA-based metabarcoding can better assess the response to heavy metal exposure in presence of subtle concentrations of the pollutant. Furthermore, our results suggest that to unlock the full potential for ecosystem assessment, eDNA and eRNA should be studied in parallel to control for potential sequence artifacts in routine ecosystem surveys.

Spatiotemporal dynamics of ammonia-oxidizing archaea and bacteria contributing to nitrification in sediments from Bohai Sea and South Yellow Sea, China

Nitrification is a central process in nitrogen cycle in the ocean. Ammonia-oxidizing archaea (AOA) and bacteria (AOB) play significant roles in ammonia oxidation which is the first and rate-limiting step in nitrification, and their differential contribution to nitrification is an important issue, attracting extensive attention. In this study, based on the quantification of archaeal and bacterial amoA gene and the measurement of potential nitrification rate (PNR), we investigated the spatiotemporal dynamics of PNRs and the amoA gene abundance and transcript abundance of aerobic ammonia oxidizers in surface sediments collected in summer and spring across ~900 km of the Bohai Sea and Yellow Sea in China. The results revealed that the contribution of AOA to nitrification was greater than that of AOB in coastal sediments, probably due to salinity and ammonia concentration. Besides, seasons had significant effect on amoA gene abundance and transcript abundance, especially for AOA, while both seasons and sea areas had significant influence on PNR of AOA and AOB. Further analysis showed complex relationships among amoA gene abundances, transcript abundances and PNRs. More importantly, both spatial (geographic distance) and environmental factors were vital in explaining the variations of ammonia-oxidizing microorganism abundances and the PNRs.

Contrasting Community Composition and Co-Occurrence Relationships of the Active Pico-Sized Haptophytes in the Surface and Subsurface Chlorophyll Maximum Layers of the Arctic Ocean in Summer

Haptophytes (Hacrobia: Haptophyta), which can perform phototrophic, phagotrophic, or mixotrophic nutritional modes, are critical for element cycling in a variety of aquatic ecosystems. However, their diversity, particularly in the changing Arctic Ocean (AO), remains largely unknown. In the present study, the biodiversity, community composition, and co-occurrence networks of pico-sized haptophytes in the surface water and subsurface chlorophyll maximum (SCM) layer of the AO were explored. Our results found higher alpha diversity estimates in the surface water compared with in the SCM based on high-throughput sequencing of haptophyte specific 18S rRNA. The community composition of the surface water was significantly different from that of the SCM, and water temperature was identified as the primary factor shaping the community compositions. Prymnesiales (mostly Chrysochromulina), uncultured Prymnesiophyceae, and Phaeocystis dominated the surface water communities, whereas Phaeocystis dominated the SCM communities, followed by Chrysochromulina, uncultured Prymnesiophyceae, and the remaining taxa. The communities of the surface water and SCM layer developed relatively independent modules in the metacommunity network. Nodes in the surface water were more closely connected to one another than those in the SCM. Network stability analysis revealed that surface water networks were more stable than SCM networks. These findings suggest that SCM communities are more susceptible to environmental fluctuations than those in surface water and that future global changes (e.g., global warming) may profoundly influence the development, persistence, and service of SCM in the AO.

Metatranscriptomic analysis reveals active microbes and genes responded to short-term Cr(VI) stress

Heavy metals have been severely polluting the environment. However, the response mechanism of microbial communities to short-term heavy metals stress remains unclear. In this study, metagenomics (MG) and metatranscriptomics (MT) was performed to observe the microbial response to short-term Cr(VI) stress. MG data showed that 99.1% of species were similar in the control and Cr(VI) treated groups. However, MT data demonstrated that 83% of the microbes were active in which 58.7% increased, while the relative abundance of 41.3% decreased after short-term Cr(VI) incubation. The MT results also revealed 9% of microbes were dormant in samples. Genes associated with oxidative stress, Cr(VI) transport, resistance, and reduction, as well as genes with unknown functions were 2-10 times upregulated after Cr(VI) treatment. To further confirm the function of unknown genes, two genes (314 and 494) were selected to detect the Cr(VI) resistance and reduction ability. The results showed that these genes significantly increased the Cr(VI) remediation ability of Escherichia coli. MT results also revealed an increase in the expression of some rare genera (at least two times) after Cr(VI) treatment, indicating these rare species played a crucial role in microbial response to short-term Cr(VI) stress. In summary, MT is an efficient way to understand the role of active and dormant microbes in specific environmental conditions.

Active bacterial and archaeal communities in coastal sediments: Biogeography pattern, assembly process and co-occurrence relationship

The biogeography of active microbial communities and the underlying mechanisms in marine sediments are important in microbial ecology but remain unclear. Here, using qPCR and high-throughput sequencing, we investigated bacterial and archaeal community abundances and activities by quantifying the abundance and expression of the 16S rRNA gene respectively, RNA-derived bacterial and archaeal community biogeography, assembly mechanisms and co-occurrence relationships in surface sediment samples from the Bohai Sea (BS), South Yellow Sea (SYS) and the north East China Sea (NECS) of the eastern Chinese marginal seas. The results revealed a higher heterogeneity of bacterial and archaeal community activities than of abundances and heterogeneous ecological functions among areas reflected by community compositions. Furthermore, clear geographic groups (i.e., the BS, SYS and NECS groups) were observed for all, abundant and rare active bacterial and archaeal communities, accompanied by significant distance-decay patterns. However, the abundant and rare taxa showed inconsistent geographic patterns. More importantly, deterministic processes played a greater role than stochastic processes in active bacterial and archaeal community assembly. The rare taxa had weaker abilities to disperse and/or adapt and more complex ecological processes than the abundant taxa. In addition, this study also showed that intertaxa competition was the dominant interaction between active bacterial and archaeal members, which could greatly contribute to dispersal limitation. Moreover, active bacterial and archaeal co-occurrence patterns showed significant distance-decay patterns, which were consistent with the community compositions.

Contrasting Community Composition of Active Microbial Eukaryotes in Melt Ponds and Sea Water of the Arctic Ocean Revealed by High Throughput Sequencing

Melt ponds (MPs), form as the result of thawing of snow and sea ice in the summer, have lower albedo than the sea ice and are thus partly responsible for the polar amplification of global warming. Knowing the community composition of MP organisms is key to understanding their roles in the biogeochemical cycles of nutrients and elements. However, the community composition of MP microbial eukaryotes has rarely been studied. In the present study, we assessed the microbial eukaryote biodiversity, community composition, and assembly processes in MPs and surface sea water (SW) using high throughput sequencing of 18S rRNA of size-fractionated samples. Alpha diversity estimates were lower in the MPs than SW across all size fractions. The community composition of MPs was significantly different from that of SW. The MP communities were dominated by members from Chrysophyceae, the ciliate classes Litostomatea and Spirotrichea, and the cercozoan groups Filosa-Thecofilosea. One open MP community was similar to SW communities, which was probably due to the advanced stage of development of the MP enabling the exchange of species between it and adjacent SW. High portions of shared species between MPs and SW may indicate the vigorous exchange of species between these two major types of environments in the Arctic Ocean. SW microbial eukaryote communities are mainly controlled by dispersal limitation whereas those of MP are mainly controlled by ecological drift.

Picoeukaryotic Diversity And Activity in the Northwestern Pacific Ocean Based on rDNA and rRNA High-Throughput Sequencing

Picoeukaryotes play an important role in the biogenic element cycle and energy flow in oligotrophic ecosystems. However, their biodiversity and activity are poorly studied in open ocean systems, such as the northwestern Pacific Ocean, which is characterized by a complex hydrological setting. Here, we investigated the diversity and activity of picoeukaryotes in the northwestern Pacific Ocean using high-throughput sequencing targeting the V9 region of 18S rDNA and rRNA. Our results showed that the DNA picoeukaryotic communities were mainly represented by Mamiellophyceae, MAST, MALV-II, Spirotrichea, Prymnesiophyceae, and MALV-I (69.33% of the total DNA reads), and the RNA communities were dominated by Spirotrichea, Mamiellophyceae, MAST, Pelagophyceae, and MALV-II (67.46% of the total RNA reads). The number of operational taxonomic units (OTUs) was significantly affected by temperature and salinity, and was decreased with the increasing nutrient concentration both in the DNA and RNA surveys. Significant differences were observed in the community composition between DNA-based and RNA-based molecular approaches, and these differences were mainly attributed to Mamiellophyceae, Spirotrichea, and Pelagophyceae. The RNA: DNA ratio was used as a proxy for relative metabolic activity of the individual OTUs. We found that the relative metabolic activities of Mamiellophyceae, Spirotrichea, and Pelagophyceae species in the northwestern Pacific Ocean were highly affected by the nutrient concentration, i.e., the NO₃ + NO₂ and SiO₃ concentration. Overall, our study shed light on picoeukaryotic diversity and distribution in the northwestern Pacific Ocean and revealed the correlation between the diversity, relative metabolic activities of marine picoeukaryotes, and the environmental factors.

Spatial Variability of Picoeukaryotic Communities in the Mariana Trench

Picoeukaryotes play prominent roles in the biogeochemical cycles in marine ecosystems. However, their molecular diversity studies have been confined in marine surface waters or shallow coastal sediments. Here, we investigated the diversity and metabolic activity of picoeukaryotic communities at depths ranging from the surface to the abyssopelagic zone in the western Pacific Ocean above the north and south slopes of the Mariana Trench. This was achieved by amplifying and sequencing the V4 region of both 18S ribosomal DNA and cDNA using Illumina HiSeq sequencing. Our study revealed: (1) Four super-groups (i.e., Alveolata, Opisthokonta, Rhizaria and Stramenopiles) dominated the picoeukaryote assemblages through the water column, although they accounted for different proportions at DNA and cDNA levels. Our data expand the deep-sea assemblages from current bathypelagic to abyssopelagic zones. (2) Using the cDNA-DNA ratio as a proxy of relative metabolic activity, the highest activity for most subgroups was usually found in the mesopelagic zone; and (3) Population shift along the vertical scale was more prominent than that on the horizontal differences, which might be explained by the sharp physicochemical gradients along the water depths. Overall, our study provides a better understanding of the diversity and metabolic activity of picoeukaryotes in water columns of the deep ocean in response to varying environmental conditions.

Considerations for incorporating real-time PCR assays into routine marine biosecurity surveillance programmes: a case study targeting the Mediterranean fanworm (Sabella spallanzanii) and club tunicate (Styela clava) 1

Molecular techniques may provide effective tools to enhance marine biosecurity surveillance. Prior to routine implementation, evidence-based consideration of their benefits and limitations is needed. In this study, we assessed the efficiency and practicality of visual diver surveys and real-time PCR assays (targeting DNA and RNA) for detecting two marine invasive species whose infestation levels varied between species and location: Sabella spallanzanii and Styela clava. Filtered water samples (n = 171) were collected in parallel with dive surveys at two locations as part of the New Zealand Marine High Risk Site Surveillance programme: Nelson Harbour (27 sites) and Waitemata Harbour (30 sites). Diver surveys resulted in a greater number of detections compared to real-time PCR: S. clava - 21 versus 5 sites in Nelson, 6 versus 1 in Auckland; S. spallanzanii - 18 versus 10 in Auckland, no detections in Nelson. Occupancy modelling derived detection probabilities for the real-time PCR for S. clava were low (14%), compared to S. spallanzanii (66%). This could be related to abundances, or species-specific differences in DNA shedding. Only one RNA sample was positive, suggesting that most detections were from extracellular DNA or non-viable fragments. While molecular methods cannot yet replace visual observations, this study shows they provide useful complementary information.

Diversity and potential activity patterns of planktonic eukaryotic microbes in a mesoeutrophic coastal area (eastern English Channel)

The diversity of planktonic eukaryotic microbes was studied at a coastal station of the eastern English Channel (EEC) from March 2011 to July 2015 (77 samples) using high throughput sequencing (454-pyrosequencing and Illumina) of the V2-V3 hypervariable region of the 18S SSU rDNA gene. Similar estimations of OTU relative abundance and taxonomic distribution for the dominant higher taxonomic groups (contributing >1% of the total number of OTUs) were observed with the two methods (Kolmogorov-Smirnov p-value = 0.22). Eight super-groups were identified throughout all samples: Alveolata, Stramenopiles, Opisthokonta, Hacrobia, Archeaplastida, Apusozoa, Rhizaria, and Amoebozoa (ordered by decreasing OTU richness). To gain further insight into microbial activity in the EEC, ribosomal RNA was extracted for samples from 2013–2015 (30 samples). Analysis of 18S rDNA and rRNA sequences led to the detection of 696 and 700 OTUs, respectively. Cluster analysis based on OTUs’ abundance indicated three major seasonal groups that were associated to spring, winter/autumn, and summer conditions. The clusters inferred from rRNA data showed a clearer seasonal representation of the community succession than the one based on rDNA. The rRNA/rDNA ratio was used as a proxy for relative cell activity. When all OTUs were considered, the average rRNA:rDNA ratio showed a linear trend around the 1:1 line, suggesting a linear relation between OTU abundance (rDNA) and activity (rRNA). However, this ratio was highly variable over time when considering individual OTUs. Interestingly, the OTU affiliated with P. globosa displayed rRNA:rDNA ratio that allowed to delimit high vs low abundance and high vs low activity periods. It unveiled quite well the Phaeocystis bloom dynamic regarding cell proliferation and activity, and could even be used as early indicator of an upcoming bloom.

Benthic Bacterial Community Composition in the Oligohaline-Marine Transition of Surface Sediments in the Baltic Sea Based on rRNA Analysis

Salinity has a strong impact on bacterial community composition such that freshwater bacterial communities are very different from those in seawater. By contrast, little is known about the composition and diversity of the bacterial community in the sediments (bacteriobenthos) at the freshwater-seawater transition (mesohaline conditions). In this study, partial 16S-rRNA sequences were used to investigate the bacterial community at five stations, representing almost freshwater (oligohaline) to marine conditions, in the Baltic Sea. Samples were obtained from the silty, top-layer (0-2.5 cm) sediments with mostly oxygenated conditions. The long water residence time characteristic of the Baltic Sea, was predicted to enable the development of autochthonous bacteriobenthos at mesohaline conditions. Our results showed that, similar to the water column, salinity is a major factor in structuring the bacteriobenthos and that there is no loss of bacterial richness at intermediate salinities. The bacterial communities of marine, mesohaline, and oligohaline sediments differed in terms of the relative rRNA abundances of the major bacterial phyla/classes. At mesohaline conditions typical marine and oligohaline operational taxonomic units (OTUs) were abundant. Putative unique OTUs in mesohaline sediments were present only at low abundances, suggesting that the mesohaline environment consists mainly of marine and oligohaline bacteria with a broad salinity tolerance. Our study provides a first overview of the diversity patterns and composition of bacteria in the sediments along the Baltic Sea salinity gradient as well as new insights into the bacteriobenthos at mesohaline conditions.

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.

Patterns and Drivers of Vertical Distribution of the Ciliate Community from the Surface to the Abyssopelagic Zone in the Western Pacific Ocean

The deep sea is one of the largest but least understood ecosystems on earth. Knowledge about the diversity and distribution patterns as well as drivers of microbial eukaryote (including ciliates) along the water column, particularly below the photic zone, is scarce. In this study, we investigated the diversity of pelagic ciliates, the main group of marine microeukaryotes, their vertical distribution from the surface to the abyssopelagic zone, as well as their horizontal distribution over a distance of 1,300 km in the Western Pacific Ocean, using high-throughput DNA and cDNA (complementary DNA) sequencing. No distance-decay relationship could be detected along the horizontal scale; instead, a distinct vertical distribution within the ciliate communities was revealed. The alpha diversity of the ciliate communities in the deep chlorophyll maximum (DCM) and the 200 m layer turned out to be significantly higher compared with the other water layers. The ciliate communities in the 200 m water layer appeared to be more similar to those in deeper layers from 1,000 m to about 5,000 m than to the surface and DCM ciliate communities. Dominant species in the bathypelagic and abyssopelagic zone, particularly some parasites, were also detected in the 200 m layer, but were almost absent in the surface layer. The 200 m layer, therefore, seems to be an important “species bank” for deep ocean layers. Statistical analyses further revealed significant effects of temperature and chlorophyll a on the partitioning of ciliate diversity, indicating that environmental factors are a stronger force in shaping marine pelagic ciliate communities than the geographic distance.

Wanted dead or alive? Using metabarcoding of environmental DNA and RNA to distinguish living assemblages for biosecurity applications

High-throughput sequencing metabarcoding studies in marine biosecurity have largely focused on targeting environmental DNA (eDNA). DNA can persist extracellularly in the environment, making discrimination of living organisms difficult. In this study, bilge water samples (i.e., water accumulating on-board a vessel during transit) were collected from 15 small recreational and commercial vessels. eDNA and eRNA molecules were co-extracted and the V4 region of the 18S ribosomal RNA gene targeted for metabarcoding. In total, 62.7% of the Operational Taxonomic Units (OTUs) were identified at least once in the corresponding eDNA and eRNA reads, with 19.5% unique to eDNA and 17.7% to eRNA. There were substantial differences in diversity between molecular compartments; 57% of sequences from eDNA-only OTUs belonged to fungi, likely originating from legacy DNA. In contrast, there was a higher percentage of metazoan (50.2%) and ciliate (31.7%) sequences in the eRNA-only OTUs. Our data suggest that the presence of eRNA-only OTUs could be due to increased cellular activities of some rare taxa that were not identified in the eDNA datasets, unusually high numbers of rRNA transcripts in ciliates, and/or artefacts produced during the reverse transcriptase, PCR and sequencing steps. The proportions of eDNA/eRNA shared and unshared OTUs were highly heterogeneous within individual bilge water samples. Multiple factors including boat type and the activities performed on-board, such as washing of scientific equipment, may play a major role in contributing to this variability. For some marine biosecurity applications analysis, eDNA-only data may be sufficient, however there are an increasing number of instances where distinguishing the living portion of a community is essential. For these circumstances, we suggest only including OTUs that are present in both eDNA and eRNA data. OTUs found only in the eRNA data need to be interpreted with caution until further research provides conclusive evidence for their origin.

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.

Microbial Eukaryote Diversity and Activity in the Water Column of the South China Sea Based on DNA and RNA High Throughput Sequencing

To study the diversity and metabolic activity of microbial eukaryotes in the water column of the South China Sea, genomic DNA and RNA were co-extracted from samples collected down to bathyal depth at two sites. V9 regions of both SSU rRNA gene and its transcript (cDNA) were amplified and sequenced using high throughput sequencing. Our study revealed: (1) DNA and RNA datasets showed significant differences in microbial eukaryote community composition, with the variability between the two datasets for the same sample exceeding that between samples within each dataset, indicating that nucleic acid source overrode environmental factors in determining the composition of microeukaryotes; (2) despite the differences in community composition between the two datasets, both DNA and RNA revealed similar depth-related distribution patterns of microbial eukaryotes; (3) using the ratio of RNA: DNA as a proxy of relative metabolic activity, a depth-related pattern was found for the relative metabolic activity of some but not all groups of microbial eukaryotes, with the highest activity for the groups with depth-related pattern usually found in the middle water layers; and (4) the presence of live and active photoautotrophic microbial eukaryotes in the deep ocean was confirmed, indicating that they play an important role in controlling the deep-sea organic carbon pool. Overall, our study sheds light on the diversity and activity of microbial eukaryotes in the water column of a tropical oligotrophic ocean and their potential contributions in the downward transportation of organic material from the surface ocean to the deep via the biological pump.

Single Cell Analysis Linking Ribosomal (r)DNA and rRNA Copy Numbers to Cell Size and Growth Rate Provides Insights into Molecular Protistan Ecology

Ribosomal (r)RNA and rDNA have been golden molecular markers in microbial ecology. However, it remains poorly understood how ribotype copy number (CN)‐based characteristics are linked with diversity, abundance, and activity of protist populations and communities observed at organismal levels. Here, we applied a single‐cell approach to quantify ribotype CNs in two ciliate species reared at different temperatures. We found that in actively growing cells, the per‐cell rDNA and rRNA CNs scaled with cell volume (CV) to 0.44 and 0.58 powers, respectively. The modeled rDNA and rRNA concentrations thus appear to be much higher in smaller than in larger cells. The observed rRNA:rDNA ratio scaled with CV^0.14. The maximum growth rate could be well predicted by a combination of per‐cell ribotype CN and temperature. Our empirical data and modeling on single‐cell ribotype scaling are in agreement with both the metabolic theory of ecology and the growth rate hypothesis, providing a quantitative framework for linking cellular rDNA and rRNA CNs with body size, growth (activity), and biomass stoichiometry. This study also demonstrates that the expression rate of rRNA genes is constrained by cell size, and favors biomass rather than abundance‐based interpretation of quantitative ribotype data in population and community ecology of protists.

Lineage-Specific and Highly Derived Gene Sequences Among Amoebozoa, Revealed by the Comparative Analysis of Transcriptomes from Twelve Amoebozoan Species

Amoebozoa represent a difficult group for traditional morphology-based taxonomy. Molecular approaches, such as gene sequencing and DNA barcoding have greatly enhanced our knowledge of the diversity of these organisms. However, metagenomic studies of Amoebozoa still did not provide as impressive results as they did among some other groups of protists. In environmental DNA surveys done on fragments of SSU rDNA gene and other traditional DNA barcodes, Amoebozoa genes normally constitute a minor part of the total gene diversity and represent only the most abundant lineages. A potential way to resolve this problem is the usage of DNA barcodes based on genes, which are unique or highly derived in this group of organisms. In the present study, we attempted to find such genes and gene families with a low level of paralogy, potentially appropriate as Amoebozoa-specific DNA barcodes. For this we re-assembled transcriptomes of 12 amoebozoan species available from the public databases and performed gene annotation and identification of orthologous genes. In our analysis Amoebozoa-specific and highly derived sequences formed 53,182 clusters of orthologs, containing from 2 to 299 proteins each. Some of these genes may be a potential target for DNA barcoding of Amoebozoa.

Haptophyte Diversity and Vertical Distribution Explored by 18S and 28S Ribosomal RNA Gene Metabarcoding and Scanning Electron Microscopy

Haptophyta encompasses more than 300 species of mostly marine pico‐ and nanoplanktonic flagellates. Our aims were to investigate the Oslofjorden haptophyte diversity and vertical distribution by metabarcoding, and to improve the approach to study haptophyte community composition, richness and proportional abundance by comparing two rRNA markers and scanning electron microscopy (SEM). Samples were collected in August 2013 at the Outer Oslofjorden, Norway. Total RNA/cDNA was amplified by haptophyte‐specific primers targeting the V4 region of the 18S, and the D1‐D2 region of the 28S rRNA. Taxonomy was assigned using curated haptophyte reference databases and phylogenetic analyses. Both marker genes showed Chrysochromulinaceae and Prymnesiaceae to be the families with highest number of Operational Taxonomic Units (OTUs), as well as proportional abundance. The 18S rRNA data set also contained OTUs assigned to eight supported and defined clades consisting of environmental sequences only, possibly representing novel lineages from family to class. We also recorded new species for the area. Comparing coccolithophores by SEM with metabarcoding shows a good correspondence with the 18S rRNA gene proportional abundances. Our results contribute to link morphological and molecular data and 28S to 18S rRNA gene sequences of haptophytes without cultured representatives, and to improve metabarcoding methodology.

First evaluation of foraminiferal metabarcoding for monitoring environmental impact from an offshore oil drilling site

At present, environmental impacts from offshore oil and gas activities are partly determined by measuring changes in macrofauna diversity. Morphological identification of macrofauna is time-consuming, expensive and dependent on taxonomic expertise. In this study, we evaluated the applicability of using foraminiferal-specific metabarcoding for routine monitoring. Sediment samples were collected along distance gradients from two oil platforms off Taranaki (New Zealand) and their physico-chemical properties, foraminiferal environmental DNA/RNA, and macrofaunal composition analyzed. Macrofaunal and foraminiferal assemblages showed similar shifts along impact gradients, but responded differently to environmental perturbations. Macrofauna were affected by hypoxia, whereas sediment grain size appeared to drive shifts in foraminifera. We identified eight foraminiferal molecular operational taxonomic units that have potential to be used as bioindicator taxa. Our results show that metabarcoding represents an effective tool for assessing foraminiferal communities near offshore oil and gas platforms, and that it can be used to complement current monitoring techniques.

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.

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.

Fungal communities from the calcareous deep-sea sediments in the Southwest India Ridge revealed by Illumina sequencing technology

The diversity and ecological significance of bacteria and archaea in deep-sea environments have been thoroughly investigated, but eukaryotic microorganisms in these areas, such as fungi, are poorly understood. To elucidate fungal diversity in calcareous deep-sea sediments in the Southwest India Ridge (SWIR), the internal transcribed spacer (ITS) regions of rRNA genes from two sediment metagenomic DNA samples were amplified and sequenced using the Illumina sequencing platform. The results revealed that 58-63 % and 36-42 % of the ITS sequences (97 % similarity) belonged to Basidiomycota and Ascomycota, respectively. These findings suggest that Basidiomycota and Ascomycota are the predominant fungal phyla in the two samples. We also found that Agaricomycetes, Leotiomycetes, and Pezizomycetes were the major fungal classes in the two samples. At the species level, Thelephoraceae sp. and Phialocephala fortinii were major fungal species in the two samples. Despite the low relative abundance, unidentified fungal sequences were also observed in the two samples. Furthermore, we found that there were slight differences in fungal diversity between the two sediment samples, although both were collected from the SWIR. Thus, our results demonstrate that calcareous deep-sea sediments in the SWIR harbor diverse fungi, which augment the fungal groups in deep-sea sediments. This is the first report of fungal communities in calcareous deep-sea sediments in the SWIR revealed by Illumina sequencing.

The influence of environmental factors on protistan microorganisms in grassland soils along a land-use gradient

In this study, we investigated the effect of land use intensity, soil parameters and vegetation on protistan communities in grassland soils. We performed qualitative (T-RFLP) and quantitative (qPCR) analyses using primers specifically targeting the 18S rRNA gene for all Eukarya and for two common flagellate groups, i.e. the Chrysophyceae and the Kinetoplastea. Both approaches were applied to extracted soil DNA and RNA, in order to distinguish between the potentially active protists (i.e. RNA pool) and the total protistan communities, including potentially inactive and encysted cells (i.e. DNA pool). Several environmental determinants such as site, soil parameters and vegetation had an impact on the T-RFLP community profiles and the abundance of the quantified 18S rRNA genes. Correlating factors often differed between quantitative (qPCR) and qualitative (T-RFLP) approaches. For instance the Chrysophyceae/Eukarya 18S rDNA ratio as determined by qPCR correlated with the C/N ratio, whereas the community composition based on T-RLFP analysis was not affected indicating that both methods taken together provide a more complete picture of the parameters driving protist diversity. Moreover, distinct T-RFs were obtained, which could serve as potential indicators for either active organisms or environmental conditions like water content. While site was the main determinant across all investigated exploratories, land use seemed to be of minor importance for structuring protist communities. The impact of other parameters differed between the target groups, e.g. Kinetoplastea reacted on changes to water content on all sites, whereas Chrysophyceae were only affected in the Schorfheide. Finally, in most cases different responses were observed on RNA- and DNA-level, respectively. Vegetation for instance influenced the two flagellate groups only at the DNA-level across all sites. Future studies should thus include different protistan groups and also distinguish between active and inactive cells, in order to reveal causal shifts in community composition and abundance in agriculturally used systems.

Combined Culture-Based and Culture-Independent Approaches Provide Insights into Diversity of Jakobids, an Extremely Plesiomorphic Eukaryotic Lineage

We used culture-based and culture-independent approaches to discover diversity and ecology of anaerobic jakobids (Excavata: Jakobida), an overlooked, deep-branching lineage of free-living nanoflagellates related to Euglenozoa. Jakobids are among a few lineages of nanoflagellates frequently detected in anoxic habitats by PCR-based studies, however only two strains of a single jakobid species have been isolated from those habitats. We recovered 712 environmental sequences and cultured 21 new isolates of anaerobic jakobids that collectively represent at least ten different species in total, from which four are uncultured. Two cultured species have never been detected by environmental, PCR-based methods. Surprisingly, culture-based and culture-independent approaches were able to reveal a relatively high proportion of overall species diversity of anaerobic jakobids—60 or 80%, respectively. Our phylogenetic analyses based on SSU rDNA and six protein-coding genes showed that anaerobic jakobids constitute a clade of morphologically similar, but genetically and ecologically diverse protists—Stygiellidae fam. nov. Our investigation combines culture-based and environmental molecular-based approaches to capture a wider extent of species diversity and shows Stygiellidae as a group that ordinarily inhabits anoxic, sulfide- and ammonium-rich marine habitats worldwide.

Large variability of bathypelagic microbial eukaryotic communities across the world's oceans

In this work, we study the diversity of bathypelagic microbial eukaryotes (0.8-20 μm) in the global ocean. Seawater samples from 3000 to 4000 m depth from 27 stations in the Atlantic, Pacific and Indian Oceans were analyzed by pyrosequencing the V4 region of the 18S ribosomal DNA. The relative abundance of the most abundant operational taxonomic units agreed with the results of a parallel metagenomic analysis, suggesting limited PCR biases in the tag approach. Although rarefaction curves for single stations were seldom saturated, the global analysis of all sequences together suggested an adequate recovery of bathypelagic diversity. Community composition presented a large variability among samples, which was poorly explained by linear geographic distance. In fact, the similarity between communities was better explained by water mass composition (26% of the variability) and the ratio in cell abundance between prokaryotes and microbial eukaryotes (21%). Deep diversity appeared dominated by four taxonomic groups (Collodaria, Chrysophytes, Basidiomycota and MALV-II) appearing in different proportions in each sample. Novel diversity amounted to 1% of the pyrotags and was lower than expected. Our study represents an essential step in the investigation of bathypelagic microbial eukaryotes, indicating dominating taxonomic groups and suggesting idiosyncratic assemblages in distinct oceanic regions.

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.

Diverse protist grazers select for virulence-related traits in Legionella

It is generally accepted that selection for resistance to grazing by protists has contributed to the evolution of Legionella pneumophila as a pathogen. Grazing resistance is becoming more generally recognized as having an important role in the ecology and evolution of bacterial pathogenesis. However, selection for grazing resistance presupposes the existence of protist grazers that provide the selective pressure. To determine whether there are protists that graze on pathogenic Legionella species, we investigated the existence of such organisms in a variety of environmental samples. We isolated and characterized diverse protists that graze on L. pneumophila and determined the effects of adding L. pneumophila on the protist community structures in microcosms made from these environmental samples. Several unrelated organisms were able to graze efficiently on L. pneumophila. The community structures of all samples were markedly altered by the addition of L. pneumophila. Surprisingly, some of the Legionella grazers were closely related to species that are known hosts for L. pneumophila, indicating the presence of unknown specificity determinants for this interaction. These results provide the first direct support for the hypothesis that protist grazers exert selective pressure on Legionella to acquire and retain adaptations that contribute to survival, and that these properties are relevant to the ability of the bacteria to cause disease in people. We also report a novel mechanism of killing of amoebae by one Legionella species that requires an intact Type IV secretion system but does not involve intracellular replication. We refer to this phenomenon as 'food poisoning'.

Protistan diversity in a permanently stratified meromictic lake (Lake Alatsee, SW Germany)

Protists play a crucial role for ecosystem function(ing) and oxygen is one of the strongest barriers against their local dispersal. However, protistan diversity in freshwater habitats with oxygen gradients received very little attention. We applied high-throughput sequencing of the V9 region (18S rRNA gene) to provide a hitherto unique spatiotemporal analysis of protistan diversity along the oxygen gradient of a freshwater meromictic lake (Lake Alatsee, SW Germany). In the mixolimnion, the communities experienced most seasonal structural changes, with Stramenopiles dominating in autumn and Dinoflagellata in summer. The suboxic interface supported the highest diversity, but only 23 OTUs95% (mainly Euglenozoa, after quality check and removal of operational taxonomic units (OTUs) with less than three sequences) were exclusively associated with this habitat. Eukaryotic communities in the anoxic monimolimnion showed the most stable seasonal pattern, with Chrysophyta and Bicosoecida being the dominant taxa. Our data pinpoint to the ecological role of the interface as a short-term 'meeting point' for protists, contributing to the coupling of the mixolimnion and the monimolimnion. Our analyses of divergent genetic diversity suggest a high degree of previously undescribed OTUs. Future research will have to reveal if this result actually points to a high number of undescribed species in such freshwater habitats.

Next-generation environmental diversity surveys of foraminifera: preparing the future

Foraminifera are commonly defined as marine testate protists, and their diversity is mainly assessed on the basis of the morphology of their agglutinated or mineralized tests. Diversity surveys based on environmental DNA (eDNA) have dramatically changed this view by revealing an unexpected diversity of naked and organic-walled lineages as well as detecting foraminiferal lineages in soil and freshwater environments. Moreover, single-cell analyses have allowed discrimination among genetically distinctive types within almost every described morphospecies. In view of these studies, the foraminiferal diversity appeared to be largely underestimated, but its accurate estimation was impeded by the low speed and coverage of a cloning-based eDNA approach. With the advent of high-throughput sequencing (HTS) technologies, these limitations disappeared in favor of exhaustive descriptions of foraminiferal diversity in numerous samples. Yet, the biases and errors identified in early HTS studies raised some questions about the accuracy of HTS data and their biological interpretation. Among the most controversial issues affecting the reliability of HTS diversity estimates are (1) the impact of technical and biological biases, (2) the sensitivity and specificity of taxonomic sequence assignment, (3) the ability to distinguish rare species, and (4) the quantitative interpretation of HTS data. Here, we document the lessons learned from previous HTS surveys and present the current advances and applications focusing on foraminiferal eDNA. We discuss the problems associated with HTS approaches and predict the future trends and avenues that hold promises for surveying foraminiferal diversity accurately and efficiently.

DNA barcoding in amoebozoa and challenges: the example of Cochliopodium

The diversity of microbial eukaryotes in general and amoeboid lineages in particular is poorly documented. Even though amoeboid lineages are among the most abundant microbes, taxonomic progress in the group has been hindered by the limitations of traditional taxonomy and technical difficultly in studying them. Studies using molecular approaches such as DNA barcoding with cytochrome oxidase I (COI) gene are slowly trickling in for Amoebozoa, and they hopefully will aid in unveiling the true diversity of the group. In this study a retrospective approach is used to test the utility of COI gene in a scale-bearing amoeba, Cochliopodium, which is morphologically well defined. A total of 126 COI sequences and 62 unique haplotypes were generated from 9 Cochliopodium species. Extensive analyses exploring effects of sequence evolution models and length of sequence on genetic diversity computations were conducted. The findings show that COI is a promising marker for Cochliopodium, except in one case where it failed to delineate two morphologically well-defined cochliopodiums. Two species delimitation approaches also recognize 8 genetic lineages out of 9 species examined. The taxonomic implications of these findings and factors that may confound COI as a barcode marker in Cochliopodium and other amoebae are discussed.

Living at the Limits: Evidence for Microbial Eukaryotes Thriving under Pressure in Deep Anoxic, Hypersaline Habitats

The advent of molecular tools in microbial ecology paved the way to exploit the diversity of microbes in extreme environments. Here, we review these tools as applied in one of the most polyextreme habitats known on our planet, namely, deep hypersaline anoxic basins (DHABs), located at ca. 3000–3500 m depth in the Eastern Mediterranean Sea. Molecular gene signatures amplified from environmental DHAB samples identified a high degree of genetic novelty, as well as distinct communities in the DHABs. Canonical correspondence analyses provided strong evidence that salinity, ion composition, and anoxia were the strongest selection factors shaping protistan community structures, largely preventing cross-colonization among the individual basins. Thus, each investigated basin represents a unique habitat (“isolated islands of evolution”), making DHABs ideal model sites to test evolutionary hypotheses. Fluorescence in situ hybridization assays using specifically designed probes revealed that the obtained genetic signatures indeed originated from indigenous polyextremophiles. Electron microscopy imaging revealed unknown ciliates densely covered with prokaryote ectosymbionts, which may enable adaptations of eukaryotes to DHAB conditions. The research reviewed here significantly advanced our knowledge on polyextremophile eukaryotes, which are excellent models for a number of biological research areas, including ecology, diversity, biotechnology, evolutionary research, physiology, and astrobiology.

Placing environmental next-generation sequencing amplicons from microbial eukaryotes into a phylogenetic context

Nucleotide positions in the hypervariable V4 and V9 regions of the small subunit (SSU)-rDNA locus are normally difficult to align and are usually removed before standard phylogenetic analyses. Yet, with next-generation sequencing data, amplicons of these regions are all that are available to answer ecological and evolutionary questions that rely on phylogenetic inferences. With ciliates, we asked how inclusion of the V4 or V9 regions, regardless of alignment quality, affects tree topologies using distinct phylogenetic methods (including PairDist that is introduced here). Results show that the best approach is to place V4 amplicons into an alignment of full-length Sanger SSU-rDNA sequences and to infer the phylogenetic tree with RAxML. A sliding window algorithm as implemented in RAxML shows, though, that not all nucleotide positions in the V4 region are better than V9 at inferring the ciliate tree. With this approach and an ancestral-state reconstruction, we use V4 amplicons from European nearshore sampling sites to infer that rather than being primarily terrestrial and freshwater, colpodean ciliates may have repeatedly transitioned from terrestrial/freshwater to marine environments.

Capillary electrophoresis finger print technique (CE-SSCP): an alternative tool for the monitoring activities of HAB species in Baja California Sur Costal

In Mexican waters, there is no a formal and well-established monitoring program of harmful algal blooms (HAB) events. Until now, most of the work has been focused on the characterization of organisms present in certain communities. Therefore, the development of new techniques for the rapid detection of HAB species is necessary. Capillary electrophoresis finger print technique (CE-SSCP) is a fingerprinting technique based on the identification of different conformers dependent of its base composition. This technique, coupled with capillary electrophoresis, has been used to compare and identify different conformers. The aim of this study was to determine if CE-SSCP analysis of ribosomal RNA (rRNA) gene fragments could be used for a rapid identification of toxic and harmful HAB species to improve monitoring activities along the coasts of Baja California Sur, Mexico.Three different highly variable regions of the 18S and 28S rRNA genes were chosen and their suitability for the discrimination of different dinoflagellate species was assessed by CE-SSCP.The CE-SSCP results obtained for the LSU D7 fragment has demonstrated that this technique with this gene region could be useful for the identification of the ten dinoflagellates species of different genera.We have shown that this method can be used to discriminate species and the next step will be to apply it to natural samples to achieve our goal of molecular monitoring for toxic algae in Mexican waters. This strategy will offer an option to improve an early warning system of HAB events for coastal BCS, allowing the possible implementation of mitigation strategies. A monitoring program of HAB species using molecular methods will permit the analysis of several samples in a short period of time, without the pressure of counting with a taxonomic expert in phytoplankton taxonomy.

Fungal diversity from various marine habitats deduced through culture-independent studies

Studies on the molecular diversity of the micro-eukaryotic community have shown that fungi occupy a central position in a large number of marine habitats. Environmental surveys using molecular tools have shown the presence of fungi from a large number of marine habitats such as deep-sea habitats, pelagic waters, coastal regions, hydrothermal vent ecosystem, anoxic habitats, and ice-cold regions. This is of interest to a variety of research disciplines like ecology, evolution, biogeochemistry, and biotechnology. In this review, we have summarized how molecular tools have helped to broaden our understanding of the fungal diversity in various marine habitats. Majority of the environmental phylotypes could be grouped as novel clades within Ascomycota, Basidiomycota, and Chytridiomycota or as basal fungal lineages. Deep-branching novel environmental clusters could be grouped within Ascomycota as the Pezizomycotina clone group, deep-sea fungal group-I, and soil clone group-I, within Basidiomycota as the hydrothermal and/or anaerobic fungal group, and within Chytridiomycota as Cryptomycota or the Rozella clade. However, a basal true marine environmental cluster is still to be identified as most of the clusters include representatives from terrestrial regions. The challenge for future research is to explore the true marine fungi using molecular techniques.

Reticulamoeba is a long-branched Granofilosean (Cercozoa) that is missing from sequence databases

We sequenced the 18S ribosomal RNA gene of seven isolates of the enigmatic marine amoeboflagellate Reticulamoeba Grell, which resolved into four genetically distinct Reticulamoeba lineages, two of which correspond to R. gemmipara Grell and R. minor Grell, another with a relatively large cell body forming lacunae, and another that has similarities to both R. minor and R. gemmipara but with a greater propensity to form cell clusters. These lineages together form a long-branched clade that branches within the cercozoan class Granofilosea (phylum Cercozoa), showing phylogenetic affinities with the genus Mesofila. The basic morphology of Reticulamoeba is a roundish or ovoid cell with a more or less irregular outline. Long and branched reticulopodia radiate from the cell. The reticulopodia bear granules that are bidirectionally motile. There is also a biflagellate dispersal stage. Reticulamoeba is frequently observed in coastal marine environmental samples. PCR primers specific to the Reticulamoeba clade confirm that it is a frequent member of benthic marine microbial communities, and is also found in brackish water sediments and freshwater biofilm. However, so far it has not been found in large molecular datasets such as the nucleotide database in NCBI GenBank, metagenomic datasets in Camera, and the marine microbial eukaryote sampling and sequencing consortium BioMarKs, although closely related lineages can be found in some of these datasets using a highly targeted approach. Therefore, although such datasets are very powerful tools in microbial ecology, they may, for several methodological reasons, fail to detect ecologically and evolutionary key lineages.

Evaluation of DNA extraction methods for freshwater eukaryotic microalgae

The use of molecular methods to investigate microalgal communities of natural and engineered freshwater resources is in its infancy, with the majority of previous studies carried out by microscopy. Inefficient or differential DNA extraction of microalgal community members can lead to bias in downstream community analysis. Three commercially available DNA extraction kits have been tested on a range of pure culture freshwater algal species with diverse cell walls and mixed algal cultures taken from eutrophic waste stabilization ponds (WSP). DNA yield and quality were evaluated, along with DNA suitability for amplification of 18S rRNA gene fragments by polymerase chain reaction (PCR). QiagenDNeasy(®) Blood and Tissue kit (QBT), was found to give the highest DNA yields and quality. Denaturant Gradient Gel Electrophoresis (DGGE) was used to assess the diversity of communities from which DNA was extracted. No significant differences were found among kits when assessing diversity. QBT is recommended for use with WSP samples, a conclusion confirmed by further testing on communities from two tropical WSP systems. The fixation of microalgal samples with ethanol prior to DNA extraction was found to reduce yields as well as diversity and is not recommended.

Unveiling trophic functions of uncultured protist taxa by incubation experiments in the brackish Baltic Sea

BACKGROUND

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

METHODOLOGY/PRINCIPAL FINDINGS

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

CONCLUSIONS/SIGNIFICANCE

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