Extensive variation between inbred mouse strains due to endogenous L1 retrotransposition

Locus-resolution analysis of L1 regulation and retrotransposition potential in mouse embryonic development

Mice harbor ∼2800 intact copies of the retrotransposon Long Interspersed Element 1 (L1). The in vivo retrotransposition capacity of an L1 copy is defined by both its sequence integrity and epigenetic status, including DNA methylation of the monomeric units constituting young mouse L1 promoters. Locus-specific L1 methylation dynamics during development may therefore elucidate and explain spatiotemporal niches of endogenous retrotransposition but remain unresolved. Here, we interrogate the retrotransposition efficiency and epigenetic fate of source (donor) L1s, identified as mobile in vivo. We show that promoter monomer loss consistently attenuates the relative retrotransposition potential of their offspring (daughter) L1 insertions. We also observe that most donor/daughter L1 pairs are efficiently methylated upon differentiation in vivo and in vitro. We use Oxford Nanopore Technologies (ONT) long-read sequencing to resolve L1 methylation genome-wide and at individual L1 loci, revealing a distinctive "smile" pattern in methylation levels across the L1 promoter region. Using Pacific Biosciences (PacBio) SMRT sequencing of L1 5' RACE products, we then examine DNA methylation dynamics at the mouse L1 promoter in parallel with transcription start site (TSS) distribution at locus-specific resolution. Together, our results offer a novel perspective on the interplay between epigenetic repression, L1 evolution, and genome stability.

Drug discovery oncology in a mouse: concepts, models and limitations

The utilization of suitable mouse models is a critical step in the drug discovery oncology workflow as their generation and use are important for target identification and validation as well as toxicity and efficacy assessments. Current murine models have been instrumental in furthering insights into the mode of action of drugs before transitioning into the clinic. Recent advancements in genome editing with the development of the CRISPR/Cas9 system and the possibility of applying such technology directly in vivo have expanded the toolkit of preclinical models available. In this review, a brief presentation of the current models used in drug discovery will be provided with a particular emphasis on the novel CRISPR/Cas9 models.

Organ-, sex- and age-dependent patterns of endogenous L1 mRNA expression at a single locus resolution

Expression of L1 mRNA, the first step in the L1 copy-and-paste amplification cycle, is a prerequisite for L1-associated genomic instability. We used a reported stringent bioinformatics method to parse L1 mRNA transcripts and measure the level of L1 mRNA expressed in mouse and rat organs at a locus-specific resolution. This analysis determined that mRNA expression of L1 loci in rodents exhibits striking organ specificity with less than 0.8% of loci shared between organs of the same organism. This organ specificity in L1 mRNA expression is preserved in male and female mice and across age groups. We discovered notable differences in L1 mRNA expression between sexes with only 5% of expressed L1 loci shared between male and female mice. Moreover, we report that the levels of total L1 mRNA expression and the number and spectrum of expressed L1 loci fluctuate with age as independent variables, demonstrating different patterns in different organs and sexes. Overall, our comparisons between organs and sexes and across ages ranging from 2 to 22 months establish previously unforeseen dynamic changes in L1 mRNA expression in vivo. These findings establish the beginning of an atlas of endogenous L1 mRNA expression across a broad range of biological variables that will guide future studies.

Distinct contributions of DNA methylation and histone acetylation to the genomic occupancy of transcription factors

Epigenetic modifications on chromatin play important roles in regulating gene expression. Although chromatin states are often governed by multilayered structure, how individual pathways contribute to gene expression remains poorly understood. For example, DNA methylation is known to regulate transcription factor binding but also to recruit methyl-CpG binding proteins that affect chromatin structure through the activity of histone deacetylase complexes (HDACs). Both of these mechanisms can potentially affect gene expression, but the importance of each, and whether these activities are integrated to achieve appropriate gene regulation, remains largely unknown. To address this important question, we measured gene expression, chromatin accessibility, and transcription factor occupancy in wild-type or DNA methylation-deficient mouse embryonic stem cells following HDAC inhibition. We observe widespread increases in chromatin accessibility at retrotransposons when HDACs are inhibited, and this is magnified when cells also lack DNA methylation. A subset of these elements has elevated binding of the YY1 and GABPA transcription factors and increased expression. The pronounced additive effect of HDAC inhibition in DNA methylation-deficient cells demonstrates that DNA methylation and histone deacetylation act largely independently to suppress transcription factor binding and gene expression.

Subfamily-specific quantification of endogenous mouse L1 retrotransposons by droplet digital PCR

Long interspersed element type 1 (LINE-1; L1) mobilizes during early embryogenesis, neurogenesis, and germ cell development, accounting for 25% of disease-causing heritable insertions and 98% of somatic insertions in cancer. To better understand the regulation and impact of L1 mobilization in the genome, reliable methods for measuring L1 copy number variation (CNV) are needed. Here we present a comprehensive analysis of a droplet digital PCR (ddPCR) based method for quantifying endogenous mouse L1. We provide experimental evidence that ddPCR assays can be designed to target specific L1 subfamilies using diagnostic single nucleotide polymorphisms (SNPs). The target and off-target L1 subfamilies form distinct droplet clusters, which were experimentally verified using both synthetic gene fragments and endogenous L1 derived plasmid clones. We further provide a roadmap for in silico assay design and evaluation of target specificity, ddPCR testing, and optimization for L1 CNV quantification. The assay can achieve a sensitivity of 5% CNV with 8 technical replicates. With 24 technical replicates, it can detect 2% CNV because of the increased precision. The same approach will serve as a guide for the development of ddPCR based assays for quantifying human L1 copy number and any other high copy genomic target sequences.

The Parkinson's disease-associated kinase LRRK2 regulates genes required for cell adhesion, polarization, and chemotaxis in activated murine macrophages

Leucine-rich repeat kinase 2 (LRRK2) encodes a complex protein that includes kinase and GTPase domains. Genome-wide association studies have identified dominant LRRK2 alleles that predispose their carriers to late-onset idiotypic Parkinson's disease (PD) and also to autoimmune disorders such as Crohn's disease. Considerable evidence indicates that PD initiation and progression involve activation of innate immune functions in microglia, which are brain-resident macrophages. Here we asked whether LRRK2 modifies inflammatory signaling and how this modification might contribute to PD and Crohn's disease. We used RNA-Seq-based high-resolution transcriptomics to compare gene expression in activated primary macrophages derived from WT and Lrrk2 knockout mice. Remarkably, expression of a single gene, Rap guanine nucleotide exchange factor 3 (Rapgef3), was strongly up-regulated in the absence of LRRK2 and down-regulated in its presence. We observed similar regulation of Rapgef3 expression in cells treated with a highly specific inhibitor of LRRK2 protein kinase activity. Rapgef3 encodes an exchange protein, activated by cAMP 1 (EPAC-1), a guanine nucleotide exchange factor that activates the small GTPase Rap-1. Rap-1 mediates cell adhesion, polarization, and directional motility, and our results indicate that LRRK2 modulates chemotaxis of microglia and macrophages. Dominant PD-associated LRRK2 alleles may suppress EPAC-1 activity, further restricting motility and preventing efficient migration of microglia to sites of neuronal damage. Functional analysis in vivo in a subclinical infection model also indicated that Lrrk2 subtly modifies the inflammatory response. These results indicate that LRRK2 modulates the expression of genes involved in murine immune cell chemotaxis.

