Comparative genomics reveals insight into the evolutionary origin of massively scrambled genomes

Abstract
Editor's evaluation
Introduction
Results
Discussion
Methods
Data availability
References
Article and author information
Metrics

Abstract

Ciliates are microbial eukaryotes that undergo extensive programmed genome rearrangement, a natural genome editing process that converts long germline chromosomes into smaller gene-rich somatic chromosomes. Three well-studied ciliates include Oxytricha trifallax, Tetrahymena thermophila, and Paramecium tetraurelia, but only the Oxytricha lineage has a massively scrambled genome, whose assembly during development requires hundreds of thousands of precisely programmed DNA joining events, representing the most complex genome dynamics of any known organism. Here we study the emergence of such complex genomes by examining the origin and evolution of discontinuous and scrambled genes in the Oxytricha lineage. This study compares six genomes from three species, the germline and somatic genomes for Euplotes woodruffi, Tetmemena sp., and the model ciliate O. trifallax. We sequenced, assembled, and annotated the germline and somatic genomes of E. woodruffi, which provides an outgroup, and the germline genome of Tetmemena sp. We find that the germline genome of Tetmemena is as massively scrambled and interrupted as Oxytricha’s: 13.6% of its gene loci require programmed translocations and/or inversions, with some genes requiring hundreds of precise gene editing events during development. This study revealed that the earlier diverged spirotrich, E. woodruffi, also has a scrambled genome, but only roughly half as many loci (7.3%) are scrambled. Furthermore, its scrambled genes are less complex, together supporting the position of Euplotes as a possible evolutionary intermediate in this lineage, in the process of accumulating complex evolutionary genome rearrangements, all of which require extensive repair to assemble functional coding regions. Comparative analysis also reveals that scrambled loci are often associated with local duplications, supporting a gradual model for the origin of complex, scrambled genomes via many small events of DNA duplication and decay.

Editor's evaluation

The study marks a significant advance in the field of evolutionary genomics of ciliates, an ancient and highly diverse eukaryotic phylum with many idiosyncrasies that teach us valuable lessons, inter alia, about sex and the plasticity of genomes. By focusing on two species from the same family, plus a more distant outgroup within the same class, this valuable study provides new and compelling information on evolutionary trends of genome rearrangement among different species of this interesting group of organisms. The work will be of interest to anyone interested in genome dynamics.

https://doi.org/10.7554/eLife.82979.sa0

Introduction

Organisms do not always contain a single, static genome. Programmed genome editing is a naturally occurring and essential part of development in many organisms, including ciliates (Chen et al., 2014), nematodes (Mitreva et al., 2005), lampreys (Smith et al., 2012), and zebra finches (Biederman et al., 2018). Most of these events involve precise removal and rejoining of large regions of DNA during postzygotic differentiation of a somatic genome from a germline genome. Ciliates are microbial eukaryotes with two types of nuclei: a somatic macronucleus (MAC) that differentiates from a germline micronucleus (MIC). In the model ciliate Oxytricha, the MAC is entirely active chromatin (Beh et al., 2019) and the hub of transcription. The three species that we compare are all spirotrichs, which have gene-sized ‘nanochromosomes’ in the MAC, present at high copy number (Swart et al., 2013; Chen et al., 2015; Wang et al., 2016; Chen et al., 2019; Lindblad et al., 2019; Vinogradov et al., 2012). The diploid MIC participates in sexual reproduction, but its megabase-sized chromosomes are mostly transcriptionally silent.

Gene loci are often arranged discontinuously in the MIC, with short genic segments called macronuclear destined sequences (MDSs), interrupted by stretches of non-coding DNA called internally eliminated sequences (IESs) (Figure 1A). During sexual development, a new MAC genome rearranges from a copy of the zygotic MIC genome. MDSs join in the correct order and orientation, whereas MIC-limited genomic regions undergo programmed deletion, including repetitive elements, intergenic regions, and IESs (Figure 1A). Though analogous to intron splicing, these events occur on DNA. The MDSs for some MAC chromosomes are scrambled if they require translocation or inversion during MAC development (Figure 1A). Pairs of short repeats, called pointers, are present at MDS-IES junctions in both scrambled and nonscrambled loci (Mitcham et al., 1992; Prescott, 1994). Pointer sequences are present twice in the MIC, at the end of MDS n and the beginning of MDS n+1. One copy of the repeat is retained at each MDS-MDS junction in a mature MAC chromosome (Figure 1A). These microhomologous regions help guide MDS recombination, but most are non-unique, and the shortest pointers are just 2 bp. Thousands of long, noncoding template RNAs collectively program MDS joining (Nowacki et al., 2008; Lindblad et al., 2017; Yerlici and Landweber, 2014).

Figure 1

Download asset Open asset

Genome rearrangements in representative ciliate species.

(A) Diagram of genome rearrangement in *Oxytricha*. Each ciliate cell contains a somatic macronucleus (MAC) and a germline micronucleus (MIC). During development, the MAC genome rearranges from a copy of the MIC genome. (1) Nonscrambled genes rearrange simply by joining consecutive macronuclear destined sequences (MDSs, blue boxes) and removing internal eliminated sequences (IESs, thin lines). (2) Rearrangement of scrambled genes requires MDS translocation and/or inversion. Pointers are microhomologous sequences (colored vertical bars) present in two copies in the MIC and only one copy in the MAC where consecutive MDSs recombine. (B) Comparison of genome rearrangement features of representative ciliates and the non-ciliate *Plasmodium falciparum* as an outgroup (phylogenetic information is based on Parfrey et al., 2011; Bracht et al., 2013). Conclusions from this study are shown in bold. * indicates that some scrambled pointers in *Euplotes woodruffi* are much longer, as discussed in the results. Statistics for pointers ≤30 bp in *E. woodruffi* are shown. Table information derives from the following sources: 1 - Swart et al., 2013; 2 - Lindblad et al., 2019; 3 - Chen et al., 2014; 4 - Chen et al., 2015; 5 - Sheng et al., 2020; 6 - Eisen et al., 2006; 7 - Hamilton et al., 2016; 8 - Aury et al., 2006; 9 - Guérin et al., 2017; 10 - Arnaiz et al., 2012; 11 - Riley and Katz, 2001; 12 - Maurer-Alcalá et al., 2018a; 13 - Katz and Kovner, 2010; 14 - Gao et al., 2014.

Numerous studies have inferred the possible scope of genome rearrangement in different ciliate species using partial genome surveys. In Paramecium, PiggyMac-depleted cells fail to remove MIC-limited regions properly, which provided a resource to annotate ~45,000 IESs prior to assembly of a draft MIC genome (Arnaiz et al., 2012). The use of single-cell sequencing has allowed pilot studies to sample partial MIC genomes of diverse species (Chen et al., 2019; Maurer-Alcalá et al., 2018a; Maurer-Alcalá et al., 2018b; Smith et al., 2020). Alignment of tentative MIC reads to either assembled MAC genomes or single-cell transcriptome data predicts over 20 candidate scrambled loci in two basal ciliates, Loxodes sp. and Blepharisma americanum (Maurer-Alcalá et al., 2018b) and hundreds of candidate loci in the tintinnid Schmidingerella arcuata (Smith et al., 2020). Nearly one-third (31%) of approximately 5000 surveyed transcripts may be scrambled in Chilodonella uncinata (Maurer-Alcalá et al., 2018a, Figure 1B), which has four confirmed cases of scrambled genes (Katz and Kovner, 2010; Gao et al., 2014). Transcriptome-based surveys offer less precise estimates and cannot distinguish RNA splicing. Several computational pipelines have been developed to facilitate the inference of genome rearrangement features by split-read mapping in the absence of complete MIC or MAC reference genomes (Denby Wilkes et al., 2016; Zheng et al., 2020; Feng et al., 2020; Seah et al., 2021). By surveying lighter genome coverage prior to full sequencing, these tools provide partial insight into germline architecture. This helps guide selection of species for full genome sequencing and subsequent construction of complete rearrangement maps between the MIC and MAC genomes. High-quality MIC genome reference assemblies are only currently available for three ciliate genera: Oxytricha (Chen et al., 2014), Tetrahymena (Hamilton et al., 2016), and Paramecium (Guérin et al., 2017; Sellis et al., 2021).

Programmed genome rearrangements in Oxytricha exhibit the highest accuracy and largest scale of any known natural gene-editing system, with exquisite control over hundreds of thousands of precise DNA cleavage/joining events. Accordingly, its germline genome structure is arguably the most complex of any model organism (Chen et al., 2014), requiring programmed deletion of over 90% of the germline DNA during development and massive descrambling of the resulting fragments to construct a new MAC genome of over 18,000 chromosomes (Lindblad et al., 2019). This differs from the distantly related Tetrahymena and Paramecium that both eliminate ~30% of the germline genome (Hamilton et al., 2016; Guérin et al., 2017). Paramecium uses exclusively 2 bp pointers and lacks evidence of any scrambled loci. A small number of scrambled loci (4 confirmed out of 2711 candidates) have been reported in Tetrahymena (Sheng et al., 2020, Figure 1B). Tetrahymena and Paramecium diverged from Oxytricha over 1 billion years ago (Parfrey et al., 2011; Bracht et al., 2013), which leaves a large gap in our understanding of the emergence of complex DNA rearrangements in the Oxytricha lineage.

Open questions include how did the Oxytricha germline genome acquire its high number of IES insertions and how do scrambled loci arise and evolve. Three previous studies tackled these questions at the level of single genes and orthologs, including DNA polymerase α, actin I, and Telomere end-binding protein subunit α (Hogan et al., 2001; Chang et al., 2005; Wong and Landweber, 2006; DuBois and Prescott, 1995). Here, we provide the first comparative genomic analysis of Oxytricha trifallax and two other spirotrichous ciliates, Tetmemena sp. and Euplotes woodruffi. Tetmemena sp. is a hypotrich similar to Tetmemena pustulata, formerly Stylonychia pustulata (Chen et al., 2015), in the same family as O. trifallax (Figure 1B; Chen et al., 2014; Chen et al., 2015). Hypotrichs are noted for the presence of scrambled genes, based on previous ortholog comparisons (Chen et al., 2015; Hogan et al., 2001; Chang et al., 2005; DuBois and Prescott, 1995; Figure 1B). E. woodruffi, together with the hypotrichous ciliates, belong to the class Spirotrichea (Figure 1B). Like hypotrichs, Euplotes also has gene-sized nanochromosomes in the MAC genome (Wang et al., 2016; Chen et al., 2019; Chen et al., 2021), but this outgroup uses a different genetic code (UGA is reassigned to cysteine, Meyer et al., 1991), and little is known about its MIC genome. A partial MIC genome of Euplotes vannus was previously assembled, and it contains highly conserved TA pointers (Chen et al., 2019), consistent with previous observations in Euplotes crassus (Klobutcher and Herrick, 1995). This differs from O. trifallax, which uses longer pointers of varying lengths, with scrambled pointers typically longer than nonscrambled ones (Chen et al., 2014, Figure 1B). This observation suggests that longer pointers may supply more information to facilitate MDS descrambling, sometimes over great distances. Therefore, the preponderance of 2 bp pointers in the other Euplotes species could indicate limited capacity to support scrambled genes, and a partial genome survey of E. vannus concluded that at least 97% of loci are nonscrambled (Chen et al., 2019). Early studies of Euplotes octocarinatus, on the other hand, demonstrated its use of longer pointers (that usually contain TA) (Tan et al., 1999; Wang et al., 2005), suggesting that some members of the Euplotes genus may have the capacity to support complex genome reorganization. To investigate the origin of scrambled genomes, we choose E. woodruffi as an outgroup, because it is closely related to E. octocarinatus (Syberg-Olsen et al., 2016) and feasible to culture in the lab.

This study includes the de novo assemblies of the micronuclear genome of Tetmemena sp. and both genomes of E. woodruffi. The availability of MIC and MAC genomes for both species allows us to annotate and compare their genome rearrangement maps and other key features to each other and to O. trifallax. The MIC genome of Tetmemena is extremely interrupted, like Oxytricha. While the E. woodruffi MIC genome is much more IES-sparse, it contains thousands of scrambled genes, whose architecture we compare to orthologous loci in the other species. We infer that the evolutionary origin of scrambled genes is associated with local duplications, providing strong support for a previously proposed simple evolutionary model requiring only duplication and decay (Gao et al., 2015) that allows for the evolutionary expansion of extremely rearranged chromosome architectures.

Results

Germline genome expansion via repetitive elements

Tetmemena sp. and E. woodruffi were both propagated in laboratory culture from single cells. The E. woodruffi MAC genome was sequenced and assembled from paired-end Illumina reads from whole cell DNA, which is mostly MAC-derived. For comparative analysis, the MAC genome of E. woodruffi was assembled using the same pipeline previously used for Tetmemena sp. (Chen et al., 2015). Because MIC DNA is significantly more sparse than MAC DNA in individual cells (Prescott, 1994), MIC DNA was enriched before sequencing (see Methods); however, this leads to much lower sequence coverage of the MIC than the MAC. Third-generation long reads (Pacific Biosciences and Oxford Nanopore Technologies) were combined with Illumina paired-end reads (Methods, see genome coverage in Supplementary file 1) to construct hybrid genome assemblies for Tetmemena sp. and E. woodruffi. Though the final genome assemblies are still fragmented, often due to transposon or other repetitive insertions at boundaries (Figure 2—figure supplement 1), the current draft assemblies cover most (>90%) MDSs for 89.1% of MAC nanochromosomes in Tetmemena, and for 90.0% of MAC nanochromosomes in E. woodruffi. This allowed us to establish near-complete rearrangement maps for the newly assembled genomes of Tetmemena and E. woodruffi, at a level comparable to the published reference for O. trifallax (Chen et al., 2014), which is appropriate for comparative analysis.

Table 1 shows a comparison of genome features for the three species. The three MAC genomes are similar in size, with most nanochromosomes bearing only one gene. The size distributions of MAC chromosomes are similar for the three species, though slightly shorter for E. woodruffi, consistent with prior observation via gel electrophoresis (Prescott, 1994, Figure 2—figure supplement 2). Like O. trifallax (Swart et al., 2013), the maximum number of genes encoded on one chromosome is 7–8 (Table 1). Surprisingly, the MIC genome sizes differ substantially: the Tetmemena MIC genome assembly is 237 Mbp, nearly half that of Oxytricha. The E. woodruffi MIC genome assembly is even smaller, approximately 172 Mbp (Table 1).

