Mammalian genomes contain numerous evolutionary harbored mobile elements, a part of which are still active and may cause genomic instability. Their movement and positional diversity occasionally result in phenotypic changes and variation by causing altered expression or disruption of neighboring host genes. Here, we describe a novel microarray-based method by which dispersed genomic locations of a type of retrotransposon in a mammalian genome can be identified. Using this method, we mapped the DNA elements for a mouse retrotransposon, intracisternal A-particle (IAP), within genomes of C3H/He and C57BL/6J inbred mouse strains; consequently we detected hundreds of probable IAP cDNA-integrated genomic regions, in which a considerable number of strain-specific putative insertions were included. In addition, by comparing genomic DNAs from radiation-induced myeloid leukemia cells and its reference normal tissue, we detected three genomic regions around which an IAP element was integrated. These results demonstrate the first successful genome-wide mapping of a retrotransposon type in a mammalian genome.

Stabilin2 (Stab2) encodes a large transmembrane protein which is predominantly expressed in the liver sinusoidal endothelial cells (LSECs) and functions as a scavenger receptor for various macromolecules including hyaluronans (HA). In DBA/2J mice, plasma HA concentration is ten times higher than in 129S6 or C57BL/6J mice, and this phenotype is genetically linked to the Stab2 locus. Stab2 mRNA in the LSECs was significantly lower in DBA/2J than in 129S6, leading to reduced STAB2 proteins in the DBA/2J LSECs. We found a retrovirus-derived transposable element, intracisternal A particle (IAP), in the promoter region of Stab2DBA which likely interferes with normal expression in the LSECs. In contrast, in other tissues of DBA/2J mice, the IAP drives high ectopic Stab2DBA transcription starting within the 5' long terminal repeat of IAP in a reverse orientation and continuing through the downstream Stab2DBA. Ectopic transcription requires the Stab2-IAP element but is dominantly suppressed by the presence of loci on 59.7-73.0 Mb of chromosome (Chr) 13 from C57BL/6J, while the same region in 129S6 requires additional loci for complete suppression. Chr13:59.9-73 Mb contains a large number of genes encoding Krüppel-associated box-domain zinc-finger proteins that target transposable elements-derived sequences and repress their expression. Despite the high amount of ectopic Stab2DBA transcript in tissues other than liver, STAB2 protein was undetectable and unlikely to contribute to the plasma HA levels of DBA/2J mice. Nevertheless, the IAP insertion and its effects on the transcription of the downstream Stab2DBA exemplify that stochastic evolutional events could significantly influence susceptibility to complex but common diseases.

Throughout evolution primate genomes have been modified by waves of retrotransposon insertions. For each wave, the host eventually finds a way to repress retrotransposon transcription and prevent further insertions. In mouse embryonic stem cells, transcriptional silencing of retrotransposons requires KAP1 (also known as TRIM28) and its repressive complex, which can be recruited to target sites by KRAB zinc-finger (KZNF) proteins such as murine-specific ZFP809 which binds to integrated murine leukaemia virus DNA elements and recruits KAP1 to repress them. KZNF genes are one of the fastest growing gene families in primates and this expansion is hypothesized to enable primates to respond to newly emerged retrotransposons. However, the identity of KZNF genes battling retrotransposons currently active in the human genome, such as SINE-VNTR-Alu (SVA) and long interspersed nuclear element 1 (L1), is unknown. Here we show that two primate-specific KZNF genes rapidly evolved to repress these two distinct retrotransposon families shortly after they began to spread in our ancestral genome. ZNF91 underwent a series of structural changes 8-12 million years ago that enabled it to repress SVA elements. ZNF93 evolved earlier to repress the primate L1 lineage until ∼12.5 million years ago when the L1PA3-subfamily of retrotransposons escaped ZNF93's restriction through the removal of the ZNF93-binding site. Our data support a model where KZNF gene expansion limits the activity of newly emerged retrotransposon classes, and this is followed by mutations in these retrotransposons to evade repression, a cycle of events that could explain the rapid expansion of lineage-specific KZNF genes.

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.

