Programmed DNA under-amplification in Paramecium primaurelia

Species-Specific Duplication of Surface Antigen Genes in Paramecium

Paramecium is a free-living ciliate that undergoes antigenic variation and still the functions of these variable surface antigen coats in this non-pathogenic ciliate remain elusive. Only a few surface antigen genes have been described, mainly in the two model species P. tetraurelia strain 51 and P. primaurelia strain 156. Given the lack of suitable sequence data to allow for phylogenetics and deeper sequence comparisons, we screened the genomes of six different Paramecium species for serotype genes and isolated 548 candidates. Our approach identified the subfamilies of the isogenes of individual serotypes that were mostly represented by intrachromosomal gene duplicates. These showed different duplication levels, and chromosome synteny suggested rather young duplication events after the emergence of the P. aurelia species complex, indicating a rapid evolution of surface antigen genes. We were able to identify the different subfamilies of the surface antigen genes with internal tandem repeats, which showed consensus motifs across species. The individual isogene families showed additional consensus motifs, indicating that the selection pressure holds individual amino acids constant in these repeats. This may be a hint of the receptor function of these antigens rather than a presentation of random epitopes, generating the variability of these surface molecules.

Analysis of sequence variability in the macronuclear DNA of Paramecium tetraurelia: a somatic view of the germline

Ciliates are the only unicellular eukaryotes known to separate germinal and somatic functions. Diploid but silent micronuclei transmit the genetic information to the next sexual generation. Polyploid macronuclei express the genetic information from a streamlined version of the genome but are replaced at each sexual generation. The macronuclear genome of Paramecium tetraurelia was recently sequenced by a shotgun approach, providing access to the gene repertoire. The 72-Mb assembly represents a consensus sequence for the somatic DNA, which is produced after sexual events by reproducible rearrangements of the zygotic genome involving elimination of repeated sequences, precise excision of unique-copy internal eliminated sequences (IES), and amplification of the cellular genes to high copy number. We report use of the shotgun sequencing data (>10(6) reads representing 13 x coverage of a completely homozygous clone) to evaluate variability in the somatic DNA produced by these developmental genome rearrangements. Although DNA amplification appears uniform, both of the DNA elimination processes produce sequence heterogeneity. The variability that arises from IES excision allowed identification of hundreds of putative new IESs, compared to 42 that were previously known, and revealed cases of erroneous excision of segments of coding sequences. We demonstrate that IESs in coding regions are under selective pressure to introduce premature termination of translation in case of excision failure.

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.

Genetic structure and evolution of the Vps25 family, a yeast ESCRT-II component

Background

Vps25p is the product of yeast gene VPS25 and is found in an endosomal sorting complex required for transport (ESCRT)-II, along with Vps22p and Vps36p. This complex is essential for sorting of ubiquitinated biosynthetic and endosomal cargoes into endosomes.

Results

We found that VPS25 is a highly conserved and widely expressed eukaryotic gene, with single orthologs in chromalveolate, excavate, amoebozoan, plant, fungal and metazoan species. Two paralogs were found in Trichomonas vaginalis. An ortholog was strikingly absent from the Encephalitozoon cuniculi genome. Intron positions were analyzed in VPS25 from 36 species. We found evidence for five ancestral VPS25 introns, intron loss, and single instances of intron gain (a Paramecium species) and intron slippage (Theileria species). Processed pseudogenes were identified in four mammalian genomes, with a notable absence in the mouse genome. Two retropseudogenes were found in the chimpanzee genome, one more recently inserted, and one evolving from a common primate ancestor. The amino acid sequences of 119 Vps25 orthologs are aligned, compared with the known secondary structure of yeast Vps25p, and used to carry out phylogenetic analysis. Residues in two amino-terminal PPXY motifs (motif I and II), involved in dimerization of Vps25p and interaction with Vps22p and Vps36p, were closely, but not absolutely conserved. Specifically, motif I was absent in Vps25 homologs of chromalveolates, euglenozoa, and diplomonads. A highly conserved carboxy-terminal lysine was identified, which suggests Vps25 is ubiquitinated. Arginine-83 of yeast Vps25p involved in Vps22p interaction was highly, but not absolutely, conserved. Human tissue expression analysis showed universal expression.

Conclusion

We have identified 119 orthologs of yeast Vps25p. Expression of mammalian VPS25 in a wide range of tissues, and the presence in a broad range of eukaryotic species, indicates a basic role in eukaryotic cell function. Intron splice site positions were highly conserved across all major eukaryotic species, suggesting an ancestral origin. Amino acid sequence analysis showed the consensus for the amino-terminal proline-rich motifs is P- [WP]-X-[YF] for motif I (when present) and P-P-[FYL]-[FY] for motif II, and that Vps25 may be ubiquitinated.

Epigenetic control of chromosome breakage at the 5' end of Paramecium tetraurelia gene A

Macronuclei and micronuclei of ciliates have related genomes, with macronuclei developing from zygotic micronuclei through programmed DNA rearrangements. While Paramecium tetraurelia wild-type strain 51 and mutant strain d48 have the same micronuclear genome, qualitative differences between their macronuclear genomes have been described, demonstrating that programmed DNA rearrangements could be epigenetically controlled in ciliates. Macronuclear chromosomes end downstream of gene A (A51 Mac ends) and at the 5' end of gene A (Ad48 Mac ends) in strains 51 and d48, respectively. To gain further insight into the process of chromosome end formation, we performed an extensive analysis of locus A rearrangement in strains d48 and 51, in strain d12, which harbors a gene A deletion, and in interstrain cross progeny. We show that (i) allele Ad12 harbors a deletion of >16 kb, (ii) A51 Mac ends distribute over four rather than three DNA regions, (iii) strains d48 and 51 display only quantitative differences (rare Ad48 and A51 Mac ends do form in strains 51 and d48, respectively), (iv) the level of A51 Mac ends is severalfold enhanced in d12- and d48-derived progeny, and (v) this level inversely correlates with the level of Ad48 Mac ends in the d48 parent. Together, these data lead to a model in which the formation of Ad48 Mac ends is epigenetically controlled by a d48 factor(s). We propose that the d48 factor(s) may be derived from RNA molecules transcribed from the Ad48 Mac ends and encompassing the truncated A gene and telomeric repeats.

