Vanadium Induces VL30 Retrotransposition at an Unusually High Level: A Possible Carcinogenesis Mechanism

A Systems Biology Approach on the Regulatory Footprint of Human Endogenous Retroviruses (HERVs)

Human endogenous retroviruses (HERVs) are a family of endogenous retroviruses that comprise the ~8.93% of the human genome sequence, with a high proportion being human specific. The recent expansion of repeated HERV sequences has offered a framework for genetic and epigenetic innovation. In the current report, a systematic approach is implemented to catalogue regulatory elements within HERVs, as a roadmap to potential functions of HERV sequences in gene networks. ENCODE Project has offered a wealth of epigenetic data based on omics technologies. I analyzed the presence of HERV sequences on consensus cis-regulatory elements (cCREs) from ENCODE data. On the one side, HERVs are in 1 out of 9 cCREs (>100.000 cCREs in total), dispersed within the genome and present in cis-regulatory regions of ~81% of human genes, as calculated following gene enrichment analysis. On the other side, promoter-associated HERV cCREs are present adjacent to (in a 200 bp window) the transcription start sites of 256 human genes. Regulatory network production, followed by centrality analysis led to the discovery of 90 core genes containing HERV-associated promoters. Pathway analysis on the core network genes and their immediate neighbors revealed a regulatory footprint that, among others, is associated with inflammation, chemokine signaling and response to viral infection. Collectively, these results support the concept that the expansion of regulatory sequences derived from HERVs is critical for epigenetic innovation that may have wired together genes into novel transcriptional networks with critical roles in cellular physiology and pathology.

Origins and Function of VL30 lncRNA Packaging in Small Extracellular Vesicles: Implications for Cellular Physiology and Pathology

Long non-coding RNAs (lncRNAs) have emerged during the post-genomic era as significant epigenetic regulators. Viral-like 30 elements (VL30s) are a family of mouse retrotransposons that are transcribed into functional lncRNAs. Recent data suggest that VL30 RNAs are efficiently packaged in small extracellular vesicles (SEVs) through an SEV enrichment sequence. We analysed VL30 elements for the presence of the distinct 26 nt SEV enrichment motif and found that SEV enrichment is an inherent hallmark of the VL30 family, contained in 36 full-length elements, with a widespread chromosomal distribution. Among them, 25 elements represent active, present-day integrations and contain an abundance of regulatory sequences. Phylogenetic analysis revealed a recent spread of SEV-VL30s from 4.4 million years ago till today. Importantly, 39 elements contain an SFPQ-binding motif, associated with the transcriptional induction of oncogenes. Most SEV-VL30s reside in transcriptionally active regions, as characterised by their distribution adjacent to candidate cis-regulatory elements (cCREs). Network analysis of SEV-VL30-associated genes suggests a distinct transcriptional footprint associated with embryonal abnormalities and neoplasia. Given the established role of VL30s in oncogenesis, we conclude that their potential to spread through SEVs represents a novel mechanism for non-coding RNA biology with numerous implications for cellular homeostasis and disease.

Tinzaparin inhibits VL30 retrotransposition induced by oxidative stress and/or VEGF in HC11 mouse progenitor mammary cells: Association between inhibition of cancer stem cell proliferation and mammosphere disaggregation

