The Paramecium Germline Genome Provides a Niche for Intragenic Parasitic DNA: Evolutionary Dynamics of Internal Eliminated Sequences

Insertions of parasitic DNA within coding sequences are usually deleterious and are generally counter-selected during evolution. Thanks to nuclear dimorphism, ciliates provide unique models to study the fate of such insertions. Their germline genome undergoes extensive rearrangements during development of a new somatic macronucleus from the germline micronucleus following sexual events. In Paramecium, these rearrangements include precise excision of unique-copy Internal Eliminated Sequences (IES) from the somatic DNA, requiring the activity of a domesticated piggyBac transposase, PiggyMac. We have sequenced Paramecium tetraurelia germline DNA, establishing a genome-wide catalogue of ∼45,000 IESs, in order to gain insight into their evolutionary origin and excision mechanism. We obtained direct evidence that PiggyMac is required for excision of all IESs. Homology with known P. tetraurelia Tc1/mariner transposons, described here, indicates that at least a fraction of IESs derive from these elements. Most IES insertions occurred before a recent whole-genome duplication that preceded diversification of the P. aurelia species complex, but IES invasion of the Paramecium genome appears to be an ongoing process. Once inserted, IESs decay rapidly by accumulation of deletions and point substitutions. Over 90% of the IESs are shorter than 150 bp and present a remarkable size distribution with a ∼10 bp periodicity, corresponding to the helical repeat of double-stranded DNA and suggesting DNA loop formation during assembly of a transpososome-like excision complex. IESs are equally frequent within and between coding sequences; however, excision is not 100% efficient and there is selective pressure against IES insertions, in particular within highly expressed genes. We discuss the possibility that ancient domestication of a piggyBac transposase favored subsequent propagation of transposons throughout the germline by allowing insertions in coding sequences, a fraction of the genome in which parasitic DNA is not usually tolerated.     

Excision of Helitron Transposons in Maize

Helitrons are novel transposons discovered by bioinformatic analysis of eukaryotic genome sequences. They are believed to move by rolling circle (RC) replication because their predicted transposases are homologous to those of bacterial RC transposons. We report here evidence of somatic Helitron excision in maize, an unexpected finding suggesting that Helitrons can exhibit an excisive mode of transposition.     

Genome gymnastics: unique modes of DNA evolution and processing in ciliates

In some ciliates, the DNA sequences of the germline genomes have been profoundly modified during evolution, providing unprecedented examples of germline DNA malleability. Although the significance of the modifications and malleability is unclear, they may reflect the evolution of mechanisms that facilitate evolution. Because of the modifications, these ciliates must perform remarkable feats of cutting, splicing, rearrangement and elimination of DNA sequences to convert the chromosomal DNA in the germline genome (micronuclear genome) into gene-sized DNA molecules in the somatic genome (macronuclear genome). How these manipulations of DNA are guided and carried out is largely unknown. However, the organization and manipulation of ciliate DNA sequences are new phenomena that expand a general appreciation for the flexibility of DNA in evolution and development.     

Comparative analysis of single-cell genome sequencing techniques toward the characterization of germline and somatic genomes in ciliated protists

Ciliated protists contain both germline micronucleus (MIC) and somatic macronucleus (MAC) in a single cytoplasm. Programmed genome rearrangements occur in ciliates during sexual processes, and the extent of rearrangements varies dramatically among species, which lead to significant differences in genomic architectures. However, genomic sequences remain largely unknown for most ciliates due to the difficulty in culturing and in separating the germline from the somatic genome in a single cell. Single-cell whole genome amplification (WGA) has emerged as a powerful technology to characterize the genomic heterogeneity at the single-cell level. In this study, we compared two single-cell WGA, multiple displacement amplification (MDA) and multiple annealing and looping-based amplification cycles (MALBAC) in characterizing the germline and somatic genomes in ciliates with different genomic architectures. Our results showed that: 1) MALBAC exhibits strong amplification bias towards MAC genome while MDA shows bias towards MIC genome of ciliates with extensively fragmented MAC genome; 2) both MDA and MALBAC could amplify MAC genome more efficiently in ciliates with moderately fragmented MAC genome. Moreover, we found that more sample replicates could help to obtain more genomic data. Our work provides a reference for selecting the appropriate method to characterize germline and somatic genomes of ciliates.     

