Transcriptome-wide expansion of non-coding regulatory switches: evidence from co-occurrence of Alu exonization, antisense and editing

Recent insights into crosstalk between genetic parasites and their host genome

The bulk of higher order organismal genomes is comprised of transposable element (TE) copies, i.e. genetic parasites. The host-parasite relation is multi-faceted, varying across genomic region (genic versus intergenic), life-cycle stages, tissue-type and of course in health versus pathological state. The reach of functional genomics though, in investigating genotype-to-phenotype relations, has been limited when TEs are involved. The aim of this review is to highlight recent progress made in understanding how TE origin biochemical activity interacts with the central dogma stages of the host genome. Such interaction can also bring about modulation of the immune context and this could have important repercussions in disease state where immunity has a role to play. Thus, the review is to instigate ideas and action points around identifying evolutionary adaptations that the host genome and the genetic parasite have evolved and why they could be relevant.

Alu RNA fold links splicing with signal recognition particle proteins

Transcriptomic diversity in primates was considerably expanded by exonizations of intronic Alu elements. To better understand their cellular mechanisms we have used structure-based mutagenesis coupled with functional and proteomic assays to study the impact of successive primate mutations and their combinations on inclusion of a sense-oriented AluJ exon in the human F8 gene. We show that the splicing outcome was better predicted by consecutive RNA conformation changes than by computationally derived splicing regulatory motifs. We also demonstrate an involvement of SRP9/14 (signal recognition particle) heterodimer in splicing regulation of Alu-derived exons. Nucleotide substitutions that accumulated during primate evolution relaxed the conserved left-arm AluJ structure including helix H1 and reduced the capacity of SRP9/14 to stabilize the closed Alu conformation. RNA secondary structure-constrained mutations that promoted open Y-shaped conformations of the Alu made the Alu exon inclusion reliant on DHX9. Finally, we identified additional SRP9/14 sensitive Alu exons and predicted their functional roles in the cell. Together, these results provide unique insights into architectural elements required for sense Alu exonization, identify conserved pre-mRNA structures involved in exon selection and point to a possible chaperone activity of SRP9/14 outside the mammalian signal recognition particle.

Exonized Alu repeats in the 3'UTR of a CYP20A1_Alu-LT transcript act as a miRNA sponge

OBJECTIVE

Alu repeats have gained huge importance in the creation and modification of regulatory networks. We previously reported a unique isoform of human CYP20A1 i.e. CYP20A1_Alu-LT with 23 Alu repeats exonized in its 9 kb long 3'UTR with 4742 potential binding sites for 994 miRNAs. The role of this transcript was hypothesized as a potential miRNA sponge in primary neurons as its expression correlated with that of 380 genes having shared miRNA sites and enriched in neuro-coagulopathy. This study provides experimental evidence for the miRNA sponge activity of CYP20A1_Alu-LT in neuronal cell lines.

RESULTS

We studied the Alu-rich fragment of the CYP20A1_Alu-LT extended 3'UTR with > 10 binding sites for miR-619-5p and miR-3677-3p. Enrichment of the Alu-rich fragment with Ago2 confirmed miRNA association of this transcript. Cloning the fragment downstream of a reporter gene led to a 90% decrease in luciferase activity. Overexpression and knockdown studies revealed a positive correlation between the expression of CYP20A1_Alu-LT and miR-619-5p / miR-3677-3p target genes. GAP43, one of the key modulators of nerve regeneration, was significantly altered by the expression of CYP20A1_Alu-LT. This study, for the first time, provides evidence for a unique regulatory function of exonized Alu repeats as miRNA sponges.

The Role of Transposable Elements of the Human Genome in Neuronal Function and Pathology

