The regulated retrotransposon transcriptome of mammalian cells

Measurement of telomere length using PNA probe by cytometry

Peptide nucleic acid (PNA) probes hybridize to denatured telomeric sequences in cells permeabilized in hot formamide. In reported protocols, the hybridization was conducted in solutions with high formamide concentrations to avoid the DNA renaturation that can hamper binding of the oligo-PNA probe to specific sequences. We postulated that telomeric DNA, confined in the nuclear microvolume, is not able to properly renature after hot formamide denaturation. Therefore, to improve hybridization conditions between the probe and the target sequences, it might be possible to add probe to sample after the complete removal of formamide.

Whole transcriptome analysis: what are we still missing?

New technologies such as tag-based sequencing and tiling arrays have provided unique insights into the transcriptional output of cells. Many new RNA classes have been uncovered in the past decade, despite limitations in current technologies. Even as the repertoire of known functional elements of the transcriptome increases and contemporary technologies become mainstream, inadequacies in conventional protocols for library preparation, sequencing and mapping continue to hamper revelation of the entire transcriptome of cells. In this article, we review current protocols and outline their deficiencies. We also provide our view on what we may be overlooking in the transcriptome, despite exhaustive investigations, and indicate future areas of technological development and research.

Interrogating the transcriptome of oocytes and preimplantation embryos

During its growth phase, a mouse oocyte accumulates RNA that is the sole template for new protein synthesis in the transcriptionally silent interval between growth completion and transcriptional activation of the embryonic genome. Over this transcriptionally silent interval, almost half the quantity of RNA accumulated in the full-grown oocyte is degraded, while stable messages undergo major transcript-specific polyadenylation fluctuations associated with timely translation of new proteins. These processes, in the background of substantial RNA degradation, create unique pitfalls for transcriptome analysis. Three particular challenges are discussed herein. (1) Systematic errors of relative quantification occur if standard approaches are used, wherein samples are normalized to a constant quantity of RNA, or when computational analyses are normalized to an apparent "constant" endogenous to the sample. We show that use of a fixed quantity of exogenous RNA per oocyte or embryo alleviates this problem. (2) Comparison of large-scale expression analyses from widely disparate platforms highlights how the differing protocols produce correspondingly different lists of genes with significant changes in transcript abundance. Only with careful attention to the differences among experiments can such discrepancies be understood. (3) The complete assessment of changes in expression requires correspondingly comprehensive assessment of the role of isoform-specific changes.

The reverse transcription inhibitor abacavir shows anticancer activity in prostate cancer cell lines

BACKGROUND

Transposable Elements (TEs) comprise nearly 45% of the entire genome and are part of sophisticated regulatory network systems that control developmental processes in normal and pathological conditions. The retroviral/retrotransposon gene machinery consists mainly of Long Interspersed Nuclear Elements (LINEs-1) and Human Endogenous Retroviruses (HERVs) that code for their own endogenous reverse transcriptase (RT). Interestingly, RT is typically expressed at high levels in cancer cells. Recent studies report that RT inhibition by non-nucleoside reverse transcriptase inhibitors (NNRTIs) induces growth arrest and cell differentiation in vitro and antagonizes growth of human tumors in animal model. In the present study we analyze the anticancer activity of Abacavir (ABC), a nucleoside reverse transcription inhibitor (NRTI), on PC3 and LNCaP prostate cancer cell lines.

PRINCIPAL FINDINGS

ABC significantly reduces cell growth, migration and invasion processes, considerably slows S phase progression, induces senescence and cell death in prostate cancer cells. Consistent with these observations, microarray analysis on PC3 cells shows that ABC induces specific and dose-dependent changes in gene expression, involving multiple cellular pathways. Notably, by quantitative Real-Time PCR we found that LINE-1 ORF1 and ORF2 mRNA levels were significantly up-regulated by ABC treatment.

CONCLUSIONS

Our results demonstrate the potential of ABC as anticancer agent able to induce antiproliferative activity and trigger senescence in prostate cancer cells. Noteworthy, we show that ABC elicits up-regulation of LINE-1 expression, suggesting the involvement of these elements in the observed cellular modifications.

Improved procedure for the measurement of telomere length in whole cells by PNA probe and flow cytometry

OBJECTIVES

Peptide nucleic acid (PNA) probes hybridize to denatured telomeric sequences in cells permeabilized in hot formamide. In reported protocols, the hybridization was conducted in solutions with high formamide concentrations to avoid the DNA renaturation that can hamper binding of the oligo-PNA probe to specific sequences. We postulated that telomeric DNA, confined in the nuclear microvolume, is not able to properly renature after hot formamide denaturation. Therefore, to improve hybridization conditions between the probe and the target sequences, it might be possible to add probe to sample after the complete removal of formamide.

MATERIALS AND METHODS

After telomeric DNA denaturation in hot formamide solution and several washes to remove the ionic solvent, cells were hybridized overnight at room temperature with human telomere-specific PNA probe conjugated with Cy5 fluorochrome, Cy5-OO-(CCCTAA)(3) . After stringency washes and staining with ethidium bromide, the cells were analysed by flow cytometry and by using a confocal microscope.

RESULTS

Using three continuous cell lines, different in DNA content and telomere length, and resting human peripheral blood T and B lymphocytes, we demonstrated that the oligo-PNA probe hybridized to telomeric sequences after complete removal of formamide and that in the preserved nucleus, telomeric sequence denaturation is irreversible.

CONCLUSION

According to our experience, oligo-PNA binding results is efficient, specific and proportional to telomere length. These, our original findings, can form the technological basis of actual in situ hybridization on preserved whole cells.

The central role of RNA in the genetic programming of complex organisms

Notwithstanding lineage-specific variations, the number and type of protein-coding genes remain relatively static across the animal kingdom. By contrast there has been a massive expansion in the extent of genomic non-proteincoding sequences with increasing developmental complexity. These non-coding sequences are, in fact, transcribed in a regulated manner to produce large numbers of large and small non-protein-coding RNAs that control gene expression at many levels including chromatin architecture, post-transcriptional processing and translation. Moreover, many RNAs are edited, especially in the nervous system, which may be the basis of epigenome-environment interactions and the function of the brain.

Repeat-element driven activation of proto-oncogenes in human malignancies

Recent data demonstrated that the aberrant activity of endogenous repetitive elements of the DNA in humans can drive the expression of proto-oncogenes. This article summarizes these results and gives an outlook on the impact of these findings on the pathogenesis and therapy of human cancer.

Update of the FANTOM web resource: from mammalian transcriptional landscape to its dynamic regulation

