Characterization and genomic mapping of chimeric ERV9 endogenous retroviruses-host gene transcripts

Utility of next-generation RNA-sequencing in identifying chimeric transcription involving human endogenous retroviruses

Several studies have shown that human endogenous retroviruses and endogenous retrovirus-like repeats (here collectively HERVs) impose direct regulation on human genes through enhancer and promoter motifs present in their long terminal repeats (LTRs). Although chimeric transcription in which novel gene isoforms containing retroviral and human sequence are transcribed from viral promoters are commonly associated with disease, regulation by HERVs is beneficial in other settings; for example, in human testis chimeric isoforms of TP63 induced by an ERV9 LTR protect the male germ line upon DNA damage by inducing apoptosis, whereas in the human globin locus the γ- and β-globin switch during normal hematopoiesis is mediated by complex interactions of an ERV9 LTR and surrounding human sequence. The advent of deep sequencing or next-generation sequencing (NGS) has revolutionized the way researchers solve important scientific questions and develop novel hypotheses in relation to human genome regulation. We recently applied next-generation paired-end RNA-sequencing (RNA-seq) together with chromatin immunoprecipitation with sequencing (ChIP-seq) to examine ERV9 chimeric transcription in human reference cell lines from Encyclopedia of DNA Elements (ENCODE). This led to the discovery of advanced regulation mechanisms by ERV9s and other HERVs across numerous human loci including transcription of large gene-unannotated genomic regions, as well as cooperative regulation by multiple HERVs and non-LTR repeats such as Alu elements. In this article, well-established examples of human gene regulation by HERVs are reviewed followed by a description of paired-end RNA-seq, and its application in identifying chimeric transcription genome-widely. Based on integrative analyses of RNA-seq and ChIP-seq, data we then present novel examples of regulation by ERV9s of tumor suppressor genes CADM2 and SEMA3A, as well as transcription of an unannotated region. Taken together, this article highlights the high suitability of contemporary sequencing methods in future analyses of human biology in relation to evolutionary acquired retroviruses in the human genome.

Comprehensive identification of genes driven by ERV9-LTRs reveals TNFRSF10B as a re-activatable mediator of testicular cancer cell death

The long terminal repeat (LTR) of human endogenous retrovirus type 9 (ERV9) acts as a germline-specific promoter that induces the expression of a proapoptotic isoform of the tumor suppressor homologue p63, GTAp63, in male germline cells. Testicular cancer cells silence this promoter, but inhibitors of histone deacetylases (HDACs) restore GTAp63 expression and give rise to apoptosis. We show here that numerous additional transcripts throughout the genome are driven by related ERV9-LTRs. 3' Rapid amplification of cDNA ends (3'RACE) was combined with next-generation sequencing to establish a large set of such mRNAs. HDAC inhibitors induce these ERV9-LTR-driven genes but not the LTRs from other ERVs. In particular, a transcript encoding the death receptor DR5 originates from an ERV9-LTR inserted upstream of the protein coding regions of the TNFRSF10B gene, and it shows an expression pattern similar to GTAp63. When treating testicular cancer cells with HDAC inhibitors as well as the death ligand TNF-related apoptosis-inducing ligand (TRAIL), rapid cell death was observed, which depended on TNFRSF10B expression. HDAC inhibitors also cooperate with cisplatin (cDDP) to promote apoptosis in testicular cancer cells. ERV9-LTRs not only drive a large set of human transcripts, but a subset of them acts in a proapoptotic manner. We propose that this avoids the survival of damaged germ cells. HDAC inhibition represents a strategy of restoring the expression of a class of ERV9-LTR-mediated genes in testicular cancer cells, thereby re-enabling tumor suppression.

The role of human endogenous retroviral long terminal repeat sequences in human cancer (Review)

