Foraminifera as a model of eukaryotic genome dynamism

In contrast to the canonical view that genomes cycle only between haploid and diploid states, many eukaryotes have dynamic genomes that change content throughout an individual's life cycle. However, the few detailed studies of microeukaryotic life cycles render our understanding of eukaryotic genome dynamism incomplete. Foraminifera (Rhizaria) are an ecologically important, yet understudied, clade of microbial eukaryotes with complex life cycles that include changes in ploidy and genome organization. Here, we apply fluorescence microscopy and image analysis techniques to over 2,800 nuclei in 110 cells to characterize the life cycle of Allogromia laticollaris strain Cold Spring Harbor (CSH), one of few cultivable foraminifera species. We show that haploidy and diploidy are brief moments in the A. laticollaris life cycle and that A. laticollaris nuclei endoreplicate up to 12,000 times the haploid genome size. We find that A. laticollaris reorganizes a highly endoreplicated nucleus into thousands of haploid genomes through a non-canonical mechanism called Zerfall, in which the nuclear envelope degrades and extrudes chromatin into the cytoplasm. Based on these findings, along with changes in nuclear architecture across the life cycle, we believe that A. laticollaris uses spatio-temporal mechanisms to delineate germline and somatic DNA within a single nucleus. The analyses here extend our understanding of the genome dynamics across the eukaryotic tree of life.IMPORTANCEIn traditional depictions of eukaryotes (i.e., cells with nuclei), life cycles alternate only between haploid and diploid phases, overlooking studies of diverse microeukaryotic lineages (e.g., amoebae, ciliates, and flagellates) that show dramatic variation in DNA content throughout their life cycles. Endoreplication of genomes enables cells to grow to large sizes and perhaps to also respond to changes in their environments. Few microeukaryotic life cycles have been studied in detail, which limits our understanding of how eukaryotes regulate and transmit their DNA across generations. Here, we use microscopy to study the life cycle of Allogromia laticollaris strain CSH, an early-diverging lineage within the Foraminifera (an ancient clade of predominantly marine amoebae). We show that DNA content changes significantly throughout their life cycle and further describe an unusual process called Zerfall, by which this species reorganizes a large nucleus with up to 12,000 genome copies into hundreds of small gametic nuclei, each with a single haploid genome. Our results are consistent with the idea that all eukaryotes demarcate germline DNA to pass on to offspring amidst more flexible somatic DNA and extend the known diversity of eukaryotic life cycles.

Somatic genome architecture and molecular evolution are decoupled in "young" linage-specific gene families in ciliates

The evolution of lineage-specific gene families remains poorly studied across the eukaryotic tree of life, with most analyses focusing on the recent evolution of de novo genes in model species. Here we explore the origins of lineage-specific genes in ciliates, a ~1 billion year old clade of microeukaryotes that are defined by their division of somatic and germline functions into distinct nuclei. Previous analyses on conserved gene families have shown the effect of ciliates' unusual genome architecture on gene family evolution: extensive genome processing-the generation of thousands of gene-sized somatic chromosomes from canonical germline chromosomes-is associated with larger and more diverse gene families. To further study the relationship between ciliate genome architecture and gene family evolution, we analyzed lineage specific gene families from a set of 46 transcriptomes and 12 genomes representing x species from eight ciliate classes. We assess how the evolution lineage-specific gene families occurs among four groups of ciliates: extensive fragmenters with gene-size somatic chromosomes, non-extensive fragmenters with "large'' multi-gene somatic chromosomes, Heterotrichea with highly polyploid somatic genomes and Karyorelictea with 'paradiploid' somatic genomes. Our analyses demonstrate that: 1) most lineage-specific gene families are found at shallow taxonomic scales; 2) extensive genome processing (i.e., gene unscrambling) during development likely influences the size and number of young lineage-specific gene families; and 3) the influence of somatic genome architecture on molecular evolution is increasingly apparent in older gene families. Altogether, these data highlight the influences of genome architecture on the evolution of lineage-specific gene families in eukaryotes.

