Single cell genome analysis of an uncultured heterotrophic stramenopile

MarFERReT, an open-source, version-controlled reference library of marine microbial eukaryote functional genes

Metatranscriptomics generates large volumes of sequence data about transcribed genes in natural environments. Taxonomic annotation of these datasets depends on availability of curated reference sequences. For marine microbial eukaryotes, current reference libraries are limited by gaps in sequenced organism diversity and barriers to updating libraries with new sequence data, resulting in taxonomic annotation of about half of eukaryotic environmental transcripts. Here, we introduce Marine Functional EukaRyotic Reference Taxa (MarFERReT), a marine microbial eukaryotic sequence library designed for use with taxonomic annotation of eukaryotic metatranscriptomes. We gathered 902 publicly accessible marine eukaryote genomes and transcriptomes and assessed their sequence quality and cross-contamination issues, selecting 800 validated entries for inclusion in MarFERReT. Version 1.1 of MarFERReT contains reference sequences from 800 marine eukaryotic genomes and transcriptomes, covering 453 species- and strain-level taxa, totaling nearly 28 million protein sequences with associated NCBI and PR2 Taxonomy identifiers and Pfam functional annotations. The MarFERReT project repository hosts containerized build scripts, documentation on installation and use case examples, and information on new versions of MarFERReT.

A k-mer-Based Approach for Phylogenetic Classification of Taxa in Environmental Genomic Data

In the age of genome sequencing, whole-genome data is readily and frequently generated, leading to a wealth of new information that can be used to advance various fields of research. New approaches, such as alignment-free phylogenetic methods that utilize k-mer-based distance scoring, are becoming increasingly popular given their ability to rapidly generate phylogenetic information from whole-genome data. However, these methods have not yet been tested using environmental data, which often tends to be highly fragmented and incomplete. Here, we compare the results of one alignment-free approach (which utilizes the D2 statistic) to traditional multi-gene maximum likelihood trees in 3 algal groups that have high-quality genome data available. In addition, we simulate lower-quality, fragmented genome data using these algae to test method robustness to genome quality and completeness. Finally, we apply the alignment-free approach to environmental metagenome assembled genome data of unclassified Saccharibacteria and Trebouxiophyte algae, and single-cell amplified data from uncultured marine stramenopiles to demonstrate its utility with real datasets. We find that in all instances, the alignment-free method produces phylogenies that are comparable, and often more informative, than those created using the traditional multi-gene approach. The k-mer-based method performs well even when there are significant missing data that include marker genes traditionally used for tree reconstruction. Our results demonstrate the value of alignment-free approaches for classifying novel, often cryptic or rare, species, that may not be culturable or are difficult to access using single-cell methods, but fill important gaps in the tree of life.

Single-cell genotyping of phytoplankton from ocean water by gel-based cell manipulation

A comprehensive understanding of phytoplankton diversity is valuable for assessing an environment of interest as phytoplankton are primary producers in various aquatic food webs. Microscopic analyses are useful for diversity assessment based on characteristic cell morphologies. However, phylogenetic classification based solely on morphology requires an extremely high level of expertise. The genetic approach is another option for evaluating phytoplankton diversity; however, it cannot reveal morphological information. To integrate these two approaches, we developed an original technology that is referred to as microcavity array (MCA)/gel-based cell manipulation (GCM). The model experiments using monocultures of various phytoplankton indicated that the efficiencies of cell recovery and isolation of single-cell plankton were dependent on cell size and shape. Cells with widths larger than the cavity width showed high level of recovery and isolation efficiency. Subsequent whole-genome amplification (WGA) of isolated single-cell plankton provided a sufficient amount (approximately 30 μg) of WGA products for genetic analyses. Furthermore, we showed that MCA/GCM could directly analyze phytoplankton in ocean water obtained from Suruga Bay, Japan, without any cumbersome pretreatment. These results indicate that MCA/GCM technology is a powerful tool for elucidating the phytoplankton diversity in marine environment. This article is protected by copyright. All rights reserved.

Genome sequencing and de novo assembly of the giant unicellular alga Acetabularia acetabulum using droplet MDA

The macroscopic single-celled green alga Acetabularia acetabulum has been a model system in cell biology for more than a century. However, no genomic information is available from this species. Since the alga has a long life cycle, is difficult to grow in dense cultures, and has an estimated diploid genome size of almost 2 Gb, obtaining sufficient genomic material for genome sequencing is challenging. Here, we have attempted to overcome these challenges by amplifying genomic DNA using multiple displacement amplification (MDA) combined with microfluidics technology to distribute the amplification reactions across thousands of microscopic droplets. By amplifying and sequencing DNA from five single cells we were able to recover an estimated ~ 7–11% of the total genome, providing the first draft of the A. acetabulum genome. We highlight challenges associated with genome recovery and assembly of MDA data due to biases arising during genome amplification, and hope that our study can serve as a reference for future attempts on sequencing the genome from non-model eukaryotes.

