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

Visualized computational predictions of transcriptional effects by intronic endogenous retroviruses

When endogenous retroviruses (ERVs) or other transposable elements (TEs) insert into an intron, the consequence on gene transcription can range from negligible to a complete ablation of normal transcripts. With the advance of sequencing technology, more and more insertionally polymorphic or private TE insertions are being identified in humans and mice, of which some could have a significant impact on host gene expression. Nevertheless, an efficient and low cost approach to prioritize their potential effect on gene transcription has been lacking. By building a computational model based on artificial neural networks (ANN), we demonstrate the feasibility of using machine-learning approaches to predict the likelihood that intronic ERV insertions will have major effects on gene transcription, focusing on the two ERV families, namely Intracisternal A-type Particle (IAP) and Early Transposon (ETn)/MusD elements, which are responsible for the majority of ERV-induced mutations in mice. We trained the ANN model using properties associated with these ERVs known to cause germ-line mutations (positive cases) and properties associated with likely neutral ERVs of the same families (negative cases), and derived a set of prediction plots that can visualize the likelihood of affecting gene transcription by ERV insertions. Our results show a highly reliable prediction power of our model, and offer a potential approach to computationally screen for other types of TE insertions that may affect gene transcription or even cause disease.

Identification and association study with lung cancer for novel insertion polymorphisms of human endogenous retrovirus

Sequences of human endogenous retroviruses (HERVs) are members of the long terminal repeat (LTR) retrotransposon family. Although the expression of HERV has long been a topic of investigation, HERV-insertion polymorphisms are not well known, and a genetic association between HERV-insertion polymorphisms and cancer has never been reported. To identify novel HERV loci in the genome from cancer tissues, we carried out the inverse PCR method targeting a conserved LTR region of HML-2, which is the most recently acquired HERV group. Novel two insertions, HML-2_sLTR(1p13.2) and HML-2_sLTR(19q12), were identified as insertionally polymorphic solo LTRs. Furthermore, a significant prevalence of HML-2_sLTR(1p13.2) homozygosity was detected in female never-smoking patients aged 60 years and over who had lung adenocarcinoma [versus the other genotyping; odds ratio (OR): 1.97; 95% confidence interval (CI): 1.01-3.81]. In another cohort consisting of female never-smoking patients with lung adenocarcinoma, a prevalence of HML-2_sLTR(1p13.2) homozygosity tended to be high in patients aged 60 years and over (versus the other genotyping; OR: 2.03; 95% CI: 0.96-4.29), whereas a low prevalence of HML-2_sLTR(1p13.2) homozygosity was detected in patients <60 years old (versus the other genotyping; OR: 0.31; 95% CI: 0.11-0.94). Our results suggest that HML-2_sLTR(1p13.2) is involved with the susceptibility to lung adenocarcinoma in female never-smokers in an age-dependent manner and that other HERV polymorphisms related to human diseases might remain to be identified in the human genome.

Are human endogenous retroviruses pathogenic? An approach to testing the hypothesis

A number of observations have led researchers to postulate that, despite being replication-defective, human endogenous retroviruses (HERVs) may have retained the potential to cause or contribute to disease. However, mechanisms of HERV pathogenicity might differ substantially from those of modern infectious retroviruses or of the infectious precursors of HERVs. Therefore, novel pathways of HERV involvement in disease pathogenesis should be investigated. Recent technological advances in sequencing and bioinformatics are making this task increasingly feasible. The accumulating knowledge of HERV biology may also facilitate the definition and general acceptance of criteria that establish HERV pathogenicity. Here, we explore possible mechanisms whereby HERVs may cause disease and examine the evidence that either has been or should be obtained in order to decisively address the pathogenic potential of HERVs.

Endogenous gammaretrovirus acquisition in Mus musculus subspecies carrying functional variants of the XPR1 virus receptor

The xenotropic and polytropic mouse leukemia viruses (X-MLVs and P-MLVs, respectively) have different host ranges but use the same functionally polymorphic receptor, XPR1, for entry. Endogenous retroviruses (ERVs) of these 2 gammaretrovirus subtypes are largely segregated in different house mouse subspecies, but both MLV types are found in the classical strains of laboratory mice, which are genetic mosaics of 3 wild mouse subspecies. To describe the subspecies origins of laboratory mouse XP-MLV ERVs and their coevolutionary trajectory with their XPR1 receptor, we screened the house mouse subspecies for known and novel Xpr1 variants and for the individual full-length XP-MLV ERVs found in the sequenced C57BL mouse genome. The 12 X-MLV ERVs predate the origins of laboratory mice; they were all traced to Japanese wild mice and are embedded in the 5% of the laboratory mouse genome derived from the Asian Mus musculus musculus and, in one case, in the <1% derived from M. m. castaneus. While all 31 P-MLV ERVs map to the 95% of the laboratory mouse genome derived from P-MLV-infected M. m. domesticus, no C57BL P-MLV ERVs were found in wild M. m. domesticus. All M. m. domesticus mice carry the fully permissive XPR1 receptor allele, but all of the various restrictive XPR1 receptors, including the X-MLV-restricting laboratory mouse Xpr1(n) and a novel M. m. castaneus allele, originated in X-MLV-infected Asian mice. Thus, P-MLV ERVs show more insertional polymorphism than X-MLVs, and these differences in ERV acquisition and fixation are linked to subspecies-specific and functionally distinct XPR1 receptor variants.

Footprint of APOBEC3 on the genome of human retroelements

Almost half of the human genome is composed of transposable elements. The genomic structures and life cycles of some of these elements suggest they are a result of waves of retroviral infection and transposition over millions of years. The reduction of retrotransposition activity in primates compared to that in nonprimates, such as mice, has been attributed to the positive selection of several antiretroviral factors, such as apolipoprotein B mRNA editing enzymes. Among these, APOBEC3G is known to mutate G to A within the context of GG in the genome of endogenous as well as several exogenous retroelements (the underlining marks the G that is mutated). On the other hand, APOBEC3F and to a lesser extent other APOBEC3 members induce G-to-A changes within the nucleotide GA. It is known that these enzymes can induce deleterious mutations in the genome of retroviral sequences, but the evolution and/or inactivation of retroelements as a result of mutation by these proteins is not clear. Here, we analyze the mutation signatures of these proteins on large populations of long interspersed nuclear element (LINE), short interspersed nuclear element (SINE), and endogenous retrovirus (ERV) families in the human genome to infer possible evolutionary pressure and/or hypermutation events. Sequence context dependency of mutation by APOBEC3 allows investigation of the changes in the genome of retroelements by inspecting the depletion of G and enrichment of A within the APOBEC3 target and product motifs, respectively. Analysis of approximately 22,000 LINE-1 (L1), 24,000 SINE Alu, and 3,000 ERV sequences showed a footprint of GG→AG mutation by APOBEC3G and GA→AA mutation by other members of the APOBEC3 family (e.g., APOBEC3F) on the genome of ERV-K and ERV-1 elements but not on those of ERV-L, LINE, or SINE.