The international Functional Annotation Of the Mammalian Genomes 4 (FANTOM4) research collaboration set out to better understand the transcriptional network that regulates macrophage differentiation and to uncover novel components of the transcriptome employing a series of high-throughput experiments. The primary and unique technique is cap analysis of gene expression (CAGE), sequencing mRNA 5′-ends with a second-generation sequencer to quantify promoter activities even in the absence of gene annotation. Additional genome-wide experiments complement the setup including short RNA sequencing, microarray gene expression profiling on large-scale perturbation experiments and ChIP–chip for epigenetic marks and transcription factors. All the experiments are performed in a differentiation time course of the THP-1 human leukemic cell line. Furthermore, we performed a large-scale mammalian two-hybrid (M2H) assay between transcription factors and monitored their expression profile across human and mouse tissues with qRT-PCR to address combinatorial effects of regulation by transcription factors. These interdependent data have been analyzed individually and in combination with each other and are published in related but distinct papers. We provide all data together with systematic annotation in an integrated view as resource for the scientific community (http://fantom.gsc.riken.jp/4/). Additionally, we assembled a rich set of derived analysis results including published predicted and validated regulatory interactions. Here we introduce the resource and its update after the initial release.

Epigenetic control of retrotransposon expression in human embryonic stem cells

Long interspersed element 1s (LINE-1s or L1s) are a family of non-long-terminal-repeat retrotransposons that predominate in the human genome. Active LINE-1 elements encode proteins required for their mobilization. L1-encoded proteins also act in trans to mobilize short interspersed elements (SINEs), such as Alu elements. L1 and Alu insertions have been implicated in many human diseases, and their retrotransposition provides an ongoing source of human genetic diversity. L1/Alu elements are expected to ensure their transmission to subsequent generations by retrotransposing in germ cells or during early embryonic development. Here, we determined that several subfamilies of Alu elements are expressed in undifferentiated human embryonic stem cells (hESCs) and that most expressed Alu elements are active elements. We also exploited expression from the L1 antisense promoter to map expressed elements in hESCs. Remarkably, we found that expressed Alu elements are enriched in the youngest subfamily, Y, and that expressed L1s are mostly located within genes, suggesting an epigenetic control of retrotransposon expression in hESCs. Together, these data suggest that distinct subsets of active L1/Alu elements are expressed in hESCs and that the degree of somatic mosaicism attributable to L1 insertions during early development may be higher than previously anticipated.

Unique signatures of long noncoding RNA expression in response to virus infection and altered innate immune signaling

Studies of the host response to virus infection typically focus on protein-coding genes. However, non-protein-coding RNAs (ncRNAs) are transcribed in mammalian cells, and the roles of many of these ncRNAs remain enigmas. Using next-generation sequencing, we performed a whole-transcriptome analysis of the host response to severe acute respiratory syndrome coronavirus (SARS-CoV) infection across four founder mouse strains of the Collaborative Cross. We observed differential expression of approximately 500 annotated, long ncRNAs and 1,000 nonannotated genomic regions during infection. Moreover, studies of a subset of these ncRNAs and genomic regions showed the following. (i) Most were similarly regulated in response to influenza virus infection. (ii) They had distinctive kinetic expression profiles in type I interferon receptor and STAT1 knockout mice during SARS-CoV infection, including unique signatures of ncRNA expression associated with lethal infection. (iii) Over 40% were similarly regulated in vitro in response to both influenza virus infection and interferon treatment. These findings represent the first discovery of the widespread differential expression of long ncRNAs in response to virus infection and suggest that ncRNAs are involved in regulating the host response, including innate immunity. At the same time, virus infection models provide a unique platform for studying the biology and regulation of ncRNAs.

Most studies examining the host transcriptional response to infection focus only on protein-coding genes. However, there is growing evidence that thousands of non-protein-coding RNAs (ncRNAs) are transcribed from mammalian genomes. While most attention to the involvement of ncRNAs in virus-host interactions has been on small ncRNAs such as microRNAs, it is becoming apparent that many long ncRNAs (>200 nucleotides [nt]) are also biologically important. These long ncRNAs have been found to have widespread functionality, including chromatin modification and transcriptional regulation and serving as the precursors of small RNAs. With the advent of next-generation sequencing technologies, whole-transcriptome analysis of the host response, including long ncRNAs, is now possible. Using this approach, we demonstrated that virus infection alters the expression of numerous long ncRNAs, suggesting that these RNAs may be a new class of regulatory molecules that play a role in determining the outcome of infection.

Genome architecture marked by retrotransposons modulates predisposition to DNA methylation in cancer

Epigenetic silencing plays an important role in cancer development. An attractive hypothesis is that local DNA features may participate in differential predisposition to gene hypermethylation. We found that, compared with methylation-resistant genes, methylation-prone genes have a lower frequency of SINE and LINE retrotransposons near their transcription start site. In several large testing sets, this distribution was highly predictive of promoter methylation. Genome-wide analysis showed that 22% of human genes were predicted to be methylation-prone in cancer; these tended to be genes that are down-regulated in cancer and that function in developmental processes. Moreover, retrotransposon distribution marks a larger fraction of methylation-prone genes compared to Polycomb group protein (PcG) marking in embryonic stem cells; indeed, PcG marking and our predictive model based on retrotransposon frequency appear to be correlated but also complementary. In summary, our data indicate that retrotransposon elements, which are widespread in our genome, are strongly associated with gene promoter DNA methylation in cancer and may in fact play a role in influencing epigenetic regulation in normal and abnormal physiological states.

Orphan CpG islands identify numerous conserved promoters in the mammalian genome

CpG islands (CGIs) are vertebrate genomic landmarks that encompass the promoters of most genes and often lack DNA methylation. Querying their apparent importance, the number of CGIs is reported to vary widely in different species and many do not co-localise with annotated promoters. We set out to quantify the number of CGIs in mouse and human genomes using CXXC Affinity Purification plus deep sequencing (CAP-seq). We also asked whether CGIs not associated with annotated transcripts share properties with those at known promoters. We found that, contrary to previous estimates, CGI abundance in humans and mice is very similar and many are at conserved locations relative to genes. In each species CpG density correlates positively with the degree of H3K4 trimethylation, supporting the hypothesis that these two properties are mechanistically interdependent. Approximately half of mammalian CGIs (>10,000) are "orphans" that are not associated with annotated promoters. Many orphan CGIs show evidence of transcriptional initiation and dynamic expression during development. Unlike CGIs at known promoters, orphan CGIs are frequently subject to DNA methylation during development, and this is accompanied by loss of their active promoter features. In colorectal tumors, however, orphan CGIs are not preferentially methylated, suggesting that cancer does not recapitulate a developmental program. Human and mouse genomes have similar numbers of CGIs, over half of which are remote from known promoters. Orphan CGIs nevertheless have the characteristics of functional promoters, though they are much more likely than promoter CGIs to become methylated during development and hence lose these properties. The data indicate that orphan CGIs correspond to previously undetected promoters whose transcriptional activity may play a functional role during development.

Impact of nuclear organization and dynamics on epigenetic regulation in the central nervous system: implications for neurological disease states

Epigenetic mechanisms that are highly responsive to interoceptive and environmental stimuli mediate the proper execution of complex genomic programs, such as cell type-specific gene transcription and posttranscriptional RNA processing, and are increasingly thought to be important for modulating the development, homeostasis, and plasticity of the central nervous system (CNS). These epigenetic processes include DNA methylation, histone modifications, and chromatin remodeling, all of which play roles in neural cellular diversity, connectivity, and plasticity. Further, large-scale transcriptomic analyses have revealed that the eukaryotic genome is pervasively transcribed, forming interleaved protein-coding RNAs and regulatory nonprotein-coding RNAs (ncRNAs), which act through a broad array of molecular mechanisms. Most of these ncRNAs are transcribed in a cell type- and developmental stage-specific manner in the CNS. A broad array of posttranscriptional processes, such as RNA editing and transport, can modulate the functions of both protein-coding RNAs and ncRNAs. Additional studies implicate nuclear organization and dynamics in mediating epigenetic regulation. The compartmentalization of DNA sequences and other molecular machinery into functional nuclear domains, such as transcription factories, Cajal bodies, promyelocytic leukemia nuclear bodies, nuclear speckles, and paraspeckles, some of which are found prominently in neural cells, is associated with regulation of transcriptional activity and posttranscriptional RNA processing. These observations suggest that genomic architecture and RNA biology in the CNS are much more complex and nuanced than previously appreciated. Increasing evidence now suggests that most, if not all, human CNS diseases are associated with either primary or secondary perturbations in one or more aspects of the epigenome. In this review, we provide an update of our emerging understanding of genomic architecture, RNA biology, and nuclear organization and highlight the interconnected roles that deregulation of these factors may play in diverse CNS disorders.

The enemy within: an epigenetic role of retrotransposons in cancer initiation

This article proposes that cancers can be initiated by retrotransposon (RTN) activation through changes in the transcriptional regulation of nearby genes. I first detail the hypothesis and then discuss the nature of physiological stress(es) in RTN activation; the role of DNA demethylation in the initiation and propagation of new RTN states; the connection between ageing and cancer incidence and the involvement of activated RTNs in the chromosomal aberrations that feature in cancer progression. The hypothesis neither replaces nor invalidates other theories of cancer, in particular the somatic mutation theory, but helps clarify and unify much of the hitherto poorly integrated, complex phenomenology of cancer.

Transposable elements: insertion pattern and impact on gene expression evolution in hominids

Transposable elements (TEs) can affect the regulation of nearby genes through several mechanisms. Here, we examine to what extent recent TE insertions have contributed to the evolution of gene expression in hominids. We compare expression levels of human and chimpanzee orthologs and detect a weak increase in expression divergence (ED) for genes with species-specific TE insertions compared with unaffected genes. However, we show that genes with TE insertions predating the human-chimpanzee split also exhibit a similar increase in ED and therefore conclude that the increase is not due to the transcriptional influence of the TEs. These results are further confirmed by lineage-specific analysis of ED, using rhesus macaque as an outgroup: Human-chimpanzee ortholog pairs, where one ortholog has suffered TE insertion but not the other, do not show increased ED along the lineage where the insertion occurred, relative to the other lineage. We also show that genes with recent TE insertions tend to produce more alternative transcripts but find no evidence that the TEs themselves promote transcript diversity. Finally, we observe that TEs are enriched upstream relative to downstream of genes and show that this is due to insertional bias, rather than selection, because this bias is only observed in genes expressed in the germ line. This provides an alternative neutral explanation for the accumulation of TEs in upstream sequences.

SVA retrotransposons: Evolution and genetic instability

SINE-VNTR-Alus (SVA) are non-autonomous hominid specific retrotransposons that are associated with disease in humans. SVAs are evolutionarily young and presumably mobilized by the LINE-1 reverse transcriptase in trans. SVAs are currently active and may impact the host through a variety of mechanisms including insertional mutagenesis, exon shuffling, alternative splicing, and the generation of differentially methylated regions (DMR). Here we review SVA biology, including SVA insertions associated with known diseases. Further, we discuss a model describing the initial formation of SVA and the mechanisms by which SVA may impact the host.

LINE-1 retrotransposons: mediators of somatic variation in neuronal genomes?

LINE-1 (L1) elements are retrotransposons that insert extra copies of themselves throughout the genome using a 'copy and paste' mechanism. L1s comprise nearly approximately 20% of the human genome and are able to influence chromosome integrity and gene expression upon reinsertion. Recent studies show that L1 elements are active and 'jumping' during neuronal differentiation. New somatic L1 insertions could generate 'genomic plasticity' in neurons by causing variation in genomic DNA sequences and by altering the transcriptome of individual cells. Thus, L1-induced variation could affect neuronal plasticity and behavior. We discuss potential consequences of L1-induced neuronal diversity and propose that a mechanism for generating diversity in the brain could broaden the spectrum of behavioral phenotypes that can originate from any single genome.

Histone modifications at the blastocyst Axin1(Fu) locus mark the heritability of in vitro culture-induced epigenetic alterations in mice

For epigenetic phenotypes to be passed on from one generation to the next, it is required that epigenetic marks between generations are not cleared during the two stages of epigenetic reprogramming: mammalian gametogenesis and preimplantation development. The molecular nature of the chromatin marks involved in these events is unknown. Using the epigenetically inherited allele Axin1(Fu) (the result of a retrotransposon insertion upstream of the Axin1 gene) we sought to establish the heritable mark during early embryonic development that determines transgenerational epigenetic inheritance and to examine a possible shift in the expression of this epiallele in future progeny induced by in vitro culture (IVC). To identify the heritable mark we analyzed 1) the level of DNA methylation shown by the Axin1(Fu) allele in sperm and embryos at blastocysts stage and 2) the histone marks (H3K4 me2, H3K9 me3, H3K9 ac, and H4K20 me3) present at the blastocyst stage at the specific Axin1(Fu) locus. According to our data, histone H3K4 me2 and H3K9 ac mark the differences between the Axin1(Fu) penetrant and the silent locus during the first period of demethylation of the preimplantation development. Moreover, suboptimal IVC (reported to produce epigenetic alterations in embryos) and the histone deacetylase inhibitor trichostatin A affect the postnatal expression of this epigenetically sensitive allele through histone modifications during early development. This finding indicates that altered histone modifications during preimplantation can drive altered gene expression later on in development.

