Dynamics of Gene Loss following Ancient Whole-Genome Duplication in the Cryptic Paramecium Complex.

Whole-genome duplications (WGDs) have shaped the gene repertoire of many eukaryotic lineages. The redundancy created by WGDs typically results in a phase of massive gene loss. However, some WGD-derived paralogs are maintained over long evolutionary periods, and the relative contributions of different selective pressures to their maintenance are still debated. Previous studies have revealed a history of three successive WGDs in the lineage of the ciliate Paramecium tetraurelia and two of its sister species from the Paramecium aurelia complex. Here, we report the genome sequence and analysis of 10 additional P. aurelia species and 1 additional out group, revealing aspects of post-WGD evolution in 13 species sharing a common ancestral WGD. Contrary to the morphological radiation of vertebrates that putatively followed two WGD events, members of the cryptic P. aurelia complex have remained morphologically indistinguishable after hundreds of millions of years. Biases in gene retention compatible with dosage constraints appear to play a major role opposing post-WGD gene loss across all 13 species. In addition, post-WGD gene loss has been slower in Paramecium than in other species having experienced genome duplication, suggesting that the selective pressures against post-WGD gene loss are especially strong in Paramecium. A near complete lack of recent single-gene duplications in Paramecium provides additional evidence for strong selective pressures against gene dosage changes. This exceptional data set of 13 species sharing an ancestral WGD and 2 closely related out group species will be a useful resource for future studies on Paramecium as a major model organism in the evolutionary cell biology.

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.

Genome-defence small RNAs exapted for epigenetic mating-type inheritance.

In the ciliate Paramecium, transposable elements and their single-copy remnants are deleted during the development of somatic macronuclei from germline micronuclei, at each sexual generation. Deletions are targeted by scnRNAs, small RNAs produced from the germ line during meiosis that first scan the maternal macronuclear genome to identify missing sequences, and then allow the zygotic macronucleus to reproduce the same deletions. Here we show that this process accounts for the maternal inheritance of mating types in Paramecium tetraurelia, a long-standing problem in epigenetics. Mating type E depends on expression of the transmembrane protein mtA, and the default type O is determined during development by scnRNA-dependent excision of the mtA promoter. In the sibling species Paramecium septaurelia, mating type O is determined by coding-sequence deletions in a different gene, mtB, which is specifically required for mtA expression. These independently evolved mechanisms suggest frequent exaptation of the scnRNA pathway to regulate cellular genes and mediate transgenerational epigenetic inheritance of essential phenotypic polymorphisms.

A mating-type mutagenesis screen identifies a zinc-finger protein required for specific DNA excision events in Paramecium.

In the ciliate Paramecium tetraurelia, functional genes are reconstituted during development of the somatic macronucleus through the precise excision of ∼45 000 single-copy Internal Eliminated Sequences (IESs), thought to be the degenerate remnants of ancient transposon insertions. Like introns, IESs are marked only by a weak consensus at their ends. How such a diverse set of sequences is faithfully recognized and precisely excised remains unclear: specialized small RNAs have been implicated, but in their absence up to ∼60% of IESs are still correctly excised. To get further insight, we designed a mutagenesis screen based on the hypersensitivity of a specific excision event in the mtA gene, which determines mating types. Unlike most IES-containing genes, the active form of mtA is the unexcised one, allowing the recovery of hypomorphic alleles of essential IES recognition/excision factors. Such is the case of one mutation recovered in the Piwi gene PTIWI09, a key player in small RNA-mediated IES recognition. Another mutation identified a novel protein with a C2H2 zinc finger, mtGa, which is required for excision of a small subset of IESs characterized by enrichment in a 5-bp motif. The unexpected implication of a sequence-specific factor establishes a new paradigm for IES recognition and/or excision.

Loss of a Fragile Chromosome Region leads to the Screwy Phenotype in Paramecium tetraurelia.

A conspicuous cell-shape phenotype known as "screwy" was reported to result from mutations at two or three uncharacterized loci in the ciliate Paramecium tetraurelia. Here, we describe a new screwy mutation, Spinning Top, which appeared spontaneously in the cross of an unrelated mutant with reference strain 51. The macronuclear (MAC) genome of the Spinning Top mutant is shown to lack a ~28.5-kb segment containing 18 genes at the end of one chromosome, which appears to result from a collinear deletion in the micronuclear (MIC) genome. We tested several candidate genes from the deleted locus by dsRNA-induced silencing in wild-type cells, and identified a single gene responsible for the phenotype. This gene, named Spade, encodes a 566-aa glutamine-rich protein with a C2HC zinc finger. Its silencing leads to a fast phenotype switch during vegetative growth, but cells recover a wild-type phenotype only 5-6 divisions after silencing is stopped. We analyzed 5 independently-obtained mutant alleles of the Sc1 locus, and concluded that all of them also lack the Spade gene and a number of neighboring genes in the MAC and MIC genomes. Mapping of the MAC deletion breakpoints revealed two different positions among the 5 alleles, both of which differ from the Spinning Top breakpoint. These results suggest that this MIC chromosome region is intrinsically unstable in strain 51.

Six domesticated PiggyBac transposases together carry out programmed DNA elimination in Paramecium.