Nowa1p and Nowa2p: novel putative RNA binding proteins involved in trans-nuclear crosstalk in Paramecium tetraurelia

BACKGROUND: The germline genome of ciliates is extensively rearranged during development of a new somatic macronucleus from the germline micronucleus, a process that follows sexual events. In Paramecium tetraurelia, single-copy internal eliminated sequences (IESs) and multicopy transposons are eliminated, whereas cellular genes are amplified to approximately 800 n. For a subset of IESs, introduction of the IES sequence into the maternal (prezygotic) macronucleus specifically inhibits excision of the homologous IES in the developing zygotic macronucleus. This and other homology-dependent maternal effects have suggested that rearrangement patterns are epigenetically determined by an RNA-mediated, trans-nuclear comparison, involving the RNA interference pathway, of germline and somatic genomes. RESULTS: We report the identification of novel developmentally regulated RNA binding proteins, Nowa1p and Nowa2p, which are required for the survival of sexual progeny. Green fluorescent protein (GFP) fusions show that Nowa1p accumulates into the maternal macronucleus shortly before meiosis of germline micronuclei and is later transported to developing macronuclei. Nowa1p/2p depletion impairs the elimination of transposons and of those IESs that are controlled by maternal effects, confirming the existence of distinct IES classes. CONCLUSIONS: The results indicate that Nowa proteins are essential components of the trans-nuclear-crosstalk mechanism that is responsible for epigenetic programming of genome rearrangements. We discuss implications for the current models of genome scanning in ciliates, a process related to the formation of heterochromatin by RNA interference in other eukaryotes.     

Early developmental,meiosis-specific proteins — Spo11, Msh4-1, and Msh5 — Affect subsequent genome reorganization in Paramecium tetraurelia

Developmental DNA elimination in Paramecium tetraurelia occurs through a trans-nuclear comparison of the genomes of two distinct types of nuclei: the germline micronucleus (MIC) and the somatic macronucleus (MAC). During sexual reproduction, which starts with meiosis of the germline nuclei, MIC-limited sequences including Internal Eliminated Sequences (IESs) and transposons are eliminated from the developing MAC in a process guided by noncoding RNAs (scnRNAs and iesRNAs). However, our current understanding of this mechanism is still very limited. Therefore, studying both genetic and epigenetic aspects of these processes is a crucial step to understand this phenomenon in more detail. Here, we describe the involvement of homologs of classical meiotic proteins, Spo11, Msh4-1, and Msh5 in this phenomenon. Based on our analyses, we propose that proper functioning of Spo11, Msh4-1, and Msh5 during Paramecium sexual reproduction are necessary for genome reorganization and viable progeny. Also, we show that double-strand breaks (DSBs) in DNA induced during meiosis by Spo11 are crucial for proper IESs excision. In summary, our investigations show that early sexual reproduction processes may significantly influence later somatic genome integrity.     

Ku-Mediated Coupling of DNA Cleavage and Repair during Programmed Genome Rearrangements in the Ciliate Paramecium tetraurelia