Comparative genomics reveals new functional insights in uncultured MAST species

Heterotrophic lineages of stramenopiles exhibit enormous diversity in morphology, lifestyle, and habitat. Among them, the marine stramenopiles (MASTs) represent numerous independent lineages that are only known from environmental sequences retrieved from marine samples. The core energy metabolism characterizing these unicellular eukaryotes is poorly understood. Here, we used single-cell genomics to retrieve, annotate, and compare the genomes of 15 MAST species, obtained by coassembling sequences from 140 individual cells sampled from the marine surface plankton. Functional annotations from their gene repertoires are compatible with all of them being phagocytotic. The unique presence of rhodopsin genes in MAST species, together with their widespread expression in oceanic waters, supports the idea that MASTs may be capable of using sunlight to thrive in the photic ocean. Additional subsets of genes used in phagocytosis, such as proton pumps for vacuole acidification and peptidases for prey digestion, did not reveal particular trends in MAST genomes as compared with nonphagocytotic stramenopiles, except a larger presence and diversity of V-PPase genes. Our analysis reflects the complexity of phagocytosis machinery in microbial eukaryotes, which contrasts with the well-defined set of genes for photosynthesis. These new genomic data provide the essential framework to study ecophysiology of uncultured species and to gain better understanding of the function of rhodopsins and related carotenoids in stramenopiles.

A single-cell genomics pipeline for environmental microbial eukaryotes

Single-cell sequencing of environmental microorganisms is an essential component of the microbial ecology toolkit. However, large-scale targeted single-cell sequencing for the whole-genome recovery of uncultivated eukaryotes is lagging. The key challenges are low abundance in environmental communities, large complex genomes, and cell walls that are difficult to break. We describe a pipeline composed of state-of-the art single-cell genomics tools and protocols optimized for poorly studied and uncultivated eukaryotic microorganisms that are found at low abundance. This pipeline consists of seven distinct steps, beginning with sample collection and ending with genome annotation, each equipped with quality review steps to ensure high genome quality at low cost. We tested and evaluated each step on environmental samples and cultures of early-diverging lineages of fungi and Chromista/SAR. We show that genomes produced using this pipeline are almost as good as complete reference genomes for functional and comparative genomics for environmental microbial eukaryotes.

Single Cell Genomics Reveals Viruses Consumed by Marine Protists

The predominant model of the role of viruses in the marine trophic web is that of the “viral shunt,” where viral infection funnels a substantial fraction of the microbial primary and secondary production back to the pool of dissolved organic matter. Here, we analyzed the composition of non-eukaryotic DNA associated with individual cells of small, planktonic protists in the Gulf of Maine (GoM) and the Mediterranean Sea. We found viral DNA associated with a substantial fraction cells from the GoM (51%) and the Mediterranean Sea (35%). While Mediterranean SAGs contained a larger proportion of cells containing bacterial sequences (49%), a smaller fraction of cells contained bacterial sequences in the GoM (19%). In GoM cells, nearly identical bacteriophage and ssDNA virus sequences where found across diverse lineages of protists, suggesting many of these viruses are non-infective. The fraction of cells containing viral DNA varied among protistan lineages and reached 100% in Picozoa and Choanozoa. These two groups also contained significantly higher numbers of viral sequences than other identified taxa. We consider mechanisms that may explain the presence of viral DNA in protistan cells and conclude that protistan predation on free viral particles contributed to the observed patterns. These findings confirm prior experiments with protistan isolates and indicate that the viral shunt is complemented by a viral link in the marine microbial food web. This link may constitute a sink of viral particles in the ocean and has implications for the flow of carbon through the microbial food web.

SCGid: a consensus approach to contig filtering and genome prediction from single-cell sequencing libraries of uncultured eukaryotes

MOTIVATION

Whole-genome sequencing of uncultured eukaryotic genomes is complicated by difficulties in acquiring sufficient amounts of tissue. Single-cell genomics (SCG) by multiple displacement amplification provides a technical workaround, yielding whole-genome libraries which can be assembled de novo. Downsides of multiple displacement amplification include coverage biases and exacerbation of contamination. These factors affect assembly continuity and fidelity, complicating discrimination of genomes from contamination and noise by available tools. Uncultured eukaryotes and their relatives are often underrepresented in large sequence data repositories, further impairing identification and separation.

RESULTS

We compare the ability of filtering approaches to remove contamination and resolve eukaryotic draft genomes from SCG metagenomes, finding significant variation in outcomes. To address these inconsistencies, we introduce a consensus approach that is codified in the SCGid software package. SCGid parallelly filters assemblies using different approaches, yielding three intermediate drafts from which consensus is drawn. Using genuine and mock SCG metagenomes, we show that our approach corrects for variation among draft genomes predicted by individual approaches and outperforms them in recapitulating published drafts in a fast and repeatable way, providing a useful alternative to available methods and manual curation.