Table 1

Statistics of somatic macronucleus (MAC) and germline micronucleus (MIC) genomes in three species.

	Oxytricha trifallax		Tetmemena sp.		Euplotes woodruffi
	MAC^a,*	MIC^b	MAC^c	MIC^†	MAC^†	MIC^†
Genome size (Mbp)	67.1	496	60.6	237	72.2	172
N50 (bp)	3745	27,807	3339	14,722	2702	44,656
GC%	31.36	28.44	37.05	32.17	36.56	35.31
Number of contigs^‡	22,426	25,720	25,206	28,446	35,099	17,655
Two-telomere contigs	14,225	-	15,802	-	19,061	-
Telomeric contigs	20,336	-	21,165	-	28,294	-
Single-gene telomeric contigs	76.1%	-	75.5%	-	68.5%	-
Maximum number of genes on a telomeric contig	8	-	7	-	8	-

a - Swart et al., 2013; b - Chen et al., 2014; c - Chen et al., 2015.
*

This study used the MAC genome of Oxytricha from Swart et al., 2013 instead of the long-read assembly in Lindblad et al., 2019, because the short MAC genomes in the present study were primarily assembled from Illumina reads, as in Swart et al., 2013. Lindblad et al., 2019 updated Swart et al., 2013 by including nanochromosomes captured in single long reads, which are currently not available for the other two species. The MIC genomes of Tetmemena and E. woodruffi were assembled to a similar N50 as the reference O. trifallax genome (Chen et al., 2014) for comparative analysis.
†

Data from this study.
‡

Telomere-bearing element (TBE) transposon contaminants in MAC contigs were removed (Methods). Therefore, 24 Oxytricha MAC contigs and 13 Tetmemena MAC contigs were removed from the published versions.

The expansion of repetitive elements in the Oxytricha lineage may contribute to the difference in MIC genome sizes (Figure 2A–C). Oxytricha has a variety of tranposable elements (TEs) in the MIC, with telomere-bearing elements (TBEs) of the Tc1/mariner family the most abundant (Chen et al., 2014; Chen and Landweber, 2016, Supplementary file 2). A complete TBE transposon contains three open reading frames (ORFs). ORF1 encodes a 42kD transposase with a DDE-catalytic motif. Though present only in the germline, TBEs are so abundant in hypotrichs that some were partially recovered and assembled from whole cell DNA (Chen and Landweber, 2016). The Oxytricha MIC genome contains ~10,000 complete TBEs and ~24,000 partial TBEs, which occupy approximately 15.20% (75 Mbp) of the genome (Figure 2A, Supplementary file 3; Chen et al., 2014; Chen and Landweber, 2016). Tetmemena, on the other hand, has many fewer TBE ORFs and only 48 complete TBEs (Supplementary file 3), comprising 1.83% (4.3 Mbp) of its MIC genome (Figure 2B). E. crassus has also been reported to have an abundant transposon family called Tec elements (Transposon of Euplotes crassus). Like TBEs, each Tec consists of three ORFs, and ORF1 also encodes a transposase from the Tc1/mariner family (Baird et al., 1989; Krikau and Jahn, 1991; Jahn et al., 1993; Jahn et al., 1989; Klobutcher and Herrick, 1997). The ~57 kD ORF2 encodes a tyrosine-type recombinase (Doak et al., 2003), and the 20kD ORF3 has unknown function (Jahn et al., 1993). Using the three ORFs of Tec1 and Tec2 as queries for search, we identified 74 complete Tec elements in E. woodruffi. Collectively, Tec ORFs occupy 3.6 Mbp, corresponding to only 2.1% of the MIC genome (Figure 2C). Notably, the transposase-encoding ORF1 is more abundant than the other two TBE/Tec ORFs in all three ciliates (Supplementary file 3), consistent with its proposed role in DNA cleavage during genome rearrangement in Oxytricha (Nowacki et al., 2009).

Figure 2 with 2 supplements see all

Download asset Open asset

The three germline micronucleus (MIC) genomes differ in repeat content, especially transposable elements.

(**A–C**) MIC genome categories for (A) *Oxytricha trifallax*, (B) *Tetmemena sp*., and (C) *Euplotes woodruffi. Oxytricha* displays the greatest proportion of repetitive elements (telomere-bearing elements [TBE], other repeats, and tandem repeats) relative to the other species. *Oxytricha* MIC-specific genes were annotated in Chen et al., 2014; Miller et al., 2021. (**D–F**) Phylogenetic analysis of the three TBE open reading frames (ORFs) in *Oxytricha* and *Tetmemena*: (D) 42 kD, (E) 22 kD, and (F) 57 kD, suggest that TBE3 (green) is the ancestral transposon family in *Oxytricha*. For each ORF, 30 protein sequences from each species were randomly subsampled and maximum likelihood trees constructed using PhyML (Guindon et al., 2010).

Oxytricha contains three families of TBEs. TBE3 appears to be the most ancient among hypotrichs, based on previous analysis of limited MIC genome data (Chen and Landweber, 2016). We constructed phylogenetic trees using randomly subsampled TBE sequences for all three ORFs from Oxytricha and Tetmemena (Figure 2D–F). This confirmed that only TBE3 is present in the Tetmemena MIC genome, as proposed in Chen and Landweber, 2016. This also suggests that TBE1 and TBE2 expanded in Oxytricha after its divergence from other hypotrichous ciliates. As illustrated in Figure 2—figure supplement 1, the MIC genome contexts of TBEs in Oxytricha and Tetmemena are similar, with many TE insertions within IESs, consistent with either IESs as hotspots for TE insertion or with the model (Klobutcher and Herrick, 1997) that some TE insertions may have generated IESs, as demonstrated in Paramecium (Sellis et al., 2021; Feng and Landweber, 2021). Subsequent sequence evolution at the edges of IES/MDS pointers (DuBois and Prescott, 1995) can give rise to boundaries that no longer correspond precisely to TBE ends. For further discussion of the conservation of TBE locations, see the section, ‘Oxytricha and Tetmemena share conserved rearrangement junctions’ below.

Additionally, Repeatmodeler/Repeatmasker identified that Oxytricha has more MIC repeats in the ‘Other’ category than Tetmemena or E. woodruffi (Figure 2, subcategories of repeat content in Supplementary file 2). 214 Mbp of the Oxytricha MIC genome (43%, which is greater than 35.9% reported in Chen et al., 2014 that used earlier versions of the software) is considered repetitive (including TBEs, tandem repeats, and other repeats in Figure 2), versus 31.7 Mbp for Tetmemena (13.4%) and 28.5 Mbp (16.8%) for E. woodruffi. Oxytricha’s additional ~180 Mbp in repeat content partially explains the significantly larger MIC genome size of Oxytricha versus the other spirotrich ciliates.

The E. woodruffi genome has fewer IESs

We used the genome rearrangement annotation tool, Scrambled DNA Rearrangement Annotation Protocol (SDRAP, Braun et al., 2022) to annotate the MIC genomes of Oxytricha, Tetmemena, and E. woodruffi (Methods). Consistent with their close genetic distance, the genomes of O. trifallax and Tetmemena have similarly high levels of discontinuity (Figure 3A). We annotated over 215,299 MDSs in Oxytricha and over 215,624 in Tetmemena with similar MDS length distributions (Figure 3A). By contrast, E. woodruffi MDSs are typically longer, which indicates a less interrupted genome (Figure 3A). We compared the number of MDSs between single-copy orthologs for single-gene MAC chromosomes across the three species and found that the orthologs have similar coding sequence (CDS) lengths (Figure 3—figure supplement 1A–B). There is a strong positive correlation between number of MDSs for orthologous genes in Oxytricha and Tetmemena (R²=0.75, Figure 3B). There is no correlation among number of MDSs between orthologs of E. woodruffi and Oxytricha (R²=0.003, Figure 3C), since E. woodruffi orthologs typically contain fewer MDSs.

Figure 3 with 2 supplements see all

Download asset Open asset

The three germline micronucleus genomes are interrupted by internally eliminated sequences (IESs) at different levels.

(A) Macronuclear destined sequences (MDSs) of *Euplotes woodruffi* are longer compared to *Oxytricha* or *Tetmemena*. (B) Positive correlation between the numbers of MDSs for orthologous genes in *Tetmemena* and in *Oxytricha* for 903 single-gene orthologs. Black line is the function of linear regression (R²=0.75). Red line is y=x. (C) Orthologs in *E. woodruffi* have fewer MDSs compared to *Oxytricha*, with no correlation (R²=0.003). Note that many highly discontinuous genes in *Oxytricha* are IES-less in *E. woodruffi* (present on one MDS). 917 single-gene orthologs are shown. (D) Distribution of pointers on single-gene somatic macronucleus (MAC) chromosomes in *Oxytricha vs*. (E) *E. woodruffi*, with MAC chromosomes oriented in gene direction. Pointers significantly accumulate at the 5’ end of single-gene MAC chromosomes in *E. woodruffi*. (F) Pointer positions on 3684 two-MDS MAC chromosomes demonstrate a preference upstream of the start codon.

The E. woodruffi genome is generally much less interrupted than that of Oxytricha or Tetmemena. 39.9% of MAC nanochromosomes in E. woodruffi lack IESs (IES-less nanochromosomes) compared to only 4.1 and 4.4% in Oxytricha and Tetmemena, respectively. The sparse IES distribution (as measured by plotting pointer distributions) in E. woodruffi displays a curious 5’ end bias on single-gene MAC chromosomes, oriented in gene direction (Figure 3E). A weak 5’ bias is also present in Oxytricha (Figure 3D) and Tetmemena (Figure 3—figure supplement 1C). In addition, E. woodruffi IESs preferentially accumulate in the 5’ UTR, a short distance upstream of start codons (Figure 3F). Notably, the median distance between the 5’ telomere addition site and the start codon in E. woodruffi is just 54 bp for single-gene chromosomes, approximately half that of Oxytricha (Swart et al., 2013).

E. woodruffi has an intermediate level of genome scrambling

Scrambled genome rearrangements exist in all three species, which we report here for the first time in Tetmemena and the early diverged E. woodruffi. Previous studies have described scrambled genes with confirmed MIC-MAC rearrangement maps for a limited species of hypotrichs (Chen et al., 2014; Chen et al., 2015; Hogan et al., 2001; Chang et al., 2005; Wong and Landweber, 2006; DuBois and Prescott, 1995) and Chilodonella (Katz and Kovner, 2010; Gao et al., 2014) but not in Euplotes. Consistent with the phylogenetic placement of Euplotes as an earlier diverged outgroup to hypotrichs (Lynn, 2008; Gao et al., 2016), the E. woodruffi genome is scrambled, but it contains approximately half as many scrambled genes (2429 genes encoded on 1913 chromosomes, or 7.3% of genes), versus 15.6% scrambled in O. trifallax (3613 genes encoded on 2852 chromosomes) and 13.6% in Tetmemena (3371 genes encoded on 2556 chromosomes). The E. woodruffi lineage may therefore reflect an evolutionary intermediate stage between ancestral genomes with only modest levels of genome scrambling and the more massively scrambled genomes of hypotrichs.

We infer that many genes were likely scrambled in the last common ancestor of Oxytricha and Tetmemena, because these two species share approximately half of their scrambled genes (Supplementary file 4). Furthermore, most scrambled genes are not new genes, since they possess at least one ortholog in other ciliate species (Supplementary file 4, Supplementary file 5).

Scrambled genes are associated with local paralogy

Notably, scrambled genes in all three species generally have more paralogs (Figure 4). We identified orthogroups containing genes derived from the same gene in the last common ancestor of the three species (Methods). For each species, orthogroups with at least one scrambled gene are significantly larger than those containing no scrambled genes (p-value <1e−5, Mann-Whitney U test, Figure 4A–C). This association suggests a possible role of gene duplication in the origin of scrambled genes.

Figure 4 with 3 supplements see all

Download asset Open asset

Scrambled genes have more paralogs than nonscrambled genes in the three species.

Orthogroups containing at least one scrambled gene (‘scrambled’) are larger than orthogroups that lack scrambled genes (‘nonscrambled’) in (A) *Oxytricha*, (B) *Tetmemena,* and (C) *Euplotes woodruffi*.

