Macronucleus Research Articles

Two paralogous PHD finger proteins participate in natural genome editing in Paramecium tetraurelia.

The unicellular eukaryote Paramecium tetraurelia contains functionally distinct nuclei: germline micronuclei (MICs) and a somatic macronucleus (MAC). During sex, the MIC genome is reorganized into a new MAC genome and the old MAC is lost. Almost 45,000 unique internal eliminated sequences (IESs) distributed throughout the genome require precise excision to guarantee a functional new MAC genome. Here, we characterize a pair of paralogous PHD finger proteins involved in DNA elimination. DevPF1, the early-expressed paralog, is present in only some of the gametic and post-zygotic nuclei during meiosis. Both DevPF1 and DevPF2 localize in the new developing MACs, where IES excision occurs. Upon DevPF2 knockdown (KD), long IESs are preferentially retained and late-expressed small RNAs decrease; no length preference for retained IESs was observed in DevPF1-KD and development-specific small RNAs were abolished. The expression of at least two genes from the new MAC with roles in genome reorganization seems to be influenced by DevPF1- and DevPF2-KD. Thus, both PHD fingers are crucial for new MAC genome development, with distinct functions, potentially via regulation of non-coding and coding transcription in the MICs and new MACs.

RebL1 is required for macronuclear structure stability and gametogenesis in Tetrahymena thermophila

Histone modification and nucleosome assembly play important roles in chromatin-related processes. Histone chaperones form different complexes and coordinate histone transportation and assembly. Various histone chaperone complexes have been identified in different organisms. The ciliate protozoa (ciliates) have various chromatin structures and different nuclear morphology. However, histone chaperone components and functions of different subunits remain unclear in ciliates. Tetrahymema thermophila contains a transcriptionally active macronucleus (MAC) and a transcriptionally inactive micronucleus (MIC) which exhibit multiple replication and various chromatin remodeling progresses during vegetative growth and sexual developmental stages. Here, we found histone chaperone RebL1 not only localized evenly in the transcriptionally active MAC but also dynamically changed in the MIC during vegetative growth and sexual developmental stages. REBL1 knockdown inhibited cellular proliferation. The macronuclear morphology became bigger in growing mutants. The abnormal macronuclear structure also occurred in the starvation stage. Furthermore, micronuclear meiosis was disturbed during sexual development, leading to a failure to generate new gametic nuclei. RebL1 potentially interacted with various factors involved in histone-modifying complexes and chromatin remodeling complexes in different developmental stages. REBL1 knockdown affected expression levels of the genes involved in chromatin organization and transcription. Taken together, RebL1 plays a vital role in maintaining macronuclear structure stability and gametogenesis in T. thermophila.

Mismatch Repair Protein Msh6Tt Is Necessary for Nuclear Division and Gametogenesis in Tetrahymena thermophila.

DNA mismatch repair (MMR) improves replication accuracy by up to three orders of magnitude. The MutS protein in E. coli or its eukaryotic homolog, the MutSα (Msh2-Msh6) complex, recognizes base mismatches and initiates the mismatch repair mechanism. Msh6 is an essential protein for assembling the heterodimeric complex. However, the function of the Msh6 subunit remains elusive. Tetrahymena undergoes multiple DNA replication and nuclear division processes, including mitosis, amitosis, and meiosis. Here, we found that Msh6Tt localized in the macronucleus (MAC) and the micronucleus (MIC) during the vegetative growth stage and starvation. During the conjugation stage, Msh6Tt only localized in MICs and newly developing MACs. MSH6Tt knockout led to aberrant nuclear division during vegetative growth. The MSH6TtKO mutants were resistant to treatment with the DNA alkylating agent methyl methanesulfonate (MMS) compared to wild type cells. MSH6Tt knockout affected micronuclear meiosis and gametogenesis during the conjugation stage. Furthermore, Msh6Tt interacted with Msh2Tt and MMR-independent factors. Downregulation of MSH2Tt expression affected the stability of Msh6Tt. In addition, MSH6Tt knockout led to the upregulated expression of several MSH6Tt homologs at different developmental stages. Msh6Tt is involved in macronuclear amitosis, micronuclear mitosis, micronuclear meiosis, and gametogenesis in Tetrahymena.