AVAILABILITY AND IMPLEMENTATION

The SCGid package is implemented in python and R. Source code is available at http://www.github.com/amsesk/SCGid under the GNU GPL 3.0 license.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

The need to account for cell biology in characterizing predatory mixotrophs in aquatic environments

Photosynthesis in eukaryotes first arose through phagocytotic processes wherein an engulfed cyanobacterium was not digested, but instead became a permanent organelle. Other photosynthetic lineages then arose when eukaryotic cells engulfed other already photosynthetic eukaryotic cells. Some of the resulting lineages subsequently lost their ability for phagocytosis, while many others maintained the ability to do both processes. These mixotrophic taxa have more complicated ecological roles, in that they are both primary producers and consumers that can shift more towards producing the organic matter that forms the base of aquatic food chains, or towards respiring and releasing CO2. We still have much to learn about which taxa are predatory mixotrophs as well as about the physiological consequences of this lifestyle, in part, because much of the diversity of unicellular eukaryotes in aquatic ecosystems remains uncultured. Here, we discuss existing methods for studying predatory mixotrophs, their individual biases, and how single-cell approaches can enhance knowledge of these important taxa. The question remains what the gold standard should be for assigning a mixotrophic status to ill-characterized or uncultured taxa-a status that dictates how organisms are incorporated into carbon cycle models and how their ecosystem roles may shift in future lakes and oceans. This article is part of a discussion meeting issue 'Single cell ecology'.

Unexpected mitochondrial genome diversity revealed by targeted single-cell genomics of heterotrophic flagellated protists

Most eukaryotic microbial diversity is uncultivated, under-studied and lacks nuclear genome data. Mitochondrial genome sampling is more comprehensive, but many phylogenetically important groups remain unsampled. Here, using a single-cell sorting approach combining tubulin-specific labelling with photopigment exclusion, we sorted flagellated heterotrophic unicellular eukaryotes from Pacific Ocean samples. We recovered 206 single amplified genomes, predominantly from underrepresented branches on the tree of life. Seventy single amplified genomes contained unique mitochondrial contigs, including 21 complete or near-complete mitochondrial genomes from formerly under-sampled phylogenetic branches, including telonemids, katablepharids, cercozoans and marine stramenopiles, effectively doubling the number of available samples of heterotrophic flagellate mitochondrial genomes. Collectively, these data identify a dynamic history of mitochondrial genome evolution including intron gain and loss, extensive patterns of genetic code variation and complex patterns of gene loss. Surprisingly, we found that stramenopile mitochondrial content is highly plastic, resembling patterns of variation previously observed only in plants.

Annual Protist Community Dynamics in a Freshwater Ecosystem Undergoing Contrasted Climatic Conditions: The Saint-Charles River (Canada)

Protists are key stone components of aquatic ecosystems, sustaining primary productivity and aquatic food webs. However, their diversity, ecology and structuring factors shaping their temporal distribution remain strongly misunderstood in freshwaters. Using high-throughput sequencing on water samples collected over 16 different months (including two summer and two winter periods), combined with geochemical measurements and climate monitoring, we comprehensively determined the pico- and nanoeukaryotic community composition and dynamics in a Canadian river undergoing prolonged ice-cover winters. Our analysis revealed a large protist diversity in this fluctuating ecosystem and clear seasonal patterns demonstrating a direct and/or indirect selective role of abiotic factors, such as water temperature or nitrogen concentrations, in structuring the eukaryotic microbial community. Nonetheless, our results also revealed that primary productivity, predatory as well as parasitism lifestyles, inferred from fine phylogenetic placements, remained potentially present over the annual cycle, despite the large seasonal fluctuations and the remodeling of the community composition under ice. In addition, potential interplays with the bacterial community composition were identified supporting a possible contribution of the bacterial community to the temporal dynamics of the protist community structure. Our results illustrate the complexity of the eukaryotic microbial community and provide a substantive and useful dataset to better understand the global freshwater ecosystem functioning.

The New Tree of Eukaryotes

For 15 years, the eukaryote Tree of Life (eToL) has been divided into five to eight major groupings, known as 'supergroups'. However, the tree has been profoundly rearranged during this time. The new eToL results from the widespread application of phylogenomics and numerous discoveries of major lineages of eukaryotes, mostly free-living heterotrophic protists. The evidence that supports the tree has transitioned from a synthesis of molecular phylogenetics and biological characters to purely molecular phylogenetics. Most current supergroups lack defining morphological or cell-biological characteristics, making the supergroup label even more arbitrary than before. Going forward, the combination of traditional culturing with maturing culture-free approaches and phylogenomics should accelerate the process of completing and resolving the eToL at its deepest levels.

