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

Spatial and diel variations of bacterioplankton and pico-nanoeukaryote communities and potential biotic interactions during macroalgal blooms

We investigated spatial heterogeneity and diel variations in bacterioplankton and pico-nanoeukaryote communities, and potential biotic interactions at the extinction stage of the Ulva prolifera bloom in the Jiaozhou Bay, Yellow Sea. It was found that the presence of Ulva canopies significantly promoted the cell abundance of heterotrophic bacteria, raised evenness, and altered the community structure of bacterioplankton. A diel pattern was solely significant for pico-nanoeukaryote community structure. >50 % of variation in the heterotrophic bacterial abundance was accounted for by the ratio of Bacteroidota to Firmicutes, and dissolved organic nitrogen effectively explained the variations in cell abundances of phytoplankton populations. The factors representing biotic interactions frequently contributed substantially more than environmental factors in explaining the variations in diversity and community structure of both bacterioplankton and pico-nanoeukaryotes. There were higher proportions of eukaryotic pathogens compared to other marine systems, suggesting a higher ecological risk associated with the Ulva blooms.

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.

Spatiotemporal assembly and functional composition of planktonic microeukaryotic communities along productivity gradients in a subtropical lake

Microeukaryotes play crucial roles in the microbial loop of freshwater ecosystems, functioning both as primary producers and bacterivorous consumers. However, understanding the assembly of microeukaryotic communities and their functional composition in freshwater lake ecosystems across diverse environmental gradients remains limited. Here, we utilized amplicon sequencing of 18S rRNA gene and multivariate statistical analyses to examine the spatiotemporal and biogeographical patterns of microeukaryotes in water columns (at depths of 0.5, 5, and 10 m) within a subtropical lake in eastern China, covering a 40 km distance during spring and autumn of 2022. Our results revealed that complex and diverse microeukaryotic communities were dominated by Chlorophyta (mainly Chlorophyceae), Fungi, Alveolata, Stramenopiles, and Cryptophyta lineages. Species richness was higher in autumn than in spring, forming significant hump-shaped relationships with chlorophyll a concentration (Chl-a, an indicator of phytoplankton biomass). Microeukaryotic communities exhibited significant seasonality and distance-decay patterns. By contrast, the effect of vertical depth was negligible. Stochastic processes mainly influenced the assembly of microeukaryotic communities, explaining 63, 67, and 55% of community variation for spring, autumn, and both seasons combined, respectively. Trait-based functional analysis revealed the prevalence of heterotrophic and phototrophic microeukaryotic plankton with a trade-off along N:P ratio, Chl-a, and dissolved oxygen (DO) gradients. Similarly, the mixotrophic proportions were significantly and positively correlated with Chl-a and DO concentrations. Overall, our findings may provide useful insights into the assembly patterns of microeukaryotes in lake ecosystem and how their functions respond to environmental changes.

Using eRNA/eDNA metabarcoding to detect community-level impacts of nanoplastic exposure to benthic estuarine ecosystems

Plastic particles are ubiquitous in marine systems and fragment into smaller pieces, such as nanoplastics (NPs). The effects of NPs on marine organisms are of growing concern but are not well understood. Marine sediments act as a sink for many contaminants, like microplastics, and are rich habitats for benthic micro- and meiofauna which are ecologically-important components of marine food webs; however, little is known about the sensitivities of specific organisms to NPs or the effects on community diversity and composition. Utilizing molecular methods, such as metabarcoding of environmental DNA/RNA, allows for the rapid and comprehensive detection of microscopic organisms via high-throughput sequencing to assess adverse effects at the community level. The objective of this study was to use a metabarcoding approach to investigate the effects of NPs on benthic micro- and meiofaunal community diversity. Mesocosms were created with sediment cores collected from the Narrow River estuary (Rhode Island, USA) and exposed to 900 nm diameter weathered polystyrene beads at concentrations of 0.1, 1, 10, or 100 mg/kg dry weight in sediment for two weeks. Following exposure, RNA and DNA were co-extracted from the sediment, RNA was reverse-transcribed, 18S and COI markers were PCR-amplified, and amplicons were sequenced on an Illumina MiSeq. Using the 18S marker and eRNA template, increases to α-diversity and significant differences to β-diversity were observed in the highest NP exposures relative to the control. Observed differences in community composition were driven by the differential abundance of several types of protists and arthropods. Significant dose-dependent shifts in composition were observed in β-diversity Jaccard and Unweighted-Unifrac metrics with the 18S marker using the RNA template. To our knowledge, this is the first demonstration of a dose-response relationship for NPs at a community level, and it highlights the value of using community-level endpoints to assess environmental impacts of nanoparticles.