Prion protein PrP nucleic acid binding and mobilization implicates retroelements as the replicative component of transmissible spongiform encephalopathy

The existence of more than 30 strains of transmissible spongiform encephalopathy (TSE) and the paucity of infectivity of purified PrP^Sc, as well as considerations of PrP structure, are inconsistent with the protein-only (prion) theory of TSE. Nucleic acid is a strong contender as a second component. We juxtapose two key findings: (i) PrP is a nucleic-acid-binding antimicrobial protein that is similar to retroviral Gag proteins in its ability to trigger reverse transcription. (ii) Retroelement mobilization is widely seen in TSE disease. Given further evidence that PrP also mediates nucleic acid transport into and out of the cell, a strong case is to be made that a second element – retroelement nucleic acid – bound to PrP constitutes the second component necessary to explain the multiple strains of TSE.

Insertion of a chimeric retrotransposon sequence in mouse Axin1 locus causes metastable kinky tail phenotype

Background

Transposable elements (TEs) make up > 50% of the human genome, and the majority of retrotransposon insertions are truncated and many are located in introns. However, the effects of retrotransposition on the host genes remain incompletely known.

Results

We report here that insertion of a chimeric L1 (cL1), but not IAP solo LTR, into intron 6 of Axin1 using CRIPSR/Cas9 induced the kinky tail phenotype with ~ 80% penetrance in heterozygous Axin ^cL1 mice. Both penetrant (with kinky tails) and silent (without kinky tails) Axin ^cL1 mice, regardless of sex, could transmit the phenotype to subsequent generations with similar penetrance (~ 80%). Further analyses revealed that a longer Axin1 transcript isoform containing partial cL1-targeted intron was present in penetrant, but absent in silent and wild type mice, and the production of this unique Axin1 transcript appeared to correlate with altered levels of an activating histone modification, H3K9ac.

Conclusions

The mechanism for Axin ^cL1 mice is different from those previously identified in mice with spontaneous retrotransposition of IAP, e.g., Axin ^Fu and A ^vy , both of which have been associated with DNA methylation changes. Our data suggest that Axin1 locus is sensitive to genetic and epigenetic alteration by retrotransposons and thus, ideally suited for studying the effects of new retrotransposition events on target gene function in mice.

Mouse germ line mutations due to retrotransposon insertions

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

MobiSeq: De novo SNP discovery in model and non-model species through sequencing the flanking region of transposable elements

In recent years, the availability of reduced representation library (RRL) methods has catalysed an expansion of genome-scale studies to characterize both model and non-model organisms. Most of these methods rely on the use of restriction enzymes to obtain DNA sequences at a genome-wide level. These approaches have been widely used to sequence thousands of markers across individuals for many organisms at a reasonable cost, revolutionizing the field of population genomics. However, there are still some limitations associated with these methods, in particular the high molecular weight DNA required as starting material, the reduced number of common loci among investigated samples, and the short length of the sequenced site-associated DNA. Here, we present MobiSeq, a RRL protocol exploiting simple laboratory techniques, that generates genomic data based on PCR targeted enrichment of transposable elements and the sequencing of the associated flanking region. We validate its performance across 103 DNA extracts derived from three mammalian species: grey wolf (Canis lupus), red deer complex (Cervus sp.) and brown rat (Rattus norvegicus). MobiSeq enables the sequencing of hundreds of thousands loci across the genome and performs SNP discovery with relatively low rates of clonality. Given the ease and flexibility of MobiSeq protocol, the method has the potential to be implemented for marker discovery and population genomics across a wide range of organisms-enabling the exploration of diverse evolutionary and conservation questions.

L1 retrotransposition in the soma: a field jumping ahead

Retrotransposons are transposable elements (TEs) capable of "jumping" in germ, embryonic and tumor cells and, as is now clearly established, in the neuronal lineage. Mosaic TE insertions form part of a broader landscape of somatic genome variation and hold significant potential to generate phenotypic diversity, in the brain and elsewhere. At present, the LINE-1 (L1) retrotransposon family appears to be the most active autonomous TE in most mammals, based on experimental data obtained from disease-causing L1 mutations, engineered L1 reporter systems tested in cultured cells and transgenic rodents, and single-cell genomic analyses. However, the biological consequences of almost all somatic L1 insertions identified thus far remain unknown. In this review, we briefly summarize the current state-of-the-art in the field, including estimates of L1 retrotransposition rate in neurons. We bring forward the hypothesis that an extensive subset of retrotransposition-competent L1s may be de-repressed and mobile in the soma but largely inactive in the germline. We discuss recent reports of non-canonical L1-associated sequence variants in the brain and propose that the elevated L1 DNA content reported in several neurological disorders may predominantly comprise accumulated, unintegrated L1 nucleic acids, rather than somatic L1 insertions. Finally, we consider the main objectives and obstacles going forward in elucidating the biological impact of somatic retrotransposition.

Heritable L1 Retrotransposition Events During Development: Understanding Their Origins: Examination of heritable, endogenous L1 retrotransposition in mice opens up exciting new questions and research directions

The retrotransposon Long Interspersed Element 1 (LINE-1 or L1) has played a major role in shaping the sequence composition of the mammalian genome. In our recent publication, "Heritable L1 retrotransposition in the mouse primordial germline and early embryo," we systematically assessed the rate and developmental timing of de novo, heritable endogenous L1 insertions in mice. Such heritable retrotransposition events allow L1 to exert an ongoing influence upon genome evolution. Here, we place our findings in the context of earlier studies, and highlight how our results corroborate, and depart from, previous research based on human patient samples and transgenic mouse models harboring engineered L1 reporter genes. In parallel, we outline outstanding questions regarding the stage-specificity, regulation, and functional impact of embryonic and germline L1 retrotransposition, and propose avenues for future research in this field.

