Zhenfeng Liu

Investigating Microbial Eukaryotic Diversity from a Global Census: Insights from a Comparison of Pyrotag and Full-Length Sequences of 18S rRNA Genes

ABSTRACT Next-generation DNA sequencing (NGS) approaches are rapidly surpassing Sanger sequencing for characterizing the diversity of natural microbial communities. Despite this rapid transition, few comparisons exist between Sanger sequences and the generally much shorter reads of NGS. Operational taxonomic units (OTUs) derived from full-length (Sanger sequencing) and pyrotag (454 sequencing of the V9 hypervariable region) sequences of 18S rRNA genes from 10 global samples were analyzed in order to compare the resulting protistan community structures and species richness. Pyrotag OTUs called at 98% sequence similarity yielded numbers of OTUs that were similar overall to those for full-length sequences when the latter were called at 97% similarity. Singleton OTUs strongly influenced estimates of species richness but not the higher-level taxonomic composition of the community. The pyrotag and full-length sequence data sets had slightly different taxonomic compositions of rhizarians, stramenopiles, cryptophytes, and haptophytes, but the two data sets had similarly high compositions of alveolates. Pyrotag-based OTUs were often derived from sequences that mapped to multiple full-length OTUs at 100% similarity. Thus, pyrotags sequenced from a single hypervariable region might not be appropriate for establishing protistan species-level OTUs. However, nonmetric multidimensional scaling plots constructed with the two data sets yielded similar clusters, indicating that beta diversity analysis results were similar for the Sanger and NGS sequences. Short pyrotag sequences can provide holistic assessments of protistan communities, although care must be taken in interpreting the results. The longer reads (>500 bp) that are now becoming available through NGS should provide powerful tools for assessing the diversity of microbial eukaryotic assemblages.

Protistan diversity and activity inferred from RNA and DNA at a coastal ocean site in the eastern North Pacific.

Microbial eukaryotes fulfill key ecological positions in marine food webs. Molecular approaches that connect protistan diversity and biogeography to their diverse metabolisms will greatly improve our understanding of marine ecosystem function. The majority of molecular-based studies to date use 18S rRNA gene sequencing to characterize natural microbial assemblages, but this approach does not necessarily discriminate between active and non-active cells. We incorporated RNA sequencing into standard 18S rRNA gene sequence surveys with the purpose of assessing those members of the protistan community contributing to biogeochemical cycling (active organisms), using the ratio of cDNA (reverse transcribed from total RNA) to 18S rRNA gene sequences within major protistan taxonomic groups. Trophically important phytoplankton, such as diatoms and chlorophytes exhibited seasonal trends in relative activity. Additionally, both radiolaria and ciliates displayed previously unreported high relative activities below the euphotic zone. This study sheds new light on the relative metabolic activity of specific protistan groups and how microbial communities respond to changing environmental conditions.

Changes in gene expression of Prymnesium parvum induced by nitrogen and phosphorus limitation

Prymnesium parvum is a globally distributed prymnesiophyte alga commonly found in brackish water marine ecosystems and lakes. It possesses a suite of toxins with ichthyotoxic, cytotoxic and hemolytic effects which, along with its mixotrophic nutritional capabilities, allows it to form massive Ecosystem Disruptive Algal Blooms (EDABs). While blooms of high abundance coincide with high levels of nitrogen (N) and phosphorus (P), reports of field and laboratory studies have noted that P. parvum toxicity appears to be augmented at high N:P ratios or P-limiting conditions. Here we present the results of a comparative analysis of P. parvum RNA-Seq transcriptomes under nutrient replete conditions, and N or P deficiency to understand how this organism responds at the transcriptional level to varying nutrient conditions. In nutrient limited conditions we found diverse transcriptional responses for genes involved in nutrient uptake, protein synthesis and degradation, photosynthesis, and toxin production. As anticipated, when either N or P was limiting, transcription levels of genes encoding transporters for the respective nutrient were higher than those under replete condition. Ribosomal and lysosomal protein genes were expressed at higher levels under either nutrient-limited condition compared to the replete condition. Photosynthesis genes and polyketide synthase genes were more highly expressed under P-limitation but not under N-limitation. These results highlight the ability of P. parvum to mount a coordinated and varied cellular and physiological response to nutrient limitation. Results also provide potential marker genes for further evaluating the physiological response and toxin production of P. parvum populations during bloom formation or to changing environmental conditions.

