Patterns of Transposable Element Expression and Insertion in Cancer

Anti-HERV-K Drugs and Vaccines, Possible Therapies against Tumors

The footprint of human endogenous retroviruses (HERV), specifically HERV-K, has been found in malignancies, such as melanoma, teratocarcinoma, osteosarcoma, breast cancer, lymphoma, and ovary and prostate cancers. HERV-K is characterized as the most biologically active HERV due to possession of open reading frames (ORF) for all Gag, Pol, and Env genes, which enables it to be more infective and obstructive towards specific cell lines and other exogenous viruses, respectively. Some factors might contribute to carcinogenicity and at least one of them has been recognized in various tumors, including overexpression/methylation of long interspersed nuclear element 1 (LINE-1), HERV-K Gag, and Env genes themselves plus their transcripts and protein products, and HERV-K reverse transcriptase (RT). Therapies effective for HERV-K-associated tumors mostly target invasive autoimmune responses or growth of tumors through suppression of HERV-K Gag or Env protein and RT. To design new therapeutic options, more studies are needed to better understand whether HERV-K and its products (Gag/Env transcripts and HERV-K proteins/RT) are the initiators of tumor formation or just the disorder’s developers. Accordingly, this review aims to present evidence that highlights the association between HERV-K and tumorigenicity and introduces some of the available or potential therapies against HERV-K-induced tumors.

Repetitive Sequence Transcription in Breast Cancer

Repetitive sequences represent about half of the human genome. They are actively transcribed and play a role during development and in epigenetic regulation. The altered activity of repetitive sequences can lead to genomic instability and they can contribute to the establishment or the progression of degenerative diseases and cancer transformation. In this work, we analyzed the expression profiles of DNA repetitive sequences in the breast cancer specimens of the HMUCC cohort. Satellite expression is generally upregulated in breast cancers, with specific families upregulated per histotype: in HER2-enriched cancers, they are the human satellite II (HSATII), in luminal A and B, they are part of the ALR family and in triple-negative, they are part of SAR and GSAT families, together with a perturbation in the transcription from endogenous retroviruses and their LTR sequences. We report that the background expression of repetitive sequences in healthy tissues of cancer patients differs from the tissues of non-cancerous controls. To conclude, peculiar patterns of expression of repetitive sequences are reported in each specimen, especially in the case of transcripts arising from satellite repeats.

Hsa-miR-422a Originated from Short Interspersed Nuclear Element Increases ARID5B Expression by Collaborating with NF-E2

MicroRNAs (miRNAs) are a class of small non-coding RNAs that regulate the expression of target messenger RNA (mRNA) complementary to the 3' untranslated region (UTR) at the post-transcriptional level. Hsa-miR-422a, which is commonly known as miRNA derived from transposable element (MDTE), was derived from short interspersed nuclear element (SINE). Through expression analysis, hsa-miR-422a was found to be highly expressed in both the small intestine and liver of crab-eating monkey. AT-Rich Interaction Domain 5 B (ARID5B) was selected as the target gene of hsa-miR-422a, which has two binding sites in both the exon and 3'UTR of ARID5B. To identify the interaction between hsa-miR-422a and ARID5B, a dual luciferase assay was conducted in HepG2 cell line. The luciferase activity of cells treated with the hsa-miR-422a mimic was upregulated and inversely downregulated when both the hsa-miR-422a mimic and inhibitor were administered. Nuclear factor erythroid-2 (NF-E2) was selected as the core transcription factor (TF) via feed forward loop analysis. The luciferase expression was downregulated when both the hsa-miR-422a mimic and siRNA of NF-E2 were treated, compared to the treatment of the hsa-miR-422a mimic alone. The present study suggests that hsa-miR-422a derived from SINE could bind to the exon region as well as the 3'UTR of ARID5B. Additionally, hsa-miR-422a was found to share binding sites in ARID5Bwith several TFs, including NF-E2. The hsa-miR-422a might thus interact with TF to regulate the expression of ARID5B, as demonstrated experimentally. Altogether, hsa-miR-422a acts as a super enhancer miRNA of ARID5Bby collaborating with TF and NF-E2.