L1 retrotransposition is a common feature of mammalian hepatocarcinogenesis

The retrotransposon Long Interspersed Element 1 (LINE-1 or L1) is a continuing source of germline and somatic mutagenesis in mammals. Deregulated L1 activity is a hallmark of cancer, and L1 mutagenesis has been described in numerous human malignancies. We previously employed retrotransposon capture sequencing (RC-seq) to analyze hepatocellular carcinoma (HCC) samples from patients infected with hepatitis B or hepatitis C virus and identified L1 variants responsible for activating oncogenic pathways. Here, we have applied RC-seq and whole-genome sequencing (WGS) to an Abcb4 (Mdr2)-/- mouse model of hepatic carcinogenesis and demonstrated for the first time that L1 mobilization occurs in murine tumors. In 12 HCC nodules obtained from 10 animals, we validated four somatic L1 insertions by PCR and capillary sequencing, including TF subfamily elements, and one GF subfamily example. One of the TF insertions carried a 3' transduction, allowing us to identify its donor L1 and to demonstrate that this full-length TF element retained retrotransposition capacity in cultured cancer cells. Using RC-seq, we also identified eight tumor-specific L1 insertions from 25 HCC patients with a history of alcohol abuse. Finally, we used RC-seq and WGS to identify three tumor-specific L1 insertions among 10 intra-hepatic cholangiocarcinoma (ICC) patients, including one insertion traced to a donor L1 on Chromosome 22 known to be highly active in other cancers. This study reveals L1 mobilization as a common feature of hepatocarcinogenesis in mammals, demonstrating that the phenomenon is not restricted to human viral HCC etiologies and is encountered in murine liver tumors.

A fine LINE-1 in mouse embryonic chromatin regulation

The functional role of repetitive elements in mammalian genomes is still largely unexplored. A new study provides evidence that LINE-1 retrotransposons regulate chromatin dynamics and are essential for normal embryonic development in mice.

Expression dynamics and relations with nearby genes of rat transposable elements across 11 organs, 4 developmental stages and both sexes

BACKGROUND

TEs pervade mammalian genomes. However, compared with mice, fewer studies have focused on the TE expression patterns in rat, particularly the comparisons across different organs, developmental stages and sexes. In addition, TEs can influence the expression of nearby genes. The temporal and spatial influences of TEs remain unclear yet.

RESULTS

To evaluate the TEs transcription patterns, we profiled their transcript levels in 11 organs for both sexes across four developmental stages of rat. The results show that most short interspersed elements (SINEs) are commonly expressed in all conditions, which are also the major TE types with commonly expression patterns. In contrast, long terminal repeats (LTRs) are more likely to exhibit specific expression patterns. The expression tendency of TEs and genes are similar in most cases. For example, few specific genes and TEs are in the liver, muscle and heart. However, TEs perform superior over genes on classing organ, which imply their higher organ specificity than genes. By associating the TEs with the closest genes in genome, we find their expression levels are correlated, independent of their distance in some cases.

CONCLUSIONS

TEs sex-dependently associate with nearest genes. A gene would be associated with more than one TE. Our works can help to functionally annotate the genome and further understand the role of TEs in gene regulation.

LINE-1 activation after fertilization regulates global chromatin accessibility in the early mouse embryo

After fertilization, to initiate development, gametes are reprogramed to become totipotent. Approximately half of the mammalian genome consists of repetitive elements, including retrotransposons, some of which are transcribed after fertilization. Retrotransposon activation is generally assumed to be a side effect of the extensive chromatin remodeling underlying the epigenetic reprogramming of gametes. Here, we used a targeted epigenomic approach to address whether specific retrotransposon families play a direct role in chromatin organization and developmental progression. We demonstrate that premature silencing of LINE-1 elements decreases chromatin accessibility, whereas prolonged activation prevents the gradual chromatin compaction that occurs naturally in developmental progression. Preventing LINE-1 activation and interfering with its silencing decreases developmental rates independently of the coding nature of the LINE-1 transcript, thus suggesting that LINE-1 functions primarily at the chromatin level. Our data suggest that activation of LINE-1 regulates global chromatin accessibility at the beginning of development and indicate that retrotransposon activation is integral to the developmental program.

Dynamic silencing of somatic L1 retrotransposon insertions reflects the developmental and cellular contexts of their genomic integration

Background

The ongoing mobilization of mammalian transposable elements (TEs) contributes to natural genetic variation. To survey the epigenetic control and expression of reporter genes inserted by L1 retrotransposition in diverse cellular and genomic contexts, we engineered highly sensitive, real-time L1 retrotransposon reporter constructs.

Results

Here we describe different patterns of expression and epigenetic controls of newly inserted sequences retrotransposed by L1 in various somatic cells and tissues including cultured human cancer cells, mouse embryonic stem cells, and tissues of pseudofounder transgenic mice and their progeny. In cancer cell lines, the newly inserted sequences typically underwent rapid transcriptional gene silencing, but they lacked cytosine methylation even after many cell divisions. L1 reporter expression was reversible and oscillated frequently. Silenced or variegated reporter expression was strongly and uniformly reactivated by treatment with inhibitors of histone deacetylation, revealing the mechanism for their silencing. By contrast, de novo integrants retrotransposed by L1 in pluripotent mouse embryonic stem (ES) cells underwent rapid silencing by dense cytosine methylation. Similarly, de novo cytosine methylation also was identified at new integrants when studied in several distinct somatic tissues of adult founder mice. Pre-existing L1 elements in cultured human cancer cells were stably silenced by dense cytosine methylation, whereas their transcription modestly increased when cytosine methylation was experimentally reduced in cells lacking DNA methyltransferases DNMT1 and DNMT3b. As a control, reporter genes mobilized by piggyBac (PB), a DNA transposon, revealed relatively stable and robust expression without apparent silencing in both cultured cancer cells and ES cells.

Conclusions

We hypothesize that the de novo methylation marks at newly inserted sequences retrotransposed by L1 in early pre-implantation development are maintained or re-established in adult somatic tissues. By contrast, histone deacetylation reversibly silences L1 reporter insertions that had mobilized at later timepoints in somatic development and differentiation, e.g., in cancer cell lines. We conclude that the cellular contexts of L1 retrotransposition can determine expression or silencing of newly integrated sequences. We propose a model whereby reporter expression from somatic TE insertions reflects the timing, molecular mechanism, epigenetic controls and the genomic, cellular and developmental contexts of their integration.

Electronic supplementary material

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

Heritable L1 retrotransposition in the mouse primordial germline and early embryo