During somatic differentiation, physiological DNA double-strand breaks (DSB) can drive programmed genome rearrangements (PGR), during which DSB repair pathways are mobilized to safeguard genome integrity. Because of their unique nuclear dimorphism, ciliates are powerful unicellular eukaryotic models to study the mechanisms involved in PGR. At each sexual cycle, the germline nucleus is transmitted to the progeny, but the somatic nucleus, essential for gene expression, is destroyed and a new somatic nucleus differentiates from a copy of the germline nucleus. In Paramecium tetraurelia, the development of the somatic nucleus involves massive PGR, including the precise elimination of at least 45,000 germline sequences (Internal Eliminated Sequences, IES). IES excision proceeds through a cut-and-close mechanism: a domesticated transposase, PiggyMac, is essential for DNA cleavage, and DSB repair at excision sites involves the Ligase IV, a specific component of the non-homologous end-joining (NHEJ) pathway. At the genome-wide level, a huge number of programmed DSBs must be repaired during this process to allow the assembly of functional somatic chromosomes. To understand how DNA cleavage and DSB repair are coordinated during PGR, we have focused on Ku, the earliest actor of NHEJ-mediated repair. Two Ku70 and three Ku80 paralogs are encoded in the genome of P. tetraurelia: Ku70a and Ku80c are produced during sexual processes and localize specifically in the developing new somatic nucleus. Using RNA interference, we show that the development-specific Ku70/Ku80c heterodimer is essential for the recovery of a functional somatic nucleus. Strikingly, at the molecular level, PiggyMac-dependent DNA cleavage is abolished at IES boundaries in cells depleted for Ku80c, resulting in IES retention in the somatic genome. PiggyMac and Ku70a/Ku80c co-purify as a complex when overproduced in a heterologous system. We conclude that Ku has been integrated in the Paramecium DNA cleavage factory, enabling tight coupling between DSB introduction and repair during PGR.     

Characterization of a Tc1-like transposable element in zebrafish (Danio rerio)

We have characterized Tdr1, a family of Tc1-like transposable elements found in the genome of zebrafish (Danio rerio). The copy number and distribution of the sequence in the zebrafish genome have been determined, and by these criteria Tdr1 can be classified as a moderately repetitive, interspersed element. Examination of the sequences and structures of several copies of Tdr1 revealed that a particular deletion derivative, 1250 by long, of the transposon has been amplified to become the dominant form of Tdr1. The deletion in these elements encompasses sequences encoding the N-terminal portion of the putative Tdr1 transposase. Sequences corresponding to the deleted region were also detected, and thus allowed prediction of the nucleotide sequence of a hypothetical full-length element. Well conserved segments of Tc1-like transposons were found in the flanking regions of known fish genes, suggesting that these elements have a long evolutionary history in piscine genomes. Tdr1 elements have long, 208 by inverted repeats, with a short DNA motif repeated four times at the termini of the inverted repeats. Although different from that of the prototype C. elegans transposon Tc1, this inverted repeat structure is shared by transposable elements from salmonid fish species and two Drosophila species. We propose that these transposons form a subgroup within the Tc1-like family. Comparison of Tc1-like transposons supports the hypothesis that the transposase genes and their flanking sequences have been shaped by independent evolutionary constraints. Although Tc1-like sequences are present in the genomes of several strains of zebrafish and in salmonid fishes, these sequences are not conserved in the genus Danio, thus raising the possibility that these elements can be exploited for gene tagging and genome mapping.     

A family of Tc1-like transposons from the genomes of fishes and frogs: evidence for horizontal transmission.

Tc1-like transposons are very widely distributed within the genomes of animal species. They consist of an inverted repeat sequence flanking a transposase gene with homology to the mobile DNA element, Tc1 of the nematode Caenorhabditis elegans. These elements seem particularly to infest the genomes of fish and amphibian species where they can account for 1% of the total genome. However, all vertebrate Tc1-like elements isolated so far are non-functional in that they contain multiple frameshifts within their transposase coding regions. Here I describe a Tc1-like transposon (PPTN) from the genome of a marine flatfish species (Pleuronectes platessa) which bears conserved inverted repeats flanking an apparently intact transposase gene. Closely related, although degenerate, Tc1-like transposons were also isolated from the genomes of Atlantic salmon (SSTN, Salmo salar) and frog (RTTN, Rana temporaria). Consensual nucleic acid sequences were derived by comparing several individual isolates from each species and conceptual amino acid sequences were thence derived for their transposases. Phylogenetic analysis of these sequences with previously isolated Tc1-like transposases shows that the elements from plaice, salmon and frog comprise a new subfamily of Tc1-like transposons. Each member is distinct in that it is not found in the genomes of the other species tested. Plaice genomes contain about 300 copies of PPTN, salmon 1200 copies of SSTN and frog genomes about 500 copies of RTTN. The presence of these closely related elements in the genomes of fish and frog species, representing evolutionary lines, which diverged more than 400 million years ago, is not consistent with a vertical transmission model for their distributions.     