TE Density: a tool to investigate the biology of transposable elements

BACKGROUND

Transposable elements (TEs) are powerful creators of genotypic and phenotypic diversity due to their inherent mutagenic capabilities and in this way they serve as a deep reservoir of sequences for genomic variation. As agents of genetic disruption, a TE's potential to impact phenotype is partially a factor of its location in the genome. Previous research has shown TEs' ability to impact the expression of neighboring genes, however our understanding of this trend is hampered by the exceptional amount of diversity in the TE world, and a lack of publicly available computational methods that quantify the presence of TEs relative to genes.

RESULTS

Here, we have developed a tool to more easily quantify TE presence relative to genes through the use of only a gene and TE annotation, yielding a new metric we call TE Density. Briefly defined as the proportion of TE-occupied base-pairs relative to a window-size of the genome. This new pipeline reports TE density for each gene in the genome, for each type descriptor of TE (order and superfamily), and for multiple positions and distances relative to the gene (upstream, intragenic, and downstream) over sliding, user-defined windows. In this way, we overcome previous limitations to the study of TE-gene relationships by focusing on all TE types present in the genome, utilizing flexible genomic distances for measurement, and reporting a TE presence metric for every gene in the genome.

CONCLUSIONS

Together, this new tool opens up new avenues for studying TE-gene relationships, genome architecture, comparative genomics, and the tremendous diversity present of the TE world. TE Density is open-source and freely available at: https://github.com/sjteresi/TE_Density .

Transposable Elements and Human Diseases: Mechanisms and Implication in the Response to Environmental Pollutants

Transposable elements (TEs) are recognized as major players in genome plasticity and evolution. The high abundance of TEs in the human genome, especially the Alu and Long Interspersed Nuclear Element-1 (LINE-1) repeats, makes them responsible for the molecular origin of several diseases. This involves several molecular mechanisms that are presented in this review: insertional mutation, DNA recombination and chromosomal rearrangements, modification of gene expression, as well as alteration of epigenetic regulations. This literature review also presents some of the more recent and/or more classical examples of human diseases in which TEs are involved. Whether through insertion of LINE-1 or Alu elements that cause chromosomal rearrangements, or through epigenetic modifications, TEs are widely implicated in the origin of human cancers. Many other human diseases can have a molecular origin in TE-mediated chromosomal recombination or alteration of gene structure and/or expression. These diseases are very diverse and include hemoglobinopathies, metabolic and neurological diseases, and common diseases. Moreover, TEs can also have an impact on aging. Finally, the exposure of individuals to stresses and environmental contaminants seems to have a non-negligible impact on the epigenetic derepression and mobility of TEs, which can lead to the development of diseases. Thus, improving our knowledge of TEs may lead to new potential diagnostic markers of diseases.

The insertion and dysregulation of transposable elements in osteosarcoma and their association with patient event-free survival

The dysregulation of transposable elements (TEs) has been explored in a variety of cancers. However, TE activities in osteosarcoma (OS) have not been extensively studied yet. By integrative analysis of RNA-seq, whole-genome sequencing (WGS), and methylation data, we showed aberrant TE activities associated with dysregulations of TEs in OS tumors. Specifically, expression levels of LINE-1 and Alu of different evolutionary ages, as well as subfamilies of SVA and HERV-K, were significantly up-regulated in OS tumors, accompanied by enhanced DNA repair responses. We verified the characteristics of LINE-1 mediated TE insertions, including target site duplication (TSD) length (centered around 15 bp) and preferential insertions into intergenic and AT-rich regions as well as intronic regions of longer genes. By filtering polymorphic TE insertions reported in 1000 genome project (1KGP), besides 148 tumor-specific somatic TE insertions, we found most OS patient-specific TE insertions (3175 out of 3326) are germline insertions, which are associated with genes involved in neuronal processes or with transcription factors important for cancer development. In addition to 68 TE-affected cancer genes, we found recurrent germline TE insertions in 72 non-cancer genes with high frequencies among patients. We also found that +/− 500 bps flanking regions of transcription start sites (TSS) of LINE-1 (young) and Alu showed lower methylation levels in OS tumor samples than controls. Interestingly, by incorporating patient clinical data and focusing on TE activities in OS tumors, our data analysis suggested that higher TE insertions in OS tumors are associated with a longer event-free survival time.