We have previously identified a novel atherosclerosis quantitative trait locus (QTL), Arch atherosclerosis 5 (Aath5), on mouse chromosome 10 by three-way QTL analyses between Apoe -/- mice on a DBA/2J, 129S6 and C57BL/6J background. The DBA/2J haplotype at the Aath5 locus was associated with smaller plaque size. One of the candidate genes underlying Aath5 was Stabilin-2 (Stab2), which encodes a clearance receptor for hyaluronan (HA) predominantly expressed in liver sinusoidal endothelial cells (LSECs). However, the role of Stab2 in atherosclerosis is unknown. A congenic line of Apoe -/- mice carrying Aath5 covering the Stab2 DBA allele on a background of 129S6 confirmed the small reductions of atherosclerotic plaque development. To further determine whether Stab2 is an underlying gene for Aath5, we generated Stab2 -/- Apoe -/- mice on a C57BL/6J background. When fed with a Western diet for 8 weeks, Stab2 -/- Apoe -/- males developed approximately 30% smaller plaques than Stab2 +/+ Apoe -/- mice. HA was accumulated in circulation but not in major organs in the Stab2 deficient mice. STAB2-binding molecules that are involved in atherosclerosis, including acLDL, apoptotic cells, heparin and vWF were not likely the direct cause of the protection in the Stab2 -/- Apoe -/- males. These data indicate that reduction of Stab2 is protective against atherosclerotic plaque development, and that Stab2 is a contributing gene underlying Aath5, although its effect is small. To test whether non-synonymous amino acid changes unique to DBA/2J affect the function of STAB2 protein, we made HEK293 cell lines expressing STAB2129 or STAB2DBA proteins, as well as STAB2129 proteins carrying each of five DBA-unique replacements that have been predicted to be deleterious. These mutant cells were capable of internalizing 125I -HA and DiI-acLDL similarly to the control cells. These results indicate that the amino acid changes unique to DBA/2J are not affecting the function of STAB2 protein, and support our previous observation that the reduced transcription of Stab2 in the liver sinusoid as a consequence of the insertion of a viral-derived sequence, intracisternal A particle, is the primary contributor to the athero-protection conferred by the DBA/2J allele.

DNA methylation is a heritable epigenetic modification involved in gene silencing, imprinting, and the suppression of retrotransposons. Global DNA demethylation occurs in the early embryo and the germ line, and may be mediated by Tet (ten eleven translocation) enzymes, which convert 5-methylcytosine (5mC) to 5-hydroxymethylcytosine (5hmC). Tet enzymes have been studied extensively in mouse embryonic stem (ES) cells, which are generally cultured in the absence of vitamin C, a potential cofactor for Fe(II) 2-oxoglutarate dioxygenase enzymes such as Tet enzymes. Here we report that addition of vitamin C to mouse ES cells promotes Tet activity, leading to a rapid and global increase in 5hmC. This is followed by DNA demethylation of many gene promoters and upregulation of demethylated germline genes. Tet1 binding is enriched near the transcription start site of genes affected by vitamin C treatment. Importantly, vitamin C, but not other antioxidants, enhances the activity of recombinant Tet1 in a biochemical assay, and the vitamin-C-induced changes in 5hmC and 5mC are entirely suppressed in Tet1 and Tet2 double knockout ES cells. Vitamin C has a stronger effect on regions that gain methylation in cultured ES cells compared to blastocysts, and in vivo are methylated only after implantation. In contrast, imprinted regions and intracisternal A particle retroelements, which are resistant to demethylation in the early embryo, are resistant to vitamin-C-induced DNA demethylation. Collectively, the results of this study establish vitamin C as a direct regulator of Tet activity and DNA methylation fidelity in ES cells.

