RNA-mediated epigenetic programming of a genome-rearrangement pathway

Genes and proteins expressed at different life cycle stages in the model protist Euplotes vannus revealed by both transcriptomic and proteomic approaches

Sexual reproduction first appeared in unicellular protists and has continued to be an essential biological process in almost all eukaryotes. Ciliated protists, which contain both germline and somatic genomes within a single cell, have evolved a special form of sexual reproduction called conjugation that involves mitosis, meiosis, fertilization, nuclear differentiation, genome rearrangement, and the development of unique cellular structures. The molecular basis and mechanisms of conjugation vary dramatically among ciliates, and many details of the process and its regulation are still largely unknown. In order to better comprehend these processes and mechanisms from an evolutionary perspective, this study provides the first comprehensive overview of the transcriptome and proteome profiles during the entire life cycle of the newly-established marine model ciliate Euplotes vannus. Transcriptome analyses from 14 life cycle stages (three vegetative stages and 11 sexual stages) revealed over 26,000 genes that are specifically expressed at different stages, many of which are related to DNA replication, transcription, translation, mitosis, meiosis, nuclear differentiation, and/or genome rearrangement. Quantitative proteomic analyses identified 338 proteins with homologs associated with conjugation and/or somatic nuclear development in other ciliates, including dicer-like proteins, Hsp90 proteins, RNA polymerase II and transcription elongation factors, ribosomal-associated proteins, and ubiquitin-related proteins. Four of these homologs belong to the PIWI family, each with different expression patterns identified and confirmed by RT-qPCR, which may function in small RNA-mediated genome rearrangement. Proteins involved in the nonhomologous end-joining pathway are induced early during meiosis and accumulate in the developing new somatic nucleus, where more than 80% of the germline sequences are eliminated from the somatic genome. A number of new candidate genes and proteins likely to play roles in conjugation and its related genome rearrangements have also been revealed. The gene expression profiles reported here will be valuable resources for further studies of the origin and evolution of sexual reproduction in this new model species.

Conserved features and diversity attributes of chimeric RNAs across accessions in four plants

As a non-collinear expression form of genetic information, chimeric RNAs increase the complexity of transcriptome in diverse organisms. Although chimeric RNAs have been identified in plants, few common features have been revealed. Here, we systemically explored the landscape of chimeric RNAs across multi-accession and multi-tissue using pan-genome and transcriptome data of four plants: rice, maize, soybean, and Arabidopsis. Among the four species, conserved characteristics of breakpoints and parental genes were discovered. In each species, chimeric RNAs displayed a high level of diversity among accessions, and the clustering of accessions using chimeric events was generally concordant with clustering based on genomic variants, implying a general relationship between genetic variations and chimeric RNAs. Through mass spectrometry, we confirmed a fusion protein OsNDC1-OsGID1L2 and observed its subcellular localization, which differed from the original proteins. Phenotypic cues in transgenic rice suggest the potential functions of OsNDC1-OsGID1L2. Moreover, an intriguing chimeric event Os01g0216500-Os01g0216900, generated by a large deletion in basmati rice, also exists in another accession without the deletion, demonstrating its convergence in evolution. Our results illuminate the characteristics and hint at the evolutionary implications of plant chimeric RNAs, which serve as a supplement to genetic variations, thus expanding our understanding of genetic diversity.

Dynamic DNA N 6-adenine methylation (6mA) governs the encystment process, showcased in the unicellular eukaryote Pseudocohnilembus persalinus

The formation of resting cysts commonly found in unicellular eukaryotes is a complex and highly regulated survival strategy against environmental stress that involves drastic physiological and biochemical changes. Although most studies have focused on the morphology and structure of cysts, little is known about the molecular mechanisms that control this process. Recent studies indicate that DNA N 6-adenine methylation (6mA) could be dynamically changing in response to external stimuli; however, its potential role in the regulation of cyst formation remains unknown. We used the ciliate Pseudocohnilembus persalinus, which can be easily induced to form cysts to investigate the dynamic pattern of 6mA in trophonts and cysts. Single-molecule real-time (SMRT) sequencing reveals high levels of 6mA in trophonts that decrease in cysts, along with a conversion of symmetric 6mA to asymmetric 6mA. Further analysis shows that 6mA, a mark of active transcription, is involved in altering the expression of encystment-related genes through changes in 6mA levels and 6mA symmetric-to-asymmetric conversion. Most importantly, we show that reducing 6mA levels by knocking down the DNA 6mA methyltransferase PpAMT1 accelerates cyst formation. Taken together, we characterize the genome-wide 6mA landscape in P. persalinus and provide insights into the role of 6mA in gene regulation under environmental stress in eukaryotes. We propose that 6mA acts as a mark of active transcription to regulate the encystment process along with symmetric-to-asymmetric conversion, providing important information for understanding the molecular response to environmental cues from the perspective of 6mA modification.

Application of RNA interference and protein localization to investigate housekeeping and developmentally regulated genes in the emerging model protozoan Paramecium caudatum

Unicellular eukaryotes represent tremendous evolutionary diversity. However, the molecular mechanisms underlying this diversity remain largely unexplored, partly due to a limitation of genetic tools to only a few model species. Paramecium caudatum is a well-known unicellular eukaryote with an unexpectedly large germline genome, of which only two percent is retained in the somatic genome following sexual processes, revealing extensive DNA elimination. However, further progress in understanding the molecular mechanisms governing this process is hampered by a lack of suitable genetic tools. Here, we report the successful application of gene knockdown and protein localization methods to interrogate the function of both housekeeping and developmentally regulated genes in P. caudatum. Using these methods, we achieved the expected phenotypes upon RNAi by feeding, and determined the localization of these proteins by microinjection of fusion constructs containing fluorescent protein or antibody tags. Lastly, we used these methods to reveal that P. caudatum PiggyMac, a domesticated piggyBac transposase, is essential for sexual development, and is likely to be an active transposase directly involved in DNA cleavage. The application of these methods lays the groundwork for future studies of gene function in P. caudatum and can be used to answer important biological questions in the future.

Programmed chromosome fragmentation in ciliated protozoa: multiple means to chromosome ends