Transposable elements (TEs) have been extensively studied for decades. In recent years, the introduction of whole-genome and whole-transcriptome approaches, as well as single-cell resolution techniques, provided a breakthrough that uncovered TE involvement in host gene expression regulation underlying multiple normal and pathological processes. Of particular interest is increased TE activity in neuronal tissue, and specifically in the hippocampus, that was repeatedly demonstrated in multiple experiments. On the other hand, numerous neuropathologies are associated with TE dysregulation. Here, we provide a comprehensive review of literature about the role of TEs in neurons published over the last three decades. The first chapter of the present review describes known mechanisms of TE interaction with host genomes in general, with the focus on mammalian and human TEs; the second chapter provides examples of TE exaptation in normal neuronal tissue, including TE involvement in neuronal differentiation and plasticity; and the last chapter lists TE-related neuropathologies. We sought to provide specific molecular mechanisms of TE involvement in neuron-specific processes whenever possible; however, in many cases, only phenomenological reports were available. This underscores the importance of further studies in this area.

The 'Alu-ome' shapes the epigenetic environment of regulatory elements controlling cellular defense

Promoters and enhancers are sites of transcription initiation (TSSs) and carry specific histone modifications, including H3K4me1, H3K4me3, and H3K27ac. Yet, the principles governing the boundaries of such regulatory elements are still poorly characterized. Alu elements are good candidates for a boundary function, being highly abundant in gene-rich regions, while essentially excluded from regulatory elements. Here, we show that the interval ranging from TSS to first upstream Alu, accommodates all H3K4me3 and most H3K27ac marks, while excluding DNA methylation. Remarkably, the average length of these intervals greatly varies in-between tissues, being longer in stem- and shorter in immune-cells. The very shortest TSS-to-first-Alu intervals were observed at promoters active in T-cells, particularly at immune genes, where first-Alus were traversed by RNA polymerase II transcription, while accumulating H3K4me1 signal. Finally, DNA methylation at first-Alus was found to evolve with age, regressing from young to middle-aged, then recovering later in life. Thus, the first-Alus upstream of TSSs appear as dynamic boundaries marking the transition from DNA methylation to active histone modifications at regulatory elements, while also participating in the recording of immune gene transcriptional events by positioning H3K4me1-modified nucleosomes.

SINE jumping contributes to large-scale polymorphisms in the pig genomes

BACKGROUND

Molecular markers based on retrotransposon insertion polymorphisms (RIPs) have been developed and are widely used in plants and animals. Short interspersed nuclear elements (SINEs) exert wide impacts on gene activity and even on phenotypes. However, SINE RIP profiles in livestock remain largely unknown, and not be revealed in pigs.

RESULTS

Our data revealed that SINEA1 displayed the most polymorphic insertions (22.5 % intragenic and 26.5 % intergenic), followed by SINEA2 (10.5 % intragenic and 9 % intergenic) and SINEA3 (12.5 % intragenic and 5.0 % intergenic). We developed a genome-wide SINE RIP mining protocol and obtained a large number of SINE RIPs (36,284), with over 80 % accuracy and an even distribution in chromosomes (14.5/Mb), and 74.34 % of SINE RIPs generated by SINEA1 element. Over 65 % of pig SINE RIPs overlap with genes, most of them (> 95 %) are in introns. Overall, about one forth (23.09 %) of the total genes contain SINE RIPs. Significant biases of SINE RIPs in the transcripts of protein coding genes were observed. Nearly half of the RIPs are common in these pig breeds. Sixteen SINE RIPs were applied for population genetic analysis in 23 pig breeds, the phylogeny tree and cluster analysis were generally consistent with the geographical distributions of native pig breeds in China.

CONCLUSIONS

Our analysis revealed that SINEA1-3 elements, particularly SINEA1, are high polymorphic across different pig breeds, and generate large-scale structural variations in the pig genomes. And over 35,000 SINE RIP markers were obtained. These data indicate that young SINE elements play important roles in creating new genetic variations and shaping the evolution of pig genome, and also provide strong evidences to support the great potential of SINE RIPs as genetic markers, which can be used for population genetic analysis and quantitative trait locus (QTL) mapping in pig.

