Evolution and assembly of an extremely scrambled gene

Combined Genome and Transcriptome Analyses of the Ciliate Schmidingerella arcuata (Spirotrichea) Reveal Patterns of DNA Elimination, Scrambling, and Inversion

Schmidingerella arcuata is an ecologically important tintinnid ciliate that has long served as a model species in plankton trophic ecology. We present a partial micronuclear genome and macronuclear transcriptome resource for S. arcuata, acquired using single-cell techniques, and we report on pilot analyses including functional annotation and genome architecture. Our analysis shows major fragmentation, elimination, and scrambling in the micronuclear genome of S. arcuata. This work introduces a new nonmodel genome resource for the study of ciliate ecology and genomic biology and provides a detailed functional counterpart to ecological research on S. arcuata.

Roles of Noncoding RNAs in Ciliate Genome Architecture

Ciliates are an interesting model system for investigating diverse functions of noncoding RNAs, especially in genome defence pathways. During sexual development, the ciliate somatic genome undergoes massive rearrangement and reduction through removal of transposable elements and other repetitive DNA. This is guided by a multitude of noncoding RNAs of different sizes and functions, the extent of which is only recently becoming clear. The genome rearrangement pathways evolved as a defence against parasitic DNA, but interestingly also use the transposable elements and transposases to execute their own removal. Thus, ciliates are also a good model for the coevolution of host and transposable element, and the mutual dependence between the two. In this review, we summarise the genome rearrangement pathways in three diverse species of ciliate, with focus on recent discoveries and the roles of noncoding RNAs.

•

Ciliate genomes undergo massive rearrangement and reduction during development.

•

Transposon elimination is guided by small RNAs and carried out by transposases.

•

New pathways for noncoding RNA production have recently been discovered in ciliates.

•

Diverse ciliate species have different mechanisms for RNA-guided genome remodeling.

Insights into an Extensively Fragmented Eukaryotic Genome: De Novo Genome Sequencing of the Multinuclear Ciliate Uroleptopsis citrina

Ciliated protists are a large group of single-celled eukaryotes with separate germline and somatic nuclei in each cell. The somatic genome is developed from the zygotic nucleus through a series of chromosomal rearrangements, including fragmentation, DNA elimination, de novo telomere addition, and DNA amplification. This unique feature makes them perfect models for research in genome biology and evolution. However, genomic research of ciliates has been limited to a few species, owing to problems with DNA contamination and obstacles in cultivation. Here, we introduce a method combining telomere-primer PCR amplification and high-throughput sequencing, which can reduce DNA contamination and obtain genomic data efficiently. Based on this method, we report a draft somatic genome of a multimacronuclear ciliate, Uroleptopsis citrina. 1) The telomeric sequence in U. citrina is confirmed to be C4A4C4A4C4 by directly blunt-end cloning. 2) Genomic analysis of the resulting chromosomes shows a "one-gene one-chromosome" pattern, with a small number of multiple-gene chromosomes. 3) Amino acid usage is analyzed, and reassignment of stop codons is confirmed. 4) Chromosomal analysis shows an obvious asymmetrical GC skew and high bias between A and T in the subtelomeric regions of the sense-strand, with the detection of an 11-bp high AT motif region in the 3' subtelomeric region. 5) The subtelomeric sequence also has an obvious 40 nt strand oscillation of nucleotide ratio. 6) In the 5' subtelomeric region of the coding strand, the distribution of potential TATA-box regions is illustrated, which accumulate between 30 and 50 nt. This work provides a valuable reference for genomic research and furthers our understanding of the dynamic nature of unicellular eukaryotic genomes.

Twisted Tales: Insights into Genome Diversity of Ciliates Using Single-Cell 'Omics