LINE-1 (L1) retrotransposons are a noted source of genetic diversity and disease in mammals. To expand its genomic footprint, L1 must mobilize in cells that will contribute their genetic material to subsequent generations. Heritable L1 insertions may therefore arise in germ cells and in pluripotent embryonic cells, prior to germline specification, yet the frequency and predominant developmental timing of such events remain unclear. Here, we applied mouse retrotransposon capture sequencing (mRC-seq) and whole-genome sequencing (WGS) to pedigrees of C57BL/6J animals, and uncovered an L1 insertion rate of ≥1 event per eight births. We traced heritable L1 insertions to pluripotent embryonic cells and, strikingly, to early primordial germ cells (PGCs). New L1 insertions bore structural hallmarks of target-site primed reverse transcription (TPRT) and mobilized efficiently in a cultured cell retrotransposition assay. Together, our results highlight the rate and evolutionary impact of heritable L1 retrotransposition and reveal retrotransposition-mediated genomic diversification as a fundamental property of pluripotent embryonic cells in vivo.

Evolution and Diversity of Transposable Elements in Vertebrate Genomes

Transposable elements (TEs) are selfish genetic elements that mobilize in genomes via transposition or retrotransposition and often make up large fractions of vertebrate genomes. Here, we review the current understanding of vertebrate TE diversity and evolution in the context of recent advances in genome sequencing and assembly techniques. TEs make up 4-60% of assembled vertebrate genomes, and deeply branching lineages such as ray-finned fishes and amphibians generally exhibit a higher TE diversity than the more recent radiations of birds and mammals. Furthermore, the list of taxa with exceptional TE landscapes is growing. We emphasize that the current bottleneck in genome analyses lies in the proper annotation of TEs and provide examples where superficial analyses led to misleading conclusions about genome evolution. Finally, recent advances in long-read sequencing will soon permit access to TE-rich genomic regions that previously resisted assembly including the gigantic, TE-rich genomes of salamanders and lungfishes.

Distinct behavioral and epileptic phenotype differences in 129/P mice compared to C57BL/6 mice subject to intraamygdala kainic acid-induced status epilepticus

Animal models of status epilepticus are important tools to understand the pathogenesis of epileptic brain injury and evaluate potential seizure-suppressive, neuroprotective, and antiepileptogenic treatments. Focal elicitation of status epilepticus by intraamygdala kainic acid in mice produces unilateral hippocampal damage and the emergence of spontaneous recurrent seizures after a short latent period. The model has been characterized in C57BL/6, BALB/c, and SJL mice where strain-specific differences were found in the extent of hippocampal damage. 129/P mice are a common background strain for genetic models and may display unique characteristics in this model. We therefore compared responses to intraamygdala kainic acid between 129/P and C57BL/6 mice. Racine scale-scored convulsive behavior during status epilepticus was substantially lower in 129/P mice compared with that in C57BL/6 mice. Analysis of surface-recorded electroencephalogram (EEG) showed differences between strains in several frequency bands; EEG total power was greater during ictal episodes while duration of seizures was slightly shorter in 129/P mice. Histological analysis revealed similar hippocampal injury between strains, with neuronal death mainly confined to the ipsilateral CA3 subfield. Expression of genes associated with gliosis and neuroinflammatory responses was also similar between strains after seizures. Video-EEG telemetry recordings showed that 129/P mice first display spontaneous seizures within a few days of status epilepticus similar to C57BL/6 mice. However, high mortality in 129/P mice prevented a quantitative comparison of the epileptic seizure phenotypes between strains. This study defined behavioral, EEG, and histopathologic features of this mouse strain in a model increasingly useful for the study of the genetic contribution to acquired epilepsy. Intraamygdala kainic acid in 129/P mice could serve as a model of nonconvulsive status epilepticus, but long-term assessments will require model adjustment to mitigate the severity of the emergent epileptic phenotype.

LINEs in mice: features, families, and potential roles in early development

Approximately half of the mammalian genome is composed of repetitive elements, including LINE-1 (L1) elements. Because of their potential ability to transpose and integrate into other regions of the genome, their activation represents a threat to genome stability. Molecular pathways have emerged to tightly regulate and repress their transcriptional activity, including DNA methylation, histone modifications, and RNA pathways. It has become evident that Line-L1 elements are evolutionary diverse and dedicated repression pathways have been recently uncovered that discriminate between evolutionary old and young elements, with RNA-directed silencing mechanisms playing a prominent role. During periods of epigenetic reprogramming in development, specific classes of repetitive elements are upregulated, presumably due to the loss of most heterochromatic marks in this process. While we have learnt a lot on the molecular mechanisms that regulate Line-L1 expression over the last years, it is still unclear whether reactivation of Line-L1 after fertilization serves a functional purpose or it is a simple side effect of reprogramming.

Effect of crossing C57BL/6 and FVB mouse strains on basal cytokine expression

C57BL/6 is the most often used laboratory mouse strain. However, sometimes it is beneficial to cross the transgenic mice on the C57BL/6 background to the other strain, such as FVB. Although this is a common strategy, the influence of crossing these different strains on homeostatic expression of cytokines is not known. Here we have investigated the differences in the expression of selected cytokines between C57BL/6J and C57BL/6JxFVB mice in serum and skeletal muscle. We have found that only few cytokines were altered by crossing of the strains. Concentrations of IL5, IL7, LIF, MIP-2, and IP-10 were higher in serum of C57BL/6J mice than in C57BL/6JxFVB mice, whereas concentration of G-CSF was lower in C57BL/6J. In the skeletal muscle only the concentration of VEGF was higher in C57BL/6J mice than in C57BL/6JxFVB mice. Concluding, the differences in cytokine expression upon crossing C57BL/6 and FVB strain in basal conditions are not profound.

LINE-1 retrotransposons: from 'parasite' sequences to functional elements

Long interspersed nuclear elements-1 (LINE-1) are the most abundant and active retrotransposons in the mammalian genomes. Traditionally, the occurrence of LINE-1 sequences in the genome of mammals has been explained by the selfish DNA hypothesis. Nevertheless, recently, it has also been argued that these sequences could play important roles in these genomes, as in the regulation of gene expression, genome modelling and X-chromosome inactivation. The non-random chromosomal distribution is a striking feature of these retroelements that somehow reflects its functionality. In the present study, we have isolated and analysed a fraction of the open reading frame 2 (ORF2) LINE-1 sequence from three rodent species, Cricetus cricetus, Peromyscus eremicus and Praomys tullbergi. Physical mapping of the isolated sequences revealed an interspersed longitudinal AT pattern of distribution along all the chromosomes of the complement in the three genomes. A detailed analysis shows that these sequences are preferentially located in the euchromatic regions, although some signals could be detected in the heterochromatin. In addition, a coincidence between the location of imprinted gene regions (as Xist and Tsix gene regions) and the LINE-1 retroelements was also observed. According to these results, we propose an involvement of LINE-1 sequences in different genomic events as gene imprinting, X-chromosome inactivation and evolution of repetitive sequences located at the heterochromatic regions (e.g. satellite DNA sequences) of the rodents' genomes analysed.