Tinzaparin is an anticoagulant and antiangiogenic drug with inhibitory properties against tumor growth. VEGF stimulates angiogenesis, while an association between reactive oxygen species (ROS) and angiogenesis is involved in tumor progression. The present study aimed to investigate the effect of tinzaparin on VL30 retrotransposition-positive mouse HC11 mammary stem-like epithelial cells, previously reported to be associated with induced mammosphere/cancer stem cell (CSC) generation and tumorigenesis. Under 24 h serum starvation, 15.2% nominal retrotransposition frequency was increased to 29%. Additionally, while treatment with 3–12 ng/ml VEGF further induced retrotransposition frequency in a dose-dependent manner (up to 40.3%), pre-incubation with tinzaparin (2 IU/ml) for 0.5–4 h reduced this frequency to 18.3% in a time-dependent manner, confirmed by analogous results in NIH3T3 fibroblasts. Treatment with 10–40 pg/ml glucose oxidase (GO) for 24 h induced HC11 cell retrotransposition in a dose-dependent manner (up to 82.5%), while a 3 h pre-incubation with tinzaparin (1 or 2 IU/ml) elicited a 13.5 or 25.5% reduction in retrotransposition, respectively. Regarding tumorigenic VL30 retrotransposition-positive HC11 cells, treatment with 2 IU/ml tinzaparin for 5 days reduced proliferation rate in a time-dependent manner (up to ~55%), and after 3 weeks, disaggregated soft agar-formed foci, as well as low-adherent mammospheres, producing single mesenchymal-like cells with a ~50% reduced retrotransposition. With respect to the VL30 retrotransposition mechanism: While 12 ng/ml VEGF increased the level of VL30 and endogenous reverse transcriptase (enRT) transcripts ~1.41- and ~1.16-fold, respectively, subsequent tinzaparin treatment reduced both endogenous/ROS- and VEGF-induced levels 1.15- and 0.40-fold (VL30) and 0.60- and 0.52-fold (enRT), respectively. To the best of our knowledge, these data demonstrate for the first time, the novel inhibition activity of tinzaparin against ROS- and VEGF-induced VL30 retrotransposition, and the proliferation and/or aggregation of mouse HC11 mammosphere/tumor-initiating CSCs, thus contributing to the inhibition of VL30 retrotransposition-induced primary tumor growth.

VL30 retrotransposition is associated with induced EMT, CSC generation and tumorigenesis in HC11 mouse mammary stem‑like epithelial cells

Retrotransposons copy their sequences via an RNA intermediate, followed by reverse transcription into cDNA and random insertion, into a new genomic locus. New retrotransposon copies may lead to cell transformation and/or tumorigenesis through insertional mutagenesis. Methylation is a major defense mechanism against retrotransposon RNA expression and retrotransposition in differentiated cells, whereas stem cells are relatively hypo-methylated. Epithelial-to-mesenchymal transition (EMT), which transforms normal epithelial cells into mesenchymal-like cells, also contributes to tumor progression and tumor metastasis. Cancer stem cells (CSCs), a fraction of undifferentiated tumor-initiating cancer cells, are reciprocally related to EMT. In the present study, the outcome of long terminal repeat (LTR)-Viral-Like 30 (VL30) retrotransposition was examined in mouse mammary stem-like/progenitor HC11 epithelial cells. The transfection of HC11 cells with a VL30 retrotransposon, engineered with an EGFP-based retrotransposition cassette, elicited a higher retrotransposition frequency in comparison to differentiated J3B1A and C127 mouse mammary cells. Fluorescence microscopy and PCR analysis confirmed the specificity of retrotransposition events. The differentiated retrotransposition-positive cells retained their epithelial morphology, while the respective HC11 cells acquired mesenchymal features associated with the loss of E-cadherin, the induction of N-cadherin, and fibronectin and vimentin protein expression, as well as an increased transforming growth factor (TGF)-β1, Slug, Snail-1 and Twist mRNA expression. In addition, they were characterized by cell proliferation in low serum, and the acquisition of CSC-like properties indicated by mammosphere formation under anchorage-independent conditions. Mammospheres exhibited an increased Nanog and Oct4 mRNA expression and a CD44⁺/CD24^−/low antigenic phenotype, as well as self-renewal and differentiation capacity, forming mammary acini-like structures. DNA sequencing analysis of retrotransposition-positive HC11 cells revealed retrotransposed VL30 copies integrated at the vicinity of EMT-, cancer type- and breast cancer-related genes. The inoculation of these cells into Balb/c mice produced cytokeratin-positive tumors containing pancytokeratin-positive cells, indicative of cell invasion features. On the whole, the findings of the present study demonstrate, for the first time, to the best of our knowledge, that stem-like epithelial HC11 cells are amenable to VL30 retrotransposition associated with the induction of EMT and CSC generation, leading to tumorigenesis.

