Inter-Strain Epigenomic Profiling Reveals a Candidate IAP Master Copy in C3H Mice

Insertions of endogenous retroviruses cause a significant fraction of mutations in inbred mice but not all strains are equally susceptible. Notably, most new Intracisternal A particle (IAP) ERV mutagenic insertions have occurred in C3H mice. We show here that strain-specific insertional polymorphic IAPs accumulate faster in C3H/HeJ mice, relative to other sequenced strains, and that IAP transcript levels are higher in C3H/HeJ embryonic stem (ES) cells compared to other ES cells. To investigate the mechanism for high IAP activity in C3H mice, we identified 61 IAP copies in C3H/HeJ ES cells enriched with H3K4me3 (a mark of active promoters) and, among those tested, all are unmethylated in C3H/HeJ ES cells. Notably, 13 of the 61 are specific to C3H/HeJ and are members of the non-autonomous 1Δ1 IAP subfamily that is responsible for nearly all new insertions in C3H. One copy is full length with intact open reading frames and hence potentially capable of providing proteins in trans to other 1Δ1 elements. This potential “master copy” is present in other strains, including 129, but its 5’ long terminal repeat (LTR) is methylated in 129 ES cells. Thus, the unusual IAP activity in C3H may be due to reduced epigenetic repression coupled with the presence of a master copy.

Endogenous Retrovirus Transcript Levels Are Associated with Immunogenic Signatures in Multiple Metastatic Cancer Types

Next-generation sequencing of solid tumors has revealed variable signatures of immunogenicity across tumors, but underlying molecular characteristics driving such variation are not fully understood. Although expression of endogenous retrovirus (ERV)-containing transcripts can provide a source of tumor-specific neoantigen in some cancer models, associations between ERV levels and immunogenicity across different types of metastatic cancer are not well established. We performed bioinformatics analysis of genomic, transcriptomic, and clinical data across an integrated cohort of 199 patients with metastatic breast, colorectal, and pancreatic ductal adenocarcinoma tumors. Within each cancer type, we identified a subgroup of viral mimicry tumors in which increased ERV levels were coupled with transcriptional signatures of autonomous antiviral response and immunogenicity. In addition, viral mimicry colorectal and pancreatic tumors showed increased expression of DNA demethylation gene TET2 Taken together, these data demonstrate the existence of an ERV-associated viral mimicry phenotype across three distinct metastatic cancer types, while indicating links between ERV abundance, epigenetic dysregulation, and immunogenicity.

Abstract B56: Endogenous retrovirus transcript levels are associated with immunogenic signatures in multiple metastatic cancer types

Variability in the immunogenic landscape of metastatic cancer lesions has revealed insight into detection of potential immunotherapy-responsive patients, though underlying mechanisms driving such variation are not fully understood. Endogenous retrovirus (ERV)-containing transcripts have recently emerged as a potential source of tumor-associated antigen that are orthogonal to somatic mutation-derived neoantigens. To characterize the intersection between ERV levels and predicted immunogenicity in metastatic cancer, we comprehensively profiled the transcript abundance of 702,533 ERV loci in 199 metastatic tumors from breast, colorectal, and pancreatic cancer patients. In all three cancer types, overall ERV transcript load was associated with upregulation of genes involved in innate antiviral response pathways as well as genes involved in both adaptive and innate immune signaling. In colorectal and pancreatic tumors, samples with concomitant increases in ERV load and antiviral response gene expression, termed viral mimicry tumors, showed high expression of the DNA demethylation gene TET2, a gene previously described to promote transcription of ERV-containing transcripts. Collectively, these data are compatible with the notion that an increased level of ERV-containing transcripts may account for increased immunogenicity in a subset of metastatic tumors and support the relationship between DNA demethylation and ERV load in colorectal and pancreatic tumors.

Citation Format: James T. Topham, Emma Titmuss, Erin Pleasance, Laura M. Williamson, Joanna M. Karasinska, Luka Culibrk, Michael K.C. Lee, Steve E. Kalloger, Shehara Mendis, Richard A. Moore, Andrew J. Mungall, Janessa Laskin, Jonathan M. Loree, Dixie L. Mager, Marco A. Marra, Steven J.M. Jones, David F. Schaeffer, Daniel J. Renouf. Endogenous retrovirus transcript levels are associated with immunogenic signatures in multiple metastatic cancer types [abstract]. In: Proceedings of the AACR Special Conference on Pancreatic Cancer: Advances in Science and Clinical Care; 2019 Sept 6-9; Boston, MA. Philadelphia (PA): AACR; Cancer Res 2019;79(24 Suppl):Abstract nr B56.

Reality check for transposon enhancers

Hundreds of retrovirus-like sequences have features that suggest they might be gene enhancers, but only a small fraction displays gene-regulating activity in experiments on mouse stem cells.

Mouse germ line mutations due to retrotransposon insertions

Transposable element (TE) insertions are responsible for a significant fraction of spontaneous germ line mutations reported in inbred mouse strains. This major contribution of TEs to the mutational landscape in mouse contrasts with the situation in human, where their relative contribution as germ line insertional mutagens is much lower. In this focussed review, we provide comprehensive lists of TE-induced mouse mutations, discuss the different TE types involved in these insertional mutations and elaborate on particularly interesting cases. We also discuss differences and similarities between the mutational role of TEs in mice and humans.

LIONS: analysis suite for detecting and quantifying transposable element initiated transcription from RNA-seq

Summary

Transposable elements (TEs) influence the evolution of novel transcriptional networks yet the specific and meaningful interpretation of how TE-derived transcriptional initiation contributes to the transcriptome has been marred by computational and methodological deficiencies. We developed LIONS for the analysis of RNA-seq data to specifically detect and quantify TE-initiated transcripts.

Availability and implementation

Source code, container, test data and instruction manual are freely available at www.github.com/ababaian/LIONS.

Supplementary information

Supplementary data are available at Bioinformatics online.

Ten things you should know about transposable elements

Transposable elements (TEs) are major components of eukaryotic genomes. However, the extent of their impact on genome evolution, function, and disease remain a matter of intense interrogation. The rise of genomics and large-scale functional assays has shed new light on the multi-faceted activities of TEs and implies that they should no longer be marginalized. Here, we introduce the fundamental properties of TEs and their complex interactions with their cellular environment, which are crucial to understanding their impact and manifold consequences for organismal biology. While we draw examples primarily from mammalian systems, the core concepts outlined here are relevant to a broad range of organisms.

