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

To ERV Is Human: A Phenotype-Wide Scan Linking Polymorphic Human Endogenous Retrovirus-K Insertions to Complex Phenotypes

Approximately 8% of the human genome is comprised of endogenous retroviral insertions (ERVs) originating from historic retroviral integration into germ cells. The function of ERVs as regulators of gene expression is well established. Less well studied are insertional polymorphisms of ERVs and their contribution to the heritability of complex phenotypes. The most recent integration of ERV, HERV-K, is expressed in a range of complex human conditions from cancer to neurologic diseases. Using an in-house computational pipeline and whole-genome sequencing data from the diverse 1,000 Genomes Phase 3 population (n = 2,504), we identified 46 polymorphic HERV-K insertions that are tagged by adjacent single nucleotide polymorphisms (SNPs). To test the potential role of polymorphic HERV-K in the heritability of complex diseases, existing databases were queried for enrichment of established relationships between the HERV-K insertion-associated SNPs (hiSNPs), and tissue specific gene expression and disease phenotypes. Overall, hiSNPs for the 46 polymorphic HERV-K sites were statistically enriched (p < 1.0E^-16) for eQTLs across 44 human tissues. Fifteen of the 46 HERV-K insertions had hiSNPs annotated in the EMBL-EBI GWAS Catalog and cumulatively associated with >100 phenotypes. Experimental factor ontology enrichment analysis suggests that polymorphic HERV-K specifically contribute to neurologic and immunologic disease phenotypes, including traits related to intra cranial volume (FDR 2.00E-09), Parkinson's disease (FDR 1.80E-09), and autoimmune diseases (FDR 1.80E-09). These results provide strong candidates for context-specific study of polymorphic HERV-K insertions in disease-related traits, serving as a roadmap for future studies of the heritability of complex disease.

Hypermethylated LTR retrotransposon exhibits enhancer activity

LTR retrotransposons are repetitive DNA elements comprising ∼10% of the human genome. They are silenced by hypermethylation of cytosines in CpG dinucleotides and are considered parasitic DNA serving no useful function for the host genome. However, hypermethylated LTRs contain enhancer and promoter sequences and can promote tissue-specific transcription of cis-linked genes. To resolve the apparent paradox of hypermethylated LTRs possessing transcriptional activities, we studied the ERV-9 LTR retrotransposon located at the 5' border of the transcriptionally active β-globin gene locus in human erythroid progenitor and erythroleukemia K562 cells. We found that the ERV-9 LTR, containing 65 CpGs in 1.7 kb DNA, was hypermethylated (with > 90% methylated CpGs). Hypermethylated LTR possessed transcriptional enhancer activity, since in vivo deletion of the LTR by CRISPR-cas9 suppressed transcription of the globin genes by > 50%. ChIP-qPCR and ChIP-seq studies showed that the hypermethylated LTR enhancer spanning recurrent CCAATCG and GATA motifs associated respectively with key transcription factors (TFs) NF-Y and GATA-1 and -2 at reduced levels, compared with the unmethylated LTR in transfected LTR-reporter gene plasmids. Electrophoretic mobility shift assays with methylated LTR enhancer probe showed that the methylated probe bound both NF-Y and GATA-1 and -2 with lower affinities than the unmethylated enhancer probe. Thus, hypermethylation drastically reduced, but did not totally abolish, the binding affinities of the enhancer motifs to the key TFs to assemble the LTR-pol II transcription complex that activated transcription of cis-linked genes at reduced efficiency.

New bioinformatic tool for quick identification of functionally relevant endogenous retroviral inserts in human genome