Repetitive elements (REs) comprise 40-60% of the mammalian genome and have been shown to epigenetically influence the expression of genes through the formation of fusion transcript (FTs). We previously showed that an intracisternal A particle forms an FT with the agouti gene in mice, causing obesity/type 2 diabetes. To determine the frequency of FTs genome-wide, we developed a TopHat-Fusion-based analytical pipeline to identify FTs with high specificity. We applied it to an RNA-seq dataset from the nucleus accumbens (NAc) of mice repeatedly exposed to cocaine. Cocaine was previously shown to increase the expression of certain REs in this brain region. Using this pipeline that can be applied to single- or paired-end reads, we identified 438 genes expressing 813 different FTs in the NAc. Although all types of studied repeats were present in FTs, simple sequence repeats were underrepresented. Most importantly, reverse-transcription and quantitative PCR validated the expression of selected FTs in an independent cohort of animals, which also revealed that some FTs are the prominent isoforms expressed in the NAc by some genes. In other RNA-seq datasets, developmental expression as well as tissue specificity of some FTs differed from their corresponding non-fusion counterparts. Finally, in silico analysis predicted changes in the structure of proteins encoded by some FTs, potentially resulting in gain or loss of function. Collectively, these results indicate the robustness of our pipeline in detecting these new isoforms of genes, which we believe provides a valuable tool to aid in better understanding the broad role of REs in mammalian cellular biology.

In order to understand whether early epigenetic mechanisms instruct the long-term behavior of neural stem cells (NSCs) and their progeny, we examined Uhrf1 (ubiquitin-like PHD ring finger-1; also known as Np95), as it is highly expressed in NSCs of the developing brain and rapidly down-regulated upon differentiation. Conditional deletion of Uhrf1 in the developing cerebral cortex resulted in rather normal proliferation and neurogenesis but severe postnatal neurodegeneration. During development, deletion of Uhrf1 lead to global DNA hypomethylation with a strong activation of the intracisternal A particle (IAP) family of endogenous retroviral elements, accompanied by an increase in 5-hydroxymethylcytosine. Down-regulation of Tet enzymes rescued the IAP activation in Uhrf1 conditional knockout (cKO) cells, suggesting an antagonistic interplay between Uhrf1 and Tet on IAP regulation. As IAP up-regulation persists into postnatal stages in the Uhrf1 cKO mice, our data show the lack of means to repress IAPs in differentiating neurons that normally never express Uhrf1 The high load of viral proteins and other transcriptional deregulation ultimately led to postnatal neurodegeneration. Taken together, these data show that early developmental NSC factors can have long-term effects in neuronal differentiation and survival. Moreover, they highlight how specific the consequences of widespread changes in DNA methylation are for certain classes of retroviral elements.

Defective germline reprogramming in Piwil4 (Miwi2)- and Dnmt3l-deficient mice results in the failure to reestablish transposon silencing, meiotic arrest and progressive loss of spermatogonia. Here we sought to understand the molecular basis for this spermatogonial dysfunction. Through a combination of imaging, conditional genetics and transcriptome analysis, we demonstrate that germ cell elimination in the respective mutants arises as a result of defective de novo genome methylation during reprogramming rather than because of a function for the respective factors within spermatogonia. In both Miwi2-/- and Dnmt3l-/- spermatogonia, the intracisternal-A particle (IAP) family of endogenous retroviruses is derepressed, but, in contrast to meiotic cells, DNA damage is not observed. Instead, we find that unmethylated IAP promoters rewire the spermatogonial transcriptome by driving expression of neighboring genes. Finally, spermatogonial numbers, proliferation and differentiation are altered in Miwi2-/- and Dnmt3l-/- mice. In summary, defective reprogramming deregulates the spermatogonial transcriptome and may underlie spermatogonial dysfunction.