A single-cell genome reveals diplonemid-like ancestry of kinetoplastid mitochondrial gene structure

Euglenozoa comprises euglenids, kinetoplastids, and diplonemids, with each group exhibiting different and highly unusual mitochondrial genome organizations. Although they are sister groups, kinetoplastids and diplonemids have very distinct mitochondrial genome architectures, requiring widespread insertion/deletion RNA editing and extensive trans-splicing, respectively, in order to generate functional transcripts. The evolutionary history by which these differing processes arose remains unclear. Using single-cell genomics, followed by small sub unit ribosomal DNA and multigene phylogenies, we identified an isolated marine cell that branches on phylogenetic trees as a sister to known kinetoplastids. Analysis of single-cell amplified genomic material identified multiple mitochondrial genome contigs. These revealed a gene architecture resembling that of diplonemid mitochondria, with small fragments of genes encoded out of order and or on different contigs, indicating that these genes require extensive trans-splicing. Conversely, no requirement for kinetoplastid-like insertion/deletion RNA-editing was detected. Additionally, while we identified some proteins so far only found in kinetoplastids, we could not unequivocally identify mitochondrial RNA editing proteins. These data invite the hypothesis that extensive genome fragmentation and trans-splicing were the ancestral states for the kinetoplastid-diplonemid clade but were lost during the kinetoplastid radiation. This study demonstrates that single-cell approaches can successfully retrieve lineages that represent important new branches on the tree of life, and thus can illuminate major evolutionary and functional transitions in eukaryotes.

This article is part of a discussion meeting issue ‘Single cell ecology’.

A single-cell genome perspective on studying intracellular associations in unicellular eukaryotes

Single-cell genomics (SCG) methods provide a unique opportunity to analyse whole genome information at the resolution of an individual cell. While SCG has been extensively used to investigate bacterial and archaeal genomes, the technique has been rarely used to access the genetic makeup of uncultivated microbial eukaryotes. In this regard, the use of SCG can provide a wealth of information; not only do the methods allow exploration of the genome, they can also help elucidate the relationship between the cell and intracellular entities extant in nearly all eukaryotes. SCG enables the study of total eukaryotic cellular DNA, which in turn allows us to better understand the evolutionary history and diversity of life, and the physiological interactions that define complex organisms.

This article is part of a discussion meeting issue ‘Single cell ecology’.

Single cell ecology

Cells are the building blocks of life, from single-celled microbes through to multi-cellular organisms. To understand a multitude of biological processes we need to understand how cells behave, how they interact with each other and how they respond to their environment. The use of new methodologies is changing the way we study cells allowing us to study them on minute scales and in unprecedented detail. These same methods are allowing researchers to begin to sample the vast diversity of microbes that dominate natural environments. The aim of this special issue is to bring together research and perspectives on the application of new approaches to understand the biological properties of cells, including how they interact with other biological entities. This article is part of a discussion meeting issue 'Single cell ecology'.

Combining morphology, behaviour and genomics to understand the evolution and ecology of microbial eukaryotes

Microbial eukaryotes (protists) are structurally, developmentally and behaviourally more complex than their prokaryotic cousins. This complexity makes it more difficult to translate genomic and metagenomic data into accurate functional inferences about systems ranging all the way from molecular and cellular levels to global ecological networks. This problem can be traced back to the advent of the cytoskeleton and endomembrane systems at the origin of eukaryotes, which endowed them with a range of complex structures and behaviours that still largely dominate how they evolve and interact within microbial communities. But unlike the diverse metabolic properties that evolved within prokaryotes, the structural and behavioural characteristics that strongly define how protists function in the environment cannot readily be inferred from genomic data, since there is generally no simple correlation between a gene and a discrete activity or function. A deeper understanding of protists at both cellular and ecological levels, therefore, requires not only high-throughput genomics but also linking such data to direct observations of natural history and cell biology. This is challenging since these observations typically require cultivation, which is lacking for most protists. Potential remedies with current technology include developing a more phylogenetically diverse range of model systems to better represent the diversity, as well as combining high-throughput, single-cell genomics with microscopic documentation of the subject cells to link sequence with structure and behaviour. This article is part of a discussion meeting issue 'Single cell ecology'.

Leveraging single-cell genomics to expand the fungal tree of life

Environmental DNA surveys reveal that most fungal diversity represents uncultured species. We sequenced the genomes of eight uncultured species across the fungal tree of life using a new single-cell genomics pipeline. We show that, despite a large variation in genome and gene space recovery from each single amplified genome (SAG), ≥90% can be recovered by combining multiple SAGs. SAGs provide robust placement for early-diverging lineages and infer a diploid ancestor of fungi. Early-diverging fungi share metabolic deficiencies and show unique gene expansions correlated with parasitism and unculturability. Single-cell genomics holds great promise in exploring fungal diversity, life cycles and metabolic potential.