Distinct roles of KAP1, HP1 and G9a/GLP in silencing of the two-cell-specific retrotransposon MERVL in mouse ES cells

BACKGROUND

In mouse embryonic stem cells (mESCs), transcriptional silencing of numerous class I and II endogenous retroviruses (ERVs), including IAP, ETn and MMERVK10C, is dependent upon the H3K9 methyltransferase (KMTase) SETDB1/ESET and its binding partner KAP1/TRIM28. In contrast, the H3K9 KMTases G9a and GLP and HP1 proteins are dispensable for this process. Intriguingly, MERVL retroelements are actively transcribed exclusively in the two-cell (2C) embryo, but the molecular basis of silencing of these class III ERVs at later developmental stages has not been systematically addressed.

RESULTS

Here, we characterized the roles of these chromatin factors in MERVL silencing in mESCs. While MMERVK10C and IAP ERVs are bound by SETDB1 and KAP1 and are induced following their deletion, MERVL ERVs show relatively low levels of SETDB1 and KAP1 binding and are upregulated exclusively following KAP1 depletion, indicating that KAP1 influences MERVL expression independent of SETDB1. In contrast to class I and class II ERVs, MERVL and MERVL LTR-driven genic transcripts are also upregulated following depletion of G9a or GLP, and G9a binds directly to these ERVs. Consistent with a direct role for H3K9me2 in MERVL repression, these elements are highly enriched for G9a-dependent H3K9me2, and catalytically active G9a is required for silencing of MERVL LTR-driven transcripts. MERVL is also derepressed in HP1α and HP1β KO ESCs. However, like KAP1, HP1α and HP1β are only modestly enriched at MERVL relative to IAP LTRs. Intriguingly, as recently shown for KAP1, RYBP, LSD1 and G9a-deficient mESCs, many genes normally expressed in the 2C embryo are also induced in HP1 KO mESCs, revealing that aberrant expression of a subset of 2C-specific genes is a common feature in each of these KO lines.

CONCLUSIONS

Our results indicate that G9a and GLP, which are not required for silencing of class I and II ERVs, are recruited to MERVL elements and play a direct role in silencing of these class III ERVs, dependent upon G9a catalytic activity. In contrast, induction of MERVL expression in KAP1, HP1α and HP1β KO ESCs may occur predominantly as a consequence of indirect effects, in association with activation of a subset of 2C-specific genes.

The genome-defence gene Tex19.1 suppresses LINE-1 retrotransposons in the placenta and prevents intra-uterine growth retardation in mice

DNA methylation plays an important role in suppressing retrotransposon activity in mammalian genomes, yet there are stages of mammalian development where global hypomethylation puts the genome at risk of retrotransposition-mediated genetic instability. Hypomethylated primordial germ cells appear to limit this risk by expressing a cohort of retrotransposon-suppressing genome-defence genes whose silencing depends on promoter DNA methylation. Here, we investigate whether similar mechanisms operate in hypomethylated trophectoderm-derived components of the mammalian placenta to couple expression of genome-defence genes to the potential for retrotransposon activity. We show that the hypomethylated state of the mouse placenta results in activation of only one of the hypomethylation-sensitive germline genome-defence genes: Tex19.1. Tex19.1 appears to play an important role in placenta function as Tex19.1(-/-) mouse embryos exhibit intra-uterine growth retardation and have small placentas due to a reduction in the number of spongiotrophoblast, glycogen trophoblast and sinusoidal trophoblast giant cells. Furthermore, we show that retrotransposon mRNAs are derepressed in Tex19.1(-/-) placentas and that protein encoded by the LINE-1 retrotransposon is upregulated in hypomethylated trophectoderm-derived cells that normally express Tex19.1. This study suggests that post-transcriptional genome-defence mechanisms are operating in the placenta to protect the hypomethylated cells in this tissue from retrotransposons and suggests that imbalances between retrotransposon activity and genome-defence mechanisms could contribute to placenta dysfunction and disease.

De novo DNA methylation of endogenous retroviruses is shaped by KRAB-ZFPs/KAP1 and ESET

Endogenous retroviruses (ERVs) undergo de novo DNA methylation during the first few days of mammalian embryogenesis, although the factors that control the targeting of this process are largely unknown. We asked whether KAP1 (KRAB-associated protein 1) is involved in this mechanism because of its previously defined role in maintaining the silencing of ERVs through the histone methyltransferase ESET and histone H3 lysine 9 trimethylation. Here, we demonstrate that introduced ERV sequences are sufficient to direct rapid de novo methylation of a flanked promoter in embryonic stem (ES) cells. This mechanism requires the presence of an ERV sequence-recognizing KRAB zinc-finger protein (ZFP) and both KAP1 and ESET. Furthermore, this process can also take place on a strong cellular promoter and leads to methylation signatures that are subsequently maintained in vivo throughout embryogenesis. Finally, we show that methylation of ERVs residing in the genome is affected by knockout of KAP1 in early embryos. KRAB-ZFPs, KAP1 and ESET are thus likely to be responsible for the early embryonic instatement of stable epigenetic marks at ERV-containing loci.

Next-generation sequencing identifies the Danforth's short tail mouse mutation as a retrotransposon insertion affecting Ptf1a expression

The semidominant Danforth's short tail (Sd) mutation arose spontaneously in the 1920s. The homozygous Sd phenotype includes severe malformations of the axial skeleton with an absent tail, kidney agenesis, anal atresia, and persistent cloaca. The Sd mutant phenotype mirrors features seen in human caudal malformation syndromes including urorectal septum malformation, caudal regression, VACTERL association, and persistent cloaca. The Sd mutation was previously mapped to a 0.9 cM region on mouse chromosome 2qA3. We performed Sanger sequencing of exons and intron/exon boundaries mapping to the Sd critical region and did not identify any mutations. We then performed DNA enrichment/capture followed by next-generation sequencing (NGS) of the critical genomic region. Standard bioinformatic analysis of paired-end sequence data did not reveal any causative mutations. Interrogation of reads that had been discarded because only a single end mapped correctly to the Sd locus identified an early transposon (ETn) retroviral insertion at the Sd locus, located 12.5 kb upstream of the Ptf1a gene. We show that Ptf1a expression is significantly upregulated in Sd mutant embryos at E9.5. The identification of the Sd mutation will lead to improved understanding of the developmental pathways that are misregulated in human caudal malformation syndromes.