The emergence of robust single-cell 'omics techniques enables studies of uncultivable species, allowing for the (re)discovery of diverse genomic features. In this study, we combine single-cell genomics and transcriptomics to explore genome evolution in ciliates (a > 1 Gy old clade). Analysis of the data resulting from these single-cell 'omics approaches show: 1) the description of the ciliates in the class Karyorelictea as "primitive" is inaccurate because their somatic macronuclei contain loci of varying copy number (i.e., they have been processed by genome rearrangements from the zygotic nucleus); 2) gene-sized somatic chromosomes exist in the class Litostomatea, consistent with Balbiani's (1890) observation of giant chromosomes in this lineage; and 3) gene scrambling exists in the underexplored Postciliodesmatophora (the classes Heterotrichea and Karyorelictea, abbreviated here as the Po-clade), one of two major clades of ciliates. Together these data highlight the complex evolutionary patterns underlying germline genome architectures in ciliates and provide a basis for further exploration of principles of genome evolution in diverse microbial lineages.

Hyper-eccentric structural genes in the mitochondrial genome of the algal parasite Hemistasia phaeocysticola

Diplonemid mitochondria are considered to have very eccentric structural genes. Coding regions of individual diplonemid mitochondrial genes are fragmented into small pieces and found on different circular DNAs. Short RNAs transcribed from each DNA molecule mature through a unique RNA maturation process involving assembly and three types of RNA editing (i.e., U insertion and A-to-I and C-to-U substitutions), although the molecular mechanism(s) of RNA maturation and the evolutionary history of these eccentric structural genes still remain to be understood. Since the gene fragmentation pattern is generally conserved among the diplonemid species studied to date, it was considered that their structural complexity has plateaued and further gene fragmentation could not occur. Here, we show the mitochondrial gene structure of Hemistasia phaeocysticola, which was recently identified as a member of a novel lineage in diplonemids, by comparison of the mitochondrial DNA sequences with cDNA sequences synthesized from mature mRNA. The genes of H. phaeocysticola are fragmented much more finely than those of other diplonemids studied to date. Furthermore, in addition to all known types of RNA editing, it is suggested that a novel processing step (i.e., secondary RNA insertion) is involved in the RNA maturation in the mitochondria of H. phaeocysticola. Our findings demonstrate the tremendous plasticity of mitochondrial gene structures.

<mds_ies_db>: a database of ciliate genome rearrangements