Genomic analysis of mouse VL30 retrotransposons

BACKGROUND

Retrotransposons are mobile elements that have a high impact on shaping the mammalian genomes. Since the availability of whole genomes, genomic analyses have provided novel insights into retrotransposon biology. However, many retrotransposon families and their possible genomic impact have not yet been analysed.

RESULTS

Here, we analysed the structural features, the genomic distribution and the evolutionary history of mouse VL30 LTR-retrotransposons. In total, we identified 372 VL30 sequences categorized as 86 full-length and 49 truncated copies as well as 237 solo LTRs, with non-random chromosomal distribution. Full-length VL30s were highly conserved elements with intact retroviral replication signals, but with no protein-coding capacity. Analysis of LTRs revealed a high number of common transcription factor binding sites, possibly explaining the known inducible and tissue-specific expression of individual elements. The overwhelming majority of full-length and truncated elements (82/86 and 40/49, respectively) contained one or two specific motifs required for binding of the VL30 RNA to the poly-pyrimidine tract-binding protein-associated splicing factor (PSF). Phylogenetic analysis revealed three VL30 groups with the oldest emerging ~17.5 Myrs ago, while the other two were characterized mostly by new genomic integrations. Most VL30 sequences were found integrated either near, adjacent or inside transcription start sites, or into introns or at the 3' end of genes. In addition, a significant number of VL30s were found near Krueppel-associated box (KRAB) genes functioning as potent transcriptional repressors.

CONCLUSION

Collectively, our study provides data on VL30s related to their: (a) number and structural features involved in their transcription that play a role in steroidogenesis and oncogenesis; (b) evolutionary history and potential for retrotransposition; and (c) unique genomic distribution and impact on gene expression.

Altering Genomic Integrity: Heavy Metal Exposure Promotes Transposable Element-Mediated Damage

Maintenance of genomic integrity is critical for cellular homeostasis and survival. The active transposable elements (TEs) composed primarily of three mobile element lineages LINE-1, Alu, and SVA comprise approximately 30% of the mass of the human genome. For the past 2 decades, studies have shown that TEs significantly contribute to genetic instability and that TE-caused damages are associated with genetic diseases and cancer. Different environmental exposures, including several heavy metals, influence how TEs interact with its host genome increasing their negative impact. This mini-review provides some basic knowledge on TEs, their contribution to disease, and an overview of the current knowledge on how heavy metals influence TE-mediated damage.

Arsenic induces VL30 retrotransposition: the involvement of oxidative stress and heat-shock protein 70

Arsenic is an environmental contaminant with known cytotoxic and carcinogenic properties, but the cellular mechanisms of its action are not fully known. As retrotransposition consists a potent mutagenic factor affecting genome stability, we investigated the effect of arsenic on retrotransposition of an enhanced green fluorescent protein (EGFP)-tagged nonautonomous long terminal repeat (LTR)-retrotransposon viral-like 30 (VL30) in a mouse NIH3T3 cell culture-retrotransposition assay. Flow cytometry analysis of assay cells treated with 2.5-20μM sodium arsenite revealed induction of retrotransposition events in a dose- and time-dependent manner, which was further confirmed as genomic integrations by PCR analysis and appearance of EGFP-positive cells by UV microscopy. Specifically, 20μM sodium arsenite strongly induced the VL30 retrotransposition frequency, which was ~90,000-fold higher than the natural one and also VL30 RNA expression was ~6.6-fold. Inhibition of the activity of endogenous reverse transcriptases by efavirenz at 15μM or nevirapine at 375μM suppressed the arsenite-induced VL30 retrotransposition by 71.16 or 79.88%, respectively. In addition, the antioxidant N-acetyl-cysteine reduced the level of arsenite-induced retrotransposition, which correlated with the rescue of arsenite-induced G2/M cell cycle arrest and cell toxicity. Treatment of assay cells ectopically overexpressing the human heat-shock protein 70 (Hsp70) with 15μM sodium arsenite resulted in an additional ~4.5-fold induction of retrotransposition compared with normal assay cells, whereas treatment with 20μM produced a massive cell death. Our results show for the first time that arsenic both as an oxidative and heat-shock mimicking agent is a potent inducer of VL30 retrotransposition in mouse cells. The impact of arsenic-induced retrotransposition, as a cellular response, on contribution to or explanation of the arsenic-associated toxicity and carcinogenicity is discussed.

