Developmental timing of programmed DNA elimination in Paramecium tetraurelia recapitulates germline transposon evolutionary dynamics.

With its nuclear dualism, the ciliate Paramecium constitutes a unique model to study how host genomes cope with transposable elements (TEs). P. tetraurelia harbors two germline micronuclei (MICs) and a polyploid somatic macronucleus (MAC) that develops from one MIC at each sexual cycle. Throughout evolution, the MIC genome has been continuously colonized by TEs and related sequences that are removed from the somatic genome during MAC development. Whereas TE elimination is generally imprecise, excision of approximately 45,000 TE-derived internal eliminated sequences (IESs) is precise, allowing for functional gene assembly. Programmed DNA elimination is concomitant with genome amplification. It is guided by noncoding RNAs and repressive chromatin marks. A subset of IESs is excised independently of this epigenetic control, raising the question of how IESs are targeted for elimination. To gain insight into the determinants of IES excision, we established the developmental timing of DNA elimination genome-wide by combining fluorescence-assisted nuclear sorting with high-throughput sequencing. Essentially all IESs are excised within only one endoreplication round (32C to 64C), whereas TEs are eliminated at a later stage. We show that DNA elimination proceeds independently of replication. We defined four IES classes according to excision timing. The earliest excised IESs tend to be independent of epigenetic factors, display strong sequence signals at their ends, and originate from the most ancient integration events. We conclude that old IESs have been optimized during evolution for early and accurate excision by acquiring stronger sequence determinants and escaping epigenetic control.

GC content, but not nucleosome positioning, directly contributes to intron splicing efficiency in Paramecium.

Eukaryotic genes are interrupted by introns that must be accurately spliced from mRNA precursors. With an average length of 25 nt, the more than 90,000 introns of Paramecium tetraurelia stand among the shortest introns reported in eukaryotes. The mechanisms specifying the correct recognition of these tiny introns remain poorly understood. Splicing can occur cotranscriptionally, and it has been proposed that chromatin structure might influence splice site recognition. To investigate the roles of nucleosome positioning in intron recognition, we determined the nucleosome occupancy along the P. tetraurelia genome. We show that P. tetraurelia displays a regular nucleosome array with a nucleosome repeat length of ∼151 bp, among the smallest periodicities reported. Our analysis has revealed that introns are frequently associated with inter-nucleosomal DNA, pointing to an evolutionary constraint favoring introns at the AT-rich nucleosome edge sequences. Using accurate splicing efficiency data from cells depleted for nonsense-mediated decay effectors, we show that introns located at the edge of nucleosomes display higher splicing efficiency than those at the center. However, multiple regression analysis indicates that the low GC content of introns, rather than nucleosome positioning, is associated with high splicing efficiency. Our data reveal a complex link between GC content, nucleosome positioning, and intron evolution in Paramecium.

Inter-generational nuclear crosstalk links the control of gene expression to programmed genome rearrangement during the Paramecium sexual cycle.

Multinucleate cells are found in many eukaryotes, but how multiple nuclei coordinate their functions is still poorly understood. In the cytoplasm of the ciliate Paramecium tetraurelia, two micronuclei (MIC) serving sexual reproduction coexist with a somatic macronucleus (MAC) dedicated to gene expression. During sexual processes, the MAC is progressively destroyed while still ensuring transcription, and new MACs develop from copies of the zygotic MIC. Several gene clusters are successively induced and switched off before vegetative growth resumes. Concomitantly, programmed genome rearrangement (PGR) removes transposons and their relics from the new MACs. Development of the new MACs is controlled by the old MAC, since the latter expresses genes involved in PGR, including the PGM gene encoding the essential PiggyMac endonuclease that cleaves the ends of eliminated sequences. Using RNA deep sequencing and transcriptome analysis, we show that impairing PGR upregulates key known PGR genes, together with ∼600 other genes possibly also involved in PGR. Among these genes, 42% are no longer induced when no new MACs are formed, including 180 genes that are co-expressed with PGM under all tested conditions. We propose that bi-directional crosstalk between the two coexisting generations of MACs links gene expression to the progression of MAC development.

Functional diversification of Paramecium Ku80 paralogs safeguards genome integrity during precise programmed DNA elimination.