A novel isoform of IL-33 revealed by screening for transposable element promoted genes in human colorectal cancer

Remnants of ancient transposable elements (TEs) are abundant in mammalian genomes. These sequences contain multiple regulatory motifs and hence are capable of influencing expression of host genes. TEs are known to be released from epigenetic repression and can become transcriptionally active in cancer. Such activation could also lead to lineage-inappropriate activation of oncogenes, as previously described in lymphomas. However, there are few reports of this mechanism occurring in non-blood cancers. Here, we re-analyzed whole transcriptome data from a large cohort of patients with colon cancer, compared to matched normal colon control samples, to detect genes or transcripts ectopically expressed through activation of TE promoters. Among many such transcripts, we identified six where the affected gene has described role in cancer and where the TE-driven gene mRNA is expressed in primary colon cancer, but not normal matched tissue, and confirmed expression in colon cancer-derived cell lines. We further characterized a TE-gene chimeric transcript involving the Interleukin 33 (IL-33) gene (termed LTR-IL-33), that is ectopically expressed in a subset of colon cancer samples through the use of an endogenous retroviral long terminal repeat (LTR) promoter of the MSTD family. The LTR-IL-33 chimeric transcript encodes a novel shorter isoform of the protein, which is missing the initial N-terminus (including many conserved residues) of Native IL-33. In vitro studies showed that LTR-IL-33 expression is required for optimal CRC cell line growth as 3D colonospheres. Taken together, these data demonstrate the significance of TEs as regulators of aberrant gene expression in colon cancer.

Induction of endogenous retroelements as a potential mechanism for mouse-specific drug-induced carcinogenicity

A number of chemical compounds have been shown to induce liver tumors in mice but not in other species. While several mechanisms for this species-specific tumorigenicity have been proposed, no definitive mechanism has been established. We examined the effects of the nongenotoxic rodent hepatic carcinogen, WY-14,643, in male mice from a high liver tumor susceptible strain (C3H/HeJ), and from a low tumor susceptible strain (C57BL/6). WY-14,643, a PPARα activator induced widespread increases in the expression of some endogenous retroelements, namely members of LTR and LINE elements in both strains. The expression of a number of known retroviral defense genes was also elevated. We also demonstrated that basal immune-mediated viral defense was elevated in C57BL/6 mice (the resistant strain) and that WY-14,643 further activated those immuno-defense processes. We propose that the previously reported >100X activity of retroelements in mice drives mouse-specific tumorigenicity. We also propose that C57BL/6's competent immune to retroviral activation allows it to remove cells before the activation of these elements can result in significant chromosomal insertions and mutation. Finally, we showed that WY-14,643 treatment induced gene signatures of DNA recombination in the sensitive C3H/HeJ strain.

Pluripotency and the endogenous retrovirus HERVH: Conflict or serendipity?

Remnants of ancient retroviral infections during evolution litter all mammalian genomes. In modern humans, such endogenous retroviral (ERV) sequences comprise at least 8% of the genome. While ERVs and other types of transposable elements undoubtedly contribute to the genomic "junk yard", functions for some ERV sequences have been demonstrated, with growing evidence that ERVs can be important players in gene regulatory processes. Here we focus on one particular large family of human ERVs, termed HERVH, which several recent studies suggest has a key regulatory role in human pluripotent stem cells. Remarkably, this is not the first instance of an ERV controlling pluripotency. We speculate as to why this convergent evolution might have come about, suggesting that it may reflect selection on the virus to extend the time available for transposition. Alternatively it may reflect serendipity alone.

Endogenous retroviral promoter exaptation in human cancer

Cancer arises from a series of genetic and epigenetic changes, which result in abnormal expression or mutational activation of oncogenes, as well as suppression/inactivation of tumor suppressor genes. Aberrant expression of coding genes or long non-coding RNAs (lncRNAs) with oncogenic properties can be caused by translocations, gene amplifications, point mutations or other less characterized mechanisms. One such mechanism is the inappropriate usage of normally dormant, tissue-restricted or cryptic enhancers or promoters that serve to drive oncogenic gene expression. Dispersed across the human genome, endogenous retroviruses (ERVs) provide an enormous reservoir of autonomous gene regulatory modules, some of which have been co-opted by the host during evolution to play important roles in normal regulation of genes and gene networks. This review focuses on the “dark side” of such ERV regulatory capacity. Specifically, we discuss a growing number of examples of normally dormant or epigenetically repressed ERVs that have been harnessed to drive oncogenes in human cancer, a process we term onco-exaptation, and we propose potential mechanisms that may underlie this phenomenon.

Transposable elements in the mammalian embryo: pioneers surviving through stealth and service

Transposable elements (TEs) are notable drivers of genetic innovation. Over evolutionary time, TE insertions can supply new promoter, enhancer, and insulator elements to protein-coding genes and establish novel, species-specific gene regulatory networks. Conversely, ongoing TE-driven insertional mutagenesis, nonhomologous recombination, and other potentially deleterious processes can cause sporadic disease by disrupting genome integrity or inducing abrupt gene expression changes. Here, we discuss recent evidence suggesting that TEs may contribute regulatory innovation to mammalian embryonic and pluripotent states as a means to ward off complete repression by their host genome.

Methylated DNA Immunoprecipitation Analysis of Mammalian Endogenous Retroviruses

Endogenous retroviruses are repetitive sequences found abundantly in mammalian genomes which are capable of modulating host gene expression. Nevertheless, most endogenous retrovirus copies are under tight epigenetic control via histone-repressive modifications and DNA methylation. Here we describe a common method used in our laboratory to detect, quantify, and compare mammalian endogenous retrovirus DNA methylation. More specifically we describe methylated DNA immunoprecipitation (MeDIP) followed by quantitative PCR.

Genistein versus ICI 182, 780: an ally or enemy in metastatic progression of prostate cancer

BACKGROUND