SUMMARYCiliated protozoa undergo large-scale developmental rearrangement of their somatic genomes when forming a new transcriptionally active macronucleus during conjugation. This process includes the fragmentation of chromosomes derived from the germline, coupled with the efficient healing of the broken ends by de novo telomere addition. Here, we review what is known of developmental chromosome fragmentation in ciliates that have been well-studied at the molecular level (Tetrahymena, Paramecium, Euplotes, Stylonychia, and Oxytricha). These organisms differ substantially in the fidelity and precision of their fragmentation systems, as well as in the presence or absence of well-defined sequence elements that direct excision, suggesting that chromosome fragmentation systems have evolved multiple times and/or have been significantly altered during ciliate evolution. We propose a two-stage model for the evolution of the current ciliate systems, with both stages involving repetitive or transposable elements in the genome. The ancestral form of chromosome fragmentation is proposed to have been derived from the ciliate small RNA/chromatin modification process that removes transposons and other repetitive elements from the macronuclear genome during development. The evolution of this ancestral system is suggested to have potentiated its replacement in some ciliate lineages by subsequent fragmentation systems derived from mobile genetic elements.

C-value paradox: Genesis in misconception that natural selection follows anthropocentric parameters of 'economy' and 'optimum'

C-value paradox refers to the lack of correlation between biological complexity and the intuitively expected protein-coding genomic information or DNA content. Here I discuss five questions about this paradox: i) Do biologically complex organisms carry more protein-coding genes? ii) Does variable accumulation of selfish/ junk/ parasitic DNA underlie the c-value paradox? iii) Can nucleoskeletal or nucleotypic function of DNA explain the enigma of orders of magnitude high levels of DNA in some 'lower' taxa or in taxonomically related species? iv) Can the newly understood noncoding but functional DNA explain the c-value paradox? and, v) Does natural selection uniformly apply the anthropocentric parameters for 'optimum' and 'economy'? Answers to Q.1-5 are largely negative. Biology presents numerous 'anomalous' examples where the same end function/ phenotype is attained in different organisms through astoundingly diverse ways that appear 'illogical' in our perceptions. Such evolutionary oddities exist because natural selection, unlike a designer, exploits random and stochastic events to modulate the existing system. Consequently, persistence of the new-found 'solution/s' often appear bizarre, uneconomic, and therefore, paradoxical to human logic. The unexpectedly high c-values in diverse organisms are irreversible evolutionary accidents that persisted, and the additional DNA often got repurposed over the evolutionary time scale. Therefore, the c-value paradox is a redundant issue. Future integrative biological studies should address evolutionary mechanisms and processes underlying sporadic DNA expansions/ contractions, and how the newly acquired DNA content has been repurposed in diverse groups.

Adversary of DNA integrity: A long non-coding RNA stimulates driver oncogenic chromosomal rearrangement in human thyroid cells

The flurry of publications devoted to the functions of long non-coding RNAs (lncRNAs) published in the last decade leaves no doubt about the exceptional importance of lncRNAs in various areas including tumor biology. However, contribution of lncRNAs to the early stages of oncogenesis remains poorly understood. In this study we explored a new role for lncRNAs: stimulation of specific chromosomal rearrangements upon DNA damage. We demonstrated that lncRNA CASTL1 (ENSG00000269945) stimulates the formation of the CCDC6-RET inversion (RET/PTC1) in human thyroid cells subjected to radiation or chemical DNA damage. Facilitation of chromosomal rearrangement requires lncRNA to contain regions complementary to the introns of both CCDC6 and RET genes as deletion of these regions deprives CASTL1 of the ability to stimulate the gene fusion. We found that CASTL1 expression is elevated in tumors with CCDC6-RET fusion which is the most frequent rearrangement in papillary thyroid carcinoma. Our results open a new venue for the studies of early oncogenesis in various tumor types, especially those associated with physical or chemical DNA damage.

Reversing the Central Dogma: RNA-guided control of DNA in epigenetics and genome editing

The Central Dogma of the flow of genetic information is arguably the crowning achievement of 20th century molecular biology. Reversing the flow of information from RNA to DNA or chromatin has come to the fore in recent years, from the convergence of fundamental discoveries and synthetic biology. Inspired by the example of long noncoding RNAs (lncRNAs) in mammalian genomes that direct chromatin modifications and gene expression, synthetic biologists have repurposed prokaryotic RNA-guided genome defense systems such as CRISPR to edit eukaryotic genomes and epigenomes. Here we explore the parallels of these two fields and highlight opportunities for synergy and future breakthroughs.

Employing template-directed CRISPR-based editing of the OsALS gene to create herbicide tolerance in Basmati rice

Rice (Oryza sativa) is one of the primary food crops which contributes major portion of daily calorie intake. It is used as model crop for various genome editing studies. Basmati rice was also explored for establishing non-homologous end joining-based genome editing. But it was not clear whether homology-directed repair (HDR)-based genome editing can be done in Basmati rice. The current study was designed to establish HDR-based genome editing in Basmati rice to develop herbicide tolerance. There is severe weed spread when rice is grown via direct planted rice method in various countries to save labour and water resources. Therefore, the use of herbicides is necessary to control weeds. These herbicides can also affect cultivated rice which creates the need to develop herbicide-tolerant rice. In current study, we introduced a point mutation in Acetolactate Synthase gene to convert tryptophan to leucine at position 548. For this purpose, different constructs for HDR were tested with different RNA scaffold and orientation of repair templates. Out of four different architectures, the one having repair template identical to the target DNA strand precisely edited the target site. We successfully established template-directed CRISPR-Cas9 system in Super Basmati rice by detecting desired substitutions at the target site in Acetolactate Synthase locus. Moreover, this editing of Acetolactate Synthase gene resulted in the production of herbicide tolerance in Super Basmati rice. This study suggests that such type of HDR system can be used to precisely edit other genes for crop improvement.

Impact of Salmonella genome rearrangement on gene expression

