SDRAP for annotating scrambled or rearranged genomes.

Genomes sometimes undergo large-scale rearrangements. Programmed genome rearrangements in ciliates offer an extreme example, making them a compelling model system to study DNA rearrangements. Currently, available methods for genome annotation are not adequate for highly scrambled genomes. We present a theoretical framework and software implementation for the systematic extraction and analysis of DNA rearrangement annotations from pairs of genome assemblies corresponding to precursor and product versions. The software makes no assumptions about the structure of the rearrangements, and permits the user to select parameters to suit the data. Compared to previous approaches, this work achieves more complete precursor-product mappings, allows for full transparency and reproducibility, and can be adapted to genomic data from different sources.

Exploration of the Nuclear Proteomes in the Ciliate Oxytricha trifallax.

Nuclear dimorphism is a fundamental feature of ciliated protozoa, which have separate somatic and germline genomes in two distinct organelles within a single cell. The transcriptionally active somatic genome, contained within the physically larger macronucleus, is both structurally and functionally different from the silent germline genome housed in the smaller micronucleus. This difference in genome architecture is particularly exaggerated in Oxytricha trifallax, in which the somatic genome comprises tens of thousands of gene-sized nanochromosomes maintained at a high and variable ploidy, while the germline has a diploid set of megabase-scale chromosomes. To examine the compositional differences between the nuclear structures housing the genomes, we performed a proteomic survey of both types of nuclei and of macronuclear histones using quantitative mass spectrometry. We note distinct differences between the somatic and germline nuclei, with many functional proteins being highly enriched in one of the two nuclei. To validate our conclusions and the efficacy of nuclear separation, we used protein localization through a combination of transformations and immunofluorescence. We also note that the macronuclear histones strikingly display only activating marks, consistent with the conclusion that the macronucleus is the hub of transcription. These observations suggest that the compartmentalization of different genome features into separate structures has been accompanied by a similar specialization of nuclear components that maintain and facilitate the functions of the genomes specific to each nucleus.

Comparative genomics reveals insight into the evolutionary origin of massively scrambled genomes.

Ciliates are microbial eukaryotes that undergo extensive programmed genome rearrangement, a natural genome editing process that converts long germline chromosomes into smaller gene-rich somatic chromosomes. Three well-studied ciliates include Oxytricha trifallax, Tetrahymena thermophila, and Paramecium tetraurelia, but only the Oxytricha lineage has a massively scrambled genome, whose assembly during development requires hundreds of thousands of precisely programmed DNA joining events, representing the most complex genome dynamics of any known organism. Here we study the emergence of such complex genomes by examining the origin and evolution of discontinuous and scrambled genes in the Oxytricha lineage. This study compares six genomes from three species, the germline and somatic genomes for Euplotes woodruffi, Tetmemena sp., and the model ciliate O. trifallax. We sequenced, assembled, and annotated the germline and somatic genomes of E. woodruffi, which provides an outgroup, and the germline genome of Tetmemena sp. We find that the germline genome of Tetmemena is as massively scrambled and interrupted as Oxytricha's: 13.6% of its gene loci require programmed translocations and/or inversions, with some genes requiring hundreds of precise gene editing events during development. This study revealed that the earlier diverged spirotrich, E. woodruffi, also has a scrambled genome, but only roughly half as many loci (7.3%) are scrambled. Furthermore, its scrambled genes are less complex, together supporting the position of Euplotes as a possible evolutionary intermediate in this lineage, in the process of accumulating complex evolutionary genome rearrangements, all of which require extensive repair to assemble functional coding regions. Comparative analysis also reveals that scrambled loci are often associated with local duplications, supporting a gradual model for the origin of complex, scrambled genomes via many small events of DNA duplication and decay.

NT-seq: a chemical-based sequencing method for genomic methylome profiling.

DNA methylation plays vital roles in both prokaryotes and eukaryotes. There are three forms of DNA methylation in prokaryotes: N6-methyladenine (6mA), N4-methylcytosine (4mC), and 5-methylcytosine (5mC). Although many sequencing methods have been developed to sequence specific types of methylation, few technologies can be used for efficiently mapping multiple types of methylation. Here, we present NT-seq for mapping all three types of methylation simultaneously. NT-seq reliably detects all known methylation motifs in two bacterial genomes and can be used for identifying de novo methylation motifs. NT-seq provides a simple and efficient solution for detecting multiple types of DNA methylation.