Hidden aliens: Application of digital PCR to track an exotic foraminifer across the Skagerrak (North Sea) correlates well with traditional morphospecies analysis

The problem of invasive species is a well-studied one, but knowledge of free-living unicellular eukaryotic invasive species is lacking. A potentially invasive foraminifer (Rhizaria), Nonionella sp. T1, was recently discovered in the Skagerrak and its fjords. Digital polymerase chain reaction (dPCR) was applied to track the spread of this non-indigenous species using a new dPCR assay (T1-1). The use of dPCR appears highly complementary to traditional hand picking of foraminiferal shells from the sediment, and is far less time-consuming. This study indicates that Nonionella sp. T1 has bypassed the outer Skagerrak strait, instead becoming established in Swedish west coast fjords, constituting up to half of the living foraminiferal community in fjord mouth areas. The ecology of Nonionella sp. T1 and its potential invasive impacts are still largely unknown, but it appears to be an opportunist using several energy sources such as nitrate respiration and kleptoplasty along with a possibly more efficient reproductive strategy to gain an advantage over the native foraminiferal species. Future ecological studies of Nonionella sp. T1 could be aided by dPCR and the novel Nonionella sp. T1-specific T1-1 assay.

Response of microalgae size-class structure to nutrients differences in northern Yellow Sea, China

Marine phytoplankton size-class structure affects ecological functions and shellfish culture. We use high-throughput sequencing and size-fractioned grading techniques to identify and analyze responses of phytoplankton differences in environmental variables at Donggang, northern Yellow Sea (high inorganic nitrogen (DIN)) and Changhai (low DIN) for 2021. The main environmental variables that correlate with differences in the proportional contributions of pico-, nano-, and microphytoplankton to the total phytoplankton community are inorganic phosphorus (DIP), nitrite to inorganic nitrogen ratio (NO2/dissolved inorganic nitrogen (DIN)), and ammonia nitrogen to inorganic nitrogen ratio (NH4/DIN), respectively. DIN, which contributes most to environmental differences, mainly positively correlates with changes in picophytoplankton biomass in high DIN waters. Nitrite (NO2) correlates mostly with changes in the proportional contribution of microphytoplankton in high DIN waters and nanophytoplankton in low DIN waters, and negatively correlates with changes in the biomass and proportional representation of microphytoplankton in low DIN waters. For near-shore phosphorus-limited waters, an increase in DIN may increase total microalgal biomass, but proportions of microphytoplankton may not increase; for high DIN waters, an increase in DIP may increase proportions of microphytoplankton, while for low DIN waters, an increase in DIP may preferentially increase proportions of picophytoplankton and nanophytoplankton. Picophytoplankton contributed little to the growth of two commercially cultured filter-feeding shellfish, Ruditapes philippinarum and Mizuhopecten yessoensis.

Diversity, community structure, and quantity of eukaryotic phytoplankton revealed using 18S rRNA and plastid 16S rRNA genes and pigment markers: a case study of the Pearl River Estuary

Understanding consistencies and discrepancies in characterizing diversity and quantity of phytoplankton is essential for better modeling ecosystem change. In this study, eukaryotic phytoplankton in the Pearl River Estuary, South China Sea were investigated using nuclear 18S rRNA and plastid 16S or 23S rRNA genes and pigment analysis. It was found that 18S abundance poorly explained the variations in total chlorophyll a (Chl-a). However, the ratios of log-transformed 18S abundance to Chl-a in the major phytoplankton groups were generally environment dependent, suggesting that the ratio has potential as an indicator of the physiological state of phytoplankton. The richness of 18S-based operational taxonomic units was positively correlated with the richness of 16S-based amplicon sequence variants of the whole phytoplankton community, but insignificant or weak for individual phytoplankton groups. Overall, the 18S based, rather than the 16S based, community structure had a greater similarity to pigment-based estimations. Relative to the pigment data, the proportion of haptophytes in the 18S dataset, and diatoms and cryptophytes in the 16S dataset, were underestimated. This study highlights that 18S metabarcoding tends to reflect biomass-based community organization of eukaryotic phytoplankton. Because there were lower copy numbers of plastid 16S than 18S per genome, metabarcoding of 16S probably approximates cell abundance-based community organization. Changes in biomass organization of the pigment-based community were sensitive to environmental changes. Taken together, multiple methodologies are recommended to be applied to more accurately profile the diversity and community composition of phytoplankton in natural ecosystems.