Reduced pachytene piRNAs and translation underlie spermiogenic arrest in Maelstrom mutant mice

Pachytene piRNAs are a class of Piwi-interacting small RNAs abundant in spermatids of the adult mouse testis. They are processed from piRNA primary transcripts by a poorly understood mechanism and, unlike fetal transposon-derived piRNAs, lack complementary targets in the spermatid transcriptome. We report that immunopurified complexes of a conserved piRNA pathway protein Maelstrom (MAEL) are enriched in MIWI (Piwi partner of pachytene piRNAs), Tudor-domain proteins and processing intermediates of pachytene piRNA primary transcripts. We provide evidence of functional significance of these complexes in Mael129 knockout mice that exhibit spermiogenic arrest with acrosome and flagellum malformation. Mael129-null mutant testes possess low levels of piRNAs derived from MAEL-associated piRNA precursors and exhibit reduced translation of numerous spermiogenic mRNAs including those encoding acrosome and flagellum proteins. These translation defects in haploid round spermatids are likely indirect, as neither MAEL nor piRNA precursors associate with polyribosomes, and they may arise from an imbalance between pachytene piRNAs and MIWI.

Identification of structural variation in mouse genomes

Structural variation is variation in structure of DNA regions affecting DNA sequence length and/or orientation. It generally includes deletions, insertions, copy-number gains, inversions, and transposable elements. Traditionally, the identification of structural variation in genomes has been challenging. However, with the recent advances in high-throughput DNA sequencing and paired-end mapping (PEM) methods, the ability to identify structural variation and their respective association to human diseases has improved considerably. In this review, we describe our current knowledge of structural variation in the mouse, one of the prime model systems for studying human diseases and mammalian biology. We further present the evolutionary implications of structural variation on transposable elements. We conclude with future directions on the study of structural variation in mouse genomes that will increase our understanding of molecular architecture and functional consequences of structural variation.

Lizards and LINEs: selection and demography affect the fate of L1 retrotransposons in the genome of the green anole (Anolis carolinensis)

Autonomous retrotransposons lacking long terminal repeats (LTR) account for much of the variation in genome size and structure among vertebrates. Mammalian genomes contain hundreds of thousands of non-LTR retrotransposon copies, mostly resulting from the amplification of a single clade known as L1. The genomes of teleost fish and squamate reptiles contain a much more diverse array of non-LTR retrotransposon families, whereas copy number is relatively low. The majority of non-LTR retrotransposon insertions in nonmammalian vertebrates also appear to be very recent, suggesting strong purifying selection limits the accumulation of non-LTR retrotransposon copies. It is however unclear whether this turnover model, originally proposed in Drosophila, applies to nonmammalian vertebrates. Here, we studied the population dynamics of L1 in the green anole lizard (Anolis carolinensis). We found that although most L1 elements are recent in this genome, truncated insertions accumulate readily, and many are fixed at both the population and species level. In contrast, full-length L1 insertions are found at lower population frequencies, suggesting that the turnover model only applies to longer L1 elements in Anolis. We also found that full-length L1 inserts are more likely to be fixed in populations of small effective size, suggesting that the strength of purifying selection against deleterious alleles is highly dependent on host demographic history. Similar mechanisms seem to be controlling the fate of non-LTR retrotransposons in both Anolis and teleostean fish, which suggests that mammals have considerably diverged from the ancestral vertebrate in terms of how they interact with their intragenomic parasites.

A single whole-body low dose X-irradiation does not affect L1, B1 and IAP repeat element DNA methylation longitudinally

The low dose radioadaptive response has been shown to be protective against high doses of radiation as well as aging-induced genomic instability. We hypothesised that a single whole-body exposure of low dose radiation would induce a radioadaptive response thereby reducing or abrogating aging-related changes in repeat element DNA methylation in mice. Following sham or 10 mGy X-irradiation, serial peripheral blood sampling was performed and differences in Long Interspersed Nucleic Element 1 (L1), B1 and Intracisternal-A-Particle (IAP) repeat element methylation between samples were assessed using high resolution melt analysis of PCR amplicons. By 420 days post-irradiation, neither radiation- or aging-related changes in the methylation of peripheral blood, spleen or liver L1, B1 and IAP elements were observed. Analysis of the spleen and liver tissues of cohorts of untreated aging mice showed that the 17-19 month age group exhibited higher repeat element methylation than younger or older mice, with no overall decline in methylation detected with age. This is the first temporal analysis of the effect of low dose radiation on repeat element methylation in mouse peripheral blood and the first to examine the long term effect of this dose on repeat element methylation in a radiosensitive tissue (spleen) and a tissue fundamental to the aging process (liver). Our data indicate that the methylation of murine DNA repeat elements can fluctuate with age, but unlike human studies, do not demonstrate an overall aging-related decline. Furthermore, our results indicate that a low dose of ionising radiation does not induce detectable changes to murine repeat element DNA methylation in the tissues and at the time-points examined in this study. This radiation dose is relevant to human diagnostic radiation exposures and suggests that a dose of 10 mGy X-rays, unlike high dose radiation, does not cause significant short or long term changes to repeat element or global DNA methylation.

An antisense promoter in mouse L1 retrotransposon open reading frame-1 initiates expression of diverse fusion transcripts and limits retrotransposition

Between 6 and 30% of human and mouse transcripts are initiated from transposable elements. However, the promoters driving such transcriptional activity are mostly unknown. We experimentally characterized an antisense (AS) promoter in mouse L1 retrotransposons for the first time, oriented antiparallel to the coding strand of L1 open reading frame-1. We found that AS transcription is mediated by RNA polymerase II. Rapid amplification of cDNA ends cloning mapped transcription start sites adjacent to the AS promoter. We identified >100 novel fusion transcripts, of which many were conserved across divergent mouse lineages, suggesting conservation of potential functions. To evaluate whether AS L1 transcription could regulate L1 retrotransposition, we replaced portions of native open reading frame-1 in donor elements by synonymously recoded sequences. The resulting L1 elements lacked AS promoter activity and retrotransposed more frequently than endogenous L1s. Overexpression of AS L1 transcripts also reduced L1 retrotransposition. This suppression of retrotransposition was largely independent of Dicer. Our experiments shed new light on how AS fusion transcripts are initiated from endogenous L1 elements across the mouse genome. Such AS transcription can contribute substantially both to natural transcriptional variation and to endogenous regulation of L1 retrotransposition.