An assessment of bioinformatics tools for the detection of human endogenous retroviral insertions in short-read genome sequencing data

There is a growing interest in the study of human endogenous retroviruses (HERVs) given the substantial body of evidence that implicates them in many human diseases. Although their genomic characterization presents numerous technical challenges, next-generation sequencing (NGS) has shown potential to detect HERV insertions and their polymorphisms in humans. Currently, a number of computational tools to detect them in short-read NGS data exist. In order to design optimal analysis pipelines, an independent evaluation of the available tools is required. We evaluated the performance of a set of such tools using a variety of experimental designs and datasets. These included 50 human short-read whole-genome sequencing samples, matching long and short-read sequencing data, and simulated short-read NGS data. Our results highlight a great performance variability of the tools across the datasets and suggest that different tools might be suitable for different study designs. However, specialized tools designed to detect exclusively human endogenous retroviruses consistently outperformed generalist tools that detect a wider range of transposable elements. We suggest that, if sufficient computing resources are available, using multiple HERV detection tools to obtain a consensus set of insertion loci may be ideal. Furthermore, given that the false positive discovery rate of the tools varied between 8% and 55% across tools and datasets, we recommend the wet lab validation of predicted insertions if DNA samples are available.

Association of HERV-K and LINE-1 hypomethylation with reduced disease-free survival in melanoma patients

Aim: To evaluate CpG methylation of long interspersed nuclear elements 1 (LINE-1) and human endogenous retrovirus K (HERV-K) retroelements as potential prognostic biomarkers in cutaneous melanoma. Materials & methods: Methylation of HERV-K and LINE-1 retroelements was assessed in resected melanoma tissues from 82 patients ranging in age from 14 to 88 years. In addition, nevi from eight patients were included for comparison with nonmalignant melanocytic lesions. Results: Methylation levels were lower in melanomas than in nevi. HERV-K and LINE-1 methylation were decreased in melanoma patients with clinical parameters associated with adverse prognosis, while they were independent of age and gender. Hypomethylation of HERV-K (but not LINE-1) was an independent predictor of reduced disease-free survival. Conclusion: HERV-K hypomethylation can be a potential independent biomarker of melanoma recurrence.

The Sophisticated Transcriptional Response Governed by Transposable Elements in Human Health and Disease

Transposable elements (TEs), which cover ~45% of the human genome, although firstly considered as “selfish” DNA, are nowadays recognized as driving forces in eukaryotic genome evolution. This capability resides in generating a plethora of sophisticated RNA regulatory networks that influence the cell type specific transcriptome in health and disease. Indeed, TEs are transcribed and their RNAs mediate multi-layered transcriptional regulatory functions in cellular identity establishment, but also in the regulation of cellular plasticity and adaptability to environmental cues, as occurs in the immune response. Moreover, TEs transcriptional deregulation also evolved to promote pathogenesis, as in autoimmune and inflammatory diseases and cancers. Importantly, many of these findings have been achieved through the employment of Next Generation Sequencing (NGS) technologies and bioinformatic tools that are in continuous improvement to overcome the limitations of analyzing TEs sequences. However, they are highly homologous, and their annotation is still ambiguous. Here, we will review some of the most recent findings, questions and improvements to study at high resolution this intriguing portion of the human genome in health and diseases, opening the scenario to novel therapeutic opportunities.