Long-Read-Based Genome Assembly Reveals Numerous Endogenous Viral Elements in the Green Algal Bacterivore Cymbomonas tetramitiformis

The marine tetraflagellate Cymbomonas tetramitiformis has drawn attention as an early diverging green alga that uses a phago-mixotrophic mode of nutrition (i.e., the ability to derive nourishment from both photosynthesis and bacterial prey). The Cymbomonas nuclear genome was sequenced previously, but due to the exclusive use of short-read (Illumina) data, the assembly suffered from missing a large proportion of the genome's repeat regions. For this study, we generated Oxford Nanopore long-read and additional short-read Illumina data and performed a hybrid assembly that significantly improved the total assembly size and contiguity. Numerous endogenous viral elements were identified in the repeat regions of the new assembly. These include the complete genome of a giant Algavirales virus along with many genomes of integrated Polinton-like viruses (PLVs) from two groups: Gezel-like PLVs and a novel group of prasinophyte-specific PLVs. The integrated ∼400 kb genome of the giant Algavirales virus is the first account of the association of the uncultured viral family AG_03 with green algae. The complete PLV genomes from C. tetramitiformis ranged between 15 and 25 kb in length and showed a diverse gene content. In addition, heliorhodopsin gene-containing repeat elements of putative mirusvirus origin were identified. These results illustrate past (and possibly ongoing) multiple alga-virus interactions that accompanied the genome evolution of C. tetramitiformis.

A small RNA-guided PRC2 complex eliminates DNA as an extreme form of transposon silencing

In animal germlines, transposons are silenced at the transcriptional or post-transcriptional level to prevent deleterious expression. Ciliates employ a more direct approach by physically eliminating transposons from their soma, utilizing piRNAs to recognize transposons and imprecisely excise them. Ancient, mutated transposons often do not require piRNAs and are precisely eliminated. Here, we characterize the Polycomb Repressive Complex 2 (PRC2) in Paramecium and demonstrate its involvement in the removal of transposons and transposon-derived DNA. Our results reveal a striking difference between the elimination of new and ancient transposons at the chromatin level and show that the complex may be guided by Piwi-bound small RNAs (sRNAs). We propose that imprecise elimination in ciliates originates from an ancient transposon silencing mechanism, much like in plants and metazoans, through sRNAs, repressive methylation marks, and heterochromatin formation. However, it is taken a step further by eliminating DNA as an extreme form of transposon silencing.

Identification of novel, functional, long noncoding RNAs involved in programmed, large-scale genome rearrangements

Noncoding RNAs (ncRNAs) make up to ∼98% percent of the transcriptome of a given organism. In recent years, one relatively new class of ncRNAs, long noncoding RNAs (lncRNAs), were shown to be more than mere by-products of gene expression and regulation. The unicellular eukaryote Paramecium tetraurelia is a member of the ciliate phylum, an extremely heterogeneous group of organisms found in most bodies of water across the globe. A hallmark of ciliate genetics is nuclear dimorphism and programmed elimination of transposons and transposon-derived DNA elements, the latter of which is essential for the maintenance of the somatic genome. Paramecium and ciliates in general harbor a plethora of different ncRNA species, some of which drive the process of large-scale genome rearrangements, including DNA elimination, during sexual development. Here, we identify and validate the first known functional lncRNAs in ciliates to date. Using deep-sequencing and subsequent bioinformatic processing and experimental validation, we show that Paramecium expresses at least 15 lncRNAs. These candidates were predicted by a highly conservative pipeline, and informatic analyses hint at differential expression during development. Depletion of two lncRNAs, lnc1 and lnc15, resulted in clear phenotypes, decreased survival, morphological impairment, and a global effect on DNA elimination.

Old genes in new places: A taxon-rich analysis of interdomain lateral gene transfer events