Multigene phylogeny and cell evolution of chromist infrakingdom Rhizaria: contrasting cell organisation of sister phyla Cercozoa and Retaria

Infrakingdom Rhizaria is one of four major subgroups with distinct cell body plans that comprise eukaryotic kingdom Chromista. Unlike other chromists, Rhizaria are mostly heterotrophic flagellates, amoebae or amoeboflagellates, commonly with reticulose (net-like) or filose (thread-like) feeding pseudopodia; uniquely for eukaryotes, cilia have proximal ciliary transition-zone hub-lattices. They comprise predominantly flagellate phylum Cercozoa and reticulopodial phylum Retaria, whose exact phylogenetic relationship has been uncertain. Given even less clear relationships amongst cercozoan classes, we sequenced partial transcriptomes of seven Cercozoa representing five classes and endomyxan retarian Filoreta marina to establish 187-gene multiprotein phylogenies. Ectoreta (retarian infraphyla Foraminifera, Radiozoa) branch within classical Cercozoa as sister to reticulose Endomyxa. This supports recent transfer of subphylum Endomyxa from Cercozoa to Retaria alongside subphylum Ectoreta which embraces classical retarians where capsules or tests subdivide cells into organelle-containing endoplasm and anastomosing pseudopodial net-like ectoplasm. Cercozoa are more homogeneously filose, often with filose pseudopodia and/or posterior ciliary gliding motility: zooflagellate Helkesimastix and amoeboid Guttulinopsis form a strongly supported clade, order Helkesida. Cercomonads are polyphyletic (Cercomonadida sister to glissomonads; Paracercomonadida deeper). Thecofilosea are a clade, whereas Imbricatea may not be; Sarcomonadea may be paraphyletic. Helkesea and Metromonadea are successively deeper outgroups within cercozoan subphylum Monadofilosa; subphylum Reticulofilosa (paraphyletic on site-heterogeneous trees) branches earliest, Granofilosea before Chlorarachnea. Our multiprotein trees confirm that Rhizaria are sisters of infrakingdom Halvaria (Alveolata, Heterokonta) within chromist subkingdom Harosa (= SAR); they further support holophyly of chromist subkingdom Hacrobia, and are consistent with holophyly of Chromista as sister of kingdom Plantae. Site-heterogeneous rDNA trees group Kraken with environmental DNA clade 'eSarcomonad', not Paracercomonadida. Ectoretan fossil dates evidence ultrarapid episodic stem sequence evolution. We discuss early rhizarian cell evolution and multigene tree coevolutionary patterns, gene-paralogue evidence for chromist monophyly, and integrate this with fossil evidence for the age of Rhizaria and eukaryote cells, and revise rhizarian classification.

Genome-reconstruction for eukaryotes from complex natural microbial communities

Microbial eukaryotes are integral components of natural microbial communities, and their inclusion is critical for many ecosystem studies, yet the majority of published metagenome analyses ignore eukaryotes. In order to include eukaryotes in environmental studies, we propose a method to recover eukaryotic genomes from complex metagenomic samples. A key step for genome recovery is separation of eukaryotic and prokaryotic fragments. We developed a k-mer-based strategy, EukRep, for eukaryotic sequence identification and applied it to environmental samples to show that it enables genome recovery, genome completeness evaluation, and prediction of metabolic potential. We used this approach to test the effect of addition of organic carbon on a geyser-associated microbial community and detected a substantial change of the community metabolism, with selection against almost all candidate phyla bacteria and archaea and for eukaryotes. Near complete genomes were reconstructed for three fungi placed within the Eurotiomycetes and an arthropod. While carbon fixation and sulfur oxidation were important functions in the geyser community prior to carbon addition, the organic carbon-impacted community showed enrichment for secreted proteases, secreted lipases, cellulose targeting CAZymes, and methanol oxidation. We demonstrate the broader utility of EukRep by reconstructing and evaluating relatively high-quality fungal, protist, and rotifer genomes from complex environmental samples. This approach opens the way for cultivation-independent analyses of whole microbial communities.

Single-cell genomics of multiple uncultured stramenopiles reveals underestimated functional diversity across oceans

Single-celled eukaryotes (protists) are critical players in global biogeochemical cycling of nutrients and energy in the oceans. While their roles as primary producers and grazers are well appreciated, other aspects of their life histories remain obscure due to challenges in culturing and sequencing their natural diversity. Here, we exploit single-cell genomics and metagenomics data from the circumglobal Tara Oceans expedition to analyze the genome content and apparent oceanic distribution of seven prevalent lineages of uncultured heterotrophic stramenopiles. Based on the available data, each sequenced genome or genotype appears to have a specific oceanic distribution, principally correlated with water temperature and depth. The genome content provides hypotheses for specialization in terms of cell motility, food spectra, and trophic stages, including the potential impact on their lifestyles of horizontal gene transfer from prokaryotes. Our results support the idea that prominent heterotrophic marine protists perform diverse functions in ocean ecology.

The biology of many marine protists, such as stramenopiles, remains obscure. Here, the authors exploit single-cell genomics and metagenomics to analyze the genome content and apparent oceanic distribution of seven prevalent lineages of uncultured heterotrophic stramenopiles.

Exploration of the Germline Genome of the Ciliate Chilodonella uncinata through Single-Cell Omics (Transcriptomics and Genomics)

Separate germline and somatic genomes are found in numerous lineages across the eukaryotic tree of life, often separated into distinct tissues (e.g., in plants, animals, and fungi) or distinct nuclei sharing a common cytoplasm (e.g., in ciliates and some foraminifera). In ciliates, germline-limited (i.e., micronuclear-specific) DNA is eliminated during the development of a new somatic (i.e., macronuclear) genome in a process that is tightly linked to large-scale genome rearrangements, such as deletions and reordering of protein-coding sequences. Most studies of germline genome architecture in ciliates have focused on the model ciliates Oxytricha trifallax, Paramecium tetraurelia, and Tetrahymena thermophila, for which the complete germline genome sequences are known. Outside of these model taxa, only a few dozen germline loci have been characterized from a limited number of cultivable species, which is likely due to difficulties in obtaining sufficient quantities of “purified” germline DNA in these taxa. Combining single-cell transcriptomics and genomics, we have overcome these limitations and provide the first insights into the structure of the germline genome of the ciliate Chilodonella uncinata, a member of the understudied class Phyllopharyngea. Our analyses reveal the following: (i) large gene families contain a disproportionate number of genes from scrambled germline loci; (ii) germline-soma boundaries in the germline genome are demarcated by substantial shifts in GC content; (iii) single-cell omics techniques provide large-scale quality germline genome data with limited effort, at least for ciliates with extensively fragmented somatic genomes. Our approach provides an efficient means to understand better the evolution of genome rearrangements between germline and soma in ciliates.

Our understanding of the distinctions between germline and somatic genomes in ciliates has largely relied on studies of a few model genera (e.g., Oxytricha, Paramecium, Tetrahymena). We have used single-cell omics to explore germline-soma distinctions in the ciliate Chilodonella uncinata, which likely diverged from the better-studied ciliates ~700 million years ago. The analyses presented here indicate that developmentally regulated genome rearrangements between germline and soma are demarcated by rapid transitions in local GC composition and lead to diversification of protein families. The approaches used here provide the basis for future work aimed at discerning the evolutionary impacts of germline-soma distinctions among diverse ciliates.

Evaluation of single-cell genomics to address evolutionary questions using three SAGs of the choanoflagellate Monosiga brevicollis

Single-cell genomics (SCG) appeared as a powerful technique to get genomic information from uncultured organisms. However, SCG techniques suffer from biases at the whole genome amplification step that can lead to extremely variable numbers of genome recovery (5-100%). Thus, it is unclear how useful can SCG be to address evolutionary questions on uncultured microbial eukaryotes. To provide some insights into this, we here analysed 3 single-cell amplified genomes (SAGs) of the choanoflagellate Monosiga brevicollis, whose genome is known. Our results show that each SAG has a different, independent bias, yielding different levels of genome recovery for each cell (6-36%). Genes often appear fragmented and are split into more genes during annotation. Thus, analyses of gene gain and losses, gene architectures, synteny and other genomic features can not be addressed with a single SAG. However, the recovery of phylogenetically-informative protein domains can be up to 55%. This means SAG data can be used to perform accurate phylogenomic analyses. Finally, we also confirm that the co-assembly of several SAGs improves the general genomic recovery. Overall, our data show that, besides important current limitations, SAGs can still provide interesting and novel insights from poorly-known, uncultured organisms.

Improved genome recovery and integrated cell-size analyses of individual uncultured microbial cells and viral particles

Microbial single-cell genomics can be used to provide insights into the metabolic potential, interactions, and evolution of uncultured microorganisms. Here we present WGA-X, a method based on multiple displacement amplification of DNA that utilizes a thermostable mutant of the phi29 polymerase. WGA-X enhances genome recovery from individual microbial cells and viral particles while maintaining ease of use and scalability. The greatest improvements are observed when amplifying high G+C content templates, such as those belonging to the predominant bacteria in agricultural soils. By integrating WGA-X with calibrated index-cell sorting and high-throughput genomic sequencing, we are able to analyze genomic sequences and cell sizes of hundreds of individual, uncultured bacteria, archaea, protists, and viral particles, obtained directly from marine and soil samples, in a single experiment. This approach may find diverse applications in microbiology and in biomedical and forensic studies of humans and other multicellular organisms.

Single-cell genomics can be used to study uncultured microorganisms. Here, Stepanauskas et al. present a method combining improved multiple displacement amplification and FACS, to obtain genomic sequences and cell size information from uncultivated microbial cells and viral particles in environmental samples.

Recent progress in diatom genomics and epigenomics

Diatoms are one of the most diverse and successful groups of phytoplankton at the base of the food chain, sustaining life in the ocean and performing vital biogeochemical functions. The last fifteen years have witnessed the comprehensive analysis of several diatom genomes, revealing that they bear traces of their endosymbiotic origins from algal and heterotrophic ancestors, as well as significant gene transfer from bacteria. Their chimeric genomes are further regulated by a range of chromatin-based processes that are characteristic of both plant and animal genomes. We discuss the conservation of gene regulatory mechanisms in diatoms and propose that epigenetic processes may have a significant role in mediating responses to a highly dynamic and unpredictable environment in these organisms.

Accessing the genomic information of unculturable oceanic picoeukaryotes by combining multiple single cells

Pico-sized eukaryotes play key roles in the functioning of marine ecosystems, but we still have a limited knowledge on their ecology and evolution. The MAST-4 lineage is of particular interest, since it is widespread in surface oceans, presents ecotypic differentiation and has defied culturing efforts so far. Single cell genomics (SCG) are promising tools to retrieve genomic information from these uncultured organisms. However, SCG are based on whole genome amplification, which normally introduces amplification biases that limit the amount of genomic data retrieved from a single cell. Here, we increase the recovery of genomic information from two MAST-4 lineages by co-assembling short reads from multiple Single Amplified Genomes (SAGs) belonging to evolutionary closely related cells. We found that complementary genomic information is retrieved from different SAGs, generating co-assembly that features >74% of genome recovery, against about 20% when assembled individually. Even though this approach is not aimed at generating high-quality draft genomes, it allows accessing to the genomic information of microbes that would otherwise remain unreachable. Since most of the picoeukaryotes still remain uncultured, our work serves as a proof-of-concept that can be applied to other taxa in order to extract genomic data and address new ecological and evolutionary questions.

Survey of the green picoalga Bathycoccus genomes in the global ocean

Bathycoccus is a cosmopolitan green micro-alga belonging to the Mamiellophyceae, a class of picophytoplankton that contains important contributors to oceanic primary production. A single species of Bathycoccus has been described while the existence of two ecotypes has been proposed based on metagenomic data. A genome is available for one strain corresponding to the described phenotype. We report a second genome assembly obtained by a single cell genomics approach corresponding to the second ecotype. The two Bathycoccus genomes are divergent enough to be unambiguously distinguishable in whole DNA metagenomic data although they possess identical sequence of the 18S rRNA gene including in the V9 region. Analysis of 122 global ocean whole DNA metagenome samples from the Tara-Oceans expedition reveals that populations of Bathycoccus that were previously identified by 18S rRNA V9 metabarcodes are only composed of these two genomes. Bathycoccus is relatively abundant and widely distributed in nutrient rich waters. The two genomes rarely co-occur and occupy distinct oceanic niches in particular with respect to depth. Metatranscriptomic data provide evidence for gain or loss of highly expressed genes in some samples, suggesting that the gene repertoire is modulated by environmental conditions.

A Phylogenomic Framework to Study the Diversity and Evolution of Stramenopiles (=Heterokonts)

Stramenopiles or heterokonts constitute one of the most speciose and diverse clades of protists. It includes ecologically important algae (such as diatoms or large multicellular brown seaweeds), as well as heterotrophic (e.g., bicosoecids, MAST groups) and parasitic (e.g., Blastocystis, oomycetes) species. Despite their evolutionary and ecological relevance, deep phylogenetic relationships among stramenopile groups, inferred mostly from small-subunit rDNA phylogenies, remain unresolved, especially for the heterotrophic taxa. Taking advantage of recently released stramenopile transcriptome and genome sequences, as well as data from the genomic assembly of the MAST-3 species Incisomonas marina generated in our laboratory, we have carried out the first extensive phylogenomic analysis of stramenopiles, including representatives of most major lineages. Our analyses, based on a large data set of 339 widely distributed proteins, strongly support a root of stramenopiles lying between two clades, Bigyra and Gyrista (Pseudofungi plus Ochrophyta). Additionally, our analyses challenge the Phaeista-Khakista dichotomy of photosynthetic stramenopiles (ochrophytes) as two groups previously considered to be part of the Phaeista (Pelagophyceae and Dictyochophyceae), branch with strong support with the Khakista (Bolidophyceae and Diatomeae). We propose a new classification of ochrophytes within the two groups Chrysista and Diatomista to reflect the new phylogenomic results. Our stramenopile phylogeny provides a robust phylogenetic framework to investigate the evolution and diversification of this group of ecologically relevant protists.