Androgen signalling through the androgen receptor (AR) plays a critical role in prostate cancer (PCa) initiation and progression. Estrogen in synergy with androgen is essential for cell growth of the normal and malignant prostate. However, the exact role that estrogen and the estrogen receptor play in prostate carcinogenesis remains unclear. We have previously demonstrated the metastasis-promoting effect of an estrogen receptor beta (ERβ) agonist (genistein) in a patient-derived PCa xenograft model mimicking localized and metastatic disease.

METHODS

To test the hypothesis that the tumor-promoting activity of genistein was due to its estrogenic properties, we treated the xenograft-bearing mice with genistein and an anti-estrogen compound (ICI 182, 780) and compared the differential gene expression using microarrays.

RESULTS

Using a second xenograft model which was derived from another patient, we showed that genistein promoted disease progression in vivo and ICI 182, 780 inhibited metastatic spread. The microarray analysis revealed that the metallothionein (MT) gene family was differentially expressed in tumors treated by these compounds. Using qRT-PCR, the differences in expression levels were validated in the metastatic and non-metastatic LTL313 PCa xenograft tumor lines, both of which were originally derived from the same PCa patient.

CONCLUSIONS

Together our data provide evidence that genistein stimulates and ICI 182, 780 inhibits metastatic progression, suggesting that these effects may be mediated by ERβ signalling.

Epigenetic interplay between mouse endogenous retroviruses and host genes

Background

Transposable elements are often the targets of repressive epigenetic modifications such as DNA methylation that, in theory, have the potential to spread toward nearby genes and induce epigenetic silencing. To better understand the role of DNA methylation in the relationship between transposable elements and genes, we assessed the methylation state of mouse endogenous retroviruses (ERVs) located near genes.

Results

We found that ERVs of the ETn/MusD family show decreased DNA methylation when near transcription start sites in tissues where the nearby gene is expressed. ERVs belonging to the IAP family, however, are generally heavily methylated, regardless of the genomic environment and the tissue studied. Furthermore, we found full-length ETn and IAP copies that display differential DNA methylation between their two long terminal repeats (LTRs), suggesting that the environment surrounding gene promoters can prevent methylation of the nearby LTR. Spreading from methylated ERV copies to nearby genes was rarely observed, with the regions between the ERVs and genes apparently acting as a boundary, enriched in H3K4me3 and CTCF, which possibly protects the unmethylated gene promoter. Furthermore, the flanking regions of unmethylated ERV copies harbor H3K4me3, consistent with spreading of euchromatin from the host gene toward ERV insertions.

Conclusions

We have shown that spreading of DNA methylation from ERV copies toward active gene promoters is rare. We provide evidence that genes can be protected from ERV-induced heterochromatin spreading by either blocking the invasion of repressive marks or by spreading euchromatin toward the ERV copy.

Transposable elements: an abundant and natural source of regulatory sequences for host genes

The fact that transposable elements (TEs) can influence host gene expression was first recognized more than 50 years ago. However, since that time, TEs have been widely regarded as harmful genetic parasites-selfish elements that are rarely co-opted by the genome to serve a beneficial role. Here, we survey recent findings that relate to TE impact on host genes and remind the reader that TEs, in contrast to other noncoding parts of the genome, are uniquely suited to gene regulatory functions. We review recent studies that demonstrate the role of TEs in establishing and rewiring gene regulatory networks and discuss the overall ubiquity of exaptation. We suggest that although individuals within a population can be harmed by the deleterious effects of new TE insertions, the presence of TE sequences in a genome is of overall benefit to the population.

Transposable elements: not as quiet as a mouse

In this issue of Genome Biology, Nellåker et al. show massive purging of deleterious transposable element variants, through negative selection, in 18 mouse strains.

C-GATE - catalogue of genes affected by transposable elements

Background

Functional regulatory sequences are present in many transposable element (TE) copies, resulting in TEs being frequently exapted by host genes. Today, many examples of TEs impacting host gene expression can be found in the literature and we believe a new catalogue of such exaptations would be useful for the field.

Findings

We have established the catalogue of genes affected by transposable elements (C-GATE), which can be found at https://sites.google.com/site/tecatalog/. To date, it holds 221 cases of biologically verified TE exaptations and more than 10,000 in silico TE-gene partnerships. C-GATE is interactive and allows users to include missed or new TE exaptation data. C-GATE provides a graphic representation of the entire library, which may be used for future statistical analysis of TE impact on host gene expression.

Conclusions

We hope C-GATE will be valuable for the TE community but also for others who have realized the role that TEs may have in their research.

Gene properties and chromatin state influence the accumulation of transposable elements in genes

Transposable elements (TEs) are mobile DNA sequences found in the genomes of almost all species. By measuring the normalized coverage of TE sequences within genes, we identified sets of genes with conserved extremes of high/low TE density in the genomes of human, mouse and cow and denoted them as ‘shared upper/lower outliers (SUOs/SLOs)’. By comparing these outlier genes to the genomic background, we show that a large proportion of SUOs are involved in metabolic pathways and tend to be mammal-specific, whereas many SLOs are related to developmental processes and have more ancient origins. Furthermore, the proportions of different types of TEs within human and mouse orthologous SUOs showed high similarity, even though most detectable TEs in these two genomes inserted after their divergence. Interestingly, our computational analysis of polymerase-II (Pol-II) occupancy at gene promoters in different mouse tissues showed that 60% of tissue-specific SUOs show strong Pol-II binding only in embryonic stem cells (ESCs), a proportion significantly higher than the genomic background (37%). In addition, our analysis of histone marks such as H3K4me3 and H3K27me3 in mouse ESCs also suggest a strong association between TE-rich genes and open-chromatin at promoters. Finally, two independent whole-transcriptome datasets show a positive association between TE density and gene expression level in ESCs. While this study focuses on genes with extreme TE densities, the above results clearly show that the probability of TE accumulation/fixation in mammalian genes is not random and is likely associated with different factors/gene properties and, most importantly, an association between the TE insertion/fixation rate and gene activity status in ES cells.

Role of runt-related transcription factor 3 (RUNX3) in transcription regulation of natural cytotoxicity receptor 1 (NCR1/NKp46), an activating natural killer (NK) cell receptor

