Origin of a substantial fraction of human regulatory sequences from transposable elements

Polymorphic mobile element insertions contribute to gene expression and alternative splicing in human tissues

BACKGROUND

Mobile elements are a major source of structural variants in the human genome, and some mobile elements can regulate gene expression and transcript splicing. However, the impact of polymorphic mobile element insertions (pMEIs) on gene expression and splicing in diverse human tissues has not been thoroughly studied. The multi-tissue gene expression and whole genome sequencing data generated by the Genotype-Tissue Expression (GTEx) project provide a great opportunity to systematically evaluate the role of pMEIs in regulating gene expression in human tissues.

RESULTS

Using the GTEx whole genome sequencing data, we identify 20,545 high-quality pMEIs from 639 individuals. Coupling pMEI genotypes with gene expression profiles, we identify pMEI-associated expression quantitative trait loci (eQTLs) and splicing quantitative trait loci (sQTLs) in 48 tissues. Using joint analyses of pMEIs and other genomic variants, pMEIs are predicted to be the potential causal variant for 3522 eQTLs and 3717 sQTLs. The pMEI-associated eQTLs and sQTLs show a high level of tissue specificity, and these pMEIs are enriched in the proximity of affected genes and in regulatory elements. Using reporter assays, we confirm that several pMEIs associated with eQTLs and sQTLs can alter gene expression levels and isoform proportions, respectively.

CONCLUSION

Overall, our study shows that pMEIs are associated with thousands of gene expression and splicing variations, indicating that pMEIs could have a significant role in regulating tissue-specific gene expression and transcript splicing. Detailed mechanisms for the role of pMEIs in gene regulation in different tissues will be an important direction for future studies.

DNA methylation enables transposable element-driven genome expansion

Multicellular eukaryotic genomes show enormous differences in size. A substantial part of this variation is due to the presence of transposable elements (TEs). They contribute significantly to a cell's mass of DNA and have the potential to become involved in host gene control. We argue that the suppression of their activities by methylation of the C-phosphate-G (CpG) dinucleotide in DNA is essential for their long-term accommodation in the host genome and, therefore, to its expansion. An inevitable consequence of cytosine methylation is an increase in C-to-T transition mutations via deamination, which causes CpG loss. Cytosine deamination is often needed for TEs to take on regulatory functions in the host genome. Our study of the whole-genome sequences of 53 organisms showed a positive correlation between the size of a genome and the percentage of TEs it contains, as well as a negative correlation between size and the CpG observed/expected (O/E) ratio in both TEs and the host DNA. TEs are seldom found at promoters and transcription start sites, but they are found more at enhancers, particularly after they have accumulated C-to-T and other mutations. Therefore, the methylation of TE DNA allows for genome expansion and also leads to new opportunities for gene control by TE-based regulatory sites.

Transposable elements, circular RNAs and mitochondrial transcription in age-related genomic regulation

Our understanding of the molecular regulation of aging and age-related diseases is still in its infancy, requiring in-depth characterization of the molecular landscape shaping these complex phenotypes. Emerging classes of molecules with promise as aging modulators include transposable elements, circRNAs and the mitochondrial transcriptome. Analytical complexity means that these molecules are often overlooked, even though they exhibit strong associations with aging and, in some cases, may directly contribute to its progress. Here, we review the links between these novel factors and age-related phenotypes, and we suggest tools that can be easily incorporated into existing pipelines to better understand the aging process.

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'.

Dynamic Signatures of the Epigenome: Friend or Foe?

Highly dynamic epigenetic signaling is influenced mainly by (micro)environmental stimuli and genetic factors. The exact mechanisms affecting particular epigenomic patterns differ dependently on the context. In the current review, we focus on the causes and effects of the dynamic signatures of the human epigenome as evaluated with the high-throughput profiling data and single-gene approaches. We will discuss three different aspects of phenotypic outcomes occurring as a consequence of epigenetics interplaying with genotype and environment. The first issue is related to the cases of environmental impacts on epigenetic profile, and its adverse and advantageous effects related to human health and evolutionary adaptation. The next topic will present a model of the interwoven co-evolution of genetic and epigenetic patterns exemplified with transposable elements (TEs) and their epigenetic repressors Krüppel-associated box zinc finger proteins (KRAB–ZNFs). The third aspect concentrates on the mitosis-based microevolution that takes place during carcinogenesis, leading to clonal diversity and expansion of tumor cells. The whole picture of epigenome plasticity and its role in distinct biological processes is still incomplete. However, accumulating data define epigenomic dynamics as an essential co-factor driving adaptation at the cellular and inter-species levels with a benefit or disadvantage to the host.

Transposable elements as a potent source of diverse cis-regulatory sequences in mammalian genomes

Eukaryotic gene regulation is mediated by cis-regulatory elements, which are embedded within the vast non-coding genomic space and recognized by the transcription factors in a sequence- and context-dependent manner. A large proportion of eukaryotic genomes, including at least half of the human genome, are composed of transposable elements (TEs), which in their ancestral form carried their own cis-regulatory sequences able to exploit the host trans environment to promote TE transcription and facilitate transposition. Although not all present-day TE copies have retained this regulatory function, the preexisting regulatory potential of TEs can provide a rich source of cis-regulatory innovation for the host. Here, we review recent evidence documenting diverse contributions of TE sequences to gene regulation by functioning as enhancers, promoters, silencers and boundary elements. We discuss how TE-derived enhancer sequences can rapidly facilitate changes in existing gene regulatory networks and mediate species- and cell-type-specific regulatory innovations, and we postulate a unique contribution of TEs to species-specific gene expression divergence in pluripotency and early embryogenesis. With advances in genome-wide technologies and analyses, systematic investigation of TEs' cis-regulatory potential is now possible and our understanding of the biological impact of genomic TEs is increasing. This article is part of a discussion meeting issue 'Crossroads between transposons and gene regulation'.

Viral Symbiosis in the Origins and Evolution of Life with a Particular Focus on the Placental Mammals

Advances in understanding over the last decade or so highlight the need for a reappraisal of the role of viruses in relation to the origins and evolution of cellular life, as well as in the homeostasis of the biosphere on which all of life depends. The relevant advances have, in particular, revealed an important contribution of viruses to the evolution of the placental mammals, while also contributing key roles to mammalian embryogenesis, genomic evolution, and physiology. Part of this reappraisal will include the origins of viruses, a redefinition of their quintessential nature, and a suggestion as to how we might view viruses in relation to the tree of life.

The role of transposable elements activity in aging and their possible involvement in laminopathic diseases

Eukaryotic genomes contain a large number of transposable elements, part of which are still active and able to transpose in the host genome. Mobile element activation is repressed to avoid deleterious effects, such as gene mutations or chromosome rearrangements. Control of transposable elements includes a variety of mechanisms comprising silencing pathways, which are based on the production of small non-coding RNAs. Silencing can occur either through transposable element RNA degradation or through the targeting of DNA sequences by heterochromatin formation and consequent transcriptional inhibition. Since the important role of the heterochromatin silencing, the gradual loss of heterochromatin marks in constitutive heterochromatin regions during the aging process promotes derepression of transposable elements, which is considered a cause of the progressive increase in genomic instability and of the activation of inflammatory responses. This review provides an overview of the effects of heterochromatin loss on the activity of transposable elements during the aging process and the possible impact on genome function. In this context, we discuss the possible role of the nuclear lamina, a major player in heterochromatin dynamics, in the regulation of transposable element activity and potential implications in laminopathic diseases.

The SUMO Ligase Su(var)2-10 Controls Hetero- and Euchromatic Gene Expression via Establishing H3K9 Trimethylation and Negative Feedback Regulation

Сhromatin is critical for genome compaction and gene expression. On a coarse scale, the genome is divided into euchromatin, which harbors the majority of genes and is enriched in active chromatin marks, and heterochromatin, which is gene-poor but repeat-rich. The conserved molecular hallmark of heterochromatin is the H3K9me3 modification, which is associated with gene silencing. We found that in Drosophila, deposition of most of the H3K9me3 mark depends on SUMO and the SUMO ligase Su(var)2-10, which recruits the histone methyltransferase complex SetDB1/Wde. In addition to repressing repeats, H3K9me3 influences expression of both hetero- and euchromatic host genes. High H3K9me3 levels in heterochromatin are required to suppress spurious transcription and ensure proper gene expression. In euchromatin, a set of conserved genes is repressed by Su(var)2-10/SetDB1-induced H3K9 trimethylation, ensuring tissue-specific gene expression. Several components of heterochromatin are themselves repressed by this pathway, providing a negative feedback mechanism to ensure chromatin homeostasis.

Diversification of the Caenorhabditis heat shock response by Helitron transposable elements

Heat Shock Factor 1 (HSF-1) is a key regulator of the heat shock response (HSR). Upon heat shock, HSF-1 binds well-conserved motifs, called Heat Shock Elements (HSEs), and drives expression of genes important for cellular protection during this stress. Remarkably, we found that substantial numbers of HSEs in multiple Caenorhabditis species reside within Helitrons, a type of DNA transposon. Consistent with Helitron-embedded HSEs being functional, upon heat shock they display increased HSF-1 and RNA polymerase II occupancy and up-regulation of nearby genes in C. elegans. Interestingly, we found that different genes appear to be incorporated into the HSR by species-specific Helitron insertions in C. elegans and C. briggsae and by strain-specific insertions among different wild isolates of C. elegans. Our studies uncover previously unidentified targets of HSF-1 and show that Helitron insertions are responsible for rewiring and diversifying the Caenorhabditis HSR.

TEffectR: an R package for studying the potential effects of transposable elements on gene expression with linear regression model

Introduction

Recent studies highlight the crucial regulatory roles of transposable elements (TEs) on proximal gene expression in distinct biological contexts such as disease and development. However, computational tools extracting potential TE -proximal gene expression associations from RNA-sequencing data are still missing.

Implementation

Herein, we developed a novel R package, using a linear regression model, for studying the potential influence of TE species on proximal gene expression from a given RNA-sequencing data set. Our R package, namely TEffectR, makes use of publicly available RepeatMasker TE and Ensembl gene annotations as well as several functions of other R-packages. It calculates total read counts of TEs from sorted and indexed genome aligned BAM files provided by the user, and determines statistically significant relations between TE expression and the transcription of nearby genes under diverse biological conditions.

Availability

TEffectR is freely available at https://github.com/karakulahg/TEffectR along with a handy tutorial as exemplified by the analysis of RNA-sequencing data including normal and tumour tissue specimens obtained from breast cancer patients.

Impact of polymorphic transposable elements on transcription in lymphoblastoid cell lines from public data

BACKGROUND

Transposable elements (TEs) are DNA sequences able to mobilize themselves and to increase their copy-number in the host genome. In the past, they have been considered mainly selfish DNA without evident functions. Nevertheless, currently they are believed to have been extensively involved in the evolution of primate genomes, especially from a regulatory perspective. Due to their recent activity they are also one of the primary sources of structural variants (SVs) in the human genome. By taking advantage of sequencing technologies and bioinformatics tools, recent surveys uncovered specific TE structural variants (TEVs) that gave rise to polymorphisms in human populations. When combined with RNA-seq data this information provides the opportunity to study the potential impact of TEs on gene expression in human.

RESULTS

In this work, we assessed the effects of the presence of specific TEs in cis on the expression of flanking genes by producing associations between polymorphic TEs and flanking gene expression levels in human lymphoblastoid cell lines. By using public data from the 1000 Genome Project and the Geuvadis consortium, we exploited an expression quantitative trait loci (eQTL) approach integrated with additional bioinformatics data mining analyses. We uncovered human loci enriched for common, less common and rare TEVs and identified 323 significant TEV-cis-eQTL associations. SINE-R/VNTR/Alus (SVAs) resulted the TE class with the strongest effects on gene expression. We also unveiled differential functional enrichments on genes associated to TEVs, genes associated to TEV-cis-eQTLs and genes associated to the genomic regions mostly enriched in TEV-cis-eQTLs highlighting, at multiple levels, the impact of TEVs on the host genome. Finally, we also identified polymorphic TEs putatively embedded in transcriptional units, proposing a novel mechanism in which TEVs may mediate individual-specific traits.

CONCLUSION

We contributed to unveiling the effect of polymorphic TEs on transcription in lymphoblastoid cell lines.

The Selfishness of Law-Abiding Genes

Selfish genes were once controversial, but it is now accepted that the genome contains parasitic elements in addition to a complement of conventional genes. This opinion article argues that 'law-abiding' genes also indulge in game playing to ensure their propagation, so that initially nonessential processes secure a genetic heritage. A gene-centered view of this kind can help to explain otherwise puzzling aspects of biology, including the complexity and stability of living systems.

DNA methylation of individual repetitive elements in hepatitis C virus infection-induced hepatocellular carcinoma

Background

The two most common repetitive elements (REs) in humans, long interspersed nuclear element-1 (LINE-1) and Alu element (Alu), have been linked to various cancers. Hepatitis C virus (HCV) may cause hepatocellular carcinoma (HCC) by suppressing host defenses, through DNA methylation that controls the mobilization of REs. We aimed to investigate the role of RE methylation in HCV-induced HCC (HCV-HCC).

Results

We studied methylation of over 30,000 locus-specific REs across the genome in HCC, cirrhotic, and healthy liver tissues obtained by surgical resection. Relative to normal liver tissue, we observed the largest number of differentially methylated REs in HCV-HCC followed by alcohol-induced HCC (EtOH-HCC). After excluding EtOH-HCC-associated RE methylation (FDR < 0.001) and those unable to be validated in The Cancer Genome Atlas (TCGA), we identified 13 hypomethylated REs (11 LINE-1 and 2 Alu) and 2 hypermethylated REs (1 LINE-1 and 1 Alu) in HCV-HCC (FDR < 0.001). A majority of these REs were located in non-coding regions, preferentially enriched with chromatin repressive marks H3K27me3, and positively associated with gene expression (median correlation r = 0.32 across REs). We further constructed an HCV-HCC RE methylation score that distinguished HCV-HCC (lowest score), HCV-cirrhosis, and normal liver (highest score) in a dose-responsive manner (p for trend < 0.001). HCV-cirrhosis had a lower score than EtOH-cirrhosis (p = 0.038) and HCV-HCC had a lower score than EtOH-HCC in TCGA (p = 0.024).

Conclusions

Our findings indicate that HCV infection is associated with loss of DNA methylation in specific REs, which could implicate molecular mechanisms in liver cancer development. If our findings are validated in larger sample sizes, methylation of these REs may be useful as an early detection biomarker for HCV-HCC and/or a target for prevention of HCC in HCV-positive individuals.

Electronic supplementary material

The online version of this article (10.1186/s13148-019-0733-y) contains supplementary material, which is available to authorized users.

SQuIRE reveals locus-specific regulation of interspersed repeat expression

Transposable elements (TEs) are interspersed repeat sequences that make up much of the human genome. Their expression has been implicated in development and disease. However, TE-derived RNA-seq reads are difficult to quantify. Past approaches have excluded these reads or aggregated RNA expression to subfamilies shared by similar TE copies, sacrificing quantitative accuracy or the genomic context necessary to understand the basis of TE transcription. As a result, the effects of TEs on gene expression and associated phenotypes are not well understood. Here, we present Software for Quantifying Interspersed Repeat Expression (SQuIRE), the first RNA-seq analysis pipeline that provides a quantitative and locus-specific picture of TE expression (https://github.com/wyang17/SQuIRE). SQuIRE is an accurate and user-friendly tool that can be used for a variety of species. We applied SQuIRE to RNA-seq from normal mouse tissues and a Drosophila model of amyotrophic lateral sclerosis. In both model organisms, we recapitulated previously reported TE subfamily expression levels and revealed locus-specific TE expression. We also identified differences in TE transcription patterns relating to transcript type, gene expression and RNA splicing that would be lost with other approaches using subfamily-level analyses. Altogether, our findings illustrate the importance of studying TE transcription with locus-level resolution.

Lung genotoxicity of benzo(a)pyrene in vivo involves reactivation of LINE-1 retrotransposon and early reprogramming of oncogenic regulatory networks

Several lines of evidence have implicated long interspersed nuclear element-1 (LINE-1) retroelement in the onset and progression of lung cancer. Retrotransposition-dependent mechanisms leading to DNA mobilization give rise to insertion mutations and DNA deletions, whereas retrotransposition-independent mechanisms disrupt epithelial programming and differentiation. Previous work by our group established that tobacco carcinogens such as benzo(a)pyrene (BaP) reactivate LINE-1 in bronchial epithelial cells through displacement of nucleosome remodeling and deacetylase (NuRD) corepressor complexes and interference with retinoblastoma-regulated epigenetic signaling. Whether LINE-1 in coordination with other genes within its regulatory network contributes to the in vivo genotoxic response to BaP remains largely unknown. Evidence is presented here that intratracheal instillation of ORFeus^LSL mice with BaP alone or in combination with adenovirus (adeno)-CRE recombinase is genotoxic to the lung and associated with activation of the human LINE-1 transgene present in these mice. LINE-1 reactivation modulated the expression of genes involved in oncogenic signaling, and these responses were most pronounced in female mice compared with males and synergized by adeno-CRE recombinase. This is the first report linking LINE-1 and genes within its oncogenic regulatory network with early sexually dimorphic responses of the lung in vivo.

Transposable Elements Adaptive Role in Genome Plasticity, Pathogenicity and Evolution in Fungal Phytopathogens

Transposable elements (TEs) are agents of genetic variability in phytopathogens as they are a source of adaptive evolution through genome diversification. Although many studies have uncovered information on TEs, the exact mechanism behind TE-induced changes within the genome remains poorly understood. Furthermore, convergent trends towards bigger genomes, emergence of novel genes and gain or loss of genes implicate a TE-regulated genome plasticity of fungal phytopathogens. TEs are able to alter gene expression by revamping the cis-regulatory elements or recruiting epigenetic control. Recent findings show that TEs recruit epigenetic control on the expression of effector genes as part of the coordinated infection strategy. In addition to genome plasticity and diversity, fungal pathogenicity is an area of economic concern. A survey of TE distribution suggests that their proximity to pathogenicity genes TEs may act as sites for emergence of novel pathogenicity factors via nucleotide changes and expansion or reduction of the gene family. Through a systematic survey of literature, we were able to conclude that the role of TEs in fungi is wide: ranging from genome plasticity, pathogenicity to adaptive behavior in evolution. This review also identifies the gaps in knowledge that requires further elucidation for a better understanding of TEs' contribution to genome architecture and versatility.

Human transposons are an abundant supply of transcription factor binding sites and promoter activities in breast cancer cell lines

Background

Transposable elements (TE) are commonly regarded as “junk DNA” with no apparent regulatory roles in the human genome. However, a growing body of evidence demonstrates that some TEs exhibit regulatory activities in a range of biological pathways and diseases, with notable examples in bile metabolism and innate immunity. TEs are typically suppressed by epigenetic modifications in healthy somatic tissues, which prevents both undesirable effects of insertional mutagenesis, and also unwanted gene activation. Interestingly, TEs are widely reported to be dysregulated in epithelial cancers, and while much attention has been paid to their effects on genome instability, relatively little has been reported on their effects on gene regulation. Here, we investigated the contribution of TEs to the transcriptional regulation in breast cancer cell lines.

Results

We found that a subset of TE subfamilies were enriched in oncogenic transcription factor binding sites and also harboured histone marks associated with active transcription, raising the possibility of these subfamilies playing a broad role in breast cancer transcriptional regulation. To directly assess promoter activity in triple negative breast cancer cell lines, we identified four breast cancer-associated genes with putative TE-derived promoters. TE deletion confirmed a contribution to promoter activity in all cases, and for two examples the promoter activity was almost completely contained within the TE.

Conclusions

Our findings demonstrate that TEs provide abundant oncogenic transcription factor binding sites in breast cancer and that individual TEs contain substantial promoter activity. Our findings provide further evidence for transcriptional regulation of human genes through TE exaptation by demonstrating the regulatory potential of TEs in multiple breast cancer cell lines.

Electronic supplementary material

The online version of this article (10.1186/s13100-019-0158-3) contains supplementary material, which is available to authorized users.

Chromatin Evolution-Key Innovations Underpinning Morphological Complexity

The history of life consists of a series of major evolutionary transitions, including emergence and radiation of complex multicellular eukaryotes from unicellular ancestors. The cells of multicellular organisms, with few exceptions, contain the same genome, however, their organs are composed of a variety of cell types that differ in both structure and function. This variation is largely due to the transcriptional activity of different sets of genes in different cell types. This indicates that complex transcriptional regulation played a key role in the evolution of complexity in eukaryotes. In this review, we summarize how gene duplication and subsequent evolutionary innovations, including the structural evolution of nucleosomes and chromatin-related factors, contributed to the complexity of the transcriptional system and provided a basis for morphological diversity.

Epigenetic Components of Myalgic Encephalomyelitis/Chronic Fatigue Syndrome Uncover Potential Transposable Element Activation

PURPOSE

Studies to determine epigenetic changes associated with myalgic encephalomyelitis/chronic fatigue syndrome (ME/CFS) remain scarce; however, current evidence clearly shows that methylation patterns of genomic DNA and noncoding RNA profiles of immune cells differ between patients and healthy subjects, suggesting an active role of these epigenetic mechanisms in the disease. The present study compares and contrasts the available ME/CFS epigenetic data in an effort to evidence overlapping pathways capable of explaining at least some of the dysfunctional immune parameters linked to this disease.

METHODS

A systematic search of the literature evaluating the ME/CFS epigenome landscape was performed following the Preferred Reporting Items for Systematic Reviews and Meta-Analyses criteria. Differential DNA methylation and noncoding RNA differential expression patterns associated with ME/CFS were used to screen for the presence of transposable elements using the Dfam browser, a search program nurtured with the Repbase repetitive sequence database and the RepeatMasker annotation tool.

FINDINGS

Unexpectedly, particular associations of transposable elements and ME/CFS epigenetic hallmarks were uncovered. A model for the disease emerged involving transcriptional induction of endogenous dormant transposons and structured cellular RNA interactions, triggering the activation of the innate immune system without a concomitant active infection.

IMPLICATIONS

Repetitive sequence filters (ie, RepeatMasker) should be avoided when analyzing transcriptomic data to assess the potential participation of repetitive sequences ("junk repetitive DNA"), representing >45% of the human genome, in the onset and evolution of ME/CFS. In addition, transposable element screenings aimed at designing cost-effective, focused empirical assays that can confirm or disprove the suspected involvement of transposon transcriptional activation in this disease, following the pilot strategy presented here, will require databases gathering large ME/CFS epigenetic datasets.

Transposable Elements: Classification, Identification, and Their Use As a Tool For Comparative Genomics

Most genomes are populated by hundreds of thousands of sequences originated from mobile elements. On the one hand, these sequences present a real challenge in the process of genome analysis and annotation. On the other hand, they are very interesting biological subjects involved in many cellular processes. Here we present an overview of transposable elements biodiversity, and we discuss different approaches to transposable elements detection and analyses.

Transposable element insertions in fission yeast drive adaptation to environmental stress

Cells are regularly exposed to a range of naturally occurring stress that can restrict growth or cause lethality. In response, cells activate expression networks with hundreds of genes that together increase resistance to common environmental insults. However, stress response networks can be insufficient to ensure survival, which raises the question of whether cells possess genetic programs that can promote adaptation to novel forms of stress. We found transposable element (TE) mobility in Schizosaccharomyces pombe was greatly increased when cells were exposed to unusual forms of stress such as heavy metals, caffeine, and the plasticizer phthalate. By subjecting TE-tagged cells to CoCl₂, we found the TE integration provided the major path to resistance. Groups of insertions that provided resistance were linked to TOR regulation and metal response genes. We extended our study of adaptation by analyzing TE positions in 57 genetically distinct wild strains. The genomic positions of 1048 polymorphic LTRs were strongly associated with a range of stress response genes, indicating TE integration promotes adaptation in natural conditions. These data provide strong support for the idea, first proposed by Barbara McClintock, that TEs provide a system to modify the genome in response to stress.

Structure, clustering and functional insights of repeats configurations in the upstream promoter region of the human coding genes

Background

Repetitive DNA sequences (Repeats) are significant regions in the human genome that have a specific genomic distribution, structure, and several binding sites for genome architecture and function. In consequence, the possible configurations of Repeats in specific and dynamic regions like the gene promoters could define footprints for molecular mechanisms, pathways, and cell function beyond their density in the genome. Here we explored the distribution of Repeats in the upstream promoter region of the human coding genes with the aim to identify specific configurations, clusters and functional meaning of those elements. Our method includes structural descriptions, hierarchical clustering, pathway association, and functional enrichment analysis.

Results

We report here several configurations of Repeats in the upstream promoter region (UPR), which define 2729 patterns for the 80% of the human coding genes. There are 47 types of Repeats in these configurations, where the most frequent were Alu, Low_complexity, MIR, Simple_repeat, LINE/L2, LINE/L1, hAT-Charlie, and ERV1. The distribution, length, and the high frequency of Repeats in the UPR defines several patterns and clusters, where the minimum frequency of configuration among Repeats was higher than 0.7. We found those clusters associated with cellular pathways and ontologies; thus, it was plausible to determine groups of Repeats to specific functional insights, for example, pathways for Genetic Information Processing or Metabolism shows particular groups of Repeats with specific configurations.

Conclusion

Based on these findings, we propose that specific configurations of repetitive elements describe frequent patterns in the upstream promoter for sets of human coding genes, which those correlated to specific and essential cell pathways and functions.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-5196-6) contains supplementary material, which is available to authorized users.

Human Endogenous Retroviruses: Residues of Ancient Times Are Differentially Expressed in Crohn's Disease

Background

Eight percent of the human genome consists of human endogenous retroviruses (HERV). These genetic elements are remnants of ancient retroviral germ-line infections. Altered HERV expression is associated with several chronic inflammatory diseases. A physiological role of the HERV-derived proteins syncytin-1 and -2 has been described for the integrity of the human placental cell layer in terms of maintaining feto-maternal tolerance. The aim of this project was to investigate HERV expression in Crohn's disease (CD) with a further focus on syncytins in the gut.

Material and Methods

Seventy-four ileal and colonic tissue samples of CD patients and healthy controls have been investigated for mRNA expression of major HERV groups by a comprehensive microarray screening. The most prominent differences have been validated by qRT-PCR. Immunohistochemistry (IHC), Western Blot (WB) and qRT-PCR were performed for syncytin-1 and -2.

Results

HERV microarray screening revealed a distinct expression profile in ileal and colonic tissue, as well as differential expression in CD compared to healthy controls. qRT-PCR validated differential expression of at least 3 HERV-groups in CD. qRT-PCR, IHC and WB showed a tissue-dependent diminished epithelial expression of syncytins in inflamed CD.

Conclusion

For the first time, HERV expression has been comprehensively studied in the gut. Between CD and healthy controls we could show a tissue dependent differential HERV expression profile. Notably, we could show that syncytin-1 and -2 are expressed in the epithelial layer in ileal and colonic tissue samples, whereas their diminished tissue-dependent expression in inflamed CD might modulate inflammatory processes at the gut barrier.

Transposable elements generate regulatory novelty in a tissue-specific fashion

BACKGROUND

Transposable elements (TE) are an important source of evolutionary novelty in gene regulation. However, the mechanisms by which TEs contribute to gene expression are largely uncharacterized.

RESULTS

Here, we leverage Roadmap and GTEx data to investigate the association of TEs with active and repressed chromatin in 24 tissues. We find 112 human TE families enriched in active regions of the genome across tissues. Short Interspersed Nuclear Elements (SINEs) and DNA transposons are the most frequently enriched classes, while Long Terminal Repeat Retrotransposons (LTRs) are often enriched in a tissue-specific manner. We report across-tissue variability in TE enrichment in active regions. Genes with consistent expression across tissues are less likely to be associated with TE insertions. TE presence in repressed regions similarly follows tissue-specific patterns. Moreover, different TE classes correlate with different repressive marks: LTRs and Long Interspersed Nuclear Elements (LINEs) are overrepresented in regions marked by H3K9me3, while the other TEs are more likely to overlap regions with H3K27me3. Young TEs are typically enriched in repressed regions and depleted in active regions. We detect multiple instances of TEs that are enriched in tissue-specific active regulatory regions. Such TEs contain binding sites for transcription factors that are master regulators for the given tissue. These TEs are enriched in intronic enhancers, and their tissue-specific enrichment correlates with tissue-specific variations in the expression of the nearest genes.

CONCLUSIONS

We provide an integrated overview of the contribution of TEs to human gene regulation. Expanding previous analyses, we demonstrate that TEs can potentially contribute to the turnover of regulatory sequences in a tissue-specific fashion.

PARP1 activation increases expression of modified tumor suppressors and pathways underlying development of aggressive hepatoblastoma

Hepatoblastoma (HBL) is a pediatric liver cancer that affects children under the age of three. Reduction of tumor suppressor proteins (TSPs) is commonly seen in liver cancer. However, in our studies we find that aggressive, chemo-resistant HBLs exhibit an elevation of TSPs. HBL patients with a classic phenotype have reduced TSP levels, but patients with aggressive HBL express elevated TSPs that undergo posttranslational modifications, eliminating their tumor suppression activities. Here we identify unique aggressive liver cancer domains (ALCDs) that are activated in aggressive HBL by PARP1-mediated chromatin remodeling leading to elevation of modified TSPs and activation of additional cancer pathways: WNT signaling and β-catenin. Inhibition of PARP1 blocks activation of ALCDs and normalizes expression of corresponding genes, therefore reducing cell proliferation. Our studies reveal PARP1 activation as a mechanism for the development of aggressive HBL, further suggesting FDA-approved PARP1 inhibitors might be used for treatment of patients with aggressive HBL.

Leila Valanejad et al. report increased expression of modified tumor suppressor proteins (TSPs) with loss of tumor suppressor activity in aggressive, chemotherapy-resistant hepatoblastoma. They find that TSP upregulation occurs via PARP1-mediated chromatin remodeling, leading to activation of multiple cancer-associated pathways.

Repetitive Element DNA Methylation is Associated with Menopausal Age

To investigate associations between the age of menopause and the DNA methylation levels of two repetitive elements, Alu and LINE-1, we performed plasma DNA extraction on 161 subjects and serum cell-free DNA extraction on 120 subjects. We grouped women by menopausal age as follows: ≤ 48 years (earlier menopause), ≥ 52 years (later menopause), and 48-52 years (control). The DNA methylation levels of Alu and LINE-1 were measured by MethyLight PCR. The results showed that the DNA methylation levels of both Alu and LINE-1 were inversely correlated with menopausal age in the plasma DNA cohort (r = 0.079, P < 0.001 for Alu; r = 0.045, P = 0.007 for LINE-1) as well as in the serum DNA cohort (r = 0.087, P = 0.001 for Alu; r = 0.041, P = 0.026 for LINE-1). Alu methylation levels in both the plasma and serum DNA cohorts and LINE-1 methylation levels in the plasma cohort were remarkably higher in the earlier menopause group than in the later menopause and control groups (P < 0.01 and P < 0.05, respectively). In the serum DNA cohort, the LINE-1 methylation levels in the later menopause group were significantly lower than that in the earlier menopause group and control group (P < 0.05). Therefore, methylation levels of Alu and LINE-1 were significantly associated with menopausal age. Women with earlier menopause showed hypermethylation in both repetitive elements, while women with later menopause showed hypomethylation. These findings suggest that altered DNA methylation in leukocytes and serum cell-free DNA may represent a biomarker of menopausal age.

Identification of sex-specific DNA methylation changes driven by specific chemicals in cord blood in a Faroese birth cohort

Faroe islanders consume marine foods contaminated with methylmercury (MeHg), polychlorinated biphenyls (PCBs), and other toxicants associated with chronic disease risks. Differential DNA methylation at specific CpG sites in cord blood may serve as a surrogate biomarker of health impacts from chemical exposures. We aimed to identify key environmental chemicals in cord blood associated with DNA methylation changes in a population with elevated exposure to chemical mixtures. We studied 72 participants of a Faroese birth cohort recruited between 1986 and 1987 and followed until adulthood. The cord blood DNA methylome was profiled using Infinium HumanMethylation450 BeadChips. We determined the associations of CpG site changes with concentrations of MeHg, major PCBs, other organochlorine compounds [hexachlorobenzene (HCB), p,p'-dichlorodiphenyldichloroethylene (p,p'-DDE) and p,p'-dichlorodiphenyltrichloroethane], and perfluoroalkyl substances. In a combined sex analysis, among the 16 chemicals studied, PCB congener 105 (CB-105) exposure was associated with the majority of differentially methylated CpG sites (214 out of a total of 250). In female-only analysis, only 73 CB-105 associated CpG sites were detected, 44 of which were mapped to genes in the ELAV1-associated cancer network. In males-only, methylation changes were seen for perfluorooctane sulfonate, HCB, and p,p'-DDE in 10,598, 1,238, and 1,473 CpG sites, respectively, 15% of which were enriched in cytobands of the X-chromosome associated with neurological disorders. In this multiple-pollutant and genome-wide study, we identified key epigenetic toxicants. The significant enrichment of specific X-chromosome sites in males implies potential sex-specific epigenome responses to prenatal chemical exposures.

Contribution of transposable elements in the plant's genome

Plants maintain extensive growth flexibility under different environmental conditions, allowing them to continuously and rapidly adapt to alterations in their environment. A large portion of many plant genomes consists of transposable elements (TEs) that create new genetic variations within plant species. Different types of mutations may be created by TEs in plants. Many TEs can avoid the host's defense mechanisms and survive alterations in transposition activity, internal sequence and target site. Thus, plant genomes are expected to utilize a variety of mechanisms to tolerate TEs that are near or within genes. TEs affect the expression of not only nearby genes but also unlinked inserted genes. TEs can create new promoters, leading to novel expression patterns or alternative coding regions to generate alternate transcripts in plant species. TEs can also provide novel cis-acting regulatory elements that act as enhancers or inserts within original enhancers that are required for transcription. Thus, the regulation of plant gene expression is strongly managed by the insertion of TEs into nearby genes. TEs can also lead to chromatin modifications and thereby affect gene expression in plants. TEs are able to generate new genes and modify existing gene structures by duplicating, mobilizing and recombining gene fragments. They can also facilitate cellular functions by sharing their transposase-coding regions. Hence, TE insertions can not only act as simple mutagens but can also alter the elementary functions of the plant genome. Here, we review recent discoveries concerning the contribution of TEs to gene expression in plant genomes and discuss the different mechanisms by which TEs can affect plant gene expression and reduce host defense mechanisms.

Silencing of transposable elements may not be a major driver of regulatory evolution in primate iPSCs

Transposable elements (TEs) comprise almost half of primate genomes and their aberrant regulation can result in deleterious effects. In pluripotent stem cells, rapidly evolving KRAB-ZNF genes target TEs for silencing by H3K9me3. To investigate the evolution of TE silencing, we performed H3K9me3 ChIP-seq experiments in induced pluripotent stem cells from 10 human and 7 chimpanzee individuals. We identified four million orthologous TEs and found the SVA and ERV families to be marked most frequently by H3K9me3. We found little evidence of inter-species differences in TE silencing, with as many as 82% of putatively silenced TEs marked at similar levels in humans and chimpanzees. TEs that are preferentially silenced in one species are a similar age to those silenced in both species and are not more likely to be associated with expression divergence of nearby orthologous genes. Our data suggest limited species-specificity of TE silencing across 6 million years of primate evolution.

Transposable element activity, genome regulation and human health

A convergence of novel genome analysis technologies is enabling population genomic studies of human transposable elements (TEs). Population surveys of human genome sequences have uncovered thousands of individual TE insertions that segregate as common genetic variants, i.e. TE polymorphisms. These recent TE insertions provide an important source of naturally occurring human genetic variation. Investigators are beginning to leverage population genomic data sets to execute genome-scale association studies for assessing the phenotypic impact of human TE polymorphisms. For example, the expression quantitative trait loci (eQTL) analytical paradigm has recently been used to uncover hundreds of associations between human TE insertion variants and gene expression levels. These include population-specific gene regulatory effects as well as coordinated changes to gene regulatory networks. In addition, analyses of linkage disequilibrium patterns with previously characterized genome-wide association study (GWAS) trait variants have uncovered TE insertion polymorphisms that are likely causal variants for a variety of common complex diseases. Gene regulatory mechanisms that underlie specific disease phenotypes have been proposed for a number of these trait associated TE polymorphisms. These new population genomic approaches hold great promise for understanding how ongoing TE activity contributes to functionally relevant genetic variation within and between human populations.

Transposable elements shape the human proteome landscape via formation of cis-acting upstream open reading frames

Transposons are major drivers of mammalian genome evolution. To obtain new insights into the contribution of transposons to the regulation of protein translation, we here examined how transposons affected the genesis and function of upstream open reading frames (uORFs), which serve as cis-acting elements to regulate translation from annotated ORFs (anORFs) located downstream of the uORFs in eukaryotic mRNAs. Among 39,786 human uORFs, 3,992 had ATG trinucleotides of a transposon origin, termed "transposon-derived upstream ATGs" or TuATGs. Luciferase reporter assays suggested that many TuATGs modulate translation from anORFs. Comparisons with transposon consensus sequences revealed that most TuATGs were generated by nucleotide substitutions in non-ATG trinucleotides of integrated transposons. Among these non-ATG trinucleotides, GTG and ACG were converted into TuATGs more frequently, indicating a CpG methylation-mediated process of TuATG formation. Interestingly, it is likely that this process accelerated human-specific upstream ATG formation within transposon sequences in 5' untranslated regions after divergence between human and nonhuman primates. Methylation-mediated TuATG formation seems to be ongoing in the modern human population and could alter the expression of disease-related proteins. This study shows that transposons have potentially been shaping the human proteome landscape via cis-acting uORF creation.

Transposable Element Mediated Innovation in Gene Regulatory Landscapes of Cells: Re-Visiting the "Gene-Battery" Model

Transposable elements (TEs) are no longer considered to be "junk" DNA. Here, we review how TEs can impact gene regulation systematically. TEs encode various regulatory elements that enables them to regulate gene expression. RJ Britten and EH Davidson hypothesized that TEs can integrate the function of various transcriptional regulators into gene regulatory networks. Uniquely TEs can deposit regulatory sites across the genome when they transpose, and thereby bring multiple genes under control of the same regulatory logic. Several studies together have robustly established that TEs participate in embryonic development and oncogenesis. We discuss the regulatory characteristics of TEs in context of evolution to understand the extent of their impact on gene networks. Understanding these features of TEs is central to future investigations of TEs in cellular processes and phenotypic presentations, which are applicable to development and disease studies. We re-visit the Britten-Davidson "gene-battery" model and understand the genetic and transcriptional impact of TEs in innovating gene regulatory networks.

Living Organisms Author Their Read-Write Genomes in Evolution

Evolutionary variations generating phenotypic adaptations and novel taxa resulted from complex cellular activities altering genome content and expression: (i) Symbiogenetic cell mergers producing the mitochondrion-bearing ancestor of eukaryotes and chloroplast-bearing ancestors of photosynthetic eukaryotes; (ii) interspecific hybridizations and genome doublings generating new species and adaptive radiations of higher plants and animals; and, (iii) interspecific horizontal DNA transfer encoding virtually all of the cellular functions between organisms and their viruses in all domains of life. Consequently, assuming that evolutionary processes occur in isolated genomes of individual species has become an unrealistic abstraction. Adaptive variations also involved natural genetic engineering of mobile DNA elements to rewire regulatory networks. In the most highly evolved organisms, biological complexity scales with "non-coding" DNA content more closely than with protein-coding capacity. Coincidentally, we have learned how so-called "non-coding" RNAs that are rich in repetitive mobile DNA sequences are key regulators of complex phenotypes. Both biotic and abiotic ecological challenges serve as triggers for episodes of elevated genome change. The intersections of cell activities, biosphere interactions, horizontal DNA transfers, and non-random Read-Write genome modifications by natural genetic engineering provide a rich molecular and biological foundation for understanding how ecological disruptions can stimulate productive, often abrupt, evolutionary transformations.

The "evolutionary field" hypothesis. Non-Mendelian transgenerational inheritance mediates diversification and evolution

Epigenetics is increasingly regarded as a potential contributing factor to evolution. Building on apparently unrelated results, here I propose that RNA-containing nanovesicles, predominantly small regulatory RNAs, are released from somatic tissues in the bloodstream, cross the Weismann barrier, reach the epididymis, and are eventually taken up by spermatozoa; henceforth the information is delivered to oocytes at fertilization. In the model, a LINE-1-encoded reverse transcriptase activity, present in spermatozoa and early embryos, plays a key role in amplifying and propagating these RNAs as extrachromosomal structures. It may be conceived that, over generations, the cumulative effects of sperm-delivered RNAs would cross a critical threshold and overcome the buffering capacity of embryos. As a whole, the process can promote the generation of an information-containing platform that drives the reshaping of the embryonic epigenetic landscape with the potential to generate ontogenic changes and redirect the evolutionary trajectory. Over time, evolutionary significant, stably acquired variations could be generated through the process. The interplay between these elements defines the concept of "evolutionary field", a self-consistent, comprehensive information-containing platform and a source of discontinuous evolutionary novelty.

Inevitability of the emergence and persistence of genetic parasites caused by evolutionary instability of parasite-free states

Genetic parasites, including viruses and mobile genetic elements, are ubiquitous among cellular life forms, and moreover, are the most abundant biological entities on earth that harbor the bulk of the genetic diversity. Here we examine simple thought experiments to demonstrate that both the emergence of parasites in simple replicator systems and their persistence in evolving life forms are inevitable because the putative parasite-free states are evolutionarily unstable.

REVIEWERS

This article has been reviewed by Yitzhak Pilpel, Bojan Zagrovic, and Eric van Nimwegen.

Prediction of genome-wide DNA methylation in repetitive elements

DNA methylation in repetitive elements (RE) suppresses their mobility and maintains genomic stability, and decreases in it are frequently observed in tumor and/or surrogate tissues. Averaging methylation across RE in genome is widely used to quantify global methylation. However, methylation may vary in specific RE and play diverse roles in disease development, thus averaging methylation across RE may lose significant biological information. The ambiguous mapping of short reads by and high cost of current bisulfite sequencing platforms make them impractical for quantifying locus-specific RE methylation. Although microarray-based approaches (particularly Illumina's Infinium methylation arrays) provide cost-effective and robust genome-wide methylation quantification, the number of interrogated CpGs in RE remains limited. We report a random forest-based algorithm (and corresponding R package, REMP) that can accurately predict genome-wide locus-specific RE methylation based on Infinium array profiling data. We validated its prediction performance using alternative sequencing and microarray data. Testing its clinical utility with The Cancer Genome Atlas data demonstrated that our algorithm offers more comprehensively extended locus-specific RE methylation information that can be readily applied to large human studies in a cost-effective manner. Our work has the potential to improve our understanding of the role of global methylation in human diseases, especially cancer.

Transposable elements are the primary source of novelty in primate gene regulation

Gene regulation shapes the evolution of phenotypic diversity. We investigated the evolution of liver promoters and enhancers in six primate species using ChIP-seq (H3K27ac and H3K4me1) to profile cis-regulatory elements (CREs) and using RNA-seq to characterize gene expression in the same individuals. To quantify regulatory divergence, we compared CRE activity across species by testing differential ChIP-seq read depths directly measured for orthologous sequences. We show that the primate regulatory landscape is largely conserved across the lineage, with 63% of the tested human liver CREs showing similar activity across species. Conserved CRE function is associated with sequence conservation, proximity to coding genes, cell-type specificity, and transcription factor binding. Newly evolved CREs are enriched in immune response and neurodevelopmental functions. We further demonstrate that conserved CREs bind master regulators, suggesting that while CREs contribute to species adaptation to the environment, core functions remain intact. Newly evolved CREs are enriched in young transposable elements (TEs), including Long-Terminal-Repeats (LTRs) and SINE-VNTR-Alus (SVAs), that significantly affect gene expression. Conversely, only 16% of conserved CREs overlap TEs. We tested the cis-regulatory activity of 69 TE subfamilies by luciferase reporter assays, spanning all major TE classes, and showed that 95.6% of tested TEs can function as either transcriptional activators or repressors. In conclusion, we demonstrated the critical role of TEs in primate gene regulation and illustrated potential mechanisms underlying evolutionary divergence among the primate species through the noncoding genome.

LINE-1 hypermethylation in white blood cell DNA is associated with high-grade cervical intraepithelial neoplasia

Background

Long Interspersed Nuclear Elements-1 (LINEs-1) methylation from white blood cells (WBCs) DNA has been proposed as biomarker associated with different types of cancer. The aim of the present study was to investigate the degree of WBCs LINE-1 methylation, according to high-risk Human Papilloma Virus (hrHPV) status in a healthy population, and the association with high-grade Cervical Intraepithelial Neoplasia (CIN2+) in hrHPV positive women.

Methods

Women with abnormal cervical cells were enrolled and classified by histological diagnosis and hrHPV infection. A structured questionnaire was used to obtain information on socio-demographic variables and lifestyle factors. LINE-1 methylation level in WBCs was measured by pyrosequencing-based methylation analysis after bisulfite conversion.

Results

Among 252 women diagnosed with normal cervical epithelium, with regard to LINE-1 methylation level no significant difference was observed between hrHPV positive and hrHPV negative women, also adjusting for known risk factors of infection. The association between WBCs LINE-1 methylation and CIN2+ status was analyzed in hrHPV positive women. The median value of LINE-1 methylation levels was higher in cases (CIN2+) than in controls (75.00% versus 73.17%; p = 0.002). For a one-unit increase in LINE-1 methylation level, the odds of being diagnosed with CIN2+ increased by 10%, adjusting for known factors related to LINE-1 methylation (adjOR: 1.10; 95% CI:1.01–1.20; p = 0.032). The Receiver-Operating Characteristic (ROC) curve analysis identified the cut-off value of 73.8% as the best threshold to separate cases from controls (sensitivity: 63.4% and specificity: 61.8%).

Conclusions

LINE-1 methylation status in WBCs DNA may represent a cost-effective and tissue-accessible biomarker for high-grade CIN in hrHPV positive women. However, LINE-1 hypermethylation cannot be considered specific for cervical cancer (CC) and a model based solely on LINE-1 methylation levels has limited performance. Further investigations are necessary to propose and validate a novel methylation biomarker panel, based on LINE-1 methylation and other differentially methylated regions, for the screening of women at risk of CC.

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.

Systematic identification and characterization of regulatory elements derived from human endogenous retroviruses

Human endogenous retroviruses (HERVs) and other long terminal repeat (LTR)-type retrotransposons (HERV/LTRs) have regulatory elements that possibly influence the transcription of host genes. We systematically identified and characterized these regulatory elements based on publicly available datasets of ChIP-Seq of 97 transcription factors (TFs) provided by ENCODE and Roadmap Epigenomics projects. We determined transcription factor-binding sites (TFBSs) using the ChIP-Seq datasets and identified TFBSs observed on HERV/LTR sequences (HERV-TFBSs). Overall, 794,972 HERV-TFBSs were identified. Subsequently, we identified "HERV/LTR-shared regulatory element (HSRE)," defined as a TF-binding motif in HERV-TFBSs, shared within a substantial fraction of a HERV/LTR type. HSREs could be an indication that the regulatory elements of HERV/LTRs are present before their insertions. We identified 2,201 HSREs, comprising specific associations of 354 HERV/LTRs and 84 TFs. Clustering analysis showed that HERV/LTRs can be grouped according to the TF binding patterns; HERV/LTR groups bounded to pluripotent TFs (e.g., SOX2, POU5F1, and NANOG), embryonic endoderm/mesendoderm TFs (e.g., GATA4/6, SOX17, and FOXA1/2), hematopoietic TFs (e.g., SPI1 (PU1), GATA1/2, and TAL1), and CTCF were identified. Regulatory elements of HERV/LTRs tended to locate nearby and/or interact three-dimensionally with the genes involved in immune responses, indicating that the regulatory elements play an important role in controlling the immune regulatory network. Further, we demonstrated subgroup-specific TF binding within LTR7, LTR5B, and LTR5_Hs, indicating that gains or losses of the regulatory elements occurred during genomic invasions of the HERV/LTRs. Finally, we constructed dbHERV-REs, an interactive database of HERV/LTR regulatory elements (http://herv-tfbs.com/). This study provides fundamental information in understanding the impact of HERV/LTRs on host transcription, and offers insights into the transcriptional modulation systems of HERV/LTRs and ancestral HERVs.

Useful parasites: the evolutionary biology and biotechnology applications of transposable elements

Transposable elements usually comprise the most abundant nongenic fraction of eukaryotic genomes. Because of their capacity to selfreplicate and to induce a wide range of mutations, transposable elements have long been considered as 'parasitic' or 'selfish'. Today, we recognize that the findings about genomic changes affected by transposable elements have considerably altered our view of the ways in which genomes evolve and work. Numerous studies have provided evidences that mobile elements have the potential to act as agents of evolution by increasing, rearranging and diversifying the genetic repertoire of their hosts. With large-scale sequencing becoming increasingly available, more and more scientists come across transposable element sequences in their data. I will provide examples that transposable elements, although having signatures of 'selfish' DNA, play a significant biological role in the maintainance of genome integrity and providing novel regulatoty networks. These features, along with the transpositional and mutagenic capacity to produce a raw genetic diversity, make the genome mobile fraction, a key player in species adaptation and microevolution. The last but not least, transposable elements stand as informative DNA markers that may complement other conventional DNA markers. Altogether, transposable elements represent a promising, but still largely unexplored research niche and deserve to be included into the agenda of molecular ecologists, evolutionary geneticists, conservation biologists and plant breeders.

Human population-specific gene expression and transcriptional network modification with polymorphic transposable elements

Transposable element (TE) derived sequences are known to contribute to the regulation of the human genome. The majority of known TE-derived regulatory sequences correspond to relatively ancient insertions, which are fixed across human populations. The extent to which human genetic variation caused by recent TE activity leads to regulatory polymorphisms among populations has yet to be thoroughly explored. In this study, we searched for associations between polymorphic TE (polyTE) loci and human gene expression levels using an expression quantitative trait loci (eQTL) approach. We compared locus-specific polyTE insertion genotypes to B cell gene expression levels among 445 individuals from 5 human populations. Numerous human polyTE loci correspond to both cis and trans eQTL, and their regulatory effects are directly related to cell type-specific function in the immune system. PolyTE loci are associated with differences in expression between European and African population groups, and a single polyTE loci is indirectly associated with the expression of numerous genes via the regulation of the B cell-specific transcription factor PAX5. The polyTE-gene expression associations we found indicate that human TE genetic variation can have important phenotypic consequences. Our results reveal that TE-eQTL are involved in population-specific gene regulation as well as transcriptional network modification.