Tissue-specific transcript annotation and expression profiling with complementary next-generation sequencing technologies

Next-generation sequencing is excellently suited to evaluate the abundance of mRNAs to study gene expression. Here we compare two alternative technologies, cap analysis of gene expression (CAGE) and serial analysis of gene expression (SAGE), for the same RNA samples. Along with quantifying gene expression levels, CAGE can be used to identify tissue-specific transcription start sites, while SAGE monitors 3′-end usage. We used both methods to get more insight into the transcriptional control of myogenesis, studying differential gene expression in differentiated and proliferating C2C12 myoblast cells with statistical evaluation of reproducibility and differential gene expression. Both CAGE and SAGE provided highly reproducible data (Pearson's correlations >0.92 among biological triplicates). With both methods we found around 10 000 genes expressed at levels 2 transcripts per million (0.3 copies per cell), with an overlap of 86%. We identified 4304 and 3846 genes differentially expressed between proliferating and differentiated C2C12 cells by CAGE and SAGE, respectively, with an overlap of 2144. We identified 196 novel regulatory regions with preferential use in proliferating or differentiated cells. Next-generation sequencing of CAGE and SAGE libraries provides consistent expression levels and can enrich current genome annotations with tissue-specific promoters and alternative 3′-UTR usage.

Long-range function of an intergenic retrotransposon

Retrotransposons including endogenous retroviruses and their solitary long terminal repeats (LTRs) compose >40% of the human genome. Many of them are located in intergenic regions far from genes. Whether these intergenic retrotransposons serve beneficial host functions is not known. Here we show that an LTR retrotransposon of ERV-9 human endogenous retrovirus located 40-70 kb upstream of the human fetal gamma- and adult beta-globin genes serves a long-range, host function. The ERV-9 LTR contains multiple CCAAT and GATA motifs and competitively recruits a high concentration of NF-Y and GATA-2 present in low abundance in adult erythroid cells to assemble an LTR/RNA polymerase II complex. The LTR complex transcribes intergenic RNAs unidirectionally through the intervening DNA to loop with and modulate transcription factor occupancies at the far downstream globin promoters, thereby modulating globin gene switching by a competitive mechanism.

Linking promoters to functional transcripts in small samples with nanoCAGE and CAGEscan

Large-scale sequencing projects have revealed an unexpected complexity in the origins, structures and functions of mammalian transcripts. Many loci are known to produce overlapping coding and noncoding RNAs with capped 5' ends that vary in size. Methods to identify the 5' ends of transcripts will facilitate the discovery of new promoters and 5' ends derived from secondary capping events. Such methods often require high input amounts of RNA not obtainable from highly refined samples such as tissue microdissections and subcellular fractions. Therefore, we developed nano-cap analysis of gene expression (nanoCAGE), a method that captures the 5' ends of transcripts from as little as 10 ng of total RNA, and CAGEscan, a mate-pair adaptation of nanoCAGE that captures the transcript 5' ends linked to a downstream region. Both of these methods allow further annotation-agnostic studies of the complex human transcriptome.

Intergenic and repeat transcription in human, chimpanzee and macaque brains measured by RNA-Seq

Transcription is the first step connecting genetic information with an organism's phenotype. While expression of annotated genes in the human brain has been characterized extensively, our knowledge about the scope and the conservation of transcripts located outside of the known genes' boundaries is limited. Here, we use high-throughput transcriptome sequencing (RNA-Seq) to characterize the total non-ribosomal transcriptome of human, chimpanzee, and rhesus macaque brain. In all species, only 20-28% of non-ribosomal transcripts correspond to annotated exons and 20-23% to introns. By contrast, transcripts originating within intronic and intergenic repetitive sequences constitute 40-48% of the total brain transcriptome. Notably, some repeat families show elevated transcription. In non-repetitive intergenic regions, we identify and characterize 1,093 distinct regions highly expressed in the human brain. These regions are conserved at the RNA expression level across primates studied and at the DNA sequence level across mammals. A large proportion of these transcripts (20%) represents 3'UTR extensions of known genes and may play roles in alternative microRNA-directed regulation. Finally, we show that while transcriptome divergence between species increases with evolutionary time, intergenic transcripts show more expression differences among species and exons show less. Our results show that many yet uncharacterized evolutionary conserved transcripts exist in the human brain. Some of these transcripts may play roles in transcriptional regulation and contribute to evolution of human-specific phenotypic traits.

Long non-coding RNAs in nervous system function and disease

Central nervous system (CNS) development, homeostasis, stress responses, and plasticity are all mediated by epigenetic mechanisms that modulate gene expression and promote selective deployment of functional gene networks in response to complex profiles of interoceptive and environmental signals. Thus, not surprisingly, disruptions of these epigenetic processes are implicated in the pathogenesis of a spectrum of neurological and psychiatric diseases. Epigenetic mechanisms involve chromatin remodeling by relatively generic complexes that catalyze DNA methylation and various types of histone modifications. There is increasing evidence that these complexes are directed to their sites of action by long non-protein-coding RNAs (lncRNAs), of which there are tens if not hundreds of thousands specified in the genome. LncRNAs are transcribed in complex intergenic, overlapping and antisense patterns relative to adjacent protein-coding genes, suggesting that many lncRNAs regulate the expression of these genes. LncRNAs also participate in a wide array of subcellular processes, including the formation and function of cellular organelles. Most lncRNAs are transcribed in a developmentally regulated and cell type specific manner, particularly in the CNS, wherein over half of all lncRNAs are expressed. While the numerous biological functions of lncRNAs are yet to be characterized fully, a number of recent studies suggest that lnRNAs are important for mediating cell identity. This function seems to be especially important for generating the enormous array of regional neuronal and glial cell subtypes that are present in the CNS. Further studies have also begun to elucidate additional roles played by lncRNAs in CNS processes, including homeostasis, stress responses and plasticity. Herein, we review emerging evidence that highlights the expression and function of lncRNAs in the CNS and suggests that lncRNA deregulation is an important factor in various CNS pathologies including neurodevelopmental, neurodegenerative and neuroimmunological disorders, primary brain tumors, and psychiatric diseases.

LINE-1 activity in facultative heterochromatin formation during X chromosome inactivation

During X chromosome inactivation (XCI), Xist RNA coats and silences one of the two X chromosomes in female cells. Little is known about how XCI spreads across the chromosome, although LINE-1 elements have been proposed to play a role. Here we show that LINEs participate in creating a silent nuclear compartment into which genes become recruited. A subset of young LINE-1 elements, however, is expressed during XCI, rather than being silenced. We demonstrate that such LINE expression requires the specific heterochromatic state induced by Xist. These LINEs often lie within escape-prone regions of the X chromosome, but close to genes that are subject to XCI, and are associated with putative endo-siRNAs. LINEs may thus facilitate XCI at different levels, with silent LINEs participating in assembly of a heterochromatic nuclear compartment induced by Xist, and active LINEs participating in local propagation of XCI into regions that would otherwise be prone to escape.

The genome sequence of the spontaneously hypertensive rat: Analysis and functional significance

The spontaneously hypertensive rat (SHR) is the most widely studied animal model of hypertension. Scores of SHR quantitative loci (QTLs) have been mapped for hypertension and other phenotypes. We have sequenced the SHR/OlaIpcv genome at 10.7-fold coverage by paired-end sequencing on the Illumina platform. We identified 3.6 million high-quality single nucleotide polymorphisms (SNPs) between the SHR/OlaIpcv and Brown Norway (BN) reference genome, with a high rate of validation (sensitivity 96.3%-98.0% and specificity 99%-100%). We also identified 343,243 short indels between the SHR/OlaIpcv and reference genomes. These SNPs and indels resulted in 161 gain or loss of stop codons and 629 frameshifts compared with the BN reference sequence. We also identified 13,438 larger deletions that result in complete or partial absence of 107 genes in the SHR/OlaIpcv genome compared with the BN reference and 588 copy number variants (CNVs) that overlap with the gene regions of 688 genes. Genomic regions containing genes whose expression had been previously mapped as cis-regulated expression quantitative trait loci (eQTLs) were significantly enriched with SNPs, short indels, and larger deletions, suggesting that some of these variants have functional effects on gene expression. Genes that were affected by major alterations in their coding sequence were highly enriched for genes related to ion transport, transport, and plasma membrane localization, providing insights into the likely molecular and cellular basis of hypertension and other phenotypes specific to the SHR strain. This near complete catalog of genomic differences between two extensively studied rat strains provides the starting point for complete elucidation, at the molecular level, of the physiological and pathophysiological phenotypic differences between individuals from these strains.

Nuclear stress bodies

Nuclear stress bodies (nSBs) are unique subnuclear organelles which form in response to heat shock. They are initiated through a direct interaction between heat shock transcription factor 1 (HSF1) and pericentric tandem repeats of satellite III sequences and correspond to active transcription sites for noncoding satellite III transcripts. Given their unusual features, nSBs are distinct from other known transcription sites. In stressed cells, they are thought to participate in rapid, transient, and global reprogramming of gene expression through different types of mechanisms including chromatin remodeling and trapping of transcription and splicing factors. The analysis of these atypical and intriguing structures uncovers new facets of the relationship between nuclear organization and nuclear function.

Hydroxyurea-induced global transcriptional suppression in mouse ES cells

Hydroxyurea (HU), as a therapeutic medicine, has been extensively used clinically. To further survey molecular mechanisms of HU treatment, we analyzed global transcriptomic alteration of mouse ES cells in response to the treatment using high-throughput sequencing. We show that the global transcriptional activity is significantly suppressed as cells are exposed to HU treatment and alters multiple key cellular pathways, including cell cycle, apoptosis and DNAs. HU treatment also alters alternative splicing mechanisms and suppresses non-coding RNA expression. Our result provides novel clues for the understanding of how cells respond to HU and further suggests that high-throughput sequencing technology provides a powerful tool to study mechanisms of clinical drugs at the cellular level.

Epigenomic diversity of colorectal cancer indicated by LINE-1 methylation in a database of 869 tumors

BACKGROUND

Genome-wide DNA hypomethylation plays a role in genomic instability and carcinogenesis. LINE-1 (L1 retrotransposon) constitutes a substantial portion of the human genome, and LINE-1 methylation correlates with global DNA methylation status. LINE-1 hypomethylation in colon cancer has been strongly associated with poor prognosis. However, whether LINE-1 hypomethylators constitute a distinct cancer subtype remains uncertain. Recent evidence for concordant LINE-1 hypomethylation within synchronous colorectal cancer pairs suggests the presence of a non-stochastic mechanism influencing tumor LINE-1 methylation level. Thus, it is of particular interest to examine whether its wide variation can be attributed to clinical, pathologic or molecular features.

DESIGN

Utilizing a database of 869 colorectal cancers in two prospective cohort studies, we constructed multivariate linear and logistic regression models for LINE-1 methylation (quantified by Pyrosequencing). Variables included age, sex, body mass index, family history of colorectal cancer, smoking status, tumor location, stage, grade, mucinous component, signet ring cells, tumor infiltrating lymphocytes, CpG island methylator phenotype (CIMP), microsatellite instability, expression of TP53 (p53), CDKN1A (p21), CTNNB1 (beta-catenin), PTGS2 (cyclooxygenase-2), and FASN, and mutations in KRAS, BRAF, and PIK3CA.

RESULTS

Tumoral LINE-1 methylation ranged from 23.1 to 90.3 of 0-100 scale (mean 61.4; median 62.3; standard deviation 9.6), and distributed approximately normally except for extreme hypomethylators [LINE-1 methylation < 40; N = 22 (2.5%), which were far more than what could be expected by normal distribution]. LINE-1 extreme hypomethylators were significantly associated with younger patients (p = 0.0058). Residual plot by multivariate linear regression showed that LINE-1 extreme hypomethylators clustered as one distinct group, separate from the main tumor group. The multivariate linear regression model could explain 8.4% of the total variability of LINE-1 methylation (R-square = 0.084). Multivariate logistic regression models for binary LINE-1 hypomethylation outcomes (cutoffs of 40, 50 and 60) showed at most fair predictive ability (area under receiver operator characteristics curve < 0.63).

CONCLUSIONS

LINE-1 extreme hypomethylators appear to constitute a previously-unrecognized, distinct subtype of colorectal cancers, which needs to be confirmed by additional studies. Our tumor LINE-1 methylation data indicate enormous epigenomic diversity of individual colorectal cancers.

Repressive and active histone methylation mark distinct promoters in human and mouse spermatozoa

In higher eukaryotes, histone methylation is involved in maintaining cellular identity during somatic development. As most nucleosomes are replaced by protamines during spermatogenesis, it is unclear whether histone modifications function in paternal transmission of epigenetic information. Here we show that two modifications important for Trithorax- and Polycomb-mediated gene regulation have methylation-specific distributions at regulatory regions in human spermatozoa. Histone H3 Lys4 dimethylation (H3K4me2) marks genes that are relevant in spermatogenesis and cellular homeostasis. In contrast, histone H3 Lys27 trimethylation (H3K27me3) marks developmental regulators in sperm, as in somatic cells. However, nucleosomes are only moderately retained at regulatory regions in human sperm. Nonetheless, genes with extensive H3K27me3 coverage around transcriptional start sites in particular tend not to be expressed during male and female gametogenesis or in preimplantation embryos. Promoters of orthologous genes are similarly modified in mouse spermatozoa. These data are compatible with a role for Polycomb in repressing somatic determinants across generations, potentially in a variegating manner.

Derepression of an endogenous long terminal repeat activates the CSF1R proto-oncogene in human lymphoma

Mammalian genomes contain many repetitive elements, including long terminal repeats (LTRs), which have long been suspected to have a role in tumorigenesis. Here we present evidence that aberrant LTR activation contributes to lineage-inappropriate gene expression in transformed human cells and that such gene expression is central for tumor cell survival. We show that B cell-derived Hodgkin's lymphoma cells depend on the activity of the non-B, myeloid-specific proto-oncogene colony-stimulating factor 1 receptor (CSF1R). In these cells, CSF1R transcription initiates at an aberrantly activated endogenous LTR of the MaLR family (THE1B). Derepression of the THE1 subfamily of MaLR LTRs is widespread in the genome of Hodgkin's lymphoma cells and is associated with impaired epigenetic control due to loss of expression of the corepressor CBFA2T3. Furthermore, we detect LTR-driven CSF1R transcripts in anaplastic large cell lymphoma, in which CSF1R is known to be expressed aberrantly. We conclude that LTR derepression is involved in the pathogenesis of human lymphomas, a finding that might have diagnostic, prognostic and therapeutic implications.

Hypomethylation of a LINE-1 promoter activates an alternate transcript of the MET oncogene in bladders with cancer

It was recently shown that a large portion of the human transcriptome can originate from within repetitive elements, leading to ectopic expression of protein-coding genes. However the mechanism of transcriptional activation of repetitive elements has not been definitively elucidated. For the first time, we directly demonstrate that hypomethylation of retrotransposons can cause altered gene expression in humans. We also reveal that active LINE-1s switch from a tetranucleosome to dinucleosome structure, acquiring H2A.Z- and nucleosome-free regions upstream of TSSs, previously shown only at active single-copy genes. Hypomethylation of a specific LINE-1 promoter was also found to induce an alternate transcript of the MET oncogene in bladder tumors and across the entire urothelium of tumor-bearing bladders. These data show that, in addition to contributing to chromosomal instability, hypomethylation of LINE-1s can alter the functional transcriptome and plays a role not only in human disease but also in disease predisposition.

A surprisingly large portion of our transcriptome originates within repetitive elements, most commonly LINE-1s. However, the mechanism of activation has not been definitively shown. We directly demonstrate for the first time the causal relationship between DNA hypomethylation and transcriptional activation of LINE-1 promoters. Hypomethylation of specific LINE-1 promoters can alter the transcriptome, including activating an alternate transcript of the MET oncogene, not only in primary bladder tumors but also in premalignant urothelium across entire bladders with tumors. Our study has important implications for tumor biology, cancer detection, and treatment, and it also answers the long-standing question of whether hypomethylation of retrotransposons induces ectopic gene expression and influences disease susceptibility in humans, a phenomenon first described in agouti mice.

Core promoter structure and genomic context reflect histone 3 lysine 9 acetylation patterns

BACKGROUND

Histone modifications play an important role in gene regulation. Acetylation of histone 3 lysine 9 (H3K9ac) is generally associated with transcription initiation and unfolded chromatin, thereby positively influencing gene expression. Deep sequencing of the 5' ends of gene transcripts using DeepCAGE delivers detailed information about the architecture and expression level of gene promoters. The combination of H3K9ac ChIP-chip and DeepCAGE in a myeloid leukemia cell line (THP-1) allowed us to study the spatial distribution of H3K9ac around promoters using a novel clustering approach. The promoter classes were analyzed for association with relevant genomic sequence features.

RESULTS

We performed a clustering of 4,481 promoters according to their surrounding H3K9ac signal and analyzed the clustered promoters for association with different sequence features. The clustering revealed three groups with major H3K9ac signal upstream, centered and downstream of the promoter. Narrow single peak promoters tend to have a concentrated activity of H3K9ac in the upstream region, while broad promoters tend to have a concentrated activity of H3K9ac and RNA polymerase II binding in the centered and downstream regions. A subset of promoters with high gene expression level, compared to subsets with low and medium gene expression, shows dramatic increase in H3K9ac activity in the upstream cluster only; this may indicate that promoters in the centered and downstream clusters are predominantly regulated at post-initiation steps. Furthermore, the upstream cluster is depleted in CpG islands and more likely to regulate un-annotated genes.

CONCLUSIONS

Clustering core promoters according to their surrounding acetylation signal is a promising approach for the study of histone modifications. When examining promoters clustered into groups according to their surrounding H3K9 acetylation signal, we find that the relative localization and intensity of H3K9ac is very specific depending on characteristic sequence features of the promoter. Experimental data from DeepCAGE and ChIP-chip experiments using undifferentiated (monocyte) and differentiated (macrophage) THP-1 cells leads us to the same conclusions.

Transcriptome sequencing in an ecologically important tree species: assembly, annotation, and marker discovery

Background

Massively parallel sequencing of cDNA is now an efficient route for generating enormous sequence collections that represent expressed genes. This approach provides a valuable starting point for characterizing functional genetic variation in non-model organisms, especially where whole genome sequencing efforts are currently cost and time prohibitive. The large and complex genomes of pines (Pinus spp.) have hindered the development of genomic resources, despite the ecological and economical importance of the group. While most genomic studies have focused on a single species (P. taeda), genomic level resources for other pines are insufficiently developed to facilitate ecological genomic research. Lodgepole pine (P. contorta) is an ecologically important foundation species of montane forest ecosystems and exhibits substantial adaptive variation across its range in western North America. Here we describe a sequencing study of expressed genes from P. contorta, including their assembly and annotation, and their potential for molecular marker development to support population and association genetic studies.

Results

We obtained 586,732 sequencing reads from a 454 GS XLR70 Titanium pyrosequencer (mean length: 306 base pairs). A combination of reference-based and de novo assemblies yielded 63,657 contigs, with 239,793 reads remaining as singletons. Based on sequence similarity with known proteins, these sequences represent approximately 17,000 unique genes, many of which are well covered by contig sequences. This sequence collection also included a surprisingly large number of retrotransposon sequences, suggesting that they are highly transcriptionally active in the tissues we sampled. We located and characterized thousands of simple sequence repeats and single nucleotide polymorphisms as potential molecular markers in our assembled and annotated sequences. High quality PCR primers were designed for a substantial number of the SSR loci, and a large number of these were amplified successfully in initial screening.

Conclusions

This sequence collection represents a major genomic resource for P. contorta, and the large number of genetic markers characterized should contribute to future research in this and other pines. Our results illustrate the utility of next generation sequencing as a basis for marker development and population genomics in non-model species.

Transcriptome sequencing in an ecologically important tree species: assembly, annotation, and marker discovery

BACKGROUND

Massively parallel sequencing of cDNA is now an efficient route for generating enormous sequence collections that represent expressed genes. This approach provides a valuable starting point for characterizing functional genetic variation in non-model organisms, especially where whole genome sequencing efforts are currently cost and time prohibitive. The large and complex genomes of pines (Pinus spp.) have hindered the development of genomic resources, despite the ecological and economical importance of the group. While most genomic studies have focused on a single species (P. taeda), genomic level resources for other pines are insufficiently developed to facilitate ecological genomic research. Lodgepole pine (P. contorta) is an ecologically important foundation species of montane forest ecosystems and exhibits substantial adaptive variation across its range in western North America. Here we describe a sequencing study of expressed genes from P. contorta, including their assembly and annotation, and their potential for molecular marker development to support population and association genetic studies.

RESULTS

We obtained 586,732 sequencing reads from a 454 GS XLR70 Titanium pyrosequencer (mean length: 306 base pairs). A combination of reference-based and de novo assemblies yielded 63,657 contigs, with 239,793 reads remaining as singletons. Based on sequence similarity with known proteins, these sequences represent approximately 17,000 unique genes, many of which are well covered by contig sequences. This sequence collection also included a surprisingly large number of retrotransposon sequences, suggesting that they are highly transcriptionally active in the tissues we sampled. We located and characterized thousands of simple sequence repeats and single nucleotide polymorphisms as potential molecular markers in our assembled and annotated sequences. High quality PCR primers were designed for a substantial number of the SSR loci, and a large number of these were amplified successfully in initial screening.

CONCLUSIONS

This sequence collection represents a major genomic resource for P. contorta, and the large number of genetic markers characterized should contribute to future research in this and other pines. Our results illustrate the utility of next generation sequencing as a basis for marker development and population genomics in non-model species.

Large-scale transcriptome data reveals transcriptional activity of fission yeast LTR retrotransposons

Background

Retrotransposons are transposable elements that proliferate within eukaryotic genomes through a process involving reverse transcription. The numbers of retrotransposons within genomes and differences between closely related species may yield insight into the evolutionary history of the elements. Less is known about the ongoing dynamics of retrotransposons, as analysis of genome sequences will only reveal insertions of retrotransposons that are fixed - or near fixation - in the population or strain from which genetic material has been extracted for sequencing. One pre-requisite for retrotransposition is transcription of the elements. Given their intrinsic sequence redundancy, transcriptome-level analyses of transposable elements are scarce. We have used recently published transcriptome data from the fission yeast Schizosaccharomyces pombe to assess the ability to detect and describe transcriptional activity from Long Terminal Repeat (LTR) retrotransposons. LTR retrotransposons are normally flanked by two LTR sequences. However, the majority of LTR sequences in S. pombe exist as solitary LTRs, i.e. as single terminal repeat sequences not flanking a retrotransposon. Transcriptional activity was analysed for both full-length LTR retrotransposons and solitary LTRs.

Results

Two independent sets of transcriptome data reveal the presence of full-length, polyadenylated transcripts from LTR retrotransposons in S. pombe during growth phase in rich medium. The redundancy of retrotransposon sequences makes it difficult to assess which elements are transcriptionally active, but data strongly indicates that only a subset of the LTR retrotransposons contribute significantly to the detected transcription. A considerable level of reverse strand transcription is also detected. Equal levels of transcriptional activity are observed from both strands of solitary LTR sequences. Transcriptome data collected during meiosis suggests that transcription of solitary LTRs is correlated with the transcription of nearby protein-coding genes.

Conclusions

Presumably, the host organism negatively regulates proliferation of LTR retrotransposons. The finding of considerable transcriptional activity of retrotransposons suggests that part of this regulation is likely to take place at a post-transcriptional level. Alternatively, the transcriptional activity may signify a hitherto unrecognized activity level of retrotransposon proliferation. Our findings underline the usefulness of transcriptome data in elucidating dynamics in retrotransposon transcription.

RNA dust: where are the genes?

Initial gene discovery efforts through analysis of genome sequences and identification and characterization of expressed RNAs have revealed that only a relatively small portion of the genome is transcribed into protein coding mRNAs in vertebrates. However, in contrast with this paucity of protein coding ‘genes’, there is an enormous complexity in transcription and the protein coding mRNAs contribute to a very small fraction of transcripts in comparison with the different varieties of non-coding RNAs (ncRNAs). This transcriptome complexity may be hypothesized to have a regulatory role that is required for the development and function of organisms as complex as vertebrates. At the same time, it raises the fundamental question of the unequivocal definition of a gene. It is intriguing to postulate that many ncRNAs might finely modulate gene activity by acting as regulatory elements. The emerging hypotheses suggest that the gene regulatory machinery may be deeply interconnected with the world of short RNAs. These RNAs may generally act for fine-tuning of the protein-coding transcriptome.

KAP1 controls endogenous retroviruses in embryonic stem cells

More than forty per cent of the mammalian genome is derived from retroelements, of which about one-quarter are endogenous retroviruses (ERVs). Some are still active, notably in mice the highly polymorphic early transposon (ETn)/MusD and intracisternal A-type particles (IAP). ERVs are transcriptionally silenced during early embryogenesis by histone and DNA methylation (and reviewed in ref. 7), although the initiators of this process, which is essential to protect genome integrity, remain largely unknown. KAP1 (KRAB-associated protein 1, also known as tripartite motif-containing protein 28, TRIM28) represses genes by recruiting the histone methyltransferase SETDB1, heterochromatin protein 1 (HP1) and the NuRD histone deacetylase complex, but few of its physiological targets are known. Two lines of evidence suggest that KAP1-mediated repression could contribute to the control of ERVs: first, KAP1 can trigger permanent gene silencing during early embryogenesis, and second, a KAP1 complex silences the retrovirus murine leukaemia virus in embryonic cells. Consistent with this hypothesis, here we show that KAP1 deletion leads to a marked upregulation of a range of ERVs, in particular IAP elements, in mouse embryonic stem (ES) cells and in early embryos. We further demonstrate that KAP1 acts synergistically with DNA methylation to silence IAP elements, and that it is enriched at the 5' untranslated region (5'UTR) of IAP genomes, where KAP1 deletion leads to the loss of histone 3 lysine 9 trimethylation (H3K9me3), a hallmark of KAP1-mediated repression. Correspondingly, IAP 5'UTR sequences can impose in cis KAP1-dependent repression on a heterologous promoter in ES cells. Our results establish that KAP1 controls endogenous retroelements during early embryonic development.

