The voyage of the microbial eukaryote

MetaEuk-sensitive, high-throughput gene discovery, and annotation for large-scale eukaryotic metagenomics

BACKGROUND

Metagenomics is revolutionizing the study of microorganisms and their involvement in biological, biomedical, and geochemical processes, allowing us to investigate by direct sequencing a tremendous diversity of organisms without the need for prior cultivation. Unicellular eukaryotes play essential roles in most microbial communities as chief predators, decomposers, phototrophs, bacterial hosts, symbionts, and parasites to plants and animals. Investigating their roles is therefore of great interest to ecology, biotechnology, human health, and evolution. However, the generally lower sequencing coverage, their more complex gene and genome architectures, and a lack of eukaryote-specific experimental and computational procedures have kept them on the sidelines of metagenomics.

RESULTS

MetaEuk is a toolkit for high-throughput, reference-based discovery, and annotation of protein-coding genes in eukaryotic metagenomic contigs. It performs fast searches with 6-frame-translated fragments covering all possible exons and optimally combines matches into multi-exon proteins. We used a benchmark of seven diverse, annotated genomes to show that MetaEuk is highly sensitive even under conditions of low sequence similarity to the reference database. To demonstrate MetaEuk's power to discover novel eukaryotic proteins in large-scale metagenomic data, we assembled contigs from 912 samples of the Tara Oceans project. MetaEuk predicted >12,000,000 protein-coding genes in 8 days on ten 16-core servers. Most of the discovered proteins are highly diverged from known proteins and originate from very sparsely sampled eukaryotic supergroups.

CONCLUSION

The open-source (GPLv3) MetaEuk software (https://github.com/soedinglab/metaeuk) enables large-scale eukaryotic metagenomics through reference-based, sensitive taxonomic and functional annotation. Video abstract.

Heterogeneity of interactions of microbial communities in regions of Taihu Lake with different nutrient loadings: A network analysis

To investigate the differences in the interactions of microbial communities in two regions in Taihu Lake with different nutrient loadings [Meiliang Bay (MLB) and Xukou Bay (XKB)], water samples were collected and both intra- and inter-kingdom microbial community interactions were examined with network analysis. It is demonstrated that all of the bacterioplankton, microeukaryotes and inter-kingdom communities networks in Taihu Lake were non-random. For the networks of bacterioplankton and inter-kingdom community in XKB, higher clustering coefficient and average degree but lower average path length indexes were observed, indicating the nodes in XKB were more clustered and closely connected with plenty edges than those of MLB. The bacterioplankton and inter-kingdom networks were considerably larger and more complex with more module hubs and connectors in XKB compared with those of MLB, whereas the microeukaryotes networks were comparable and had no module hubs or connectors in the two lake zones. The phyla of Acidobacteria, Cyanobacteria and Planctomycetes maintained greater cooperation with other phyla in XKB, rather than competition. The relationships between microbial communities and environmental factors in MLB were weaker. Compared with the microbial community networks of XKB, less modules in networks of MLB were significantly correlated with total phosphorous and total nitrogen.

The heterogeneity of composition and assembly processes of the microbial community between different nutrient loading lake zones in Taihu Lake

To investigate the differences in the microbial community composition and assembly process in two lake zones (Meiliang Bay (MLB) and Xukou Bay (XKB) in Taihu Lake, China) with different nutrient loadings, water samples were collected. Both the 16S ribosomal RNA (rRNA) gene for the bacterial community and the 18S rRNA gene for the microeukaryote community were investigated using the Illumina second-generation sequencing platform (2 × 250 paired-end). The results indicated that both the bacterioplankton and microeukaryote community composition derived from the two lake zones were significantly different. Significantly higher operational taxonomic unit (OTU) richness (P < 0.01) and phylogenetic diversity (P < 0.05) were found for the bacterioplankton community of MLB. However, a comparable alpha diversity was found between the microeukaryote communities of MLB and XKB (P > 0.05). Environmental factors significantly affected the community compositions in XKB for both the bacterioplankton and microeukaryotes. However, they did not significantly influence the microbial community composition in MLB, except for a weak correlation between dissolved organic carbon (DOC) and the microeukaryote community. The microbial communities tended to be more phylogenetically clustered than expected by chance in the two lake zones. Moreover, the results of the phylogenetic structure suggest that deterministic processes played overwhelming roles in driving the assembly of both the bacterioplankton and microeukaryote community in XKB.

Evidence-based green algal genomics reveals marine diversity and ancestral characteristics of land plants

Background

Prasinophytes are widespread marine green algae that are related to plants. Cellular abundance of the prasinophyte Micromonas has reportedly increased in the Arctic due to climate-induced changes. Thus, studies of these unicellular eukaryotes are important for marine ecology and for understanding Viridiplantae evolution and diversification.

Results

We generated evidence-based Micromonas gene models using proteomics and RNA-Seq to improve prasinophyte genomic resources. First, sequences of four chromosomes in the 22 Mb Micromonas pusilla (CCMP1545) genome were finished. Comparison with the finished 21 Mb genome of Micromonas commoda (RCC299; named herein) shows they share ≤8,141 of ~10,000 protein-encoding genes, depending on the analysis method. Unlike RCC299 and other sequenced eukaryotes, CCMP1545 has two abundant repetitive intron types and a high percent (26 %) GC splice donors. Micromonas has more genus-specific protein families (19 %) than other genome sequenced prasinophytes (11 %). Comparative analyses using predicted proteomes from other prasinophytes reveal proteins likely related to scale formation and ancestral photosynthesis. Our studies also indicate that peptidoglycan (PG) biosynthesis enzymes have been lost in multiple independent events in select prasinophytes and plants. However, CCMP1545, polar Micromonas CCMP2099 and prasinophytes from other classes retain the entire PG pathway, like moss and glaucophyte algae. Surprisingly, multiple vascular plants also have the PG pathway, except the Penicillin-Binding Protein, and share a unique bi-domain protein potentially associated with the pathway. Alongside Micromonas experiments using antibiotics that halt bacterial PG biosynthesis, the findings highlight unrecognized phylogenetic complexity in PG-pathway retention and implicate a role in chloroplast structure or division in several extant Viridiplantae lineages.

Conclusions

Extensive differences in gene loss and architecture between related prasinophytes underscore their divergence. PG biosynthesis genes from the cyanobacterial endosymbiont that became the plastid, have been selectively retained in multiple plants and algae, implying a biological function. Our studies provide robust genomic resources for emerging model algae, advancing knowledge of marine phytoplankton and plant evolution.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-016-2585-6) contains supplementary material, which is available to authorized users.