Transposon debris in ciliate genomes.

The germline genomes of ciliated protists are replete with "junk" DNA insertions that need to be removed for gene expression. Unlike introns, these are spliced as DNA. What is their source, and why are they so abundant? A new study in PLOS Biology supports a classic model of transposon origins.

Transcribed germline-limited coding sequences in Oxytricha 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.

SIGAR: Inferring Features of Genome Architecture and DNA Rearrangements by Split-Read Mapping.

Ciliates are microbial eukaryotes with distinct somatic and germline genomes. Postzygotic development involves extensive remodeling of the germline genome to form somatic chromosomes. Ciliates therefore offer a valuable model for studying the architecture and evolution of programed genome rearrangements. Current studies usually focus on a few model species, where rearrangement features are annotated by aligning reference germline and somatic genomes. Although many high-quality somatic genomes have been assembled, a high-quality germline genome assembly is difficult to obtain due to its smaller DNA content and abundance of repetitive sequences. To overcome these hurdles, we propose a new pipeline, SIGAR (Split-read Inference of Genome Architecture and Rearrangements) to infer germline genome architecture and rearrangement features without a germline genome assembly, requiring only short DNA sequencing reads. As a proof of principle, 93% of rearrangement junctions identified by SIGAR in the ciliate Oxytricha trifallax were validated by the existing germline assembly. We then applied SIGAR to six diverse ciliate species without germline genome assemblies, including Ichthyophthirius multifilii, a fish pathogen. Despite the high level of somatic DNA contamination in each sample, SIGAR successfully inferred rearrangement junctions, short eliminated sequences, and potential scrambled genes in each species. This pipeline enables pilot surveys or exploration of DNA rearrangements in species with limited DNA material access, thereby providing new insights into the evolution of chromosome rearrangements.

Programmed genome rearrangements in Oxytricha produce transcriptionally active extrachromosomal circular DNA.

Extrachromosomal circular DNA (eccDNA) is both a driver of eukaryotic genome instability and a product of programmed genome rearrangements, but its extent had not been surveyed in Oxytricha, a ciliate with elaborate DNA elimination and translocation during development. Here, we captured rearrangement-specific circular DNA molecules across the genome to gain insight into its processes of programmed genome rearrangement. We recovered thousands of circularly excised Tc1/mariner-type transposable elements and high confidence non-repetitive germline-limited loci. We verified their bona fide circular topology using circular DNA deep-sequencing, 2D gel electrophoresis and inverse polymerase chain reaction. In contrast to the precise circular excision of transposable elements, we report widespread heterogeneity in the circular excision of non-repetitive germline-limited loci. We also demonstrate that circular DNAs are transcribed in Oxytricha, producing rearrangement-specific long non-coding RNAs. The programmed formation of thousands of eccDNA molecules makes Oxytricha a model system for studying nucleic acid topology. It also suggests involvement of eccDNA in programmed genome rearrangement.

Programmed Chromosome Deletion in the Ciliate Oxytricha trifallax.

The ciliate Oxytricha trifallax contains two nuclei: a germline micronucleus and a somatic macronucleus. These two nuclei diverge significantly in genomic structure. The micronucleus contains approximately 100 chromosomes of megabase scale, while the macronucleus contains 16,000 gene-sized, high ploidy "nanochromosomes." During its sexual cycle, a copy of the zygotic germline micronucleus develops into a somatic macronucleus via DNA excision and rearrangement. The rearrangement process is guided by multiple RNA-based pathways that program the epigenetic inheritance of sequences in the parental macronucleus of the subsequent generation. Here, we show that the introduction of synthetic DNA molecules homologous to a complete native nanochromosome during the rearrangement process results in either loss or heavy copy number reduction of the targeted nanochromosome in the macronucleus of the subsequent generation. This phenomenon was tested on a variety of nanochromosomes with different micronuclear structures, with deletions resulting in all cases. Deletion of the targeted nanochromosome results in the loss of expression of the targeted genes, including gene knockout phenotypes that were phenocopied using alternative knockdown approaches. Further investigation of the chromosome deletion showed that, although the full length nanochromosome was lost, remnants of the targeted chromosome remain. We were also able to detect the presence of telomeres on these remnants. The chromosome deletions and remnants are epigenetically inherited when backcrossed to wild type strains, suggesting that an undiscovered mechanism programs DNA elimination and cytoplasmically transfers to both daughter cells during conjugation. Programmed deletion of targeted chromosomes provides a novel approach to investigate genome rearrangement and expands the available strategies for gene knockout in Oxytricha trifallax.