In addition to nucleotide variation, many bacteria also undergo changes at a much larger scale via rearrangement of their genome structure (GS) around long repeat sequences. These rearrangements result in genome fragments shifting position and/or orientation in the genome without necessarily affecting the underlying nucleotide sequence. To date, scalable techniques have not been applied to GS identification, so it remains unclear how extensive this variation is and the extent of its impact upon gene expression. However, the emergence of multiplexed, long-read sequencing overcomes the scale problem, as reads of several thousand bases are routinely produced that can span long repeat sequences to identify the flanking chromosomal DNA, allowing GS identification. Genome rearrangements were generated in Salmonella enterica serovar Typhi through long-term culture at ambient temperature. Colonies with rearrangements were identified via long-range PCR and subjected to long-read nanopore sequencing to confirm genome variation. Four rearrangements were investigated for differential gene expression using transcriptomics. All isolates with changes in genome arrangement relative to the parent strain were accompanied by changes in gene expression. Rearrangements with similar fragment movements demonstrated similar changes in gene expression. The most extreme rearrangement caused a large imbalance between the origin and terminus of replication and was associated with differential gene expression as a factor of distance moved toward or away from the origin of replication. Genome structure variation may provide a mechanism through which bacteria can quickly adapt to new environments and warrants routine assessment alongside traditional nucleotide-level measures of variation.

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.

Precision genome editing in plants using gene targeting and prime editing: existing and emerging strategies

Precise modification of plant genomes, such as seamless insertion, deletion, or replacement of DNA sequences at a predefined site, is a challenging task. Gene targeting (GT) and prime editing are currently the best approaches for this purpose. However, these techniques are inefficient in plants, which limits their applications for crop breeding programs. Recently, substantial developments have been made to improve the efficiency of these techniques in plants. Several strategies, such as RNA donor templating, chemically modified donor DNA template, and tandem-repeat homology-directed repair, are aimed at improving GT. Additionally, improved prime editing gRNA design, use of engineered reverse transcriptase enzymes, and splitting prime editing components have improved the efficacy of prime editing in plants. These emerging strategies and existing technologies are reviewed along with various perspectives on their future improvement and the development of robust precision genome editing technologies for plants.

Timing and characteristics of nuclear events during conjugation and genomic exclusion in Paramecium multimicronucleatum

Ciliated protists are ideal material for studying the origin and evolution of sex, because of their nuclear dimorphism (containing both germline micronucleus and somatic macronucleus in the same cytoplasm), special sexual processes (conjugation and autogamy), and high diversity of mating-type systems. However, the study of sexual process is limited to only a few species, due to the difficulties in inducing or observing conjugation. In the present study, we investigate the conjugation process in Paramecium multimicronucleatum: (1) of the three prezygotic divisions, all micronuclei undergo the first two divisions (meiosis I, II), while a variable number of nuclei undergo the third division (mitosis); (2) the synkaryon divides three times after fertilization, giving rise to eight products that differentiate into four macronuclear anlagen and four micronuclei; (3) cells restore the vegetative stage after two successive cell fissions during which the macronuclear anlagen are distributed into daughter cells without division, while micronuclei divide mitotically; (4) the parental macronucleus begins to fragment following the first meiotic division and finally degenerates completely; (5) the entire process takes about 110 h, of which about 85 h are required for macronuclear development. In addition, we describe for the first time the process of genomic exclusion occurring between amicronucleate and micronucleate cells of P. multimicronucleatum, during which the micronucleate cell contributes a pronucleus to the amicronucleate cell, resulting in both exconjugants being homozygotes. These results provide new insights into the diversity of sexual processes and lay an important cytological basis for future in-depth studies of mating systems in ciliates.

Identification of novel, functional, long noncoding RNAs involved in programmed, large-scale genome rearrangements

Noncoding RNAs (ncRNAs) make up to ∼98% percent of the transcriptome of a given organism. In recent years, one relatively new class of ncRNAs, long noncoding RNAs (lncRNAs), were shown to be more than mere by-products of gene expression and regulation. The unicellular eukaryote Paramecium tetraurelia is a member of the ciliate phylum, an extremely heterogeneous group of organisms found in most bodies of water across the globe. A hallmark of ciliate genetics is nuclear dimorphism and programmed elimination of transposons and transposon-derived DNA elements, the latter of which is essential for the maintenance of the somatic genome. Paramecium and ciliates in general harbor a plethora of different ncRNA species, some of which drive the process of large-scale genome rearrangements, including DNA elimination, during sexual development. Here, we identify and validate the first known functional lncRNAs in ciliates to date. Using deep-sequencing and subsequent bioinformatic processing and experimental validation, we show that Paramecium expresses at least 15 lncRNAs. These candidates were predicted by a highly conservative pipeline, and informatic analyses hint at differential expression during development. Depletion of two lncRNAs, lnc1 and lnc15, resulted in clear phenotypes, decreased survival, morphological impairment, and a global effect on DNA elimination.

Chimeric RNAs Discovered by RNA Sequencing and Their Roles in Cancer and Rare Genetic Diseases

Chimeric RNAs are transcripts that are generated by gene fusion and intergenic splicing events, thus comprising nucleotide sequences from different parental genes. In the past, Northern blot analysis and RT-PCR were used to detect chimeric RNAs. However, they are low-throughput and can be time-consuming, labor-intensive, and cost-prohibitive. With the development of RNA-seq and transcriptome analyses over the past decade, the number of chimeric RNAs in cancer as well as in rare inherited diseases has dramatically increased. Chimeric RNAs may be potential diagnostic biomarkers when they are specifically expressed in cancerous cells and/or tissues. Some chimeric RNAs can also play a role in cell proliferation and cancer development, acting as tools for cancer prognosis, and revealing new insights into the cell origin of tumors. Due to their abilities to characterize a whole transcriptome with a high sequencing depth and intergenically identify spliced chimeric RNAs produced with the absence of chromosomal rearrangement, RNA sequencing has not only enhanced our ability to diagnose genetic diseases, but also provided us with a deeper understanding of these diseases. Here, we reviewed the mechanisms of chimeric RNA formation and the utility of RNA sequencing for discovering chimeric RNAs in several types of cancer and rare inherited diseases. We also discussed the diagnostic, prognostic, and therapeutic values of chimeric RNAs.

Paramecium epigenetics in development and proliferation