Gene duplication and diversification drive the emergence of novel functions during evolution. Because of whole genome duplications, ciliates from the Paramecium aurelia group constitute a remarkable system to study the evolutionary fate of duplicated genes. Paramecium species harbor two types of nuclei: a germline micronucleus (MIC) and a somatic macronucleus (MAC) that forms from the MIC at each sexual cycle. During MAC development, ~45,000 germline Internal Eliminated Sequences (IES) are excised precisely from the genome through a 'cut-and-close' mechanism. Here, we have studied the P. tetraurelia paralogs of KU80, which encode a key DNA double-strand break repair factor involved in non-homologous end joining. The three KU80 genes have different transcription patterns, KU80a and KU80b being constitutively expressed, while KU80c is specifically induced during MAC development. Immunofluorescence microscopy and high-throughput DNA sequencing revealed that Ku80c stably anchors the PiggyMac (Pgm) endonuclease in the developing MAC and is essential for IES excision genome-wide, providing a molecular explanation for the previously reported Ku-dependent licensing of DNA cleavage at IES ends. Expressing Ku80a under KU80c transcription signals failed to complement a depletion of endogenous Ku80c, indicating that the two paralogous proteins have distinct properties. Domain-swap experiments identified the α/ß domain of Ku80c as the major determinant for its specialized function, while its C-terminal part is required for excision of only a small subset of IESs located in IES-dense regions. We conclude that Ku80c has acquired the ability to license Pgm-dependent DNA cleavage, securing precise DNA elimination during programmed rearrangements. The present study thus provides novel evidence for functional diversification of genes issued from a whole-genome duplication.

Sequence-specific DNA binding activity of the cross-brace zinc finger motif of the piggyBac transposase.

The piggyBac transposase (PB) is distinguished by its activity and utility in genome engineering, especially in humans where it has highly promising therapeutic potential. Little is known, however, about the structure-function relationships of the different domains of PB. Here, we demonstrate in vitro and in vivo that its C-terminal Cysteine-Rich Domain (CRD) is essential for DNA breakage, joining and transposition and that it binds to specific DNA sequences in the left and right transposon ends, and to an additional unexpectedly internal site at the left end. Using NMR, we show that the CRD adopts the specific fold of the cross-brace zinc finger protein family. We determine the interaction interfaces between the CRD and its target, the 5'-TGCGT-3'/3'-ACGCA-5' motifs found in the left, left internal and right transposon ends, and use NMR results to propose docking models for the complex, which are consistent with our site-directed mutagenesis data. Our results provide support for a model of the PB/DNA interactions in the context of the transpososome, which will be useful for the rational design of PB mutants with increased activity.

Multimerization properties of PiggyMac, a domesticated piggyBac transposase involved in programmed genome rearrangements.

During sexual processes, the ciliate Paramecium eliminates 25-30% of germline DNA from its somatic genome. DNA elimination includes excision of ∼45 000 short, single-copy internal eliminated sequences (IESs) and depends upon PiggyMac (Pgm), a domesticated piggyBac transposase that is essential for DNA cleavage at IES ends. Pgm carries a core transposase region with a putative catalytic domain containing three conserved aspartic acids, and a downstream cysteine-rich (CR) domain. A C-terminal extension of unknown function is predicted to adopt a coiled-coil (CC) structure. To address the role of the three domains, we designed an in vivo complementation assay by expressing wild-type or mutant Pgm-GFP fusions in cells depleted for their endogenous Pgm. The DDD triad and the CR domain are essential for Pgm activity and mutations in either domain have a dominant-negative effect in wild-type cells. A mutant lacking the CC domain is partially active in the presence of limiting Pgm amounts, but inactive when Pgm is completely absent, suggesting that presence of the mutant protein increases the overall number of active complexes. We conclude that IES excision involves multiple Pgm subunits, of which at least a fraction must contain the CC domain.

TFIIS-Dependent Non-coding Transcription Regulates Developmental Genome Rearrangements.