The domestication of transposable elements has repeatedly occurred during evolution and domesticated transposases have often been implicated in programmed genome rearrangements, as remarkably illustrated in ciliates. In Paramecium, PiggyMac (Pgm), a domesticated PiggyBac transposase, carries out developmentally programmed DNA elimination, including the precise excision of tens of thousands of gene-interrupting germline Internal Eliminated Sequences (IESs). Here, we report the discovery of five groups of distant Pgm-like proteins (PgmLs), all able to interact with Pgm and essential for its nuclear localization and IES excision genome-wide. Unlike Pgm, PgmLs lack a conserved catalytic site, suggesting that they rather have an architectural function within a multi-component excision complex embedding Pgm. PgmL depletion can increase erroneous targeting of residual Pgm-mediated DNA cleavage, indicating that PgmLs contribute to accurately position the complex on IES ends. DNA rearrangements in Paramecium constitute a rare example of a biological process jointly managed by six distinct domesticated transposases.

Functional analysis of cauliflower mosaic virus 35S promoter: re-evaluation of the role of subdomains B5, B4 and B2 in promoter activity.

The cauliflower mosaic virus 35S (35S) promoter is used extensively for transgene expression in plants. The promoter has been delineated into different subdomains based on deletion analysis and gain-of-function studies. However, cis-elements important for promoter activity have been identified only in the domains B1 (as-2 element), A1 (as-1 element) and minimal promoter (TATA box). No cis-elements have been described in subdomains B2-B5, although these are reported to be important for the overall activity of the 35S promoter. We have re-evaluated the contribution of three of these subdomains, namely B5, B4 and B2, to 35S promoter activity by developing several modified promoters. The analysis of beta-glucuronidase gene expression driven by the modified promoters in different tissues of primary transgenic tobacco lines, as well as in seedlings of the T(1) generation, revealed new facets about the functional organization of the 35S promoter. This study suggests that: (i) the 35S promoter truncated up to -301 functions in a similar manner to the -343 (full-length) 35S promoter; (ii) the Dof core and I-box core observed in the subdomain B4 are important for 35S promoter activity; and (iii) the subdomain B2 is essential for maintaining an appropriate distance between the proximal and distal regions of the 35S promoter. These observations will aid in the development of functional synthetic 35S promoters with decreased sequence homology. Such promoters can be used to drive multiple transgenes without evoking promoter homology-based gene silencing when attempting gene stacking.

Analysis of promoter activity in transgenic plants by normalizing expression with a reference gene: anomalies due to the influence of the test promoter on the reference promoter.

Variations in transgene expression due to position effect and copy number are normalized when analysing and comparing the strengths of different promoters. In such experiments, the promoter to be tested is placed upstream to a reporter gene and a second expression cassette is introduced in a linked fashion in the same transfer DNA (T-DNA). Normalization in the activity of the test promoter is carried out by calculating the ratio of activities of the test and reference promoters. When an appropriate number of independent transgenic events are analysed, normalization facilitates assessment of the relative strengths of the test promoters being compared. In this study, using different modified versions of the Cauliflower Mosaic Virus (CaMV) 35S promoter expressing the reporter gene beta-glucuronidase (gus) (test cassette) linked to a chloramphenicol acetyl transferase (cat) gene under the wild-type 35S promoter (reference cassette) in transgenic tobacco lines, we observed that cat gene expression varied depending upon the strength of the modified 35S promoter expressing the gus gene. The 35S promoter in the reference cassette was found to have been upregulated in cases where the modified 35S promoter was weaker than the wild-type 35S promoter. Many studies have been carried out in different organisms to study the phenomenon of transcriptional interference, which refers to the reduced expression of the downstream promoter by a closely linked upstream promoter. However, we observed a positive interaction wherein the weakened activity of a promoter led to upregulation of a contiguous promoter. These observations suggest that, in situations where the promoters of the test and reference gene share the same transcription factors, the activity of the test promoter can influence the activity of the reference promoter in a way that the test promoter's strength is underestimated when normalized by the reference promoter.

Strategies for development of functionally equivalent promoters with minimum sequence homology for transgene expression in plants: cis-elements in a novel DNA context versus domain swapping.

The cauliflower mosaic virus 35S (35S) promoter has been extensively used for the constitutive expression of transgenes in dicotyledonous plants. The repetitive use of the same promoter is known to induce transgene inactivation due to promoter homology. As a way to circumvent this problem, we tested two different strategies for the development of synthetic promoters that are functionally equivalent but have a minimum sequence homology. Such promoters can be generated by (a) introducing known cis-elements in a novel or synthetic stretch of DNA or (b) "domain swapping," wherein domains of one promoter can be replaced with functionally equivalent domains from other heterologous promoters. We evaluated the two strategies for promoter modifications using domain A (consisting of minimal promoter and subdomain A1) of the 35S promoter as a model. A set of modified 35S promoters were developed whose strength was compared with the 35S promoter per se using beta-glucuronidase as the reporter gene. Analysis of the expression of the reporter gene in transient assay system showed that domain swapping led to a significant fall in promoter activity. In contrast, promoters developed by placing cis-elements in a novel DNA context showed levels of expression comparable with that of the 35S. Two promoter constructs Mod2A1T and Mod3A1T were then designed by placing the core sequences of minimal promoter and subdomain A1 in divergent DNA sequences. Transgenics developed in tobacco (Nicotiana tabacum) with the two constructs and with 35S as control were used to assess the promoter activity in different tissues of primary transformants. Mod2A1T and Mod3A1T were found to be active in all of the tissues tested, at levels comparable with that of 35S. Further, the expression of the Mod2A1T promoter in the seedlings of the T1 generation was also similar to that of the 35S promoter. The present strategy opens up the possibility of creating a set of synthetic promoters with minimum sequence homology and with expression levels comparable with the wild-type prototype by modifying sequences present between cis-elements for transgene expression in plants.