Intronic L1 Retrotransposons and Nested Genes Cause Transcriptional Interference by Inducing Intron Retention, Exonization and Cryptic Polyadenylation

Transcriptional interference has been recently recognized as an unexpectedly complex and mostly negative regulation of genes. Despite a relatively few studies that emerged in recent years, it has been demonstrated that a readthrough transcription derived from one gene can influence the transcription of another overlapping or nested gene. However, the molecular effects resulting from this interaction are largely unknown.Using in silico chromosome walking, we searched for prematurely terminated transcripts bearing signatures of intron retention or exonization of intronic sequence at their 3' ends upstream to human L1 retrotransposons, protein-coding and noncoding nested genes. We demonstrate that transcriptional interference induced by intronic L1s (or other repeated DNAs) and nested genes could be characterized by intron retention, forced exonization and cryptic polyadenylation. These molecular effects were revealed from the analysis of endogenous transcripts derived from different cell lines and tissues and confirmed by the expression of three minigenes in cell culture. While intron retention and exonization were comparably observed in introns upstream to L1s, forced exonization was preferentially detected in nested genes. Transcriptional interference induced by L1 or nested genes was dependent on the presence or absence of cryptic splice sites, affected the inclusion or exclusion of the upstream exon and the use of cryptic polyadenylation signals.Our results suggest that transcriptional interference induced by intronic L1s and nested genes could influence the transcription of the large number of genes in normal as well as in tumor tissues. Therefore, this type of interference could have a major impact on the regulation of the host gene expression

Improvement of RNA secondary structure prediction using RNase H cleavage and randomized oligonucleotides

RNA secondary structure prediction using free energy minimization is one method to gain an approximation of structure. Constraints generated by enzymatic mapping or chemical modification can improve the accuracy of secondary structure prediction. We report a facile method that identifies single-stranded regions in RNA using short, randomized DNA oligonucleotides and RNase H cleavage. These regions are then used as constraints in secondary structure prediction. This method was used to improve the secondary structure prediction of Escherichia coli 5S rRNA. The lowest free energy structure without constraints has only 27% of the base pairs present in the phylogenetic structure. The addition of constraints from RNase H cleavage improves the prediction to 100% of base pairs. The same method was used to generate secondary structure constraints for yeast tRNAPhe, which is accurately predicted in the absence of constraints (95%). Although RNase H mapping does not improve secondary structure prediction, it does eliminate all other suboptimal structures predicted within 10% of the lowest free energy structure. The method is advantageous over other single-stranded nucleases since RNase H is functional in physiological conditions. Moreover, it can be used for any RNA to identify accessible binding sites for oligonucleotides or small molecules

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

Disease Dynamics and Bird Migration—Linking Mallards Anas platyrhynchos and Subtype Diversity of the Influenza A Virus in Time and Space

The mallard Anas platyrhynchos is a reservoir species for influenza A virus in the northern hemisphere, with particularly high prevalence rates prior to as well as during its prolonged autumn migration. It has been proposed that the virus is brought from the breeding grounds and transmitted to conspecifics during subsequent staging during migration, and so a better understanding of the natal origin of staging ducks is vital to deciphering the dynamics of viral movement pathways. Ottenby is an important stopover site in southeast Sweden almost halfway downstream in the major Northwest European flyway, and is used by millions of waterfowl each year. Here, mallards were captured and sampled for influenza A virus infection, and positive samples were subtyped in order to study possible links to the natal area, which were determined by a novel approach combining banding recovery data and isotopic measurements (δ2H) of feathers grown on breeding grounds. Geographic assignments showed that the core natal areas of studied mallards were in Estonia, southern and central Finland, and northwestern Russia. This study demonstrates a clear temporal succession of latitudes of natal origin during the course of autumn migration. We also demonstrate a corresponding and concomitant shift in virus subtypes. Acknowledging that these two different patterns were based in part upon different data, a likely interpretation worth further testing is that the early arriving birds with more proximate origins have different influenza A subtypes than the more distantly originating late autumn birds. If true, this knowledge would allow novel insight into the origins and transmission of the influenza A virus among migratory hosts previously unavailable through conventional approaches

An Amphioxus Gli Gene Reveals Conservation of Midline Patterning and the Evolution of Hedgehog Signalling Diversity in Chordates

Background. Hedgehog signalling, interpreted in receiving cells by Gli transcription factors, plays a central role in the development of vertebrate and Drosphila embryos. Many aspects of the signalling pathway are conserved between these lineages, however vertebrates have diverged in at least one key aspect: they have evolved multiple Gli genes encoding functionally-distinct proteins, increasing the complexity of the hedgehog-dependent transcriptional response. Amphioxus is one of the closest living relatives of the vertebrates, having split from the vertebrate lineage prior to the widespread gene duplication prominent in early vertebrate evolution. Principal findings. We show that amphioxus has a single Gli gene, which is deployed in tissues adjacent to sources of hedgehog signalling derived from the midline and anterior endoderm. This shows the duplication and divergence of the Gli family, and hence the origin of vertebrate Gli functional diversity, was specific to the vertebrate lineage. However we also show that the single amphioxus Gli gene produces two distinct transcripts encoding different proteins. We utilise three tests of Gli function to examine the transcription regulatory capacities of these different proteins, demonstrating one has activating activity similar to Gli2, while the other acts as a weak repressor, similar to Gli3. Conclusions. These data show that the vertebrates and amphioxus have evolved functionally-similar repertoires of Gli proteins using parallel molecular routes; vertebrates via gene duplication and divergence, and amphioxus via alternate splicing of a single gene. Our results demonstrate that similar functional complexity of intercellular signalling can be achieved via different evolutionary pathways