Endogenous retroviruses (ERVs) and LTR retrotransposons (LRs) occupy ∼8% of human genome. Deep sequencing technologies provide clues to understanding of functional relevance of individual ERVs/LRs by enabling direct identification of transcription factor binding sites (TFBS) and other landmarks of functional genomic elements. Here, we performed the genome-wide identification of human ERVs/LRs containing TFBS according to the ENCODE project. We created the first interactive ERV/LRs database that groups the individual inserts according to their familial nomenclature, number of mapped TFBS and divergence from their consensus sequence. Information on any particular element can be easily extracted by the user. We also created a genome browser tool, which enables quick mapping of any ERV/LR insert according to genomic coordinates, known human genes and TFBS. These tools can be used to easily explore functionally relevant individual ERV/LRs, and for studying their impact on the regulation of human genes. Overall, we identified ∼110,000 ERV/LR genomic elements having TFBS. We propose a hypothesis of “domestication” of ERV/LR TFBS by the genome milieu including subsequent stages of initial epigenetic repression, partial functional release, and further mutation-driven reshaping of TFBS in tight coevolution with the enclosing genomic loci.

Molecular functions of human endogenous retroviruses in health and disease

Human endogenous retroviruses (HERVs) and related genetic elements form 504 distinct families and occupy ~8% of human genome. Recent success of high-throughput experimental technologies facilitated understanding functional impact of HERVs for molecular machinery of human cells. HERVs encode active retroviral proteins, which may exert important physiological functions in the body, but also may be involved in the progression of cancer and numerous human autoimmune, neurological and infectious diseases. The spectrum of related malignancies includes, but not limits to, multiple sclerosis, psoriasis, lupus, schizophrenia, multiple cancer types and HIV. In addition, HERVs regulate expression of the neighboring host genes and modify genomic regulatory landscape, e.g., by providing regulatory modules like transcription factor binding sites (TFBS). Indeed, recent bioinformatic profiling identified ~110,000 regulatory active HERV elements, which formed at least ~320,000 human TFBS. These and other peculiarities of HERVs might have played an important role in human evolution and speciation. In this paper, we focus on the current progress in understanding of normal and pathological molecular niches of HERVs, on their implications in human evolution, normal physiology and disease. We also review the available databases dealing with various aspects of HERV genetics.

VlincRNAs controlled by retroviral elements are a hallmark of pluripotency and cancer

BACKGROUND

The function of the non-coding portion of the human genome remains one of the most important questions of our time. Its vast complexity is exemplified by the recent identification of an unusual and notable component of the transcriptome - very long intergenic non-coding RNAs, termed vlincRNAs.

RESULTS

Here we identify 2,147 vlincRNAs covering 10 percent of our genome. We show they are present not only in cancerous cells, but also in primary cells and normal human tissues, and are controlled by canonical promoters. Furthermore, vlincRNA promoters frequently originate from within endogenous retroviral sequences. Strikingly, the number of vlincRNAs expressed from endogenous retroviral promoters strongly correlates with pluripotency or the degree of malignant transformation. These results suggest a previously unknown connection between the pluripotent state and cancer via retroviral repeat-driven expression of vlincRNAs. Finally, we show that vlincRNAs can be syntenically conserved in humans and mouse and their depletion using RNAi can cause apoptosis in cancerous cells.

CONCLUSIONS

These intriguing observations suggest that vlincRNAs could create a framework that combines many existing short ESTs and lincRNAs into a landscape of very long transcripts functioning in the regulation of gene expression in the nucleus. Certain types of vlincRNAs participate at specific stages of normal development and, based on analysis of a limited set of cancerous and primary cell lines, they appear to be co-opted by cancer-associated transcriptional programs. This provides additional understanding of transcriptome regulation during the malignant state, and could lead to additional targets and options for its reversal.

Genes associated with the cis-regulatory functions of intragenic LINE-1 elements

BACKGROUND

Thousands of intragenic long interspersed element 1 sequences (LINE-1 elements or L1s) reside within genes. These intragenic L1 sequences are conserved and regulate the expression of their host genes. When L1 methylation is decreased, either through chemical induction or in cancer, the intragenic L1 transcription is increased. The resulting L1 mRNAs form RISC complexes with pre-mRNA to degrade the complementary mRNA. In this study, we screened for genes that are involved in intragenic L1 regulation networks.

RESULTS