Vertical inheritance is foundational to Darwinian evolution, but fails to explain major innovations such as the rapid spread of antibiotic resistance among bacteria and the origin of photosynthesis in eukaryotes. While lateral gene transfer (LGT) is recognized as an evolutionary force in prokaryotes, the role of LGT in eukaryotic evolution is less clear. With the exception of the transfer of genes from organelles to the nucleus, a process termed endosymbiotic gene transfer (EGT), the extent of interdomain transfer from prokaryotes to eukaryotes is highly debated. A common critique of studies of interdomain LGT is the reliance on the topology of single-gene trees that attempt to estimate more than one billion years of evolution. We take a more conservative approach by identifying cases in which a single clade of eukaryotes is found in an otherwise prokaryotic gene tree (i.e. exclusive presence). Starting with a taxon-rich dataset of over 13,600 gene families and passing data through several rounds of curation, we identify and categorize the function of 306 interdomain LGT events into diverse eukaryotes, including 189 putative EGTs, 52 LGTs into Opisthokonta (i.e. animals, fungi and their microbial relatives), and 42 LGTs nearly exclusive to anaerobic eukaryotes. To assess differential gene loss as an explanation for exclusive presence, we compare branch lengths within each LGT tree to a set of vertically-inherited genes subsampled to mimic gene loss (i.e. with the same taxonomic sampling) and consistently find shorter relative distance between eukaryotes and prokaryotes in LGT trees, a pattern inconsistent with gene loss. Our methods provide a framework for future studies of interdomain LGT and move the field closer to an understanding of how best to model the evolutionary history of eukaryotes.

How Ciliated Protists Survive by Cysts: Some Key Points During Encystment and Excystment

Forming cysts is a common and important bionomic strategy for microorganisms to persist in harsh environments. In ciliated protists, many species have been reported to form cysts when facing unfavorable conditions. Despite traditional studies on the morphological features of cysts and the chemical composition of cyst wall, recent research has focused more on the molecular mechanisms of encystment. The present work reviews studies on developmental features and molecular information of resting cysts in ciliates, and pays more attention to the following questions: what are the inducing factors of encystment and excystment? How does the cell change morphologically during these dynamic processes? And what molecular mechanisms underlie those changes? We also present and summarize the characteristics of cysts from diverse ciliate lineages in a phylogenetic framework, aiming to provide new perspectives for studies on adaptive evolution of unicellular eukaryotes.

Deciphering phylogenetic relationships in class Karyorelictea (Protista, Ciliophora) based on updated multi-gene information with establishment of a new order Wilbertomorphida n. ord

The class Karyorelictea, a unique assemblage of ciliates, is a key group in deciphering ciliate evolution history. However, the systematic relationships among members of this class remain poorly understood. Here we newly obtained eight small subunit (SSU) rDNA, 24 large subunit (LSU) rDNA, and 25 ITS1-5.8S-ITS2 sequences (covering 25 species, 10 genera and 4 out of 6 families) to analyze the phylogenetic relationships within Karyorelictea. Our results indicate that: (1) considering its unique morphology and early branching position in the SSU rDNA-based tree, the family Wilbertomorphidae represents a new taxon at order level, hence the new order Wilbertomorphida n. ord. is established; (2) all five families with available molecular information are monophyletic, as expected, and the orders Loxodida and Protostomatida show a closer relationship than with Protoheterotrichida; (3) in Trachelocercidae, the compound circumoral kineties is believed to be a plesiomorphic feature while the single circumoral kinety is synapomorphic; and (4) the freshwater genus Loxodes could be derived from the marine Remanella and both share most morphological features. Taken together, these muti-gene analyses provide further insights into the phylogeny of the diverse clades in Karyorelictea.

Combined Genome and Transcriptome Analyses of the Ciliate Schmidingerella arcuata (Spirotrichea) Reveal Patterns of DNA Elimination, Scrambling, and Inversion

Schmidingerella arcuata is an ecologically important tintinnid ciliate that has long served as a model species in plankton trophic ecology. We present a partial micronuclear genome and macronuclear transcriptome resource for S. arcuata, acquired using single-cell techniques, and we report on pilot analyses including functional annotation and genome architecture. Our analysis shows major fragmentation, elimination, and scrambling in the micronuclear genome of S. arcuata. This work introduces a new nonmodel genome resource for the study of ciliate ecology and genomic biology and provides a detailed functional counterpart to ecological research on S. arcuata.