Distinctive patterns of age-dependent hypomethylation in interspersed repetitive sequences

Interspersed repetitive sequences (IRSs) are a major contributor to genome size and may contribute to cellular functions. IRSs are subdivided according to size and functionally related structures into short interspersed elements, long interspersed elements (LINEs), DNA transposons, and LTR-retrotransposons. Many IRSs may produce RNA and regulate genes by a variety of mechanisms. The majority of DNA methylation occurs in IRSs and is believed to suppress IRS activities. Global hypomethylation, or the loss of genome-wide methylation, is a common epigenetic event not only in senescent cells but also in cancer cells. Loss of LINE-1 methylation has been characterized in many cancers. Here, we evaluated the methylation levels of peripheral blood mononuclear cells of LINE-1, Alu, and human endogenous retrovirus K (HERV-K) in 177 samples obtained from volunteers between 20 and 88 yr of age. Age was negatively associated with methylation levels of Alu (r = -0.452, P < 10(-3)) and HERV-K (r = -0.326, P < 10(-3)) but not LINE-1 (r = 0.145, P = 0.055). Loss of methylation of Alu occurred during ages 34-68 yr, and loss of methylation of HERV-K occurred during ages 40-63 yr and again during ages 64-83 yr. Interestingly, methylation of Alu and LINE-1 are directly associated, particularly at ages 49 yr and older (r = 0.49, P < 10(-3)). Therefore, only some types of IRSs lose methylation at certain ages. Moreover, Alu and HERV-K become hypomethylated differently. Finally, there may be several mechanisms of global methylation. However, not all of these mechanisms are age-dependent. This finding may lead to a better understanding of not only the biological causes and consequences of genome-wide hypomethylation but also the role of IRSs in the aging process.