Transformation with Artificial Chromosomes in Oxytricha trifallax and Their Applications.

Oxytricha trifallax, like other ciliates, has separate germline and somatic nuclei. The diploid germline genome in the micronucleus is composed of long conventional chromosomes. The macronucleus contains a somatic genome which is naturally fragmented into thousands of kilobase-sized chromosomes. Here, we develop a method to stably incorporate artificial chromosomes into the macronucleus. We report two cases of successful transformation and demonstrate the use of somatic transformation to investigate gene regulation and gene function in Oxytricha We show that the transformed artificial chromosomes are maintained through multiple asexual divisions. Furthermore, they support the transcriptional regulation of the native chromosome from which they were derived and are translated to produce functional proteins. To test if transformed chromosomes are amenable to practical applications, we generated a tagged version of a representative gene (AL1) and used it to co-precipitate associated proteins. This revealed an association with nucleic acid binding proteins, specifically RNA-binding proteins, and RNA immunoprecipitation of AL1 revealed its association with multiple RNAs. The use of artificial chromosomes in Oxytricha enables an array of genetic and molecular biological assays, as well as new avenues of inquiry into the epigenetic programming of macronuclear development and genome rearrangement.

Identification of a DNA N6-Adenine Methyltransferase Complex and Its Impact on Chromatin Organization.

DNA N6-adenine methylation (6mA) has recently been described in diverse eukaryotes, spanning unicellular organisms to metazoa. Here, we report a DNA 6mA methyltransferase complex in ciliates, termed MTA1c. It consists of two MT-A70 proteins and two homeobox-like DNA-binding proteins and specifically methylates dsDNA. Disruption of the catalytic subunit, MTA1, in the ciliate Oxytricha leads to genome-wide loss of 6mA and abolishment of the consensus ApT dimethylated motif. Mutants fail to complete the sexual cycle, which normally coincides with peak MTA1 expression. We investigate the impact of 6mA on nucleosome occupancy in vitro by reconstructing complete, full-length Oxytricha chromosomes harboring 6mA in native or ectopic positions. We show that 6mA directly disfavors nucleosomes in vitro in a local, quantitative manner, independent of DNA sequence. Furthermore, the chromatin remodeler ACF can overcome this effect. Our study identifies a diverged DNA N6-adenine methyltransferase and defines the role of 6mA in chromatin organization.

Capture of complete ciliate chromosomes in single sequencing reads reveals widespread chromosome isoforms.

BACKGROUND: Whole-genome shotgun sequencing, which stitches together millions of short sequencing reads into a single genome, ushered in the era of modern genomics and led to a rapid expansion of the number of genome sequences available. Nevertheless, assembly of short reads remains difficult, resulting in fragmented genome sequences. Ultimately, only a sequencing technology capable of capturing complete chromosomes in a single run could resolve all ambiguities. Even "third generation" sequencing technologies produce reads far shorter than most eukaryotic chromosomes. However, the ciliate Oxytricha trifallax has a somatic genome with thousands of chromosomes averaging only 3.2 kbp, making it an ideal candidate for exploring the benefits of sequencing whole chromosomes without assembly. RESULTS: We used single-molecule real-time sequencing to capture thousands of complete chromosomes in single reads and to update the published Oxytricha trifallax JRB310 genome assembly. In this version, over 50% of the completed chromosomes with two telomeres derive from single reads. The improved assembly includes over 12,000 new chromosome isoforms, and demonstrates that somatic chromosomes derive from variable rearrangements between somatic segments encoded up to 191,000 base pairs away. However, while long reads reduce the need for assembly, a hybrid approach that supplements long-read sequencing with short reads for error correction produced the most complete and accurate assembly, overall. CONCLUSIONS: This assembly provides the first example of complete eukaryotic chromosomes captured by single sequencing reads and demonstrates that traditional approaches to genome assembly can mask considerable structural variation.

Russian Doll Genes and Complex Chromosome Rearrangements in Oxytricha trifallax.