High serum levels of triglycerides (TG) and low levels of high-density lipoprotein cholesterol (HDL-C) increase the risk of coronary heart disease in humans. Herein, we first reported that the C3H/HeNSlc (C3H-S) mouse, a C3H/HeN-derived substrain, is a novel model for dyslipidemia. C3H-S showed hypertriglyceridemia and low total cholesterol (TC), HDL-C, and phospholipid (PL) concentrations. To identify the gene locus causing dyslipidemia in C3H-S, we performed genetic analysis. In F2 intercrosses between C3H-S mice and strains with normal serum lipids, the locus associated with serum lipids was identified as 163-168 Mb on chromosome 2. The phospholipid transfer protein (Pltp) gene was a candidate gene within this locus. Pltp expression and serum PLTP activity were markedly lower in C3H-S mice. Pltp expression was negatively correlated with serum TG and positively correlated with serum TC and HDL-C in F2 mice. Genome sequencing analysis revealed that an endogenous retrovirus (ERV) sequence called intracisternal A particle was inserted into intron 12 of Pltp in C3H-S. These results suggest that ERV insertion within Pltp causes aberrant splicing, leading to reduced Pltp expression in C3H-S. This study demonstrated the contribution of C3H-S to our understanding of the relationship between TG, TC, and PL metabolism via PLTP.

"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 ability of mammalian cytidine deaminases encoded by the APOBEC3 (A3) genes to restrict a broad number of endogenous retroelements and exogenous retroviruses, including murine leukemia virus and human immunodeficiency virus (HIV)-1, is now well established. The RNA editing family member apolipoprotein B (apo B)-editing catalytic subunit 1 (APOBEC1; A1) from a variety of mammalian species, a protein involved in lipid transport and which mediates C-U deamination of mRNA for apo B, has also been shown to modify a range of exogenous retroviruses, but its activity against endogenous retroelements remains unclear. Here, we show in cell culture-based retrotransposition assays that A1 family proteins from multiple mammalian species can also reduce the mobility and infectivity potential of LINE-1 (long interspersed nucleotide sequence-1, L1) and long-terminal repeats (LTRs) retrotransposons (or endogenous retroviruses), such as murine intracisternal A-particle (IAP) and MusD sequences. The anti-L1 activity of A1 was mainly mediated by a deamination-independent mechanism, and was not affected by subcellular localization of the proteins. In contrast, the inhibition of LTR-retrotransposons appeared to require the deaminase activity of A1 proteins. Thus, the AID/APOBEC family proteins including A1s employ multiple mechanisms to regulate the mobility of autonomous retrotransposons in several mammalian species.

Germline genome defense evolves to recognize and suppress retrotransposons. One of defensive mechanisms is the PIWI-associated RNA (piRNA) pathway, which employs small RNAs for sequence-specific repression. The loss of the piRNA pathway in mice causes male sterility while females remain fertile. Unlike spermatogenic cells, mouse oocytes posses also RNA interference (RNAi), another small RNA pathway capable of retrotransposon suppression. To examine whether RNAi compensates the loss of the piRNA pathway, we produced a new RNAi pathway mutant DicerSOM and crossed it with a catalytically-dead mutant of Mili, an essential piRNA gene. Normal follicular and oocyte development in double mutants showed that RNAi does not suppress a strong ovarian piRNA knock-out phenotype. However, we observed redundant and non-redundant targeting of specific retrotransposon families illustrating stochasticity of recognition and targeting of invading retrotransposons. Intracisternal A Particle retrotransposon was mainly targeted by the piRNA pathway, MaLR and RLTR10 retrotransposons were targeted mainly by RNAi. Double mutants showed accumulations of LINE-1 retrotransposon transcripts. However, we did not find strong evidence for transcriptional activation and mobilization of retrotransposition competent LINE-1 elements suggesting that while both defense pathways are simultaneously expendable for ovarian oocyte development, yet another transcriptional silencing mechanism prevents mobilization of LINE-1 elements.

Transposable elements (TEs) are mobile genetic elements that make up a large fraction of mammalian genomes. While select TEs have been co-opted in host genomes to have function, the majority of these elements are epigenetically silenced by DNA methylation in somatic cells. However, some TEs in mice, including the Intracisternal A-particle (IAP) subfamily of retrotransposons, have been shown to display interindividual variation in DNA methylation. Recent work has revealed that IAP sequence differences and strain-specific KRAB zinc finger proteins (KZFPs) may influence the methylation state of these IAPs. However, the mechanisms underlying the establishment and maintenance of interindividual variability in DNA methylation still remain unclear. Here, we report that sequence content and genomic context influence the likelihood that IAPs become variably methylated. IAPs that differ from consensus IAP sequences have altered KZFP recruitment that can lead to decreased KAP1 recruitment when in proximity of constitutively expressed genes. These variably methylated loci have a high CpG density, similar to CpG islands, and can be bound by ZF-CxxC proteins, providing a potential mechanism to maintain this permissive chromatin environment and protect from DNA methylation. These observations indicate that variably methylated IAPs escape silencing through both attenuation of KZFP binding and recognition by ZF-CxxC proteins to maintain a hypomethylated state.