An atlas of transposable element-derived alternative splicing in cancer

Transposable element (TE)-derived sequences comprise more than half of the human genome, and their presence has been documented to alter gene expression in a number of different ways, including the generation of alternatively spliced transcript isoforms. Alternative splicing has been associated with tumorigenesis for a number of different cancers. The objective of this study was to broadly characterize the role of human TEs in generating alternatively spliced transcript isoforms in cancer. To do so, we screened for the presence of TE-derived sequences co-located with alternative splice sites that are differentially used in normal versus cancer tissues. We analysed a comprehensive set of alternative splice variants characterized for 614 matched normal-tumour tissue pairs across 13 cancer types, resulting in the discovery of 4820 TE-generated alternative splice events distributed among 723 cancer-associated genes. Short interspersed nuclear elements (Alu) and long interspersed nuclear elements (L1) were found to contribute the majority of TE-generated alternative splice sites in cancer genes. A number of cancer-associated genes, including MYH11, WHSC1 and CANT1, were shown to have overexpressed TE-derived isoforms across a range of cancer types. TE-derived isoforms were also linked to cancer-specific fusion transcripts, suggesting a novel mechanism for the generation of transcriptome diversity via trans-splicing mediated by dispersed TE repeats. This article is part of a discussion meeting issue 'Crossroads between transposons and gene regulation'.

Biomarker potential of repetitive-element transcriptome in lung cancer

Since repetitive elements (REs) account for nearly 53% of the human genome, profiling its transcription after an oncogenic change might help in the search for new biomarkers. Lung cancer was selected as target since it is the most frequent cause of cancer death. A bioinformatic workflow based on well-established bioinformatic tools (such as RepEnrich, RepBase, SAMTools, edgeR and DESeq2) has been developed to identify differentially expressed RNAs from REs. It was trained and tested with public RNA-seq data from matched sequencing of tumour and healthy lung tissues from the same patient to reveal differential expression within the RE transcriptome. Healthy lung tissues express a specific set of REs whose expression, after an oncogenic process, is strictly and specifically changed. Discrete sets of differentially expressed REs were found for lung adenocarcinoma, for small-cell lung cancer, and for both cancers. Differential expression affects more HERV-than LINE-derived REs and seems biased towards down-regulation in cancer cells. REs behaving consistently in all patients were tested in a different patient cohort to validate the proposed biomarkers. Down-regulation of AluYg6 and LTR18B was confirmed as potential lung cancer biomarkers, while up-regulation of HERVK11D-Int is specific for lung adenocarcinoma and up-regulation of UCON88 is specific for small cell lung cancer. Hence, the study of RE transcriptome might be considered another research target in cancer, making REs a promising source of lung cancer biomarkers.

Transposable element expression in tumors is associated with immune infiltration and increased antigenicity

Profound global loss of DNA methylation is a hallmark of many cancers. One potential consequence of this is the reactivation of transposable elements (TEs) which could stimulate the immune system via cell-intrinsic antiviral responses. Here, we develop REdiscoverTE, a computational method for quantifying genome-wide TE expression in RNA sequencing data. Using The Cancer Genome Atlas database, we observe increased expression of over 400 TE subfamilies, of which 262 appear to result from a proximal loss of DNA methylation. The most recurrent TEs are among the evolutionarily youngest in the genome, predominantly expressed from intergenic loci, and associated with antiviral or DNA damage responses. Treatment of glioblastoma cells with a demethylation agent results in both increased TE expression and de novo presentation of TE-derived peptides on MHC class I molecules. Therapeutic reactivation of tumor-specific TEs may synergize with immunotherapy by inducing inflammation and the display of potentially immunogenic neoantigens.

Treatment with demethylation agents can reactivate transposable elements. Here in glioblastoma, the authors also show that this is accompanied by de novo presentation of TE-derived peptides on MHC class I molecules.

Transposable Element Expression in Acute Myeloid Leukemia Transcriptome and Prognosis