Ciliates have two different types of nuclei per cell, with one acting as a somatic, transcriptionally active nucleus (macronucleus; abbr. MAC) and another serving as a germline nucleus (micronucleus; abbr. MIC). Furthermore, Oxytricha trifallax undergoes extensive genome rearrangements during sexual conjugation and post-zygotic development of daughter cells. These rearrangements are necessary because the precursor MIC loci are often both fragmented and scrambled, with respect to the corresponding MAC loci. Such genome architectures are remarkably tolerant of encrypted MIC loci, because RNA-guided processes during MAC development reorganize the gene fragments in the correct order to resemble the parental MAC sequence. Here, we describe the germline organization of several nested and highly scrambled genes in Oxytricha trifallax These include cases with multiple layers of nesting, plus highly interleaved or tangled precursor loci that appear to deviate from previously described patterns. We present mathematical methods to measure the degree of nesting between precursor MIC loci, and revisit a method for a mathematical description of scrambling. After applying these methods to the chromosome rearrangement maps of O. trifallax we describe cases of nested arrangements with up to five layers of embedded genes, as well as the most scrambled loci in O. trifallax.

Small RNA-mediated regulation of DNA dosage in the ciliate Oxytricha.

Dicer-dependent small noncoding RNAs play important roles in gene regulation in a wide variety of organisms. Endogenous small interfering RNAs (siRNAs) are part of an ancient pathway of transposon control in plants and animals. The ciliate, Oxytricha trifallax, has approximately 16,000 gene-sized chromosomes in its somatic nucleus. Long noncoding RNAs establish high ploidy levels at the onset of sexual development, but the factors that regulate chromosome copy numbers during cell division and growth have been a mystery. We report a novel function of a class of Dicer (Dcl-1)- and RNA-dependent RNA polymerase (RdRP)-dependent endogenous small RNAs in regulating chromosome copy number and gene dosage in O. trifallax Asexually growing populations express an abundant class of 21-nt sRNAs that map to both coding and noncoding regions of most chromosomes. These sRNAs are bound to chromatin and their levels surprisingly do not correlate with mRNA levels. Instead, the levels of these small RNAs correlate with genomic DNA copy number. Reduced sRNA levels in dcl-1 or rdrp mutants lead to concomitant reduction in chromosome copy number. Furthermore, these cells show no signs of transposon activation, but instead display irregular nuclear architecture and signs of replication stress. In conclusion, Oxytricha Dcl-1 and RdRP-dependent small RNAs that derive from the somatic nucleus contribute to the maintenance of gene dosage, possibly via a role in DNA replication, offering a novel role for these small RNAs in eukaryotes.

Non-model model organisms.

Model organisms are widely used in research as accessible and convenient systems to study a particular area or question in biology. Traditionally only a handful of organisms have been widely studied, but modern research tools are enabling researchers to extend the set of model organisms to include less-studied and more unusual systems. This Forum highlights a range of 'non-model model organisms' as emerging systems for tackling questions across the whole spectrum of biology (and beyond), the opportunities and challenges, and the outlook for the future.

Thousands of RNA-cached copies of whole chromosomes are present in the ciliate Oxytricha during development.

The ciliate Oxytricha trifallax maintains two genomes: a germline genome that is active only during sexual conjugation and a transcriptionally active, somatic genome that derives from the germline via extensive sequence reduction and rearrangement. Previously, we found that long noncoding (lnc) RNA "templates"-telomere-containing, RNA-cached copies of mature chromosomes-provide the information to program the rearrangement process. Here we used a modified RNA-seq approach to conduct the first genome-wide search for endogenous, telomere-to-telomere RNA transcripts. We find that during development, Oxytricha produces long noncoding RNA copies for over 10,000 of its 16,000 somatic chromosomes, consistent with a model in which Oxytricha transmits an RNA-cached copy of its somatic genome to the sexual progeny. Both the primary sequence and expression profile of a somatic chromosome influence the temporal distribution and abundance of individual template RNAs. This suggests that Oxytricha may undergo multiple rounds of DNA rearrangement during development. These observations implicate a complex set of thousands of long RNA molecules in the wiring and maintenance of a highly elaborate somatic genome architecture.

Chromosome fusions triggered by noncoding RNA.