Natural cytotoxicity receptor 1 (NCR1), also known as NKp46, is a natural killer (NK) lymphocyte-activating receptor. It is involved in major aspects of NK immune function and shows a high degree of lineage specificity in blood and bone marrow. The nature of its NK-restricted expression is not well understood. In this study, we confirm that human NCR1 NK-specific expression is achieved at the mRNA level. We found two key cis-regulatory elements in the immediate vicinity upstream of the gene. One element acts as an essential promoter, whereas the other acts as a tissue-dependent enhancer/repressor. This latter regulatory element contains a runt related-transcription factor (RUNX) recognition motif that preferentially binds RUNX3. Interfering with RUNX proteins using a dominant negative form results in decreased Ncr1 expression. RUNX3 overexpression had the opposite effect. These findings shed light on the role of RUNX3 in the control of an important NK-activating receptor.

Human Th1 and Th17 cells exhibit epigenetic stability at signature cytokine and transcription factor loci

The linear model of Th cell lineage commitment is being revised due to reports that mature Th cells can trans-differentiate into alternate lineages. This ability of Th cells to reprogram is thought to be regulated by epigenetic mechanisms that control expression of transcription factors characteristic of opposing lineages. It is unclear, however, to what extent this new model of Th cell plasticity holds true in human Th cell subsets that develop under physiological conditions in vivo. We isolated in vivo-differentiated human Th1 and Th17 cells, as well as intermediate Th1/17 cells, and identified distinct epigenetic signatures at cytokine (IFNG and IL17A) and transcription factor (TBX21, RORC, and RORA) loci. We also examined the phenotypic and epigenetic stability of human Th17 cells exposed to Th1-polarizing conditions and found that although they could upregulate TBX21 and IFN-γ, this occurred without loss of IL-17 or RORC expression, and resulted in cells with a Th1/17 phenotype. Similarly, Th1 cells could upregulate IL-17 upon enforced expression of RORC2, but did not lose expression of IFN-γ or TBX21. Despite alterations in expression of these signature genes, epigenetic modifications were remarkably stable aside from the acquisition of active histone methylation marks at cytokine gene promoters. The limited capacity of human Th17 and Th1 cells to undergo complete lineage conversion suggests that the bipotent Th1/17 cells may arise from Th1 and/or Th17 cells. These data also question the broad applicability of the new model of Th cell lineage plasticity to in vivo-polarized human Th cell subsets.

Placenta-specific expression of the interleukin-2 (IL-2) receptor β subunit from an endogenous retroviral promoter

The long terminal repeat (LTR) sequences of endogenous retroviruses and retroelements contain promoter elements and are known to form chimeric transcripts with nearby cellular genes. Here we show that an LTR of the THE1D retroelement family has been domesticated as an alternative promoter of human IL2RB, the gene encoding the β subunit of the IL-2 receptor. The LTR promoter confers expression specifically in the placental trophoblast as opposed to its native transcription in the hematopoietic system. Rather than sequence-specific determinants, DNA methylation was found to regulate transcription initiation and splicing efficiency in a tissue-specific manner. Furthermore, we detected the cytoplasmic signaling domain of the IL-2Rβ protein in the placenta, suggesting that IL-2Rβ undergoes preferential proteolytic cleavage in this tissue. These findings implicate novel functions for this cytokine receptor subunit in the villous trophoblast and reveal an intriguing example of ancient LTR exaptation to drive tissue-specific gene expression.

H3K9me3-binding proteins are dispensable for SETDB1/H3K9me3-dependent retroviral silencing

Background

Endogenous retroviruses (ERVs) are parasitic sequences whose derepression is associated with cancer and genomic instability. Many ERV families are silenced in mouse embryonic stem cells (mESCs) via SETDB1-deposited trimethylated lysine 9 of histone 3 (H3K9me3), but the mechanism of H3K9me3-dependent repression remains unknown. Multiple proteins, including members of the heterochromatin protein 1 (HP1) family, bind H3K9me2/3 and are involved in transcriptional silencing in model organisms. In this work, we address the role of such H3K9me2/3 "readers" in the silencing of ERVs in mESCs.

Results

We demonstrate that despite the reported function of HP1 proteins in H3K9me-dependent gene repression and the critical role of H3K9me3 in transcriptional silencing of class I and class II ERVs, the depletion of HP1α, HP1β and HP1γ, alone or in combination, is not sufficient for derepression of these elements in mESCs. While loss of HP1α or HP1β leads to modest defects in DNA methylation of ERVs or spreading of H4K20me3 into flanking genomic sequence, respectively, neither protein affects H3K9me3 or H4K20me3 in ERV bodies. Furthermore, using novel ERV reporter constructs targeted to a specific genomic site, we demonstrate that, relative to Setdb1, knockdown of the remaining known H3K9me3 readers expressed in mESCs, including Cdyl, Cdyl2, Cbx2, Cbx7, Mpp8, Uhrf1 and Jarid1a-c, leads to only modest proviral reactivation.

Conclusion

Taken together, these results reveal that each of the known H3K9me3-binding proteins is dispensable for SETDB1-mediated ERV silencing. We speculate that H3K9me3 might maintain ERVs in a silent state in mESCs by directly inhibiting deposition of active covalent histone marks.

DNA methylation and SETDB1/H3K9me3 regulate predominantly distinct sets of genes, retroelements, and chimeric transcripts in mESCs

DNA methylation and histone H3 lysine 9 trimethylation (H3K9me3) play important roles in silencing of genes and retroelements. However, a comprehensive comparison of genes and repetitive elements repressed by these pathways has not been reported. Here we show that in mouse embryonic stem cells (mESCs), the genes upregulated after deletion of the H3K9 methyltransferase Setdb1 are distinct from those derepressed in mESC deficient in the DNA methyltransferases Dnmt1, Dnmt3a, and Dnmt3b, with the exception of a small number of primarily germline-specific genes. Numerous endogenous retroviruses (ERVs) lose H3K9me3 and are concomitantly derepressed exclusively in SETDB1 knockout mESCs. Strikingly, ~15% of upregulated genes are induced in association with derepression of promoter-proximal ERVs, half in the context of "chimeric" transcripts that initiate within these retroelements and splice to genic exons. Thus, SETDB1 plays a previously unappreciated yet critical role in inhibiting aberrant gene transcription by suppressing the expression of proximal ERVs.

Distributions of transposable elements reveal hazardous zones in mammalian introns