Antisense RNAs during early vertebrate development are divided in groups with distinct features

Long noncoding RNAs or lncRNAs are a class of non-protein-coding RNAs that are >200 nt in length. Almost 50% of lncRNAs during zebrafish development are transcribed in an antisense direction to a protein-coding gene. However, the role of these natural antisense transcripts (NATs) during development remains enigmatic. To understand NATs in early vertebrate development, we took a computational biology approach and analyzed existing as well as novel data sets. Our analysis indicates that zebrafish NATs can be divided into two major classes based on their coexpression patterns with respect to the overlapping protein-coding genes. Group 1 NATs have characteristics similar to maternally deposited RNAs in that their levels decrease as development progresses. Group 1 NAT levels are negatively correlated with that of overlapping sense-strand protein-coding genes. Conversely, Group 2 NATs are coexpressed with overlapping protein-coding genes. In contrast to Group 1, which is enriched in genes involved in developmental pathways, Group 2 protein-coding genes are enriched in housekeeping functions. Group 1 NATs also show larger overlap and higher complementarity with the sense-strand mRNAs compared to other NATs. In addition, our transcriptomics data, quantifying RNA levels from cytoplasmic and nuclear compartments, indicates that Group 1 NATs are more abundant in the cytosol. Based on their expression pattern, cytosolic nature, and their higher complementarity to the overlapping developmental mRNAs, we speculate that Group 1 NATs function post-transcriptionally to silence spurious expression of developmental genes.

Multiple Alu Exonization in 3'UTR of a Primate-Specific Isoform of CYP20A1 Creates a Potential miRNA Sponge