The term epigenetics is used for any layer of genetic information aside from the DNA base-sequence information. Mammalian epigenetic research increased our understanding of chromatin dynamics in terms of cytosine methylation and histone modification during differentiation, aging, and disease. Instead, ciliate epigenetics focused more on small RNA-mediated effects. On the one hand, these do concern the transport of RNA from parental to daughter nuclei, representing a regulated transfer of epigenetic information across generations. On the other hand, studies of Paramecium, Tetrahymena, Oxytricha, and Stylonychia revealed an almost unique function of transgenerational RNA. Rather than solely controlling chromatin dynamics, they control sexual progeny's DNA content quantitatively and qualitatively. Thus epigenetics seems to control genetics, at least genetics of the vegetative macronucleus. This combination offers ciliates, in particular, an epigenetically controlled genetic variability. This review summarizes the epigenetic mechanisms that contribute to macronuclear heterogeneity and relates these to nuclear dimorphism. This system's adaptive and evolutionary possibilities raise the critical question of whether such a system is limited to unicellular organisms or binuclear cells. We discuss here the relevance of ciliate genetics and epigenetics to multicellular organisms.

An update on precision genome editing by homology-directed repair in plants

Beneficial alleles derived from local landraces or related species, or even orthologs from other plant species, are often caused by differences of one or several single-nucleotide polymorphisms or indels in either the promoter region or the encoding region of a gene and often account for major differences in agriculturally important traits. Clustered regularly interspaced short palindromic repeats-associated endonuclease Cas9 system (CRISPR/Cas9)-mediated precision genome editing enables targeted allele replacement or insertion of flag or foreign genes at specific loci via homology-directed repair (HDR); however, HDR efficiency is low due to the intrinsic rare occurrence of HDR and insufficient DNA repair template in the proximity of a double-stranded break (DSB). Precise replacement of the targeted gene with elite alleles from landraces or relatives into a commercial variety through genome editing has been a holy grail in the crop genome editing field. In this update, we briefly summarize CRISPR/Cas-mediated HDR in plants. We describe diverse strategies to improve HDR efficiency by manipulating the DNA repair pathway, timing DSB induction, and donor delivery, and so on. Lastly, we outline open questions and challenges in HDR-mediated precision genome editing in both plant biological research and crop improvement.

Membrane Traffic and Ca2+ -Signals in Ciliates

A Paramecium cell has as many types of membrane interactions as mammalian cells, as established with monoclonal antibodies by R. Allen and A. Fok. Since then, we have identified key-players, such as SNARE-proteins, Ca²⁺ -regulating proteins, including Ca²⁺ -channels, Ca²⁺ -pumps, Ca²⁺ -binding proteins of different affinity etc. at the molecular level, probed their function and localized them at the light and electron microscopy level. SNARE-proteins, in conjunction with a synaptotagmin-like Ca²⁺ -sensor protein, mediate membrane fusion. This interaction is additionally regulated by monomeric GTPases whose spectrum in Tetrahymena and Paramecium has been established by A. Turkewitz. As known from mammalian cells, GTPases are activated on membranes in conjunction with lumenal acidification by an H⁺ -ATPase. For these complex molecules we found in Paramecium an unsurpassed number of 17 a-subunit paralogs which connect the polymeric head and basis part, V1 and V0. (This multitude may reflect different local functional requirements.) Together with plasmalemmal Ca²⁺ -influx-channels, locally enriched intracellular InsP₃ -type (InsP₃ R, mainly in osmoregulatory system) and ryanodine receptor-like Ca²⁺ -release channels (ryanodine receptor-like proteins, RyR-LP), this complexity mediates Ca²⁺ signals for most flexible local membrane-to-membrane interactions. As we found, the latter channel types miss a substantial portion of the N-terminal part. Caffeine and 4-chloro-meta-cresol (the agent used to probe mutations of RyRs in man during surgery in malignant insomnia patients) initiate trichocyst exocytosis by activating Ca²⁺ -release channels type CRC-IV in the peripheral part of alveolar sacs. This is superimposed by Ca²⁺ -influx, i.e. a mechanism called "store-operated Ca²⁺ -entry" (SOCE). For the majority of key players, we have mapped paralogs throughout the Paramecium cell, with features in common or at variance in the different organelles participating in vesicle trafficking. Local values of free Ca²⁺ -concentration, [Ca²⁺ ]_i , and their change, e.g. upon exocytosis stimulation, have been registered by flurochromes and chelator effects. In parallel we have registered release of Ca²⁺ from alveolar sacs by quenched-flow analysis combined with cryofixation and x-ray microanalysis.

Experimental Evolution in Tetrahymena

Experimental evolution has provided novel insight into a wide array of biological processes. Species in the genus Tetrahymena are proving to be a highly useful system for studying a range of questions using experimental evolution. Their unusual genomic architecture, diversity of life history traits, importance as both predator and prey, and amenability to laboratory culture allow them to be studied in a variety of contexts. In this paper, we review what we are learning from experimental evolution with Tetrahymena about mutation, adaptation, and eco-evolutionary dynamics. We predict that future experimental evolution studies using Tetrahyemena will continue to shed new light on these processes.

An RNA-Based Theory of Natural Universal Computation

Life is confronted with computation problems in a variety of domains including animal behavior, single-cell behavior, and embryonic development. Yet we currently do not know of a naturally existing biological system that is capable of universal computation, i.e., Turing-equivalent in scope. Generic finite-dimensional dynamical systems (which encompass most models of neural networks, intracellular signaling cascades, and gene regulatory networks) fall short of universal computation, but are assumed to be capable of explaining cognition and development. I present a class of models that bridge two concepts from distant fields: combinatory logic (or, equivalently, lambda calculus) and RNA molecular biology. A set of basic RNA editing rules can make it possible to compute any computable function with identical algorithmic complexity to that of Turing machines. The models do not assume extraordinarily complex molecular machinery or any processes that radically differ from what we already know to occur in cells. Distinct independent enzymes can mediate each of the rules and RNA molecules solve the problem of parenthesis matching through their secondary structure. In the most plausible of these models all of the editing rules can be implemented with merely cleavage and ligation operations at fixed positions relative to predefined motifs. This demonstrates that universal computation is well within the reach of molecular biology. It is therefore reasonable to assume that life has evolved - or possibly began with - a universal computer that yet remains to be discovered. The variety of seemingly unrelated computational problems across many scales can potentially be solved using the same RNA-based computation system. Experimental validation of this theory may immensely impact our understanding of memory, cognition, development, disease, evolution, and the early stages of life.