Because of their nuclear dimorphism, ciliates provide a unique opportunity to study the role of non-coding RNAs (ncRNAs) in the communication between germline and somatic lineages. In these unicellular eukaryotes, a new somatic nucleus develops at each sexual cycle from a copy of the zygotic (germline) nucleus, while the old somatic nucleus degenerates. In the ciliate Paramecium tetraurelia, the genome is massively rearranged during this process through the reproducible elimination of repeated sequences and the precise excision of over 45,000 short, single-copy Internal Eliminated Sequences (IESs). Different types of ncRNAs resulting from genome-wide transcription were shown to be involved in the epigenetic regulation of genome rearrangements. To understand how ncRNAs are produced from the entire genome, we have focused on a homolog of the TFIIS elongation factor, which regulates RNA polymerase II transcriptional pausing. Six TFIIS-paralogs, representing four distinct families, can be found in P. tetraurelia genome. Using RNA interference, we showed that TFIIS4, which encodes a development-specific TFIIS protein, is essential for the formation of a functional somatic genome. Molecular analyses and high-throughput DNA sequencing upon TFIIS4 RNAi demonstrated that TFIIS4 is involved in all kinds of genome rearrangements, including excision of ~48% of IESs. Localization of a GFP-TFIIS4 fusion revealed that TFIIS4 appears specifically in the new somatic nucleus at an early developmental stage, before IES excision. RT-PCR experiments showed that TFIIS4 is necessary for the synthesis of IES-containing non-coding transcripts. We propose that these IES+ transcripts originate from the developing somatic nucleus and serve as pairing substrates for germline-specific short RNAs that target elimination of their homologous sequences. Our study, therefore, connects the onset of zygotic non coding transcription to the control of genome plasticity in Paramecium, and establishes for the first time a specific role of TFIIS in non-coding transcription in eukaryotes.

Improved methods and resources for paramecium genomics: transcription units, gene annotation and gene expression.

BACKGROUND: The 15 sibling species of the Paramecium aurelia cryptic species complex emerged after a whole genome duplication that occurred tens of millions of years ago. Given extensive knowledge of the genetics and epigenetics of Paramecium acquired over the last century, this species complex offers a uniquely powerful system to investigate the consequences of whole genome duplication in a unicellular eukaryote as well as the genetic and epigenetic mechanisms that drive speciation. High quality Paramecium gene models are important for research using this system. The major aim of the work reported here was to build an improved gene annotation pipeline for the Paramecium lineage. RESULTS: We generated oriented RNA-Seq transcriptome data across the sexual process of autogamy for the model species Paramecium tetraurelia. We determined, for the first time in a ciliate, candidate P. tetraurelia transcription start sites using an adapted Cap-Seq protocol. We developed TrUC, multi-threaded Perl software that in conjunction with TopHat mapping of RNA-Seq data to a reference genome, predicts transcription units for the annotation pipeline. We used EuGene software to combine annotation evidence. The high quality gene structural annotations obtained for P. tetraurelia were used as evidence to improve published annotations for 3 other Paramecium species. The RNA-Seq data were also used for differential gene expression analysis, providing a gene expression atlas that is more sensitive than the previously established microarray resource. CONCLUSIONS: We have developed a gene annotation pipeline tailored for the compact genomes and tiny introns of Paramecium species. A novel component of this pipeline, TrUC, predicts transcription units using Cap-Seq and oriented RNA-Seq data. TrUC could prove useful beyond Paramecium, especially in the case of high gene density. Accurate predictions of 3' and 5' UTR will be particularly valuable for studies of gene expression (e.g. nucleosome positioning, identification of cis regulatory motifs). The P. tetraurelia improved transcriptome resource, gene annotations for P. tetraurelia, P. biaurelia, P. sexaurelia and P. caudatum, and Paramecium-trained EuGene configuration are available through ParameciumDB ( http://paramecium.i2bc.paris-saclay.fr ). TrUC software is freely distributed under a GNU GPL v3 licence ( https://github.com/oarnaiz/TrUC ).

Is non-homologous end-joining really an inherently error-prone process?

DNA double-strand breaks (DSBs) are harmful lesions leading to genomic instability or diversity. Non-homologous end-joining (NHEJ) is a prominent DSB repair pathway, which has long been considered to be error-prone. However, recent data have pointed to the intrinsic precision of NHEJ. Three reasons can account for the apparent fallibility of NHEJ: 1) the existence of a highly error-prone alternative end-joining process; 2) the adaptability of canonical C-NHEJ (Ku- and Xrcc4/ligase IV-dependent) to imperfect complementary ends; and 3) the requirement to first process chemically incompatible DNA ends that cannot be ligated directly. Thus, C-NHEJ is conservative but adaptable, and the accuracy of the repair is dictated by the structure of the DNA ends rather than by the C-NHEJ machinery. We present data from different organisms that describe the conservative/versatile properties of C-NHEJ. The advantages of the adaptability/versatility of C-NHEJ are discussed for the development of the immune repertoire and the resistance to ionizing radiation, especially at low doses, and for targeted genome manipulation.

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.