Ciliated protists exhibit nuclear dimorphism through the presence of somatic macronuclei (MAC) and germline micronuclei (MIC). In some ciliates, DNA from precursor segments in the MIC genome rearranges to form transcriptionally active genes in the mature MAC genome, making these ciliates model organisms to study the process of somatic genome rearrangement. Similar broad scale, somatic rearrangement events occur in many eukaryotic cells and tumors. The <mds_ies_db> (http://oxytricha.princeton.edu/mds_ies_db) is a database of genome recombination and rearrangement annotations, and it provides tools for visualization and comparative analysis of precursor and product genomes. The database currently contains annotations for two completely sequenced ciliate genomes: Oxytricha trifallax and Tetrahymena thermophila.

Analyses of alternatively processed genes in ciliates provide insights into the origins of scrambled genomes and may provide a mechanism for speciation

Chromosome rearrangements occur in a variety of eukaryotic life cycles, including during the development of the somatic macronuclear genome in ciliates. Previous work on the phyllopharyngean ciliate Chilodonella uncinata revealed that macronuclear β-tubulin and protein kinase gene families share alternatively processed germ line segments nested within divergent regions. To study genome evolution in this ciliate further, we characterized two additional alternatively processed gene families from two cryptic species of the ciliate morphospecies C. uncinata: those encoding histidine acid phosphatase protein (Hap) and leishmanolysin family protein (Lei). Analyses of the macronuclear Hap and Lei sequences reveal that each gene family consists of three members in the macronucleus that are marked by identical regions nested among highly divergent regions. Investigation of the micronuclear Hap sequences revealed a complex pattern in which the three macronuclear sequences are derived either from a single micronuclear region or from a combination of this shared region recombined with additional duplicate micronuclear copies of Hap. We propose a model whereby gene scrambling evolves by gene duplication followed by partial and reciprocal degradation of the duplicate sequences. In this model, alternative processing represents an intermediate step in the evolution of scrambled genes. Finally, we speculate on the possible role of genome architecture in speciation in ciliates by describing what might happen if changes in alternatively processed loci occur in subdivided populations.

Genome rearrangements occur in a variety of eukaryotic cells and serve as an important mechanism for generating genomic diversity. The unusual genome architecture of ciliates with separate germline and somatic nuclei in each cell, provides an ideal system to study further principles of genome evolution. Previous analyses revealed complex forms of chromosome rearrangements, including gene scrambling and alternative processing of germ line chromosomes. Here we describe more complex rearrangements between germ line and somatic chromosomes than previously seen in alternatively processed gene families. Drawing on the present and previous findings, we propose a model in which alternative processing of duplicated micronuclear regions represents an intermediate stage in the evolution of scrambled genes. Under this model, alternative processing may provide insights into a mechanism for speciation in ciliates. Our data on gene scrambling and alternative processing also enhance views on the dynamic nature of genomes across the eukaryotic tree of life.

RNA-Mediated Epigenetic Programming of Genome Rearrangements

RNA, normally thought of as a conduit in gene expression, has a novel mode of action in ciliated protozoa. Maternal RNA templates provide both an organizing guide for DNA rearrangements and a template that can transport somatic mutations to the next generation. This opportunity for RNA-mediated genome rearrangement and DNA repair is profound in the ciliate Oxytricha, which deletes 95% of its germline genome during development in a process that severely fragments its chromosomes and then sorts and reorders the hundreds of thousands of pieces remaining. Oxytricha's somatic nuclear genome is therefore an epigenome formed through RNA templates and signals arising from the previous generation. Furthermore, this mechanism of RNA-mediated epigenetic inheritance can function across multiple generations, and the discovery of maternal template RNA molecules has revealed new biological roles for RNA and has hinted at the power of RNA molecules to sculpt genomic information in cells.

RNA-dependent control of gene amplification

We exploit the unusual genome organization of the ciliate cell to analyze the control of specific gene amplification during a nuclear differentiation process. Ciliates contain two types of nuclei within one cell, the macronucleus and the micronucleus; and after sexual reproduction a new macronucleus is formed from a micronuclear derivative. During macronuclear differentiation, most extensive DNA reorganization, elimination, and fragmentation processes occur, resulting in a macronucleus containing short DNA molecules (nanochromosomes) representing individual genetic units and each being present in high copy number. It is believed that these processes are controlled by small nuclear RNAs but also by a template derived from the old macronucleus. We first describe the exact copy numbers of selected nanochromosomes in the macronucleus, and define the timing during nuclear differentiation at which copy number is determined. This led to the suggestion that DNA processing and copy number control may be closely related mechanisms. Degradation of an RNA template derived from the macronucleus leads to significant decrease in copy number, whereas injection of additional template molecules results in an increase in copy number and enhanced expression of the corresponding gene. These observations can be incorporated into a mechanistic model about an RNA-dependent epigenetic regulation of gene copy number during nuclear differentiation. This highlights that RNA, in addition to its well-known biological functions, can also be involved in the control of gene amplification.

Systematically fragmented genes in a multipartite mitochondrial genome

Arguably, the most bizarre mitochondrial DNA (mtDNA) is that of the euglenozoan eukaryote Diplonema papillatum. The genome consists of numerous small circular chromosomes none of which appears to encode a complete gene. For instance, the cox1 coding sequence is spread out over nine different chromosomes in non-overlapping pieces (modules), which are transcribed separately and joined to a contiguous mRNA by trans-splicing. Here, we examine how many genes are encoded by Diplonema mtDNA and whether all are fragmented and their transcripts trans-spliced. Module identification is challenging due to the sequence divergence of Diplonema mitochondrial genes. By employing most sensitive protein profile search algorithms and comparing genomic with cDNA sequence, we recognize a total of 11 typical mitochondrial genes. The 10 protein-coding genes are systematically chopped up into three to 12 modules of 60–350 bp length. The corresponding mRNAs are all trans-spliced. Identification of ribosomal RNAs is most difficult. So far, we only detect the 3′-module of the large subunit ribosomal RNA (rRNA); it does not trans-splice with other pieces. The small subunit rRNA gene remains elusive. Our results open new intriguing questions about the biochemistry and evolution of mitochondrial trans-splicing in Diplonema.

The unusual way to make a genetically active nucleus

During macronuclear differentiation in ciliated protozoa, extensive DNA rearrangement and DNA excision processes occur, and these are most profound in stichotrichous ciliates, such as Stylonychia or Oxytricha. This review describes the morphological and molecular events taking place during macronuclear development in stichotrichous ciliates. Various models for the regulation of macronuclear differentiation have been proposed and will be discussed here. Finally, an attempt to speculate about the biological consequences of these rearrangement and excision processes will be made. Because specific elimination of DNA sequences not required in the differentiated nucleus can be regarded as the most extreme form of gene silencing, results obtained in these cells may also be relevant for our understanding of differentiation processes in higher eukaryotic organisms.

Macronuclear genome structure of the ciliate Nyctotherus ovalis: single-gene chromosomes and tiny introns

BACKGROUND

Nyctotherus ovalis is a single-celled eukaryote that has hydrogen-producing mitochondria and lives in the hindgut of cockroaches. Like all members of the ciliate taxon, it has two types of nuclei, a micronucleus and a macronucleus. N. ovalis generates its macronuclear chromosomes by forming polytene chromosomes that subsequently develop into macronuclear chromosomes by DNA elimination and rearrangement.

RESULTS

We examined the structure of these gene-sized macronuclear chromosomes in N. ovalis. We determined the telomeres, subtelomeric regions, UTRs, coding regions and introns by sequencing a large set of macronuclear DNA sequences (4,242) and cDNAs (5,484) and comparing them with each other. The telomeres consist of repeats CCC(AAAACCCC)n, similar to those in spirotrichous ciliates such as Euplotes, Sterkiella (Oxytricha) and Stylonychia. Per sequenced chromosome we found evidence for either a single protein-coding gene, a single tRNA, or the complete ribosomal RNAs cluster. Hence the chromosomes appear to encode single transcripts. In the short subtelomeric regions we identified a few overrepresented motifs that could be involved in gene regulation, but there is no consensus polyadenylation site. The introns are short (21-29 nucleotides), and a significant fraction (1/3) of the tiny introns is conserved in the distantly related ciliate Paramecium tetraurelia. As has been observed in P. tetraurelia, the N. ovalis introns tend to contain in-frame stop codons or have a length that is not dividable by three. This pattern causes premature termination of mRNA translation in the event of intron retention, and potentially degradation of unspliced mRNAs by the nonsense-mediated mRNA decay pathway.

CONCLUSION

The combination of short leaders, tiny introns and single genes leads to very minimal macronuclear chromosomes. The smallest we identified contained only 150 nucleotides.

The pathway to detangle a scrambled gene

Background

Programmed DNA elimination and reorganization frequently occur during cellular differentiation. Development of the somatic macronucleus in some ciliates presents an extreme case, involving excision of internal eliminated sequences (IESs) that interrupt coding DNA segments (macronuclear destined sequences, MDSs), as well as removal of transposon-like elements and extensive genome fragmentation, leading to 98% genome reduction in Stylonychia lemnae. Approximately 20–30% of the genes are estimated to be scrambled in the germline micronucleus, with coding segment order permuted and present in either orientation on micronuclear chromosomes. Massive genome rearrangements are therefore critical for development.

Methodology/Principal Findings

To understand the process of DNA deletion and reorganization during macronuclear development, we examined the population of DNA molecules during assembly of different scrambled genes in two related organisms in a developmental time-course by PCR. The data suggest that removal of conventional IESs usually occurs first, accompanied by a surprising level of error at this step. The complex events of inversion and translocation seem to occur after repair and excision of all conventional IESs and via multiple pathways.

Conclusions/Significance

This study reveals a temporal order of DNA rearrangements during the processing of a scrambled gene, with simpler events usually preceding more complex ones. The surprising observation of a hidden layer of errors, absent from the mature macronucleus but present during development, also underscores the need for repair or screening of incorrectly-assembled DNA molecules.

RNA-mediated epigenetic programming of a genome-rearrangement pathway

Genome-wide DNA rearrangements occur in many eukaryotes but are most exaggerated in ciliates, making them ideal model systems for epigenetic phenomena. During development of the somatic macronucleus, Oxytricha trifallax destroys 95% of its germ line, severely fragmenting its chromosomes, and then unscrambles hundreds of thousands of remaining fragments by permutation or inversion. Here we demonstrate that DNA or RNA templates can orchestrate these genome rearrangements in Oxytricha, supporting an epigenetic model for sequence-dependent comparison between germline and somatic genomes. A complete RNA cache of the maternal somatic genome may be available at a specific stage during development to provide a template for correct and precise DNA rearrangement. We show the existence of maternal RNA templates that could guide DNA assembly, and that disruption of specific RNA molecules disables rearrangement of the corresponding gene. Injection of artificial templates reprogrammes the DNA rearrangement pathway, suggesting that RNA molecules guide genome rearrangement.

Associative learning in biochemical networks

It has been recently suggested that there are likely generic features characterizing the emergence of systems constructed from the self-organization of self-replicating agents acting under one or more selection pressures. Therefore, structures and behaviors at one length scale may be used to infer analogous structures and behaviors at other length scales. Motivated by this suggestion, we seek to characterize various "animate" behaviors in biochemical networks, and the influence that these behaviors have on genomic evolution. Specifically, in this paper, we develop a simple, chemostat-based model illustrating how a process analogous to associative learning can occur in a biochemical network. Associative learning is a form of learning whereby a system "learns" to associate two stimuli with one another. Associative learning, also known as conditioning, is believed to be a powerful learning process at work in the brain (associative learning is essentially "learning by analogy"). In our model, two types of replicating molecules, denoted as A and B, are present in some initial concentration in the chemostat. Molecules A and B are stimulated to replicate by some growth factors, denoted as G(A) and G(B), respectively. It is also assumed that A and B can covalently link, and that the conjugated molecule can be stimulated by either the G(A) or G(B) growth factors (and can be degraded). We show that, if the chemostat is stimulated by both growth factors for a certain time, followed by a time gap during which the chemostat is not stimulated at all, and if the chemostat is then stimulated again by only one of the growth factors, then there will be a transient increase in the number of molecules activated by the other growth factor. Therefore, the chemostat bears the imprint of earlier, simultaneous stimulation with both growth factors, which is indicative of associative learning. It is interesting to note that the dynamics of our model is consistent with certain aspects of Pavlov's original series of conditioning experiments in dogs. We discuss how associative learning can potentially be performed in vitro within RNA, DNA, or peptide networks. We also describe how such a mechanism could be involved in genomic evolution, and suggest relevant bioinformatics studies that could potentially resolve these issues.

RNA-guided DNA assembly

We propose molecular models for homologous DNA recombination events that are guided by either double-stranded RNA (dsRNA) or single-stranded RNA (ssRNA) templates. The models are applied to explain DNA rearrangements in some groups of ciliates, such as Stylonychia or Oxytricha, where extensive gene rearrangement occurs during differentiation of a somatic macronucleus from a germline micronucleus. We describe a model for RNA template guided DNA recombination, such that the template serves as a catalyst that remains unchanged after DNA recombination. This recombination can be seen as topological braiding of the DNA, with the template-guided alignment proceeding through DNA branch migration. We show that a virtual knot diagram can provide a physical representation of the DNA at the time of recombination. Schematically, the braiding process can be represented as a crossing in the virtual knot diagram. The homologous recombination corresponds to removal of the crossings in the knot diagram (called smoothing). We show that if all recombinations are performed at the same time (i.e., simultaneous smoothings of the crossings) then one of the resulting DNA molecules will always contain all of the gene segments in their correct, linear order, which produces the mature DNA sequence.

Intron Evolution and Information processing in the DNA polymerase alpha gene in spirotrichous ciliates: a hypothesis for interconversion between DNA and RNA deletion

BACKGROUND

The somatic DNA molecules of spirotrichous ciliates are present as linear chromosomes containing mostly single-gene coding sequences with short 5' and 3' flanking regions. Only a few conserved motifs have been found in the flanking DNA. Motifs that may play roles in promoting and/or regulating transcription have not been consistently detected. Moreover, comparing subtelomeric regions of 1,356 end-sequenced somatic chromosomes failed to identify more putatively conserved motifs.

RESULTS

We sequenced and compared DNA and RNA versions of the DNA polymerase alpha (pol alpha) gene from nine diverged spirotrichous ciliates. We identified a G-C rich motif aaTACCGC(G/C/T) upstream from transcription start sites in all nine pol alpha orthologs. Furthermore, we consistently found likely polyadenylation signals, similar to the eukaryotic consensus AAUAAA, within 35 nt upstream of the polyadenylation sites. Numbers of introns differed among orthologs, suggesting independent gain or loss of some introns during the evolution of this gene. Finally, we discuss the occurrence of short direct repeats flanking some introns in the DNA pol alpha genes. These introns flanked by direct repeats resemble a class of DNA sequences called internal eliminated sequences (IES) that are deleted from ciliate chromosomes during development.

CONCLUSION

Our results suggest that conserved motifs are present at both 5' and 3' untranscribed regions of the DNA pol alpha genes in nine spirotrichous ciliates. We also show that several independent gains and losses of introns in the DNA pol alpha genes have occurred in the spirotrichous ciliate lineage. Finally, our statistical results suggest that proven introns might also function in an IES removal pathway. This could strengthen a recent hypothesis that introns evolve into IESs, explaining the scarcity of introns in spirotrichs. Alternatively, the analysis suggests that ciliates might occasionally use intron splicing to correct, at the RNA level, failures in IES excision during developmental DNA elimination.

Interconversion of germline-limited and somatic DNA in a scrambled gene

Ciliates have a somatic and a germline nucleus; after sexual conjugation a new somatic nucleus forms from the new zygotic germline nucleus. Formation of the somatic nucleus involves precise elimination of a large portion of DNA sequences from the germline. Here we compare the architecture of the germline and somatic versions of the actin I gene in two geographically isolated strains of Stylonychia lemnae. We show that the structure of the germline gene is surprisingly mercurial, with the distinction between germline-limited and somatic sequences variable over the course of evolution. This is, to our knowledge, the first example of evolutionary swapping of retained versus deleted sequences during ciliate development, with sequences deleted during development that are specifically retained in another strain.

RNA-guided DNA deletion in Tetrahymena: an RNAi-based mechanism for programmed genome rearrangements

Ciliated protozoan are unicellular eukaryotes. Most species in this diverse group display nuclear dualism, a special feature that supports both somatic and germline nuclei in the same cell. Probably due to this unique life style, they exhibit unusual nuclear characteristics that have intrigued researchers for decades. Among them are large-scale DNA rearrangements, which restructure the somatic genome to become drastically different from its germline origin. They resemble the classical phenomenon of chromatin diminution in some nematodes discovered more than a century ago. The mechanisms of such rearrangements, their biological roles, and their evolutionary origins have been difficult to understand. Recent studies have revealed a clear link to RNA interference, and begin to shed light on these issues. Using the simple ciliate Tetrahymena as a model, this chapter summarizes the physical characterization of these processes, describes recent findings that connect them to RNA interference, and discusses the details of their mechanisms, potential roles in genome defense, and possible occurrences in other organisms.

A new scrambled gene in the ciliate Uroleptus

In the germline micronucleus of spirotrichous ciliates, the gene segments, or macronuclear destined sequences (MDSs), that give rise to the somatic macronucleus are interrupted by internal eliminated sequences (IESs). For some genes, the MDSs are not arranged sequentially, but rather are scrambled, in the micronucleus. Three scrambled genes have been extensively studied in many species: actin I, alpha-telomere binding protein, and DNA polymerase alpha. However, in the past decade, no new scrambled genes have been reported, and the prevalence of scrambled genes is still an important question. To screen for scrambled genes, we completely sequenced 11 macronuclear chromosomes in the spirotrich Uroleptus sp., and then pursued their micronuclear organization. This allowed us to identify new scrambled genes, which also display novel features. In this study we describe one of these newly discovered scrambled genes. This gene, tentatively named USG1 (Unknown Scrambled Gene 1), encodes a putative protein of 1016 aa. While the function of this protein product is not clear, dN/dS calculated from the two alleles suggests the encoded protein is under purifying selection. USG1 consists of 16 germline MDSs, of which 14 are located on one locus. The other locus, which is at least 3 kb away from the main locus, contains two scrambled MDSs separated by a nonscrambled IES. Curiously, one MDS and its outgoing (3') pointer (direct repeat) overlap intron splice sites, indicating that these DNA sequences may be under dual (or multiple) constraints. Our findings identify a new scrambled gene in the micronuclear genome of a spirotrichous ciliate, and suggest that even more complicated structures may be present.

The evolutionary origin of a complex scrambled gene

Some species of ciliates undergo massive DNA elimination and genome rearrangement to construct gene-sized "chromosomes" in their somatic nucleus. An example is the extensively scrambled DNA polymerase alpha gene that is broken into 48 pieces and distributed over two unlinked loci in Stylonychia. To understand the emergence of this complex phenomenon during evolution, we examined DNA polymerase alpha genes in several earlier diverging species, representing evolutionary intermediates. Mapping these data onto an evolutionary tree suggests that this gene became extensively fragmented and scrambled over evolutionary time through a series of steps, each leading to greater complexity. Our results also suggest a possible mechanism for intron loss by deletion of intron sequences as DNA during development of the somatic nucleus.

MDS_IES_DB: a database of macronuclear and micronuclear genes in spirotrichous ciliates

Ciliated protozoa have two kinds of nuclei: Macronuclei (MAC) and Micronuclei (MIC). In some ciliate classes, such as spirotrichs, most genes undergo several layers of DNA rearrangement during macronuclear development. Because of such processes, these organisms provide ideal systems for studying mechanisms of recombination and gene rearrangement. Here, we describe a database that contains all spirotrich genes for which both MAC and MIC versions are sequenced, with consistent annotation and easy access to all the features. An interface to query the database is available at http://oxytricha.princeton.edu/dimorphism/database.htm.

Quantum computation, non-demolition measurements, and reflective control in living systems

Internal computation underlies robust non-equilibrium living process. The smallest details of living systems are molecular devices that realize non-demolition quantum measurements. These smaller devices form larger devices (macromolecular complexes), up to living body. The quantum device possesses its own potential internal quantum state (IQS), which is maintained for a prolonged time via reflective error-correction. Decoherence-free IQS can exhibit itself by a creative generation of iteration limits in the real world. It resembles the properties of a quasi-particle, which interacts with the surround, applying decoherence commands to it. In this framework, enzymes are molecular automata of the extremal quantum computer, the set of which maintains highly ordered robust coherent state, and genome represents a concatenation of error-correcting codes into a single reflective set. The biological evolution can be viewed as a functional evolution of measurement constraints in which limits of iteration are established, possessing criteria of perfection and having selective values.

Arithmetic inside the universal genetic code

The first information system emerged on the earth as primordial version of the genetic code and genetic texts. The natural appearance of arithmetic power in such a linguistic milieu is theoretically possible and practical for producing information systems of extremely high efficiency. In this case, the arithmetic symbols should be incorporated into an alphabet, i.e. the genetic code. A number is the fundamental arithmetic symbol produced by the system of numeration. If the system of numeration were detected inside the genetic code, it would be natural to expect that its purpose is arithmetic calculation e.g., for the sake of control, safety, and precise alteration of the genetic texts. The nucleons of amino acids and the bases of nucleic acids seem most suitable for embodiments of digits. These assumptions were used for the analyzing the genetic code. The compressed, life-size, and split representation of the Escherichia coli and Euplotes octocarinatus code versions were considered simultaneously. An exact equilibration of the nucleon sums of the amino acid standard blocks and/or side chains was found repeatedly within specified sets of the genetic code. Moreover, the digital notations of the balanced sums acquired, in decimal representation, the unique form 111, 222...., 999. This form is a consequence of the criterion of divisibility by 037. The criterion could simplify some computing mechanism of a cell if any and facilitate its computational procedure. The cooperative symmetry of the genetic code demonstrates that possibly a zero was invented and used by this mechanism. Such organization of the genetic code could be explained by activities of some hypothetical molecular organelles working as natural biocomputers of digital genetic texts. It is well known that if mutation replaces an amino acid, the change of hydrophobicity is generally weak, while that of size is strong. The antisymmetrical correlation between the amino acid size and the degeneracy number is known as well. It is shown that these and some other familiar properties may be a physicochemical effect of arithmetic inside the genetic code. The "frozen accident" model, giving unlimited freedom to the mapping function, could optimally support the appearance of both arithmetic symbols and physicochemical protection inside the genetic code.

Template-guided recombination for IES elimination and unscrambling of genes in stichotrichous ciliates

The micronuclear versions of genes in stichotrichous ciliates are interrupted by multiple, short, non-coding DNA segments called internal eliminated segments, or IESs. IESs divide a gene into macronuclear destined segments, or MDSs. In some micronuclear genes MDSs are in a scrambled disorder. During development of a micronucleus into a macronucleus after cell mating the IESs are excised from micronuclear genes and the MDSs are spliced in the sequentially correct order. Pairs of short repeat sequences in the ends of MDSs undergo homologous recombination to excise IESs and splice MDSs. However, the repeat sequences are too short to guide unambiguously their own alignment in preparation for recombination. Based on experiments by others on the distantly related ciliate, Paramecium, we propose a molecular model of template-guided recombination to explain the excision of the 100,000-150,000 IESs and splicing of MDSs, including unscrambling, in the genome of stichotrichous ciliates. The model solves the problem of correct pairing of pointers, precisely identifies MDS-IES junctions, and provides for irreversible recombination.

Genome remodeling in ciliated protozoa

The germline genomes of ciliated protozoa are dynamic structures, undergoing massive DNA rearrangement during the formation of a functional macronucleus. Macronuclear development involves chromosome fragmentation coupled with de novo telomere synthesis, numerous DNA splicing events that remove internal segments of DNA, and, in some ciliates, the reordering of scrambled gene segments. Despite the fact that all ciliates share similar forms of DNA rearrangement, there appears to be great diversity in both the nature of the rearranged DNA and the molecular mechanisms involved. Epigenetic effects on rearrangement have also been observed, and recent work suggests that chromatin differentiation plays a role in specifying DNA segments either for rearrangement or for elimination.