H2O2 signals via iron induction of VL30 retrotransposition correlated with cytotoxicity

The impact of oxidative stress on mobilization of endogenous retroviruses and their effects on cell fate is unknown. We investigated the action of H2O2 on retrotransposition of an EGFP-tagged mouse LTR-retrotransposon, VL30, in an NIH3T3 cell-retrotransposition assay. H2O2 treatment of assay cells caused specific retrotranspositions documented by UV microscopy and PCR analysis. Flow cytometric analysis revealed an unusually high dose- and time-dependent retrotransposition frequency induced, ∼420,000-fold at 40 μM H2O2 compared to the natural frequency, which was reduced by ectopic expression of catalase. Remarkably, H2O2 moderately induced the RNA expression of retrotransposon B2 without affecting the basal expression of VL30s and L1 and significantly induced the expression of various endogenous reverse transcriptase genes. Further, whereas treatment with 50 μM FeCl2 alone was ineffective, cotreatment with 10 μM H2O2 and 50 μM FeCl2 caused a 6-fold higher retrotransposition induction than H2O2 alone, which was associated with cytotoxicity. H2O2- or H2O2/FeCl2-induced retrotransposition was significantly reduced by the iron chelator DFO or the antioxidant NAC, respectively. Furthermore, both H2O2-induced retrotransposition and associated cytotoxicity were inhibited after pretreatment of cells with DFO or the reverse transcriptase inhibitors efavirenz and etravirine. Our data show for the first time that H2O2, acting via iron, is a potent stimulus of retrotransposition contributing to oxidative stress-induced cell damage.

DNA damage induction in human cells exposed to vanadium oxides in vitro

Vanadium and vanadium salts cause genotoxicity and elicit variable biological effects depending on several factors. In the present study, we analyzed and compared the DNA damage and repair processes induced by vanadium in three oxidation states. We used human blood leukocytes in vitro and in a single cell gel electrophoresis assay at two pH values. We observed that vanadium(III) trioxide and vanadium(V) pentoxide produced DNA single-strand breaks at all of the concentrations (1, 2, 4, or 8 μg/ml) and treatment times (2, 4, or 6 h) tested. Vanadium(IV) tetraoxide treatment significantly increased DNA damage at all concentrations for 4 or 6 h of treatment but not for 2 h of treatment. The DNA repair kinetics indicated that most of the cells exposed to vanadium III and V for 4 h recovered within the repair incubation time of 90 min; however, those exposed to vanadium(IV) repaired their DNA within 120 min. The data at pH 9 indicated that vanadium(IV) tetraoxide induced DNA double-strand breaks. Our results show that the genotoxic effect of vanadium can be produced by any of its three oxidation states. However, vanadium(IV) induces double-strand breaks, and it is known that these lesions are linked with forming structural chromosomal aberrations.

Epigenetic regulation of a murine retrotransposon by a dual histone modification mark

Large fractions of eukaryotic genomes contain repetitive sequences of which the vast majority is derived from transposable elements (TEs). In order to inactivate those potentially harmful elements, host organisms silence TEs via methylation of transposon DNA and packaging into chromatin associated with repressive histone marks. The contribution of individual histone modifications in this process is not completely resolved. Therefore, we aimed to define the role of reversible histone acetylation, a modification commonly associated with transcriptional activity, in transcriptional regulation of murine TEs. We surveyed histone acetylation patterns and expression levels of ten different murine TEs in mouse fibroblasts with altered histone acetylation levels, which was achieved via chemical HDAC inhibition with trichostatin A (TSA), or genetic inactivation of the major deacetylase HDAC1. We found that one LTR retrotransposon family encompassing virus-like 30S elements (VL30) showed significant histone H3 hyperacetylation and strong transcriptional activation in response to TSA treatment. Analysis of VL30 transcripts revealed that increased VL30 transcription is due to enhanced expression of a limited number of genomic elements, with one locus being particularly responsive to HDAC inhibition. Importantly, transcriptional induction of VL30 was entirely dependent on the activation of MAP kinase pathways, resulting in serine 10 phosphorylation at histone H3. Stimulation of MAP kinase cascades together with HDAC inhibition led to simultaneous phosphorylation and acetylation (phosphoacetylation) of histone H3 at the VL30 regulatory region. The presence of the phosphoacetylation mark at VL30 LTRs was linked with full transcriptional activation of the mobile element. Our data indicate that the activity of different TEs is controlled by distinct chromatin modifications. We show that activation of a specific mobile element is linked to a dual epigenetic mark and propose a model whereby phosphoacetylation of histone H3 is crucial for full transcriptional activation of VL30 elements.

The majority of genomic sequences in higher eukaryotes do not contain protein coding genes. Large fractions are covered by repetitive sequences, many of which are derived from transposable elements (TEs). These selfish genes, only containing sequences necessary for self-propagation, can multiply and change their location within the genome, threatening host genome integrity and provoking mutational bursts. Therefore host organisms have evolved a diverse repertoire of defence mechanisms to counteract and silence these genomic parasites. One way is to package DNA sequences containing TEs into transcriptionally inert heterochromatin, which is partly achieved via chemical modification of the packaging proteins associated with DNA, the histones. To better understand the contribution of histone acetylation in the activation of TEs, we treated mouse fibroblasts with a specific histone deacetylase inhibitor. By monitoring the expression of ten different types of murine mobile elements, we identified a defined subset of VL30 transposons specifically reactivated upon increased histone acetylation. Importantly, phosphorylation of histone H3, a modification that is triggered by stress, is required for acetylation-dependent activation of VL30 elements. We present a model where concomitant histone phosphorylation and acetylation cooperate in the transcriptional induction of VL30 elements.

Retrotransposon RNA expression and evidence for retrotransposition events in human oocytes

Although human diseases of retrotransposition-derived etiology have been documented, retrotransposon RNA expression and the occurrence of retrotransposition events in the human oocyte are not studied. We investigated the RNA expression of L1 and HERV-K10 retrotransposons in human oocytes by RT-PCR analysis with designed primers. Using denucleated germinal vesicles (GVs), we detected RT-PCR products of expressed L1, HERV-K10 and, unexpectedly, SINE-R, VNTR and Alu (SVA) retrotransposons. Their transcript specificities were identified as such following RNA-FISH and their origin by cloning and sequence alignment analyses. Assessing the expression level in comparison with somatic cells by densitometry analysis, we found that although in normal lymphocytes and transformed HeLa cells their profile was in an order of L1 > HERV-K10 > SVA, remarkably this was reversed in oocytes. To investigate whether de novo retrotransposition events occur and reverse transcriptases are expressed in the human oocyte, we introduced in GVs either a retrotransposition active human L1 or mouse reverse transcriptase deficient-VL30 retrotransposon tagged with an EGFP-based retrotransposition cassette. Interestingly, in both the cases, we observed EGFP-positive oocytes, associated with an abnormal morphology for L1 and granulation for VL30, and the retrotransposition events were confirmed by PCR. Our results: (i) show that L1, HERV-K10 and SVA retrotransposons are transcriptionally expressed and (ii) provide evidence, for the first time, for retrotransposition events occurring in the human oocyte. These findings suggest that both, network of retrotransposon transcripts and controlled retrotranspositions, might serve important functions required for oocyte development and fertilization while the uncontrolled ones might explain the onset of genetic disorders.