Tetrahymena ORC contains a ribosomal RNA fragment that participates in rDNA origin recognition

Developmentally Programmed Switches in DNA Replication: Gene Amplification and Genome-Wide Endoreplication in Tetrahymena

Locus-specific gene amplification and genome-wide endoreplication generate the elevated copy number of ribosomal DNA (rDNA, 9000 C) and non-rDNA (90 C) chromosomes in the developing macronucleus of Tetrahymena thermophila. Subsequently, all macronuclear chromosomes replicate once per cell cycle during vegetative growth. Here, we describe an unanticipated, programmed switch in the regulation of replication initiation in the rDNA minichromosome. Early in development, the 21 kb rDNA minichromosome is preferentially amplified from 2 C to ~800 C from well-defined origins, concurrent with genome-wide endoreplication (2 C to 8-16 C) in starved mating Tetrahymena (endoreplication (ER) Phase 1). Upon refeeding, rDNA and non-rDNA chromosomes achieve their final copy number through resumption of just the endoreplication program (ER Phase 2). Unconventional rDNA replication intermediates are generated primarily during ER phase 2, consistent with delocalized replication initiation and possible formation of persistent RNA-DNA hybrids. Origin usage and replication fork elongation are affected in non-rDNA chromosomes as well. Despite the developmentally programmed 10-fold reduction in the ubiquitous eukaryotic initiator, the Origin Recognition Complex (ORC), active initiation sites are more closely spaced in ER phases 1 and 2 compared to vegetative growing cells. We propose that initiation site selection is relaxed in endoreplicating macronuclear chromosomes and may be less dependent on ORC.

Transcriptome analysis of the binucleate ciliate Tetrahymena thermophila with asynchronous nuclear cell cycles

Tetrahymena thermophila harbors two functionally and physically distinct nuclei within a shared cytoplasm. During vegetative growth, the "cell cycles" of the diploid micronucleus and polyploid macronucleus are offset. Micronuclear S phase initiates just before cytokinesis and is completed in daughter cells before onset of macronuclear DNA replication. Mitotic micronuclear division occurs mid-cell cycle, while macronuclear amitosis is coupled to cell division. Here we report the first RNA-seq cell cycle analysis of a binucleated ciliated protozoan. RNA was isolated across 1.5 vegetative cell cycles, starting with a macronuclear G1 population synchronized by centrifugal elutriation. Using MetaCycle, 3244 of the 26,000+ predicted genes were shown to be cell cycle regulated. Proteins present in both nuclei exhibit a single mRNA peak that always precedes their macronuclear function. Nucleus-limited genes, including nucleoporins and importins, are expressed before their respective nucleus-specific role. Cyclin D and A/B gene family members exhibit different expression patterns that suggest nucleus-restricted roles. Periodically expressed genes cluster into seven cyclic patterns. Four clusters have known PANTHER gene ontology terms associated with G1/S and G2/M phase. We propose that these clusters encode known and novel factors that coordinate micro- and macronuclear-specific events such as mitosis, amitosis, DNA replication, and cell division.

An Optimized and Versatile Counter-Flow Centrifugal Elutriation Workflow to Obtain Synchronized Eukaryotic Cells

Cell synchronization is a powerful tool to understand cell cycle events and its regulatory mechanisms. Counter-flow centrifugal elutriation (CCE) is a more generally desirable method to synchronize cells because it does not significantly alter cell behavior and/or cell cycle progression, however, adjusting specific parameters in a cell type/equipment-dependent manner can be challenging. In this paper, we used the unicellular eukaryotic model organism, Tetrahymena thermophila as a testing system for optimizing CCE workflow. Firstly, flow cytometry conditions were identified that reduced nuclei adhesion and improved the assessment of cell cycle stage. We then systematically examined how to achieve the optimal conditions for three critical factors affecting the outcome of CCE, including loading flow rate, collection flow rate and collection volume. Using our optimized workflow, we obtained a large population of highly synchronous G1-phase Tetrahymena as measured by 5-ethynyl-2'-deoxyuridine (EdU) incorporation into nascent DNA strands, bulk DNA content changes by flow cytometry, and cell cycle progression by light microscopy. This detailed protocol can be easily adapted to synchronize other eukaryotic cells.

Time-course profiling of bovine alphaherpesvirus 1.1 transcriptome using multiplatform sequencing

Long-read sequencing (LRS) has become a standard approach for transcriptome analysis in recent years. Bovine alphaherpesvirus 1 (BoHV-1) is an important pathogen of cattle worldwide. This study reports the profiling of the dynamic lytic transcriptome of BoHV-1 using two long-read sequencing (LRS) techniques, the Oxford Nanopore Technologies MinION, and the LoopSeq synthetic LRS methods, using multiple library preparation protocols. In this work, we annotated viral mRNAs and non-coding transcripts, and a large number of transcript isoforms, including transcription start and end sites, as well as splice variants of BoHV-1. Our analysis demonstrated an extremely complex pattern of transcriptional overlaps.

The completed macronuclear genome of a model ciliate Tetrahymena thermophila and its application in genome scrambling and copy number analyses