Genes containing L1s were obtained from L1Base (http://l1base.molgen.mpg.de). The expression profiles of 205 genes in 516 gene knockdown experiments were obtained from the Gene Expression Omnibus (GEO) (http://www.ncbi.nlm.nih.gov/geo). The expression levels of the genes with and without L1s were compared using Pearson's chi-squared test. After a permutation based statistical analysis and a multiple hypothesis testing, 73 genes were found to induce significant regulatory changes (upregulation and/or downregulation) in genes with L1s. In detail, 5 genes were found to induce both the upregulation and downregulation of genes with L1s, whereas 27 and 37 genes induced the downregulation and upregulation, respectively, of genes with L1s. These regulations sometimes differed depending on the cell type and the orientation of the intragenic L1s. Moreover, the siRNA-regulating genes containing L1s possess a variety of molecular functions, are responsible for many cellular phenotypes and are associated with a number of diseases.

CONCLUSIONS

Cells use intragenic L1s as cis-regulatory elements within gene bodies to modulate gene expression. There may be several mechanisms by which L1s mediate gene expression. Intragenic L1s may be involved in the regulation of several biological processes, including DNA damage and repair, inflammation, immune function, embryogenesis, cell differentiation, cellular response to external stimuli and hormonal responses. Furthermore, in addition to cancer, intragenic L1s may alter gene expression in a variety of diseases and abnormalities.

Functional and structural divergence of an unusual LTR retrotransposon family in plants

Retrotransposons with long terminal repeats (LTRs) more than 3 kb are not frequent in most eukaryotic genomes. Rice LTR retrotransposon, Retrosat2, has LTRs greater than 3.2 kb and two open reading frames (ORF): ORF1 encodes enzymes for retrotransposition whereas no function can be assigned to ORF0 as it is not found in any other organism. A variety of experimental and in silico approaches were used to determine the origin of Retrosat2 and putative function of ORF0. Our data show that not only is Retrosat2 highly abundant in the Oryza genus, it may yet be active in rice. Homologs of Retrosat2 were identified in maize, sorghum, Arabidopsis and other plant genomes suggesting that the Retrosat2 family is of ancient origin. Several putatively cis-acting elements, some multicopy, that regulate retrotransposon replication or responsiveness to environmental factors were found in the LTRs of Retrosat2. Unlike the ORF1, the ORF0 sequences from Retrosat2 and homologs are divergent at the sequence level, 3D-structures and predicted biological functions. In contrast to other retrotransposon families, Retrosat2 and its homologs are dispersed throughout genomes and not concentrated in the specific chromosomal regions, such as centromeres. The genomic distribution of Retrosat2 homologs varies across species which likely reflects the differing evolutionary trajectories of this retrotransposon family across diverse species.

Smoking increases transcription of human endogenous retroviruses in a newly established in vitro cell model and in normal urothelium

Human endogenous retroviruses (HERVs) accounting for 9% of the human genome are considered as surrogate markers for genetic instability and as a driving force of genetic variation. Moreover, they modulate regular gene activities and give rise to expression of disease-associated peptides that may serve as diagnostic markers or even targets for T cell-based immune responses. To date, no data are available on the potential link between urothelial carcinogenesis, HERV activity, and tobacco smoking, the main risk for bladder cancer. Here, we report on potential alterations in HERV transcription induced by smoking in a newly established in vitro model and in human urothelium. Normal human dermal fibroblasts were cultivated with urine from never (n = 6) and current smokers (n = 6) and transcription levels for the HERV subfamilies HERV-E 4-1, HERV-T S71-TK1, and HERV-K HML-6 were measured by quantitative real-time PCR. Tendencies toward increased mean transcript levels were detected for cells treated with urine from current smokers. Equally, activity measured in human urothelium supported an increase of HERV transcription in current smokers (n = 9) compared to never smokers (n = 4).

Inactive X chromosome-specific histone H3 modifications and CpG hypomethylation flank a chromatin boundary between an X-inactivated and an escape gene

In mammals, the dosage compensation of sex chromosomes between males and females is achieved by transcriptional inactivation of one of the two X chromosomes in females. However, a number of genes escape X-inactivation in humans. It remains poorly understood how the transcriptional activity of these ‘escape genes’ is maintained despite the chromosome-wide heterochromatin formation. To address this question, we analyzed a putative chromatin boundary between the inactivated RBM10 and an escape gene, UBA1/UBE1. Chromatin immunoprecipitation revealed that trimethylated histone H3 lysine 9 and H4 lysine 20 were enriched in the last exon through the proximal downstream region of RBM10, but were remarkably diminished at ∼2 kb upstream of the UBA1 transcription start site. Whereas RNA polymerase II was not loaded onto the intergenic region, CTCF (CCCTC binding factor) was enriched around the boundary, where some CpG sites were hypomethylated specifically on inactive X. These findings suggest that local DNA hypomethylation and CTCF binding are involved in the formation of a chromatin boundary, which protects the UBA1 escape gene against the chromosome-wide transcriptional silencing.

Unique functions of repetitive transcriptomes

Repetitive sequences occupy a huge fraction of essentially every eukaryotic genome. Repetitive sequences cover more than 50% of mammalian genomic DNAs, whereas gene exons and protein-coding sequences occupy only ~3% and 1%, respectively. Numerous genomic repeats include genes themselves. They generally encode "selfish" proteins necessary for the proliferation of transposable elements (TEs) in the host genome. The major part of evolutionary "older" TEs accumulated mutations over time and fails to encode functional proteins. However, repeats have important functions also on the RNA level. Repetitive transcripts may serve as multifunctional RNAs by participating in the antisense regulation of gene activity and by competing with the host-encoded transcripts for cellular factors. In addition, genomic repeats include regulatory sequences like promoters, enhancers, splice sites, polyadenylation signals, and insulators, which actively reshape cellular transcriptomes. TE expression is tightly controlled by the host cells, and some mechanisms of this regulation were recently decoded. Finally, capacity of TEs to proliferate in the host genome led to the development of multiple biotechnological applications.

Retroelements and their impact on genome evolution and functioning

Retroelements comprise a considerable fraction of eukaryotic genomes. Since their initial discovery by Barbara McClintock in maize DNA, retroelements have been found in genomes of almost all organisms. First considered as a "junk DNA" or genomic parasites, they were shown to influence genome functioning and to promote genetic innovations. For this reason, they were suggested as an important creative force in the genome evolution and adaptation of an organism to altered environmental conditions. In this review, we summarize the up-to-date knowledge of different ways of retroelement involvement in structural and functional evolution of genes and genomes, as well as the mechanisms generated by cells to control their retrotransposition.

A facilitated tracking and transcription mechanism of long-range enhancer function

In the human ε−globin gene locus, the HS2 enhancer in the Locus Control Region regulates transcription of the embryonic ε-globin gene located over 10 kb away. The mechanism of long-range HS2 enhancer function was not fully established. Here we show that the HS2 enhancer complex containing the enhancer DNA together with RNA polymerase II (pol II) and TBP tracks along the intervening DNA, synthesizing short, polyadenylated, intergenic RNAs to ultimately loop with the ε-globin promoter. Guided by this facilitated tracking and transcription mechanism, the HS2 enhancer delivers pol II and TBP to the cis-linked globin promoter to activate mRNA synthesis from the target gene. An insulator inserted in the intervening DNA between the enhancer and the promoter traps the enhancer DNA and the associated pol II and TBP at the insulator site, blocking mid-stream the facilitated tracking and transcription mechanism of the enhancer complex, thereby blocking long-range enhancer function.

Genotypic heterogeneity and correlation to intergenic haplotype within high HbF beta-thalassemia intermedia

OBJECTIVES

A molecular study was carried out of beta-thalassemia intermedia patients, compound heterozygotes for mutations usually found in beta-thalassemia major, with high levels of HbF in the absence of hereditary persistence of fetal hemoglobin (HPFH) syndrome. Our objective was to locate cis-DNA structures, DNA haplotypes, motifs, or polymorphisms that may correlate with the presence of high HbF.

METHODS

Allele-specific oligonucleotide (ASO) hybridization was used for the detection of mutations and restriction fragment length polymorphism (RFLP) analysis and automated sequencing for motifs, haplotypes, and polymorphisms. Southern blot was used for investigating alpha-thalassemia and/or alpha- or gamma-globin genes triplications. RNA extracted from burst forming unit-erythroid (BFU-e) colonies of peripheral blood mononuclear cell cultures was used in reverse transcriptase-polymerase chain reaction (RT-PCR) to investigate intergenic transcription.

RESULTS

We established that (i) the combination: T haplotype of the Agamma-delta-globin intergenic region, the motif (TA)9N10(TA)10 in the HS2 site of locus control region (LCR), and TAG pre-Ggamma haplotype is sufficient but not necessary for high HbF, (ii) the genetic determinant(s) for high HbF involves an element associated with this combination and must be present in the specific R haplotype occurring in beta-thalassemia intermedia and (iii) the genetic determinant(s) for high HbF does not involve the abolition of intergenic transcription in the Agamma-delta-globin intergenic region.

CONCLUSIONS

The genetic determinant(s) of high HbF in the absence of HPFH is linked to intergenic haplotype T and does not disrupt intergenic transcription.

The HS2 enhancer of the beta-globin locus control region initiates synthesis of non-coding, polyadenylated RNAs independent of a cis-linked globin promoter

The HS2 enhancer in the beta-globin locus control region (LCR) regulates transcription of the globin genes 10-50 kb away. Earlier studies show that a transcription mechanism initiated by the HS2 enhancer through the intervening DNA in the direction of the cis-linked promoter and gene mediates long-range enhancer function. Here, we further analyzed the enhancer-initiated RNAs and their mode of transcription from the HS2 enhancer in the endogenous genome of erythroid K562 cells, in plasmids integrated into K562 cells and in purified DNA used as template in in vitro transcription reactions. We found that the HS2 enhancer was able to initiate transcription autonomously in the absence of a cis-linked globin promoter. The enhancer-initiated, intergenic RNAs were different from the mRNA synthesized at the promoter in several aspects. The enhancer RNAs were synthesized not from a defined site but from multiple sites both within and as far as 1 kb downstream of the enhancer. The enhancer RNAs did not appear to contain a normal cap structure at the 5' ends. They were polyadenylated at multiple sites within 3 kb downstream of their initiation sites and were therefore shorter than 3 kb in lengths. The enhancer RNAs remained in discrete spots within the nucleus and were not processed into mRNA or translated into proteins. These particular features of enhancer-initiated transcription indicate that the transcriptional complex assembled by the enhancer was different from the basal transcription complex assembled at the promoter. The results suggest that in synthesizing non-coding, intergenic RNAs, the enhancer-assembled transcription complex could track through the intervening DNA to reach the basal promoter complex and activate efficient mRNA synthesis from the promoter.

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.

HS2 enhancer function is blocked by a transcriptional terminator inserted between the enhancer and the promoter

The HS2 enhancer in the beta-globin locus control region regulates transcription of the globin genes 10-50 kb away. How the HS2 enhancer acts over this distance is not clearly understood. Earlier studies show that in erythroid cells the HS2 enhancer initiates synthesis of intergenic RNAs from sites within and downstream of the enhancer, and the enhancer-initiated RNAs are transcribed through the intervening DNA into the cis-linked promoter and gene. To investigate the functional significance of the enhancer-initiated transcription, here we inserted the lac operator sequence in the intervening DNA between the HS2 enhancer and the epsilon-globin promoter in reporter plasmids and integrated the plasmids into erythroid K562 cells expressing the lac repressor protein. We found that the interposed lac operator/repressor complex blocked the elongation of enhancer-initiated transcription through the intervening DNA and drastically reduced HS2 enhancer function as measured by the level of mRNA synthesized from the epsilon-globin promoter. The results indicate that the tracking and transcription mechanism of the HS2 enhancer-assembled transcriptional machinery from the enhancer through the intervening DNA into the cis-linked promoter can mediate enhancer-promoter interaction over a long distance.