Gene silencing in microalgae: mechanisms and biological roles

Microalgae exhibit enormous diversity and can potentially contribute to the production of biofuels and high value compounds. However, for most species, our knowledge of their physiology, metabolism, and gene regulation is fairly limited. In eukaryotes, gene silencing mechanisms play important roles in both the reversible repression of genes that are required only in certain contexts and the suppression of genome invaders such at transposons. The recent sequencing of several algal genomes is providing insights into the complexity of these mechanisms in microalgae. Collectively, glaucophyte, red, and green microalgae contain the machineries involved in repressive histone H3 lysine methylation, DNA cytosine methylation, and RNA interference. However, individual species often only have subsets of these gene silencing mechanisms. Moreover, current evidence suggests that algal silencing systems function in transposon and transgene repression but their role(s) in gene regulation or other cellular processes remains virtually unexplored, hindering rational genetic engineering efforts.

The Marine Microbial Eukaryote Transcriptome Sequencing Project (MMETSP): illuminating the functional diversity of eukaryotic life in the oceans through transcriptome sequencing

Current sampling of genomic sequence data from eukaryotes is relatively poor, biased, and inadequate to address important questions about their biology, evolution, and ecology; this Community Page describes a resource of 700 transcriptomes from marine microbial eukaryotes to help understand their role in the world's oceans.

RNA-mediated silencing in Algae: biological roles and tools for analysis of gene function

Algae are a large group of aquatic, typically photosynthetic, eukaryotes that include species from very diverse phylogenetic lineages, from those similar to land plants to those related to protist parasites. The recent sequencing of several algal genomes has provided insights into the great complexity of these organisms. Genomic information has also emphasized our lack of knowledge of the functions of many predicted genes, as well as the gene regulatory mechanisms in algae. Core components of the machinery for RNA-mediated silencing show widespread distribution among algal lineages, but they also seem to have been lost entirely from several species with relatively small nuclear genomes. Complex sets of endogenous small RNAs, including candidate microRNAs and small interfering RNAs, have now been identified by high-throughput sequencing in green, red, and brown algae. However, the natural roles of RNA-mediated silencing in algal biology remain poorly understood. Limited evidence suggests that small RNAs may function, in different algae, in defense mechanisms against transposon mobilization, in responses to nutrient deprivation and, possibly, in the regulation of recently evolved developmental processes. From a practical perspective, RNA interference (RNAi) is becoming a promising tool for assessing gene function by sequence-specific knockdown. Transient gene silencing, triggered with exogenously synthesized nucleic acids, and/or stable gene repression, involving genome-integrated transgenes, have been achieved in green algae, diatoms, yellow-green algae, and euglenoids. The development of RNAi technology in conjunction with system level "omics" approaches may provide the tools needed to advance our understanding of algal physiological and metabolic processes.

Discovery of nuclear-encoded genes for the neurotoxin saxitoxin in dinoflagellates

Saxitoxin is a potent neurotoxin that occurs in aquatic environments worldwide. Ingestion of vector species can lead to paralytic shellfish poisoning, a severe human illness that may lead to paralysis and death. In freshwaters, the toxin is produced by prokaryotic cyanobacteria; in marine waters, it is associated with eukaryotic dinoflagellates. However, several studies suggest that saxitoxin is not produced by dinoflagellates themselves, but by co-cultured bacteria. Here, we show that genes required for saxitoxin synthesis are encoded in the nuclear genomes of dinoflagellates. We sequenced >1.2×10⁶ mRNA transcripts from the two saxitoxin-producing dinoflagellate strains Alexandrium fundyense CCMP1719 and A. minutum CCMP113 using high-throughput sequencing technology. In addition, we used in silico transcriptome analyses, RACE, qPCR and conventional PCR coupled with Sanger sequencing. These approaches successfully identified genes required for saxitoxin-synthesis in the two transcriptomes. We focused on sxtA, the unique starting gene of saxitoxin synthesis, and show that the dinoflagellate transcripts of sxtA have the same domain structure as the cyanobacterial sxtA genes. But, in contrast to the bacterial homologs, the dinoflagellate transcripts are monocistronic, have a higher GC content, occur in multiple copies, contain typical dinoflagellate spliced-leader sequences and eukaryotic polyA-tails. Further, we investigated 28 saxitoxin-producing and non-producing dinoflagellate strains from six different genera for the presence of genomic sxtA homologs. Our results show very good agreement between the presence of sxtA and saxitoxin-synthesis, except in three strains of A. tamarense, for which we amplified sxtA, but did not detect the toxin. Our work opens for possibilities to develop molecular tools to detect saxitoxin-producing dinoflagellates in the environment.