Comprising nearly half of the human and mouse genomes, transposable elements (TEs) are found within most genes. Although the vast majority of TEs in introns are fixed in the species and presumably exert no significant effects on the enclosing gene, some markedly perturb transcription and result in disease or a mutated phenotype. Factors determining the likelihood that an intronic TE will affect transcription are not clear. In this study, we examined intronic TE distributions in both human and mouse and found several factors that likely contribute to whether a particular TE can influence gene transcription. Specifically, we observed that TEs near exons are greatly underrepresented compared to random distributions, but the size of these “underrepresentation zones” differs between TE classes. Compared to elsewhere in introns, TEs within these zones are shorter on average and show stronger orientation biases. Moreover, TEs in extremely close proximity (<20 bp) to exons show a strong bias to be near splice-donor sites. Interestingly, disease-causing intronic TE insertions show the opposite distributional trends, and by examining expressed sequence tag (EST) databases, we found that the proportion of TEs contributing to chimeric TE-gene transcripts is significantly higher within their underrepresentation zones. In addition, an analysis of predicted splice sites within human long terminal repeat (LTR) elements showed a significantly lower total number and weaker strength for intronic LTRs near exons. Based on these factors, we selectively examined a list of polymorphic mouse LTR elements in introns and showed clear evidence of transcriptional disruption by LTR element insertions in the Trpc6 and Kcnh6 genes. Taken together, these studies lend insight into the potential selective forces that have shaped intronic TE distributions and enable identification of TEs most likely to exert transcriptional effects on genes.

Sequences derived from transposable elements (TEs) are major constituents of mammalian genomes and are found within introns of most genes. While nearly all TEs within introns appear harmless, some de novo intronic TE insertions do disrupt gene transcription and splicing and cause disease. It is unclear why some intronic TEs perturb gene transcription whereas most do not. Here, we examined intronic TE distributions in both human and mouse genes to gain insight into which TEs may be more likely to affect transcription. We found evidence that TEs near exons are likely subject to strong negative selection but the size of the region under selection or “underrepresentation zone” differs for different TE classes. Strikingly, all reported human disease-causing intronic TE insertions fall within these underrepresentation zones, and the proportion of TEs contributing to chimeric TE-gene transcripts is significantly higher when TEs are located in these zones. We also examined insertionally polymorphic mouse TEs located within underrepresentation zones and found evidence of transcriptional disruption in two genes. Given the growing appreciation for ongoing activity of TEs in human, our results should be of value in prioritizing insertionally polymorphic TEs for study of their potential contributions to gene expression differences and phenotypic variability.

The clf2 gene has an epigenetic role in the multifactorial etiology of cleft lip and palate in the A/WySn mouse strain

BACKGROUND

The A/WySn mouse strain with 15 to 20% penetrance of cleft lip and palate (CLP) is an animal model for human multifactorial CLP. The CLP is due to two unlinked genes that interact epistatically, Wnt9b(clf1) and clf2, plus a maternal effect. The Wnt9b(clf1) mutation is an IAP transposon insertion. The clf2 gene, with unknown function, was located in a 13.6 Mb region of chromosome 13 containing 145 genes.

METHODS

To reduce the clf2 candidate region, 1146 mice segregating for A/WySn and C57BL/6J alleles at clf2 were screened for recombinants by simple sequence-length polymorphism haplotypes; recombinants' testcross progeny were typed for CLP and simple-sequence length polymorphisms. To identify the function of clf2, the effect of clf2 genotype on risk of CLP was tested in Wnt9b(null/null) knockouts and in compound mutants (Wnt9b(clf1/null) ), and the methylation of the IAP at Wnt9b was assayed in the Wnt9b(clf1/null) mutants by combined bisulfite restriction analysis.

RESULTS

The location of clf2 was redefined to 3.0 Mb between Cntnap3 and AK029746 containing 48 genes, of which 30 are Zfp genes. The clf2 genotype had no detectable effect on Wnt9b(null/null) embryos, but strongly affected risk of CLP and methylation of the IAP in Wnt9b(clf1/null) embryos. CLP was associated with low levels of methylation of the IAP.

CONCLUSIONS

The clf2 gene is the first identified polymorphism that affects the epigenetic methylation and silencing of IAP retrotransposons. This CLP model raises the question of whether parallel epigenetic factors are involved in risk and environmental sensitivity of human CLP.

Potential mechanisms of endogenous retroviral-mediated genomic instability in human cancer

Malignancy results from a complex combination of genetic and epigenetic changes, the full effects of which are still largely unknown. Here we summarize current knowledge of the origin, retrotranspositional activity, epigenetic state, and transcription of human endogenous retroviruses (HERVs), and then discuss the potential effects of their deregulation in cancer. Evidence suggests that cancer-associated epigenetic changes most likely underlie potential HERV-mediated effects on genome and transcriptome instability and may play a role in malignancy. Despite our currently limited understanding of the importance of HERVs or other transposable elements in cancer development, we believe that the emerging era of high-throughput sequencing of cancer genomes, epigenomes, and transcriptomes will provide unprecedented opportunities to investigate these roles in the future.

Expression of the leukemic prognostic marker CD7 is linked to epigenetic modifications in chronic myeloid leukemia

Background

Expression levels of the cell surface glycoprotein, CD7, and the serine protease, elastase 2 (ELA2), in the leukemic cells of patients with chronic myeloid leukemia (CML) have been associated with clinical outcome. However, little is known about the mechanisms that underlie the variable expression of these genes in the leukemic cells.

Results

To address this question, we compared the level of their expression with the DNA methylation and histone acetylation status of 5' sequences of both genes in leukemic cell lines and primitive (lin^-CD34⁺) leukemic cells from chronic phase CML patients. DNA methylation of the ELA2 gene promoter did not correlate with its expression pattern in lin^-CD34⁺ cells from chronic phase CML patient samples even though there was clear differential DNA methylation of this locus in ELA2-expressing and non-expressing cell lines. In contrast, we found a strong relation between CD7 expression and transcription-permissive chromatin modifications, both at the level of DNA methylation and histone acetylation with evidence of hypomethylation of the CD7 promoter region in the lin^-CD34⁺ cells from CML patients with high CD7 expression.

Conclusion

These findings indicate a link between epigenetic modifications and CD7 expression in primitive CML cells.