Human endogenous retrovirus (HERV) and solitary long terminal repeats (LTRs) constitute 8% of the human genome. Although most HERV genes are partially deleted and not intact, HERV LTRs comprise features including promoters, enhancers, selective splicer sites and polyadenylation sites in order to regulate the expression of neighboring genes. Owing to the genetic instability of LTRs, their wide distributions along human chromosomes are not only non-random, but are also correlated with gene density. Considerable evidence indicates that HERV LTRs regulate the expression of their adjacent viral and cellular genes in placental development and tumorigenesis. However, the regulatory mechanism of HERV LTRs on the expression of its neighboring cancer-associated genes in human cancers remains to be elucidated. Insertional mutagenesis, recombination and polymorphism are three principal factors of LTR that contribute to its genetic instability. Moreover, genetic instability, hypomethylation, transactivation and the antisense transcript of LTRs enhance the activity of LTRs and regulate the expression of their adjacent genes in human cancers. Therefore, in the present review, we examined the mechanism of HERV LTRs in tumorigenesis in combination with the structure and function of LTRs.

The long terminal repeat (LTR) of ERV-9 human endogenous retrovirus binds to NF-Y in the assembly of an active LTR enhancer complex NF-Y/MZF1/GATA-2

The solitary ERV-9 long terminal repeat (LTR) located upstream of the HS5 site in the human beta-globin locus control region exhibits prominent enhancer activity in embryonic and erythroid cells. The LTR enhancer contains 14 tandemly repeated subunits with recurrent CCAAT, GTGGGGA, and GATA motifs. Here we showed that in erythroid K562 cells these DNA motifs bound the following three transcription factors: ubiquitous NF-Y and hematopoietic MZF1 and GATA-2. These factors and their target DNA motifs exhibited a hierarchy of DNA/protein and protein/protein binding affinities: NF-Y/CCAAT > NF-Y/GATA-2 > NF-Y/MZF1 > MZF1/GTGGGGA; GATA-2/GATA. Through protein/protein interactions, NF-Y bound at the CCAAT motif recruited MZF1 and GATA-2, but not Sp1 and GATA-1, and stabilized their binding to the neighboring GTGGGGA and GATA sites to assemble a novel LTR enhancer complex, NF-Y/MZF1/GATA-2. In the LTR-HS5-epsilonp-GFP plasmid integrated into K562 cells, mutation of the CCAAT motif in the LTR enhancer to abolish NF-Y binding inactivated the enhancer, closed down the chromatin structure of the epsilon-globin promoter, and silenced transcription of the green fluorescent protein gene. The results indicated that NF-Y bound at the CCAAT motifs assembled a robust LTR enhancer complex, which could act over the intervening DNA to remodel the chromatin structure and to stimulate the transcription of the downstream gene locus.

The ERV-9 LTR enhancer is not blocked by the HS5 insulator and synthesizes through the HS5 site non-coding, long RNAs that regulate LTR enhancer function

A solitary long terminal repeat (LTR) of ERV-9 human endogenous retrovirus is located upstream of the HS5 site in the human beta-globin locus control region and possesses unique enhancer activity in erythroid K562 cells. In cells transfected with plasmid LTR-HS5-epsilonp-GFP, the LTR enhancer activates the GFP reporter gene and is not blocked by the interposed HS5 site, which has been reported to have insulator function. The LTR enhancer initiates synthesis of long RNAs from the LTR promoter through the intervening HS5 site into the epsilon-globin promoter and the GFP gene. Synthesis of the sense, long LTR RNAs is correlated with high level synthesis of GFP mRNA from the epsilon-globin promoter. Mutations of the LTR promoter and/or the epsilon-globin promoter show that (i) the LTR enhancer can autonomously initiate synthesis of LTR RNAs independent of the promoters and (ii) the LTR RNAs are not processed into GFP mRNA or translated into GFP. However, reversing the orientation of the LTR in plasmid (LTR)rev-HS5-epsilonp-GFP, thus reversing the direction of synthesis of LTR RNAs in the antisense direction away from the epsilon-globin promoter and GFP gene drastically reduces the level of GFP mRNA and thus LTR enhancer function. The results suggest that the LTR-assembled transcription machinery in synthesizing non-coding, LTR RNAs can reach the downstream epsilon-globin promoter to activate transcription of the GFP gene.

The solitary long terminal repeats of ERV-9 endogenous retrovirus are conserved during primate evolution and possess enhancer activities in embryonic and hematopoietic cells