Birth defects are the leading cause of infant mortality in the United States, accounting for 1 in 5 infant deaths annually. Birth defects that affect development of the caudal portion of the embryo can include malformations of the spine, such as spina bifida, and malformations of the kidneys and lower gastrointestinal tract. Little is known regarding the genetic causes of human caudal birth defects. The Danforth's short tail (Sd) mouse shares many similarities with these caudal birth defects that occur in humans. In this manuscript, we used next-generation sequencing to identify the genetic cause of the Sd mouse phenotype. We found that the Sd mutation is a retrotransposon insertion that inappropriately turns on a nearby gene that is normally important for pancreas development. Future studies of Sd mice will help us understand the pathogenesis of caudal birth defects in humans.

Ectopic expression of Ptf1a induces spinal defects, urogenital defects, and anorectal malformations in Danforth's short tail mice

Danforth's short tail (Sd) is a semidominant mutation on mouse chromosome 2, characterized by spinal defects, urogenital defects, and anorectal malformations. However, the gene responsible for the Sd phenotype was unknown. In this study, we identified the molecular basis of the Sd mutation. By positional cloning, we identified the insertion of an early transposon in the Sd candidate locus approximately 12-kb upstream of Ptf1a. We found that insertion of the transposon caused overexpression of three neighboring genes, Gm13344, Gm13336, and Ptf1a, in Sd mutant embryos and that the Sd phenotype was not caused by disruption of an as-yet-unknown gene in the candidate locus. Using multiple knockout and knock-in mouse models, we demonstrated that misexpression of Ptf1a, but not of Gm13344 or Gm13336, in the notochord, hindgut, cloaca, and mesonephros was sufficient to replicate the Sd phenotype. The ectopic expression of Ptf1a in the caudal embryo resulted in attenuated expression of Cdx2 and its downstream target genes T, Wnt3a, and Cyp26a1; we conclude that this is the molecular basis of the Sd phenotype. Analysis of Sd mutant mice will provide insight into the development of the spinal column, anus, and kidney.

Caudal regression syndrome (CRS) is a congenital heterogeneous constellation of caudal anomalies that includes varying degrees of agenesis of the spinal column, anorectal malformations, and genitourinary anomalies. Its pathogenesis is unclear. However, it could be the result of excessive physiologic regression of the embryonic caudal region based on analyses of the various mouse mutants carrying caudal agenesis. Among the mouse mutants, the Danforth's short tail (Sd) mouse is considered a best model for human CRS. Sd is a semidominant mutation, characterized by spinal defects, urogenital defects, and anorectal malformations, thus showing phenotypic similarity to human CRS. Although Sd is known to map to mouse chromosome 2, little is known about the molecular nature of the mutation. Here, we demonstrate an insertion of one type of retrotransposon near the Ptf1a gene. This resulted in ectopic expression of Ptf1a gene in the caudal region of the embryo and downregulation of Cdx2 and its downstream targets, leading to characteristic phenotypes in Sd mouse. Thus, Sd mutant mice will provide insight into the development of the spinal column, anus, and kidney.

A retrotransposon insertion in the 5' regulatory domain of Ptf1a results in ectopic gene expression and multiple congenital defects in Danforth's short tail mouse

Danforth's short tail mutant (Sd) mouse, first described in 1930, is a classic spontaneous mutant exhibiting defects of the axial skeleton, hindgut, and urogenital system. We used meiotic mapping in 1,497 segregants to localize the mutation to a 42.8-kb intergenic segment on chromosome 2. Resequencing of this region identified an 8.5-kb early retrotransposon (ETn) insertion within the highly conserved regulatory sequences upstream of Pancreas Specific Transcription Factor, 1a (Ptf1a). This mutation resulted in up to tenfold increased expression of Ptf1a as compared to wild-type embryos at E9.5 but no detectable changes in the expression levels of other neighboring genes. At E9.5, Sd mutants exhibit ectopic Ptf1a expression in embryonic progenitors of every organ that will manifest a developmental defect: the notochord, the hindgut, and the mesonephric ducts. Moreover, at E 8.5, Sd mutant mice exhibit ectopic Ptf1a expression in the lateral plate mesoderm, tail bud mesenchyme, and in the notochord, preceding the onset of visible defects such as notochord degeneration. The Sd heterozygote phenotype was not ameliorated by Ptf1a haploinsufficiency, further suggesting that the developmental defects result from ectopic expression of Ptf1a. These data identify disruption of the spatio-temporal pattern of Ptf1a expression as the unifying mechanism underlying the multiple congenital defects in Danforth's short tail mouse. This striking example of an enhancer mutation resulting in profound developmental defects suggests that disruption of conserved regulatory elements may also contribute to human malformation syndromes.

Birth defects are a major cause of childhood morbidity and mortality. We studied the Danforth's short tail mouse, a classic mouse model of birth defects involving the skeleton, gut, and urinary system. We precisely localized the mutation responsible for these birth defects to a 42.8-kb segment on chromosome 2 and identified the mutation as an 8.5-kb transposon that disrupts highly conserved regulatory sequences upstream of the Pancreas Specific Transcription Factor, 1a (Ptf1a). The insertion disrupts a Ptf1a regulatory domain that is highly conserved across evolution and results in spatiotemporal defects in Ptf1a expression: we detected increased expression, temporally premature expression, and (most important for elucidating the mutant phenotype) the ectopic expression of Ptf1a in the notochord, hindgut, and mesonephros—the three sites that will give rise to organ defects in Danforth's short tail mouse. Our data also provide a striking example of how a noncoding, regulatory mutation can produce transient spatio-temporal dsyregulation of gene expression and result in profound developmental defects, highlighting the critical role of noncoding elements for coordinated gene expression in the vertebrate genome. Finally, these data provide novel insight into the role of Ptf1a in embryogenesis and lay the groundwork for elucidation of novel mechanisms underlying birth defects in humans.

Latent regulatory potential of human-specific repetitive elements

At least half of the human genome is derived from repetitive elements, which are often lineage specific and silenced by a variety of genetic and epigenetic mechanisms. Using a transchromosomic mouse strain that transmits an almost complete single copy of human chromosome 21 via the female germline, we show that a heterologous regulatory environment can transcriptionally activate transposon-derived human regulatory regions. In the mouse nucleus, hundreds of locations on human chromosome 21 newly associate with activating histone modifications in both somatic and germline tissues, and influence the gene expression of nearby transcripts. These regions are enriched with primate and human lineage-specific transposable elements, and their activation corresponds to changes in DNA methylation at CpG dinucleotides. This study reveals the latent regulatory potential of the repetitive human genome and illustrates the species specificity of mechanisms that control it.