PiggyMac, a domesticated piggyBac transposase involved in programmed genome rearrangements in the ciliate Paramecium tetraurelia.

Programmed genome rearrangements drive functional gene assembly in ciliates during the development of the somatic macronucleus. The elimination of germline sequences is directed by noncoding RNAs and is initiated by DNA double-strand breaks, but the enzymes responsible for DNA cleavage have not been identified. We show here that PiggyMac (Pgm), a domesticated piggyBac transposase, is required for these rearrangements in Paramecium tetraurelia. A GFP-Pgm fusion localizes in developing macronuclei, where rearrangements take place, and RNAi-mediated silencing of PGM abolishes DNA cleavage. This is the first in vivo evidence suggesting an essential endonucleolytic function of a domesticated piggyBac transposase.

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.

Translational control of intron splicing in eukaryotes.

Most eukaryotic genes are interrupted by non-coding introns that must be accurately removed from pre-messenger RNAs to produce translatable mRNAs. Splicing is guided locally by short conserved sequences, but genes typically contain many potential splice sites, and the mechanisms specifying the correct sites remain poorly understood. In most organisms, short introns recognized by the intron definition mechanism cannot be efficiently predicted solely on the basis of sequence motifs. In multicellular eukaryotes, long introns are recognized through exon definition and most genes produce multiple mRNA variants through alternative splicing. The nonsense-mediated mRNA decay (NMD) pathway may further shape the observed sets of variants by selectively degrading those containing premature termination codons, which are frequently produced in mammals. Here we show that the tiny introns of the ciliate Paramecium tetraurelia are under strong selective pressure to cause premature termination of mRNA translation in the event of intron retention, and that the same bias is observed among the short introns of plants, fungi and animals. By knocking down the two P. tetraurelia genes encoding UPF1, a protein that is crucial in NMD, we show that the intrinsic efficiency of splicing varies widely among introns and that NMD activity can significantly reduce the fraction of unspliced mRNAs. The results suggest that, independently of alternative splicing, species with large intron numbers universally rely on NMD to compensate for suboptimal splicing efficiency and accuracy.

Highly precise and developmentally programmed genome assembly in Paramecium requires ligase IV-dependent end joining.

During the sexual cycle of the ciliate Paramecium, assembly of the somatic genome includes the precise excision of tens of thousands of short, non-coding germline sequences (Internal Eliminated Sequences or IESs), each one flanked by two TA dinucleotides. It has been reported previously that these genome rearrangements are initiated by the introduction of developmentally programmed DNA double-strand breaks (DSBs), which depend on the domesticated transposase PiggyMac. These DSBs all exhibit a characteristic geometry, with 4-base 5' overhangs centered on the conserved TA, and may readily align and undergo ligation with minimal processing. However, the molecular steps and actors involved in the final and precise assembly of somatic genes have remained unknown. We demonstrate here that Ligase IV and Xrcc4p, core components of the non-homologous end-joining pathway (NHEJ), are required both for the repair of IES excision sites and for the circularization of excised IESs. The transcription of LIG4 and XRCC4 is induced early during the sexual cycle and a Lig4p-GFP fusion protein accumulates in the developing somatic nucleus by the time IES excision takes place. RNAi-mediated silencing of either gene results in the persistence of free broken DNA ends, apparently protected against extensive resection. At the nucleotide level, controlled removal of the 5'-terminal nucleotide occurs normally in LIG4-silenced cells, while nucleotide addition to the 3' ends of the breaks is blocked, together with the final joining step, indicative of a coupling between NHEJ polymerase and ligase activities. Taken together, our data indicate that IES excision is a "cut-and-close" mechanism, which involves the introduction of initiating double-strand cleavages at both ends of each IES, followed by DSB repair via highly precise end joining. This work broadens our current view on how the cellular NHEJ pathway has cooperated with domesticated transposases for the emergence of new mechanisms involved in genome dynamics.