Developmentally programmed DNA rearrangement in Tetrahymena thermophila: Isolation and sequence characterization of three new alternative deletion systems

Extensive developmentally programmed DNA rearrangements, including thousands of internal deletions, occur in the differentiating somatic macronucleus in Tetrahymena thermophila. Some deletion systems involve the use of multiple alternative deletion sites. We report here the cloning and the sequences of three new alternative deletion systems (RR, RP and B) obtained using genomic subtraction. The RP and RR deletion systems are 2 kb apart on chromosome IR, and both involve the removal of < 2 kb of micronuclear sequences. The B deletion system is on chromosome 5 and involves a deletion of > 5 kb. All three deleted regions are very AT rich (∼ 80%) and do not appear to encode any protein. Sequences of the regions flanking the deletion junctions of all three systems revealed no sequence similarity among them nor with any previously reported deletion systems, suggesting that different cis-acting elements are involved for rearrangement. Unlike other deletion systems in ciliates, the B deletion system lacks short terminal direct repeats. Our results suggest an average of at least one alternative deletion system per 134 kb of micronuclear DNA and lead to an estimate that at least 25% of all deletion systems in Tetrahymena utilize alternative ends. The genomic subtraction method employed in this study could prove useful for the isolation of alternatively deleted DNA in special-purpose cases in Tetrahymena and other ciliates. The hybridization parameters for genomic subtraction worked out here for highly AT-rich DNA may have wider usefulness.     

Random sequencing of Paramecium somatic DNA

We report a random survey of 1 to 2% of the somatic genome of the free-living ciliate Paramecium tetraurelia by single-run sequencing of the ends of plasmid inserts. As in all ciliates, the germ line genome of Paramecium (100 to 200 Mb) is reproducibly rearranged at each sexual cycle to produce a somatic genome of expressed or potentially expressed genes, stripped of repeated sequences, transposons, and AT-rich unique sequence elements limited to the germ line. We found the somatic genome to be compact (>68% coding, estimated from the sequence of several complete library inserts) and to feature uniformly small introns (18 to 35 nucleotides). This facilitated gene discovery: 722 open reading frames (ORFs) were identified by similarity with known proteins, and 119 novel ORFs were tentatively identified by internal comparison of the data set. We determined the phylogenetic position of Paramecium with respect to eukaryotes whose genomes have been sequenced by the distance matrix neighbor-joining method by using random combined protein data from the project. The unrooted tree obtained is very robust and in excellent agreement with accepted topology, providing strong support for the quality and consistency of the data set. Our study demonstrates that a random survey of the somatic genome of Paramecium is a good strategy for gene discovery in this organism.     

Local Effect of Enhancer of Zeste-Like Reveals Cooperation of Epigenetic and cis-Acting Determinants for Zygotic Genome Rearrangements