Phylogenetic and DNA methylation analysis reveal novel regions of variable methylation in the mouse IAP class of transposons

Background

Select retrotransposons in the long terminal repeat (LTR) class exhibit interindividual variation in DNA methylation that is altered by developmental environmental exposures. Yet, neither the full extent of variability at these “metastable epialleles,” nor the phylogenetic relationship underlying variable elements is well understood. The murine metastable epialleles, A^vy and Cabp^IAP, result from independent insertions of an intracisternal A particle (IAP) mobile element, and exhibit remarkably similar sequence identity (98.5%).

Results

Utilizing the C57BL/6 genome we identified 10802 IAP LTRs overall and a subset of 1388 in a family that includes A^vy and Cabp^IAP. Phylogenetic analysis revealed two duplication and divergence events subdividing this family into three clades. To characterize interindividual variation across clades, liver DNA from 17 isogenic mice was subjected to combined bisulfite and restriction analysis (CoBRA) for 21 separate LTR transposons (7 per clade). The lowest and highest mean methylation values were 59% and 88% respectively, while methylation levels at individual LTRs varied widely, ranging from 9% to 34%. The clade with the most conserved elements had significantly higher mean methylation across LTRs than either of the two diverged clades (p = 0.040 and p = 0.017). Within each mouse, average methylation across all LTRs was not significantly different (71%-74%, p > 0.99).

Conclusions

Combined phylogenetic and DNA methylation analysis allows for the identification of novel regions of variable methylation. This approach increases the number of known metastable epialleles in the mouse, which can serve as biomarkers for environmental modifications to the epigenome.

Retrotransposon insertion in the T-cell acute lymphocytic leukemia 1 (Tal1) gene is associated with severe renal disease and patchy alopecia in Hairpatches (Hpt) mice

“Hairpatches” (Hpt) is a naturally occurring, autosomal semi-dominant mouse mutation. Hpt/Hpt homozygotes die in utero, while Hpt/+ heterozygotes exhibit progressive renal failure accompanied by patchy alopecia. This mutation is a model for the rare human disorder “glomerulonephritis with sparse hair and telangiectases" (OMIM 137940). Fine mapping localized the Hpt locus to a 6.7 Mb region of Chromosome 4 containing 62 known genes. Quantitative real time PCR revealed differential expression for only one gene in the interval, T-cell acute lymphocytic leukemia 1 (Tal1), which was highly upregulated in the kidney and skin of Hpt/+ mice. Southern blot analysis of Hpt mutant DNA indicated a new EcoRI site in the Tal1 gene. High throughput sequencing identified an endogenous retroviral class II intracisternal A particle insertion in Tal1 intron 4. Our data suggests that the IAP insertion in Tal1 underlies the histopathological changes in the kidney by three weeks of age, and that glomerulosclerosis is a consequence of an initial developmental defect, progressing in severity over time. The Hairpatches mouse model allows an investigation into the effects of Tal1, a transcription factor characterized by complex regulation patterns, and its effects on renal disease.

The APOBEC3 family of retroelement restriction factors

The ability to regulate and even target mutagenesis is an extremely valuable cellular asset. Enzyme-catalyzed DNA cytosine deamination is a molecular strategy employed by vertebrates to promote antibody diversity and defend against foreign nucleic acids. Ten years ago, a family of cellular enzymes was first described with several proving capable of deaminating DNA and inhibiting HIV-1 replication. Ensuing studies on the apolipoprotein B mRNA-editing enzyme catalytic polypeptide-like 3 (APOBEC3) restriction factors have uncovered a broad-spectrum innate defense network that suppresses the replication of numerous endogenous and exogenous DNA-based parasites. Although many viruses possess equally elaborate counter-defense mechanisms, the APOBEC3 enzymes offer a tantalizing possibility of leveraging innate immunity to fend off viral infection. Here, we focus on mechanisms of retroelement restriction by the APOBEC3 family of restriction enzymes, and we consider the therapeutic benefits, as well as the possible pathological consequences, of arming cells with active DNA deaminases.

RetroSeq: transposable element discovery from next-generation sequencing data

A significant proportion of eukaryote genomes consist of transposable element (TE)-derived sequence. These elements are known to have the capacity to modulate gene function and genome evolution. We have developed RetroSeq for detecting non-reference TE insertions from Illumina paired-end whole-genome sequencing data. We evaluate RetroSeq on a human trio from the 1000 Genomes Project, showing that it produces highly accurate TE calls.

AVAILABILTY

RetroSeq is open-source and available from https://github.com/tk2/RetroSeq.

The correlation of genome size and DNA methylation rate in metazoans

Total DNA methylation rates are well known to vary widely between different metazoans. The phylogenetic distribution of this variation, however, has not been investigated systematically. We combine here publicly available data on methylcytosine content with the analysis of nucleotide compositions of genomes and transcriptomes of 78 metazoan species to trace the evolution of abundance and distribution of DNA methylation. The depletion of CpG and the associated enrichment of TpG and CpA dinucleotides are used to infer the intensity and localization of germline CpG methylation and to estimate its evolutionary dynamics. We observe a positive correlation of the relative methylation of CpG motifs with genome size. We tested this trend successfully by measuring total DNA methylation with LC/MS in orthopteran insects with very different genome sizes: house crickets, migratory locusts and meadow grasshoppers. We hypothesize that the observed correlation between methylation rate and genome size is due to a dependence of both variables from long-term effective population size and is driven by the accumulation of repetitive sequences that are typically methylated during periods of small population sizes. This process may result in generally methylated, large genomes such as those of jawed vertebrates. In this case, the emergence of a novel demethylation pathway and of novel reader proteins for methylcytosine may have enabled the usage of cytosine methylation for promoter-based gene regulation. On the other hand, persistently large populations may lead to a compression of the genome and to the loss of the DNA methylation machinery, as observed, e.g., in nematodes.

Transposable elements and viruses as factors in adaptation and evolution: an expansion and strengthening of the TE-Thrust hypothesis