Programmed genome rearrangements in ciliates

Ciliates are a highly divergent group of unicellular eukaryotes with separate somatic and germline genomes found in distinct dimorphic nuclei. This characteristic feature is tightly linked to extremely laborious developmentally regulated genome rearrangements in the development of a new somatic genome/nuclei following sex. The transformation from germline to soma genome involves massive DNA elimination mediated by non-coding RNAs, chromosome fragmentation, as well as DNA amplification. In this review, we discuss the similarities and differences in the genome reorganization processes of the model ciliates Paramecium and Tetrahymena (class Oligohymenophorea), and the distantly related Euplotes, Stylonychia, and Oxytricha (class Spirotrichea).

PhyloToL: A Taxon/Gene-Rich Phylogenomic Pipeline to Explore Genome Evolution of Diverse Eukaryotes

Estimating multiple sequence alignments (MSAs) and inferring phylogenies are essential for many aspects of comparative biology. Yet, many bioinformatics tools for such analyses have focused on specific clades, with greatest attention paid to plants, animals, and fungi. The rapid increase in high-throughput sequencing (HTS) data from diverse lineages now provides opportunities to estimate evolutionary relationships and gene family evolution across the eukaryotic tree of life. At the same time, these types of data are known to be error-prone (e.g., substitutions, contamination). To address these opportunities and challenges, we have refined a phylogenomic pipeline, now named PhyloToL, to allow easy incorporation of data from HTS studies, to automate production of both MSAs and gene trees, and to identify and remove contaminants. PhyloToL is designed for phylogenomic analyses of diverse lineages across the tree of life (i.e., at scales of >100 My). We demonstrate the power of PhyloToL by assessing stop codon usage in Ciliophora, identifying contamination in a taxon- and gene-rich database and exploring the evolutionary history of chromosomes in the kinetoplastid parasite Trypanosoma brucei, the causative agent of African sleeping sickness. Benchmarking PhyloToL’s homology assessment against that of OrthoMCL and a published paper on superfamilies of bacterial and eukaryotic organellar outer membrane pore-forming proteins demonstrates the power of our approach for determining gene family membership and inferring gene trees. PhyloToL is highly flexible and allows users to easily explore HTS data, test hypotheses about phylogeny and gene family evolution and combine outputs with third-party tools (e.g., PhyloChromoMap, iGTP).

De novo Sequencing, Assembly, and Annotation of the Transcriptome for the Free-Living Testate Amoeba Arcella intermedia

Arcella, a diverse understudied genus of testate amoebae is a member of Tubulinea in Amoebozoa group. Transcriptomes are a powerful tool for characterization of these organisms as they are an efficient way of characterizing the protein-coding potential of the genome. In this work, we employed both single-cell and clonal populations transcriptomics to create a reference transcriptome for Arcella. We compared our results with annotations of Dictyostelium discoideum, a model Amoebozoan. We assembled a pool of 38 Arcella intermedia transcriptomes, which after filtering are composed of a total of 14,712 translated proteins. There are GO categories enriched in Arcella including mainly intracellular signal transduction pathways; we also used KEGG to annotate 11,546 contigs, which also have similar distribution to Dictyostelium. A large portion of data is still impossible to assign to a gene family, probably due to a combination of lineage-specific genes, incomplete sequences in the transcriptome and rapidly evolved genes. Some absences in pathways could also be related to low expression of these genes. We provide a reference database for Arcella, and we highlight the emergence of the need for further gene discovery in Arcella.

Single-Cell Transcriptomics Reveal a Correlation between Genome Architecture and Gene Family Evolution in Ciliates