In the ciliate Paramecium tetraurelia, differentiation of the somatic nucleus from the zygotic nucleus is characterized by massive and reproducible deletion of transposable elements and of 45,000 short, dispersed, single-copy sequences. A specific class of small RNAs produced by the germline during meiosis, the scnRNAs, are involved in the epigenetic regulation of DNA deletion but the underlying mechanisms are poorly understood. Here, we show that trimethylation of histone H3 (H3K27me3 and H3K9me3) displays a dynamic nuclear localization that is altered when the endonuclease required for DNA elimination is depleted. We identified the putative histone methyltransferase Ezl1 necessary for H3K27me3 and H3K9me3 establishment and show that it is required for correct genome rearrangements. Genome-wide analyses show that scnRNA-mediated H3 trimethylation is necessary for the elimination of long, repeated germline DNA, while single copy sequences display differential sensitivity to depletion of proteins involved in the scnRNA pathway, Ezl1- a putative histone methyltransferase and Dcl5- a protein required for iesRNA biogenesis. Our study reveals cis-acting determinants, such as DNA length, also contribute to the definition of germline sequences to delete. We further show that precise excision of single copy DNA elements, as short as 26 bp, requires Ezl1, suggesting that development specific H3K27me3 and H3K9me3 ensure specific demarcation of very short germline sequences from the adjacent somatic sequences.     

Varied truncation and clustering characterize some short repeats identified in micronucleus-specific DNA of Tetrahymena thermophila

There are over 6000 internally eliminated DNA sequences (IESs) in the Tetrahymena genome that are deleted in a programmed fashion during the development of a polyploid, somatic macronucleus from a diploid germline micronucleus. Recently, based on several results, a homology and small RNA-based mechanism has been proposed for the efficient elimination of IES elements. Since the RNAi machinery is proposed to be intimately involved in silencing potentially harmful repeats such as transposons and viruses, characterization of repeats and the conditions for their developmental elimination from the somatic genome is warranted. Three short (500–600 bp) repeat families, members of which had been experimentally identified in IESs, that is, in micronucleus-specific DNA, are examined here using the Tetrahymena genome database. Members of all three families display varied degrees of truncation and are represented in macronuclear sequences. A 200 bp segment of one of the families can appear in the genome on its own, or as part of a 600 bp repeat detected experimentally, or in association with an unrelated 1 kb sequence to form a 1.2 kb repeat that is also frequently truncated. The 1 kb sequence contains a 300 bp section similar to a repeat associated with a non-long terminal repeat-like element and is often found accompanied by several more copies of this shorter repeat. These observations indicate that transposition may have had a role in the evolution of the short repeat families.     

NAR Breakthrough Article: Paramecium tetraurelia chromatin assembly factor-1-like protein PtCAF-1 is involved in RNA-mediated control of DNA elimination

Genome-wide DNA remodelling in the ciliate Paramecium is ensured by RNA-mediated trans-nuclear crosstalk between the germline and the somatic genomes during sexual development. The rearrangements include elimination of transposable elements, minisatellites and tens of thousands non-coding elements called internally eliminated sequences (IESs). The trans-nuclear genome comparison process employs a distinct class of germline small RNAs (scnRNAs) that are compared against the parental somatic genome to select the germline-specific subset of scnRNAs that subsequently target DNA elimination in the progeny genome. Only a handful of proteins involved in this process have been identified so far and the mechanism of DNA targeting is unknown. Here we describe chromatin assembly factor-1-like protein (PtCAF-1), which we show is required for the survival of sexual progeny and localizes first in the parental and later in the newly developing macronucleus. Gene silencing shows that PtCAF-1 is required for the elimination of transposable elements and a subset of IESs. PTCAF-1 depletion also impairs the selection of germline-specific scnRNAs during development. We identify specific histone modifications appearing during Paramecium development which are strongly reduced in PTCAF-1 depleted cells. Our results demonstrate the importance of PtCAF-1 for the epigenetic trans-nuclear cross-talk mechanism.     

Evolution of Germline-Limited Sequences in Two Populations of the Ciliate Chilodonella uncinata