The ciliate Tetrahymena thermophila has been a powerful model system for molecular and cellular biology. However, some investigations have been limited due to the incomplete closure and sequencing of the macronuclear genome assembly, which for many years has been stalled at 1,158 scaffolds, with large sections of unknown sequences (available in Tetrahymena Genome Database, TGD, http://ciliate.org/ ). Here we completed the first chromosome-level Tetrahymena macronuclear genome assembly, with approximately 300× long Single Molecule, Real-Time reads of the wild-type SB210 cells-the reference strain for the initial macronuclear genome sequencing project. All 181 chromosomes were capped with two telomeres and gaps were entirely closed. The completed genome shows significant improvements over the current assembly (TGD 2014) in both chromosome structure and sequence integrity. The majority of previously identified gene models shown in TGD were retained, with the addition of 36 new genes and 883 genes with modified gene models. The new genome and annotation were incorporated into TGD. This new genome allows for pursuit in some underexplored areas that were far more challenging previously; two of them, genome scrambling and chromosomal copy number, were investigated in this study. We expect that the completed macronuclear genome will facilitate many studies in Tetrahymena biology, as well as multiple lines of research in other eukaryotes.

Sequencing and characterization of the macronuclear rDNA minichromosome of the protozoan Tetrahymena pyriformis

Tetrahymena ribosomal DNA (rDNA) is an ideal system for studying eukaryotic DNA replication and gene transcription. In this study, we developed a new method to isolate rDNA from Tetrahymena cells and used it to sequence and annotate the complete 19,670 bp macronuclear rDNA minichromosome of Tetrahymena pyriformis, a species that lacks the germ-line micronucleus and is unable to undergo sexual reproduction. The key features of T. pyriformis and Tetrahymena thermophila rDNA sequences were then compared. Our results showed (i) the short inverted repeats (M repeats) essential for formation of rDNA minichromosome palindromic structure during sexual reproduction in Tetrahymena are highly conserved in T. pyriformis; (ii) in contrast to T. thermophila, which has two tandem domains that coordinately regulate rDNA replication, T. pyriformis has only a single domain; (iii) the 35S pre-rRNA precursor has 80.25% similarity between the two species; and (iv) the G + C content is higher in the transcribed region than the non-transcribed region in both species, but the GC-skew is more stable in T. pyriformis. The new isolation method and annotated information for the T. pyriformis rDNA minichromosome will provide a useful resource for studying DNA replication and chromosome copy number control in Tetrahymena.

Origins of DNA replication

In all kingdoms of life, DNA is used to encode hereditary information. Propagation of the genetic material between generations requires timely and accurate duplication of DNA by semiconservative replication prior to cell division to ensure each daughter cell receives the full complement of chromosomes. DNA synthesis of daughter strands starts at discrete sites, termed replication origins, and proceeds in a bidirectional manner until all genomic DNA is replicated. Despite the fundamental nature of these events, organisms have evolved surprisingly divergent strategies that control replication onset. Here, we discuss commonalities and differences in replication origin organization and recognition in the three domains of life.

Novel classes of replication-associated transcripts discovered in viruses

The role of RNA molecules in the priming of DNA replication and in providing a template for telomerase extension has been known for decades. Since then, several transcripts have been discovered, which play diverse roles in governing replication, including regulation of RNA primer formation, the recruitment of replication origin (Ori) recognition complex, and the assembly of replication fork. Recent studies on viral transcriptomes have revealed novel classes of replication-associated (ra)RNAs, which are expressed from the genomic locations in close vicinity to the Ori. Many of them overlap the Ori, whereas others are terminated close to the replication origin. These novel transcripts can be both protein-coding and non-coding RNAs. The Ori-overlapping part of the mRNAs is generally either the 5ʹ-untranslated regions (UTRs), or the 3ʹ-UTRs of the longer isoforms. Several raRNAs have been identified in various viral families using primarily third-generation long-read sequencing. Hyper-editing of these transcripts has also been described.

Knock-Down of a Novel snoRNA in Tetrahymena Reveals a Dual Role in 5.8S rRNA Processing and Generation of a 26S rRNA Fragment

In eukaryotes, 18S, 5.8S, and 28S rRNAs are transcribed as precursor molecules that undergo extensive modification and nucleolytic processing to form the mature rRNA species. Central in the process are the small nucleolar RNAs (snoRNAs). The majority of snoRNAs guide site specific chemical modifications but a few are involved in defining pre-rRNA cleavages. Here, we describe an unusual snoRNA (TtnuCD32) belonging to the box C/D subgroup from the ciliate Tetrahymena thermophila. We show that TtnuCD32 is unlikely to function as a modification guide snoRNA and that it is critical for cell viability. Cell lines with genetic knock-down of TtnuCD32 were impaired in growth and displayed two novel and apparently unrelated phenotypes. The most prominent phenotype is the accumulation of processing intermediates of 5.8S rRNA. The second phenotype is the decrease in abundance of a ~100 nt 26S rRNA fragment of unknown function. Sequence analysis demonstrated that TtnuCD32 share features with the essential snoRNA U14 but an alternative candidate (TtnuCD25) was more closely related to other U14 sequences. This, together with the fact that the observed rRNA processing phenotypes were not similar to what has been observed in U14 depleted cells, suggests that TtnuCD32 is a U14 homolog that has gained novel functions.

Nonreplicative functions of the origin recognition complex

Origin recognition complex (ORC), a heteromeric six-subunit complex, is the central component of the eukaryotic pre-replication complex. Recent data from yeast, frogs, flies and mammals present compelling evidence that ORC and its individual subunits have nonreplicative functions as well. The majority of these functions, such as heterochromatin formation, chromosome condensation, and segregation are dependent on ORC-DNA interactions. Furthermore, ORC is involved in the control of cell division via its participation in centrosome duplication and cytokinesis. Recent findings have also demonstrated a direct interaction between ORC and mRNPs and highlighted an essential role of ORC in mRNA nuclear export. Along with the growth of evolutionary complexity of organisms, ORC complex functions become more elaborate and new functions of the ORC sub-complexes and individual subunits have emerged.

Transgenerational function of Tetrahymena Piwi protein Twi8p at distinctive noncoding RNA loci

Transgenerational transmission of genome-regulatory epigenetic information can determine phenotypes in the progeny of sexual reproduction. Sequence specificity of transgenerational regulation derives from small RNAs assembled into Piwi-protein complexes. Known targets of transgenerational regulation are primarily transposons and transposon-derived sequences. Here, we extend the scope of Piwi-mediated transgenerational regulation to include unique noncoding RNA loci. Ciliates such as Tetrahymena have a phenotypically silent germline micronucleus and an expressed somatic macronucleus, which is differentiated anew from a germline genome copy in sexual reproduction. We show that the nuclear-localized Tetrahymena Piwi protein Twi8p shuttles from parental to zygotic macronuclei. Genetic elimination of Twi8p has no phenotype for cells in asexual growth. On the other hand, cells lacking Twi8p arrest in sexual reproduction with zygotic nuclei that retain the germline genome structure, without the DNA elimination and fragmentation required to generate a functional macronucleus. Twi8p-bound small RNAs originate from long-noncoding RNAs with a terminal hairpin, which become detectable in the absence of Twi8p. Curiously, the loci that generate Twi8p-bound small RNAs are essential for asexual cell growth, even though Twi8 RNPs are essential only in sexual reproduction. Our findings suggest the model that Twi8 RNPs act on silent germline chromosomes to permit their conversion to expressed macronuclear chromosomes. Overall this work reveals that a Piwi protein carrying small RNAs from long-noncoding RNA loci has transgenerational function in establishing zygotic nucleus competence for gene expression.

Non-coding Y RNAs associate with early replicating euchromatin in concordance with the origin recognition complex

Non-coding Y RNAs are essential for the initiation of chromosomal DNA replication in vertebrates, yet their association with chromatin during the cell cycle is not characterised. Here, we quantify human Y RNA levels in soluble and chromatin-associated intracellular fractions and investigate, topographically, their dynamic association with chromatin during the cell cycle. We find that, on average, about a million Y RNA molecules are present in the soluble fraction of a proliferating cell, and 5-10-fold less are in association with chromatin. These levels decrease substantially during quiescence. No significant differences are apparent between cancer and non-cancer cell lines. Y RNAs associate with euchromatin throughout the cell cycle. Their levels are 2-4-fold higher in S phase than in G1 phase or mitosis. Y RNAs are not detectable at active DNA replication foci, and re-associate with replicated euchromatin during mid and late S phase. The dynamics and sites of Y1 RNA association with chromatin are in concordance with those of the origin recognition complex (ORC). Our data therefore suggest a functional role of Y RNAs in a common pathway with ORC.

ORC interacts with THSC/TREX-2 and its subunits promote Nxf1 association with mRNP and mRNA export in Drosophila

The origin recognition complex (ORC) of eukaryotes associates with the replication origins and initiates the pre-replication complex assembly. In the literature, there are several reports of interaction of ORC with different RNAs. Here, we demonstrate for the first time a direct interaction of ORC with the THSC/TREX-2 mRNA nuclear export complex. The THSC/TREX-2 was purified from the Drosophila embryonic extract and found to bind with a fraction of the ORC. This interaction occurred via several subunits and was essential for Drosophila viability. Also, ORC was associated with mRNP, which was facilitated by TREX-2. ORC subunits interacted with the Nxf1 receptor mediating the bulk mRNA export. The knockdown of Orc5 led to a drop in the Nxf1 association with mRNP, while Orc3 knockdown increased the level of mRNP-bound Nxf1. The knockdown of Orc5, Orc3 and several other ORC subunits led to an accumulation of mRNA in the nucleus, suggesting that ORC participates in the regulation of the mRNP export.

Diverged composition and regulation of the Trypanosoma brucei origin recognition complex that mediates DNA replication initiation

Initiation of DNA replication depends upon recognition of genomic sites, termed origins, by AAA+ ATPases. In prokaryotes a single factor binds each origin, whereas in eukaryotes this role is played by a six-protein origin recognition complex (ORC). Why eukaryotes evolved a multisubunit initiator, and the roles of each component, remains unclear. In Trypanosoma brucei, an ancient unicellular eukaryote, only one ORC-related initiator, TbORC1/CDC6, has been identified by sequence homology. Here we show that three TbORC1/CDC6-interacting factors also act in T. brucei nuclear DNA replication and demonstrate that TbORC1/CDC6 interacts in a high molecular complex in which a diverged Orc4 homologue and one replicative helicase subunit can also be found. Analysing the subcellular localization of four TbORC1/CDC6-interacting factors during the cell cycle reveals that one factor, TbORC1B, is not a static constituent of ORC but displays S-phase restricted nuclear localization and expression, suggesting it positively regulates replication. This work shows that ORC architecture and regulation are diverged features of DNA replication initiation in T. brucei, providing new insight into this key stage of eukaryotic genome copying.

DNA replication origin activation in space and time

DNA replication begins with the assembly of pre-replication complexes (pre-RCs) at thousands of DNA replication origins during the G1 phase of the cell cycle. At the G1-S-phase transition, pre-RCs are converted into pre-initiation complexes, in which the replicative helicase is activated, leading to DNA unwinding and initiation of DNA synthesis. However, only a subset of origins are activated during any S phase. Recent insights into the mechanisms underlying this choice reveal how flexibility in origin usage and temporal activation are linked to chromosome structure and organization, cell growth and differentiation, and replication stress.

Checkpoint Activation of an Unconventional DNA Replication Program in Tetrahymena

The intra-S phase checkpoint kinase of metazoa and yeast, ATR/MEC1, protects chromosomes from DNA damage and replication stress by phosphorylating subunits of the replicative helicase, MCM2-7. Here we describe an unprecedented ATR-dependent pathway in Tetrahymena thermophila in which the essential pre-replicative complex proteins, Orc1p, Orc2p and Mcm6p are degraded in hydroxyurea-treated S phase cells. Chromosomes undergo global changes during HU-arrest, including phosphorylation of histone H2A.X, deacetylation of histone H3, and an apparent diminution in DNA content that can be blocked by the deacetylase inhibitor sodium butyrate. Most remarkably, the cell cycle rapidly resumes upon hydroxyurea removal, and the entire genome is replicated prior to replenishment of ORC and MCMs. While stalled replication forks are elongated under these conditions, DNA fiber imaging revealed that most replicating molecules are produced by new initiation events. Furthermore, the sole origin in the ribosomal DNA minichromosome is inactive and replication appears to initiate near the rRNA promoter. The collective data raise the possibility that replication initiation occurs by an ORC-independent mechanism during the recovery from HU-induced replication stress.

DNA damage and replication stress activate cell cycle checkpoint responses that protect the integrity of eukaryotic chromosomes. A well-conserved response involves the reversible phosphorylation of the replicative helicase, MCM2-7, which together with the origin recognition complex (ORC) dictates when and where replication initiates in chromosomes. The central role of ORC and MCMs in DNA replication is illustrated by the fact that small changes in abundance of these pre-replicative complex (pre-RC) components are poorly tolerated from yeast to humans. Here we describe an unprecedented replication stress checkpoint response in the early branching eukaryote, Tetrahymena thermophila, that is triggered by the depletion of dNTP pools with hydroxyurea (HU). Instead of transiently phosphorylating MCM subunits, ORC and MCM proteins are physically degraded in HU-treated Tetrahymena. Unexpectedly, upon HU removal the genome is completely and effortlessly replicated prior to replenishment of ORC and MCM components. Using DNA fiber imaging and 2D gel electrophoresis, we show that ORC-dependent mechanisms are bypassed during the recovery phase to produce bidirectional replication forks throughout the genome. Our findings suggest that Tetrahymena enlists an alternative mechanism for replication initiation, and that the underlying process can operate on a genome-wide scale.

Non-coding stem-bulge RNAs are required for cell proliferation and embryonic development in C. elegans

Stem bulge RNAs (sbRNAs) are a family of small non-coding stem-loop RNAs present in Caenorhabditis elegans and other nematodes, the function of which is unknown. Here, we report the first functional characterisation of nematode sbRNAs. We demonstrate that sbRNAs from a range of nematode species are able to reconstitute the initiation of chromosomal DNA replication in the presence of replication proteins in vitro, and that conserved nucleotide sequence motifs are essential for this function. By functionally inactivating sbRNAs with antisense morpholino oligonucleotides, we show that sbRNAs are required for S phase progression, early embryonic development and the viability of C. elegans in vivo. Thus, we demonstrate a new and essential role for sbRNAs during the early development of C. elegans. sbRNAs show limited nucleotide sequence similarity to vertebrate Y RNAs, which are also essential for the initiation of DNA replication. Our results therefore establish that the essential function of small non-coding stem-loop RNAs during DNA replication extends beyond vertebrates.

R-loops and initiation of DNA replication in human cells: a missing link?

The unanticipated widespread occurrence of stable hybrid DNA/RNA structures (R-loops) in human cells and the increasing evidence of their involvement in several human malignancies have invigorated the research on R-loop biology in recent years. Here we propose that physiological R-loop formation at CpG island promoters can contribute to DNA replication origin specification at these regions, the most efficient replication initiation sites in mammalian cells. Quite likely, this occurs by the strand-displacement reaction activating the formation of G-quadruplex structures that target the origin recognition complex (ORC) in the single-stranded conformation. In agreement with this, we found that R-loops co-localize with the ORC within the same CpG island region in a significant fraction of these efficient replication origins, precisely at the position displaying the highest density of G4 motifs. This scenario builds on the connection between transcription and replication in human cells and suggests that R-loop dysregulation at CpG island promoter-origins might contribute to the phenotype of DNA replication abnormalities and loss of genome integrity detected in cancer cells.

Developmental regulation of the Tetrahymena thermophila origin recognition complex

The Tetrahymena thermophila DNA replication machinery faces unique demands due to the compartmentalization of two functionally distinct nuclei within a single cytoplasm, and complex developmental program. Here we present evidence for programmed changes in ORC and MCM abundance that are not consistent with conventional models for DNA replication. As a starting point, we show that ORC dosage is critical during the vegetative cell cycle and development. A moderate reduction in Orc1p induces genome instability in the diploid micronucleus, aberrant division of the polyploid macronucleus, and failure to generate a robust intra-S phase checkpoint response. In contrast to yeast ORC2 mutants, replication initiation is unaffected; instead, replication forks elongation is perturbed, as Mcm6p levels decline in parallel with Orc1p. Experimentally induced down-regulation of ORC and MCMs also impairs endoreplication and gene amplification, consistent with essential roles during development. Unexpectedly Orc1p and Mcm6p levels fluctuate dramatically in developing wild type conjugants, increasing for early cycles of conventional micronuclear DNA replication and macronuclear anlagen replication (endoreplication phase I, rDNA gene amplification). This increase does not reflect the DNA replication load, as much less DNA is synthesized during this developmental window compared to vegetative S phase. Furthermore, although Orc1p levels transiently increase prior to endoreplication phase II, Orc1p and Mcm6p levels decline when the replication load increases and unconventional DNA replication intermediates are produced. We propose that replication initiation is re-programmed to meet different requirements or challenges during the successive stages of Tetrahymena development.

The Origin Recognition Complex is required for site-specific replication initiation in eukaryotic chromosomes. Null mutations are lethal in yeast and metazoa, and hypomorphs induce genome instability, a hallmark of cancer. We exploited the unique biology of Tetrahymena to explore ORC's role in conventional and alternative replication programs. Modest experimental down-regulation of ORC1 induces genome instability in vegetative growing Tetrahymena, and diminishes the capacity to support developmentally regulated endoreplication and gene amplification, consistent with essential roles in all of these processes. ORC mutants fail to activate the ATR checkpoint response, and are compromised in their ability to elongate existing replication forks. Remarkably, ORC and MCM levels fluctuate in unexpected ways during wild type development. Most notably, programmed changes in ORC abundance do not reflect the impending DNA replication load. Relative to the vegetative cell cycle, ORC and MCM levels increase dramatically and are highest early in development, when the replication load is lowest. Conversely, ORC levels are lowest during genome-wide macronuclear endoreplication, when the replication load increases. Endocycling cells generate unconventional replication intermediates that distinguish them from vegetative ORC1 knockdown mutants. The collective data suggest that the dependence on ORC may be relaxed during late stages of macronuclear development.

Identification and molecular structure analysis of a new noncoding RNA, a sbRNA homolog, in the silkworm Bombyx mori genome

The small noncoding group of RNAs called stem-bulge RNAs (sbRNAs), first reported in Caenorhabditis elegans, is described as molecules homologous to the Y RNAs, a specific class of noncoding RNAs that is present in vertebrates. This homology indicates the possibility of the existence of sbRNAs in other invertebrate organisms. In this work, we used bioinformatic tools and conserved sequences of sbRNAs from C. Elegans and Y RNAs to search for homologous sbRNA sequences in the Bombyx mori genome. This analysis led to the discovery of one noncoding gene, which was translated into RNA segments and comparatively analysed with segments from human and hamster Y RNAs and C. elegans sbRNAs in molecular dynamic simulations. This gene represents the first evidence for a new sbRNA-like noncoding RNA, the BmsbRNA gene, in this Lepidoptera genome.

Noncoding RNAs in the regulation of DNA replication

Noncoding RNAs (ncRNAs) have crucial roles in epigenetic, transcriptional, and post-transcriptional regulation. Recent studies have begun to reveal a role of ncRNAs in DNA replication. Here, we review the roles of ncRNAs in regulating different aspects of DNA replication in prokaryotic and eukaryotic systems. We speculate that ncRNAs might function to guide the origin recognition complex (ORC) to chromosomal DNA during replication initiation in higher eukaryotes.

A knowledge-based scoring function for protein-RNA interactions derived from a statistical mechanics-based iterative method

Protein-RNA interactions play important roles in many biological processes. Given the high cost and technique difficulties in experimental methods, computationally predicting the binding complexes from individual protein and RNA structures is pressingly needed, in which a reliable scoring function is one of the critical components. Here, we have developed a knowledge-based scoring function, referred to as ITScore-PR, for protein-RNA binding mode prediction by using a statistical mechanics-based iterative method. The pairwise distance-dependent atomic interaction potentials of ITScore-PR were derived from experimentally determined protein–RNA complex structures. For validation, we have compared ITScore-PR with 10 other scoring methods on four diverse test sets. For bound docking, ITScore-PR achieved a success rate of up to 86% if the top prediction was considered and up to 94% if the top 10 predictions were considered, respectively. For truly unbound docking, the respective success rates of ITScore-PR were up to 24 and 46%. ITScore-PR can be used stand-alone or easily implemented in other docking programs for protein–RNA recognition.

RNA-dependent genome processing during nuclear differentiation: the model systems of stichotrichous ciliates

We introduce ciliated protozoa, and more specifically the stichotrichous ciliates Oxytricha and Stylonychia, as biological model systems for the analysis of programmed DNA-reorganization processes during nuclear differentiation. These include DNA excision, DNA elimination, reordering of gene segments and specific gene amplification. We show that small nuclear RNAs specify DNA sequences to be excised or retained, but also discuss the need for a RNA template molecule derived from the parental nucleus for these processes. This RNA template guides reordering of gene segments to become functional genes and determines gene copy number in the differentiated nucleus. Since the template is derived from the parental macronucleus, gene reordering and DNA amplification are inherited in a non-Mendelian epigenetic manner.

Human origin recognition complex binds preferentially to G-quadruplex-preferable RNA and single-stranded DNA

Origin recognition complex (ORC), consisting of six subunits ORC1-6, is known to bind to replication origins and function in the initiation of DNA replication in eukaryotic cells. In contrast to the fact that Saccharomyces cerevisiae ORC recognizes the replication origin in a sequence-specific manner, metazoan ORC has not exhibited strict sequence-specificity for DNA binding. Here we report that human ORC binds preferentially to G-quadruplex (G4)-preferable G-rich RNA or single-stranded DNA (ssDNA). We mapped the G-rich RNA-binding domain in the ORC1 subunit, in a region adjacent to its ATPase domain. This domain itself has an ability to preferentially recognize G4-preferable sequences of ssDNA. Furthermore, we found, by structure modeling, that the G-rich RNA-binding domain is similar to the N-terminal portion of AdoMet_MTase domain of mammalian DNA methyltransferase 1. Therefore, in contrast with the binding to double-stranded DNA, human ORC has an apparent sequence preference with respect to its RNA/ssDNA binding. Interestingly, this specificity coincides with the common signature present in most of the human replication origins. We expect that our findings provide new insights into the regulations of function and chromatin binding of metazoan ORCs.

Genome-wide mapping of human DNA-replication origins: levels of transcription at ORC1 sites regulate origin selection and replication timing

We report the genome-wide mapping of ORC1 binding sites in mammals, by chromatin immunoprecipitation and parallel sequencing (ChIP-seq). ORC1 binding sites in HeLa cells were validated as active DNA replication origins (ORIs) using Repli-seq, a method that allows identification of ORI-containing regions by parallel sequencing of temporally ordered replicating DNA. ORC1 sites were universally associated with transcription start sites (TSSs) of coding or noncoding RNAs (ncRNAs). Transcription levels at the ORC1 sites directly correlated with replication timing, suggesting the existence of two classes of ORIs: those associated with moderate/high transcription levels (≥1 RNA copy/cell), firing in early S and mapping to the TSSs of coding RNAs; and those associated with low transcription levels (<1 RNA copy/cell), firing throughout the entire S and mapping to TSSs of ncRNAs. These findings are compatible with a scenario whereby TSS expression levels influence the efficiency of ORC1 recruitment at G₁ and the probability of firing during S.

Emerging players in the initiation of eukaryotic DNA replication

Faithful duplication of the genome in eukaryotes requires ordered assembly of a multi-protein complex called the pre-replicative complex (pre-RC) prior to S phase; transition to the pre-initiation complex (pre-IC) at the beginning of DNA replication; coordinated progression of the replisome during S phase; and well-controlled regulation of replication licensing to prevent re-replication. These events are achieved by the formation of distinct protein complexes that form in a cell cycle-dependent manner. Several components of the pre-RC and pre-IC are highly conserved across all examined eukaryotic species. Many of these proteins, in addition to their bona fide roles in DNA replication are also required for other cell cycle events including heterochromatin organization, chromosome segregation and centrosome biology. As the complexity of the genome increases dramatically from yeast to human, additional proteins have been identified in higher eukaryotes that dictate replication initiation, progression and licensing. In this review, we discuss the newly discovered components and their roles in cell cycle progression.

A nonredundant structure dataset for benchmarking protein-RNA computational docking

Protein-RNA interactions play an important role in many biological processes. The ability to predict the molecular structures of protein-RNA complexes from docking would be valuable for understanding the underlying chemical mechanisms. We have developed a novel nonredundant benchmark dataset for protein-RNA docking and scoring. The diverse dataset of 72 targets consists of 52 unbound-unbound test complexes, and 20 unbound-bound test complexes. Here, unbound-unbound complexes refer to cases in which both binding partners of the cocrystallized complex are either in apo form or in a conformation taken from a different protein-RNA complex, whereas unbound-bound complexes are cases in which only one of the two binding partners has another experimentally determined conformation. The dataset is classified into three categories according to the interface root mean square deviation and the percentage of native contacts in the unbound structures: 49 easy, 16 medium, and 7 difficult targets. The bound and unbound cases of the benchmark dataset are expected to benefit the development and improvement of docking and scoring algorithms for the docking community. All the easy-to-view structures are freely available to the public at http://zoulab.dalton.missouri.edu/RNAbenchmark/.

Identification of ORC1/CDC6-interacting factors in Trypanosoma brucei reveals critical features of origin recognition complex architecture

DNA Replication initiates by formation of a pre-replication complex on sequences termed origins. In eukaryotes, the pre-replication complex is composed of the Origin Recognition Complex (ORC), Cdc6 and the MCM replicative helicase in conjunction with Cdt1. Eukaryotic ORC is considered to be composed of six subunits, named Orc1–6, and monomeric Cdc6 is closely related in sequence to Orc1. However, ORC has been little explored in protists, and only a single ORC protein, related to both Orc1 and Cdc6, has been shown to act in DNA replication in Trypanosoma brucei. Here we identify three highly diverged putative T. brucei ORC components that interact with ORC1/CDC6 and contribute to cell division. Two of these factors are so diverged that we cannot determine if they are eukaryotic ORC subunit orthologues, or are parasite-specific replication factors. The other we show to be a highly diverged Orc4 orthologue, demonstrating that this is one of the most widely conserved ORC subunits in protists and revealing it to be a key element of eukaryotic ORC architecture. Additionally, we have examined interactions amongst the T. brucei MCM subunits and show that this has the conventional eukaryotic heterohexameric structure, suggesting that divergence in the T. brucei replication machinery is limited to the earliest steps in origin licensing.

Whole genome studies of Tetrahymena

Within the past decade, genomic studies have emerged as essential and highly productive tools to explore the biology of Tetrahymena thermophila. The current major resources, which have been extensively mined by the research community, are the annotated macronuclear genome assembly, transcriptomic data and the databases that house this information. Efforts in progress will soon improve these data sources and expand their scope, including providing annotated micronuclear and comparative genomic sequences. Future studies of Tetrahymena cell and molecular biology, development, physiology, evolution and ecology will benefit greatly from these resources and the advanced genomic technologies they enable.

Different roles of the human Orc6 protein in the replication initiation process

In eukaryotes, binding of the six-subunit origin recognition complex (ORC) to DNA provides an interactive platform for the sequential assembly of pre-replicative complexes. This process licenses replication origins competent for the subsequent initiation step. Here, we analyze the contribution of human Orc6, the smallest subunit of ORC, to DNA binding and pre-replicative complex formation. We show that Orc6 not only interacts with Orc1-Orc5 but also with the initiation factor Cdc6. Biochemical and imaging experiments reveal that this interaction is required for licensing DNA replication competent. Furthermore, we demonstrate that Orc6 contributes to the interaction of ORC with the chaperone protein HMGA1a (high mobility group protein A1a). Binding of human ORC to replication origins is not specified at the level of DNA sequence and the functional organization of origins is poorly understood. We have identified HMGA1a as one factor that might direct ORC to AT-rich heterochromatic regions. The systematic analysis of the interaction between ORC and HMGA1a revealed that Orc6 interacts with the acidic C-terminus of HMGA1a and also with its AT-hooks. Both domains support autonomous replication if targeted to DNA templates. As such, Orc6 functions at different stages of the replication initiation process. Orc6 can interact with ORC chaperone proteins such as HMGA1a to facilitate chromatin binding of ORC and is also an essential factor for pre-RC formation.

Characterization of Leishmania donovani MCM4: expression patterns and interaction with PCNA

Events leading to origin firing and fork elongation in eukaryotes involve several proteins which are mostly conserved across the various eukaryotic species. Nuclear DNA replication in trypanosomatids has thus far remained a largely uninvestigated area. While several eukaryotic replication protein orthologs have been annotated, many are missing, suggesting that novel replication mechanisms may apply in this group of organisms. Here, we characterize the expression of Leishmania donovani MCM4, and find that while it broadly resembles other eukaryotes, noteworthy differences exist. MCM4 is constitutively nuclear, signifying that, unlike what is seen in S.cerevisiae, varying subcellular localization of MCM4 is not a mode of replication regulation in Leishmania. Overexpression of MCM4 in Leishmania promastigotes causes progress through S phase faster than usual, implicating a role for MCM4 in the modulation of cell cycle progression. We find for the first time in eukaryotes, an interaction between any of the proteins of the MCM2-7 (MCM4) and PCNA. MCM4 colocalizes with PCNA in S phase cells, in keeping with the MCM2-7 complex being involved not only in replication initiation, but fork elongation as well. Analysis of a LdMCM4 mutant indicates that MCM4 interacts with PCNA via the PIP box motif of MCM4 - perhaps as an integral component of the MCM2-7 complex, although we have no direct evidence that MCM4 harboring a PIP box mutation can still functionally associate with the other members of the MCM2-7 complex- and the PIP box motif is important for cell survival and viability. In Leishmania, MCM4 may possibly help in recruiting PCNA to chromatin, a role assigned to MCM10 in other eukaryotes.

The midblastula transition defines the onset of Y RNA-dependent DNA replication in Xenopus laevis

Noncoding Y RNAs are essential for the initiation of chromosomal DNA replication in mammalian cell extracts, but their role in this process during early vertebrate development is unknown. Here, we use antisense morpholino nucleotides (MOs) to investigate Y RNA function in Xenopus laevis and zebrafish embryos. We show that embryos in which Y RNA function is inhibited by MOs develop normally until the midblastula transition (MBT) but then fail to replicate their DNA and die before gastrulation. Consistent with this observation, Y RNA function is not required for DNA replication in Xenopus egg extracts but is required for replication in a post-MBT cell line. Y RNAs do not bind chromatin in karyomeres before MBT, but they associate with interphase nuclei after MBT in an origin recognition complex (ORC)-dependent manner. Y RNA-specific MOs inhibit the association of Y RNAs with ORC, Cdt1, and HMGA1a proteins, suggesting that these molecular associations are essential for Y RNA function in DNA replication. The MBT is thus a transition point between Y RNA-independent and Y RNA-dependent control of vertebrate DNA replication. Our data suggest that in vertebrates Y RNAs function as a developmentally regulated layer of control over the evolutionarily conserved eukaryotic DNA replication machinery.

Dynamic interaction of Y RNAs with chromatin and initiation proteins during human DNA replication

Non-coding Y RNAs are required for the initiation of chromosomal DNA replication in mammalian cells. It is unknown how they perform this function or if they associate with a nuclear structure during DNA replication. Here, we investigate the association of Y RNAs with chromatin and their interaction with replication proteins during DNA replication in a human cell-free system. Our results show that fluorescently labelled Y RNAs associate with unreplicated euchromatin in late G1 phase cell nuclei before the initiation of DNA replication. Following initiation, Y RNAs are displaced locally from nascent and replicated DNA present in replication foci. In intact human cells, a substantial fraction of endogenous Y RNAs are associated with G1 phase nuclei, but not with G2 phase nuclei. Y RNAs interact and colocalise with the origin recognition complex (ORC), the pre-replication complex (pre-RC) protein Cdt1, and other proteins implicated in the initiation of DNA replication. These data support a molecular 'catch and release' mechanism for Y RNA function during the initiation of chromosomal DNA replication, which is consistent with Y RNAs acting as replication licensing factors.

PRIDB: a Protein-RNA interface database

The Protein–RNA Interface Database (PRIDB) is a comprehensive database of protein–RNA interfaces extracted from complexes in the Protein Data Bank (PDB). It is designed to facilitate detailed analyses of individual protein–RNA complexes and their interfaces, in addition to automated generation of user-defined data sets of protein–RNA interfaces for statistical analyses and machine learning applications. For any chosen PDB complex or list of complexes, PRIDB rapidly displays interfacial amino acids and ribonucleotides within the primary sequences of the interacting protein and RNA chains. PRIDB also identifies ProSite motifs in protein chains and FR3D motifs in RNA chains and provides links to these external databases, as well as to structure files in the PDB. An integrated JMol applet is provided for visualization of interacting atoms and residues in the context of the 3D complex structures. The current version of PRIDB contains structural information regarding 926 protein–RNA complexes available in the PDB (as of 10 October 2010). Atomic- and residue-level contact information for the entire data set can be downloaded in a simple machine-readable format. Also, several non-redundant benchmark data sets of protein–RNA complexes are provided. The PRIDB database is freely available online at http://bindr.gdcb.iastate.edu/PRIDB.

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.

Non-coding RNAs: new players in the field of eukaryotic DNA replication

The machinery required for the replication of eukaryotic chromosomal DNA is made up of proteins whose function, structure and main interaction partners are evolutionarily conserved. Several new cases have been reported recently, however, in which non-coding RNAs play additional and specialised roles in the initiation of eukaryotic DNA replication in different classes of organisms. These non-coding RNAs include Y RNAs in vertebrate somatic cells, 26T RNA in somatic macronuclei of the ciliate Tetrahymena, and G-rich RNA in the Epstein-Barr DNA tumour virus and its human host cells. Here, I will give an overview of the experimental evidence in favour of roles for these non-coding RNAs in the regulation of eukaryotic DNA replication, and compare and contrast their biosynthesis and mechanisms of action.

Y RNA functions at the initiation step of mammalian chromosomal DNA replication

Non-coding Y RNAs have recently been identified as essential novel factors for chromosomal DNA replication in mammalian cell nuclei, but mechanistic details of their function have not been defined. Here, we identify the execution point for Y RNA function during chromosomal DNA replication in a mammalian cell-free system. We determined the effect of degradation of Y3 RNA on replication origin activation and on fork progression rates at single-molecule resolution by DNA combing and nascent-strand analysis. Degradation of Y3 RNA inhibits the establishment of new DNA replication forks at the G1- to S-phase transition and during S phase. This inhibition is negated by addition of exogenous Y1 RNA. By contrast, progression rates of DNA replication forks are not affected by degradation of Y3 RNA or supplementation with exogenous Y1 RNA. These data indicate that Y RNAs are required for the establishment, but not for the elongation, of chromosomal DNA replication forks in mammalian cell nuclei. We conclude that the execution point for non-coding Y RNA function is the activation of chromosomal DNA replication origins.

Eukaryotic DNA replication control: lock and load, then fire

The initiation of chromosomal DNA replication involves initiator proteins that recruit and load hexameric DNA helicases at replication origins. This helicase loading step is tightly regulated in bacteria and eukaryotes. In contrast to the situation in bacteria, the eukaryotic helicase is loaded in an inactive form. This extra 'lock and load' mechanism in eukaryotes allows regulation of a second step, helicase activation. The temporal separation of helicase loading and activation is crucial for the coordination of DNA replication with cell growth and extracellular signals, the prevention of re-replication and the control of origin activity in response to replication stress. Initiator proteins in bacteria and eukaryotes are structurally homologous; yet the replicative helicases they load are unrelated. Understanding how these helicases are loaded and how they act during unwinding may have important implications for understanding how DNA replication is regulated in different domains of life.

The origin recognition complex protein family

The proteins of the origin recognition complex are found throughout all eukaryotes and have roles beyond that of DNA replication.

Origin recognition complex (ORC) proteins were first discovered as a six-subunit assemblage in budding yeast that promotes the initiation of DNA replication. Orc1-5 appear to be present in all eukaryotes, and include both AAA+ and winged-helix motifs. A sixth protein, Orc6, shows no structural similarity to the other ORC proteins, and is poorly conserved between budding yeast and most other eukaryotic species. The replication factor Cdc6 has extensive sequence similarity with Orc1 and phylogenetic analysis suggests the genes that encode them may be paralogs. ORC proteins have also been found in the archaea, and the bacterial DnaA replication protein has ORC-like functional domains. In budding yeast, Orc1-6 are bound to origins of DNA replication throughout the cell cycle. Following association with Cdc6 in G1 phase, the sequential hydrolysis of Cdc6 - then ORC-bound ATP loads the Mcm2-7 helicase complex onto DNA. Localization of ORC subunits to the kinetochore and centrosome during mitosis and to the cleavage furrow during cytokinesis has been observed in metazoan cells and, along with phenotypes observed following knockdown with short interfering RNAs, point to additional roles at these cell-cycle stages. In addition, ORC proteins function in epigenetic gene silencing through interactions with heterochromatin factors such as Sir1 in budding yeast and HP1 in higher eukaryotes. Current avenues of research have identified roles for ORC proteins in the development of neuronal and muscle tissue, and are probing their relationship to genome integrity.

Differential targeting of Tetrahymena ORC to ribosomal DNA and non-rDNA replication origins

The Tetrahymena thermophila origin recognition complex (ORC) contains an integral RNA subunit, 26T RNA, which confers specificity to the amplified ribosomal DNA (rDNA) origin by base pairing with an essential cis-acting replication determinant--the type I element. Using a plasmid maintenance assay, we identified a 6.7 kb non-rDNA fragment containing two closely associated replicators, ARS1-A (0.8 kb) and ARS1-B (1.2 kb). Both replicators lack type I elements and hence complementarity to 26T RNA, suggesting that ORC is recruited to these sites by an RNA-independent mechanism. Consistent with this prediction, although ORC associated exclusively with origin sequences in the 21 kb rDNA minichromosome, the interaction between ORC and the non-rDNA ARS1 chromosome changed across the cell cycle. In G(2) phase, ORC bound to all tested sequences in a 60 kb interval spanning ARS1-A/B. Remarkably, ORC and Mcm6 associated with just the ARS1-A replicator in G(1) phase when pre-replicative complexes assemble. We propose that ORC is stochastically deposited onto newly replicated non-rDNA chromosomes and subsequently targeted to preferred initiation sites prior to the next S phase.

RNA-dependent recruitment of the origin recognition complex

The origin recognition complex (ORC) has an important function in determining the initiation sites of DNA replication. In higher eukaryotes, ORC lacks sequence-specific DNA binding, and the mechanisms of ORC recruitment and origin determination are poorly understood. ORC is recruited with high efficiency to the Epstein-Barr virus origin of plasmid replication (OriP) through a complex mechanism involving interactions with the virus-encoded EBNA1 protein. We present evidence that ORC recruitment to OriP and DNA replication function depends on RGG-like motifs, referred to as LR1 and LR2, in the EBNA1 amino-terminal domain. Moreover, we show that LR1 and LR2 recruitment of ORC is RNA dependent. HMGA1a, which can functionally substitute for LR1 and LR2 domain, can also recruit ORC in an RNA-dependent manner. EBNA1 and HMGA1a RGG motifs bound to structured G-rich RNA, as did ORC1 peptides, which interact with EBNA1. RNase A treatment of cellular chromatin released a fraction of the total ORC, suggesting that ORC association with chromatin, and possibly cellular origins, is stabilized by RNA. We propose that structural RNA molecules mediate ORC recruitment at some cellular and viral origins, similar to OriP.

Expression and subcellular localization of ORC1 in Leishmania major

The mechanism of DNA replication is highly conserved in eukaryotes, with the process being preceded by the ordered assembly of pre-replication complexes (pre-RCs). Pre-RC formation is triggered by the association of the origin replication complex (ORC) with chromatin. Leishmania major appears to have only one ORC ortholog, ORC1. ORC1 in other eukaryotes is the largest of the ORC subunits and is believed to play a significant role in modulating replication initiation. Here we report for the first time, the cloning of ORC1 from L. major, and the analysis of its expression in L. major promastigotes. In human cells ORC1 levels have been found to be upregulated in G1 and subsequently degraded, thus playing a role in controlling replication initiation. We examine the subcellular localization of L. major ORC1 in relation to the different stages of the cell cycle. Our results show that, unlike what is widely believed to be the case with ORC1 in human cells, ORC1 in L. major is nuclear at all stages of the cell cycle.