Supplementary Information

The online version contains supplementary material available at 10.1007/s42995-023-00186-x.

Evolutionary scaling of maximum growth rate with organism size

Data from nearly 1000 species reveal the upper bound to rates of biomass production achievable by natural selection across the Tree of Life. For heterotrophs, maximum growth rates scale positively with organism size in bacteria but negatively in eukaryotes, whereas for phototrophs, the scaling is negligible for cyanobacteria and weakly negative for eukaryotes. These results have significant implications for understanding the bioenergetic consequences of the transition from prokaryotes to eukaryotes, and of the expansion of some groups of the latter into multicellularity. The magnitudes of the scaling coefficients for eukaryotes are significantly lower than expected under any proposed physical-constraint model. Supported by genomic, bioenergetic, and population-genetic data and theory, an alternative hypothesis for the observed negative scaling in eukaryotes postulates that growth-diminishing mutations with small effects passively accumulate with increasing organism size as a consequence of associated increases in the power of random genetic drift. In contrast, conditional on the structural and functional features of ribosomes, natural selection has been able to promote bacteria with the fastest possible growth rates, implying minimal conflicts with both bioenergetic constraints and random genetic drift. If this extension of the drift-barrier hypothesis is correct, the interpretations of comparative studies of biological traits that have traditionally ignored differences in population-genetic environments will require revisiting.

FISHing for ciliates: Catalyzed reporter deposition fluorescence in situ hybridization for the detection of planktonic freshwater ciliates

Planktonic ciliate species form multiple trophic guilds and are central components of freshwater food webs. Progress in molecular analytical tools has opened new insight into ciliate assemblages. However, high and variable 18S rDNA copy numbers, typical for ciliates, make reliable quantification by amplicon sequencing extremely difficult. For an exact determination of abundances, the classical morphology-based quantitative protargol staining is still the method of choice. Morphotype analyses, however, are time consuming and need specific taxonomic expertise. Catalyzed reporter deposition fluorescence in situ hybridization (CARD-FISH) may represent a promising tool for the analysis of planktonic ciliates by combining molecular identification with microscopic quantification. We tested the applicability of CARD-FISH using nine cultured ciliate species. Eight species- and three genus-specific oligonucleotide probes were designed based on their 18S rRNA genes. The CARD-FISH protocol was adapted and the specificity of probes was established. We subsequently examined the precision of quantitation by CARD-FISH on single cultures and mock assemblages. Successful tests on lake water samples proved that planktonic ciliates could be identified and quantified in field samples by CARD-FISH. Double hybridizations allowed studying interspecific predator prey interactions between two ciliate species. In summary, we demonstrate that CARD-FISH with species-specific probes can facilitate studies on the population dynamics of closely related, small sized or cryptic species at high sampling frequencies.

Environmental RNA as a Tool for Marine Community Biodiversity Assessments

Microscopic organisms are often overlooked in traditional diversity assessments due to the difficulty of identifying them based on morphology. Metabarcoding is a method for rapidly identifying organisms where Environmental DNA (eDNA) is used as a template. However, legacy DNA is problematically detected from organisms no longer in the environment during sampling. Environmental RNA (eRNA), which is only produced by living organisms, can also be collected from environmental samples and used for metabarcoding. The aim of this study was to determine differences in community composition and diversity between eRNA and eDNA templates for metabarcoding. Using mesocosms containing field-collected communities from an estuary, RNA and DNA were co-extracted from sediment, libraries were prepared for two loci (18S and COI), and sequenced using an Illumina MiSeq. Results show a higher number of unique sequences detected from eRNA in both markers and higher α-diversity compared to eDNA. Significant differences between eRNA and eDNA for all β-diversity metrics were also detected. This study is the first to demonstrate community differences detected with eRNA compared to eDNA from an estuarine system and illustrates the broad applications of eRNA as a tool for assessing benthic community diversity, particularly for environmental conservation and management applications.