The solitary long terminal repeats (LTRs) of ERV-9 endogenous retrovirus contain the U3, R, and U5 regions but no internal viral genes. They are middle repetitive DNAs present at 2,000 to 4,000 copies in primate genomes. Sequence analyses of the 5" boundary area of the erythroid beta-globin locus control region (beta-LCR) and the intron of the embryonic axin gene show that a solitary ERV-9 LTR has been stably integrated in the respective loci for at least 15 million years in the higher primates from orangutan to human. Functional studies utilizing the green fluorescent protein (GFP) gene as the reporter in transfection experiments show that the U3 region of the LTRs possesses strong enhancer activity in embryonic cells of widely different tissue origins and in adult cells of blood lineages. In both the genomic LTRs of embryonic placental cells and erythroid K562 cells and transfected LTRs of recombinant GFP plasmids in K562 cells, the U3 enhancer activates synthesis of RNAs that are initiated from a specific site 25 bases downstream of the AATAAA (TATA) motif in the U3 promoter. A second AATAAA motif in the R region does not serve as the TATA box or as the polyadenylation signal. The LTR-initiated RNAs extend through the R and U5 regions into the downstream genomic DNA. The results suggest that the ERV-9 LTR-initiated transcription process may modulate transcription of the associated gene loci in embryonic and hematopoietic cells.

Endogenous retroviruses provide the primary polyadenylation signal for two new human genes (HHLA2 and HHLA3)

By screening the expressed sequence tag (EST) database, we identified transcripts of two new human genes that are polyadenylated within a long terminal repeat (LTR) of the HERV-H endogenous retrovirus family. The first gene, termed HHLA2, is represented by two EST clones and one cDNA clone, all of which have a polyadenylated LTR as their 3' end. The gene has an open reading frame (ORF) of 414 amino acids with three immunoglobulin-like domains and is expressed primarily in intestinal tissues, kidney, and lung. Seven small EST clones from several different tissues were found for the second gene, termed HHLA3. As with HHLA2, all HHLA3 ESTs utilized a HERV-H LTR as the polyadenylation signal. Three types of alternatively spliced HHLA3 transcripts that could encode proteins of 76, 121, or 153 amino acids were detected. Interestingly, the ORF for two of these transcripts continues into the LTR. For both HHLA2 and 3, no major human transcripts that utilized a non-LTR polyadenylation signal were detected. Analysis of RNA from baboon, which lacks the LTRs at these genomic loci, showed that the baboon HHLA2 and 3 genes use other polyadenylation signals. This study demonstrates that ancient retroviral insertions have assumed gene regulatory functions during the course of human evolution.

A long terminal repeat of the human endogenous retrovirus ERV-9 is located in the 5' boundary area of the human beta-globin locus control region

Transcription of the human beta-like globin genes in erythroid cells is regulated by the far-upstream locus control region (LCR). In an attempt to define the 5' border of the LCR, we have cloned and sequenced 5 kb of new upstream DNA. We found an LTR retrotransposon belonging to the ERV-9 family of human endogenous retroviruses in the apparent 5' boundary area of the LCR. This ERV-9 LTR contains an unusual U3 enhancer region composed of 14 tandem repeats with recurrent GATA, CACCC, and CCAAT motifs. This LTR is conserved in human and gorilla, indicating its evolutionary stability in the genomes of the higher primates. In both recombinant constructs and the endogenous human genome, the LTR enhancer and promoter activate the transcription of cis-linked DNA preferentially in erythroid cells. Our findings suggest the possibility that this LTR retrotransposon may serve a relevant host function in regulating the transcription of the beta-globin LCR.

Perpetually mobile footprints of ancient infections in human genome

Up to 1% of the human genome is represented by human endogenous retroviruses (HERVs) and their fragments that are likely footprints of ancient primate germ-cell infections by retroviruses that occurred 10-60 million years ago. HERV solitary long terminal repeats (LTRs) can be often met in close vicinity to functional genes. The LTRs comprise a set of regulatory sequences like promoters, enhancers, hormone responsive elements and polyadenylation signals that might come out as new regulatory signals to resident genes and thus change their regulation in evolution. Moreover, the LTRs have a potential for chromatin remodeling that can also modulate gene expression. This review describes the integration specificity and distribution of the HERVs and LTRs in the human genome and discusses possible functional consequences of their integration in the vicinity of genes.