Variable DNA methylation of transposable elements: the case study of mouse Early Transposons

Phenotypic variation stems from both genetic and epigenetic differences between individuals. In order to elucidate how phenotypes are determined, it is necessary to understand the forces that generate variation in genome sequence as well as its epigenetic state. In both contexts, transposable elements (TEs) may play an important role. It is well established that TE activity is a major generator of genetic variation, but recent research also suggests that TEs contribute to epigenetic variation. Stochastic epigenetic silencing of some TE insertions in mice has been shown to cause phenotypic variability between individuals. However, the prevalence of this phenomenon has never been evaluated. Here, we use 18 insertions of a mouse Endogenous Retrovirus (ERV) family, the Early Transposons (ETns), to detect insertion-dependent determinants of DNA methylation levels and variability between both cells and individuals. We show that the structure and age of insertions influence methylation levels and variability, resulting in a subgroup of loci that displays unexpectedly high variability in methylation and suggesting stochastic events during methylation establishment. Despite variation in methylation according to the age and structure of each locus, homologous CpG sites show similar tendencies in methylation levels across loci, emphasizing the role of the insertion's sequence in methylation determination. Our results show that differences in methylation of ETns between individuals is not a sporadic phenomenon and support the hypothesis that ERVs contribute to phenotypic variability through their stochastic silencing.

Evidence for high bi-allelic expression of activating Ly49 receptors

Stochastic expression is a hallmark of the Ly49 family that encode the main MHC class-I-recognizing receptors of mouse natural killer (NK) cells. This highly polygenic and polymorphic family includes both activating and inhibitory receptor genes and is one of genome's fastest evolving loci. The inhibitory Ly49 genes are expressed in a stochastic mono-allelic manner, possibly under the control of an upstream bi-directional early promoter and show mono-allelic DNA methylation patterns. To date, no studies have directly addressed the transcriptional regulation of the activating Ly49 receptors. Our study shows differences in DNA methylation pattern between activating and inhibitory genes in C57BL/6 and F1 hybrid mouse strains. We also show a bias towards bi-allelic expression of the activating receptors based on allele-specific single-cell RT–PCR in F1 hybrid NK cells for Ly49d and Ly49H expression in Ly49h^+/− mice. Furthermore, we have identified a region of high sequence identity with possible transcriptional regulatory capacity for the activating Ly49 genes. Our results also point to a likely difference between NK and T-cells in their ability to transcribe the activating Ly49 genes. These studies highlight the complex regulation of this rapidly evolving gene family of central importance in mouse NK cell function.

Isolation of human iPS cells using EOS lentiviral vectors to select for pluripotency

Induced pluripotent stem (iPS) cells may be of use in regenerative medicine. However, the low efficiency of reprogramming is a major impediment to the generation of patient-specific iPS cell lines. Here we report the first selection system for the isolation of human iPS cells. We developed the EOS (Early Transposon promoter and Oct-4 (Pou5f1) and Sox2 enhancers) lentiviral vector to specifically express in mouse and human embryonic stem cells but not in primary fibroblasts. The bicistronic EOS vector marked emerging mouse and human iPS cell colonies with EGFP, and we used puromycin selection to aid the isolation of iPS cell lines that expressed endogenous pluripotency markers. These lines differentiated into cell types from all three germ layers. Reporter expression was extinguished upon differentiation and therefore monitored the residual pluripotent cells that form teratomas. Finally, we used EOS selection to establish Rett syndrome-specific mouse and human iPS cell lines with known mutations in MECP2.

The promoter for the long isoform of mouse NKG2D is a B1 retrotransposon with a functional GABP site (136.24)

The mouse gene for the natural killer (NK) cell activating receptor NKG2D produces two protein isoforms, NKG2D-S and NKG2D-L, which differ by 13 amino acids at the N-terminus and have different signaling capabilities. These two isoforms are produced through the use of alternative promoters but the regulation of these promoters has not been investigated. In this study, we show that rat NKG2D has the same splicing pattern as that of the mouse and we mapped transcriptional start sites in both species. We also show that the NKG2D-L specific promoter, active in both resting and stimulated NK cells, is derived from a rodent B1 retrotransposon that inserted before mouse-rat divergence. This B1 insertion is associated with loss of a nearby splice acceptor site that subsequently allowed creation of the short NKG2D isoform found in mouse but not human. Reporter assays suggest that the B1 promoter has no inherent lymphoid tissue-specificity. Finally, we have identified different binding sites for the ETS family member GABP within both the mouse and rat B1 elements that are necessary for high promoter activity.

This work was supported by the Canadian Institutes of Health Research and the Cancer Research Society.

Creation of the two isoforms of rodent NKG2D was driven by a B1 retrotransposon insertion

The mouse gene for the natural killer (NK) cell-activating receptor Nkg2d produces two protein isoforms, NKG2D-S and NKG2D-L, which differ by 13 amino acids at the N-terminus and have different signalling capabilities. These two isoforms are produced through differential splicing, but their regulation has not been investigated. In this study, we show that rat Nkg2d has the same splicing pattern as that of the mouse, and we mapped transcriptional start sites in both species. We found that the splice forms arise from alternative promoters and that the NKG2D-L promoter is derived from a rodent B1 retrotransposon that inserted before mouse–rat divergence. This B1 insertion is associated with loss of a nearby splice acceptor site that subsequently allowed creation of the short NKG2D isoform found in mouse but not human. Transient reporter assays indicate that the B1 element is a strong promoter with no inherent lymphoid tissue-specificity. We have also identified different binding sites for the ETS family member GABP within both the mouse and rat B1 elements that are necessary for high-promoter activity and for full Nkg2d-L expression. These findings demonstrate that a retroelement insertion has led to gene-regulatory change and functional diversification of rodent NKG2D.

Keeping active endogenous retroviral-like elements in check: the epigenetic perspective

Endogenous retrovirus-like elements, or ERVs, are an abundant component of all eukaryotic genomes. Their transcriptional and retrotranspositional activities have great potential for deleterious effects on gene expression. Consequences of such activity may include germline mutagenesis and cancerous transformation. As a result, mammalian genomes have evolved means of counteracting ERV transcription and mobilization. In this review, we discuss epigenetic mechanisms of ERV and LTR retrotransposon control during mouse development, focusing on involvement of DNA methylation, histone modifications, small RNAs and their interaction with one another. We also address relevance of research performed in the mouse system to human and challenges associated with studying repetitive families. (Part of a multi-author review).