Independent and Complementary Functions of Caf1b and Hir1 for Chromatin Assembly in Tetrahymena thermophila.

Histones and DNA associate to form the nucleosomes of eukaryotic chromatin. Chromatin assembly factor 1 (CAF-1) complex and histone regulatory protein A (HIRA) complex mediate replication-couple (RC) and replication-independent (RI) nucleosome assembly, respectively. CHAF1B and HIRA share a similar domain but play different roles in nucleosome assembly by binding to the different interactors. At present, there is limited understanding for the similarities and differences in their respective functions. Tetrahymena thermophila contains transcriptionally active polyploid macronuclei (MAC) and transcriptionally silent diploid micronuclei (MIC). Here, the distribution patterns of Caf1b and Hir1 exhibited both similarities and distinctions. Both proteins localized to the MAC and MIC during growth, and to the MIC during conjugation. However, Hir1 exhibited additional signaling on parental MAC and new MAC during sexual reproduction and displayed a punctate signal on developing anlagen. Caf1b and Hir1 only co-localized in the MIC with Pcna1 during conjugation. Knockdown of CAF1B impeded cellular growth and arrested sexual reproductive development. Loss of HIR1 led to MIC chromosome defects and aborted sexual development. Co-interference of CAF1B and HIR1 led to a more severe phenotype. Moreover, CAF1B knockdown led to the up-regulation of HIR1 expression, while knockdown of HIR1 also led to an increase in CAF1B expression. Furthermore, Caf1b and Hir1 interacted with different interactors. These results showed that CAF-1 and Hir1 have independent and complementary functions for chromatin assembly in T. thermophila.

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.

Nuclear events during conjugation in the poorly studied model ciliate Paramecium jenningsi

Ciliated protists are highly differentiated unicellular eukaryotes that possess special sexual processes (conjugation and autogamy) that rely on their unique nuclear dimorphism, i.e., the presence of both a germline micronucleus (MIC) and a somatic macronucleus (MAC) in a single cell. The sexual processes show a high diversity in different ciliates. To better understand the differentiation and evolution of sexual processes in closely related species, we investigated the nuclear events during conjugation in Paramecium jenningsi strain GZ, a poorly studied close relative of the well-known P. aurelia complex. The main results include: (1) the conjugation process takes about 48 ​h, including three prezygotic divisions (meiosis I, II, and mitosis) and three postzygotic divisions; (2) the MICs are dominated by the “parachute” phase at the prophase of meiosis I; (3) after meiosis II, a variable number of nuclei undergo the third prezygotic division, and the two products near the paroral cone become the genetically identical migratory and stationary pronuclei; (4) the synkaryon divides three times to form the MIC and the MAC anlagen, but only one nuclear product of the first postzygotic division completes the next two divisions; and (5) an extra cell division is required to complete the last step of conjugation, during which two MIC anlagen undergo mitotic division while two MAC anlagen are distributed between the daughter cells without division. The comparison of the nuclear events during conjugation in representatives of the class Oligohymenophorea reveals that usually there are three prezygotic divisions and a variable number of postzygotic divisions. However, the number of selected and differentiated nuclei after each division differs among species. This study provided a new model ciliate for further investigations of nuclear selection and differentiation as well as nuclear morphology during meiosis.

Mismatch Repair Protein Msh2 Is Necessary for Macronuclear Stability and Micronuclear Division in Tetrahymena thermophila

Mismatch repair (MMR) is a conserved mechanism that is primarily responsible for the repair of DNA mismatches during DNA replication. Msh2 forms MutS heterodimer complexes that initiate the MMR in eukaryotes. The function of Msh2 is less clear under different chromatin structures. Tetrahymena thermophila contains a transcriptionally active macronucleus (MAC) and a transcriptionally silent micronucleus (MIC) in the same cytoplasm. Msh2 is localized in the MAC and MIC during vegetative growth. Msh2 is localized in the perinuclear region around the MIC and forms a spindle-like structure as the MIC divides. During the early conjugation stage, Msh2 is localized in the MIC and disappears from the parental MAC. Msh2 is localized in the new MAC and new MIC during the late conjugation stage. Msh2 also forms a spindle-like structure with a meiotic MIC and mitotic gametic nucleus. MSH2 knockdown inhibits the division of MAC and MIC during vegetative growth and affects cellular proliferation. MSH2 knockdown mutants are sensitive to cisplatin treatment. MSH2 knockdown also affects micronuclear meiosis and gametogenesis during sexual development. Furthermore, Msh2 interacts with MMR-dependent and MMR-independent factors. Therefore, Msh2 is necessary for macronuclear stability, as well as micronuclear mitosis and meiosis in Tetrahymena.

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.

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.