Ciliates are microbial eukaryotes that separate their nuclear functions into a germline micronucleus and a somatic macronucleus. During development of the macronucleus the genome undergoes a series of reorganization events that includes the precise excision of intervening DNA. Here, we determine the architecture of four loci in the micronuclear and macronuclear genomes of the ciliate Chilodonella uncinata and compare the levels of variation in micronuclear-limited sequences to macronuclear destined sequences at two of these loci. We find that within a population, germline-limited sequences are evolving at the same rate as other putatively neutral sites, but between populations germline-limited sequences are accumulating mutations at a much faster rate than other sites. We also find evidence of macronuclear recombination and incomplete elimination of intervening DNA, which result in increased diversity in the macronuclear genome. Our results support the assertion that the unusual genomic features of ciliates can result in rapid and unpredicted patterns of diversification.     

Generation of Tandem Direct Duplications by Reversed-Ends Transposition of Maize Ac Elements

Tandem direct duplications are a common feature of the genomes of eukaryotes ranging from yeast to human, where they comprise a significant fraction of copy number variations. The prevailing model for the formation of tandem direct duplications is non-allelic homologous recombination (NAHR). Here we report the isolation of a series of duplications and reciprocal deletions isolated de novo from a maize allele containing two Class II Ac/Ds transposons. The duplication/deletion structures suggest that they were generated by alternative transposition reactions involving the termini of two nearby transposable elements. The deletion/duplication breakpoint junctions contain 8 bp target site duplications characteristic of Ac/Ds transposition events, confirming their formation directly by an alternative transposition mechanism. Tandem direct duplications and reciprocal deletions were generated at a relatively high frequency (∼0.5 to 1%) in the materials examined here in which transposons are positioned nearby each other in appropriate orientation; frequencies would likely be much lower in other genotypes. To test whether this mechanism may have contributed to maize genome evolution, we analyzed sequences flanking Ac/Ds and other hAT family transposons and identified three small tandem direct duplications with the structural features predicted by the alternative transposition mechanism. Together these results show that some class II transposons are capable of directly inducing tandem sequence duplications, and that this activity has contributed to the evolution of the maize genome.     

Massive colonization of protein-coding exons by selfish genetic elements in Paramecium germline genomes

Ciliates are unicellular eukaryotes with both a germline genome and a somatic genome in the same cytoplasm. The somatic macronucleus (MAC), responsible for gene expression, is not sexually transmitted but develops from a copy of the germline micronucleus (MIC) at each sexual generation. In the MIC genome of Paramecium tetraurelia, genes are interrupted by tens of thousands of unique intervening sequences called internal eliminated sequences (IESs), which have to be precisely excised during the development of the new MAC to restore functional genes. To understand the evolutionary origin of this peculiar genomic architecture, we sequenced the MIC genomes of 9 Paramecium species (from approximately 100 Mb in Paramecium aurelia species to >1.5 Gb in Paramecium caudatum). We detected several waves of IES gains, both in ancestral and in more recent lineages. While the vast majority of IESs are single copy in present-day genomes, we identified several families of mobile IESs, including nonautonomous elements acquired via horizontal transfer, which generated tens to thousands of new copies. These observations provide the first direct evidence that transposable elements can account for the massive proliferation of IESs in Paramecium. The comparison of IESs of different evolutionary ages indicates that, over time, IESs shorten and diverge rapidly in sequence while they acquire features that allow them to be more efficiently excised. We nevertheless identified rare cases of IESs that are under strong purifying selection across the aurelia clade. The cases examined contain or overlap cellular genes that are inactivated by excision during development, suggesting conserved regulatory mechanisms. Similar to the evolution of introns in eukaryotes, the evolution of Paramecium IESs highlights the major role played by selfish genetic elements in shaping the complexity of genome architecture and gene expression.

A comparative genomics study of nine Paramecium species reveals successful invasion of genes by transposable elements in their germline genomes, showing that the internal eliminated sequences (IESs) followed an evolutionary trajectory remarkably similar to that of spliceosomal introns.