In addition to the strong divergent evolution and significant and episodic evolutionary transitions and speciation we previously attributed to TE-Thrust, we have expanded the hypothesis to more fully account for the contribution of viruses to TE-Thrust and evolution. The concept of symbiosis and holobiontic genomes is acknowledged, with particular emphasis placed on the creativity potential of the union of retroviral genomes with vertebrate genomes. Further expansions of the TE-Thrust hypothesis are proposed regarding a fuller account of horizontal transfer of TEs, the life cycle of TEs, and also, in the case of a mammalian innovation, the contributions of retroviruses to the functions of the placenta. The possibility of drift by TE families within isolated demes or disjunct populations, is acknowledged, and in addition, we suggest the possibility of horizontal transposon transfer into such subpopulations. “Adaptive potential” and “evolutionary potential” are proposed as the extremes of a continuum of “intra-genomic potential” due to TE-Thrust. Specific data is given, indicating “adaptive potential” being realized with regard to insecticide resistance, and other insect adaptations. In this regard, there is agreement between TE-Thrust and the concept of adaptation by a change in allele frequencies. Evidence on the realization of “evolutionary potential” is also presented, which is compatible with the known differential survivals, and radiations of lineages. Collectively, these data further suggest the possibility, or likelihood, of punctuated episodes of speciation events and evolutionary transitions, coinciding with, and heavily underpinned by, intermittent bursts of TE activity.

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.

Origins of albino and hooded rats: implications from molecular genetic analysis across modern laboratory rat strains

Albino and hooded (or piebald) rats are one of the most frequently used laboratory animals for the past 150 years. Despite this fact, the origin of the albino mutation as well as the genetic basis of the hooded phenotype remained unclear. Recently, the albino mutation has been identified as the Arg299His missense mutation in the Tyrosinase gene and the hooded (H) locus has been mapped to the ∼460-kb region in which only the Kit gene exists. Here, we surveyed 172 laboratory rat strains for the albino mutation and the hooded (h) mutation that we identified by positional cloning approach to investigate possible genetic roots and relationships of albino and hooded rats. All of 117 existing laboratory albino rats shared the same albino missense mutation, indicating they had only one single ancestor. Genetic fine mapping followed by de novo sequencing of BAC inserts covering the H locus revealed that an endogenous retrovirus (ERV) element was inserted into the first intron of the Kit gene where the hooded allele maps. A solitary long terminal repeat (LTR) was found at the same position to the ERV insertion in another allele of the H locus, which causes the so called Irish (hⁱ) phenotype. The ERV and the solitary LTR insertions were completely associated with the hooded and Irish coat patterns, respectively, across all colored rat strains examined. Interestingly, all 117 albino rat strains shared the ERV insertion without any exception, which strongly suggests that the albino mutation had originally occurred in hooded rats.

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.

Expression of endogenous retroviruses is negatively regulated by the pluripotency marker Rex1/Zfp42

Rex1/Zfp42 is a Yy1-related zinc-finger protein whose expression is frequently used to identify pluripotent stem cells. We show that depletion of Rex1 levels notably affected self-renewal of mouse embryonic stem (ES) cells in clonal assays, in the absence of evident differences in expression of marker genes for pluripotency or differentiation. By contrast, marked differences in expression of several endogenous retroviral elements (ERVs) were evident upon Rex1 depletion. We demonstrate association of REX1 to specific elements in chromatin-immunoprecipitation assays, most strongly to muERV-L and to a lower extent to IAP and musD elements. Rex1 regulates muERV-L expression in vivo, as we show altered levels upon transient gain-and-loss of Rex1 function in pre-implantation embryos. We also find REX1 can associate with the lysine-demethylase LSD1/KDM1A, suggesting they act in concert. Similar to REX1 binding to retrotransposable elements (REs) in ES cells, we also detected binding of the REX1 related proteins YY1 and YY2 to REs, although the binding preferences of the two proteins were slightly different. Altogether, we show that Rex1 regulates ERV expression in mouse ES cells and during pre-implantation development and suggest that Rex1 and its relatives have evolved as regulators of endogenous retroviral transcription.

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.

Alopecia in a viable phospholipase C delta 1 and phospholipase C delta 3 double mutant

BACKGROUND

Inositol 1,4,5trisphosphate (IP(3)) and diacylglycerol (DAG) are important intracellular signalling molecules in various tissues. They are generated by the phospholipase C family of enzymes, of which phospholipase C delta (PLCD) forms one class. Studies with functional inactivation of Plcd isozyme encoding genes in mice have revealed that loss of both Plcd1 and Plcd3 causes early embryonic death. Inactivation of Plcd1 alone causes loss of hair (alopecia), whereas inactivation of Plcd3 alone has no apparent phenotypic effect. To investigate a possible synergy of Plcd1 and Plcd3 in postnatal mice, novel mutations of these genes compatible with life after birth need to be found.

METHODOLOGY/PRINCIPAL FINDINGS

We characterise a novel mouse mutant with a spontaneously arisen mutation in Plcd3 (Plcd3(mNab)) that resulted from the insertion of an intracisternal A particle (IAP) into intron 2 of the Plcd3 gene. This mutation leads to the predominant expression of a truncated PLCD3 protein lacking the N-terminal PH domain. C3H mice that carry one or two mutant Plcd3(mNab) alleles are phenotypically normal. However, the presence of one Plcd3(mNab) allele exacerbates the alopecia caused by the loss of functional Plcd1 in Del(9)olt1Pas mutant mice with respect to the number of hair follicles affected and the body region involved. Mice double homozygous for both the Del(9)olt1Pas and the Plcd3(mNab) mutations survive for several weeks and exhibit total alopecia associated with fragile hair shafts showing altered expression of some structural genes and shortened phases of proliferation in hair follicle matrix cells.

CONCLUSIONS/SIGNIFICANCE

The Plcd3(mNab) mutation is a novel hypomorphic mutation of Plcd3. Our investigations suggest that Plcd1 and Plcd3 have synergistic effects on the murine hair follicle in specific regions of the body surface.

The genomic landscape shaped by selection on transposable elements across 18 mouse strains

Background

Transposable element (TE)-derived sequence dominates the landscape of mammalian genomes and can modulate gene function by dysregulating transcription and translation. Our current knowledge of TEs in laboratory mouse strains is limited primarily to those present in the C57BL/6J reference genome, with most mouse TEs being drawn from three distinct classes, namely short interspersed nuclear elements (SINEs), long interspersed nuclear elements (LINEs) and the endogenous retrovirus (ERV) superfamily. Despite their high prevalence, the different genomic and gene properties controlling whether TEs are preferentially purged from, or are retained by, genetic drift or positive selection in mammalian genomes remain poorly defined.

Results

Using whole genome sequencing data from 13 classical laboratory and 4 wild-derived mouse inbred strains, we developed a comprehensive catalogue of 103,798 polymorphic TE variants. We employ this extensive data set to characterize TE variants across the Mus lineage, and to infer neutral and selective processes that have acted over 2 million years. Our results indicate that the majority of TE variants are introduced though the male germline and that only a minority of TE variants exert detectable changes in gene expression. However, among genes with differential expression across the strains there are twice as many TE variants identified as being putative causal variants as expected.