Uncoupling programmed DNA cleavage and repair scrambles the Paramecium somatic genome.

In the ciliate Paramecium, precise excision of numerous internal eliminated sequences (IESs) from the somatic genome is essential at each sexual cycle. DNA double-strands breaks (DSBs) introduced by the PiggyMac endonuclease are repaired in a highly concerted manner by the non-homologous end joining (NHEJ) pathway, illustrated by complete inhibition of DNA cleavage when Ku70/80 proteins are missing. We show that expression of a DNA-binding-deficient Ku70 mutant (Ku70-6E) permits DNA cleavage but leads to the accumulation of unrepaired DSBs. We uncoupled DNA cleavage and repair by co-expressing wild-type and mutant Ku70. High-throughput sequencing of the developing macronucleus genome in these conditions identifies the presence of extremities healed by de novo telomere addition and numerous translocations between IES-flanking sequences. Coupling the two steps of IES excision ensures that both extremities are held together throughout the process, suggesting that DSB repair proteins are essential for assembly of a synaptic precleavage complex.

Functional specialization of Piwi proteins in Paramecium tetraurelia from post-transcriptional gene silencing to genome remodelling.

Proteins of the Argonaute family are small RNA carriers that guide regulatory complexes to their targets. The family comprises two major subclades. Members of the Ago subclade, which are present in most eukaryotic phyla, bind different classes of small RNAs and regulate gene expression at both transcriptional and post-transcriptional levels. Piwi subclade members appear to have been lost in plants and fungi and were mostly studied in metazoa, where they bind piRNAs and have essential roles in sexual reproduction. Their presence in ciliates, unicellular organisms harbouring both germline micronuclei and somatic macronuclei, offers an interesting perspective on the evolution of their functions. Here, we report phylogenetic and functional analyses of the 15 Piwi genes from Paramecium tetraurelia. We show that four constitutively expressed proteins are involved in siRNA pathways that mediate gene silencing throughout the life cycle. Two other proteins, specifically expressed during meiosis, are required for accumulation of scnRNAs during sexual reproduction and for programmed genome rearrangements during development of the somatic macronucleus. Our results indicate that Paramecium Piwi proteins have evolved to perform both vegetative and sexual functions through mechanisms ranging from post-transcriptional mRNA cleavage to epigenetic regulation of genome rearrangements.

Programmed chromosome fragmentation in ciliated protozoa: multiple means to chromosome ends.

SUMMARYCiliated protozoa undergo large-scale developmental rearrangement of their somatic genomes when forming a new transcriptionally active macronucleus during conjugation. This process includes the fragmentation of chromosomes derived from the germline, coupled with the efficient healing of the broken ends by de novo telomere addition. Here, we review what is known of developmental chromosome fragmentation in ciliates that have been well-studied at the molecular level (Tetrahymena, Paramecium, Euplotes, Stylonychia, and Oxytricha). These organisms differ substantially in the fidelity and precision of their fragmentation systems, as well as in the presence or absence of well-defined sequence elements that direct excision, suggesting that chromosome fragmentation systems have evolved multiple times and/or have been significantly altered during ciliate evolution. We propose a two-stage model for the evolution of the current ciliate systems, with both stages involving repetitive or transposable elements in the genome. The ancestral form of chromosome fragmentation is proposed to have been derived from the ciliate small RNA/chromatin modification process that removes transposons and other repetitive elements from the macronuclear genome during development. The evolution of this ancestral system is suggested to have potentiated its replacement in some ciliate lineages by subsequent fragmentation systems derived from mobile genetic elements.

Global trends of whole-genome duplications revealed by the ciliate Paramecium tetraurelia.