Highly conserved motifs in non-coding regions of Sirevirus retrotransposons: the key for their pattern of distribution within and across plants?

BACKGROUND

Retrotransposons are key players in the evolution of eukaryotic genomes. Moreover, it is now known that some retrotransposon classes, like the abundant and plant-specific Sireviruses, have intriguingly distinctive host preferences. Yet, it is largely unknown if this bias is supported by different genome structures.

RESULTS

We performed sensitive comparative analysis of the genomes of a large set of Ty1/copia retrotransposons. We discovered that Sireviruses are unique among Pseudoviridae in that they constitute an ancient genus characterized by vastly divergent members, which however contain highly conserved motifs in key non-coding regions: multiple polypurine tract (PPT) copies cluster upstream of the 3' long terminal repeat (3'LTR), of which the terminal PPT tethers to a distinctive attachment site and is flanked by a precisely positioned inverted repeat. Their LTRs possess a novel type of repeated motif (RM) defined by its exceptionally high copy number, symmetry and core CGG-CCG signature. These RM boxes form CpG islands and lie a short distance upstream of a conserved promoter region thus hinting towards regulatory functions. Intriguingly, in the envelope-containing Sireviruses additional boxes cluster at the 5' vicinity of the envelope. The 5'LTR/internal domain junction and a polyC-rich integrase signal are also highly conserved domains of the Sirevirus genome.