Over half of the human genome is comprised of transposable elements (TE). Despite large-scale studies of the transcriptome in cancer, a comprehensive look at TE expression and its relationship to various mutations or prognosis has not been performed. We characterized the expression of TE in 178 adult acute myeloid leukemia (AML) patients using transcriptome data from The Cancer Genome Atlas. We characterized mutation specific dysregulation of TE expression using a multivariate linear model. We identified distinct patterns of TE expression associated with specific mutations and transcriptional networks. Genes regulating methylation was not associated with significant change in TE expression. Using an unpenalized cox regression analysis we identified a TE expression signature that predicted prognosis in AML. We identified 14 candidate prognostic TE transcripts (TEP) that classified AML as high/low-risk and this was independent of mutation-based and coding-gene expression based risk-stratification. TEP was able to predict prognosis in independent cohorts of 284 pediatric AML patients and 19 relapsed adult AML patients. This first comprehensive study of TE expression in AML demonstrates that TE expression can serve as a biomarker for prognosis in AML, and provides novel insights into the biology of AML. Studies characterizing its role in other cancers are warranted.

Transposable element dysregulation in systemic lupus erythematosus and regulation by histone conformation and Hsp90

Systemic lupus erythematosus (SLE) represents an autoimmune disease in which activation of the type I interferon pathway leads to dysregulation of tolerance and the generation of autoantibodies directed against nuclear constituents. The mechanisms driving the activation of the interferon pathway in SLE have been the subject of intense investigation but are still incompletely understood. Transposable elements represent an enormous source of RNA that could potentially stimulate the cell intrinsic RNA-recognition pathway, leading to upregulation of interferons. We used RNA-seq to define transposable element families and subfamilies in three cell types in SLE and found diverse effects on transposable element expression in the three cell types and even within a given family of transposable elements. When potential mechanisms were examined, we found that Hsp90 inhibition could drive increased expression of multiple type of transposable elements. Both direct inhibition and the delivery of a heat shock itself, which redirects heat shock regulators (including Hsp90) off of basal expression promoters and onto heat shock-responsive promoters, led to increased transposable element expression. This effect was amplified by the concurrent delivery of a histone deacetylase inhibitor. We conclude that transposable elements are dysregulated in SLE and there are tissue-specific effects and locus-specific effects. The magnitude of RNAs attributable to transposable elements makes their dysregulation of critical interest in SLE where transposable element RNA complexed with proteins has been shown to drive interferon expression.

Enhancer reprogramming in tumor progression: a new route towards cancer cell plasticity

Cancer heterogeneity arises during tumor progression as a consequence of genetic insults, environmental cues, and reversible changes in the epigenetic state, favoring tumor cell plasticity. The role of enhancer reprogramming is emerging as a relevant field in cancer biology as it supports adaptation of cancer cells to those environmental changes encountered during tumor progression and metastasis seeding. In this review, we describe the cancer-related alterations that drive oncogenic enhancer activity, leading to dysregulated transcriptional programs. We discuss the molecular mechanisms of both cis- and trans-factors in overriding the regulatory circuits that maintain cell-type specificity and imposing an alternative, de-regulated enhancer activity in cancer cells. We further comment on the increasing evidence which implicates stress response and aging-signaling pathways in the enhancer landscape reprogramming during tumorigenesis. Finally, we focus on the potential therapeutic implications of these enhancer-mediated subverted transcriptional programs, putting particular emphasis on the lack of information regarding tumor progression and the metastatic outgrowth, which still remain the major cause of mortality related to cancer.

Identification of transposable elements fused in the exonic region of the olive flounder genome