Ciliates, a eukaryotic clade that is over 1 billion years old, are defined by division of genome function between transcriptionally inactive germline micronuclei and functional somatic macronuclei. To date, most analyses of gene family evolution have been limited to cultivable model lineages (e.g., Tetrahymena, Paramecium, Oxytricha, and Stylonychia). Here, we focus on the uncultivable Karyorelictea and its understudied sister class Heterotrichea, which represent two extremes in genome architecture. Somatic macronuclei within the Karyorelictea are described as nearly diploid, while the Heterotrichea have hyperpolyploid somatic genomes. Previous analyses indicate that genome architecture impacts ciliate gene family evolution as the most diverse and largest gene families are found in lineages with extensively processed somatic genomes (i.e., possessing thousands of gene-sized chromosomes). To further assess ciliate gene family evolution, we analyzed 43 single-cell transcriptomes from 33 ciliate species representing 10 classes. Focusing on conserved eukaryotic genes, we use estimates of transcript diversity as a proxy for the number of paralogs in gene families among four focal clades: Karyorelictea, Heterotrichea, extensive fragmenters (with gene-size somatic chromosomes), and non-extensive fragmenters (with more traditional somatic chromosomes), the latter two within the subphylum Intramacronucleata. Our results show that (i) the Karyorelictea have the lowest average transcript diversity, while Heterotrichea are highest among the four groups; (ii) proteins in Karyorelictea are under the highest functional constraints, and the patterns of selection in ciliates may reflect genome architecture; and (iii) stop codon reassignments vary among members of the Heterotrichea and Spirotrichea but are conserved in other classes.IMPORTANCE To further our understanding of genome evolution in eukaryotes, we assess the relationship between patterns of molecular evolution within gene families and variable genome structures found among ciliates. We combine single-cell transcriptomics with bioinformatic tools, focusing on understudied and uncultivable lineages selected from across the ciliate tree of life. Our analyses show that genome architecture correlates with patterns of protein evolution as lineages with more canonical somatic genomes, such as the class Karyorelictea, have more conserved patterns of molecular evolution compared to other classes. This study showcases the power of single-cell transcriptomics for investigating genome architecture and evolution in uncultivable microbial lineages and provides transcriptomic resources for further research on genome evolution.

Evolutionary origins and impacts of genome architecture in ciliates

Genome architecture is well diversified among eukaryotes in terms of size and content, with many being radically shaped by ancient and ongoing genome conflicts with transposable elements (e.g., the large transposon‐rich genomes common among plants). In ciliates, a group of microbial eukaryotes with distinct somatic and germ‐line genomes present in a single cell, the consequences of these genome conflicts are most apparent in their developmentally programmed genome rearrangements. This complicated developmental phenomenon has largely overshadowed and outpaced our understanding of how germ‐line and somatic genome architectures have influenced the evolutionary dynamism and potential in these taxa. In our review, we highlight three central concepts: how the evolution of atypical ciliate germ‐line genome architectures is linked to ancient genome conflicts; how the complex, epigenetically guided transformation of germline to soma during development can generate widespread genetic variation; and how these features, coupled with their unusual life cycle, have increased the rate of molecular evolution linked to genome architecture in these taxa.

Twisted Tales: Insights into Genome Diversity of Ciliates Using Single-Cell 'Omics

The emergence of robust single-cell 'omics techniques enables studies of uncultivable species, allowing for the (re)discovery of diverse genomic features. In this study, we combine single-cell genomics and transcriptomics to explore genome evolution in ciliates (a > 1 Gy old clade). Analysis of the data resulting from these single-cell 'omics approaches show: 1) the description of the ciliates in the class Karyorelictea as "primitive" is inaccurate because their somatic macronuclei contain loci of varying copy number (i.e., they have been processed by genome rearrangements from the zygotic nucleus); 2) gene-sized somatic chromosomes exist in the class Litostomatea, consistent with Balbiani's (1890) observation of giant chromosomes in this lineage; and 3) gene scrambling exists in the underexplored Postciliodesmatophora (the classes Heterotrichea and Karyorelictea, abbreviated here as the Po-clade), one of two major clades of ciliates. Together these data highlight the complex evolutionary patterns underlying germline genome architectures in ciliates and provide a basis for further exploration of principles of genome evolution in diverse microbial lineages.