CONCLUSIONS

Our comparative analysis of retrotransposon genomes using advanced in silico methods highlighted the unique genome organization of Sireviruses. Their structure may dictate a life cycle with different regulation and transmission strategy compared to other Pseudoviridae, which may contribute towards their pattern of distribution within and across plants.

Functions and regulation of RNA editing by ADAR deaminases

One type of RNA editing converts adenosines to inosines (A-->I editing) in double-stranded RNA (dsRNA) substrates. A-->I RNA editing is mediated by adenosine deaminase acting on RNA (ADAR) enzymes. A-->I RNA editing of protein-coding sequences of a limited number of mammalian genes results in recoding and subsequent alterations of their functions. However, A-->I RNA editing most frequently targets repetitive RNA sequences located within introns and 5' and 3' untranslated regions (UTRs). Although the biological significance of noncoding RNA editing remains largely unknown, several possibilities, including its role in the control of endogenous short interfering RNAs (esiRNAs), have been proposed. Furthermore, recent studies have revealed that the biogenesis and functions of certain microRNAs (miRNAs) are regulated by the editing of their precursors. Here, I review the recent findings that indicate new functions for A-->I editing in the regulation of noncoding RNAs and for interactions between RNA editing and RNA interference mechanisms.

deepBase: a database for deeply annotating and mining deep sequencing data