Estimating Protistan Diversity Using High‐Throughput Sequencing

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

Gene expression in the mixotrophic prymnesiophyte, Prymnesium parvum, responds to prey availability

The mixotrophic prymnesiophyte, Prymnesium parvum, is a widely distributed alga with significant ecological importance. It produces toxins and can form ecosystem disruptive blooms that result in fish kills and changes in planktonic food web structure. However, the relationship between P. parvum and its prey on the molecular level is poorly understood. In this study, we used RNA-Seq technology to study changes in gene transcription of P. parvum in three treatments with different microbial populations available as potential prey: axenic P. parvum (no prey), bacterized P. paruvm, and axenic P. parvum with ciliates added as prey. Thousands of genes were differentially expressed among the three treatments. Most notably, transcriptome data indicated that P. parvum obtained organic carbon, including fatty acids, from both bacteria and ciliate prey for energy and cellular building blocks. The data also suggested that different prey provided P. parvum with macro- and micro-nutrients, namely organic nitrogen in the form of amino acids from ciliates, and iron from bacteria. However, both transcriptomic data and growth experiments indicated that P. parvum did not grow faster in the presence of prey despite the gains in nutrients, although algal abundances attained in culture were slightly greater in the presence of prey. The relationship between phototrophy, heterotrophy and growth of P. parvum is discussed.

Comparative Transcriptome Analysis of Four Prymnesiophyte Algae

Genomic studies of bacteria, archaea and viruses have provided insights into the microbial world by unveiling potential functional capabilities and molecular pathways. However, the rate of discovery has been slower among microbial eukaryotes, whose genomes are larger and more complex. Transcriptomic approaches provide a cost-effective alternative for examining genetic potential and physiological responses of microbial eukaryotes to environmental stimuli. In this study, we generated and compared the transcriptomes of four globally-distributed, bloom-forming prymnesiophyte algae: Prymnesium parvum, Chrysochromulina brevifilum, Chrysochromulina ericina and Phaeocystis antarctica. Our results revealed that the four transcriptomes possess a set of core genes that are similar in number and shared across all four organisms. The functional classifications of these core genes using the euKaryotic Orthologous Genes (KOG) database were also similar among the four study organisms. More broadly, when the frequencies of different cellular and physiological functions were compared with other protists, the species clustered by both phylogeny and nutritional modes. Thus, these clustering patterns provide insight into genomic factors relating to both evolutionary relationships as well as trophic ecology. This paper provides a novel comparative analysis of the transcriptomes of ecologically important and closely related prymnesiophyte protists and advances an emerging field of study that uses transcriptomics to reveal ecology and function in protists.

Gene expression characterizes different nutritional strategies among three mixotrophic protists

Mixotrophic protists, i.e. protists that can carry out both phototrophy and heterotrophy, are a group of organisms with a wide range of nutritional strategies. The ecological and biogeochemical importance of these species has recently been recognized. In this study, we investigated and compared the gene expression of three mixotrophic protists, Prymnesium parvum, Dinobyron sp. and Ochromonas sp. under light and dark conditions in the presence of prey using RNA-Seq. Gene expression of the obligately phototrophic P. parvum and Dinobryon sp. changed significantly between light and dark treatments, while that of primarily heterotrophic Ochromonas sp. was largely unchanged. Gene expression of P. parvum and Dinobryon sp. shared many similarities, especially in the expression patterns of genes related to reproduction. However, key genes involved in central carbon metabolism and phagotrophy had different expression patterns between these two species, suggesting differences in prey consumption and heterotrophic nutrition in the dark. Transcriptomic data also offered clues to other physiological traits of these organisms such as preference of nitrogen sources and photo-oxidative stress. These results provide potential target genes for further exploration of the mechanisms of mixotrophic physiology and demonstrate the potential usefulness of molecular approaches in characterizing the nutritional modes of mixotrophic protists.

Single-cell transcriptomics of small microbial eukaryotes: limitations and potential