Does Selection against Transcriptional Interference Shape Retroelement-Free Regions in Mammalian Genomes?

BACKGROUND: Eukaryotic genomes are scattered with retroelements that proliferate through retrotransposition. Although retroelements make up around 40 percent of the human genome, large regions are found to be completely devoid of retroelements. This has been hypothesised to be a result of genomic regions being intolerant to insertions of retroelements. The inadvertent transcriptional activity of retroelements may affect neighbouring genes, which in turn could be detrimental to an organism. We speculate that such retroelement transcription, or transcriptional interference, is a contributing factor in generating and maintaining retroelement-free regions in the human genome. METHODOLOGY/PRINCIPAL FINDINGS: Based on the known transcriptional properties of retroelements, we expect long interspersed elements (LINEs) to be able to display a high degree of transcriptional interference. In contrast, we expect short interspersed elements (SINEs) to display very low levels of transcriptional interference. We find that genomic regions devoid of long interspersed elements (LINEs) are enriched for protein-coding genes, but that this is not the case for regions devoid of short interspersed elements (SINEs). This is expected if genes are subject to selection against transcriptional interference. We do not find microRNAs to be associated with genomic regions devoid of either SINEs or LINEs. We further observe an increased relative activity of genes overlapping LINE-free regions during early embryogenesis, where activity of LINEs has been identified previously. CONCLUSIONS/SIGNIFICANCE: Our observations are consistent with the notion that selection against transcriptional interference has contributed to the maintenance and/or generation of retroelement-free regions in the human genome

RISCI - Repeat Induced Sequence Changes Identifier: a comprehensive, comparative genomics-based, in silico subtractive hybridization pipeline to identify repeat induced sequence changes in closely related genomes

<p>Abstract</p> <p>Background -</p> <p>The availability of multiple whole genome sequences has facilitated <it>in silico </it>identification of fixed and polymorphic transposable elements (TE). Whereas polymorphic loci serve as makers for phylogenetic and forensic analysis, fixed species-specific transposon insertions, when compared to orthologous loci in other closely related species, may give insights into their evolutionary significance. Besides, TE insertions are not isolated events and are frequently associated with subtle sequence changes concurrent with insertion or post insertion. These include duplication of target site, 3' and 5' flank transduction, deletion of the target locus, 5' truncation or partial deletion and inversion of the transposon, and post insertion changes like inter or intra element recombination, disruption etc. Although such changes have been studied independently, no automated platform to identify differential transposon insertions and the associated array of sequence changes in genomes of the same or closely related species is available till date. To this end, we have designed RISCI - 'Repeat Induced Sequence Changes Identifier' - a comprehensive, comparative genomics-based, <it>in silico </it>subtractive hybridization pipeline to identify differential transposon insertions and associated sequence changes using specific alignment signatures, which may then be examined for their downstream effects.</p> <p>Results -</p> <p>We showcase the utility of RISCI by comparing full length and truncated L1HS and AluYa5 retrotransposons in the reference human genome with the chimpanzee genome and the alternate human assemblies (Celera and HuRef). Comparison of the reference human genome with alternate human assemblies using RISCI predicts 14 novel polymorphisms in full length L1HS, 24 in truncated L1HS and 140 novel polymorphisms in AluYa5 insertions, besides several insertion and post insertion changes. We present comparison with two previous studies to show that RISCI predictions are broadly in agreement with earlier reports. We also demonstrate its versatility by comparing various strains of <it>Mycobacterium tuberculosis </it>for IS 6100 insertion polymorphism.</p> <p>Conclusions -</p> <p>RISCI combines comparative genomics with subtractive hybridization, inferring changes only when exclusive to one of the two genomes being compared. The pipeline is generic and may be applied to most transposons and to any two or more genomes sharing high sequence similarity. Such comparisons, when performed on a larger scale, may pull out a few critical events, which may have seeded the divergence between the two species under comparison.</p

Long interspersed nuclear element-1 hypomethylation in cancer: biology and clinical applications