Laterally Transferred Gene Recruited as a Venom in Parasitoid Wasps

Parasitoid wasps use venom to manipulate the immunity and metabolism of their host insects in a variety of ways to provide resources for their offspring. Yet, how genes are recruited and evolve to perform venom functions remain open questions. A recently recognized source of eukaryotic genome innovation is lateral gene transfer (LGT). Glycoside hydrolase family 19 (GH19) chitinases are widespread in bacteria, microsporidia, and plants where they are used in nutrient acquisition or defense, but have previously not been known in metazoans. In this study, a GH19 chitinase LGT is described from the unicellular microsporidia/Rozella clade into parasitoid wasps of the superfamily Chalcidoidea, where it has become recruited as a venom protein. The GH19 chitinase is present in 15 species of chalcidoid wasps representing four families, and phylogenetic analysis indicates that it was laterally transferred near or before the origin of Chalcidoidea (∼95 Ma). The GH19 chitinase gene is highly expressed in the venom gland of at least seven species, indicating a role in the complex host manipulations performed by parasitoid wasp venom. RNAi knockdown in the model parasitoid Nasonia vitripennis reveals that-following envenomation-the GH19 chitinase induces fly hosts to upregulate genes involved in an immune response to fungi. A second, independent LGT of GH19 chitinase from microsporidia into mosquitoes was also found, also supported by phylogenetic reconstructions. Besides these two LGT events, GH19 chitinase is not found in any other sequenced animal genome, or in any fungi outside the microsporidia/Rozella clade.

Ecological and evolutionary significance of novel protist lineages

Environmental molecular surveys targeting protist diversity have unveiled novel and uncultured lineages in a variety of ecosystems, ranging from completely new high-rank lineages, to new taxa moderately related to previously described organisms. The ecological roles of some of these novel taxa have been studied, showing that in certain habitats they may be responsible for critical environmental processes. Moreover, from an evolutionary perspective they still need to be included in a more accurate and wider understanding of the eukaryotic tree of life. These seminal discoveries promoted the development and use of a wide range of more in-depth culture-independent approaches to access this diversity, from metabarcoding and metagenomics to single cell genomics and FISH. Nonetheless, culturing using classical or innovative approaches is also essential to better characterize this new diversity. Ecologists and evolutionary biologists now face the challenge of apprehending the significance of this new diversity within the eukaryotic tree of life.

Capturing prokaryotic dark matter genomes

Prokaryotes are the most diverse and abundant cellular life forms on Earth. Most of them, identified by indirect molecular approaches, belong to microbial dark matter. The advent of metagenomic and single-cell genomic approaches has highlighted the metabolic capabilities of numerous members of this dark matter through genome reconstruction. Thus, linking functions back to the species has revolutionized our understanding of how ecosystem function is sustained by the microbial world. This review will present discoveries acquired through the illumination of prokaryotic dark matter genomes by these innovative approaches.

Emerging concepts promising new horizons for marine biodiscovery and synthetic biology

The vast oceans of the world, which comprise a huge variety of unique ecosystems, are emerging as a rich and relatively untapped source of novel bioactive compounds with invaluable biotechnological and pharmaceutical potential. Evidence accumulated over the last decade has revealed that the diversity of marine microorganisms is enormous with many thousands of bacterial species detected that were previously unknown. Associated with this diversity is the production of diverse repertoires of bioactive compounds ranging from peptides and enzymes to more complex secondary metabolites that have significant bioactivity and thus the potential to be exploited for innovative biotechnology. Here we review the discovery and functional potential of marine bioactive peptides such as lantibiotics, nanoantibiotics and peptidomimetics, which have received particular attention in recent years in light of their broad spectrum of bioactivity. The significance of marine peptides in cell-to-cell communication and how this may be exploited in the discovery of novel bioactivity is also explored. Finally, with the recent advances in bioinformatics and synthetic biology, it is becoming clear that the integration of these disciplines with genetic and biochemical characterization of the novel marine peptides, offers the most potential in the development of the next generation of societal solutions.

Single-cell sequencing technologies: current and future

Intensively developed in the last few years, single-cell sequencing technologies now present numerous advantages over traditional sequencing methods for solving the problems of biological heterogeneity and low quantities of available biological materials. The application of single-cell sequencing technologies has profoundly changed our understanding of a series of biological phenomena, including gene transcription, embryo development, and carcinogenesis. However, before single-cell sequencing technologies can be used extensively, researchers face the serious challenge of overcoming inherent issues of high amplification bias, low accuracy and reproducibility. Here, we simply summarize the techniques used for single-cell isolation, and review the current technologies used in single-cell genomic, transcriptomic, and epigenomic sequencing. We discuss the merits, defects, and scope of application of single-cell sequencing technologies and then speculate on the direction of future developments.