The duplication of entire genomes has long been recognized as having great potential for evolutionary novelties, but the mechanisms underlying their resolution through gene loss are poorly understood. Here we show that in the unicellular eukaryote Paramecium tetraurelia, a ciliate, most of the nearly 40,000 genes arose through at least three successive whole-genome duplications. Phylogenetic analysis indicates that the most recent duplication coincides with an explosion of speciation events that gave rise to the P. aurelia complex of 15 sibling species. We observed that gene loss occurs over a long timescale, not as an initial massive event. Genes from the same metabolic pathway or protein complex have common patterns of gene loss, and highly expressed genes are over-retained after all duplications. The conclusion of this analysis is that many genes are maintained after whole-genome duplication not because of functional innovation but because of gene dosage constraints.

Functional study of genes essential for autogamy and nuclear reorganization in Paramecium.

Like all ciliates, Paramecium tetraurelia is a unicellular eukaryote that harbors two kinds of nuclei within its cytoplasm. At each sexual cycle, a new somatic macronucleus (MAC) develops from the germ line micronucleus (MIC) through a sequence of complex events, which includes meiosis, karyogamy, and assembly of the MAC genome from MIC sequences. The latter process involves developmentally programmed genome rearrangements controlled by noncoding RNAs and a specialized RNA interference machinery. We describe our first attempts to identify genes and biological processes that contribute to the progression of the sexual cycle. Given the high percentage of unknown genes annotated in the P. tetraurelia genome, we applied a global strategy to monitor gene expression profiles during autogamy, a self-fertilization process. We focused this pilot study on the genes carried by the largest somatic chromosome and designed dedicated DNA arrays covering 484 genes from this chromosome (1.2% of all genes annotated in the genome). Transcriptome analysis revealed four major patterns of gene expression, including two successive waves of gene induction. Functional analysis of 15 upregulated genes revealed four that are essential for vegetative growth, one of which is involved in the maintenance of MAC integrity and another in cell division or membrane trafficking. Two additional genes, encoding a MIC-specific protein and a putative RNA helicase localizing to the old and then to the new MAC, are specifically required during sexual processes. Our work provides a proof of principle that genes essential for meiosis and nuclear reorganization can be uncovered following genome-wide transcriptome analysis.

The unusual structure of the PiggyMac cysteine-rich domain reveals zinc finger diversity in PiggyBac-related transposases.

BACKGROUND: Transposons are mobile genetic elements that colonize genomes and drive their plasticity in all organisms. DNA transposon-encoded transposases bind to the ends of their cognate transposons and catalyze their movement. In some cases, exaptation of transposon genes has allowed novel cellular functions to emerge. The PiggyMac (Pgm) endonuclease of the ciliate Paramecium tetraurelia is a domesticated transposase from the PiggyBac family. It carries a core catalytic domain typical of PiggyBac-related transposases and a short cysteine-rich domain (CRD), flanked by N- and C-terminal extensions. During sexual processes Pgm catalyzes programmed genome rearrangements (PGR) that eliminate ~ 30% of germline DNA from the somatic genome at each generation. How Pgm recognizes its DNA cleavage sites in chromatin is unclear and the structure-function relationships of its different domains have remained elusive. RESULTS: We provide insight into Pgm structure by determining the fold adopted by its CRD, an essential domain required for PGR. Using Nuclear Magnetic Resonance, we show that the Pgm CRD binds two Zn2+ ions and forms an unusual binuclear cross-brace zinc finger, with a circularly permutated treble-clef fold flanked by two flexible arms. The Pgm CRD structure clearly differs from that of several other PiggyBac-related transposases, among which is the well-studied PB transposase from Trichoplusia ni. Instead, the arrangement of cysteines and histidines in the primary sequence of the Pgm CRD resembles that of active transposases from piggyBac-like elements found in other species and of human PiggyBac-derived domesticated transposases. We show that, unlike the PB CRD, the Pgm CRD does not bind DNA. Instead, it interacts weakly with the N-terminus of histone H3, whatever its lysine methylation state. CONCLUSIONS: The present study points to the structural diversity of the CRD among transposases from the PiggyBac family and their domesticated derivatives, and highlights the diverse interactions this domain may establish with chromatin, from sequence-specific DNA binding to contacts with histone tails. Our data suggest that the Pgm CRD fold, whose unusual arrangement of cysteines and histidines is found in all PiggyBac-related domesticated transposases from Paramecium and Tetrahymena, was already present in the ancestral active transposase that gave rise to ciliate domesticated proteins.