RNA Fluorescence in situ Hybridization (FISH) to Visualize Microbial Colonization and Infection in Caenorhabditis elegans Intestines

The intestines of wild Caenorhabditis nematodes are inhabited by a variety of microorganisms, including gut microbiome bacteria and pathogens, such as microsporidia and viruses. Because of the similarities between Caenorhabditis elegans and mammalian intestinal cells, as well as the power of the C. elegans system, this host has emerged as a model system to study host intestine-microbe interactions in vivo. While it is possible to observe some aspects of these interactions with bright-field microscopy, it is difficult to accurately classify microbes and characterize the extent of colonization or infection without more precise tools. RNA fluorescence in situ hybridization (FISH) can be used as a tool to identify and visualize microbes in nematodes from the wild or to experimentally characterize and quantify infection in nematodes infected with microbes in the lab. FISH probes, labeling the highly abundant small subunit ribosomal RNA, produce a bright signal for bacteria and microsporidian cells. Probes designed to target conserved regions of ribosomal RNA common to many species can detect a broad range of microbes, whereas targeting divergent regions of the ribosomal RNA is useful for narrower detection. Similarly, probes can be designed to label viral RNA. A protocol for RNA FISH staining with either paraformaldehyde (PFA) or acetone fixation is presented. PFA fixation is ideal for nematodes associated with bacteria, microsporidia, and viruses, whereas acetone fixation is necessary for the visualization of microsporida spores. Animals were first washed and fixed in paraformaldehyde or acetone. After fixation, FISH probes were incubated with samples to allow for the hybridization of probes to the desired target. The animals were again washed and then examined on microscope slides or using automated approaches. Overall, this FISH protocol enables detection, identification, and quantification of the microbes that inhabit the C. elegans intestine, including microbes for which there are no genetic tools available.

Seagrass Colonization Alters Diversity, Abundance, Taxonomic, and Functional Community Structure of Benthic Microbial Eukaryotes

Seagrass form high productive ecosystems in coastal environments. However, the effects of these coastal plants on the structure and function of the belowground eukaryotic microbiome remain elusive. In this study, we characterized the community of microbial eukaryotes (microeukaryotes) in both vegetated and unvegetated sediments using 18S rRNA gene amplicon sequencing and quantitative PCR. Analysis of sequencing data showed that the eelgrass (Zostera marina) colonization decreased the alpha diversity indices of benthic microeukaryotes. Apicomplexa represented an average of 83% of reads across all samples, with a higher proportion at the vegetated sites. The taxonomic community structure was significantly different between these two types of sediments, for which the concentration ofNH 4 + in sediment porewater and salinity could account. Phylogenetic analyses of long 18S rRNA genes (around 1,030 bp) indicated these apicomplexan parasites are closely related to gregarine Lecudina polymorpha. Determination of 18S rRNA gene abundances provided evidence that the eelgrass markedly promoted the biomass of the gregarine and all microeukaryotes in the seagrass-colonized sediments and confirmed that the gregarine was hosted by a polychaete species. Significantly higher gene abundances of heterotrophs and mixotrophs were found at the vegetated sites, which could be explained by the finer sediments and short supply of dissolved inorganic nitrogen, respectively. The pigmented protists were more abundant in 18S rRNA gene copies at the lower and higher pH levels than at the intermediate. Nevertheless, the fractions of heterotrophs and phototrophs in the community were significantly related to porewater N:P ratio. These results indicate that seagrass colonization significantly induces an increase in overall biomass and a decrease in diversity of benthic microeukaryotes, making them more heterotrophic. This study also highlights that the hotspot of eukaryotic parasites could be linked with the high productivity of a natural ecosystem.

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

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

Impacts of UV-C irradiation on marine biofilm community succession