Advances in high-throughput next-generation sequencing technology have reshaped the transcriptomic research landscape. However, exploration of these massive data remains a daunting challenge. In this study, we describe a novel database, deepBase, which we have developed to facilitate the comprehensive annotation and discovery of small RNAs from transcriptomic data. The current release of deepBase contains deep sequencing data from 185 small RNA libraries from diverse tissues and cell lines of seven organisms: human, mouse, chicken, Ciona intestinalis, Drosophila melanogaster, Caenhorhabditis elegans and Arabidopsis thaliana. By analyzing ∼14.6 million unique reads that perfectly mapped to more than 284 million genomic loci, we annotated and identified ∼380 000 unique ncRNA-associated small RNAs (nasRNAs), ∼1.5 million unique promoter-associated small RNAs (pasRNAs), ∼4.0 million unique exon-associated small RNAs (easRNAs) and ∼6 million unique repeat-associated small RNAs (rasRNAs). Furthermore, 2038 miRNA and 1889 snoRNA candidates were predicted by miRDeep and snoSeeker. All of the mapped reads can be grouped into about 1.2 million RNA clusters. For the purpose of comparative analysis, deepBase provides an integrative, interactive and versatile display. A convenient search option, related publications and other useful information are also provided for further investigation. deepBase is available at: http://deepbase.sysu.edu.cn/.