Scrambled pointers are generally longer than nonscrambled ones in all three species (Figure 3—figure supplement 2), consistent with prior observations (Chen et al., 2014) and the possibility that longer pointers participate in more complex rearrangements, including recombination between MDSs separated by greater distances (Landweber et al., 2000). Scrambled and nonscrambled IESs also differ in their length distribution (Figure 3—figure supplement 2). Curiously, scrambled ‘pointers’ in E. woodruffi can be as long as several hundred base pairs (median 48 bp, average 212 bp) unlike the more typical 2–20 bp canonical pointers. These long ‘pointers’ in E. woodruffi are more likely partial MDS duplications (Figure 4—figure supplement 1A). We also identified MDSs that map to two or more paralogous regions within the same MIC contig (Supplementary file 6), therefore representing MDS duplications and not alleles. Such paralogous regions could be alternatively incorporated into the rearranged MAC product. Moreover, we find that, for all three species, there are significantly more scrambled chromosomes than nonscrambled MAC chromosomes that contain at least one paralogous MDS (chi-square test, p-value <1e−10; Supplementary file 6). An example is shown in Figure 4—figure supplement 1A (MDS 7 and 7').

The presence of paralogous MDSs can contribute to the origin of scrambled rearrangements, as proposed in an elegant model by Gao et al., 2015; illustrated in Figure 4—figure supplement 1B. The model proposes that initial MDS duplications permit alternative use of either MDS copy into the mature MAC chromosome. As mutations accumulate in redundant paralogs, cells that incorporate the least decayed MDS regions into the MAC gene would have both a fitness advantage and a better match to the template RNA (Nowacki et al., 2008) that guides rearrangement, thus increasing the likelihood of incorporation into the MAC chromosome. The paralogous regions containing more mutations would gradually decay into IESs, and scrambled pointers eventually be reduced to a shorter length. The extended length ‘pointers’ that we identified in E. woodruffi may reflect an intermediate stage in the origin of scrambled genes (Figure 4—figure supplement 1B).

This model may generally explain the abundance and expansion of ‘odd-even’ patterns in ciliate scrambled genes (Landweber et al., 2000; Burns et al., 2016). As illustrated in Figure 4—figure supplement 1A, the even- and odd-numbered MDSs for many scrambled genes derive from different MIC genome clusters. The model predicts that the IES between MDS n−1 and n+1 often derives from ancestral duplication of a region containing MDS n (Figure 4—figure supplement 2A). To test this hypothesis explicitly, we extracted from all odd-even scrambled loci in the three species all sets of corresponding MDS/IES pairs that are flanked by identical pointers on both sides, i.e., all pairs of scrambled MDSs and IESs, where the IES between MDS n−1 and n+1 is directly exchanged for MDS n during DNA rearrangement (S1 and S2 in Figure 4—figure supplement 2A). To exclude the possibility of alleles confounding this analysis, MDS and IES pairs were only considered if they map to the same MIC contig. In E. woodruffi, the lengths of these MDS/IES pairs strongly correlate (Spearman correlation ρ=0.755, p<1e−5, Figure 4—figure supplement 2B). Moreover, many MDS and IES sequence pairs also share sequence similarity, consistent with paralogy: for 248 MDS-IES pairs of similar length, 90.3% share a core sequence with ~97.5% identity across 8–100% of both the IES and MDS length. The lowest end of these observations is also compatible with an alternative model (Chang et al., 2005) in which direct recombination between IESs and MDSs at short repeats can lead to expansion of odd-even patterns. For Oxytricha and Tetmemena, the MDS and IES lengths for such MDS/IES pairs also display a weakly-positive correlation (p-values and Spearman correlation ρ shown in Figure 4—figure supplement 2D–E). Remarkably, the odd-even-containing loci that are species-specific, and therefore became scrambled more recently, have the strongest length correlation (Figure 4—figure supplement 2C–E) and more pairs that display sequence similarity (Supplementary file 7) relative to older loci (scrambled in two or more species). This result is consistent with an evolutionary process in which mutations accumulate in one copy of the MDS, gradually obscuring its sequence homology and ability to be incorporated as a functional MDS, and eventually its ability to be recognized by the template RNAs that guide DNA rearrangement. This analysis also suggests that most of the odd-even scrambled loci in E. woodruffi arose recently, because there is greater sequence similarity between MDSs and the corresponding IESs that they replace. Conversely, we infer that most loci that are scrambled in both Oxytricha and Tetmemena became scrambled earlier in evolution, since they display weaker sequence similarity between exchanged MDS and IES regions.

Scrambled and nonscrambled genes display nearly identical expression support (the presence of at least one read in all three replicates) in both Oxytricha (Supplementary file 8) and Tetmemena. E. woodruffi has slightly more expression support for nonscrambled vs. scrambled genes (Figure 4—figure supplement 3), which could be explained by more recent acquisition of thousands of scrambled loci in E. woodruffi. In some of those cases the nonscrambled paralogs may still contribute the major function. The distribution of expression levels is similar for scrambled vs. nonscrambled genes in all three species, supporting their authenticity (Figure 4—figure supplement 3), although in a Mann-Whitney U test, the average expression level of three replicates is significantly higher in nonscrambled genes for Oxytricha and E. woodruffi, but not significant for Tetmemena.

Oxytricha and Tetmemena share conserved DNA rearrangement junctions

To understand the conservation of genome rearrangement patterns, we developed a pipeline guided by protein sequence alignment to compare pointer positions for orthologous genes between any two species (Methods, Figure 5A). We compared pointers for 2503 three-species single-copy orthologs. 4448 pointer locations are conserved between Oxytricha and Tetmemena on 1345 ortholog pairs (Supplementary file 9), representing 38.3% of pointers in these orthologs in Oxytricha and 30.9% in Tetmemena. For Oxytricha/E. woodruffi and Tetmemena/E. woodruffi comparisons, 56 and 58 pointer pairs are conserved, respectively. We also identified 23 pointer locations shared among all three species (Supplementary file 9, Figure 5B, Figure 5—source data 1).

Figure 5 with 1 supplement see all

Download asset Open asset

Identification and examples of conserved pointers.

(A) Pipeline for comparison of pointer positions in orthologs. Orthologs are first grouped by OrthoFinder (Emms and Kelly, 2019), and protein sequences of single-copy orthologs aligned by Clustal Omega (Sievers et al., 2011). Then the protein alignments are reverse translated to coding sequence (CDS) alignments by a modified script of pal2nal (105, Methods). Pointers are annotated on the CDS alignments for comparison between any two orthologs. (B) Two examples of pointer conservation across three species. Gray lines represent the alignment of orthologous CDS regions, and boxes show magnified regions containing conserved pointers. The top panel shows a conserved scrambled pointer (*Oxytricha*: Contig889.1.g68; *Tetmemena*: LASU02015390.1.g1; *Euplotes woodruffi*: EUPWOO_MAC_30,105 .g1). The bottom panel shows a conserved nonscrambled pointer (*Oxytricha*: Contig19750.0.g98; *Tetmemena*: LASU02002033.1.g1; *E. woodruffi*: EUPWOO_MAC_31,621 .g1). Pointer sequences are noted, and commas indicate reading frame. Protein domains detected by HMMER (Finn et al., 2011) are marked in purple. (C) Examples of telomere-bearing element (TBE) insertions in nonscrambled internally eliminated sequences. The upper pair of sequences shows an *Oxytricha* TBE pointer (orange insertion of an incomplete TBE2 transposon containing the 42-kD and 57-kD open reading frames) conserved with a *Tetmemena* non-TBE pointer (*Oxytricha*: Contig736.1.g130; *Tetmemena*: LASU02012221.1.g1). Both species have a TA pointer at this junction. The bottom pair of sequences illustrates a case of nonconserved TBE pointers (*Oxytricha*: Contig17579.0.g71; *Tetmemena*: LASU02007616.1.g1).

Figure 5—source data 1 Pointers conserved in all three species.: https://cdn.elifesciences.org/articles/82979/elife-82979-fig5-data1-v2.xlsx
Download elife-82979-fig5-data1-v2.xlsx
Figure 5—source data 2 The telomere-bearing element (TBE) pointers in Oxytricha that are conserved with non-TBE pointers in Tetmemena.: https://cdn.elifesciences.org/articles/82979/elife-82979-fig5-data2-v2.xlsx
Download elife-82979-fig5-data2-v2.xlsx

To test if these pointer locations are genuinely conserved versus coincidental matching by chance, we performed a Monte Carlo simulation, as also used to study intron conservation (Rogozin et al., 2003). We randomly shuffled pointer positions on CDS regions 1000 times and counted the number of conserved pointer pairs expected for each simulation (Methods). Of the 1000 simulations, none exceeded the observed number of conserved pointer pairs between Oxytricha and Tetmemena (p-value <0.001), suggesting evolutionary conservation of pointer positions (Supplementary file 9). A similar result was obtained for pointers conserved in all three species (Supplementary file 9). However, the numbers of pointer pairs conserved between Oxytricha/E. woodruffi and Tetmemena/E. woodruffi is similar to the expectations by chance (Supplementary file 9). The low level of pointer conservation of either hypotrichs with E. woodruffi may reflect the smaller number of IESs in E. woodruffi; hence, most pointers would have arisen in the hypotrich lineage. Furthermore, E. woodruffi is genetically more distant from the two hypotrichs; hence, the accumulation of substitutions would obscure protein sequence homology, which we used to compare pointer locations. For ortholog pairs between Oxytricha and Tetmemena, scrambled pointers are significantly more conserved than nonscrambled ones (chi-square test, p-value <1e−10, Supplementary file 10). We also find that most pointer sequences differ even if the positions are conserved (Figure 5B, Figure 5—source data 1, Supplementary file 11), suggesting that substitutions may accumulate in pointers without substantially altering rearrangement boundaries.

Oxytricha and Tetmemena both contain a high copy number of TBE transposons (Chen et al., 2014; Chen and Landweber, 2016; Supplementary file 3). We investigated the level of TBE conservation between these two species. To identify orthologous insertions, we focus on TBE insertions in nonscrambled IESs on single-copy orthologs, which include 1706 Oxytricha TBEs inserted in 1296 nonscrambled IESs (multiple TBEs can be inserted into an IES) and 180 Tetmemena TBEs inserted into 170 nonscrambled IESs. We refer to the pointer flanking a TBE-containing IES as a TBE pointer. No TBE pointer locations are conserved between two species. This suggests that TBEs might invade the genomes of Oxytricha and Tetmemena independently, or still be actively mobile in the genome. Only 27 Oxytricha TBE pointers (containing 36 TBEs) are conserved with non-TBE pointers in Tetmemena (Figure 5—source data 2, Figure 5C). No Tetmemena TBE pointer is conserved with an Oxytricha non-TBE pointer. This suggests that TBE insertions may preferentially produce new rearrangement junctions instead of inserting into an existing IES.

Intron locations sometimes coincide with DNA rearrangement junctions

Ciliate genomes are generally intron-poor. Oxytricha averages 1.7 introns/gene, Tetmemena has 1.1, and E. woodruffi has 2.2. Among three-species orthologs, intron locations sometimes map near pointer positions (within a 20-bp window, Figure 5B, Figure 5—figure supplement 1). IESs and introns are both noncoding regions that are removed from mature transcripts, though at different stages. A previous single-gene study observed that an IES in Paraurostyla overlaps the position of an intron in Uroleptus, Urostyla, and also the human homolog (Chang et al., 2005). This observation suggested an intron-IES conversion model in which the ability to eliminate non-CDS regions as either DNA or RNA provides a potential backup mechanism. Such interconversion has also been observed between two strains of Stylonychia (Möllenbeck et al., 2006). In the present study, we identified 174 potential cases of intron-IES conversion in the three species (Figure 5—figure supplement 1, Supplementary file 12): 103 (59.2%) E. woodruffi introns map near Oxytricha/Tetmemena pointers. We used a 20-bp window for this analysis, since one would only expect the boundaries of introns and IESs to coincide precisely if they were recent evolutionary conversions. A Monte Carlo simulation for these intron-IES comparisons (Supplementary file 12) revealed that p<0.001 for most three-species comparisons. For two-species comparisons, we identify 306 cases where an intron boundary in one species precisely coincides with a pointer sequence in another species, with strongest statistical support for the comparison between Oxytricha intron positions and Tetmemena IES junctions (p=0.008) (Supplementary file 13). Notably, Tetmemena intron locations rarely coincide with Oxytricha IESs (Supplementary file 13), suggesting a possible bias in the direction of intron-IES conversion during evolution.

The observation that E. woodruffi has the most introns but the smallest number of IESs per gene (Figure 3) is consistent with removal of intragenic non-CDS regions as either DNA or RNA. The intron-sparseness of ciliates is compatible with a hypothesis that it is advantageous to eliminate noncoding regions earlier at the DNA level, with intron deletion sometimes providing an opportunity for repair if they fail to be excised as IESs (Chang et al., 2005).

Evolution of complex genome rearrangements: Russian doll genes

Genome rearrangements in the Oxytricha lineage can include overlapping and nested loci, with MDSs for different MAC loci embedded in each other (Chen et al., 2014; Braun et al., 2018). When multiple gene loci are nested in each other, these have been called Russian doll loci (Braun et al., 2018). Oxytricha contains two loci with five or more layers of nested genes (Braun et al., 2018). Oxytricha and Tetmemena display a high degree of synteny and conservation in both Russian doll loci. In the first Russian doll gene cluster, one nested gene (green) is present in Oxytricha but absent in Tetmemena (Figure 6A, Figure 6—figure supplement 1, Figure 6—figure supplement 2), confirmed by PCR (Methods). Oxytricha also has a complete TBE3 insertion in the green gene (Figure 6A, Figure 6—figure supplement 1A), hinting at a possible link between transposition and new gene insertion. In addition, a two-gene chromosome in Oxytricha (orange) is present as two single-gene chromosomes in Tetmemena (Figure 6A, Figure 6—figure supplement 1). In Oxytricha, seven orange MDSs ligate across two other loci via an 18-bp pointer (TATATCTATACTAAACTT) to form a two-gene nanochromosome. However, in Tetmemena, telomeres are added to the ends of both gene loci instead, forming two independent MAC chromosomes (Figure 6A, Figure 6—figure supplement 1). The second Russian doll locus has an example of a long, conserved pointer (orange dotted line) that bridges three other loci (the green and blue scrambled loci and one nonscrambled locus, Figure 6B). Close to this region is a decayed TBE insertion (769 bp) in Oxytricha. None of the E. woodruffi orthologs of both Russian doll loci maps to the same MIC contig, which suggests that the Russian doll clusters arose after the divergence of Euplotes from the common ancestor of Oxytricha and Tetmemena.

Figure 6 with 2 supplements see all

Download asset Open asset

Synteny in ‘Russian doll’ loci in *Oxytricha* and *Tetmemena*.

(A) Schematic comparison of the Russian doll gene cluster on *Oxytricha* germline micronucleus (MIC) contig OXYTRI_MIC_87484 vs. *Tetmemena* MIC contig TMEMEN_MIC_21461. Boxes of the same color represent clusters of macronuclear destined sequences (MDSs) for orthologous genes (detailed map in Figure 6—figure supplement 1 and Figure 6—figure supplement 2). Numbers in brackets indicate the number of MDSs in each cluster, grouped by somatic macronucleus (MAC) chromosome. One nested gene (green) in *Oxytricha* is absent from *Tetmemena*. A two-gene chromosome (orange) that derives from seven MDSs in *Oxytricha* is processed as two single-gene chromosomes in *Tetmemena* instead (indicated by black border around orange boxes). The purple gene in *Oxytricha* has two paralogs in *Tetmemena*. Black triangles represent conserved, orthologous, and nonscrambled gene loci inserted between nested Russian doll genes. Empty triangle represents scrambled MDSs for other loci. Gray triangles, complete nonscrambled MAC loci embedded between gene layers in one species with no orthologous gene detected in the other species. Black star, a complete telomere-bearing element (TBE) transposon insertion. Gray star, a partial TBE insertion. (B) *Oxytricha* MIC contig OXYTRI_MIC_69233 vs. *Tetmemena* MIC contig TMEMEN_MIC_22886. Pointer sequences bridging the nested MDSs of orange and green genes are highlighted. The underlined pointer portions are conserved between species, e.g., the last 8 bp of the *Oxytricha* pointer, TAAGTTCAAAGTAG, is identical to the first 8 bp of CAAAGTAGCTCAATC in *Tetmemena*, illustrating pointer sliding (DuBois and Prescott, 1995), or gradual shifting of MDS/IES boundaries. White star indicates a decayed TBE with no open reading frame identified.

Discussion

The highly diverse ciliate clade provides a valuable resource for evolutionary studies of genome rearrangement. However, full assembly and annotation of germline MIC genomes have concentrated on the model ciliates Tetrahymena, Paramecium, and Oxytricha. To provide insight into genome evolution in this lineage, we assembled and compared germline and somatic genomes of Tetmemena sp. and an outgroup, E. woodruffi, to that of O. trifallax. This expands our knowledge of the diversity of ciliate genome structures and the evolutionary origin of complex genome rearrangements.

Dramatic variation in transposon copy number (TBE and Tec elements) from the Tc1/mariner family appears to explain most of the variation in MIC genome size. In many eukaryotic taxa, genome size can differ dramatically even for closely related species, a phenomenon known as the ‘C-value paradox’ (Thomas, 1971). Our present observations are compatible with previous reports that the repeat content of the genome, especially transposon content, positively correlates with genome size (Elliott and Gregory, 2015).

Oxytricha has three TBE families in the MIC genome, but only TBE3 is present in Tetmemena, consistent with our previous conclusion that TBE3 is ancestral to the base of the transposon lineage in hypotrichous ciliates (Chen and Landweber, 2016). Tens of thousands of TBE1/2 transposons then expanded specifically in Oxytricha. Despite a high copy number of TBEs in both Oxytricha and Tetmemena, we find no identical TBE locations in nonscrambled IESs, even among syntenic Russian doll regions. These observations suggest that TBEs may be active in these genomes and contribute to the evolution of genome structure.

In the relatively IES-poor genome of E. woodruffi, IESs accumulate upstream of start codons, similar to the 5’ bias of introns in intron-poor organisms (Mourier and Jeffares, 2003). The simplest model to explain 5’ intron bias is homologous recombination between a reverse transcript of an intron-lacking mRNA and the original DNA locus to erase introns in the coding region (Mourier and Jeffares, 2003). A similar mechanism could simultaneously erase IESs in coding regions via germline recombination between the MIC chromosome and a reverse transcript; however, they are usually in different subcellular locations. More plausibly, a source for DNA recombination could be a MAC nanochromosome, since they are already abundant at high copy number, but another source could be by capture of a reverse transcript of a long non-coding template RNA that guides DNA rearrangement (Nowacki et al., 2008; Lindblad et al., 2017). Either recombination event in the germline would lead to loss of IESs, while retaining introns, but neither would necessarily provide a bias for IES-loss in coding regions. Any of these infrequent events would be meaningful on an evolutionary time scale, even if developmentally rare. The 5’ bias of IESs could also reflect an evolutionary bias for continuous coding regions. Alternatively, upstream IESs might regulate gene expression or cell growth (Sellis et al., 2021), like some introns (Parenteau et al., 2019; Morgan et al., 2019).

This study investigated the evolution of scrambled genes by comparing Oxytricha and Tetmemena to E. woodruffi, as an earlier diverged representative of the spirotrich lineage. While E. woodruffi has approximately half as many scrambled genes as Tetmemena and Oxytricha, its genes are also much more continuous. For example, the most scrambled gene in E. woodruffi, encoding a DNA replication licensing factor (EUPWOO_MAC_28518, 3 kb), has only 20 scrambled junctions. The most scrambled gene in Tetmemena (LASU02015934.1, 14.7 kb, encoding a hydrocephalus-inducing-like protein) has 204 scrambled pointers, and the most scrambled gene in Oxytricha (Contig17454.0, 13.7 kb, encoding a dynein heavy chain family protein, Chen et al., 2014) is similarly complex, with 195 scrambled junctions. Together, these observations are consistent with our interpretation that E. woodruffi reflects an evolutionary intermediate stage, as it contains both fewer scrambled loci and fewer scrambled junctions within its scrambled loci. The observation that the most scrambled locus differs in each species is also consistent with the conclusion that complex gene architectures may continue to elaborate independently.

We observed that scrambled genes in each species tend to have more paralogs than nonscrambled genes. Similarly, in C. uncinata (Maurer-Alcalá et al., 2018a), a distantly related ciliate in the class Phyllopharyngea that also has scrambled genes, scrambled gene families (orthogroups) contain more genes (~2.9) than nonscrambled gene families (~1.3) (Maurer-Alcalá et al., 2018a). Apart from duplications at the gene level, E. woodruffi often contains partial MDS duplications at scrambled junctions, annotated as unusually long ‘pointers’ (Figure 4—figure supplement 1). We also demonstrate that odd-even scrambled patterns could readily arise from local duplications (Figure 4—figure supplement 2). These observations are most consistent with a simple model (Gao et al., 2015) in which local duplications permit combinatorial DNA recombination between paralogous germline regions, and mutation accumulation in either paralogs establishes an odd-even scrambled pattern that can propagate by weaving together segments from paralogous sources. Other proposed models include Hoffman and Prescott, 1997 IES-invasion model that suggested that pairs of IESs could invade an MDS, and then subsequently recombine with another IES to yield odd-even scrambled regions; however, a previous examination did not find support for this model (Chang et al., 2005). Prescott et al., 1998 also proposed that some odd-even scrambled loci could arise suddenly via reciprocal recombination with loops of A/T-rich DNA, but this does not exploit paralogy, only the high A/T content in the MIC. We previously proposed a gradual model (Chang et al., 2005; Landweber et al., 2000) in which MDS/IES recombination at short AT-rich repeats (precursors to pointers) could generate and propagate odd-even scrambled patterns. While limited comparisons of orthologs favored the stepwise recombination models (Hogan et al., 2001; Chang et al., 2005; Wong and Landweber, 2006; DuBois and Prescott, 1995), none of the earlier models accounted for the widespread existence of partial paralogy, revealed by genome assemblies.

Local duplications provide a buffer against mutations, allowing paralogous MDSs to repair the MAC locus during assembly of odd/even scrambled genes. Therefore, once an odd/even scrambled locus is established, a consequence is that evolution can only proceed in the direction of accumulating more scrambled junctions, as each new mutation in one paralog necessitates repair via incorporation of the other paralog (Figure 4—figure supplement 1B). This shortens the length of the respective MDSs and increases the number of recombination junctions, creating an evolutionary ratchet that drives the increase in scrambling. The lack of the presence of an error-free, continuous version of this locus in the germline reduces the possibility of losing the scrambled pattern from the MIC genome, relative to the trend toward decreasing MDS lengths as more mutations accumulate in either paralogs, with a resulting increase in the levels of scrambling and fragmentation (Landweber, 2007; Speijer, 2008). The only opportunity to repair a scrambled locus in the MIC would be a rare event that replaces the locus via recombination with a continuous version from the parental MAC, with the source being either parental MAC DNA or a reverse transcript of a template RNA (Nowacki et al., 2008; Lindblad et al., 2017), as discussed above.

Recent exciting reports have also described scrambled genomes in metazoa, including cephalopods (Schmidbaur et al., 2022; Albertin et al., 2022), but those events entail primarily evolutionary shuffling of gene order, without accompanying genome editing or repair. The ciliate lineage is remarkable in having evolved a sophisticated mechanism of RNA-guided genome editing that allows accurate and precise DNA repair of translocations and inversions. The future opportunity to harness this system to develop novel tools for genome editing outside of Oxytricha offers exciting directions.

Methods

DNA collection and sequencing of Tetmemena sp.

Tetmemena sp. (strain SeJ-2015; Chen et al., 2015) was isolated as a single cell from a stock culture and propagated as a clonal strain via vegetative (asexual) cell culture. Cells were cultured in Pringsheim media (0.11 mM Na₂HPO₄, 0.08 mM MgSO₄, 0.85 mM Ca(NO₃)₂, 0.35 mM KCl, pH 7.0) and fed with Chlamydomonas reinhardtii, together with 0.1%(v/v) of an overnight culture of non-virulent Klebsiella pneumoniae. Macronuclei and micronuclei were isolated using sucrose gradient centrifugation (Lauth et al., 1976). Genomic DNA was subsequently purified using the Nucleospin Tissue Kit (Takara Bio USA, Inc). Macronuclear DNA was sequenced and assembled in Chen et al., 2015. Micronuclear DNA was further size-selected via BluePippin (Sage Science) for PacBio sequencing, or via 0.6% (w/v) SeaKem Gold agarose electrophoresis (Lonza) for Illumina sequencing. Micronuclear DNA purification and sequencing protocols are described in Chen et al., 2014.

DNA collection and sequencing for E. woodruffi

E. woodruffi (strain Iz01) was cultured in Volvic water at room temperature and fed with green algae every 2–3 days. We fed cells with C. reinhardtii for MAC DNA collection, and switched to Chlorogonium capillatum for MIC DNA collection. In order to remove algal contamination, cells were starved for at least 2–3 days before collection. Cells were washed and concentrated as in Chen et al., 2014. Because MAC DNA is predominant in whole cell DNA, we used whole cell DNA (purified via NucleoSpin Tissue kit, Takara Bio USA, Inc) for MAC genome sequencing. Paired-end sequencing was performed on an Illumina Hiseq2000 at the Princeton University Genomics Core Facility.

MIC DNA was enriched from whole cell DNA and sequenced via three sequencing platforms (Illumina, Pacific Biosciences, and Oxford Nanopore Technologies). We used conventional and pulse-field gel electrophoresis (PFGE) to enrich MIC DNA:

High-molecular-weight DNA was separated from whole cell DNA by gel-electrophoresis (0.25% agarose gel at 4°C, 120 V for 4 hr). The top band was cut from the gel and purified with the QIAGEN QIAquick kit. The purified high-molecular-weight DNA was directly sent to the group of Dr. Robert Sebra at the Icahn School of Medicine at Mount Sinai for library construction and sequencing. BluePippin (Sage Science) separation was used before sequencing to select DNA >10 kb. DNA was sequenced on two platforms: Illumina HiSeq2500 (150 bp paired-end reads) and PacBio Sequel (SMRT reads).
High-molecular-weight DNA was also enriched by PFGE. E. woodruffi cells were mixed with 1% low-melt agarose to form plugs according to Akematsu et al., 2017, with addition of 1 hr incubation with 50 μg/ml RNase (Invitrogen AM2288) in 10 mM Tris-HCl (pH7.5) at 37°C for RNA depletion. After three washes of 1 hr with 1× TE buffer, the DNA plugs were incubated in 1 mM phenylmethylsulfonyl fluoride (PMSF) to inactivate proteinase K, followed by MspJI (New England Biolabs) digestion at ^mCNNR(9/13) sites to remove contaminant DNA (^mC indicates C5-methylation or C5- hydroxymethylation). Previous reports have shown that no methylcytosine is detectable in vegetative cells of Oxytricha (Bracht et al., 2012), Tetrahymena (Gorovsky et al., 1973), and Paramecium (Cummings et al., 1974), suggesting that C5-methylation and C5-hydroxymethylation are rarely involved in the vegetative growth of the ciliate lineage. We also validated by qPCR that the quantity of two randomly selected MIC loci is not changed after the MspJI digestion. On the contrary, algal genomic DNA is significantly digested by MspJI. Based on these results, we conclude that MspJI digestion can be used to remove bacterial and algal DNA with C5-methylation and C5-hydroxymethylation, leaving E. woodruffi MIC DNA intact. The agarose plugs containing digested DNA were then inserted into wells of 1.0% Certified Megabase agarose gel (Bio-Rad) for PFGE (CHEF-DR II System, Bio-Rad). The DNA was separated at 6 V, 14°C with 0.5× TBE buffer at a 120° angle for 24 hr with switch time of 60–120 s. We validated by qPCR that the E. woodruffi MIC chromosomes were not mobilized from the well, while the MAC DNA migrated into the gel. The MIC DNA was then extracted by phenol-chloroform purification. Library preparation and sequencing were performed at Oxford Nanopore Technologies (New York, NY).

MAC genome assembly of E. woodruffi

We assembled the MAC genome of E. woodruffi using the same pipeline for Tetmemena sp. (Chen et al., 2015) for comparative analysis: two draft genomes were assembled by SPAdes (Bankevich et al., 2012) and Trinity (Grabherr et al., 2011), and were then merged by CAP3 (Huang and Madan, 1999). Trinity, which is a software developed for de novo transcriptome assembly (Grabherr et al., 2011), has been used to assemble hypotrich MAC genomes (Chen et al., 2015) because their nanochromosome genome structure is similar to transcriptomes, including properties such as variable copy number and alternative isoforms (Lindblad et al., 2019). Telomeric reads were mapped to contigs by BLAT (Kent, 2002), and contigs were further extended and capped by telomeres when at least five reads pile up at a position near ends by custom python scripts (https://github.com/yifeng-evo/Oxytricha_Tetmemena_Euplotes/tree/main/MAC_genome_telomere_capping) (Feng, 2022a). The mitochondrial DNA was removed if the contig has a TBLASTX (Camacho et al., 2009) hit on the Oxytricha mitochondrial genome (Genbank accession JN383842.1 and JN383843.1) or two Euplotes mitochondrial genomes (Euplotes minuta GQ903130.1, E. crassus GQ903131.1). Algal contigs were removed by BLASTN to all C. reinhardtii nucleotide sequences downloaded from Genbank. Non-telomeric contigs were mapped to bacterial NR by BLASTX to remove bacterial contaminations. The genome was further compressed by CD-HIT (Fu et al., 2012) in two steps: (1) contigs <500 bp were removed if 90% of the short contig can be aligned to a contig ≥ 500 bp with 90% similarity (-c 0.9 -aS 0.9 -uS 0.1); (2) then the genome was compressed by 95% similarity (-c 0.95 -aS 0.9 -uS 0.1). Contigs shorter than 500 bp without telomeres were removed. Nine contigs, likely Tec contaminants from the MIC genome, were also excluded (Tblastn, ‘-db_gencode 10 -evalue 1e-5’), and they could be assembled due to the high copy number in the MIC genome (47, 48, Genbank accessions of Tec ORFs are AAA62601.1, AAA62602.1, AAA62603.1, AAA91339.1, AAA91340.1, AAA91341.1, AAA91342.1).

RNA sequencing of E. woodruffi and Tetmemena sp.

Three biological replicates of total RNA was isolated from asexually growing E. woodruffi and Tetmemena sp. cells using TRIzol reagent (Thermo Fisher Scientific) and enriched for the poly(A)+fraction using the NEBNext Poly(A) mRNA Magnetic Isolation Module (New England Biolabs). Stranded RNA-seq libraries were constructed using the ScriptSeq v2 RNA-seq library preparation kit (Epicentre) and sequenced on an Illumina Nextseq500 at the Columbia Genome Center. For E. woodruffi, the transcriptome was assembled by Trinity (Grabherr et al., 2011), and transcript alignments to the MAC genome were generated by PASA (Haas et al., 2003).

Gene prediction of the E. woodruffi MAC genome and validation of MAC genome completeness

We followed the gene prediction pipeline developed by the Broad institute (https://github.com/PASApipeline/PASApipeline/wiki); using EVidenceModeler (EVM, Haas et al., 2008) to generate the final gene predictions. EVM produced gene structures by weighted combination of evidence from three resources: ab initio prediction, protein alignments, and transcript alignments (the weight was 3, 3, and 10 respectively). Ab initio prediction was generated by BRAKER2 pipeline (Brůna et al., 2021). Protein alignments for EVM were generated by mapping Oxytricha proteins to the E. woodruffi MAC genome by Exonerate (Slater and Birney, 2005). EVM predicted 33,379 genes on MAC chromosomes with at least one telomere.

We assessed MAC genome completeness using three methods: (1) 28,294 (80.6%) of the 35,099 E. woodruffi MAC contigs have at least one telomere. (2) In the E. woodruffi genes predicted on telomeric contigs, 88.8% of BUSCO (Simão et al., 2015; Manni et al., 2021) genes in the lineage database alveolata_odb10 were identified as complete. Within the 171 BUSCO genes, 135 are complete and single-copy, 17 are complete and duplicated, 7 are fragmented, and 12 are missing. This represents the best Euplotes MAC genome assembly available. (3) We identified 51 tRNA genes encoding all 20 amino acids by tRNAscan-SE (Lowe and Eddy, 1997) in the MAC genome, including two suppressor tRNAs of UAA and UAG.

MIC genome assembly of Tetmemena sp.

The MIC genome of Tetmemena was assembled with a hybrid approach to combine reads from different sequencing platforms. Tetmemena Illumina reads were first assembled by SPAdes (77, parameters ‘-k 21,33,55,77,99,127 –careful’). PacBio reads were error corrected by FMLRC (Wang et al., 2018) using Illumina reads with default parameters. Corrected PacBio reads were aligned to both the MAC genome and the Illumina MIC assembly with BLASTN. Reads were removed if they start or end with telomeres or are aligned better to the MAC. The remaining reads were assembled with wtdbg2 (Ruan and Li, 2020, parameters ‘-x rs’). The PacBio assembly was polished by Pilon (Walker et al., 2014) with the ‘--diploid’ option. The Illumina and PacBio assemblies were merged by quickmerge (Chakraborty et al., 2016) with the ‘-l 5000’ option.

MIC genome assembly of E. woodruffi

The MIC genome of E. woodruffi was assembled using a similar procedure as described above for Tetmemena. E. woodruffi reads were filtered to remove bacterial contamination, including abundant high-GC-content contaminants, possibly endosymbionts (Boscaro et al., 2019). Nanopore reads with GC content ≥55% were assembled by Flye (Kolmogorov et al., 2019) with the parameter ‘--meta’ for metagenomic assembly of bacterial contigs. We used kaiju (Menzel et al., 2016) to identify bacteria taxa for these contigs. 9 of 10 top-covered contigs derive from Proteobacteria, from which many Euplotes symbionts derive (Boscaro et al., 2019). Bacterial contamination was removed from Illumina reads if perfectly mapping to these metagenomic contigs by Bowtie2 (Langmead and Salzberg, 2012). The cleaned Illumina reads were then assembled by SPAdes with ‘-k 21,33,55,77,99,127’ (Bankevich et al., 2012). Pacbio raw reads and Nanopore raw reads with GC content <55% were aligned to a concatenated database containing both the MAC genome and the Illumina MIC assembly with BLASTN. Reads were removed if they start or end with telomeres or align better to the MAC. Remaining PacBio/Nanopore reads were assembled by Flye with ‘--meta’ mode. The PacBio-Nanopore assembly was polished by Pilon with the ‘--diploid’ option. Illumina and PacBio-Nanopore assemblies were merged by quickmerge with the ‘-l 10000’ option. Contigs shorter than 1 kb were removed.

MIC genome decontamination

The draft MIC genome of Tetmemena was first mapped to telomeric MAC contigs by BLASTN. MIC contigs containing MDSs were included in the final assembly. The rest of the MIC contigs were filtered by a decontamination pipeline: (1) contigs were aligned to the K. pneumoniae genome, C. reinhardtii genome, and the Oxytricha mitochondrial genome by BLASTN to remove contaminants; (2) the remaining contigs were then searched against the bacteria NR database and a ciliate protein database (including protein sequences annotated in Tetrahymena thermophila: http://www.ciliate.org/system/downloads/tet-latest/4-Protein%20fasta.fasta; Paramecium tetraurelia: http://paramecium.cgm.cnrs-gif.fr; and O. trifallax: https://oxy.ciliate.org) by BLASTX. Contigs with higher bit score to bacteria NR or G+C >45% were removed. The E. woodruffi MIC genome was decontaminated, similarly, with addition of all Chlorogonium sequences (the algal food source) on NCBI and the two Euplotes mitochondrial genomes (E. minuta GQ903130.1, E. crassus GQ903131.1) to filter contaminants.

Repeat identification

The repeat content in the MIC genomes was identified by RepeatModeler 1.0.10 (Smit and Hubley, 2008) and RepeatMasker 4.0.7 (Smit et al., 2013) with default parameters.

TBE/Tec detection

Representative Oxytricha TBE ORFs (Genbank accession AAB42034.1, AAB42016.1, and AAB42018.1) were used as queries to search TBEs in the Oxytricha and Tetmemena MIC genomes by TBLASTN (-db_gencode 6 -evalue 1e-7 -max_target_seqs 30000). Tec ORFs were similarly detected by using E. crassus Tec1 and Tec2 ORFs as queries (-db_gencode 10 -evalue 1e-5 -max_target_seqs 30000, Genbank accessions of Tec ORFs are AAA62601.1, AAA62602.1, AAA62603.1, AAA91339.1, AAA91340.1, AAA91341.1, AAA91342.1). Complete TBEs/Tecs were determined by custom python scripts when three ORFs are within 2000 bp from each other and in correct orientation (https://github.com/yifeng-evo/Oxytricha_Tetmemena_Euplotes/tree/main/TBE_ORFs/TBE_to_oxy_genome_tblastn_parse.py, Chen and Landweber, 2016). 30 TBE ORFs with >70% completeness were subsampled from each species for phylogenetic analysis (except for the 57 kD ORF in Tetmemena, for which 21 were subsampled). The subsampled TBE ORFs were aligned using MUSCLE (Edgar, 2004), and the alignments were trimmed by trimAl ‘-automated1’ (Capella-Gutiérrez et al., 2009). Phylogenetic trees were constructed using PhyML 3.3 (Guindon et al., 2010).

Rearrangement annotations

SDRAP (Braun et al., 2022) was used to annotate MDSs, pointers, and MIC-specific regions (minimum percent identity for preliminary match annotation = 95, minimum percent identity for additional match annotation = 90, minimum length of pointer annotation = 2). SDRAP requires MAC and MIC genomes as input. For the SDRAP annotation of Oxytricha, we used the MAC genome from Swart et al., 2013 instead of the latest hybrid assembly that incorporated PacBio reads (Lindblad et al., 2019), because the former version was primarily based on Illumina reads, similar to the MAC genomes of Tetmemena (7, Genbank GCA_001273295.2) and E. woodruffi which are also Illumina assemblies. Oxytricha and Tetmemena MAC genomes were preprocessed by removing MAC contigs with TBE ORFs, considered MIC contaminants (Chen and Landweber, 2016). SDRAP is a new program that can output the rearrangement annotations with minor differences from Chen et al., 2014, but most annotations are robust (Figure 3—figure supplement 2). Scrambled and nonscrambled junctions/IESs were annotated by custom python scripts (https://github.com/yifeng-evo/Oxytricha_Tetmemena_Euplotes/tree/main/scrambled_nonscrambled_IES_pointer).

MIC genome categories

Each MIC genome region is assigned to only one category in Figure 2A–C, even if it belongs to more than one category. The assignment is based on the following priority: MDS, TBE/Tec, MIC genes (only available for Oxytricha, which has developmental RNA-seq data), IES, tandem repeats, other repeats, and non-coding non-repetitive regions. For example, an MIC region can be a TBE in an IES, and it is only considered as TBE in Figure 2A–C.

Ortholog comparison pipeline and Monte Carlo simulations

Orthogroups of genes on telomeric MAC contigs were detected by OrthoFinder with ‘-S blast’ (Emms and Kelly, 2019). Single-copy orthologs were aligned by Clustal Omega (Sievers et al., 2011). Protein alignments were reversely translated to CDS alignments by a modified script of pal2nal (Suyama et al., 2006, https://github.com/yifeng-evo/Oxytricha_Tetmemena_Euplotes/tree/main/Ortholog_comparison/pal2nal.pl). Two modifications were made in the script: (1) the modified script allows pal2nal to take different genetic codes for three sequences (-codontable 6,6,10); (2) the script also fixed an error in the original pal2nal script in which codontable 10 for the Euplotid nuclear code was the same as the universal code. Visualization of pointer positions and intron locations on orthologs was implemented by a custom python script (https://github.com/yifeng-evo/Oxytricha_Tetmemena_Euplotes/blob/main/Ortholog_comparison/visualization_of_ortholog_comparison.py). Pointer positions or intron locations are considered conserved if they are within a 20-bp alignment window on the CDS alignment. Protein domains were annotated by HMMER (Finn et al., 2011). We performed Monte Carlo simulations by randomly shuffling pointer locations on the CDS but keeping their original position distribution. This was implemented by a custom python script, which transforms the CDS to a circle, rotates pointer positions on the circle, and outputs the shuffled position on the re-linearized CDS (https://github.com/yifeng-evo/Oxytricha_Tetmemena_Euplotes/blob/main/Ortholog_comparison/shuffle_simulation.py). The null hypothesis of the Monte Carlo test is that pointer positions are conserved by chance. p-Value of Monte Carlo test is given by N_{expected>observed}/N_total (N_{expected>observed} is the number of simulations when there are more conserved pointers in the simulation than the observation from real data, N_total = 1000 in this study).

PCR validation of Russian doll locus

The complex Russian doll locus on MIC contig TMEMEN_MIC_21461 in Tetmemena was validated by PCR to confirm the Tetmemena MIC genome assembly. Tetmemena micronuclear DNA was purified as described previously and used as template for PCR using PrimeSTAR Max DNA polymerase (Takara Bio). 11 primer sets (Supplementary file 14) were designed to amplify products between 3 kb and 6 kb in length, with overlapping regions between consecutive primer pairs. The resulting PCR products were visualized through agarose gel electrophoresis, and bands of the expected size were extracted using a Monarch DNA Gel Extraction Kit (New England Biolabs). The purified gel bands were cloned using a TOPO XL-2 Complete PCR Cloning Kit (Invitrogen), transformed into One Shot OmniMAX 2 T1R E. coli cells (Invitrogen), and individual clones were grown and their plasmids harvested with a QIAprep Spin Miniprep Kit (QIAGEN). The plasmid ends were Sanger sequenced, as well as the region where the Oxytricha MIC assembly contains inserted MDSs (Genewiz). Sanger sequencing reads were mapped to the Tetmemena MIC contig TMEMEN_MIC_21461 and visualized using Geneious Prime 2021.1.1 (https://www.geneious.com).

Availability of data and materials

Custom scripts are public on https://github.com/yifeng-evo/Oxytricha_Tetmemena_Euplotes, (Feng, 2022b copy archived at swh:1:rev:fd66a0efeaf9feb2d79e183313192d641b4e5400). DNA-seq reads and genome assemblies are available at GenBank under Bioprojects PRJNA694964 (Tetmemena sp.) and PRJNA781979 (E. woodruffi). Genbank accession numbers for genomes are JAJKFJ000000000 (Tetmemena sp. Micronucleus genome), JAJLLS000000000 (E. woodruffi Micronucleus genome), and JAJLLT000000000 (E. woodruffi Macronucleus genome).

Three replicates of RNA-seq reads for vegetative cells are available at GenBank under accession numbers of SRR21815378, SRR21815379, and SRR21815380 for E. woodruffi and SRR21817702, SRR21817703, and SRR21817704 for Tetmemena sp.

MDS annotations for three species are available at https://doi.org/10.5061/dryad.5dv41ns96 and https://knot.math.usf.edu/mds_ies_db/2022/downloads.html (please select species from the drop-down menu).

Data availability

Custom scripts are public on https://github.com/yifeng-evo/Oxytricha_Tetmemena_Euplotes, (copy archived at swh:1:rev:fd66a0efeaf9feb2d79e183313192d641b4e5400). DNA-seq reads and genome assemblies are available at GenBank under Bioprojects PRJNA694964 (Tetmemena sp.) and PRJNA781979 (Euplotes woodruffi). Genbank accession numbers for genomes are JAJKFJ000000000 (Tetmemena sp. Micronucleus genome), JAJLLS000000000 (Euplotes woodruffi Micronucleus genome), and JAJLLT000000000 (Euplotes woodruffi Macronucleus genome). Three replicates of RNA-seq reads for vegetative cells are available at GenBank under accession numbers of SRR21815378, SRR21815379, SRR21815380 for E. woodruffi and SRR21817702, SRR21817703 and SRR21817704 for Tetmemena sp. MDS annotations for three species are available at https://doi.org/10.5061/dryad.5dv41ns96 and https://knot.math.usf.edu/mds_ies_db/2022/downloads.html (please select species from the drop-down menu).

The following data sets were generated

1. Feng Y
2. Neme R
3. Beh LY
4. Chen X
5. Braun J
6. Lu MW
7. Landweber LF
(2022) NCBI BioProject
ID PRJNA781979. Euplotes woodruffi genome sequencing and assembly.

https://www.ncbi.nlm.nih.gov/bioproject/PRJNA781979
1. Feng Y
2. Neme R
3. Beh LY
4. Chen X
5. Braun J
6. Lu MW
7. Landweber LF
(2022) NCBI BioProject
ID PRJNA694964. Tetmemena sp. micronucleus genome sequencing and assembly.

https://www.ncbi.nlm.nih.gov/bioproject/PRJNA694964
1. Feng Y
2. Neme R
3. Beh LY
4. Chen X
5. Braun J
6. Lu MW
7. Landweber LF
(2022) NCBI BioProject
ID PRJNA781602. Euplotes woodruffi strain:Iz01.

https://www.ncbi.nlm.nih.gov/bioproject/PRJNA781602
1. Feng Y
2. Neme R
3. Beh LY
4. Chen X
5. Braun J
6. Lu MW
7. Landweber LF
(2022) NCBI BioProject
ID PRJNA887426. RNA-seq of Tetmemena sp.

https://www.ncbi.nlm.nih.gov/bioproject/PRJNA887426
(2022) MDSIESDB
ID db/2022/downloads. MDS-IES database.

https://knot.math.usf.edu/mds_ies_db/2022/downloads.html
1. Feng Y
2. Neme R
3. Beh L
4. Chen X
5. Braun J
6. Lu M
7. Landweber L
(2022) Dryad Digital Repository
MDS and IES annotations for Euplotes woodruff, Tetmemena sp. and Oxytricha trifallax.

https://doi.org/10.5061/dryad.5dv41ns96
1. Feng Y
2. Neme R
3. Beh LY
4. Chen X
5. Braun J
6. MW Lu
7. Landweber LF
(2022) NCBI GenBank
ID JAJKFJ000000000. Tetmemena sp. Micronucleus genome.

https://www.ncbi.nlm.nih.gov/nuccore/JAJKFJ000000000
1. Feng Y
2. Neme R
3. Beh LY
4. Chen X
5. Braun J
6. MW Lu
7. Landweber LF
(2022) NCBI GenBank
ID JAJLLS000000000. Euplotes woodruffi Micronucleus genome.

https://www.ncbi.nlm.nih.gov/nuccore/JAJLLS000000000
1. Feng Y
2. Neme R
3. Beh LY
4. Chen X
5. Braun J
6. MW Lu
7. Landweber LF
(2022)
NCBI GenBank

Euplotes woodruffi Macronucleus genome.

The following previously published data sets were used

1. Chen X
2. Bracht JR
3. Goldman AD
4. Dolzhenko E
5. Clay DM
6. Swart EC
7. Perlman DH
8. Doak TG
9. Stuart A
10. Amemiya CT
11. Sebra RP
12. Landweber LF
(2014) NCBI Assembly
ID GCA_000711775.1. Oxytricha trifallax micronucleus genome.

https://www.ncbi.nlm.nih.gov/assembly/GCA_000711775.1
1. Swart EC
2. Bracht JR
3. Magrini V
4. Minx P
5. Chen X
6. Zhou Y
7. Khurana JS
8. Goldman AD
9. Nowacki M
10. Schotanus K
11. Jung S
12. Ly A
13. McGrath S
14. Haub K
15. Wiggins JL
16. Storton D
17. Matese JC
18. Parsons L
19. Chang WJ
20. Bowen MS
21. Stover NA
22. Jones TA
23. Eddy SR
24. Herrick GA
25. Doak TG
26. Wilson RK
27. Mardis ER
28. Landweber LF
(2013) NCBI Assembly
ID GCA_000295675.1. Oxytricha trifallax macronucleus genome.

https://www.ncbi.nlm.nih.gov/assembly/GCA_000295675.1/
1. Chen X
2. Jung S
3. Beh LY
4. Eddy SR
5. Landweber LF
(2015) NCBI Assembly
ID GCA_001273295.2. Tetmemena sp. macronucleus genome.

https://www.ncbi.nlm.nih.gov/assembly/GCA_001273295.2
1. Beh LY
2. Debelouchina GT
3. Clay DM
4. Thompson RE
5. Lindblad KA
6. Hutton ER
7. Bracht JR
8. Sebra RP
9. Muir TW
10. Landweber LF
(2019) NCBI Gene Expression Omnibus
ID GSE94421. Genome-wide analysis of chromatin and transcription in the ciliates Oxytricha trifallax and Tetrahymena thermophila.

https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE94421

References

1. Akematsu T
2. Fukuda Y
3. Garg J
4. Fillingham JS
5. Pearlman RE
6. Loidl J
(2017) Post-meiotic DNA double-strand breaks occur in Tetrahymena, and require topoisomerase II and SPO11
eLife 6:e26176.

https://doi.org/10.7554/eLife.26176
- PubMed
- Google Scholar
1. Albertin CB
2. Medina-Ruiz S
3. Mitros T
4. Schmidbaur H
5. Sanchez G
6. Wang ZY
7. Grimwood J
8. Rosenthal JJC
9. Ragsdale CW
10. Simakov O
11. Rokhsar DS
(2022) Genome and transcriptome mechanisms driving cephalopod evolution
Nature Communications 13:2427.

https://doi.org/10.1038/s41467-022-29748-w
- PubMed
- Google Scholar
1. Arnaiz O
2. Mathy N
3. Baudry C
4. Malinsky S
5. Aury JM
6. Denby Wilkes C
7. Garnier O
8. Labadie K
9. Lauderdale BE
10. Le Mouël A
11. Marmignon A
12. Nowacki M
13. Poulain J
14. Prajer M
15. Wincker P
16. Meyer E
17. Duharcourt S
18. Duret L
19. Bétermier M
20. Sperling L
(2012) The Paramecium germline genome provides a niche for intragenic parasitic DNA: evolutionary dynamics of internal eliminated sequences
PLOS Genetics 8:e1002984.

https://doi.org/10.1371/journal.pgen.1002984
- PubMed
- Google Scholar
1. Aury JM
2. Jaillon O
3. Duret L
4. Noel B
5. Jubin C
6. Porcel BM
7. Ségurens B
8. Daubin V
9. Anthouard V
10. Aiach N
11. Arnaiz O
12. Billaut A
13. Beisson J
14. Blanc I
15. Bouhouche K
16. Câmara F
17. Duharcourt S
18. Guigo R
19. Gogendeau D
20. Katinka M
21. Keller AM
22. Kissmehl R
23. Klotz C
24. Koll F
25. Le Mouël A
26. Lepère G
27. Malinsky S
28. Nowacki M
29. Nowak JK
30. Plattner H
31. Poulain J
32. Ruiz F
33. Serrano V
34. Zagulski M
35. Dessen P
36. Bétermier M
37. Weissenbach J
38. Scarpelli C
39. Schächter V
40. Sperling L
41. Meyer E
42. Cohen J
43. Wincker P
(2006) Global trends of whole-genome duplications revealed by the ciliate Paramecium tetraurelia
Nature 444:171–178.

https://doi.org/10.1038/nature05230
- PubMed
- Google Scholar
(1989) Micronuclear genome organization in Euplotes crassus: a transposonlike element is removed during macronuclear development
Molecular and Cellular Biology 9:3793–3807.

https://doi.org/10.1128/mcb.9.9.3793-3807.1989
- PubMed
- Google Scholar
1. Bankevich A
2. Nurk S
3. Antipov D
4. Gurevich AA
5. Dvorkin M
6. Kulikov AS
7. Lesin VM
8. Nikolenko SI
9. Pham S
10. Prjibelski AD
11. Pyshkin AV
12. Sirotkin AV
13. Vyahhi N
14. Tesler G
15. Alekseyev MA
16. Pevzner PA
(2012) SPAdes: a new genome assembly algorithm and its applications to single-cell sequencing
Journal of Computational Biology 19:455–477.

https://doi.org/10.1089/cmb.2012.0021
- PubMed
- Google Scholar
1. Beh LY
2. Debelouchina GT
3. Clay DM
4. Thompson RE
5. Lindblad KA
6. Hutton ER
7. Bracht JR
8. Sebra RP
9. Muir TW
10. Landweber LF
(2019) Identification of a DNA N6-adenine methyltransferase complex and its impact on chromatin organization
Cell 177:1781–1796.

https://doi.org/10.1016/j.cell.2019.04.028
- PubMed
- Google Scholar
(2018) Discovery of the first germline-restricted gene by subtractive transcriptomic analysis in the zebra finch, Taeniopygia guttata
Current Biology 28:1620–1627.

https://doi.org/10.1016/j.cub.2018.03.067
- PubMed
- Google Scholar
(2019) Symbionts of the ciliate Euplotes: diversity, patterns and potential as models for bacteria-eukaryote endosymbioses
Proceedings of the Royal Society B: Biological Sciences 286:20190693.

https://doi.org/10.1098/rspb.2019.0693
- PubMed
- Google Scholar
(2012) Cytosine methylation and hydroxymethylation mark DNA for elimination in Oxytricha trifallax
Genome Biology 13:1–23.

https://doi.org/10.1186/gb-2012-13-10-r99
- PubMed
- Google Scholar
1. Bracht JR
2. Fang W
3. Goldman AD
4. Dolzhenko E
5. Stein EM
6. Landweber LF
(2013) Genomes on the edge: programmed genome instability in ciliates
Cell 152:406–416.

https://doi.org/10.1016/j.cell.2013.01.005
- PubMed
- Google Scholar
1. Braun J
2. Nabergall L
3. Neme R
4. Landweber LF
5. Saito M
6. Jonoska N
(2018) Russian doll genes and complex chromosome rearrangements in Oxytricha trifallax
G3: Genes, Genomes, Genetics 8:1669–1674.

https://doi.org/10.1534/g3.118.200176
- PubMed
- Google Scholar
Preprint
1. Braun J
2. Neme R
3. Feng Y
4. Landweber LF
5. Jonoska N
(2022) SDRAP for Annotating Scrambled or Rearranged Genomes
bioRxiv.

https://doi.org/10.1101/2022.10.24.513505
- Google Scholar
1. Brůna T
2. Hoff KJ
3. Lomsadze A
4. Stanke M
5. Borodovsky M
(2021) BRAKER2: automatic eukaryotic genome annotation with genemark-EP+ and AUGUSTUS supported by a protein database
NAR Genomics and Bioinformatics 3:lqaa108.

https://doi.org/10.1093/nargab/lqaa108
- PubMed
- Google Scholar
1. Burns J
2. Kukushkin D
3. Chen X
4. Landweber LF
5. Saito M
6. Jonoska N
(2016) Recurring patterns among scrambled genes in the encrypted genome of the ciliate Oxytricha trifallax
Journal of Theoretical Biology 410:171–180.

https://doi.org/10.1016/j.jtbi.2016.08.038
- PubMed
- Google Scholar
1. Camacho C
2. Coulouris G
3. Avagyan V
4. Ma N
5. Papadopoulos J
6. Bealer K
7. Madden TL
(2009) BLAST+: architecture and applications
BMC Bioinformatics 10:1–9.

https://doi.org/10.1186/1471-2105-10-421
- PubMed
- Google Scholar
(2009) TrimAl: a tool for automated alignment trimming in large-scale phylogenetic analyses
Bioinformatics 25:1972–1973.

https://doi.org/10.1093/bioinformatics/btp348
- PubMed
- Google Scholar
(2016) Contiguous and accurate de novo assembly of metazoan genomes with modest long read coverage
Nucleic Acids Research 44:e147.

https://doi.org/10.1093/nar/gkw654
- PubMed
- Google Scholar
1. Chang WJ
2. Bryson PD
3. Liang H
4. Shin MK
5. Landweber LF
(2005) The evolutionary origin of a complex scrambled gene
PNAS 102:15149–15154.

https://doi.org/10.1073/pnas.0507682102
- PubMed
- Google Scholar
1. Chen X
2. Bracht JR
3. Goldman AD
4. Dolzhenko E
5. Clay DM
6. Swart EC
7. Perlman DH
8. Doak TG
9. Stuart A
10. Amemiya CT
11. Sebra RP
12. Landweber LF
(2014) The architecture of a scrambled genome reveals massive levels of genomic rearrangement during development
Cell 158:1187–1198.

https://doi.org/10.1016/j.cell.2014.07.034
- PubMed
- Google Scholar
1. Chen X
2. Jung S
3. Beh LY
4. Eddy SR
5. Landweber LF
(2015) Combinatorial DNA rearrangement facilitates the origin of new genes in ciliates
Genome Biology and Evolution 7:2859–2870.

https://doi.org/10.1093/gbe/evv172
- PubMed
- Google Scholar
1. Chen X
2. Landweber LF
(2016) Phylogenomic analysis reveals genome-wide purifying selection on TBE transposons in the ciliate Oxytricha
Mobile DNA 7:2.

https://doi.org/10.1186/s13100-016-0057-9
- PubMed
- Google Scholar
1. Chen X
2. Jiang Y
3. Gao F
4. Zheng W
5. Krock TJ
6. Stover NA
7. Lu C
8. Katz LA
9. Song W
(2019) Genome analyses of the new model protist Euplotes vannus focusing on genome rearrangement and resistance to environmental stressors
Molecular Ecology Resources 19:1292–1308.

https://doi.org/10.1111/1755-0998.13023
- PubMed
- Google Scholar
1. Chen W
2. Zuo C
3. Wang C
4. Zhang T
5. Lyu L
6. Qiao Y
7. Zhao F
8. Miao M
(2021) The hidden genomic diversity of ciliated protists revealed by single-cell genome sequencing
BMC Biology 19:264.

https://doi.org/10.1186/s12915-021-01202-1
- PubMed
- Google Scholar
(1974) Methylated bases in DNA from Paramecium aurelia
Biochimica et Biophysica Acta - Nucleic Acids and Protein Synthesis 374:1–11.

https://doi.org/10.1016/0005-2787(74)90194-4
- Google Scholar
(2016) ParTIES: a toolbox for Paramecium interspersed DNA elimination studies
Bioinformatics 32:599–601.

https://doi.org/10.1093/bioinformatics/btv691
- PubMed
- Google Scholar
(2003) Selection on the genes of Euplotes crassus Tec1 and Tec2 transposons: evolutionary appearance of a programmed frameshift in a Tec2 gene encoding a tyrosine family site-specific recombinase
Eukaryotic Cell 2:95–102.

https://doi.org/10.1128/EC.2.1.95-102.2003
- PubMed
- Google Scholar
1. DuBois MI
2. Prescott DM
(1995) Scrambling of the actin I gene in two Oxytricha species
PNAS 92:3888–3892.

https://doi.org/10.1073/pnas.92.9.3888
- PubMed
- Google Scholar
1. Edgar RC
(2004) MUSCLE: multiple sequence alignment with high accuracy and high throughput
Nucleic Acids Research 32:1792–1797.

https://doi.org/10.1093/nar/gkh340
- PubMed
- Google Scholar
1. Eisen JA
2. Coyne RS
3. Wu M
4. Wu D
5. Thiagarajan M
6. Wortman JR
7. Badger JH
8. Ren Q
9. Amedeo P
10. Jones KM
11. Tallon LJ
12. Delcher AL
13. Salzberg SL
14. Silva JC
15. Haas BJ
16. Majoros WH
17. Farzad M
18. Carlton JM
19. Smith RK
20. Garg J
21. Pearlman RE
22. Karrer KM
23. Sun L
24. Manning G
25. Elde NC
26. Turkewitz AP
27. Asai DJ
28. Wilkes DE
29. Wang Y
30. Cai H
31. Collins K
32. Stewart BA
33. Lee SR
34. Wilamowska K
35. Weinberg Z
36. Ruzzo WL
37. Wloga D
38. Gaertig J
39. Frankel J
40. Tsao CC
41. Gorovsky MA
42. Keeling PJ
43. Waller RF
44. Patron NJ
45. Cherry JM
46. Stover NA
47. Krieger CJ
48. del Toro C
49. Ryder HF
50. Williamson SC
51. Barbeau RA
52. Hamilton EP
53. Orias E
(2006) Macronuclear genome sequence of the ciliate Tetrahymena thermophila, a model eukaryote
PLOS Biology 4:e286.

https://doi.org/10.1371/journal.pbio.0040286
- PubMed
- Google Scholar
1. Elliott TA
2. Gregory TR
(2015) What’s in a genome? the C-value enigma and the evolution of eukaryotic genome content
Philosophical Transactions of the Royal Society of London. Series B, Biological Sciences 370:20140331.

https://doi.org/10.1098/rstb.2014.0331
- PubMed
- Google Scholar
1. Emms DM
2. Kelly S
(2019) OrthoFinder: phylogenetic orthology inference for comparative genomics
Genome Biology 20:238.

https://doi.org/10.1186/s13059-019-1832-y
- PubMed
- Google Scholar
1. Feng Y
2. Beh LY
3. Chang WJ
4. Landweber LF
(2020) SIGAR: inferring features of genome architecture and DNA rearrangements by split-read mapping
Genome Biology and Evolution 12:1711–1718.

https://doi.org/10.1093/gbe/evaa147
- PubMed
- Google Scholar
1. Feng Y
2. Landweber LF
(2021) Transposon debris in ciliate genomes
PLOS Biology 19:e3001354.

https://doi.org/10.1371/journal.pbio.3001354
- PubMed
- Google Scholar
Software
1. Feng Y
(2022a) MAC genome telomere capping script, version 871eb00
Github.

https://github.com/yifeng-evo/Oxytricha_Tetmemena_Euplotes/tree/main/MAC_genome_telomere_capping
Software
1. Feng Y
(2022b) Oxytricha_Tetmemena_Euplotes, version swh:1:rev:fd66a0efeaf9feb2d79e183313192d641b4e5400
Software Heritage.

https://archive.softwareheritage.org/swh:1:dir:8ad132d58c3073da701bdde6700a37e2cdc01509;origin=https://github.com/yifeng-evo/Oxytricha_Tetmemena_Euplotes;visit=swh:1:snp:3e53ca9f9f0b0bc48a5c56d379e0def68cce596f;anchor=swh:1:rev:fd66a0efeaf9feb2d79e183313192d641b4e5400
(2011) HMMER web server: interactive sequence similarity searching
Nucleic Acids Research 39:W29–W37.

https://doi.org/10.1093/nar/gkr367
- PubMed
- Google Scholar
1. Fu L
2. Niu B
3. Zhu Z
4. Wu S
5. Li W
(2012) CD-HIT: accelerated for clustering the next-generation sequencing data
Bioinformatics 28:3150–3152.

https://doi.org/10.1093/bioinformatics/bts565
- PubMed
- Google Scholar
1. Gao F
2. Song W
3. Katz LA
(2014) Genome structure drives patterns of gene family evolution in ciliates, a case study using Chilodonella uncinata (Protista, Ciliophora, Phyllopharyngea)
Evolution; International Journal of Organic Evolution 68:2287–2295.

https://doi.org/10.1111/evo.12430
- PubMed
- Google Scholar
1. Gao F
2. Roy SW
3. Katz LA
(2015) Analyses of alternatively processed genes in ciliates provide insights into the origins of scrambled genomes and may provide a mechanism for speciation
MBio 6:e01998-14.

https://doi.org/10.1128/mBio.01998-14
- PubMed
- Google Scholar
1. Gao F
2. Warren A
3. Zhang Q
4. Gong J
5. Miao M
6. Sun P
7. Xu D
8. Huang J
9. Yi Z
10. Song W
(2016) The all-data-based evolutionary hypothesis of ciliated protists with a revised classification of the phylum Ciliophora (eukaryota, alveolata)
Scientific Reports 6:24874.

https://doi.org/10.1038/srep24874
- PubMed
- Google Scholar
(1973) (6 N) methyl adenine in the nuclear DNA of a eucaryote, Tetrahymena pyriformis
The Journal of Cell Biology 56:697–701.

https://doi.org/10.1083/jcb.56.3.697
- PubMed
- Google Scholar
1. Grabherr MG
2. Haas BJ
3. Yassour M
4. Levin JZ
5. Thompson DA
6. Amit I
7. Adiconis X
8. Fan L
9. Raychowdhury R
10. Zeng Q
11. Chen Z
12. Mauceli E
13. Hacohen N
14. Gnirke A
15. Rhind N
16. di Palma F
17. Birren BW
18. Nusbaum C
19. Lindblad-Toh K
20. Friedman N
21. Regev A
(2011) Full-Length transcriptome assembly from RNA-Seq data without a reference genome
Nature Biotechnology 29:644–652.

https://doi.org/10.1038/nbt.1883
- PubMed
- Google Scholar
(2017) Flow cytometry sorting of nuclei enables the first global characterization of Paramecium germline DNA and transposable elements
BMC Genomics 18:327.

https://doi.org/10.1186/s12864-017-3713-7
- PubMed
- Google Scholar
(2010) New algorithms and methods to estimate maximum-likelihood phylogenies: assessing the performance of PhyML 3.0
Systematic Biology 59:307–321.

https://doi.org/10.1093/sysbio/syq010
- PubMed
- Google Scholar
1. Haas BJ
2. Delcher AL
3. Mount SM
4. Wortman JR
5. Smith RK
6. Hannick LI
7. Maiti R
8. Ronning CM
9. Rusch DB
10. Town CD
11. Salzberg SL
12. White O
(2003) Improving the Arabidopsis genome annotation using maximal transcript alignment assemblies
Nucleic Acids Research 31:5654–5666.

https://doi.org/10.1093/nar/gkg770
- PubMed
- Google Scholar
1. Haas BJ
2. Salzberg SL
3. Zhu W
4. Pertea M
5. Allen JE
6. Orvis J
7. White O
8. Buell CR
9. Wortman JR
(2008) Automated eukaryotic gene structure annotation using EVidenceModeler and the program to assemble spliced alignments
Genome Biology 9:1–22.

https://doi.org/10.1186/gb-2008-9-1-r7
- PubMed
- Google Scholar
1. Hamilton EP
2. Kapusta A
3. Huvos PE
4. Bidwell SL
5. Zafar N
6. Tang H
7. Hadjithomas M
8. Krishnakumar V
9. Badger JH
10. Caler EV
11. Russ C
(2016) Structure of the germline genome of Tetrahymena thermophila and relationship to the massively rearranged somatic genome
eLife 5:e19090.

https://doi.org/10.7554/eLife.19090
- Google Scholar
1. Hoffman DC
2. Prescott DM
(1997) Evolution of internal eliminated segments and scrambling in the micronuclear gene encoding DNA polymerase alpha in two Oxytricha species
Nucleic Acids Research 25:1883–1889.

https://doi.org/10.1093/nar/25.10.1883
- PubMed
- Google Scholar
1. Hogan DJ
2. Hewitt EA
3. Orr KE
4. Prescott DM
5. Müller KM
(2001) Evolution of IESs and scrambling in the actin I gene in hypotrichous ciliates
PNAS 98:15101–15106.

https://doi.org/10.1073/pnas.011578598
- Google Scholar
1. Huang X
2. Madan A
(1999) CAP3: a DNA sequence assembly program
Genome Research 9:868–877.

https://doi.org/10.1101/gr.9.9.868
- PubMed
- Google Scholar
(1989) Developmentally coordinated en masse excision of a highly repetitive element in E. crassus
Cell 59:1009–1018.

https://doi.org/10.1016/0092-8674(89)90757-5
- PubMed
- Google Scholar
(1993) Structures of the Euplotes crassus Tec1 and Tec2 elements: identification of putative transposase coding regions
Gene 133:71–78.

https://doi.org/10.1016/0378-1119(93)90226-s
- PubMed
- Google Scholar
1. Katz LA
2. Kovner AM
(2010) Alternative processing of scrambled genes generates protein diversity in the ciliate Chilodonella uncinata
Journal of Experimental Zoology. Part B, Molecular and Developmental Evolution 314:480–488.

https://doi.org/10.1002/jez.b.21354
- PubMed
- Google Scholar
1. Kent WJ
(2002) BLAT -- the BLAST-like alignment tool
Genome Research 12:656–664.

https://doi.org/10.1101/gr.229202
- PubMed
- Google Scholar
1. Klobutcher LA
2. Herrick G
(1995) Consensus inverted terminal repeat sequence of Paramecium IESs: resemblance to termini of Tc1-related and Euplotes Tec transposons
Nucleic Acids Research 23:2006–2013.

https://doi.org/10.1093/nar/23.11.2006
- PubMed
- Google Scholar
1. Klobutcher LA
2. Herrick GL
(1997) Developmental genome reorganization in ciliated protozoa: the transposon link
Progress in Nucleic Acid Research and Molecular Biology 56:1–62.

https://doi.org/10.1016/s0079-6603(08)61001-6
- PubMed
- Google Scholar
1. Kolmogorov M
2. Yuan J
3. Lin Y
4. Pevzner PA
(2019) Assembly of long, error-prone reads using repeat graphs
Nature Biotechnology 37:540–546.

https://doi.org/10.1038/s41587-019-0072-8
- PubMed
- Google Scholar
1. Krikau MF
2. Jahn CL
(1991) Tec2, a second transposon-like element demonstrating developmentally programmed excision in Euplotes crassus
Molecular and Cellular Biology 11:4751–4759.

https://doi.org/10.1128/mcb.11.9.4751-4759.1991
- PubMed
- Google Scholar
(2000) Evolution and assembly of an extremely scrambled gene
PNAS 97:3298–3303.

https://doi.org/10.1073/pnas.97.7.3298
- Google Scholar
1. Landweber LF
(2007) Why genomes in pieces?
Science 318:405–407.

https://doi.org/10.1126/science.1150280
- PubMed
- Google Scholar
1. Langmead B
2. Salzberg SL
(2012) Fast gapped-read alignment with Bowtie 2
Nature Methods 9:357–359.

https://doi.org/10.1038/nmeth.1923
- PubMed
- Google Scholar
(1976) DNA of ciliated protozoa: DNA sequence diminution during macronuclear development of Oxytricha
Cell 7:67–74.

https://doi.org/10.1016/0092-8674(76)90256-7
- PubMed
- Google Scholar
(2017) Thousands of RNA-cached copies of whole chromosomes are present in the ciliate Oxytricha during development
RNA 23:1200–1208.

https://doi.org/10.1261/rna.058511.116
- PubMed
- Google Scholar
(2019) Capture of complete ciliate chromosomes in single sequencing reads reveals widespread chromosome isoforms
BMC Genomics 20:1037.

https://doi.org/10.1186/s12864-019-6189-9
- PubMed
- Google Scholar
1. Lowe TM
2. Eddy SR
(1997) tRNAscan-SE: a program for improved detection of transfer RNA genes in genomic sequence
Nucleic Acids Research 25:955–964.

https://doi.org/10.1093/nar/25.5.955
- PubMed
- Google Scholar
Book
1. Lynn D
(2008)
The Ciliated Protozoa: Characterization, Classification, and Guide to the Literature

Springer Science & Business Media.
- Google Scholar
1. Manni M
2. Berkeley MR
3. Seppey M
4. Simão FA
5. Zdobnov EM
(2021) BUSCO update: novel and streamlined workflows along with broader and deeper phylogenetic coverage for scoring of eukaryotic, prokaryotic, and viral genomes
Molecular Biology and Evolution 38:4647–4654.

https://doi.org/10.1093/molbev/msab199
- PubMed
- Google Scholar
(2018a) Exploration of the germline genome of the ciliate Chilodonella uncinata through single-cell omics (transcriptomics and genomics)
MBio 9:e01836-17.

https://doi.org/10.1128/mBio.01836-17
- PubMed
- Google Scholar
(2018b) Twisted tales: insights into genome diversity of ciliates using single-cell ’omics
Genome Biology and Evolution 10:1927–1939.

https://doi.org/10.1093/gbe/evy133
- PubMed
- Google Scholar
1. Menzel P
2. Ng KL
3. Krogh A
(2016) Fast and sensitive taxonomic classification for metagenomics with Kaiju
Nature Communications 7:11257.

https://doi.org/10.1038/ncomms11257
- PubMed
- Google Scholar
1. Meyer F
2. Schmidt HJ
3. Plümper E
4. Hasilik A
5. Mersmann G
6. Meyer HE
7. Engström A
8. Heckmann K
(1991) UGA is translated as cysteine in pheromone 3 of Euplotes octocarinatus
PNAS 88:3758–3761.

https://doi.org/10.1073/pnas.88.9.3758
- PubMed
- Google Scholar
1. Miller RV
2. Neme R
3. Clay DM
4. Pathmanathan JS
5. Lu MW
6. Yerlici VT
7. Khurana JS
8. Landweber LF
(2021) Transcribed germline-limited coding sequences in Oxytricha trifallax
G3 11:jkab092.

https://doi.org/10.1093/g3journal/jkab092
- PubMed
- Google Scholar
(1992) Analysis of a scrambled gene: the gene encoding alpha-telomere-binding protein in Oxytricha nova
Genes & Development 6:788–800.

https://doi.org/10.1101/gad.6.5.788
- PubMed
- Google Scholar
(2005) Comparative genomics of nematodes
Trends in Genetics 21:573–581.

https://doi.org/10.1016/j.tig.2005.08.003
- PubMed
- Google Scholar
(2006) Interconversion of germline-limited and somatic DNA in a scrambled gene
Journal of Molecular Evolution 63:69–73.

https://doi.org/10.1007/s00239-005-0166-4
- PubMed
- Google Scholar
(2019) Excised linear introns regulate growth in yeast
Nature 565:606–611.

https://doi.org/10.1038/s41586-018-0828-1
- PubMed
- Google Scholar
1. Mourier T
2. Jeffares DC
(2003) Eukaryotic intron loss
Science 300:1393.

https://doi.org/10.1126/science.1080559
- PubMed
- Google Scholar
1. Nowacki M
2. Vijayan V
3. Zhou Y
4. Schotanus K
5. Doak TG
6. Landweber LF
(2008) RNA-Mediated epigenetic programming of a genome-rearrangement pathway
Nature 451:153–158.

https://doi.org/10.1038/nature06452
- PubMed
- Google Scholar
(2009) A functional role for transposases in a large eukaryotic genome
Science 324:935–938.

https://doi.org/10.1126/science.1170023
- PubMed
- Google Scholar
(2019) Introns are mediators of cell response to starvation
Nature 565:612–617.

https://doi.org/10.1038/s41586-018-0859-7
- PubMed
- Google Scholar
1. Parfrey LW
2. Lahr DJG
3. Knoll AH
4. Katz LA
(2011) Estimating the timing of early eukaryotic diversification with multigene molecular clocks
PNAS 108:13624–13629.

https://doi.org/10.1073/pnas.1110633108
- PubMed
- Google Scholar
1. Prescott DM
(1994) The DNA of ciliated protozoa
Microbiological Reviews 58:233–267.

https://doi.org/10.1128/mr.58.2.233-267.1994
- PubMed
- Google Scholar
(1998) Evolution of the scrambled germline gene encoding alpha-telomere binding protein in three hypotrichous ciliates
Chromosoma 107:293–303.

https://doi.org/10.1007/s004120050311
- PubMed
- Google Scholar
1. Riley JL
2. Katz LA
(2001) Widespread distribution of extensive chromosomal fragmentation in ciliates
Molecular Biology and Evolution 18:1372–1377.

https://doi.org/10.1093/oxfordjournals.molbev.a003921
- PubMed
- Google Scholar
1. Rogozin IB
2. Wolf YI
3. Sorokin AV
4. Mirkin BG
5. Koonin EV
(2003) Remarkable interkingdom conservation of intron positions and massive, lineage-specific intron loss and gain in eukaryotic evolution
Current Biology 13:1512–1517.

https://doi.org/10.1016/s0960-9822(03)00558-x
- PubMed
- Google Scholar
1. Ruan J
2. Li H
(2020) Fast and accurate long-read assembly with wtdbg2
Nature Methods 17:155–158.

https://doi.org/10.1038/s41592-019-0669-3
- PubMed
- Google Scholar
1. Schmidbaur H
2. Kawaguchi A
3. Clarence T
4. Fu X
5. Hoang OP
6. Zimmermann B
7. Ritschard EA
8. Weissenbacher A
9. Foster JS
10. Nyholm SV
11. Bates PA
12. Albertin CB
13. Tanaka E
14. Simakov O
(2022) Emergence of novel cephalopod gene regulation and expression through large-scale genome reorganization
Nature Communications 13:2172.

https://doi.org/10.1038/s41467-022-29694-7
- PubMed
- Google Scholar
(2021) BleTIES: annotation of natural genome editing in ciliates using long read sequencing
Bioinformatics 37:3929–3931.

https://doi.org/10.1093/bioinformatics/btab613
- Google Scholar
1. Sellis D
2. Guérin F
3. Arnaiz O
4. Pett W
5. Lerat E
6. Boggetto N
7. Krenek S
8. Berendonk T
9. Couloux A
10. Aury JM
11. Labadie K
12. Malinsky S
13. Bhullar S
14. Meyer E
15. Sperling L
16. Duret L
17. Duharcourt S
(2021) Massive colonization of protein-coding exons by selfish genetic elements in Paramecium germline genomes
PLOS Biology 19:e3001309.

https://doi.org/10.1371/journal.pbio.3001309
- Google Scholar
1. Sheng Y
2. Duan L
3. Cheng T
4. Qiao Y
5. Stover NA
6. Gao S
(2020) The completed macronuclear genome of a model ciliate Tetrahymena thermophila and its application in genome scrambling and copy number analyses
Science China. Life Sciences 63:1534–1542.

https://doi.org/10.1007/s11427-020-1689-4
- PubMed
- Google Scholar
1. Sievers F
2. Wilm A
3. Dineen D
4. Gibson TJ
5. Karplus K
6. Li W
7. Lopez R
8. McWilliam H
9. Remmert M
10. Söding J
11. Thompson JD
12. Higgins DG
(2011) Fast, scalable generation of high-quality protein multiple sequence alignments using Clustal Omega
Molecular Systems Biology 7:539.

https://doi.org/10.1038/msb.2011.75
- PubMed
- Google Scholar
(2015) BUSCO: assessing genome assembly and annotation completeness with single-copy orthologs
Bioinformatics 31:3210–3212.

https://doi.org/10.1093/bioinformatics/btv351
- PubMed
- Google Scholar
1. Slater GSC
2. Birney E
(2005) Automated generation of heuristics for biological sequence comparison
BMC Bioinformatics 6:1.

https://doi.org/10.1186/1471-2105-6-31
- PubMed
- Google Scholar
Website
1. Smit AF
2. Hubley R
(2008) RepeatModeler Open-1.0
Accessed September 23, 2020.

http://www.repeatmasker.org
Website
1. Smit AF
2. Hubley R
3. Green P
(2013) RepeatMasker Open-4.0
Accessed September 23, 2020.

http://www.repeatmasker.org
(2012) Genetic consequences of programmed genome rearrangement
Current Biology 22:1524–1529.

https://doi.org/10.1016/j.cub.2012.06.028
- PubMed
- Google Scholar
(2020) Combined genome and transcriptome analyses of the ciliate schmidingerella arcuata (Spirotrichea) reveal patterns of DNA elimination, scrambling, and inversion
Genome Biology and Evolution 12:1616–1622.

https://doi.org/10.1093/gbe/evaa185
- PubMed
- Google Scholar
1. Speijer D
(2008) Making sense of scrambled genomes
Science 319:901–902.

https://doi.org/10.1126/science.319.5865.901a
- PubMed
- Google Scholar
(2006) PAL2NAL: robust conversion of protein sequence alignments into the corresponding codon alignments
Nucleic Acids Research 34:W609–W612.

https://doi.org/10.1093/nar/gkl315
- PubMed
- Google Scholar
1. Swart EC
2. Bracht JR
3. Magrini V
4. Minx P
5. Chen X
6. Zhou Y
7. Khurana JS
8. Goldman AD
9. Nowacki M
10. Schotanus K
11. Jung S
12. Fulton RS
13. Ly A
14. McGrath S
15. Haub K
16. Wiggins JL
17. Storton D
18. Matese JC
19. Parsons L
20. Chang WJ
21. Bowen MS
22. Stover NA
23. Jones TA
24. Eddy SR
25. Herrick GA
26. Doak TG
27. Wilson RK
28. Mardis ER
29. Landweber LF
(2013) The Oxytricha trifallax macronuclear genome: a complex eukaryotic genome with 16,000 tiny chromosomes
PLOS Biology 11:e1001473.

https://doi.org/10.1371/journal.pbio.1001473
- PubMed
- Google Scholar
(2016) Biogeography and character evolution of the ciliate genus euplotes (spirotrichea, euplotia), with description of euplotes curdsi sp. nov
PLOS ONE 11:e0165442.

https://doi.org/10.1371/journal.pone.0165442
- PubMed
- Google Scholar
(1999) Isolation of micronuclei from euplotes octocarinatus and identification of an internal eliminated sequence in the micronuclear gene encoding γ-tubulin 2
European Journal of Protistology 35:208–216.

https://doi.org/10.1016/S0932-4739(99)80039-X
- Google Scholar
1. Thomas CA
(1971) The genetic organization of chromosomes
Annual Review of Genetics 5:237–256.

https://doi.org/10.1146/annurev.ge.05.120171.001321
- PubMed
- Google Scholar
(2012) Draft macronuclear genome of a ciliate Euplotes crassus
Molekuliarnaia Biologiia 46:361–366.

https://doi.org/10.1134/S0026893312020197
- PubMed
- Google Scholar
1. Walker BJ
2. Abeel T
3. Shea T
4. Priest M
5. Abouelliel A
6. Sakthikumar S
7. Cuomo CA
8. Zeng Q
9. Wortman J
10. Young SK
11. Earl AM
(2014) Pilon: an integrated tool for comprehensive microbial variant detection and genome assembly improvement
PLOS ONE 9:e112963.

https://doi.org/10.1371/journal.pone.0112963
- PubMed
- Google Scholar
1. Wang W
2. Zhi H
3. Chai B
4. Liang A
(2005) Cloning and sequence analysis of the micronuclear and macronuclear gene encoding Rab protein of Euplotes octocarinatus
Bioscience, Biotechnology, and Biochemistry 69:649–652.

https://doi.org/10.1271/bbb.69.649
- PubMed
- Google Scholar
1. Wang R
2. Xiong J
3. Wang W
4. Miao W
5. Liang A
(2016) High frequency of +1 programmed ribosomal frameshifting in Euplotes octocarinatus
Scientific Reports 6:21139.

https://doi.org/10.1038/srep21139
- PubMed
- Google Scholar
1. Wang JR
2. Holt J
3. McMillan L
4. Jones CD
(2018) FMLRC: hybrid long read error correction using an FM-index
BMC Bioinformatics 19:50.

https://doi.org/10.1186/s12859-018-2051-3
- PubMed
- Google Scholar
1. Wong LC
2. Landweber LF
(2006) Evolution of programmed DNA rearrangements in a scrambled gene
Molecular Biology and Evolution 23:756–763.

https://doi.org/10.1093/molbev/msj089
- PubMed
- Google Scholar
1. Yerlici VT
2. Landweber LF
(2014) Programmed genome rearrangements in the ciliate Oxytricha
Microbiology Spectrum 2:2–6.

https://doi.org/10.1128/microbiolspec.MDNA3-0025-2014
- PubMed
- Google Scholar
1. Zheng W
2. Chen J
3. Doak TG
4. Song W
5. Yan Y
(2020) ADFinder: accurate detection of programmed DNA elimination using NGS high-throughput sequencing data
Bioinformatics 36:3632–3636.

https://doi.org/10.1093/bioinformatics/btaa226
- PubMed
- Google Scholar

Article and author information

Author details

Yi Feng

Departments of Biochemistry and Molecular Biophysics and Biological Sciences, Columbia University, New York, United States

Contribution
Conceptualization, Data curation, Formal analysis, Investigation, Visualization, Methodology, Writing - original draft, Writing – review and editing

Competing interests
No competing interests declared

"This ORCID iD identifies the author of this article:" 0000-0002-2393-1700
Rafik Neme
1. Departments of Biochemistry and Molecular Biophysics and Biological Sciences, Columbia University, New York, United States
2. Department of Chemistry and Biology, Universidad del Norte, Barranquilla, Colombia
Contribution
Conceptualization, Methodology

Competing interests
No competing interests declared

"This ORCID iD identifies the author of this article:" 0000-0001-8462-5291
Leslie Y Beh

Departments of Biochemistry and Molecular Biophysics and Biological Sciences, Columbia University, New York, United States

Contribution
Investigation

Competing interests
The author is currently employed by Illumina
Xiao Chen

Pacific Biosciences, Menlo Park, United States

Contribution
Investigation

Competing interests
employed by Pacific Biosciences

"This ORCID iD identifies the author of this article:" 0000-0002-1432-268X
Jasper Braun

Department of Mathematics and Statistics, University of South Florida, Tampa, United States

Present address
Department of Pathology, Beth Israel Deaconess Medical Center, Boston, United States

Contribution
Software

Competing interests
No competing interests declared

"This ORCID iD identifies the author of this article:" 0000-0003-1250-4399
Michael W Lu

Departments of Biochemistry and Molecular Biophysics and Biological Sciences, Columbia University, New York, United States

Contribution
Investigation, Validation

Competing interests
No competing interests declared

"This ORCID iD identifies the author of this article:" 0000-0002-4926-8839
Laura F Landweber

Departments of Biochemistry and Molecular Biophysics and Biological Sciences, Columbia University, New York, United States

Contribution
Conceptualization, Supervision, Funding acquisition, Investigation, Project administration, Writing – review and editing

For correspondence
Laura.Landweber@columbia.edu

Competing interests
No competing interests declared

"This ORCID iD identifies the author of this article:" 0000-0002-7030-8540

Funding

National Institutes of Health (R35GM122555)

Laura F Landweber
Yi Feng
Leslie Y Beh
Michael W Lu

National Science Foundation (DMS1764366)

Yi Feng

Pew Latin American Fellows Program

Rafik Neme

National Science Foundation (DBI1062432)

Yi Feng

National Science Foundation (ABI1458641)

Yi Feng

National Science Foundation (ABI1759906)

Yi Feng

The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.

Acknowledgements

We thank Toshinobu Suzaki (Kobe University) for the gift of E. woodruffi (strain Iz01 from Shizuoka Prefecture) and Chlorogonium capillatum. We thank Sheela George for laboratory support and help with cell collection. We thank David Dai, Eoghan Harrington, John Beaulaurier, and Sissel Juul at Oxford Nanopore Technologies in New York for providing sequencing and advice. We thank Robert Sebra and Melissa Smith for advice and PacBio sequencing. We thank Takahiko Akematsu, Lorraine Symington, and Lea Marie for helping with PFGE. We thank Kaiyi Zhu, Shaojie He, Molly Przeworski, Harmen Bussemaker, and Nataša Jonoska for advice on Monte Carlo simulations. We also thank Scott Roy, Samuel Sternberg, Bill Jack, and all current and past Landweber lab members for discussion about the origin of scrambled genes, as well as David Prescott and Klaus Heckmann for inspiration, and Margarita T Angelova, Sindhuja Devanapally, Danylo Villano and Kehan Bao for comments on the manuscript. This work was supported by the National Institutes of Health, R35GM122555, and National Science Foundation, DMS1764366, and the National Center for Genome Analysis Support computing resources (supported by National Science Foundation DBI1062432, ABI1458641, and ABI1759906 to Indiana University). Rafik Neme was supported by the Pew Latin American Fellows Program.

Copyright

This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited.