Regulatory activities of transposable elements: from conflicts to benefits

Genome-wide mapping of transcriptional enhancer candidates using DNA and chromatin features in maize

Background

While most cells in multicellular organisms carry the same genetic information, in each cell type only a subset of genes is being transcribed. Such differentiation in gene expression depends, for a large part, on the activation and repression of regulatory sequences, including transcriptional enhancers. Transcriptional enhancers can be located tens of kilobases from their target genes, but display characteristic chromatin and DNA features, allowing their identification by genome-wide profiling. Here we show that integration of chromatin characteristics can be applied to predict distal enhancer candidates in Zea mays, thereby providing a basis for a better understanding of gene regulation in this important crop plant.

Result

To predict transcriptional enhancers in the crop plant maize (Zea mays L. ssp. mays), we integrated available genome-wide DNA methylation data with newly generated maps for chromatin accessibility and histone 3 lysine 9 acetylation (H3K9ac) enrichment in young seedling and husk tissue. Approximately 1500 intergenic regions, displaying low DNA methylation, high chromatin accessibility and H3K9ac enrichment, were classified as enhancer candidates. Based on their chromatin profiles, candidate sequences can be classified into four subcategories. Tissue-specificity of enhancer candidates is defined based on the tissues in which they are identified and putative target genes are assigned based on tissue-specific expression patterns of flanking genes.

Conclusions

Our method identifies three previously identified distal enhancers in maize, validating the new set of enhancer candidates and enlarging the toolbox for the functional characterization of gene regulation in the highly repetitive maize genome.

Electronic supplementary material

The online version of this article (doi:10.1186/s13059-017-1273-4) contains supplementary material, which is available to authorized users.

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.

Friends-Enemies: Endogenous Retroviruses Are Major Transcriptional Regulators of Human DNA

Endogenous retroviruses are mobile genetic elements hardly distinguishable from infectious, or “exogenous,” retroviruses at the time of insertion in the host DNA. Human endogenous retroviruses (HERVs) are not rare. They gave rise to multiple families of closely related mobile elements that occupy ~8% of the human genome. Together, they shape genomic regulatory landscape by providing at least ~320,000 human transcription factor binding sites (TFBS) located on ~110,000 individual HERV elements. The HERVs host as many as 155,000 mapped DNaseI hypersensitivity sites, which denote loci active in the regulation of gene expression or chromatin structure. The contemporary view of the HERVs evolutionary dynamics suggests that at the early stages after insertion, the HERV is treated by the host cells as a foreign genetic element, and is likely to be suppressed by the targeted methylation and mutations. However, at the later stages, when significant number of mutations has been already accumulated and when the retroviral genes are broken, the regulatory potential of a HERV may be released and recruited to modify the genomic balance of transcription factor binding sites. This process goes together with further accumulation and selection of mutations, which reshape the regulatory landscape of the human DNA. However, developmental reprogramming, stress or pathological conditions like cancer, inflammation and infectious diseases, can remove the blocks limiting expression and HERV-mediated host gene regulation. This, in turn, can dramatically alter the gene expression equilibrium and shift it to a newer state, thus further amplifying instability and exacerbating the stressful situation.

EnHERV: Enrichment analysis of specific human endogenous retrovirus patterns and their neighboring genes