RNA-driven JAZF1-SUZ12 gene fusion in human endometrial stromal cells

Oncogenic fusion genes as the result of chromosomal rearrangements are important for understanding genome instability in cancer cells and developing useful cancer therapies. To date, the mechanisms that create such oncogenic fusion genes are poorly understood. Previously we reported an unappreciated RNA-driven mechanism in human prostate cells in which the expression of chimeric RNA induces specified gene fusions in a sequence-dependent manner. One fundamental question yet to be addressed is whether such RNA-driven gene fusion mechanism is generalizable, or rather, a special case restricted to prostate cells. In this report, we demonstrated that the expression of designed chimeric RNAs in human endometrial stromal cells leads to the formation of JAZF1-SUZ12, a cancer fusion gene commonly found in low-grade endometrial stromal sarcomas. The process is specified by the sequence of chimeric RNA involved and inhibited by estrogen or progesterone. Furthermore, it is the antisense rather than sense chimeric RNAs that effectively drive JAZF1-SUZ12 gene fusion. The induced fusion gene is validated both at the RNA and the genomic DNA level. The ability of designed chimeric RNAs to drive and recapitulate the formation of JAZF1-SUZ12 gene fusion in endometrial cells represents another independent case of RNA-driven gene fusion, suggesting that RNA-driven genomic recombination is a permissible mechanism in mammalian cells. The results could have fundamental implications in the role of RNA in genome stability, and provide important insight in early disease mechanisms related to the formation of cancer fusion genes.

Fusion genes resulting from chromosomal translocations are important for understanding cancer mechanisms and developing anti-cancer therapies. Fusion gene are presumed to occur prior to fusion RNA expression. However, studies have reported the presence of fusion RNAs in individuals who were negative for chromosomal translocations. The observation, that fusion RNA could be present prior to fusion gene formation, raises the provocative hypothesis that fusion RNA, or any cellular RNA with sequence compositions resembling that of fusion RNA, could act as a template to mediate genomic rearrangement which leads to the final gene fusion. In this report, we demonstrated that the expression of designed chimeric RNAs in human endometrial stromal cells leads to the formation of JAZF1-SUZ12, a cancer fusion gene found in endometrial stromal sarcomas. The process is specified by the sequence of chimeric RNA involved and inhibited by estrogen or progesterone. Furthermore, it is the antisense rather than sense chimeric RNAs that effectively drive JAZF1-SUZ12 gene fusion. The results could have fundamental implications in the role of RNA in mammalian genome stability, provide important insight in early disease mechanism, as well as developing gene editing technology via mechanisms native to mammalian cells.

Transcription-coupled DNA double-strand break repair

The genomic DNA is constantly under attack by cellular and/or environmental factors. Fortunately, the cell is armed to safeguard its genome by various mechanisms such as nucleotide excision, base excision, mismatch and DNA double-strand break repairs. While these processes maintain the integrity of the genome throughout, DNA repair occurs preferentially faster at the transcriptionally active genes. Such transcription-coupled repair phenomenon plays important roles to maintain active genome integrity, failure of which would interfere with transcription, leading to an altered gene expression (and hence cellular pathologies/diseases). Among the various DNA damages, DNA double-strand breaks are quite toxic to the cells. If DNA double-strand break occurs at the active gene, it would interfere with transcription/gene expression, thus threatening cellular viability. Such DNA double-strand breaks are found to be repaired faster at the active gene in comparison to its inactive state or the inactive gene, thus supporting the existence of a new phenomenon of transcription-coupled DNA double-strand break repair. Here, we describe the advances of this repair process.

m6A regulators are associated with osteosarcoma metastasis and have prognostic significance: A study based on public databases

BACKGROUND

Osteosarcoma represents the most common malignant bone tumor with high metastatic potential and inferior prognosis. RNA methylation (N6-methyladenosine [m6A]) is a prevalent RNA modification that epigenetically influences numerous biological processes including tumorigenesis. This study aims to determine that m6A regulators are significant biomarkers for osteosarcoma, and establish a prognostic model to predict the survival of patients.

METHODS

In this study, we comprehensively analyzed the underlying associations between m6A regulators' mRNA expressions and metastasis as well as prognosis of osteosarcoma patients in the Cancer Genome Atlas. Multivariate Cox-regression analysis was used to screen regulators that were significantly associated with overall survival of osteosarcoma patients. Least absolute shrinkage and selection operator (LASSO) Cox-regression analysis was used for constructing m6A regulator-based osteosarcoma prognostic signature.

RESULTS

Some of the regulators exhibited aberrant mRNA levels between osteosarcoma samples with and without metastasis. Multivariate Cox-regression analysis identified several regulators with potential prognostic significance. A risk score formula consisted of methyltransferase-like 3, YTH domains of Homo sapiens, and fat mass and obesity-associated protein was obtained through which patients could be prognostically stratified independently of potential confounding factors. The signature was also significantly associated with the metastatic potential of osteosarcoma. All the analyses could be well reproduced in another independent osteosarcoma cohort from the Gene Expression Omnibus.

CONCLUSIONS

In conclusion, this study first revealed potential roles of m6A regulators in osteosarcoma metastasis and prognosis, which should be helpful for its clinical decision-making.

The Role of RNA in DNA Breaks, Repair and Chromosomal Rearrangements

Incorrect reparation of DNA double-strand breaks (DSB) leading to chromosomal rearrangements is one of oncogenesis's primary causes. Recently published data elucidate the key role of various types of RNA in DSB formation, recognition and repair. With growing interest in RNA biology, increasing RNAs are classified as crucial at the different stages of the main pathways of DSB repair in eukaryotic cells: nonhomologous end joining (NHEJ) and homology-directed repair (HDR). Gene mutations or variation in expression levels of such RNAs can lead to local DNA repair defects, increasing the chromosome aberration frequency. Moreover, it was demonstrated that some RNAs could stimulate long-range chromosomal rearrangements. In this review, we discuss recent evidence demonstrating the role of various RNAs in DSB formation and repair. We also consider how RNA may mediate certain chromosomal rearrangements in a sequence-specific manner.

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.