Marine biofilms are diverse microbial communities and important ecological habitats forming on surfaces submerged in the ocean. Biofilm communities resist environmental disturbance, making them a nuisance to some human activities ('biofouling'). Anti-fouling solutions rarely address the underlying stability or compositional responses of these biofilms. Using bulk measurements and molecular analyses, we examined temporal and UV-C antifouling-based shifts in marine biofilms in the coastal Western North Atlantic Ocean during early fall. Over a 24-d period, bacterial communities shifted from early dominance of Gammaproteobacteria to increased proportions of Alphaproteobacteria, Bacteroidia and Acidimicrobiia. In a network analysis based on temporal covariance, Rhodobacteraceae (Alphaproteobacteria) nodes were abundant and densely connected with generally positive correlations. In the eukaryotic community, persistent algal, protistan, and invertebrate groups were observed, although consistent temporal succession was not detected. Biofilm UV-C treatment at 13 and 20 days resulted in losses of chlorophyll a and transparent exopolymer particles, indicating biomass disruption. Bacterial community shifts suggested that UV-C treatment decreased biofilm maturation rate and was associated with proportional shifts among diverse bacterial taxa. UV-C treatment was also associated with increased proportions of protists potentially involved in detritivory and parasitism. Older biofilm communities had increased resistance to UV-C, suggesting that early biofilms are more susceptible to UV-C based antifouling. The results suggest that UV-C irradiation is potentially an effective antifouling method in marine environments in terms of biomass removal and in slowing maturation. However, as they mature, biofilm communities may accumulate microbial members that are tolerant or resilient under UV-treatment. Importance Marine biofilms regulate processes from organic matter and pollutant turnover to eukaryotic settlement and growth. Biofilm growth and eukaryotic settlement interfering with human activities via growth on ship hulls, aquaculture operations, or other marine infrastructure are called 'biofouling'. There is a need to develop sustainable anti-fouling techniques by minimizing impacts to surrounding biota. We use the biofouling-antifouling framework to test hypotheses about marine biofilm succession and stability in response to disturbance, using a novel UV-C LED device. We demonstrate strong bacterial biofilm successional patterns and detect taxa potentially contributing to stability under UV-C stress. Despite UV-C-associated biomass losses and varying UV susceptibility of microbial taxa, we detected high compositional resistance among biofilm bacterial communities, suggesting decoupling of disruption in biomass and community composition following UV-C irradiation. We also report microbial covariance patterns over 24 days of biofilm growth, pointing to areas for study of microbial interactions and targeted antifouling.

Ecology of planktonic ciliates in a changing world: Concepts, methods, and challenges

Plankton ecologists ultimately focus on forecasting, both applied and environmental outcomes. We review how appreciating planktonic ciliates has become central to these predictions. We explore the 350‐year‐old canon on planktonic ciliates and examine its steady progression, which has been punctuated by conceptual insights and technological breakthroughs. By reflecting on this process, we offer suggestions as to where future leaps are needed, with an emphasis on predicting outcomes of global warming. We conclude that in terms of climate change research: (i) climatic hotspots (e.g. polar oceans) require attention; (ii) simply adding ciliate measurements to zooplankton/phytoplankton‐based sampling programs is inappropriate; (iii) elucidating the rare biosphere's functional ecology requires culture‐independent genetic methods; (iv) evaluating genetic adaptation (microevolution) and population composition shifts is required; (v) contrasting marine and freshwaters needs attention; (vi) mixotrophy needs attention; (vii) laboratory and field studies must couple automated measurements and molecular assessment of functional gene expression; (viii) ciliate trophic diversity requires appreciation; and (ix) marrying gene expression and function, coupled with climate change scenarios is needed. In short, continued academic efforts and financial support are essential to achieve the above; these will lead to understanding how ciliates will respond to climate change, providing tools for forecasting.

Dependence of genome size and copy number of rRNA gene on cell volume in dinoflagellates