Conclusions

Most TE variants that cause gene expression changes appear to be purged rapidly by purifying selection. Our findings demonstrate that past TE insertions have often been highly deleterious, and help to prioritize TE variants according to their likely contribution to gene expression or phenotype variation.

Studies of endogenous retroviruses reveal a continuing evolutionary saga

Retroviral replication involves the formation of a DNA provirus integrated into the host genome. Through this process, retroviruses can colonize the germ line to form endogenous retroviruses (ERVs). ERV inheritance can have multiple adverse consequences for the host, some resembling those resulting from exogenous retrovirus infection but others arising by mechanisms unique to ERVs. Inherited retroviruses can also confer benefits on the host. To meet the different threats posed by endogenous and exogenous retroviruses, various host defences have arisen during evolution, acting at various stages on the retrovirus life cycle. In this Review, I describe our current understanding of the distribution and architecture of ERVs, the consequences of their acquisition for the host and the emerging details of the intimate evolutionary relationship between virus and vertebrate host.

Microarray analysis of LTR retrotransposon silencing identifies Hdac1 as a regulator of retrotransposon expression in mouse embryonic stem cells

Retrotransposons are highly prevalent in mammalian genomes due to their ability to amplify in pluripotent cells or developing germ cells. Host mechanisms that silence retrotransposons in germ cells and pluripotent cells are important for limiting the accumulation of the repetitive elements in the genome during evolution. However, although silencing of selected individual retrotransposons can be relatively well-studied, many mammalian retrotransposons are seldom analysed and their silencing in germ cells, pluripotent cells or somatic cells remains poorly understood. Here we show, and experimentally verify, that cryptic repetitive element probes present in Illumina and Affymetrix gene expression microarray platforms can accurately and sensitively monitor repetitive element expression data. This computational approach to genome-wide retrotransposon expression has allowed us to identify the histone deacetylase Hdac1 as a component of the retrotransposon silencing machinery in mouse embryonic stem cells, and to determine the retrotransposon targets of Hdac1 in these cells. We also identify retrotransposons that are targets of other retrotransposon silencing mechanisms such as DNA methylation, Eset-mediated histone modification, and Ring1B/Eed-containing polycomb repressive complexes in mouse embryonic stem cells. Furthermore, our computational analysis of retrotransposon silencing suggests that multiple silencing mechanisms are independently targeted to retrotransposons in embryonic stem cells, that different genomic copies of the same retrotransposon can be differentially sensitive to these silencing mechanisms, and helps define retrotransposon sequence elements that are targeted by silencing machineries. Thus repeat annotation of gene expression microarray data suggests that a complex interplay between silencing mechanisms represses retrotransposon loci in germ cells and embryonic stem cells.

Env-less endogenous retroviruses are genomic superspreaders

Endogenous retroviruses (ERVs) differ from typical retroviruses in being inherited through the host germline and therefore are a unique combination of pathogen and selfish genetic element. Some ERV lineages proliferate by infecting germline cells, as do typical retroviruses, whereas others lack the env gene required for virions to enter cells and thus behave like retrotransposons. We wished to know what factors determined the relative abundance of different ERV lineages, so we analyzed ERV loci recovered from 38 mammal genomes by in silico screening. By modeling the relationship between proliferation and replication mechanism in detail within one group, the intracisternal A-type particles (IAPs), and performing simple correlations across all ERV lineages, we show that when ERVs lose the env gene their proliferation within that genome is boosted by a factor of ∼30. We also show that ERV abundance follows the Pareto principle or 20/80 rule, with ∼20% of lineages containing 80% of the loci. This rule is observed in many biological systems, including infectious disease epidemics, where commonly ∼20% of the infected individuals are responsible for 80% of onward infection. We thus borrow simple epidemiological and ecological models and show that retrotransposition and loss of env is the trait that leads endogenous retroviruses to becoming genomic superspreaders that take over a significant proportion of their host's genome.

Endogenous viruses: insights into viral evolution and impact on host biology

Recent studies have uncovered myriad viral sequences that are integrated or 'endogenized' in the genomes of various eukaryotes. Surprisingly, it appears that not just retroviruses but almost all types of viruses can become endogenous. We review how these genomic 'fossils' offer fresh insights into the origin, evolutionary dynamics and structural evolution of viruses, which are giving rise to the burgeoning field of palaeovirology. We also examine the multitude of ways through which endogenous viruses have influenced, for better or worse, the biology of their hosts. We argue that the conflict between hosts and viruses has led to the invention and diversification of molecular arsenals, which, in turn, promote the cellular co-option of endogenous viruses.

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.

Bioinformatics and genomic analysis of transposable elements in eukaryotic genomes

A major portion of most eukaryotic genomes are transposable elements (TEs). During evolution, TEs have introduced profound changes to genome size, structure, and function. As integral parts of genomes, the dynamic presence of TEs will continue to be a major force in reshaping genomes. Early computational analyses of TEs in genome sequences focused on filtering out "junk" sequences to facilitate gene annotation. When the high abundance and diversity of TEs in eukaryotic genomes were recognized, these early efforts transformed into the systematic genome-wide categorization and classification of TEs. The availability of genomic sequence data reversed the classical genetic approaches to discovering new TE families and superfamilies. Curated TE databases and their accurate annotation of genome sequences in turn facilitated the studies on TEs in a number of frontiers including: (1) TE-mediated changes of genome size and structure, (2) the influence of TEs on genome and gene functions, (3) TE regulation by host, (4) the evolution of TEs and their population dynamics, and (5) genomic scale studies of TE activity. Bioinformatics and genomic approaches have become an integral part of large-scale studies on TEs to extract information with pure in silico analyses or to assist wet lab experimental studies. The current revolution in genome sequencing technology facilitates further progress in the existing frontiers of research and emergence of new initiatives. The rapid generation of large-sequence datasets at record low costs on a routine basis is challenging the computing industry on storage capacity and manipulation speed and the bioinformatics community for improvement in algorithms and their implementations.

Mutation of Rubie, a novel long non-coding RNA located upstream of Bmp4, causes vestibular malformation in mice

Background

The vestibular apparatus of the vertebrate inner ear uses three fluid-filled semicircular canals to sense angular acceleration of the head. Malformation of these canals disrupts the sense of balance and frequently causes circling behavior in mice. The Epistatic circler (Ecl) is a complex mutant derived from wildtype SWR/J and C57L/J mice. Ecl circling has been shown to result from the epistatic interaction of an SWR-derived locus on chromosome 14 and a C57L-derived locus on chromosome 4, but the causative genes have not been previously identified.