RNA Interference by Cyanobacterial Feeding Demonstrates the SCSG1 Gene Is Essential for Ciliogenesis during Oral Apparatus Regeneration in Stentor

In the giant ciliate Stentor coeruleus, oral apparatus (OA) regeneration is an experimentally tractable regeneration paradigm that occurs via a series of morphological steps. OA regeneration is thought to be driven by a complex regulatory system that orchestrates the temporal expression of conserved and specific genes. We previously identified a S. coeruleus-specific gene (named SCSG1) that was significantly upregulated during the ciliogenesis stages of OA regeneration, with an expression peak at the stage of the first OA cilia appearance. We established a novel RNA interference (RNAi) method through cyanobacteria Synechocystis sp. PCC6803 feeding in S. coeruleus. The expression of SCSG1 gene was significantly knocked down by using this method and induced abnormal ciliogenesis of OA regeneration in S. coeruleus, suggesting that SCSG1 is essential for OA regeneration in S. coeruleus. This novel RNAi method by cyanobacterial feeding has potential utility for studying other ciliates.

Programmed genome rearrangements in ciliates

Ciliates are a highly divergent group of unicellular eukaryotes with separate somatic and germline genomes found in distinct dimorphic nuclei. This characteristic feature is tightly linked to extremely laborious developmentally regulated genome rearrangements in the development of a new somatic genome/nuclei following sex. The transformation from germline to soma genome involves massive DNA elimination mediated by non-coding RNAs, chromosome fragmentation, as well as DNA amplification. In this review, we discuss the similarities and differences in the genome reorganization processes of the model ciliates Paramecium and Tetrahymena (class Oligohymenophorea), and the distantly related Euplotes, Stylonychia, and Oxytricha (class Spirotrichea).

Functional chimeric genes in ciliates: An instructive case from Euplotes raikovi

In ciliates, with every sexual event the transcriptionally active genes of the sub-chromosomic somatic genome that resides in the cell macronucleus are lost. They are de novo assembled starting from 'Macronuclear Destined Sequences' that arise from the fragmentation of transcriptionally silent DNA sequences of the germline chromosomic genome enclosed in the cell micronucleus. The RNA-mediated epigenetic mechanism that drives the assembly of these sequences is subject to errors which result in the formation of chimeric genes. Studying a gene family that in Euplotes raikovi controls the synthesis of protein signal pheromones responsible for a self/not-self recognition mechanism, we identified the chimeric structure of an 851-bp macronuclear gene previously known to specify soluble and membrane-bound pheromone molecules through an intron-splicing mechanism. This chimeric gene, designated mac-er-1*, conserved the native pheromone-gene structure throughout its coding and 3' regions. Instead, its 5' region is completely unrelated to the pheromone gene structure at the level of a 360-bp sequence, which derives from the assembly with a MDS destined to compound a 2417-bp gene encoding a 696-amino acid protein with unknown function. This mac-er-1* gene characterization provides further evidence that ciliates rely on functional chimeric genes that originate in non-programmed phenomena of somatic MDS recombination to increase the species genetic variability independently of gene reshuffling phenomena of the germline genome.

Genetic Characterization of Three Distinct Mechanisms Supporting RNA-Driven DNA Repair and Modification Reveals Major Role of DNA Polymerase ζ

DNA double-stranded breaks (DSBs) are dangerous lesions threatening genomic stability. Fidelity of DSB repair is best achieved by recombination with a homologous template sequence. In yeast, transcript RNA was shown to template DSB repair of DNA. However, molecular pathways of RNA-driven repair processes remain obscure. Utilizing assays of RNA-DNA recombination with and without an induced DSB in yeast DNA, we characterize three forms of RNA-mediated genomic modifications: RNA- and cDNA-templated DSB repair (R-TDR and c-TDR) using an RNA transcript or a DNA copy of the RNA transcript for DSB repair, respectively, and a new mechanism of RNA-templated DNA modification (R-TDM) induced by spontaneous or mutagen-induced breaks. While c-TDR requires reverse transcriptase, translesion DNA polymerase ζ (Pol ζ) plays a major role in R-TDR, and it is essential for R-TDM. This study characterizes mechanisms of RNA-DNA recombination, uncovering a role of Pol ζ in transferring genetic information from transcript RNA to DNA.

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.

Relationship among DNA double-strand break (DSB), DSB repair, and transcription prevents genome instability and cancer

DNA double‐strand break (DSB) is a serious type of DNA damage and is known to trigger multiple responses within cells. In these responses, novel relationships among DSB, DSB repair, and transcription machineries are created. First, transcription is repressed if DSB occurs near or at the transcription site, termed DSB‐induced transcriptional repression, which contributes to DSB repair with the aid of DNA damage‐signaling pathways, ATM‐ or DNA‐PKcs‐signaling pathways. DSB‐induced transcriptional repression is also regulated by transcriptional factors TLP1, NELF, and ENL, as well as chromatin remodeling and organizing factors ZMYND8, CDYL1, PBAF, and cohesin. Second, transcription and RNA promote DSB repair for genome integrity. Transcription factors such as LEDGF, SETD2, and transcriptionally active histone modification, H3K36, facilitate homologous recombination to overcome DSB. At transcriptional active sites, DNA:RNA hybrids, termed R‐loops, which are formed by DSB, are processed by RAD52 and XPG leading to an activation of the homologous recombination pathway. Even in a transcriptionally inactive non‐genic sites, noncoding RNAs that are produced by RNA polymerase II, DICER, and DROSHA, help to recruit DSB repair proteins at the DSB sites. Third, transcriptional activation itself, however, can induce DSB. Transcriptional activation often generates specific DNA structures such as R‐loops and topoisomerase‐induced DSBs, which cause genotoxic stress and may lead to genome instability and consequently to cancer. Thus, transcription and DSB repair machineries interact and cooperate to prevent genome instability and cancer.

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.

Conjugation in Euplotes raikovi (Protista, Ciliophora): New Insights into Nuclear Events and Macronuclear Development from Micronucleate and Amicronucleate Cells