Epigenetic changes in long interspersed nuclear element-1s (LINE-1s or L1s) occur early during the process of carcinogenesis. A lower methylation level (hypomethylation) of LINE-1 is common in most cancers, and the methylation level is further decreased in more advanced cancers. Consequently, several previous studies have suggested the use of LINE-1 hypomethylation levels in cancer screening, risk assessment, tumor staging, and prognostic prediction. Epigenomic changes are complex, and global hypomethylation influences LINE-1s in a generalized fashion. However, the methylation levels of some loci are dependent on their locations. The consequences of LINE-1 hypomethylation are genomic instability and alteration of gene expression. There are several mechanisms that promote both of these consequences in cis. Therefore, the methylation levels of different sets of LINE-1s may represent certain phenotypes. Furthermore, the methylation levels of specific sets of LINE-1s may indicate carcinogenesis-dependent hypomethylation. LINE-1 methylation pattern analysis can classify LINE-1s into one of three classes based on the number of methylated CpG dinucleotides. These classes include hypermethylation, partial methylation, and hypomethylation. The number of partial and hypermethylated loci, but not hypomethylated LINE-1s, is different among normal cell types. Consequently, the number of hypomethylated loci is a more promising marker than methylation level in the detection of cancer DNA. Further genome-wide studies to measure the methylation level of each LINE-1 locus may improve PCR-based methylation analysis to allow for a more specific and sensitive detection of cancer DNA or for an analysis of certain cancer phenotypes

Hypomethylation of Intragenic LINE-1 Represses Transcription in Cancer Cells through AGO2

In human cancers, the methylation of long interspersed nuclear element -1 (LINE-1 or L1) retrotransposons is reduced. This occurs within the context of genome wide hypomethylation, and although it is common, its role is poorly understood. L1s are widely distributed both inside and outside of genes, intragenic and intergenic, respectively. Interestingly, the insertion of active full-length L1 sequences into host gene introns disrupts gene expression. Here, we evaluated if intragenic L1 hypomethylation influences their host gene expression in cancer. First, we extracted data from L1base (http://l1base.molgen.mpg.de), a database containing putatively active L1 insertions, and compared intragenic and intergenic L1 characters. We found that intragenic L1 sequences have been conserved across evolutionary time with respect to transcriptional activity and CpG dinucleotide sites for mammalian DNA methylation. Then, we compared regulated mRNA levels of cells from two different experiments available from Gene Expression Omnibus (GEO), a database repository of high throughput gene expression data, (http://www.ncbi.nlm.nih.gov/geo) by chi-square. The odds ratio of down-regulated genes between demethylated normal bronchial epithelium and lung cancer was high (p<1E−27; OR = 3.14; 95% CI = 2.54–3.88), suggesting cancer genome wide hypomethylation down-regulating gene expression. Comprehensive analysis between L1 locations and gene expression showed that expression of genes containing L1s had a significantly higher likelihood to be repressed in cancer and hypomethylated normal cells. In contrast, many mRNAs derived from genes containing L1s are elevated in Argonaute 2 (AGO2 or EIF2C2)-depleted cells. Hypomethylated L1s increase L1 mRNA levels. Finally, we found that AGO2 targets intronic L1 pre-mRNA complexes and represses cancer genes. These findings represent one of the mechanisms of cancer genome wide hypomethylation altering gene expression. Hypomethylated intragenic L1s are a nuclear siRNA mediated cis-regulatory element that can repress genes. This epigenetic regulation of retrotransposons likely influences many aspects of genomic biology

GLI1 Confers Profound Phenotypic Changes upon LNCaP Prostate Cancer Cells That Include the Acquisition of a Hormone Independent State

The GLI (GLI1/GLI2) transcription factors have been implicated in the development and progression of prostate cancer although our understanding of how they actually contribute to the biology of these common tumours is limited. We observed that GLI reporter activity was higher in normal (PNT-2) and tumourigenic (DU145 and PC-3) androgen-independent cells compared to androgen-dependent LNCaP prostate cancer cells and, accordingly, GLI mRNA levels were also elevated. Ectopic expression of GLI1 or the constitutively active ΔNGLI2 mutant induced a distinct cobblestone-like morphology in LNCaP cells that, regarding the former, correlated with increased GLI2 as well as expression of the basal/stem-like markers CD44, β1-integrin, ΔNp63 and BMI1, and decreased expression of the luminal marker AR (androgen receptor). LNCaP-GLI1 cells were viable in the presence of the AR inhibitor bicalutamide and gene expression profiling revealed that the transcriptome of LNCaP-GLI1 cells was significantly closer to DU145 and PC-3 cells than to control LNCaP-pBP (empty vector) cells, as well as identifying LCN2/NGAL as a highly induced transcript which is associated with hormone independence in breast and prostate cancer. Functionally, LNCaP-GLI1 cells displayed greater clonal growth and were more invasive than control cells but they did not form colonies in soft agar or prostaspheres in suspension suggesting that they do not possess inherent stem cell properties. Moreover, targeted suppression of GLI1 or GLI2 with siRNA did not reverse the transformed phenotype of LNCaP-GLI1 cells nor did double GLI1/GLI2 knockdowns activate AR expression in DU145 or PC-3 cells. As such, early targeting of the GLI oncoproteins may hinder progression to a hormone independent state but a more detailed understanding of the mechanisms that maintain this phenotype is required to determine if their inhibition will enhance the efficacy of anti-hormonal therapy through the induction of a luminal phenotype and increased dependency upon AR function