Human endogenous retroviruses (HERVs) are flanked by long terminal repeats (LTRs), which contain the regulation part of the retrovirus. Remaining HERVs constitute 7% to 8% of the present day human genome, and most have been identified as solo LTRs. The HERV sequences have been associated with several molecular functions as well as certain diseases in human, but their roles in human diseases are yet to be established. We designed EnHERV to make accessible the identified endogenous retrovirus repetitive sequences from Repbase Update (a database of eukaryotic repetitive elements) that are present in the human genome. Defragmentation process was done to improve the RepeatMasker annotation output. The defragmented elements were used as core database in EnHERV. EnHERV is available at http://sysbio.chula.ac.th/enherv and can be searched using either gene lists of user interest or HERV characteristics. Besides the search function, EnHERV also provides an enrichment analysis function that allows users to perform enrichment analysis between selected HERV characteristics and user-input gene lists, especially genes with the expression profile of a certain disease. EnHERV will facilitate exploratory studies of specific HERV characteristics that control gene expression patterns related to various disease conditions. Here we analyzed 25 selected HERV groups/names from all four HERV superfamilies, using the sense and anti-sense directions of the HERV and gene expression profiles from 49 specific tissue and disease conditions. We found that intragenic HERVs were associated with down-regulated genes in most cancer conditions and in psoriatic skin tissues and associated with up-regulated genes in immune cells particularly from systemic lupus erythematosus (SLE) patients. EnHERV allowed the analysis of how different types of LTRs were differentially associated with specific gene expression profiles in particular disease conditions for further studies into their mechanisms and functions.

Somatizing the transposons action

The somatic mobilization of transposable elements is more common than previously thought. In this review we discuss how the intensity and the biologic consequences of somatic mobilization are dependent on the transposable elements landscapes of each genome, and on the "momentum" of each particular TE with respect to the mechanisms that control its transposition and the possibility to escape this control. Additionally, the biologic consequences of somatic mobilization vary among organisms that show an early separation between the germline and somatic cells and those organisms that do not exhibit this separation or that reproduce asexually. In the former, somatic transposition can be involved in phenotypic plasticity, detrimental conditions such as disease, or processes such as aging. For the organisms without separation between the germ and soma, somatic mobilization can be a source of genetic variability.

Genetic and epigenetic variation of transposable elements in Arabidopsis

Transposable elements are mobile genetic elements that are prevalent in plant genomes and are silenced by epigenetic modification. Different epigenetic modification pathways play distinct roles in the control of transposable element transcription, replication and recombination. The Arabidopsis genome contains families of all of the major transposable element classes, which are differentially enriched in particular genomic regions. Whole genome sequencing and DNA methylation profiling of hundreds of natural Arabidopsis accessions has revealed that transposable elements exhibit significant intraspecific genetic and epigenetic variation, and that genetic variation often underlies epigenetic variation. Together, epigenetic modification and the forces of selection define the scope within which transposable elements can contribute to, and control, genome evolution.

The double-edged sword of (re)expression of genes by hypomethylating agents: from viral mimicry to exploitation as priming agents for targeted immune checkpoint modulation

Hypomethylating agents (HMAs) have been widely used over the last decade, approved for use in myelodysplastic syndrome (MDS), chronic myelomonocytic leukemia (CMML) and acute myeloid leukemia (AML). The proposed central mechanism of action of HMAs, is the reversal of aberrant methylation in tumor cells, thus reactivating CpG-island promoters and leading to (re)expression of tumor suppressor genes. Recent investigations into the mode of action of azacitidine (AZA) and decitabine (DAC) have revealed new molecular mechanisms that impinge on tumor immunity via induction of an interferon response, through activation of endogenous retroviral elements (ERVs) that are normally epigenetically silenced. Although the global demethylation of DNA by HMAs can induce anti-tumor effects, it can also upregulate the expression of inhibitory immune checkpoint receptors and their ligands, resulting in secondary resistance to HMAs. Recent studies have, however, suggested that this could be exploited to prime or (re)sensitize tumors to immune checkpoint inhibitor therapies. In recent years, immune checkpoints have been targeted by novel therapies, with the aim of (re)activating the host immune system to specifically eliminate malignant cells. Antibodies blocking checkpoint receptors have been FDA-approved for some solid tumors and a plethora of clinical trials testing these and other checkpoint inhibitors are under way. This review will discuss AZA and DAC novel mechanisms of action resulting from the re-expression of pathologically hypermethylated promoters of gene sets that are related to interferon signaling, antigen presentation and inflammation. We also review new insights into the molecular mechanisms of action of transient, low-dose HMAs on various tumor types and discuss the potential of new treatment options and combinations.

Integration site selection by retroviruses and transposable elements in eukaryotes