The methylation of DNA repeat elements is sex-dependent and temporally different in response to X radiation in radiosensitive and radioresistant mouse strains

The effects of ionizing radiation on DNA methylation are of importance due to the role that DNA methylation plays in maintaining genome stability, and the presence of aberrant DNA methylation in many cancers. There is limited evidence that radiation-sensitivity may influence the modulation of DNA methylation by ionizing radiation, resulting in a loss of methylation. The BALB/c, CBA and C57Bl/6 strains are the most commonly utilized mouse strains in radiation research and are classified as radiation sensitive (BALB/c and CBA) or radiation resistant (C57Bl/6). We present here the first direct comparison of changes in repeat element DNA methylation (L1, B1 and Intracisternal A Particle; IAP) over time in these three mouse strains after high-dose radiation exposure. Using a high-resolution melt assay, methylation of the spleen repeat elements was investigated between 1 and 14 days after whole-body irradiation with 1 Gy X rays. Our study demonstrated that rather than a loss of methylation at the elements, all strains exhibited an early increase in L1 methylation one day after irradiation. In the most radiosensitive strain (BALB/c) the increase was also detected at 6 days postirradiation. The radioresistant C57Bl/6 strain exhibited a loss of L1 methylation at 14 days postirradiation. Less extensive changes to the B1 and IAP elements were detected at various time points, and pyrosequencing revealed that the responses of the strains were influenced by sex, with the male BALB/c and CBA mice exhibiting a greater response to the irradiation. The results of our study do not support the hypothesis that the most radiosensitive strains exhibit the greatest loss of repeat element DNA methylation after exposure to high-dose radiation. While the exact mechanism and biological outcome of the changes in DNA methylation observed here are still to be elucidated, this study provides the first evidence that radiation exposure elicits time-dependent changes in the methylation of repeat elements that are influenced by the genetic background, gender and the type of repeat element investigated. Furthermore, it suggest that any induced changes may not be persistent.

Non-LTR retrotransposons and microsatellites: Partners in genomic variation

The human genome is laden with both non-LTR (long-terminal repeat) retrotransposons and microsatellite repeats. Both types of sequences are able to, either actively or passively, mutagenize the genomes of human individuals and are therefore poised to dynamically alter the human genomic landscape across generations. Non-LTR retrotransposons, such as L1 and Alu, are a major source of new microsatellites, which are born both concurrently and subsequently to L1 and Alu integration into the genome. Likewise, the mutation dynamics of microsatellite repeats have a direct impact on the fitness of their non-LTR retrotransposon parent owing to microsatellite expansion and contraction. This review explores the interactions and dynamics between non-LTR retrotransposons and microsatellites in the context of genomic variation and evolution.

The specificity and flexibility of l1 reverse transcription priming at imperfect T-tracts

L1 retrotransposons have a prominent role in reshaping mammalian genomes. To replicate, the L1 ribonucleoprotein particle (RNP) first uses its endonuclease (EN) to nick the genomic DNA. The newly generated DNA end is subsequently used as a primer to initiate reverse transcription within the L1 RNA poly(A) tail, a process known as target-primed reverse transcription (TPRT). Prior studies demonstrated that most L1 insertions occur into sequences related to the L1 EN consensus sequence (degenerate 5′-TTTT/A-3′ sites) and frequently preceded by imperfect T-tracts. However, it is currently unclear whether—and to which degree—the liberated 3′-hydroxyl extremity on the genomic DNA needs to be accessible and complementary to the poly(A) tail of the L1 RNA for efficient priming of reverse transcription. Here, we employed a direct assay for the initiation of L1 reverse transcription to define the molecular rules that guide this process. First, efficient priming is detected with as few as 4 matching nucleotides at the primer 3′ end. Second, L1 RNP can tolerate terminal mismatches if they are compensated within the 10 last bases of the primer by an increased number of matching nucleotides. All terminal mismatches are not equally detrimental to DNA extension, a C being extended at higher levels than an A or a G. Third, efficient priming in the context of duplex DNA requires a 3′ overhang. This suggests the possible existence of additional DNA processing steps, which generate a single-stranded 3′ end to allow L1 reverse transcription. Based on these data we propose that the specificity of L1 reverse transcription initiation contributes, together with the specificity of the initial EN cleavage, to the distribution of new L1 insertions within the human genome.

Jumping genes are DNA sequences present in the genome of most living organisms. They contribute to genome dynamics and occasionally result in hereditary genetic diseases or cancer. L1 elements are the only autonomously active jumping genes in the human genome. They replicate through an RNA–mediated copy-and-paste mechanism by cleaving the host genome and then using this new DNA end as a primer to reverse transcribe its own RNA, generating a new L1 DNA copy. The molecular determinants that influence L1 target site choice are not fully understood. Here we present a quantitative assay to measure the influence of DNA target site sequence and structure on the reverse transcription step. By testing more than 65 potential DNA primers, we observe that not all sites are equally extended by the L1 machinery, and we define the rules guiding this process. In particular, we highlight the importance of partial sequence complementarity between the target site and the L1 RNA extremity, but also the high level of flexibility of this process, since detrimental terminal mismatches can be compensated by an increasing number of interacting nucleotides. We propose that this mechanism contributes to the distribution of new L1 insertions within the human genome.

Genome-Wide Detection of Gene Coexpression Domains Showing Linkage to Regions Enriched with Polymorphic Retrotransposons in Recombinant Inbred Mouse Strains

Although gene coexpression domains have been reported in most eukaryotic organisms, data available to date suggest that coexpression rarely concerns more than doublets or triplets of adjacent genes in mammals. Using expression data from hearts of mice from the panel of AxB/BxA recombinant inbred mice, we detected (according to window sizes) 42−53 loci linked to the expression levels of clusters of three or more neighboring genes. These loci thus formed “cis-expression quantitative trait loci (eQTL) clusters” because their position matched that of the genes whose expression was linked to the loci. Compared with matching control regions, genes contained within cis-eQTL clusters showed much greater levels of coexpression. Corresponding regions showed: (1) a greater abundance of polymorphic elements (mostly short interspersed element retrotransposons), and (2) significant enrichment for the motifs of binding sites for various transcription factors, with binding sites for the chromatin-organizing CCCTC-binding factor showing the greatest levels of enrichment in polymorphic short interspersed elements. Similar cis-eQTL clusters also were detected when we used data obtained with several tissues from BxD recombinant inbred mice. In addition to strengthening the evidence for gene expression domains in mammalian genomes, our data suggest a possible mechanism whereby noncoding polymorphisms could affect the coordinate expression of several neighboring genes.