Retrotransposons compose a staggering 40% of the mammalian genome. Among them, endogenous retroviruses (ERV) represent sequences that closely resemble the proviruses created from exogenous retroviral infection. ERVs make up 8 to 10% of human and mouse genomes and range from evolutionarily ancient sequences to recent acquisitions. Studies in Drosophila have provided a causal link between genomic retroviral elements and cognitive decline; however, in mammals, the role of ERVs in learning and memory remains unclear. Here we studied 2 independent murine models for ERV activation: muMT strain (lacking B cells and antibody production) and intracerebroventricular injection of streptozotocin (ICVI-STZ). We conducted behavioral assessments (contextual fear memory and spatial learning), as well as gene and protein analysis (RNA sequencing, PCR, immunohistochemistry, and western blot assays). Mice lacking mitochondrial antiviral-signaling protein (MAVS) and mice lacking stimulator of IFN genes protein (STING), 2 downstream sensors of ERV activation, provided confirmation of ERV impact. We found that muMT mice and ICVI-STZ mice induced hippocampal ERV activation, as shown by increased gene and protein expression of the Gag sequence of the transposable element intracisternal A-particle. ERV activation was accompanied by significant hippocampus-related memory impairment in both models. Notably, the deficiency of the MAVS pathway was protective against ICVI-STZ-induced cognitive pathology. Overall, our results demonstrate that ERV activation is associated with cognitive impairment in mice. Moreover, they provide a molecular target for strategies aimed at attenuating retroviral element sensing, via MAVS, to treat dementia and neuropsychiatric disorders.

Xenotropic murine leukemia virus-related virus (XMRV) is a newly identified gamma retrovirus and may be associated with prostate cancer- (PC) and chronic fatigue syndrome (CFS). Since its identification in 2006 and detection of polytropic murine lenkemia virus (MLV)-like sequences in CFS patients in 2010, several test methods including nucleic acid testing methods and serological assays have been developed for detection of XMRV and/or MLV-like sequences. However, these research assays have not yet been validated and evaluated due to the lack of well-characterized reference materials. Mouse DNA contamination should be carefully checked when testing human specimens in order to avoid false-positive detection of XMRV or MLV-like sequences.

Obesity is one of the largest health problems facing the world today. Although twin and family studies suggest about two-thirds of obesity is caused by genetic factors, only a small fraction of this variance has been unraveled. There are still large numbers of genes to be identified that cause variations in body fatness and the associated diseases encompassed in the metabolic syndrome (MetS). A locus near a sequence tagged site (STS) marker D6S1009 has been linked to obesity or body mass index (BMI). However, its genetic entity is unknown. D6S1009 is located in the intergenic region between SLC35D3 and NHEG1. Here we report that the ros mutant mice harboring a recessive mutation in the Slc35d3 gene show obesity and MetS and reduced membrane dopamine receptor D1 (D1R) with impaired dopamine signaling in striatal neurons. SLC35D3 is localized to both endoplasmic reticulum (ER) and early endosomes and interacts with D1R. In ros striatal D1 neurons, lack of SLC35D3 causes the accumulation of D1R on the ER to impair its ER exit. The MetS phenotype is reversible by the administration of D1R agonist to the ros mutant. In addition, we identified two mutations in the SLC35D3 gene in patients with MetS, which alter the subcellular localization of SLC35D3. Our results suggest that the SLC35D3 gene, close to the D6S1009 locus, is a candidate gene for MetS, which is involved in metabolic control in the central nervous system by regulating dopamine signaling.