Absolute quantification of chromosome copy numbers in the polyploid macronucleus of Tetrahymena thermophila at the single-cell level.

Amitosis is widespread among eukaryotes, but the underlying mechanisms are poorly understood. The polyploid macronucleus (MAC) of unicellular ciliates divides by amitosis, making ciliates a potentially valuable model system to study this process. However, a method to accurately quantify the copy number of MAC chromosomes has not yet been established. Here, we used droplet digital PCR (ddPCR) to quantify the absolute copy number of the MAC chromosomes in Tetrahymena thermophila. We first confirmed that ddPCR is a sensitive and reproducible method to determine accurate chromosome copy numbers at the single-cell level. We then used ddPCR to determine the copy number of different MAC chromosomes by analyzing individual T.thermophila cells in the G1 and the amitotic (AM) phases. The average copy number of MAC chromosomes was 90.9 at G1 phase, approximately half the number at AM phase (189.8). The copy number of each MAC chromosome varied among individual cells in G1 phase and correlated with cell size, suggesting that amitosis accompanied by unequal cytokinesis causes copy number variability. Furthermore, the fact that MAC chromosome copy number is less variable among AM-phase cells suggests that the copy number is standardized by regulating DNA replication. We also demonstrated that copy numbers differ among different MAC chromosomes and that interchromosomal variations in copy number are consistent across individual cells. Our findings demonstrate that ddPCR can be used to model amitosis in T.thermophila and possibly in other ciliates.

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.

Histone Chaperone Nrp1 Mutation Affects the Acetylation of H3K56 in Tetrahymena thermophila.

Histone modification and nucleosome assembly are mainly regulated by various histone-modifying enzymes and chaperones. The roles of histone-modification enzymes have been well analyzed, but the molecular mechanism of histone chaperones in histone modification and nucleosome assembly is incompletely understood. We previously found that the histone chaperone Nrp1 is localized in the micronucleus (MIC) and the macronucleus (MAC) and involved in the chromatin stability and nuclear division of Tetrahymena thermophila. In the present work, we found that truncated C-terminal mutant HA-Nrp1TrC abnormally localizes in the cytoplasm. The truncated-signal-peptide mutants HA-Nrp1TrNLS1 and HA-Nrp1TrNLS2 are localized in the MIC and MAC. Overexpression of Nrp1TrNLS1 inhibited cellular proliferation and disrupted micronuclear mitosis during the vegetative growth stage. During sexual development, Nrp1TrNLS1 overexpression led to abnormal bouquet structures and meiosis arrest. Furthermore, Histone H3 was not transported into the nucleus; instead, it formed an abnormal speckled cytoplastic distribution in the Nrp1TrNLS1 mutants. The acetylation level of H3K56 in the mutants also decreased, leading to significant changes in the transcription of the genome of the Nrp1TrNLS1 mutants. The histone chaperone Nrp1 regulates the H3 nuclear import and acetylation modification of H3K56 and affects chromatin stability and genome transcription in Tetrahymena.

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.

The histone chaperone Nrp1 is required for chromatin stability and nuclear division in Tetrahymena thermophila