DNA methylation in ES cells requires the lysine methyltransferase G9a but not its catalytic activity

Histone H3K9 methylation is required for DNA methylation and silencing of repetitive elements in plants and filamentous fungi. In mammalian cells however, deletion of the H3K9 histone methyltransferases (HMTases) Suv39h1 and Suv39h2 does not affect DNA methylation of the endogenous retrovirus murine leukaemia virus, indicating that H3K9 methylation is dispensable for DNA methylation of retrotransposons, or that a different HMTase is involved. We demonstrate that embryonic stem (ES) cells lacking the H3K9 HMTase G9a show a significant reduction in DNA methylation of retrotransposons, major satellite repeats and densely methylated CpG-rich promoters. Surprisingly, demethylated retrotransposons remain transcriptionally silent in G9a(-/-) cells, and show only a modest decrease in H3K9me2 and no decrease in H3K9me3 or HP1alpha binding, indicating that H3K9 methylation per se is not the relevant trigger for DNA methylation. Indeed, introduction of catalytically inactive G9a transgenes partially 'rescues' the DNA methylation defect observed in G9a(-/-) cells. Taken together, these observations reveal that H3K9me3 and HP1alpha recruitment to retrotransposons occurs independent of DNA methylation in ES cells and that G9a promotes DNA methylation independent of its HMTase activity.

MiRNAs, epigenetics, and cancer

By virtue of having multiple targets, a microRNA (miRNA) can have variable effects on oncogenesis by acting as tumor suppressor or oncogene in a context-dependent manner. Genome-wide epigenetic changes that occur in various cancers affect the transcription of many genes. Since the transcriptional regulation of miRNAs remains an unexplored field, it is still unknown how epigenetic changes will affect the regulation of miRNAs. Many miRNAs are intron-bound within the body of a protein-coding gene. Any change to the transcription of the "host" gene affects the transcription and genesis of the resident miRNA. It is therefore reasonable to deduce that epigenetic changes brought on by transformation can potentially affect miRNA expression in both direct and indirect ways. We have reviewed the literature pertaining to the epigenetic regulation of miRNA genes in the context of various cancers and have speculated on the potential role of epigenetic modifications on the transcriptional regulation and expression of these genes.

Widely variable endogenous retroviral methylation levels in human placenta

It is generally assumed that transposable elements, including endogenous retroviruses (ERVs), are silenced by DNA methylation/chromatin structure in mammalian cells. However, there have been very few experimental studies to examine the methylation status of human ERVs. In this study, we determined and compared the methylation status of the 5′ long terminal repeats (LTRs) of different copies of the human endogenous retrovirus (HERV) family HERV-E, which are inserted in various genomic contexts. We found that three HERV-E LTRs which function as alternative gene promoters in placenta are unmethylated in that tissue but heavily methylated in blood cells, where these LTRs are not active promoters. This difference is not solely due to global hypomethylation in placenta, since two general measures of methylation levels of HERV-E and HERV-K LTRs suggest only 10–15% lower overall HERV methylation in placenta compared to blood. Comparisons between methylation levels of the LTR-derived gene promoters and six random HERV-E LTRs in placenta showed that the former display significantly lower methylation levels than random LTRs. Moreover, the differences in methylation between LTRs cannot always be explained by their genomic environment, since methylation of flanking sequences can be very different from methylation of the LTR itself.

Expression of murine killer immunoglobulin-like receptor KIRL1 on CD1d-independent NK1.1(+) T cells

Three mouse killer immunoglobulin-like receptors (KIRs), namely, KIR3DL1, KIRL1, and KIRL2, have recently been identified in C56BL/6 (B6) mice. However, only two Kir genes are found in the B6 mouse genome sequence data base. To clarify this discrepancy, we cloned Kir cDNAs from multiple strains of mice. Sequencing of the cDNA clones showed that the Kir3dl1 gene is found in C3H/HeJ and CBA/J but not in B6 mice. Analysis of the single nucleotide polymorphism data base suggested that Kir3dl1 is the C3H/HeJ and CBA/J allele of Kirl1. We generated mAb to the recombinant KIRL1 protein to investigate its expression pattern. The anti-KIRL1 mAb bound to NK1.1(+) T cells but only very weakly or at undetectable levels to other lymphocytes including natural killer (NK) cells and conventional T cells. Among NK1.1(+) T cells, conventional NK T cells stained with CD1d tetramer did not significantly bind anti-KIRL1 mAb, whereas CD1d-tetramer-negative subset was KIRL1-positive. Furthermore, the expression of KIRL1 is readily detected on NK1.1(+) T cells from beta(2)-microglobulin-deficient B6 mice. Thus, KIRL1 is predominantly expressed on CD1d-independent NK1.1(+) T cells.

Multiple effects govern endogenous retrovirus survival patterns in human gene introns

BACKGROUND

Endogenous retroviruses (ERVs) and solitary long terminal repeats (LTRs) have a significant antisense bias when located in gene introns, suggesting strong negative selective pressure on such elements oriented in the same transcriptional direction as the enclosing gene. It has been assumed that this bias reflects the presence of strong transcriptional regulatory signals within LTRs but little work has been done to investigate this phenomenon further.

RESULTS

In the analysis reported here, we found significant differences between individual human ERV families in their prevalence within genes and degree of antisense bias and show that, regardless of orientation, ERVs of most families are less likely to be found in introns than in intergenic regions. Examination of density profiles of ERVs across transcriptional units and the transcription signals present in the consensus ERVs suggests the importance of splice acceptor sites, in conjunction with splice donor and polyadenylation signals, as the major targets for selection against most families of ERVs/LTRs. Furthermore, analysis of annotated human mRNA splicing events involving ERV sequence revealed that the relatively young human ERVs (HERVs), HERV9 and HERV-K (HML-2), are involved in no human mRNA splicing events at all when oriented antisense to gene transcription, while elements in the sense direction in transcribed regions show considerable bias for use of strong splice sites.

CONCLUSION

Our observations suggest suppression of splicing among young intronic ERVs oriented antisense to gene transcription, which may account for their reduced mutagenicity and higher fixation rate in gene introns.