Transposable elements (TEs) are mobile genetic sequences that comprise a large portion of vertebrate genomes. The olive flounder (Paralichthys olivaceus) is a valuable marine resource in East Asia. The scope of most genomic studies on the olive flounder is limited to its immunology as their focus is the prevention of mass mortality of this species. Thus, for a broader understanding of the species, its genomic information is consistently in demand. Transcripts sequences were acquired from transcriptome analysis using gill tissues of 12 olive flounders. Distribution of TEs inserted in exonic region of the olive flounder genome was analyzed using RepeatMasker ( http://www.repeatmasker.org/ ). We found 1140 TEs in the exonic region of the genome and long interspersed nuclear elements (LINEs) and long terminal repeats (LTRs) insertions occurred with forward orientation preferences. Transposons belonging to the hAt, Gypsy, and LINE 1 (L1) subfamilies were the most abundant DNA transposons, LTRs, and long interspersed elements (LINEs), respectively. Finally, we carried out a gene ontology analysis to determine the function of TE-fused genes. These results provide some genomic information about TEs that is useful for future research on changes in properties and functions of genes by TEs in the olive flounder genome.

The epigenetic alterations of endogenous retroelements in aging

Endogenous retroelements, transposons that mobilize through RNA intermediates, include some of the most abundant repetitive sequences of the human genome, such as Alu and LINE-1 sequences, and human endogenous retroviruses. Recent discoveries demonstrate that these mobile genetic elements not only act as intragenomic parasites, but also exert regulatory roles in living cells. The risk of genomic instability represented by endogenous retroelements is normally counteracted by a series of epigenetic control mechanisms which include, among the most important, CpG DNA methylation. Indeed, most of the genomic CpG sites subjected to DNA methylation in the nuclear DNA are carried by these repetitive elements. As other parts of the genome, endogenous retroelements and other transposable elements are subjected to deep epigenetic alterations during aging, repeatedly observed in the context of organismal and cellular senescence, in human and other species. This review summarizes the current status of knowledge about the epigenetic alterations occurring in this large, non-genic portion of the genome in aging and age-related conditions, with a focus on the causes and the possible functional consequences of these alterations.

Human Retrotransposon Insertion Polymorphisms Are Associated with Health and Disease via Gene Regulatory Phenotypes

The human genome hosts several active families of transposable elements (TEs), including the Alu, LINE-1, and SVA retrotransposons that are mobilized via reverse transcription of RNA intermediates. We evaluated how insertion polymorphisms generated by human retrotransposon activity may be related to common health and disease phenotypes that have been previously interrogated through genome-wide association studies (GWAS). To address this question, we performed a genome-wide screen for retrotransposon polymorphism disease associations that are linked to TE induced gene regulatory changes. Our screen first identified polymorphic retrotransposon insertions found in linkage disequilibrium (LD) with single nucleotide polymorphisms that were previously associated with common complex diseases by GWAS. We further narrowed this set of candidate disease associated retrotransposon polymorphisms by identifying insertions that are located within tissue-specific enhancer elements. We then performed expression quantitative trait loci analysis on the remaining set of candidates in order to identify polymorphic retrotransposon insertions that are associated with gene expression changes in B-cells of the human immune system. This progressive and stringent screen yielded a list of six retrotransposon insertions as the strongest candidates for TE polymorphisms that lead to disease via enhancer-mediated changes in gene regulation. For example, we found an SVA insertion within a cell-type specific enhancer located in the second intron of the B4GALT1 gene. B4GALT1 encodes a glycosyltransferase that functions in the glycosylation of the Immunoglobulin G (IgG) antibody in such a way as to convert its activity from pro- to anti-inflammatory. The disruption of the B4GALT1 enhancer by the SVA insertion is associated with down-regulation of the gene in B-cells, which would serve to keep the IgG molecule in a pro-inflammatory state. Consistent with this idea, the B4GALT1 enhancer SVA insertion is linked to a genomic region implicated by GWAS in both inflammatory conditions and autoimmune diseases, such as systemic lupus erythematosus and Crohn’s disease. We explore this example and the other cases uncovered by our genome-wide screen in an effort to illuminate how retrotransposon insertion polymorphisms can impact human health and disease by causing changes in gene expression.