Histone chaperones facilitate DNA replication and repair by promoting chromatin assembly, disassembly and histone exchange. Following histones synthesis and nucleosome assembly, the histones undergo posttranslational modification by different enzymes and are deposited onto chromatins by various histone chaperones. In Tetrahymena thermophila, histones from macronucleus (MAC) and micronucleus (MIC) have been comprehensively investigated, but the function of histone chaperones remains unclear. Histone chaperone Nrp1 in Tetrahymena contains four conserved tetratricopepeptide repeat (TPR) domains and one C-terminal nuclear localization signal. TPR2 is typically interrupted by a large acidic motif. Immunofluorescence staining showed that Nrp1 is located in the MAC and MICs, but disappeared in the apoptotic parental MAC and the degraded MICs during the conjugation stage. Nrp1 was also colocalized with α-tubulin around the spindle structure. NRP1 knockdown inhibited cellular proliferation and led to the loss of chromosome, abnormal macronuclear amitosis, and disorganized micronuclear mitosis during the vegetative growth stage. During sexual developmental stage, the gametic nuclei failed to be selected and abnormally degraded in NRP1 knockdown mutants. Affinity purification combined with mass spectrometry analysis indicated that Nrp1 is co-purified with core histones, heat shock proteins, histone chaperones, and DNA damage repair proteins. The physical direct interaction of Nrp1 and Asf1 was also confirmed by pull-down analysis in vitro. The results show that histone chaperone Nrp1 is involved in micronuclear mitosis and macronuclear amitosis in the vegetative growth stage and maintains gametic nuclei formation during the sexual developmental stage. Nrp1 is required for chromatin stability and nuclear division in Tetrahymena thermophila.

New contribution to epigenetic studies: Isolation of micronuclei with high purity and DNA integrity in the model ciliated protist, Tetrahymena thermophila

The ciliated protist Tetrahymena thermophila is a well-known model organism with typical nuclear dimorphism containing a somatic macronucleus (MAC) and a germline micronucleus (MIC). The presence in the same cell compartment of two nuclei with distinctly different structural and functional properties provides an ideal model system to explore mechanisms of genome maintenance. Although methods for the isolation of MIC have been available for many years, cross-contamination and DNA degradation remain unresolved. Here, we describe a reliable and quick method to isolate MIC with high purity and DNA integrity in T. thermophila. Different factors are examined to optimize the MIC purification. The MAC contamination ratio in purified MIC is about 0.19% and DNA integrity of purified MIC is maintained. We also establish a more accurate method to detect the contamination rate of nuclei including microscopic observation and PCR detection. This study will facilitate further epigenetic research in Tetrahymena.

Transcribed germline-limited coding sequences inOxytricha trifallax

The germline-soma divide is a fundamental distinction in developmental biology, and different genes are expressed in germline and somatic cells throughout metazoan life cycles. Ciliates, a group of microbial eukaryotes, exhibit germline-somatic nuclear dimorphism within a single cell with two different genomes. The ciliate Oxytricha trifallax undergoes massive RNA-guided DNA elimination and genome rearrangement to produce a new somatic macronucleus (MAC) from a copy of the germline micronucleus (MIC). This process eliminates noncoding DNA sequences that interrupt genes and also deletes hundreds of germline-limited open reading frames (ORFs) that are transcribed during genome rearrangement. Here, we update the set of transcribed germline-limited ORFs (TGLOs) in O. trifallax. We show that TGLOs tend to be expressed during nuclear development and then are absent from the somatic MAC. We also demonstrate that exposure to synthetic RNA can reprogram TGLO retention in the somatic MAC and that TGLO retention leads to transcription outside the normal developmental program. These data suggest that TGLOs represent a group of developmentally regulated protein-coding sequences whose gene expression is terminated by DNA elimination.

The macronuclear genome of anaerobic ciliate Entodinium caudatum reveals its biological features adapted to the distinct rumen environment

Entodinium caudatum is an anaerobic binucleated ciliate representing the most dominant protozoal species in the rumen. However, its biological features are largely unknown due to the inability to establish an axenic culture. In this study, we primally sequenced its macronucleus (MAC) genome to aid the understanding of its metabolism, physiology, ecology. We isolated the MAC of E. caudatum strain MZG-1 and sequenced the MAC genome using Illumina MiSeq, MinION, and PacBio RSII systems. De novo assembly of the MiSeq sequence reads followed with subsequent scaffolding with MinION and PacBio reads resulted in a draft MAC genome about 117 Mbp. A large number of carbohydrate-active enzymes were likely acquired through horizontal gene transfer. About 8.74% of the E. caudatum predicted proteome was predicted as proteases. The MAC genome of E. caudatum will help better understand its important roles in rumen carbohydrate metabolism, and interaction with other members of the rumen microbiome.