How do mammalian transposons induce genetic variation? A conceptual framework: the age, structure, allele frequency, and genome context of transposable elements may define their wide-ranging biological impacts

In this essay, we discuss new insights into the wide-ranging impacts of mammalian transposable elements (TE) on gene expression and function. Nearly half of each mammalian genome is comprised of these mobile, repetitive elements. While most TEs are ancient relics, certain classes can move from one chromosomal location to another even now. Indeed, striking recent data show that extensive transposition occurs not only in the germline over evolutionary time, but also in developing somatic tissues and particular human cancers. While occasional germline TE insertions may contribute to genetic variation, many other, similar TEs appear to have little or no impact on neighboring genes. However, the effects of somatic insertions on gene expression and function remain almost completely unknown. We present a conceptual framework to understand how the ages, allele frequencies, molecular structures, and especially the genomic context of mammalian TEs each can influence their various possible functional consequences.

Revisiting the evolution of mouse LINE-1 in the genomic era

Background

LINE-1 (L1) is the dominant category of transposable elements in placental mammals. L1 has significantly affected the size and structure of all mammalian genomes and understanding the nature of the interactions between L1 and its mammalian host remains a question of crucial importance in comparative genomics. For this reason, much attention has been dedicated to the evolution of L1. Among the most studied elements is the mouse L1 which has been the subject of a number of studies in the 1980s and 1990s. These seminal studies, performed in the pre-genomic era when only a limited number of L1 sequences were available, have significantly improved our understanding of L1 evolution. Yet, no comprehensive study on the evolution of L1 in mouse has been performed since the completion of this genome sequence.

Results

Using the Genome Parsing Suite we performed the first evolutionary analysis of mouse L1 over the entire length of the element. This analysis indicates that the mouse L1 has recruited novel 5’UTR sequences more frequently than previously thought and that the simultaneous activity of non-homologous promoters seems to be one of the conditions for the co-existence of multiple L1 families or lineages. In addition the exchange of genetic information between L1 families is not limited to the 5’UTR as evidence of inter-family recombination was observed in ORF1, ORF2, and the 3’UTR. In contrast to the human L1, there was little evidence of rapid amino-acid replacement in the coiled-coil of ORF1, although this region is structurally unstable. We propose that the structural instability of the coiled-coil domain might be adaptive and that structural changes in this region are selectively equivalent to the rapid evolution at the amino-acid level reported in the human lineage.

Conclusions

The pattern of evolution of L1 in mouse shows some similarity with human suggesting that the nature of the interactions between L1 and its host might be similar in these two species. Yet, some notable differences, particularly in the evolution of ORF1, suggest that the molecular mechanisms involved in host-L1 interactions might be different in these two species.

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.

Cocaine produces conditioned place aversion in mice with a cocaine-insensitive dopamine transporter

Cocaine is an inhibitor of the dopamine, norepinephrine and serotonin reuptake transporters. Because its administration would elevate signaling of all these three neurotransmitters, many studies have been aimed at attributing individual effects of cocaine to specific transmitter systems. Using mice with a cocaine-insensitive dopamine transporter (DAT-CI mice), we previously showed that cocaine-induced dopamine elevations were necessary for its rewarding and stimulating effects. In this study, we observe that DAT-CI mice exhibit cocaine-conditioned place aversion (CPA), and that its expression depends on their genetic background. Specifically, DAT-CI mice backcrossed to the C57Bl/6J strain background did not display a preference or an aversion to cocaine, whereas DAT-CI mice that were on a mixed 129S1/SvImJ × C57Bl/6J (129B6) background had a robust CPA to cocaine. These results indicate that while inhibition of the DAT is necessary for cocaine reward, other cocaine targets and neurotransmitter systems may mediate the aversive properties of cocaine. Furthermore, the aversive effect of cocaine can be observed in the absence of a DAT-mediated rewarding effect, and it is affected by genomic differences between these two mouse strains.

A streamlined method for detecting structural variants in cancer genomes by short read paired-end sequencing

Defining the architecture of a specific cancer genome, including its structural variants, is essential for understanding tumor biology, mechanisms of oncogenesis, and for designing effective personalized therapies. Short read paired-end sequencing is currently the most sensitive method for detecting somatic mutations that arise during tumor development. However, mapping structural variants using this method leads to a large number of false positive calls, mostly due to the repetitive nature of the genome and the difficulty of assigning correct mapping positions to short reads. This study describes a method to efficiently identify large tumor-specific deletions, inversions, duplications and translocations from low coverage data using SVDetect or BreakDancer software and a set of novel filtering procedures designed to reduce false positive calls. Applying our method to a spontaneous T cell lymphoma arising in a core RAG2/p53-deficient mouse, we identified 40 validated tumor-specific structural rearrangements supported by as few as 2 independent read pairs.

Transposable element small RNAs as regulators of gene expression

Transposable elements (TEs) are a source of endogenous small RNAs in animals and plants. These TE-derived small RNAs have been traditionally treated as functionally distinct from gene-regulating small RNAs, such as miRNAs. Two recent reports in Drosophila and Arabidopsis have blurred the lines of this distinction. In both examples, epigenetically and developmentally regulated bursts in TE expression produce gene-regulating small RNAs. In the Drosophila early embryo, maternally deposited TE-derived PIWI-interacting small RNAs (piRNAs) play a role in regulating the nanos mRNA through small RNA binding sites in the nanos 3' untranslated region (UTR). In Arabidopsis, when Athila retrotransposons are epigenetically activated, their transcripts are processed into small RNAs, which directly target the 3'UTR of the genic oligouridylate binding protein 1B (UBP1b) mRNA. Based on these two examples, we suggest that other TE-derived small RNAs regulate additional genes and propose that, through small RNAs, the epigenetic status of TEs could widely influence the genic transcriptome.

Circular retrotransposition products generated by a LINE retrotransposon