Dinoflagellates are an ecologically important group of protists in aquatic environment and have evolved many unusual and enigmatic genomic features such as immense genome sizes, high repeated genes, and a large portion of hydroxymethyluracil in DNA. Although previous studies have observed positive correlations between the large subunit (LSU) rRNA gene copy number and genome size of a variety of eukaryotic organisms (e.g. higher plants and animals), or between cell volume and LSU rRNA gene copy number, and/or between genome size and cell size, which suggests a possible co-evolution among these three features in different lineages of life, it remains an open question regarding the relationships among these three parameters in dinoflagellates. For the first time, we estimated the copy numbers of the LSU rRNA gene, the genome sizes, and cell volumes within a broad range of dinoflagellates (covering 15 species of 11 genera) using single-cell qPCR-based assay (determining LSU rRNA gene copy number), FlowCAM (cell volume measurement), and ultraviolet spectrophotometry (genome size estimation). The measured copy number of LSU rRNA gene ranged from 398 ± 184 (Prorocentrum minimum) to 152,078 ± 33,555 copies•cell^-1 (Alexandrium pacificum), while the genome size and the cell volume ranged from 5.6 ± 0.2 (Karlodinium veneficum) to 853 ± 19.9 pg•cell^-1 (Pseliodinium pirum), and from 1,070 ± 225 (Kar. veneficum) to 168,474 ± 124,180 μm³ (Ps. pirum), respectively. Together with the three parameters measured in literature, there are significant positive linear correlations between LSU rRNA gene copy numbers and genome sizes, cell volumes and LSU rRNA gene copy numbers, and between genome sizes and cell volumes via comparisons of multi-model regression analyses, suggesting a dependence of genome size and rRNA gene copy number on the cell volumes of dinoflagellates. Validation of the measurement methods was conducted via comparisons between reported data in the literature and that predicted using the linear equations we obtained, and between genome size measured by flow cytometry (FCM) and ultraviolet spectrophotometry (Nanodrop). These results provide insightful understandings of dinoflagellate evolution in terms of the relationships among genomes, gene copy number, and cell volume, and of rRNA gene-based studies in intra-populational and intra-individual genetic diversity, taxonomy, and diversity assessment in the environment of dinoflagellates. The results also provide a dataset useful for reads calibration in environmental metabarcoding studies of dinoflagellates and selection of candidate species for whole genome sequencing.

A decadal perspective on north water microbial eukaryotes as Arctic Ocean sentinels

The North Water region, between Greenland and Ellesmere Island, with high populations of marine birds and mammals, is an Arctic icon. Due to climate related changes, seasonal patterns in water column primary production are changing but the implications for the planktonic microbial eukaryote communities that support the ecosystem are unknown. Here we report microbial community phenology in samples collected over 12 years (2005–2018) from July to October and analysed using high throughput 18S rRNA V4 amplicon sequencing. Community composition was tied to seasonality with summer communities more variable than distinct October communities. In summer, sentinel pan-Arctic species, including a diatom in the Chaetoceros socialis-gelidus complex and the picochlorophyte Micromonas polaris dominated phytoplankton and were summer specialists. In autumn, uncultured undescribed open water dinoflagellates were favored, and their ubiquity suggests they are sentinels of arctic autumn conditions. Despite the input of nutrients into surface waters, autumn chlorophyll concentrations remained low, refuting projected scenarios that longer ice-free seasons are synonymous with high autumn production and a diatom dominated bloom. Overall, the summer sentinel microbial taxa are persisting, and a subset oceanic dinoflagellate should be monitored for possible ecosystem shifts as later autumn ice formation becomes prevalent elsewhere.

Changes in the Community Structure of Under-Ice and Open-Water Microbiomes in Urban Lakes Exposed to Road Salts

Salinization of freshwater is increasingly observed in regions where chloride de-icing salts are applied to the roads in winter, but little is known about the effects on microbial communities. In this study, we analyzed the planktonic microbiomes of four lakes that differed in degree of urbanization, eutrophication and salinization, from an oligotrophic reference lake with no surrounding roads, to a eutrophic, salinized lake receiving runoff from a highway. We tested the hypothesis that an influence of road salts would be superimposed on the effects of season and trophic status. We evaluated the microbial community structure by 16S rRNA sequencing for Bacteria, and by four methods for eukaryotes: 16S rRNA chloroplast analysis, 18S rRNA sequencing, photosynthetic pigment analysis and microscopy. Consistent with our hypothesis, chloride and total nitrogen concentrations were among the most important statistical factors explaining the differences in taxonomic composition. These factors were positively correlated with the abundance of cryptophytes, haptophytes, and cyanobacteria. Ice-cover was also a major structuring factor, with clear differences between the winter communities and those of the open-water period. Nitrifying and methane oxidizing bacteria were more abundant in winter, suggesting the importance of anaerobic sediment processes and release of reduced compounds into the ice-covered water columns. The four methods for eukaryotic analysis provided complementary information. The 18S rRNA observations were strongly influenced by the presence of ribosome-rich ciliates, but revealed a much higher degree of taxonomic richness and greater separation of lakes, seasonal changes and potential salinity effects than the other methods.