Nuclear Features of the Heterotrich Ciliate Blepharisma americanum: Genomic Amplification, Life Cycle, and Nuclear Inclusion

Blepharisma americanum, a member of the understudied ciliate class Heterotrichea, has a moniliform somatic macronucleus that resembles beads on a string. Blepharisma americanum is distinguishable by its pink coloration derived from the autofluorescent pigment blepharismin and tends to have a single somatic macronucleus with 3-6 nodes and multiple germline micronuclei. We used fluorescence confocal microscopy to explore the DNA content and amplification between the somatic and germline nuclei of B. americanum through its life cycle. We estimate that the DNA content of the macronucleus and micronucleus are 43 ± 8 Gbp and 83 ± 16 Mbp respectively. This correlates with an approximate DNA content difference of 500-fold from micronucleus to macronucleus and a macronuclear ploidy of ~1,100 N as compared to the presumably diploid micronucleus. We also investigate a previously reported macronuclear inclusion, which is present sporadically across all life cycle stages; this inclusion looks as if it contains blepharismin based on its fluorescent properties, but its function remains unknown. We also provide additional detail to our understanding of life cycles changes in B. americanum by analyses of fluorescent images. Overall, the data analyzed here contribute to our understanding of the diversity of nuclear architecture in ciliates by providing details on the highly polyploid somatic macronucleus of B. americanum.

Nuclear Architecture and Patterns of Molecular Evolution Are Correlated in the Ciliate Chilodonella uncinata

The relationship between nuclear architecture and patterns of molecular evolution in lineages across the eukaryotic tree of life is not well understood, partly because molecular evolution is traditionally explored as changes in base pairs along a linear sequence without considering the context of nuclear position of chromosomes. The ciliate Chilodonella uncinata is an ideal system to address the relationship between nuclear architecture and patterns of molecular evolution as the somatic macronucleus of this ciliate is composed of a peripheral DNA-rich area (orthomere) and a DNA-poor central region (paramere) to form a “heteromeric” macronucleus. Moreover, because the somatic chromosomes of C. uncinata are highly processed into “gene-sized” chromosomes (i.e., nanochromosomes), we can assess fine-scale relationships between location and sequence evolution. By combining fluorescence microscopy and analyses of transcriptome data from C. uncinata, we find that highly expressed genes have the greatest codon usage bias and are enriched in DNA-poor regions. In contrast, genes with less biased sequences tend to be concentrated in DNA abundant areas, at least during vegetative growth. Our analyses are consistent with recent work in plants and animals where nuclear architecture plays a role in gene expression. At the same time, the unusual localization of nanochromosomes suggests that the highly structured nucleus in C. uncinata may create a “gene bank” that facilitates rapid changes in expression of genes required only in specific life history stages. By using “nonmodel” organisms like C. uncinata, we can explore the universality of eukaryotic features while also providing examples of novel properties (i.e., the presence of a gene bank) that build from these features.

An epigenetic toolkit allows for diverse genome architectures in eukaryotes

Genome architecture varies considerably among eukaryotes in terms of both size and structure (e.g. distribution of sequences within the genome, elimination of DNA during formation of somatic nuclei). The diversity in eukaryotic genome architectures and the dynamic processes are only possible due to the well-developed epigenetic toolkit, which probably existed in the Last Eukaryotic Common Ancestor (LECA). This toolkit may have arisen as a means of navigating the genomic conflict that arose from the expansion of transposable elements within the ancestral eukaryotic genome. This toolkit has been coopted to support the dynamic nature of genomes in lineages across the eukaryotic tree of life. Here we highlight how the changes in genome architecture in diverse eukaryotes are regulated by epigenetic processes, such as DNA elimination, genome rearrangements, and adaptive changes to genome architecture. The ability to epigenetically modify and regulate genomes has contributed greatly to the diversity of eukaryotes observed today.