Non-long terminal repeat (non-LTR) retrotransposons are highly abundant elements that are present in chromosomes throughout the eukaryotic domain of life. The long interspersed nuclear element (LINE-1) (L1) clade of non-LTR retrotransposons has been particularly successful in mammals, accounting for 30–40% of human genome sequence. The current model of LINE retrotransposition, target-primed reverse transcription, culminates in a chromosomally integrated end product. Using a budding yeast model of non-LTR retrotransposition, we show that in addition to producing these ‘classical’, chromosomally integrated products, a fungal L1 clade member (Zorro3) can generate abundant, RNA-derived episomal products. Genetic evidence suggests that these products are likely to be formed via a variation of target-primed reverse transcription. These episomal products are a previously unseen alternative fate of LINE retrotransposition, and may represent an unexpected source for de novo retrotransposition.

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.

Next-generation sequencing of experimental mouse strains

Since the turn of the century the complete genome sequence of just one mouse strain, C57BL/6J, has been available. Knowing the sequence of this strain has enabled large-scale forward genetic screens to be performed, the creation of an almost complete set of embryonic stem (ES) cell lines with targeted alleles for protein-coding genes, and the generation of a rich catalog of mouse genomic variation. However, many experiments that use other common laboratory mouse strains have been hindered by a lack of whole-genome sequence data for these strains. The last 5 years has witnessed a revolution in DNA sequencing technologies. Recently, these technologies have been used to expand the repertoire of fully sequenced mouse genomes. In this article we review the main findings of these studies and discuss how the sequence of mouse genomes is helping pave the way from sequence to phenotype. Finally, we discuss the prospects for using de novo assembly techniques to obtain high-quality assembled genome sequences of these laboratory mouse strains, and what advances in sequencing technologies may be required to achieve this goal.

The bonobo genome compared with the chimpanzee and human genomes

Two African apes are the closest living relatives of humans: the chimpanzee (Pan troglodytes) and the bonobo (Pan paniscus). Although they are similar in many respects, bonobos and chimpanzees differ strikingly in key social and sexual behaviours, and for some of these traits they show more similarity with humans than with each other. Here we report the sequencing and assembly of the bonobo genome to study its evolutionary relationship with the chimpanzee and human genomes. We find that more than three per cent of the human genome is more closely related to either the bonobo or the chimpanzee genome than these are to each other. These regions allow various aspects of the ancestry of the two ape species to be reconstructed. In addition, many of the regions that overlap genes may eventually help us understand the genetic basis of phenotypes that humans share with one of the two apes to the exclusion of the other.

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.

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.

Mouse endogenous retroviruses can trigger premature transcriptional termination at a distance

Endogenous retrotransposons have caused extensive genomic variation within mammalian species, but the functional implications of such mobilization are mostly unknown. We mapped thousands of endogenous retrovirus (ERV) germline integrants in highly divergent, previously unsequenced mouse lineages, facilitating a comparison of gene expression in the presence or absence of local insertions. Polymorphic ERVs occur relatively infrequently in gene introns and are particularly depleted from genes involved in embryogenesis or that are highly expressed in embryonic stem cells. Their genomic distribution implies ongoing negative selection due to deleterious effects on gene expression and function. A polymorphic, intronic ERV at Slc15a2 triggers up to 49-fold increases in premature transcriptional termination and up to 39-fold reductions in full-length transcripts in adult mouse tissues, thereby disrupting protein expression and functional activity. Prematurely truncated transcripts also occur at Polr1a, Spon1, and up to ∼5% of other genes when intronic ERV polymorphisms are present. Analysis of expression quantitative trait loci (eQTLs) in recombinant BxD mouse strains demonstrated very strong genetic associations between the polymorphic ERV in cis and disrupted transcript levels. Premature polyadenylation is triggered at genomic distances up to >12.5 kb upstream of the ERV, both in cis and between alleles. The parent of origin of the ERV is associated with variable expression of nonterminated transcripts and differential DNA methylation at its 5'-long terminal repeat. This study defines an unexpectedly strong functional impact of ERVs in disrupting gene transcription at a distance and demonstrates that ongoing retrotransposition can contribute significantly to natural phenotypic diversity.

Accumulation and rapid decay of non-LTR retrotransposons in the genome of the three-spine stickleback

The diversity and abundance of non–long terminal repeat (LTR) retrotransposons (nLTR-RT) differ drastically among vertebrate genomes. At one extreme, the genome of placental mammals is littered with hundreds of thousands of copies resulting from the activity of a single clade of nLTR-RT, the L1 clade. In contrast, fish genomes contain a much more diverse repertoire of nLTR-RT, represented by numerous active clades and families. Yet, the number of nLTR-RT copies in teleostean fish is two orders of magnitude smaller than in mammals. The vast majority of insertions appear to be very recent, suggesting that nLTR-RT do not accumulate in fish genomes. This pattern had previously been explained by a high rate of turnover, in which the insertion of new elements is offset by the selective loss of deleterious inserts. The turnover model was proposed because of the similarity between fish and Drosophila genomes with regard to their nLTR-RT profile. However, it is unclear if this model applies to fish. In fact, a previous study performed on the puffer fish suggested that transposable element insertions behave as neutral alleles. Here we examined the dynamics of amplification of nLTR-RT in the three-spine stickleback (Gasterosteus aculeatus). In this species, the vast majority of nLTR-RT insertions are relatively young, as suggested by their low level of divergence. Contrary to expectations, a majority of these insertions are fixed in lake and oceanic populations; thus, nLTR-RT do indeed accumulate in the genome of their fish host. This is not to say that nLTR-RTs are fully neutral, as the lack of fixed long elements in this genome suggests a deleterious effect related to their length. This analysis does not support the turnover model and strongly suggests that a much higher rate of DNA loss in fish than in mammals is responsible for the relatively small number of nLTR-RT copies and for the scarcity of ancient elements in fish genomes. We further demonstrate that nLTR-RT decay in fish occurs mostly through large deletions and not by the accumulation of small deletions.

Genome sequencing of mouse induced pluripotent stem cells reveals retroelement stability and infrequent DNA rearrangement during reprogramming

The biomedical utility of induced pluripotent stem cells (iPSCs) will be diminished if most iPSC lines harbor deleterious genetic mutations. Recent microarray studies have shown that human iPSCs carry elevated levels of DNA copy number variation compared with those in embryonic stem cells, suggesting that these and other classes of genomic structural variation (SV), including inversions, smaller duplications and deletions, complex rearrangements, and retroelement transpositions, may frequently arise as a consequence of reprogramming. Here we employ whole-genome paired-end DNA sequencing and sensitive mapping algorithms to identify all classes of SV in three fully pluripotent mouse iPSC lines. Despite the improved scope and resolution of this study, we find few spontaneous mutations per line (one or two) and no evidence for endogenous retroelement transposition. These results show that genome stability can persist throughout reprogramming, and argue that it is possible to generate iPSCs lacking gene-disrupting mutations using current reprogramming methods.