Repeated recruitment of LTR retrotransposons as promoters by the anti-apoptotic locus NAIP during mammalian evolution

Neuronal apoptosis inhibitory protein (NAIP, also known as BIRC1) is a member of the conserved inhibitor of apoptosis protein (IAP) family. Lineage-specific rearrangements and expansions of this locus have yielded different copy numbers among primates and rodents, with human retaining a single functional copy and mouse possessing several copies, depending on the strain. Roles for this gene in disease have been documented, but little is known about transcriptional regulation of NAIP. We show here that NAIP has multiple promoters sharing no similarity between human and rodents. Moreover, we demonstrate that multiple, domesticated long terminal repeats (LTRs) of endogenous retroviral elements provide NAIP promoter function in human, mouse, and rat. In human, an LTR serves as a tissue-specific promoter, active primarily in testis. However, in rodents, our evidence indicates that an ancestral LTR common to all rodent genes is the major, constitutive promoter for these genes, and that a second LTR found in two of the mouse genes is a minor promoter. Thus, independently acquired LTRs have assumed regulatory roles for orthologous genes, a remarkable evolutionary scenario. We also demonstrate that 5′ flanking regions of IAP family genes as a group, in both human and mouse are enriched for LTR insertions compared to average genes. We propose several potential explanations for these findings, including a hypothesis that recruitment of LTRs near NAIP or other IAP genes may represent a host-cell adaptation to modulate apoptotic responses.

When retroviruses infect cells, the viral DNA inserts into the cellular genome. If this happens in gametes (egg or sperm), the viral DNA will be transmitted from parent to offspring, like all chromosomal DNA. Through evolutionary time, such infections of gametes have been so prevalent that 8%–10% of the normal human and mouse genomes are now composed of ancient viral DNA, termed endogenous retroviruses (ERVs). In human, these ERVs are mutated or “dead” but it has been shown that ERV regulatory regions can be employed by the host to help control expression of cellular genes. Here, we report on a remarkable example of this phenomenon. We demonstrate that both the human and rodent neuronal apoptosis inhibitory protein (NAIP) genes, involved in preventing cell death, use different ERV sequences to drive gene expression. Moreover, in each of the primate and rodent lineages, two separate ERVs contribute to NAIP gene expression. This repeated ERV recruitment by NAIP genes throughout evolution is very unlikely to have occurred by chance. We offer a number of potential explanations, including the intriguing possibility that it may be advantageous for anti-cell death genes like NAIP to use ERVs to control their expression. These results support the view that not all retroviral remnants in our genome are simply junk DNA.

Prolactin in man: a tale of two promoters

The pituitary hormone prolactin (PRL) is best known for its role in the regulation of lactation. Recent evidence furthermore indicates PRL is required for normal reproduction in rodents. Here, we report on the insertion of two transposon-like DNA sequences in the human prolactin gene, which together function as an alternative promoter directing extrapituitary PRL expression. Indeed, the transposable elements contain transcription factor binding sites that have been shown to mediate PRL transcription in human uterine decidualised endometrial cells and lymphocytes. We hypothesize that the transposon insertion event has resulted in divergent (pituitary versus extrapituitary) expression of prolactin in primates, and in differential actions of pituitary versus extrapituitary prolactin in lactation versus pregnancy respectively. Importantly, the TE insertion might provide a context for some of the conflicting results obtained in studies of PRL function in mice and man.

Stochastic epigenetic silencing of retrotransposons: does stability come with age?

It has been hypothesized that phenotypic variation in mammals could in part be due to incomplete and variable silencing of retrotransposons in somatic cells. This theory is based on the fact that some recent endogenous retroviral (ERV) insertions in the mouse exert variable effects on genes in isogenic animals, depending on the variable state of ERV methylation. In this article, we review the evidence for this and related phenomena and suggest that such stochastic epigenetic silencing is restricted to very recent insertions. We also present a model to explain the acquisition of a more stable epigenetic state for transposable element insertions through time.

Retroviral elements and their hosts: insertional mutagenesis in the mouse germ line

The inbred mouse is an invaluable model for human biology and disease. Nevertheless, when considering genetic mechanisms of variation and disease, it is important to appreciate the significant differences in the spectra of spontaneous mutations that distinguish these species. While insertions of transposable elements are responsible for only approximately 0.1% of de novo mutations in humans, the figure is 100-fold higher in the laboratory mouse. This striking difference is largely due to the ongoing activity of mouse endogenous retroviral elements. Here we briefly review mouse endogenous retroviruses (ERVs) and their influence on gene expression, analyze mechanisms of interaction between ERVs and the host cell, and summarize the variety of mutations caused by ERV insertions. The prevalence of mouse ERV activity indicates that the genome of the laboratory mouse is presently behind in the "arms race" against invasion.

A Role for DNA Hypomethylation and Histone Acetylation in Maintaining Allele-Specific Expression of Mouse NKG2A in Developing and Mature NK Cells

The repertoire of receptors that is expressed by NK cells is critical for their ability to kill virally infected or transformed cells. However, the molecular mechanisms that determine whether and when NK receptor genes are transcribed during hemopoiesis remain unclear. In this study, we show that hypomethylation of a CpG-rich region in the mouse NKG2A gene is associated with transcription of NKG2A in ex vivo NK cells and NK cell lines. This observation was extended to various developmental stages of NK cells sorted from bone marrow, in which we demonstrate that the CpGs are methylated in the NKG2A-negative stages (hemopoietic stem cells, NK progenitors, and NKG2A-negative NK cells), and hypomethylated specifically in the NKG2A-positive NK cells. Furthermore, we provide evidence that DNA methylation is important in maintaining the allele-specific expression of NKG2A. Finally, we show that acetylated histones are associated with the CpG-rich region in NKG2A positive, but not negative, cell lines, and that treatment with the histone deacetylase inhibitor trichostatin A alone is sufficient to induce NKG2A expression. Treatment with the methyltransferase inhibitor 5-azacytidine only is insufficient to induce transcription, but cotreatment with both drugs resulted in a significantly greater induction, suggesting a cooperative role for DNA methylation and histone acetylation status in regulating gene expression. These results enhance our understanding of the formation and maintenance of NK receptor repertoires in developing and mature NK cells.