Ciliates form a distinct group of single-celled eukaryotes that host two types of nuclei (micro and macronucleus) in the same cytoplasm and have a special sexual process known as conjugation, which involves mitosis, meiosis, fertilization, nuclear differentiation, and development. Due to their high species diversity, ciliates have evolved different patterns of nuclear events during conjugation. In the present study, we investigate these events in detail in the marine species Euplotes raikovi. Our results indicate that: (i) conjugation lasts for about 50 h, the longest stage being the development of the new macronucleus (ca. 36 h); (ii) there are three prezygotic micronuclear divisions (mitosis and meiosis I and II) and two postzygotic synkaryon divisions; and (iii) a fragment of the parental macronucleus fuses with the new developing macronucleus. In addition, we describe for the first time conjugation in amicronucleate E. raikovi cells. When two amicronucleate cells mate, they separate after about 4 h without evident nuclear changes; when one amicronucleate cell mates with a micronucleate cell, the micronucleus undergoes regular prezygotic divisions to form migratory and stationary pronuclei, but the two pronuclei fuse in the same cell. In the amicronucleate cell, the parental macronucleus breaks into fragments, which are then recovered to form a new functional macronucleus. These results add new information on the process of conjugation in both micronucleate and amicronucleate Euplotes cells.

Precise gene replacement in plants through CRISPR/Cas genome editing technology: current status and future perspectives

CRISPR/Cas, as a simple, versatile, robust and cost-effective system for genome manipulation, has dominated the genome editing field over the past few years. The application of CRISPR/Cas in crop improvement is particularly important in the context of global climate change, as well as diverse agricultural, environmental and ecological challenges. Various CRISPR/Cas toolboxes have been developed and allow for targeted mutagenesis at specific genome loci, transcriptome regulation and epigenome editing, base editing, and precise targeted gene/allele replacement or tagging in plants. In particular, precise replacement of an existing allele with an elite allele in a commercial variety through homology-directed repair (HDR) is a holy grail in genome editing for crop improvement as it has been very difficult, laborious and time-consuming to introgress the elite alleles into commercial varieties without any linkage drag from parental lines within a few generations in crop breeding practice. However, it still remains very challenging in crop plants. This review intends to provide an informative summary of the latest development and breakthroughs in gene replacement using CRISPR/Cas technology, with a focus on achievements, potential mechanisms and future perspectives in plant biological science as well as crop improvement.

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.

Should dsRNA treatments applied in outdoor environments be regulated?

The New Zealand Environmental Protection Authority (EPA) issued a Decision that makes the use of externally applied double-stranded (ds)RNA molecules on eukaryotic cells or organisms technically out of scope of legislation on new organisms, making risk assessments of such treatments in the open environment unnecessary. The Decision was based on its view that the treatment does not create new or genetically modified organisms and rests on the EPA's conclusions that dsRNA is not heritable and is not a mutagen. For these reasons EPA decided that treatments using dsRNA do not modify genes or other genetic material. I found from an independent review of the literature on the topic indicated, however, that each of the major scientific justifications relied upon by the EPA was based on either an inaccurate interpretation of evidence or failure to consult the research literature pertaining to additional types of eukaryotes. The Decision also did not take into account the unknown and unique eukaryotic biodiversity of New Zealand. The safe use of RNA-based technology holds promise for addressing complex and persistent challenges in public health, agriculture and conservation. However, by failing to restrict the source or means of modifying the dsRNA, the EPA removed regulatory oversight that could prevent unintended consequences of this new technology such as suppression of genes other than those selected for suppression or the release of viral genes or genomes by failing to restrict the source or means of modifying the dsRNA.

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.

The Role of Noncoding RNAs in Double-Strand Break Repair

Genome stability is constantly threatened by DNA lesions generated by different environmental factors as well as endogenous processes. If not properly and timely repaired, damaged DNA can lead to mutations or chromosomal rearrangements, well-known reasons for genetic diseases or cancer in mammals, or growth abnormalities and/or sterility in plants. To prevent deleterious consequences of DNA damage, a sophisticated system termed DNA damage response (DDR) detects DNA lesions and initiates DNA repair processes. In addition to many well-studied canonical proteins involved in this process, noncoding RNA (ncRNA) molecules have recently been discovered as important regulators of the DDR pathway, extending the broad functional repertoire of ncRNAs to the maintenance of genome stability. These ncRNAs are mainly connected with double-strand breaks (DSBs), the most dangerous type of DNA lesions. The possibility to intentionally generate site-specific DSBs in the genome with endonucleases constitutes a powerful tool to study, in vivo, how DSBs are processed and how ncRNAs participate in this crucial event. In this review, we will summarize studies reporting the different roles of ncRNAs in DSB repair and discuss how genome editing approaches, especially CRISPR/Cas systems, can assist DNA repair studies. We will summarize knowledge concerning the functional significance of ncRNAs in DNA repair and their contribution to genome stability and integrity, with a focus on plants.

Gene fusions and chimeric RNAs, and their implications in cancer

Gene fusions are appreciated as ideal cancer biomarkers and therapeutic targets. Chimeric RNAs are traditionally thought to be products of gene fusions, and thus, also cancer-specific. Recent research has demonstrated that chimeric RNAs can be generated by intergenic splicing in the absence of gene fusion, and such chimeric RNAs are also found in normal physiology. These new findings challenge the traditional theory of chimeric RNAs exclusivity to cancer, and complicates use of chimeric RNAs in cancer detection. Here, we provide an overview of gene fusions and chimeric RNAs, and emphasize their differences. We note that gene fusions are able to generate chimeric RNAs in accordance with the central dogma of biology, and that chimeric RNAs may also be able to influence the generation of the gene fusions per the “horse before the cart” hypothesis. We further expand upon the “horse before the cart” hypothesis, summarizing current evidence in support of the theory and exploring its potential impact on the field.

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.

Evolutionary origins and impacts of genome architecture in ciliates