Discovery of short pseudogenes derived from messenger RNAs

More than 40% of the human genome is generated by retrotransposition, a series of in vivo processes involving reverse transcription of RNA molecules and integration of the transcripts into the genomic sequence. The mechanism of retrotransposition, however, is not fully understood, and additional genomic elements generated by retrotransposition may remain to be discovered. Here, we report that the human genome contains many previously unidentified short pseudogenes generated by retrotransposition of mRNAs. Genomic elements generated by non-long terminal repeat retrotransposition have specific sequence signatures: a poly-A tract that is immediately downstream and a pair of duplicated sequences, called target site duplications (TSDs), at either end. Using a new computer program, TSDscan, that can accurately detect pseudogenes based on the presence of the poly-A tract and TSDs, we found 654 short (≤300 bp), previously unknown pseudogenes derived from mRNAs. Comprehensive analyses of the pseudogenes that we identified and their parent mRNAs revealed that the pseudogene length depends on the parent mRNA length: long mRNAs generate more short pseudogenes than do short mRNAs. To explain this phenomenon, we hypothesize that most long mRNAs are truncated before they are reverse transcribed. Truncated mRNAs would be rapidly degraded during reverse transcription, resulting in the generation of short pseudogenes.

Multiple evolutionary rate classes in animal genome evolution

The proportion of functional sequence in the human genome is currently a subject of debate. The most widely accepted figure is that approximately 5% is under purifying selection. In Drosophila, estimates are an order of magnitude higher, though this corresponds to a similar quantity of sequence. These estimates depend on the difference between the distribution of genomewide evolutionary rates and that observed in a subset of sequences presumed to be neutrally evolving. Motivated by the widening gap between these estimates and experimental evidence of genome function, especially in mammals, we developed a sensitive technique for evaluating such distributions and found that they are much more complex than previously apparent. We found strong evidence for at least nine well-resolved evolutionary rate classes in an alignment of four Drosophila species and at least seven classes in an alignment of four mammals, including human. We also identified at least three rate classes in human ancestral repeats. By positing that the largest of these ancestral repeat classes is neutrally evolving, we estimate that the proportion of nonneutrally evolving sequence is 30% of human ancestral repeats and 45% of the aligned portion of the genome. However, we also question whether any of the classes represent neutrally evolving sequences and argue that a plausible alternative is that they reflect variable structure-function constraints operating throughout the genomes of complex organisms.

A global view of genomic information--moving beyond the gene and the master regulator

The current view of gene regulation in complex organisms holds that gene expression is largely controlled by the combinatoric actions of transcription factors and other regulatory proteins, some of which powerfully influence cell type. Recent large-scale studies have confirmed that cellular differentiation involves many different regulatory factors. However, other studies indicate that the genome is pervasively transcribed to produce a variety of short and long non-protein-coding RNAs, including those derived from retrotransposed sequences, which also play important roles in the epigenetic regulation of gene expression. The evidence suggests that ontogenesis requires interplay between state-specific regulatory proteins, multitasked effector complexes and target-specific RNAs that recruit these complexes to their sites of action. Moreover, the semi-continuous nature of the transcriptome prompts the reassessment of 'genes' as discrete entities and indicates that the mammalian genome might be more accurately viewed as islands of protein-coding information in a sea of cis- and trans-acting regulatory sequences.

The impact of retrotransposons on human genome evolution

Their ability to move within genomes gives transposable elements an intrinsic propensity to affect genome evolution. Non-long terminal repeat (LTR) retrotransposons--including LINE-1, Alu and SVA elements--have proliferated over the past 80 million years of primate evolution and now account for approximately one-third of the human genome. In this Review, we focus on this major class of elements and discuss the many ways that they affect the human genome: from generating insertion mutations and genomic instability to altering gene expression and contributing to genetic innovation. Increasingly detailed analyses of human and other primate genomes are revealing the scale and complexity of the past and current contributions of non-LTR retrotransposons to genomic change in the human lineage.

Many LINE1 elements contribute to the transcriptome of human somatic cells

Over 600 LINE 1 elements are shown to be transcribed in humans; 400 of these are full-length elements in the reference genome.

Background

While LINE1 (L1) retroelements comprise nearly 20% of the human genome, the majority are thought to have been rendered transcriptionally inactive, due to either mutation or epigenetic suppression. How many L1 elements 'escape' these forms of repression and contribute to the transcriptome of human somatic cells? We have cloned out expressed sequence tags corresponding to the 5' and 3' flanks of L1 elements in order to characterize the population of elements that are being actively transcribed. We also examined expression of a select number of elements in different individuals.

Results

We isolated expressed sequence tags from human lymphoblastoid cell lines corresponding to 692 distinct L1 element sites, including 410 full-length elements. Four of the expression tagged sites corresponding to full-length elements from the human specific L1Hs subfamily were examined in European-American individuals and found to be differentially expressed in different family members.

Conclusions

A large number of different L1 element sites are expressed in human somatic tissues, and this expression varies among different individuals. Paradoxically, few elements were tagged at high frequency, indicating that the majority of expressed L1s are transcribed at low levels. Based on our preliminary expression studies of a limited number of elements in a single family, we predict a significant degree of inter-individual transcript-level polymorphism in this class of sequence.

Pervasive transcription of the eukaryotic genome: functional indices and conceptual implications

Genome-wide analyses of the eukaryotic transcriptome have revealed that the majority of the genome is transcribed, producing large numbers of non-protein-coding RNAs (ncRNAs). This surprising observation challenges many assumptions about the genetic programming of higher organisms and how information is stored and organized within the genome. Moreover, the rapid advances in genomics have given little opportunity for biologists to integrate these emerging findings into their intellectual and experimental frameworks. This problem has been compounded by the perception that genome-wide studies often generate more questions than answers, which in turn has led to confusion and controversy. In this article, we address common questions associated with the phenomenon of pervasive transcription and consider the indices that can be used to evaluate the function (or lack thereof) of the resulting ncRNAs. We suggest that many lines of evidence, including expression profiles, conservation signatures, chromatin modification patterns and examination of increasing numbers of individual cases, argue in favour of the widespread functionality of non-coding transcription. We also discuss how informatic and experimental approaches used to analyse protein-coding genes may not be applicable to ncRNAs and how the general perception that protein-coding genes form the main informational output of the genome has resulted in much of the misunderstanding surrounding pervasive transcription and its potential significance. Finally, we present the conceptual implications of the majority of the eukaryotic genome being functional and describe how appreciating this perspective will provide considerable opportunity to further understand the molecular basis of development and complex diseases.