Alu repeats contribute to phylogenetic novelties in conserved regulatory networks in primates. Our study highlights how exonized Alus could nucleate large-scale mRNA-miRNA interactions. Using a functional genomics approach, we characterize a transcript isoform of an orphan gene, CYP20A1 (CYP20A1_Alu-LT) that has exonization of 23 Alus in its 3'UTR. CYP20A1_Alu-LT, confirmed by 3'RACE, is an outlier in length (9 kb 3'UTR) and widely expressed. Using publically available data sets, we demonstrate its expression in higher primates and presence in single nucleus RNA-seq of 15,928 human cortical neurons. miRanda predicts ∼4,700 miRNA recognition elements (MREs) for ∼1,000 miRNAs, primarily originated within these 3'UTR-Alus. CYP20A1_Alu-LT could be a potential multi-miRNA sponge as it harbors ≥10 MREs for 140 miRNAs and has cytosolic localization. We further tested whether expression of CYP20A1_Alu-LT correlates with mRNAs harboring similar MRE targets. RNA-seq with conjoint miRNA-seq analysis was done in primary human neurons where we observed CYP20A1_Alu-LT to be downregulated during heat shock response and upregulated in HIV1-Tat treatment. In total, 380 genes were positively correlated with its expression (significantly downregulated in heat shock and upregulated in Tat) and they harbored MREs for nine expressed miRNAs which were also enriched in CYP20A1_Alu-LT. MREs were significantly enriched in these 380 genes compared with random sets of differentially expressed genes (P = 8.134e-12). Gene ontology suggested involvement of these genes in neuronal development and hemostasis pathways thus proposing a novel component of Alu-miRNA-mediated transcriptional modulation that could govern specific physiological outcomes in higher primates.

New Insights on the Evolution of Genome Content: Population Dynamics of Transposable Elements in Flies and Humans

Understanding the abundance, diversity, and distribution of TEs in genomes is crucial to understand genome structure, function, and evolution. Advances in whole-genome sequencing techniques, as well as in bioinformatics tools, have increased our ability to detect and analyze the transposable element content in genomes. In addition to reference genomes, we now have access to population datasets in which multiple individuals within a species are sequenced. In this chapter, we highlight the recent advances in the study of TE population dynamics focusing on fruit flies and humans, which represent two extremes in terms of TE abundance, diversity, and activity. We review the most recent methodological approaches applied to the study of TE dynamics as well as the new knowledge on host factors involved in the regulation of TE activity. In addition to transposition rates, we also focus on TE deletion rates and on the selective forces that affect the dynamics of TEs in genomes.

Exaptation at the molecular genetic level

The realization that body parts of animals and plants can be recruited or coopted for novel functions dates back to, or even predates the observations of Darwin. S.J. Gould and E.S. Vrba recognized a mode of evolution of characters that differs from adaptation. The umbrella term aptation was supplemented with the concept of exaptation. Unlike adaptations, which are restricted to features built by selection for their current role, exaptations are features that currently enhance fitness, even though their present role was not a result of natural selection. Exaptations can also arise from nonaptations; these are characters which had previously been evolving neutrally. All nonaptations are potential exaptations. The concept of exaptation was expanded to the molecular genetic level which aided greatly in understanding the enormous potential of neutrally evolving repetitive DNA-including transposed elements, formerly considered junk DNA-for the evolution of genes and genomes. The distinction between adaptations and exaptations is outlined in this review and examples are given. Also elaborated on is the fact that such distinctions are sometimes more difficult to determine; this is a widespread phenomenon in biology, where continua abound and clear borders between states and definitions are rare.

The angiotensin-I-converting enzyme insertion/deletion in polymorphic element codes for an AluYa5 RNA that downregulates gene expression

Angiotensin-I-converting enzyme (ACE) is involved in the synthesis and degradation of important bioactive peptides. The ACE gene has a 287-bp insertion/deletion polymorphism that controls ACE expression through a mechanism that remains elusive. In this study, we found that the 287-bp polymorphic element of the ACE gene, a member of the AluYa5 sub-family of Alu elements, codes for an RNA molecule that controls the levels of ACE mRNA. Transient transfection of a plasmid containing a CMV promoter upstream of the ACE polymorphic element resulted in significant expression of an AluYa5 RNA and reduced ACE mRNA expression as well as ACE enzymatic activity in AD 293 cells. The AluYa5 element also independently reduced the expression of other genes, regardless of whether these genes harbored Alu elements within their genomic context. Interestingly, the CMV promoter was not required for the expression of the AluYa5 element in AD 293 cells. The 287-bp sequence was sufficient to produce AluYa5 RNA and led to a significant reduction in ACE gene expression. Moreover, the removal of an 11-bp fragment of the 3' end of the ACE polymorphic sequence, which is specific to this particular AluYa5 element, did not prevent this element from being expressed but did affect its ability to target ACE expression. Thus, the expression of the AluYa5 polymorphic element within the ACE gene could explain why patients carrying the ACE insertion polymorphism have reduced risk of developing several chronic diseases.

Natural antisense transcripts drive a regulatory cascade controlling c-MYC transcription

Cis-natural antisense transcripts (cis-NATs) are long noncoding RNAs transcribed from the opposite strand and overlapping coding and noncoding genes on the sense strand. cis-NATs are widely present in the human genome and can be involved in multiple mechanisms of gene regulation. Here, we describe the presence of cis-NATs in the 3′ distal region of the c-MYC locus and investigate their impact on transcriptional regulation of this key oncogene in human cancers. We found that cis-NATs are produced as consequence of the activation of cryptic transcription initiation sites in the 3′ distal region downstream of the c-MYC 3′UTR. The process is tightly regulated and leads to the formation of two main transcripts, NAT6531 and NAT6558, which differ in their ability to fold into stem-loop secondary structures. NAT6531 acts as a substrate for DICER and as a source of small RNAs capable of modulating c-MYC transcription. This complex system, based on the interplay between cis-NATs and NAT-derived small RNAs, may represent an important layer of epigenetic regulation of the expression of c-MYC and other genes in human cells.

Frequency of Alu insertions within the ACE and PR loci in Northwestern Mexicans

OBJECTIVE

Presently, non-LTR retrotransposons are the most active mobile elements in the human genome. Among these, Alu elements are highly represented in the modern population. Worldwide, distribution of Alu polymorphisms (insertion/deletion; I/D) shows variability between different populations. Two Alu insertion loci, ACE and PR, are significant biomarkers that have served in several genotype-phenotype association studies. In Mexico, studies concerning the frequency of these biomarkers have been conducted mainly in subpopulations from central and southern regions. Here, we screened a population sample of the northwestern region to gain further knowledge regarding the prevalence of Alu polymorphisms within ACE and PR loci.

RESULTS

For ACE locus, the observed genotype frequencies were 26.5, 51.0 and 22.5% for II, ID, and DD, respectively; and allelic frequencies for I and D were 52 and 48%. Whereas respective genotype frequencies for PR locus were 2.7, 26.5 and 70.8%, and the corresponding allele frequencies were 16 and 84%. Furthermore, the insertion frequency within ACE locus was similar between central, western and northwestern subpopulations, and rather higher in southeastern subpopulation (p < 0.05). Although the occurrence of Alu polymorphisms within PR locus has not been widely examined, the insertion frequency was higher in northwestern subpopulation, as compared with western and southeastern subpopulations (p < 0.05). Based on the frequency of Alu insertions found in ACE and PR loci, subpopulations from the northwestern, western and central regions share a common genetic origin, but apparently not with the subpopulation from the southeastern region, in accordance with the notion that assumes the existence of a broad genomic diversity in the Mexican population. In addition, the high prevalence of Alu insertions reveals their potential application as biomarkers with prognostic value for the associated diseases; e.g., as part of the standard protocols for clinical diagnosis.

Identification of RNA Polymerase III-Transcribed SINEs at Single-Locus Resolution from RNA Sequencing Data

Short Interspersed Element (SINE) retrotransposons are one of the most abundant DNA repeat elements in the human genome. They have been found to impact the expression of protein-coding genes, but the possible roles in cell physiology of their noncoding RNAs, generated by RNA polymerase (Pol) III, are just starting to be elucidated. For this reason, Short Interspersed Element (SINE) expression profiling is becoming mandatory to obtain a comprehensive picture of their regulatory roles. However, their repeated nature and frequent location within Pol II-transcribed genes represent a serious obstacle to the identification and quantification of genuine, Pol III-derived SINE transcripts at single-locus resolution on a genomic scale. Among the recent Next Generation Sequencing technologies, only RNA sequencing (RNA-Seq) holds the potential to solve these issues, even though both technical and biological matters need to be taken into account. A bioinformatic pipeline has been recently set up that, by exploiting RNA-seq features and knowledge of SINE transcription mechanisms, allows for easy identification and profiling of transcriptionally active genomic loci which are a source of genuine Pol III SINE transcripts.

Alu-miRNA interactions modulate transcript isoform diversity in stress response and reveal signatures of positive selection

Primate-specific Alus harbor different regulatory features, including miRNA targets. In this study, we provide evidence for miRNA-mediated modulation of transcript isoform levels during heat-shock response through exaptation of Alu-miRNA sites in mature mRNA. We performed genome-wide expression profiling coupled with functional validation of miRNA target sites within exonized Alus, and analyzed conservation of these targets across primates. We observed that two miRNAs (miR-15a-3p and miR-302d-3p) elevated in stress response, target RAD1, GTSE1, NR2C1, FKBP9 and UBE2I exclusively within Alu. These genes map onto the p53 regulatory network. Ectopic overexpression of miR-15a-3p downregulates GTSE1 and RAD1 at the protein level and enhances cell survival. This Alu-mediated fine-tuning seems to be unique to humans as evident from the absence of orthologous sites in other primate lineages. We further analyzed signatures of selection on Alu-miRNA targets in the genome, using 1000 Genomes Phase-I data. We found that 198 out of 3177 Alu-exonized genes exhibit signatures of selection within Alu-miRNA sites, with 60 of them containing SNPs supported by multiple evidences (global-F_ST > 0.3, pair-wise-F_ST > 0.5, Fay-Wu’s H < −20, iHS > 2.0, high ΔDAF) and implicated in p53 network. We propose that by affecting multiple genes, Alu-miRNA interactions have the potential to facilitate population-level adaptations in response to environmental challenges.

Alu Elements as Novel Regulators of Gene Expression in Type 1 Diabetes Susceptibility Genes?

Despite numerous studies implicating Alu repeat elements in various diseases, there is sparse information available with respect to the potential functional and biological roles of the repeat elements in Type 1 diabetes (T1D). Therefore, we performed a genome-wide sequence analysis of T1D candidate genes to identify embedded Alu elements within these genes. We observed significant enrichment of Alu elements within the T1D genes (p-value < 10e−16), which highlights their importance in T1D. Functional annotation of T1D genes harboring Alus revealed significant enrichment for immune-mediated processes (p-value < 10e−6). We also identified eight T1D genes harboring inverted Alus (IRAlus) within their 3' untranslated regions (UTRs) that are known to regulate the expression of host mRNAs by generating double stranded RNA duplexes. Our in silico analysis predicted the formation of duplex structures by IRAlus within the 3'UTRs of T1D genes. We propose that IRAlus might be involved in regulating the expression levels of the host T1D genes.

The RNA-centred view of the synapse: non-coding RNAs and synaptic plasticity

If mRNAs were the only RNAs made by a neuron, there would be a simple mapping of mRNAs to proteins. However, microRNAs and other non-coding RNAs (ncRNAs; endo-siRNAs, piRNAs, BC1, BC200, antisense and long ncRNAs, repeat-related transcripts, etc.) regulate mRNAs via effects on protein translation as well as transcriptional and epigenetic mechanisms. Not only are genes ON or OFF, but their ability to be translated can be turned ON or OFF at the level of synapses, supporting an enormous increase in information capacity. Here, I review evidence that ncRNAs are expressed pervasively within dendrites in mammalian brain; that some are activity-dependent and highly enriched near synapses; and that synaptic ncRNAs participate in plasticity responses including learning and memory. Ultimately, ncRNAs can be viewed as the post-it notes of the neuron. They have no literal meaning of their own, but derive their functions from where (and to what) they are stuck. This may explain, in part, why ncRNAs differ so dramatically from protein-coding genes, both in terms of the usual indicators of functionality and in terms of evolutionary constraints. ncRNAs do not appear to be direct mediators of synaptic transmission in the manner of neurotransmitters or receptors, yet they orchestrate synaptic plasticity—and may drive species-specific changes in cognition.

Dynamic Alu methylation during normal development, aging, and tumorigenesis

DNA methylation primarily occurs on CpG dinucleotides and plays an important role in transcriptional regulations during tissue development and cell differentiation. Over 25% of CpG dinucleotides in the human genome reside within Alu elements, the most abundant human repeats. The methylation of Alu elements is an important mechanism to suppress Alu transcription and subsequent retrotransposition. Decades of studies revealed that Alu methylation is highly dynamic during early development and aging. Recently, many environmental factors were shown to have a great impact on Alu methylation. In addition, aberrant Alu methylation has been documented to be an early event in many tumors and Alu methylation levels have been associated with tumor aggressiveness. The assessment of the Alu methylation has become an important approach for early diagnosis and/or prognosis of cancer. This review focuses on the dynamic Alu methylation during development, aging, and tumor genesis. The cause and consequence of Alu methylation changes will be discussed.

One hundred million adenosine-to-inosine RNA editing sites: hearing through the noise

The most recent work toward compiling a comprehensive database of adenosine-to-inosine RNA editing events suggests that the potential for RNA editing is much more pervasive than previously thought; indeed, it is manifest in more than 100 million potential editing events located primarily within Alu repeat elements of the human transcriptome. Pairs of inverted Alu repeats are found in a substantial number of human genes, and when transcribed, they form long double-stranded RNA structures that serve as optimal substrates for RNA editing enzymes. A small subset of edited Alu elements has been shown to exhibit diverse functional roles in the regulation of alternative splicing, miRNA repression, and cis-regulation of distant RNA editing sites. The low level of editing for the remaining majority may be non-functional, yet their persistence in the primate genome provides enhanced genomic flexibility that may be required for adaptive evolution.

The persistent contributions of RNA to eukaryotic gen(om)e architecture and cellular function

Currently, the best scenario for earliest forms of life is based on RNA molecules as they have the proven ability to catalyze enzymatic reactions and harbor genetic information. Evolutionary principles valid today become apparent in such models already. Furthermore, many features of eukaryotic genome architecture might have their origins in an RNA or RNA/protein (RNP) world, including the onset of a further transition, when DNA replaced RNA as the genetic bookkeeper of the cell. Chromosome maintenance, splicing, and regulatory function via RNA may be deeply rooted in the RNA/RNP worlds. Mostly in eukaryotes, conversion from RNA to DNA is still ongoing, which greatly impacts the plasticity of extant genomes. Raw material for novel genes encoding protein or RNA, or parts of genes including regulatory elements that selection can act on, continues to enter the evolutionary lottery.

Physiology of the read-write genome

Discoveries in cytogenetics, molecular biology, and genomics have revealed that genome change is an active cell-mediated physiological process. This is distinctly at variance with the pre-DNA assumption that genetic changes arise accidentally and sporadically. The discovery that DNA changes arise as the result of regulated cell biochemistry means that the genome is best modelled as a read-write (RW) data storage system rather than a read-only memory (ROM). The evidence behind this change in thinking and a consideration of some of its implications are the subjects of this article. Specific points include the following: cells protect themselves from accidental genome change with proofreading and DNA damage repair systems; localized point mutations result from the action of specialized trans-lesion mutator DNA polymerases; cells can join broken chromosomes and generate genome rearrangements by non-homologous end-joining (NHEJ) processes in specialized subnuclear repair centres; cells have a broad variety of natural genetic engineering (NGE) functions for transporting, diversifying and reorganizing DNA sequences in ways that generate many classes of genomic novelties; natural genetic engineering functions are regulated and subject to activation by a range of challenging life history events; cells can target the action of natural genetic engineering functions to particular genome locations by a range of well-established molecular interactions, including protein binding with regulatory factors and linkage to transcription; and genome changes in cancer can usefully be considered as consequences of the loss of homeostatic control over natural genetic engineering functions.

Vespucci: a system for building annotated databases of nascent transcripts

Global run-on sequencing (GRO-seq) is a recent addition to the series of high-throughput sequencing methods that enables new insights into transcriptional dynamics within a cell. However, GRO-sequencing presents new algorithmic challenges, as existing analysis platforms for ChIP-seq and RNA-seq do not address the unique problem of identifying transcriptional units de novo from short reads located all across the genome. Here, we present a novel algorithm for de novo transcript identification from GRO-sequencing data, along with a system that determines transcript regions, stores them in a relational database and associates them with known reference annotations. We use this method to analyze GRO-sequencing data from primary mouse macrophages and derive novel quantitative insights into the extent and characteristics of non-coding transcription in mammalian cells. In doing so, we demonstrate that Vespucci expands existing annotations for mRNAs and lincRNAs by defining the primary transcript beyond the polyadenylation site. In addition, Vespucci generates assemblies for un-annotated non-coding RNAs such as those transcribed from enhancer-like elements. Vespucci thereby provides a robust system for defining, storing and analyzing diverse classes of primary RNA transcripts that are of increasing biological interest.