Genome architecture is well diversified among eukaryotes in terms of size and content, with many being radically shaped by ancient and ongoing genome conflicts with transposable elements (e.g., the large transposon‐rich genomes common among plants). In ciliates, a group of microbial eukaryotes with distinct somatic and germ‐line genomes present in a single cell, the consequences of these genome conflicts are most apparent in their developmentally programmed genome rearrangements. This complicated developmental phenomenon has largely overshadowed and outpaced our understanding of how germ‐line and somatic genome architectures have influenced the evolutionary dynamism and potential in these taxa. In our review, we highlight three central concepts: how the evolution of atypical ciliate germ‐line genome architectures is linked to ancient genome conflicts; how the complex, epigenetically guided transformation of germline to soma during development can generate widespread genetic variation; and how these features, coupled with their unusual life cycle, have increased the rate of molecular evolution linked to genome architecture in these taxa.

Estimating Strengths of Individual Hydrogen Bonds in RNA Base Pairs: Toward a Consensus between Different Computational Approaches

Noncoding RNA molecules are composed of a large variety of noncanonical base pairs that shape up their functionally competent folded structures. Each base pair is composed of at least two interbase hydrogen bonds (H-bonds). It is expected that the characteristic geometry and stability of different noncanonical base pairs are determined collectively by the properties of these interbase H-bonds. We have studied the ground-state electronic properties [using density functional theory (DFT) and DFT-D3-based methods] of all the 118 normal base pairs and 36 modified base pairs, belonging to 12 different geometric families (cis and trans of WW, WH, HH, WS, HS, and SS) that occur in a nonredundant set of high-resolution RNA crystal structures. Having addressed some of the limitations of the earlier approaches, we provide here a comprehensive compilation of the average energies of different types of interbase H-bonds (E _HB). We have also characterized each interbase H-bond using 13 different parameters that describe its geometry, charge distribution at its bond critical point (BCP), and n → σ*-type charge transfer from filled π orbitals of the H-bond acceptor to the empty antibonding orbital of the H-bond donor. On the basis of the extent of their linear correlation with the H-bonding energy, we have shortlisted five parameters to model linear equations for predicting E _HB values. They are (i) electron density at the BCP: ρ, (ii) its Laplacian: ∇²ρ, (iii) stabilization energy due to n → σ*-type charge transfer: E(2), (iv) donor-hydrogen distance, and (v) hydrogen-acceptor distance. We have performed single variable and multivariable linear regression analysis over the normal base pairs and have modeled sets of linear relationships between these five parameters and E _HB. Performance testing of our model over the set of modified base pairs shows promising results, at least for the moderately strong H-bonds.

Precise gene replacement in rice by RNA transcript-templated homologous recombination

One of the main obstacles to gene replacement in plants is efficient delivery of a donor repair template (DRT) into the nucleus for homology-directed DNA repair (HDR) of double-stranded DNA breaks. Production of RNA templates in vivo for transcript-templated HDR (TT-HDR) could overcome this problem, but primary transcripts are often processed and transported to the cytosol, rendering them unavailable for HDR. We show that coupling CRISPR-Cpf1 (CRISPR from Prevotella and Francisella 1) to a CRISPR RNA (crRNA) array flanked with ribozymes, along with a DRT flanked with either ribozymes or crRNA targets, produces primary transcripts that self-process to release the crRNAs and DRT inside the nucleus. We replaced the rice acetolactate synthase gene (ALS) with a mutated version using a DNA-free ribonucleoprotein complex that contains the recombinant Cpf1, crRNAs, and DRT transcripts. We also produced stable lines with two desired mutations in the ALS gene using TT-HDR.

The prognostic potential of long noncoding RNA HOTAIR expression in human digestive system carcinomas: A meta-analysis

Homeobox transcript antisense intergenic RNA (HOTAIR), one of the well-known long noncoding RNAs (lncRNAs), plays an important role in initiation and development of various tumors. Elevated level of HOTAIR is associated with metastatic behavior of primary tumor and poor outcome in several cancers. Therefore, we conducted a meta-analysis to clearly measure the prognostic impact of HOTAIR in patients with digestive system carcinomas. Fourteen studies including 2,666 patients with five different type of digestive system cancers were selected to be entered in meta-analysis. Finding demonstrated that HOTAIR overexpression could predict unfavorable outcome in digestive system carcinomas (hazard ratio [HR] = 2.4, 95% confidence interval [CI]: 2.0-2.9; p < 0.001; fixed-effect model). In stratified analysis, increased level of HOTAIR predicted poor overall survival in gastric cancer (HR = 2.1, 95% CI: 1.6-2.9; p < 0.001), colorectal cancer (HR = 4.1, 95% CI: 1.6-10.2; p = 0.002), esophageal squamous cell carcinoma (HR = 2.3, 95% CI: 1.7-3.0; p < 0.001), and hepatocellular carcinoma (HR = 3.4, 95% CI: 1.9-6.1; p < 0.001). Our meta-analysis results clearly support the prognostic value of HOTAIR to predict unfavorable prognostic outcomes in diverse digestive system carcinomas.

RNA-mediated gene fusion in mammalian cells

One of the hallmarks of cancer is the formation of oncogenic fusion genes as a result of chromosomal translocations. Fusion genes are presumed to form before fusion RNA expression. However, studies have reported the presence of fusion RNAs in individuals who were negative for chromosomal translocations. These observations give rise to "the cart before the horse" hypothesis, in which the genesis of a fusion RNA precedes the fusion gene. The fusion RNA then guides the genomic rearrangements that ultimately result in a gene fusion. However, RNA-mediated genomic rearrangements in mammalian cells have never been demonstrated. Here we provide evidence that expression of a chimeric RNA drives formation of a specified gene fusion via genomic rearrangement in mammalian cells. The process is: (i) specified by the sequence of chimeric RNA involved, (ii) facilitated by physiological hormone levels, (iii) permissible regardless of intrachromosomal (TMPRSS2-ERG) or interchromosomal (TMPRSS2-ETV1) fusion, and (iv) can occur in normal cells before malignant transformation. We demonstrate that, contrary to "the cart before the horse" model, it is the antisense rather than sense chimeric RNAs that effectively drive gene fusion, and that this disparity can be explained by transcriptional conflict. Furthermore, we identified an endogenous RNA AZI1 that functions as the "initiator" RNA to induce TMPRSS2-ERG fusion. RNA-driven gene fusion demonstrated in this report provides important insight in early disease mechanisms, and could have fundamental implications in the biology of mammalian genome stability, as well as gene-editing technology via mechanisms native to mammalian cells.