Transposable elements and retroviruses are found in most genomes, can be pathogenic and are widely used as gene-delivery and functional genomics tools. Exploring whether these genetic elements target specific genomic sites for integration and how this preference is achieved is crucial to our understanding of genome evolution, somatic genome plasticity in cancer and ageing, host-parasite interactions and genome engineering applications. High-throughput profiling of integration sites by next-generation sequencing, combined with large-scale genomic data mining and cellular or biochemical approaches, has revealed that the insertions are usually non-random. The DNA sequence, chromatin and nuclear context, and cellular proteins cooperate in guiding integration in eukaryotic genomes, leading to a remarkable diversity of insertion site distribution and evolutionary strategies.

KRAB zinc-finger proteins contribute to the evolution of gene regulatory networks

The human genome encodes some 350 Krüppel-associated box (KRAB) domain-containing zinc-finger proteins (KZFPs), the products of a rapidly evolving gene family that has been traced back to early tetrapods. The function of most KZFPs is unknown, but a few have been demonstrated to repress transposable elements in embryonic stem (ES) cells by recruiting the transcriptional regulator TRIM28 and associated mediators of histone H3 Lys9 trimethylation (H3K9me3)-dependent heterochromatin formation and DNA methylation. Depletion of TRIM28 in human or mouse ES cells triggers the upregulation of a broad range of transposable elements, and recent data based on a few specific examples have pointed to an arms race between hosts and transposable elements as an important driver of KZFP gene selection. Here, to obtain a global view of this phenomenon, we combined phylogenetic and genomic studies to investigate the evolutionary emergence of KZFP genes in vertebrates and to identify their targets in the human genome. First, we unexpectedly reassigned the root of the family to a common ancestor of coelacanths and tetrapods. Second, although we confirmed that the majority of KZFPs bind transposable elements and pinpoint cases of ongoing co-evolution, we found that most of their transposable element targets have lost all transposition potential. Third, by examining the interplay between human KZFPs and other transcriptional modulators, we obtained evidence that KZFPs exploit evolutionarily conserved fragments of transposable elements as regulatory platforms long after the arms race against these genetic invaders has ended. Together, our results demonstrate that KZFPs partner with transposable elements to build a largely species-restricted layer of epigenetic regulation.

Ecological networks to unravel the routes to horizontal transposon transfers

Transposable elements (TEs) represent the single largest component of numerous eukaryotic genomes, and their activity and dispersal constitute an important force fostering evolutionary innovation. The horizontal transfer of TEs (HTT) between eukaryotic species is a common and widespread phenomenon that has had a profound impact on TE dynamics and, consequently, on the evolutionary trajectory of many species' lineages. However, the mechanisms promoting HTT remain largely unknown. In this article, we argue that network theory combined with functional ecology provides a robust conceptual framework and tools to delineate how complex interactions between diverse organisms may act in synergy to promote HTTs.

Exchange of Genetic Sequences Between Viruses and Hosts

Although genetic transfer between viruses and vertebrate hosts occurs less frequently than gene flow between bacteriophages and prokaryotes, it is extensive and has affected the evolution of both parties. With retroviruses, the integration of proviral DNA into chromosomal DNA can result in the activation of adjacent host gene expression and in the transduction of host transcripts into retroviral genomes as oncogenes. Yet in contrast to lysogenic phage, there is little evidence that viral oncogenes persist in a chain of natural transmission or that retroviral transduction is a significant driver of the horizontal spread of host genes. Conversely, integration of proviruses into the host germ line has generated endogenous retroviral genomes (ERV) in all vertebrate genomes sequenced to date. Some of these genomes retain potential infectivity and upon reactivation may transmit to other host species. During mammalian evolution, sequences of retroviral origin have been repurposed to serve host functions, such as the viral envelope glycoproteins crucial to the development of the placenta. Beyond retroviruses, DNA viruses with complex genomes have acquired numerous genes of host origin which influence replication, pathogenesis and immune evasion, while host species have accumulated germline sequences of both DNA and RNA viruses. A codicil is added on lateral transmission of cancer cells between hosts and on migration of host mitochondria into cancer cells.