Integrated Space-Time Dataset Reveals High Diversity and Distinct Community Structure of Ciliates in Mesopelagic Waters of the Northern South China Sea

Little is known about diversity distribution and community structure of ciliates in mesopelagic waters, especially how they are related to spatial and temporal changes. Here, an integrative approach, combining high-throughput cDNA sequencing and quantitative protargol stain, was used to analyze ciliate communities collected temporally along a transect from coastal to oceanic regions at depths ranging from the surface to 1000 m. The mesopelagic zone exhibited comparable alpha diversity to surface water which was consistent over temporal variation, with high diversity occurring at the interface with the euphotic zone. Comparison with the northeastern and the western Pacific Ocean revealed consistency of this vertical distribution of ciliates across oceanic basins. Mesopelagic ciliates harbored distinct community structure without significant seasonal differences, with the vertical variations driven largely by members of the classes Spirotrichea and Oligohymenophorea. Operational taxonomic units (OTUs) affiliated with Scuticociliatia, Astomatida and Apostomatida, members of which are known to be bacterivorous and/or commensal/parasitic species, were more abundant in mesopelagic waters than above, implying they are an important component of food webs in the mesopelagic zone. A combination of depth, geographic distance and environment shaped the ciliate communities, with depth being the most influential factor. Phylogenetic null modeling analysis further indicated that 57.1 and 33.3% of mesopelagic community variation was governed by dispersal limitation and heterogeneous selection, respectively, probably due to the marked biochemical and physical gradients down the water column. This suggests that ciliate community structure in the mesopelagic zone is mainly controlled by stochastic processes. Collectively, this study reports mesopelagic ciliates exhibited high diversity and distinct community structure across spatiotemporal scales and informs the processes mediating ciliate assembly in the mesopelagic zone. These should be fully considered in future studies to build a more comprehensive understanding of mesopelagic microbial assemblages.

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.

Further analyses of variation of ribosome DNA copy number and polymorphism in ciliates provide insights relevant to studies of both molecular ecology and phylogeny

Sequence-based approaches, such as analyses of ribosome DNA (rDNA) clone libraries and high-throughput amplicon sequencing, have been used extensively to infer evolutionary relationships and elucidate the biodiversity in microbial communities. However, recent studies demonstrate both rDNA copy number variation and intra-individual (intra-genomic) sequence variation in many organisms, which challenges the application of the rDNA-based surveys. In ciliates, an ecologically important clade of microbial eukaryotes, rDNA copy number and sequence variation are rarely studied. In the present study, we estimate the intraindividual small subunit rDNA (SSU rDNA) copy number and sequence variation in a wide range of taxa covering nine classes and 18 orders of the phylum Ciliophora. Our studies reveal that: (i) intra-individual sequence variation of SSU rDNA is ubiquitous in all groups of ciliates detected and the polymorphic level varies among taxa; (ii) there is a most common version of SSU rDNA sequence in each cell that is highly predominant and may represent the germline micronuclear template; (iii) compared with the most common version, other variant sequences differ in only 1-3 nucleotides, likely generated during macronuclear (somatic) amplification; (iv) the intra-cell sequence variation is unlikely to impact phylogenetic analyses; (v) the rDNA copy number in ciliates is highly variable, ranging from 10³ to 10⁶, with the highest record in Stentor roeselii. Overall, these analyses indicate the need for careful consideration of SSU rDNA variation in analyses of the role of ciliates in ecosystems.

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.

Vertical Distribution of Microbial Eukaryotes From Surface to the Hadal Zone of the Mariana Trench