Methodology/Principal Findings

We developed a mouse chromosome substitution strain (CSS-14) that carries an SWR/J chromosome 14 on a C57BL/10J genetic background and, like Ecl, exhibits circling behavior due to lateral semicircular canal malformation. We utilized CSS-14 to identify the chromosome 14 Ecl gene by positional cloning. Our candidate interval is located upstream of bone morphogenetic protein 4 (Bmp4) and contains an inner ear-specific, long non-coding RNA that we have designated Rubie (RNA upstream of Bmp4 expressed in inner ear). Rubie is spliced and polyadenylated, and is expressed in developing semicircular canals. However, we discovered that the SWR/J allele of Rubie is disrupted by an intronic endogenous retrovirus that causes aberrant splicing and premature polyadenylation of the transcript. Rubie lies in the conserved gene desert upstream of Bmp4, within a region previously shown to be important for inner ear expression of Bmp4. We found that the expression patterns of Bmp4 and Rubie are nearly identical in developing inner ears.

Conclusions/Significance

Based on these results and previous studies showing that Bmp4 is essential for proper vestibular development, we propose that Rubie is the gene mutated in Ecl mice, that it is involved in regulating inner ear expression of Bmp4, and that aberrant Bmp4 expression contributes to the Ecl phenotype.

Expressing genes do not forget their LINEs: transposable elements and gene expression

Historically the accumulated mass of mammalian transposable elements (TEs), particularly those located within gene boundaries, was viewed as a genetic burden potentially detrimental to the genomic landscape. This notion has been strengthened by the discovery that transposable sequences can alter the architecture of the transcriptome, not only through insertion, but also long after the integration process is completed. Insertions previously considered harmless are now known to impact the expression of host genes via modification of the transcript quality or quantity, transcriptional interference, or by the control of pathways that affect the mRNA life-cycle. Conversely, several examples of the evolutionary advantageous impact of TEs on the host gene structure that diversified the cellular transcriptome are reported. TE-induced changes in gene expression can be tissue- or disease-specific, raising the possibility that the impact of TE sequences may vary during development, among normal cell types, and between normal and disease-affected tissues. The understanding of the rules and abundance of TE-interference with gene expression is in its infancy, and its contribution to human disease and/or evolution remains largely unexplored.

Active transposition in genomes

Transposons are DNA sequences capable of moving in genomes. Early evidence showed their accumulation in many species and suggested their continued activity in at least isolated organisms. In the past decade, with the development of various genomic technologies, it has become abundantly clear that ongoing activity is the rule rather than the exception. Active transposons of various classes are observed throughout plants and animals, including humans. They continue to create new insertions, have an enormous variety of structural and functional impact on genes and genomes, and play important roles in genome evolution. Transposon activities have been identified and measured by employing various strategies. Here, we summarize evidence of current transposon activity in various plant and animal genomes.

Polymorphic integrations of an endogenous gammaretrovirus in the mule deer genome

Endogenous retroviruses constitute a significant genomic fraction in all mammalian species. Typically they are evolutionarily old and fixed in the host species population. Here we report on a novel endogenous gammaretrovirus (CrERVγ; for cervid endogenous gammaretrovirus) in the mule deer (Odocoileus hemionus) that is insertionally polymorphic among individuals from the same geographical location, suggesting that it has a more recent evolutionary origin. Using PCR-based methods, we identified seven CrERVγ proviruses and demonstrated that they show various levels of insertional polymorphism in mule deer individuals. One CrERVγ provirus was detected in all mule deer sampled but was absent from white-tailed deer, indicating that this virus originally integrated after the split of the two species, which occurred approximately one million years ago. There are, on average, 100 CrERVγ copies in the mule deer genome based on quantitative PCR analysis. A CrERVγ provirus was sequenced and contained intact open reading frames (ORFs) for three virus genes. Transcripts were identified covering the entire provirus. CrERVγ forms a distinct branch of the gammaretrovirus phylogeny, with the closest relatives of CrERVγ being endogenous gammaretroviruses from sheep and pig. We demonstrated that white-tailed deer (Odocoileus virginianus) and elk (Cervus canadensis) DNA contain proviruses that are closely related to mule deer CrERVγ in a conserved region of pol; more distantly related sequences can be identified in the genome of another member of the Cervidae, the muntjac (Muntiacus muntjak). The discovery of a novel transcriptionally active and insertionally polymorphic retrovirus in mammals could provide a useful model system to study the dynamic interaction between the host genome and an invading retrovirus.

Global methylation in exposure biology and translational medical science

BACKGROUND

Many groups are actively investigating how the epigenetic state relates to environmental exposures and development of disease, including cancer. There are myriad choices for capturing and measuring the epigenetic state of a tissue, ranging from assessing the total methyl-CpG content to array-based platforms that simultaneously probe hundreds of thousands of CpG loci. There is an emerging literature that uses CpG methylation at repetitive sequences, including LINE-1 (long interspersed nuclear element-1) elements, to capture the epigenomic state.

OBJECTIVES

We explored the complexity of using CpG methylation at repetitive sequences in epidemiology and translational medical research and suggest needed avenues of research to clarify its meaning and utility.

CONCLUSIONS

Among the most urgent avenues of research is the need for prospective studies to eliminate the possibilities of reverse causality, and development of new LINE-1 assays that capture both class of LINE-1 element and copy number.

LINE-1 elements in structural variation and disease

The completion of the human genome reference sequence ushered in a new era for the study and discovery of human transposable elements. It now is undeniable that transposable elements, historically dismissed as junk DNA, have had an instrumental role in sculpting the structure and function of our genomes. In particular, long interspersed element-1 (LINE-1 or L1) and short interspersed elements (SINEs) continue to affect our genome, and their movement can lead to sporadic cases of disease. Here, we briefly review the types of transposable elements present in the human genome and their mechanisms of mobility. We next highlight how advances in DNA sequencing and genomic technologies have enabled the discovery of novel retrotransposons in individual genomes. Finally, we discuss how L1-mediated retrotransposition events impact human genomes.

Characterization and potential functional significance of human-chimpanzee large INDEL variation

BACKGROUND

Although humans and chimpanzees have accumulated significant differences in a number of phenotypic traits since diverging from a common ancestor about six million years ago, their genomes are more than 98.5% identical at protein-coding loci. This modest degree of nucleotide divergence is not sufficient to explain the extensive phenotypic differences between the two species. It has been hypothesized that the genetic basis of the phenotypic differences lies at the level of gene regulation and is associated with the extensive insertion and deletion (INDEL) variation between the two species. To test the hypothesis that large INDELs (80 to 12,000 bp) may have contributed significantly to differences in gene regulation between the two species, we categorized human-chimpanzee INDEL variation mapping in or around genes and determined whether this variation is significantly correlated with previously determined differences in gene expression.