Single-cell transcriptomics is an emerging research tool that has huge untapped potential in the study of microbial eukaryotes. Its application has been tested in microbial eukaryotes 50u2009μm or larger, and it generated transcriptomes similar to those obtained from culture-based RNA-seq. However, microbial eukaryotes have a wide range of sizes and can be as small as 1u2009μm. Single-cell RNA-seq was tested in two smaller protists (8 and 15u2009μm). Transcript recovery rate was much lower and randomness in observed gene expression levels was much higher in single-cell transcriptomes than those derived from bulk cultures of cells. We found that the reason of such observation is that the smaller organisms had much lower mRNA copy numbers. We discuss the application of single-cell RNA-seq in studying smaller microbial eukaryotes in the context of these limitations.

Effect of light and prey availability on gene expression of the mixotrophic chrysophyte, Ochromonas sp.

BackgroundOchromonas is a genus of mixotrophic chrysophytes that is found ubiquitously in many aquatic environments. Species in this genus can be important consumers of bacteria but vary in their ability to perform photosynthesis. We studied the effect of light and bacteria on growth and gene expression of a predominantly phagotrophic Ochromonas species. Axenic cultures of Ochromonas sp. were fed with heat-killed bacteria (HKB) and grown in constant light or darkness. RNA was extracted from cultures in the light or in the dark with HKB present (Lightu2009+u2009HKB; Darku2009+u2009HKB), and in the light after HKB were depleted (Lightu2009+u2009depleted HKB).ResultsThere were no significant differences in the growth or bacterial ingestion rates between algae grown in light or dark conditions. The availability of light led to a differential expression of only 8% of genes in the transcriptome. A number of genes associated with photosynthesis, phagotrophy, and tetrapyrrole synthesis was upregulated in the Lightu2009+u2009HKB treatment compared to Darku2009+u2009HKB. Conversely, the comparison between the Lightu2009+u2009HKB and Lightu2009+u2009depleted HKB treatments revealed that the presence of HKB led to differential expression of 59% of genes, including the majority of genes involved in major carbon and nitrogen metabolic pathways. Genes coding for unidirectional enzymes for the utilization of glucose were upregulated in the presence of HKB, implying increased glycolytic activities during phagotrophy. Algae without HKB upregulated their expression of genes coding for ammonium transporters, implying uptake of inorganic nitrogen from the culture medium when prey were unavailable.ConclusionsTranscriptomic results agreed with previous observations that light had minimal effect on the population growth of Ochromonas sp. However, light led to the upregulation of a number of phototrophy- and phagotrophy-related genes, while the availability of bacterial prey led to prominent changes in major carbon and nitrogen metabolic pathways. Our study demonstrated the potential of transcriptomic approaches to improve our understanding of the trophic physiologies of complex mixotrophs, and revealed responses in Ochromonas sp. not apparent from traditional culture studies.

A tale of two mixotrophic chrysophytes: Insights into the metabolisms of two Ochromonas species (Chrysophyceae) through a comparison of gene expression

Ochromonas spp. strains CCMP1393 and BG-1 are phagotrophic phytoflagellates with different nutritional strategies. Strain CCMP1393 is an obligate phototroph while strain BG-1 readily grows in continuous darkness in the presence of bacterial prey. Growth and gene expression of strain CCMP1393 were investigated under conditions allowing phagotrophic, mixotrophic, or phototrophic nutrition. The availability of light and bacterial prey led to the differential expression of 42% or 45–59% of all genes, respectively. Data from strain CCMP1393 were compared to those from a study conducted previously on strain BG-1, and revealed notable differences in carbon and nitrogen metabolism between the 2 congeners under similar environmental conditions. Strain BG-1 utilized bacterial carbon and amino acids through glycolysis and the tricarboxylic acid cycle, while downregulating light harvesting and carbon fixation in the Calvin cycle when both light and bacteria were available. In contrast, the upregulation of genes related to photosynthesis, light harvesting, chlorophyll synthesis, and carbon fixation in the presence of light and prey for strain CCMP1393 implied that this species is more phototrophic than strain BG-1, and that phagotrophy may have enhanced phototrophy. Cellular chlorophyll a content was also significantly higher in strain CCMP1393 supplied with bacteria compared to those without prey. Our results thus point to very different physiological strategies for mixotrophic nutrition in these closely related chrysophyte species.