Marine microbial eukaryotes are ubiquitous, comprised of phylogenetically diverse groups and play key roles in microbial food webs and global biogeochemical cycling. However, their vertical distribution in the deep sea has received little attention. In this study, we investigated the composition and diversity of the eukaryotes of both 0.2–3 μm and >3 μm size fractions from the surface to the hadal zone (8727 m) of the Mariana Trench using Illumina MiSeq sequencing for the 18S rDNA. The microbial eukaryotic community structure differed substantially across size fractions and depths. Operational taxonomic unit (OTU) richness in the >3 μm fraction was higher than that in the 0.2–3 μm fraction at the same depth. For the 0.2–3 μm fraction, sequences of Retaria (Rhizaria) were most abundant in the surface water (53.5%). Chrysophyceae (Stramenopiles) sequences dominated mostly in the samples from water depths below 1795 m. For the >3 μm fraction, sequences of Dinophyceae (Alveolata) were most abundant in surface waters (49.3%) and remained a significant proportion of total sequences at greater depths (9.8%, on average). Retaria sequences were abundant in samples of depths ≥1000 m. Amoebozoa and Apusozoa sequences were enriched in the hadal sample, comprising 38 and 20.4% of total sequences, respectively. Fungi (Opisthokonta) sequences were most abundant at 1759 m in both size fractions. Strong positive associations were found between Syndiniales (mainly MALV-I and MALV-II) and Retaria while negative associations were shown between MALV-II and Fungi in a co-occurrence analysis. This study compared the community structure of microbial eukaryotes in different zones in the deep sea and identified a distinct hadal community in the larger size fraction, suggesting the uniqueness of the eukaryotes in the biosphere in the Mariana Trench.

Evaluation of the Relationship Between the 18S rRNA/rDNA Ratio and Population Growth in the Marine Diatom Skeletonema tropicum via the Application of an Exogenous Nucleic Acid Standard

Ribosomal RNA (rRNA) has been regarded as a proxy for metabolic activity and population growth in microbes, but the limitations and assumptions of this approach should be better defined, particularly in eukaryotic microalgae. In this study, the 18S rRNA/rDNA ratio of a marine diatom, Skeletonema tropicum, was examined in batch and semi-continuous cultures subjected to low nitrogen and phosphorus treatments at a temperature of 20 °C. In the semi-continuous cultures, the measured 18S rRNA/rDNA ratio ranged from 4.0 × 10² to 5.0 × 10³ , and the logarithmic form of this ratio increased linearly with the population growth rate under both low nitrogen and low phosphorus conditions. In batch cultures grown under low nitrogen or low phosphorus conditions, log (rRNA/rDNA) also increased linearly with growth rate when the latter ranged between -0.4 and 1.5 day^-1 . The 18S rRNA/rDNA ratios of Skeletonema sampled from in the southern East China Sea were substantially lower than measured from laboratory cultures. Among the field samples, ratios obtained at a coastal station were higher than those obtained farther offshore. These results imply higher growth rate at the coastal station, but the influences of other factors, such as cell size and temperature, cannot be ruled out.

Sample Preservation, DNA or RNA Extraction and Data Analysis for High-Throughput Phytoplankton Community Sequencing

Phytoplankton is the basis for aquatic food webs and mirrors the water quality. Conventionally, phytoplankton analysis has been done using time consuming and partly subjective microscopic observations, but next generation sequencing (NGS) technologies provide promising potential for rapid automated examination of environmental samples. Because many phytoplankton species have tough cell walls, methods for cell lysis and DNA or RNA isolation need to be efficient to allow unbiased nucleic acid retrieval. Here, we analyzed how two phytoplankton preservation methods, three commercial DNA extraction kits and their improvements, three RNA extraction methods, and two data analysis procedures affected the results of the NGS analysis. A mock community was pooled from phytoplankton species with variation in nucleus size and cell wall hardness. Although the study showed potential for studying Lugol-preserved sample collections, it demonstrated critical challenges in the DNA-based phytoplankton analysis in overall. The 18S rRNA gene sequencing output was highly affected by the variation in the rRNA gene copy numbers per cell, while sample preservation and nucleic acid extraction methods formed another source of variation. At the top, sequence-specific variation in the data quality introduced unexpected bioinformatics bias when the sliding-window method was used for the quality trimming of the Ion Torrent data. While DNA-based analyses did not correlate with biomasses or cell numbers of the mock community, rRNA-based analyses were less affected by different RNA extraction procedures and had better match with the biomasses, dry weight and carbon contents, and are therefore recommended for quantitative phytoplankton analyses.

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.