A strategy for complete telomere-to-telomere assembly of ciliate macronuclear genome using ultra-high coverage Nanopore data

Ciliates contain two kinds of nuclei: the germline micronucleus (MIC) and the somatic macronucleus (MAC) in a single cell. The MAC usually have fragmented chromosomes. These fragmented chromosomes, capped with telomeres at both ends, could be gene size to several megabases in length among different ciliate species. So far, no telomere-to-telomere assembly of entire MAC genome in ciliate species has been finished. Development of the third generation sequencing technologies allows to generate sequencing reads up to megabases in length that could possibly span an entire MAC chromosome. Taking advantage of the ultra-long Nanopore reads, we established a simple strategy for the complete assembly of ciliate MAC genomes. Using this strategy, we assembled the complete MAC genomes of two ciliate species Tetrahymena thermophila and Tetrahymena shanghaiensis, composed of 181 and 214 chromosomes telomere-to-telomere respectively. The established strategy as well as the high-quality genome data will provide a useful approach for ciliate genome assembly, and a valuable community resource for further biological, evolutionary and population genomic studies.

Genome plasticity in Paramecium bursaria revealed by population genomics

BackgroundCiliates are an ancient and diverse eukaryotic group found in various environments. A unique feature of ciliates is their nuclear dimorphism, by which two types of nuclei, the diploid germline micronucleus (MIC) and polyploidy somatic macronucleus (MAC), are present in the same cytoplasm and serve different functions. During each sexual cycle, ciliates develop a new macronucleus in which newly fused genomes are extensively rearranged to generate functional minichromosomes. Interestingly, each ciliate species seems to have its way of processing genomes, providing a diversity of resources for studying genome plasticity and its regulation. Here, we sequenced and analyzed the macronuclear genome of different strains of Paramecium bursaria, a highly divergent species of the genus Paramecium which can stably establish endosymbioses with green algae.ResultsWe assembled a high-quality macronuclear genome of P. bursaria and further refined genome annotation by comparing population genomic data. We identified several species-specific expansions in protein families and gene lineages that are potentially associated with endosymbiosis. Moreover, we observed an intensive chromosome breakage pattern that occurred during or shortly after sexual reproduction and contributed to highly variable gene dosage throughout the genome. However, patterns of copy number variation were highly correlated among genetically divergent strains, suggesting that copy number is adjusted by some regulatory mechanisms or natural selection. Further analysis showed that genes with low copy number variation among populations tended to function in basic cellular pathways, whereas highly variable genes were enriched in environmental response pathways.ConclusionsWe report programmed DNA rearrangements in the P. bursaria macronuclear genome that allow cells to adjust gene copy number globally according to individual gene functions. Our results suggest that large-scale gene copy number variation may represent an ancient mechanism for cells to adapt to different environments.

Identification and Characterization of Base-Substitution Mutations in the Macronuclear Genome of the Ciliate Tetrahymena thermophila.

Polyploidy can provide adaptive advantages and drive evolution. Amitotic division of the polyploid macronucleus (MAC) in ciliates acts as a nonsexual genetic mechanism to enhance adaptation to stress conditions and thus provides a unique model to investigate the evolutionary role of polyploidy. Mutation is the primary source of the variation responsible for evolution and adaptation; however, to date, de novo mutations that occur in ciliate MAC genomes during these processes have not been characterized and their biological impacts are undefined. Here, we carried out long-term evolution experiments to directly explore de novo MAC mutations and their molecular features in the model ciliate, Tetrahymena thermophila. A simple but effective method was established to detect base-substitution mutations in evolving populations whereas filtering out most of the false positive base-substitutions caused by repetitive sequences and the programmed genome rearrangements. The detected mutations were rigorously validated using the MassARRAY system. Validated mutations showed a strong G/C→A/T bias, consistent with observations in other species. Moreover, a progressive increase in growth rate of the evolving populations suggested that some of these mutations might be responsible for cell fitness. The established mutation identification and validation methods will be an invaluable resource to make ciliates an important model system to study the role of polyploidy in evolution.