RESULTS

Extensive, large INDEL variation exists between the human and chimpanzee genomes. This variation is primarily attributable to retrotransposon insertions within the human lineage. There is a significant correlation between differences in gene expression and large human-chimpanzee INDEL variation mapping in genes or in proximity to them.

CONCLUSIONS

The results presented herein are consistent with the hypothesis that large INDELs, particularly those associated with retrotransposons, have played a significant role in human-chimpanzee regulatory evolution.

Retrotransposon-induced heterochromatin spreading in the mouse revealed by insertional polymorphisms

The "arms race" relationship between transposable elements (TEs) and their host has promoted a series of epigenetic silencing mechanisms directed against TEs. Retrotransposons, a class of TEs, are often located in repressed regions and are thought to induce heterochromatin formation and spreading. However, direct evidence for TE-induced local heterochromatin in mammals is surprisingly scarce. To examine this phenomenon, we chose two mouse embryonic stem (ES) cell lines that possess insertionally polymorphic retrotransposons (IAP, ETn/MusD, and LINE elements) at specific loci in one cell line but not the other. Employing ChIP-seq data for these cell lines, we show that IAP elements robustly induce H3K9me3 and H4K20me3 marks in flanking genomic DNA. In contrast, such heterochromatin is not induced by LINE copies and only by a minority of polymorphic ETn/MusD copies. DNA methylation is independent of the presence of IAP copies, since it is present in flanking regions of both full and empty sites. Finally, such spreading into genes appears to be rare, since the transcriptional start sites of very few genes are less than one Kb from an IAP. However, the B3galtl gene is subject to transcriptional silencing via IAP-induced heterochromatin. Hence, although rare, IAP-induced local heterochromatin spreading into nearby genes may influence expression and, in turn, host fitness.

Retrotransposon-centered analysis of piRNA targeting shows a shift from active to passive retrotransposon transcription in developing mouse testes

Background

Piwi-associated RNAs (piRNAs) bind transcripts from retrotransposable elements (RTE) in mouse germline cells and seemingly act as guides for genomic methylation, thereby repressing the activity of RTEs. It is currently unknown if and how Piwi proteins distinguish RTE transcripts from other cellular RNAs. During germline development, the main target of piRNAs switch between different types of RTEs. Using the piRNA targeting of RTEs as an indicator of RTE activity, and considering the entire population of genomic RTE loci along with their age and location, this study aims at further elucidating the dynamics of RTE activity during mouse germline development.

Results

Due to the inherent sequence redundancy between RTE loci, assigning piRNA targeting to specific loci is problematic. This limits the analysis, although certain features of piRNA targeting of RTE loci are apparent. As expected, young RTEs display a much higher level of piRNA targeting than old RTEs. Further, irrespective of age, RTE loci near protein-coding coding genes are targeted to a greater extent than RTE loci far from genes. During development, a shift in piRNA targeting is observed, with a clear increase in the relative piRNA targeting of RTEs residing within boundaries of protein-coding gene transcripts.

Conclusions

Reanalyzing published piRNA sequences and taking into account the features of individual RTE loci provide novel insight into the activity of RTEs during development. The obtained results are consistent with some degree of proportionality between what transcripts become substrates for Piwi protein complexes and the level by which the transcripts are present in the cell. A transition from active transcription of RTEs to passive co-transcription of RTE sequences residing within protein-coding transcripts appears to take place in postnatal development. Hence, the previously reported increase in piRNA targeting of SINEs in postnatal testis development does not necessitate widespread active transcription of SINEs, but may simply be explained by the prevalence of SINEs residing in introns.

Dynamic interactions between transposable elements and their hosts

Transposable elements (TEs) have a unique ability to mobilize to new genomic locations, and the major advance of second-generation DNA sequencing has provided insights into the dynamic relationship between TEs and their hosts. It now is clear that TEs have adopted diverse strategies - such as specific integration sites or patterns of activity - to thrive in host environments that are replete with mechanisms, such as small RNAs or epigenetic marks, that combat TE amplification. Emerging evidence suggests that TE mobilization might sometimes benefit host genomes by enhancing genetic diversity, although TEs are also implicated in diseases such as cancer. Here, we discuss recent findings about how, where and when TEs insert in diverse organisms.

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.

Restless genomes humans as a model organism for understanding host-retrotransposable element dynamics

Since their initial discovery in maize, there have been various attempts to categorize the relationship between transposable elements (TEs) and their host organisms. These have ranged from TEs being selfish parasites to their role as essential, functional components of organismal biology. Research over the past several decades has, in many respects, only served to complicate the issue even further. On the one hand, investigators have amassed substantial evidence concerning the negative effects that TE-mutagenic activity can have on host genomes and organismal fitness. On the other hand, we find an increasing number of examples, across several taxa, of TEs being incorporated into functional biological roles for their host organism. Some 45% of our own genomes are comprised of TE copies. While many of these copies are dormant, having lost their ability to mobilize, several lineages continue to actively proliferate in modern human populations. With its complement of ancestral and active TEs, the human genome exhibits key aspects of the host-TE dynamic that has played out since early on in organismal evolution. In this review, we examine what insights the particularly well-characterized human system can provide regarding the nature of the host-TE interaction.

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.

IAP display: a simple method to identify mouse strain specific IAP insertions

Intracisternal A-type particle (IAP) elements are high copy number long terminal repeat (LTR) rodent retrotransposons. Some IAP elements can transpose, and are responsible for ~12% of spontaneous mouse mutations. Inbred mouse strains show variation in genomic IAP distribution, contributing to inter-strain genetic variability. Additionally IAP elements can influence the transcriptional regulation of neighbouring genes through their strong LTR promoter, effecting phenotypic variation. This genetic and phenotypic variability can translate into experimental variability between mouse strains. For example, it has been demonstrated that strain-specific genetic/epigenetic factors can interact to yield variable responses to drugs. Therefore, in experimental contexts it is essential to unequivocally identify mouse strains. Recently it was estimated that any two inbred strains share only ~40% of their IAP insertions. Of the remaining 60%, some insertions will be strain specific, fixed during inbreeding. These fixed insertions can be exploited as genetic markers to identify inbred strains, if they can be identified simply and efficiently. Here, we report the development of a PCR-based system allowing direct acquisition of strain-specific IAP insertions. In a pilot study, we identified 21 IAP loci, genotyped IAP insertions at 9 loci, and discovered two strain-specific insertions that could reliably identify these strains.

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.