High coding density on the largest Paramecium tetraurelia somatic chromosome

Paramecium, like other ciliates, remodels its entire germline genome at each sexual generation to produce a somatic genome stripped of transposons and other multicopy elements. The germline chromosomes are fragmented by a DNA elimination process that targets heterochromatin to give a reproducible set of some 200 linear molecules 50 kb to 1 Mb in size. These chromosomes are maintained at a ploidy of 800n in the somatic macronucleus and assure all gene expression. We isolated and sequenced the largest megabase somatic chromosome in order to explore its organization and gene content. The AT-rich (72%) chromosome is compact, with very small introns (average size 25 nt), short intergenic regions (median size 202 nt), and a coding density of at least 74%, higher than that reported for budding yeast (70%) or any other free-living eukaryote. Similarity to known proteins could be detected for 57% of the 460 potential protein coding genes. Thirty-two of the proteins are shared with vertebrates but absent from yeast, consistent with the morphogenetic complexity of Paramecium, a long-standing model for differentiated functions shared with metazoans but often absent from simpler eukaryotes. Extrapolation to the whole genome suggests that Paramecium has at least 30,000 genes.

PAK paradox: Paramecium appears to have more K(+)-channel genes than humans

K(+)-selective ion channels (K(+) channels) have been found in bacteria, archaea, eucarya, and viruses. In Paramecium and other ciliates, K(+) currents play an essential role in cilia-based motility. We have retrieved and sequenced seven closely related Paramecium K(+)-channel gene (PAK) sequences by using previously reported fragments. An additional eight unique K(+)-channel sequences were retrieved from an indexed library recently used in a pilot genome sequencing project. Alignments of these protein translations indicate that while these 15 genes have diverged at different times, they all maintain many characteristics associated with just one subclass of metazoan K(+) channels (CNG/ERG type). Our results indicate that most of the genes are expressed, because all predicted frameshifts and several gaps in the homolog alignments contain Paramecium intron sequences deleted from reverse transcription-PCR products. Some of the variations in the 15 genomic nucleotide sequences involve an absence of introns, even between very closely related sequences, suggesting a potential occurrence of reverse transcription in the past. Extrapolation from the available genome sequence indicates that Paramecium harbors as many as several hundred of this one type of K(+)-channel gene. This quantity is far more numerous than those of K(+)-channel genes of all types known in any metazoan (e.g., approximately 80 in humans, approximately 30 in flies, and approximately 15 in Arabidopsis). In an effort to understand this plurality, we discuss several possible reasons for their maintenance, including variations in expression levels in response to changes in the freshwater environment, like that seen with other major plasma membrane proteins in Paramecium.

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.

Control of DNA excision efficiency in Paramecium

Programmed excision of internal eliminated sequences (IESs) occurs at thousands of sites in ciliate genomes. How this is controlled is largely unknown. Here, we report the characterization of the non-efficiently excised 156psiG-11 IES from Paramecium primaurelia strain 156 and that of the efficiently excised 168psiG-11 IES, an allelic variant from strain 168. Then, we report a genetic and molecular analysis of IES excision efficiency in F(1) progeny derived from interstrain crosses and in F(2) homozygous progeny derived from F(1) autogamy. IES 168psiG-11 excision efficiency was approximately 100% in all cases. IES 156psiG-11 excision efficiency was 19 +/- 13% in F(1) progeny and 0.6 +/- 1.1% in F(2) progeny. No trans-excision event between IESs 156psiG-11 and 168psiG-11 was detected within the F(1) progeny. These data demonstrate that the excision efficiency of this IES is variable and controlled by a cis-acting element. This should encompass positions 8 and/or 9 of the right IES end, which display allele differences. Finally, the 30-fold stimulation of IES 156psiG-11 excision efficiency within F(1) progeny relative to F(2) progeny demonstrates that Paramecium IES excision efficiency is sensitive either to a conjugation-specific trans-acting factor provided by the zygotic genome, or to homologous chromosome cross-talk.

Developmentally programmed excision of internal DNA sequences in Paramecium aurelia

The development of a new somatic nucleus (macronucleus) during sexual reproduction of the ciliate Paramecium aurelia involves reproducible chromosomal rearrangements that affect the entire germline genome. Macronuclear development can be induced experimentally, which makes P. aurelia an attractive model for the study of the mechanism and the regulation of DNA rearrangements. Two major types of rearrangements have been identified: the fragmentation of the germline chromosomes, followed by the formation of the new macronuclear chromosome ends in association with imprecise DNA elimination, and the precise excision of internal eliminated sequences (IESs). All IESs identified so far are short, A/T rich and non-coding elements. They are flanked by a direct repeat of a 5'-TA-3' dinucleotide, a single copy of which remains at the macronuclear junction after excision. The number of these single-copy sequences has been estimated to be around 60,000 per haploid genome. This review focuses on the current knowledge about the genetic and epigenetic determinants of IES elimination in P. aurelia, the analysis of excision products, and the tightly regulated timing of excision throughout macronuclear development. Several models for the molecular mechanism of IES excision will be discussed in relation to those proposed for DNA elimination in other ciliates.