Chromosomal fusions are common in normal and cancer cells and can produce aberrant gene products that promote transformation. The mechanisms driving these fusions are poorly understood, but recurrent fusions are widespread. This suggests an underlying mechanism, and some authors have proposed a possible role for RNA in this process. The unicellular eukaryote Oxytricha trifallax displays an exorbitant capacity for natural genome editing, when it rewrites its germline genome to form a somatic epigenome. This developmental process provides a powerful model system to directly test the influence of small noncoding RNAs on chromosome fusion events during somatic differentiation. Here we show that small RNAs are capable of inducing chromosome fusions in 4 distinct cases (out of 4 tested), including one fusion of 3 chromosomes. We further show that these RNA-mediated chromosome fusions are heritable over multiple sexual generations and that transmission of the acquired fusion is associated with endogenous production of novel piRNA molecules that target the fused junction. We also demonstrate the capacity of a long noncoding RNA (lncRNA) to induce chromosome fusion of 2 distal germline loci. These results underscore the ability of short-lived, aberrant RNAs to act as drivers of chromosome fusion events that can be stably transmitted to future generations.

Recurring patterns among scrambled genes in the encrypted genome of the ciliate Oxytricha trifallax.

Some genera of ciliates, such as Oxytricha and Stylonychia, undergo massive genome reorganization during development and provide model organisms to study DNA rearrangement. A common feature of these ciliates is the presence of two types of nuclei: a germline micronucleus and a transcriptionally-active somatic macronucleus containing over 16,000 gene sized "nano-chromosomes". During conjugation the old parental macronucleus disintegrates and a new macronucleus forms from a copy of the zygotic micronucleus. During this process, macronuclear chromosomes assemble through DNA processing events that delete 90-98% of the DNA content of the micronucleus. This includes the deletion of noncoding DNA segments that interrupt precursor DNA regions in the micronucleus, as well as transposons and other germline-limited DNA. Each macronuclear locus may be present in the micronucleus as several nonconsecutive, permuted, and/or inverted DNA segments. Here we investigate the genome-wide range of scrambled gene architectures that describe all precursor-product relationships in Oxytricha trifallax, the first completely sequenced scrambled genome. We find that five general, recurrent patterns in the sets of scrambled micronuclear precursor pieces can describe over 80% of Oxytricha's scrambled genes. These include instances of translocations and inversions, and other specific patterns characterized by alternating stretches of consecutive odd and even DNA segments. Moreover, we find that iterating patterns of alternating odd-even segments up to four times can describe over 96% of the scrambled precursor loci. Recurrence of these highly structured genetic architectures within scrambled genes presumably reflects recurrent evolutionary events that gave rise to over 3000 of scrambled loci in the germline genome.

What Is a Genome?

The genome is often described as the information repository of an organism. Whether millions or billions of letters of DNA, its transmission across generations confers the principal medium for inheritance of organismal traits. Several emerging areas of research demonstrate that this definition is an oversimplification. Here, we explore ways in which a deeper understanding of genomic diversity and cell physiology is challenging the concepts of physical permanence attached to the genome as well as its role as the sole information source for an organism.

Phylogenomic analysis reveals genome-wide purifying selection on TBE transposons in the ciliate Oxytricha.

BACKGROUND: Transposable elements are a major player contributing to genetic variation and shaping genome evolution. Multiple independent transposon domestication events have occurred in ciliates, recruiting transposases to key roles in cellular processes. In the ciliate Oxytricha trifallax, the telomere-bearing elements (TBE), a Tc1/mariner transposon, occupy a significant portion of the germline genome and are involved in programmed genome rearrangements that produce a transcriptionally active somatic nucleus from a copy of the germline nucleus during development. RESULTS: Here we provide a thorough characterization of the distribution and sequences of TBE transposons in the Oxytricha germline genome. We annotate more than 10,000 complete and 24,000 partial TBE sequences. TBEs cluster into four major families and display a preference for either insertion into DNA segments that are retained in the somatic genome or their maintenance at such sites. The three TBE-encoded genes in all four families display dN/dS ratios much lower than 1, suggesting genome-wide purifying selection. We also identify TBE homologs in other ciliate species for phylogenomic analysis. CONCLUSIONS: This paper provides genome-wide characterization of a major class of ciliate transposons. Phylogenomic analysis reveals selective constraints on transposon-encoded genes, shedding light on the evolution and domesticated functions of these transposons.