Whole-genome resequencing allows detection of many rare LINE-1 insertion alleles in humans

Evolutionary insights from profiling LINE-1 activity at allelic resolution in a single human genome

Transposable elements have created the majority of the sequence in many genomes. In mammals, LINE-1 retrotransposons have been expanding for more than 100 million years as distinct, consecutive lineages; however, the drivers of this recurrent lineage emergence and disappearance are unknown. Most human genome assemblies provide a record of this ancient evolution, but fail to resolve ongoing LINE-1 retrotranspositions. Utilizing the human CHM1 long-read-based haploid assembly, we identified and cloned all full-length, intact LINE-1s, and found 29 LINE-1s with measurable in vitro retrotransposition activity. Among individuals, these LINE-1s varied in their presence, their allelic sequences, and their activity. We found that recently retrotransposed LINE-1s tend to be active in vitro and polymorphic in the population relative to more ancient LINE-1s. However, some rare allelic forms of old LINE-1s retain activity, suggesting older lineages can persist longer than expected. Finally, in LINE-1s with in vitro activity and in vivo fitness, we identified mutations that may have increased replication in ancient genomes and may prove promising candidates for mechanistic investigations of the drivers of LINE-1 evolution and which LINE-1 sequences contribute to human disease.

Translational Significance of the LINE-1 Jumping Gene in Skeletal Muscle

Retrotransposons are gene segments that proliferate in the genome, and the Long INterspersed Element 1 (LINE-1 or L1) retrotransposon is active in humans. Although older mammals show enhanced skeletal muscle L1 expression, exercise generally reverses this trend. We hypothesize skeletal muscle L1 expression influences muscle physiology, and additional innovative investigations are needed to confirm this hypothesis.

A comprehensive analysis of gorilla-specific LINE-1 retrotransposons

BACKGROUND

Long interspersed element-1 (LINE-1 or L1) is the most abundant retrotransposons in the primate genome. They have approximately 520,000 copies and make up ~ 17% of the primate genome. Full-length L1s can mobilize to a new genomic location using their enzymatic machinery. Gorilla is the second closest species to humans after the chimpanzee, and human-gorilla split 7-12 million years ago. The gorilla genome provides an opportunity to explore primate origins and evolution.

OBJECTIVE

L1s have contributed to genome diversity and variations during primate evolution. This study aimed to identify gorilla-specific L1s using a more recent version of the gorilla reference genome (Mar. 2016 GSMRT3/gorGor5).

METHODS

We collected gorilla-specific L1 candidates through computational analysis and manual inspection. L1Xplorer was used to identify whether full-length gorilla-specific L1s were intact. In addition, to determine the level of sequence conservation between intact fulllength gorilla-specific L1s, two ORFs of intact L1s were aligned with the L1PA2 consensus sequence.

RESULTS

2002 gorilla-specific L1 candidates were identified through computational analysis. Among them, we manually inspected 1,883 gorilla-specific L1s, among which most of them belong to the L1PA2 subfamily and 12 were intact L1s that could influence genomic variations in the gorilla genome. Interestingly, the 12 intact full-length gorilla-specific L1s have 14 highly conserved nonsynonymous mutations, including 6 mutations and 8 mutations in ORF1 and ORF2, respectively. In comparison to the intact full-length chimpanzee-specific L1s and human-specific hot-L1s, two of these in ORF1 (L256F and E293G) were shown as gorilla-specific nonsynonymous mutations.

CONCLUSION

The gorilla-specific L1s may have had significantly affected the gorilla genome to compose a genome different form that of other primates during primate evolution.

Human L1 Transposition Dynamics Unraveled with Functional Data Analysis

Long INterspersed Elements-1 (L1s) constitute >17% of the human genome and still actively transpose in it. Characterizing L1 transposition across the genome is critical for understanding genome evolution and somatic mutations. However, to date, L1 insertion and fixation patterns have not been studied comprehensively. To fill this gap, we investigated three genome-wide data sets of L1s that integrated at different evolutionary times: 17,037 de novo L1s (from an L1 insertion cell-line experiment conducted in-house), and 1,212 polymorphic and 1,205 human-specific L1s (from public databases). We characterized 49 genomic features-proxying chromatin accessibility, transcriptional activity, replication, recombination, etc.-in the ±50 kb flanks of these elements. These features were contrasted between the three L1 data sets and L1-free regions using state-of-the-art Functional Data Analysis statistical methods, which treat high-resolution data as mathematical functions. Our results indicate that de novo, polymorphic, and human-specific L1s are surrounded by different genomic features acting at specific locations and scales. This led to an integrative model of L1 transposition, according to which L1s preferentially integrate into open-chromatin regions enriched in non-B DNA motifs, whereas they are fixed in regions largely free of purifying selection-depleted of genes and noncoding most conserved elements. Intriguingly, our results suggest that L1 insertions modify local genomic landscape by extending CpG methylation and increasing mononucleotide microsatellite density. Altogether, our findings substantially facilitate understanding of L1 integration and fixation preferences, pave the way for uncovering their role in aging and cancer, and inform their use as mutagenesis tools in genetic studies.

Inherited L1 Retrotransposon Insertions Associated With Risk for Schizophrenia and Bipolar Disorder

Studies of the genetic heritability of schizophrenia and bipolar disorder examining single nucleotide polymorphisms (SNPs) and copy number variations have failed to explain a large portion of the genetic liability, resulting in substantial missing heritability. Long interspersed element 1 (L1) retrotransposons are a type of inherited polymorphic variant that may be associated with risk for schizophrenia and bipolar disorder. We performed REBELseq, a genome wide assay for L1 sequences, on DNA from male and female persons with schizophrenia and controls (n = 63 each) to identify inherited L1 insertions and validated priority insertions. L1 insertions of interest were genotyped in DNA from a replication cohort of persons with schizophrenia, bipolar disorder, and controls (n = 2268 each) to examine differences in carrier frequencies. We identified an inherited L1 insertion in ARHGAP24 and a quadallelic SNP (rs74169643) inside an L1 insertion in SNTG2 that are associated with risk for developing schizophrenia and bipolar disorder (all odds ratios ~1.2). Pathway analysis identified 15 gene ontologies that were differentially affected by L1 burden, including multiple ontologies related to glutamatergic signaling and immune function, which have been previously associated with schizophrenia. These findings provide further evidence supporting the role of inherited repetitive genetic elements in the heritability of psychiatric disorders.

Nanopore Sequencing Enables Comprehensive Transposable Element Epigenomic Profiling

Transposable elements (TEs) drive genome evolution and are a notable source of pathogenesis, including cancer. While CpG methylation regulates TE activity, the locus-specific methylation landscape of mobile human TEs has to date proven largely inaccessible. Here, we apply new computational tools and long-read nanopore sequencing to directly infer CpG methylation of novel and extant TE insertions in hippocampus, heart, and liver, as well as paired tumor and non-tumor liver. As opposed to an indiscriminate stochastic process, we find pronounced demethylation of young long interspersed element 1 (LINE-1) retrotransposons in cancer, often distinct to the adjacent genome and other TEs. SINE-VNTR-Alu (SVA) retrotransposons, including their internal tandem repeat-associated CpG island, are near-universally methylated. We encounter allele-specific TE methylation and demethylation of aberrantly expressed young LINE-1s in normal tissues. Finally, we recover the complete sequences of tumor-specific LINE-1 insertions and their retrotransposition hallmarks, demonstrating how long-read sequencing can simultaneously survey the epigenome and detect somatic TE mobilization.

Non-full-length Water-Soluble CXCR4QTY and CCR5QTY Chemokine Receptors: Implication for Overlooked Truncated but Functional Membrane Receptors

It was posited that functionalities of GPCRs require full-length sequences that are negated by residue deletions. Here we report that significantly truncated nfCCR5^QTY and nfCXCR4^QTY still bind native ligands. Receptor-ligand interactions were discovered from yeast 2-hybrid screening and confirmed by mating selection. Two nfCCR5^QTY (SZ218a, SZ190b) and two nfCXCR4^QTY (SZ158a, SZ146a) were expressed in E. coli. Synthesized receptors exhibited α-helical structures and bound respective ligands with reduced affinities. SZ190b and SZ158a were reconverted into non-QTY forms and expressed in HEK293T cells. Reconverted receptors localized on cell membranes and functioned as negative regulators for ligand-induced signaling when co-expressed with full-length receptors. CCR5-SZ190b individually can perform signaling at a reduced level with higher ligand concentration. Our findings provide insight into essential structural components for CCR5 and CXCR4 functionality, while raising the possibility that non-full-length receptors may be resulted from alternative splicing and that pseudo-genes in genomes may be present and functional in living organisms.

•

Y2H screening reveals ligand interaction from truncated CXCR4 and CCR5 in QTY form

•

Truncated CCR5^QTY and CXCR4^QTY can be produced in E. coli and bind native ligands

•

Reconverted receptors localize on membranes and regulate cell signaling in HEK293

•

Our finding indicates potential presence and function for truncated receptors

Restriction Enzyme Based Enriched L1Hs Sequencing (REBELseq): A Scalable Technique for Detection of Ta Subfamily L1Hs in the Human Genome

Long interspersed element-1 retrotransposons (LINE-1 or L1) are ∼6 kb mobile DNA elements implicated in the origins of many Mendelian and complex diseases. The actively retrotransposing L1s are mostly limited to the L1 human specific (L1Hs) transcriptional active (Ta) subfamily. In this manuscript, we present REBELseq as a method for the construction of Ta subfamily L1Hs-enriched next-generation sequencing libraries and bioinformatic identification. REBELseq was performed on DNA isolated from NeuN+ neuronal nuclei from postmortem brain samples of 177 individuals and empirically-driven bioinformatic and experimental cutoffs were established. Putative L1Hs insertions passing bioinformatics cutoffs were experimentally validated. REBELseq reliably identified both known and novel Ta subfamily L1Hs insertions distributed throughout the genome. Differences in the proportion of individuals possessing a given reference or non-reference retrotransposon insertion were identified. We conclude that REBELseq is an unbiased, whole genome approach to the amplification and detection of Ta subfamily L1Hs retrotransposons.

Human transposon insertion profiling by sequencing (TIPseq) to map LINE-1 insertions in single cells

Long interspersed element-1 (LINE-1, L1) sequences, which comprise about 17% of human genome, are the product of one of the most active types of mobile DNAs in modern humans. LINE-1 insertion alleles can cause inherited and de novo genetic diseases, and LINE-1-encoded proteins are highly expressed in some cancers. Genome-wide LINE-1 mapping in single cells could be useful for defining somatic and germline retrotransposition rates, and for enabling studies to characterize tumour heterogeneity, relate insertions to transcriptional and epigenetic effects at the cellular level, or describe cellular phylogenies in development. Our laboratories have reported a genome-wide LINE-1 insertion site mapping method for bulk DNA, named transposon insertion profiling by sequencing (TIPseq). There have been significant barriers applying LINE-1 mapping to single cells, owing to the chimeric artefacts and features of repetitive sequences. Here, we optimize a modified TIPseq protocol and show its utility for LINE-1 mapping in single lymphoblastoid cells. Results from single-cell TIPseq experiments compare well to known LINE-1 insertions found by whole-genome sequencing and TIPseq on bulk DNA. Among the several approaches we tested, whole-genome amplification by multiple displacement amplification followed by restriction enzyme digestion, vectorette ligation and LINE-1-targeted PCR had the best assay performance. This article is part of a discussion meeting issue 'Crossroads between transposons and gene regulation'.

A potential new mechanism for pregnancy loss: considering the role of LINE-1 retrotransposons in early spontaneous miscarriage

LINE1 retrotransposons are mobile DNA elements that copy and paste themselves into new sites in the genome. To ensure their evolutionary success, heritable new LINE-1 insertions accumulate in cells that can transmit genetic information to the next generation (i.e., germ cells and embryonic stem cells). It is our hypothesis that LINE1 retrotransposons, insertional mutagens that affect expression of genes, may be causal agents of early miscarriage in humans. The cell has evolved various defenses restricting retrotransposition-caused mutation, but these are occasionally relaxed in certain somatic cell types, including those of the early embryo. We predict that reduced suppression of L1s in germ cells or early-stage embryos may lead to excessive genome mutation by retrotransposon insertion, or to the induction of an inflammatory response or apoptosis due to increased expression of L1-derived nucleic acids and proteins, and so disrupt gene function important for embryogenesis. If correct, a novel threat to normal human development is revealed, and reverse transcriptase therapy could be one future strategy for controlling this cause of embryonic damage in patients with recurrent miscarriages.

L1 Retrotransposon Heterogeneity in Ovarian Tumor Cell Evolution

LINE-1 (L1) retrotransposons are a source of insertional mutagenesis in tumor cells. However, the clinical significance of L1 mobilization during tumorigenesis remains unclear. Here, we applied retrotransposon capture sequencing (RC-seq) to multiple single-cell clones isolated from five ovarian cancer cell lines and HeLa cells and detected endogenous L1 retrotransposition in vitro. We then applied RC-seq to ovarian tumor and matched blood samples from 19 patients and identified 88 tumor-specific L1 insertions. In one tumor, an intronic de novo L1 insertion supplied a novel cis-enhancer to the putative chemoresistance gene STC1. Notably, the tumor subclone carrying the STC1 L1 mutation increased in prevalence after chemotherapy, further increasing STC1 expression. We also identified hypomethylated donor L1s responsible for new L1 insertions in tumors and cultivated cancer cells. These congruent in vitro and in vivo results highlight L1 insertional mutagenesis as a common component of ovarian tumorigenesis and cancer genome heterogeneity.

Comparative Genomics Analysis Reveals High Levels of Differential Retrotransposition among Primates from the Hominidae and the Cercopithecidae Families

Mobile elements (MEs), making ∼50% of primate genomes, are known to be responsible for generating inter- and intra-species genomic variations and play important roles in genome evolution and gene function. Using a bioinformatics comparative genomics approach, we performed analyses of species-specific MEs (SS-MEs) in eight primate genomes from the families of Hominidae and Cercopithecidae, focusing on retrotransposons. We identified a total of 230,855 SS-MEs, with which we performed normalization based on evolutionary distances, and we also analyzed the most recent SS-MEs in these genomes. Comparative analysis of SS-MEs reveals striking differences in ME transposition among these primate genomes. Interesting highlights of our results include: 1) the baboon genome has the highest number of SS-MEs with a strong bias for SINEs, while the crab-eating macaque genome has a sustained extremely low transposition for all ME classes, suggesting the existence of a genome-wide mechanism suppressing ME transposition; 2) while SS-SINEs represent the dominant class in general, the orangutan genome stands out by having SS-LINEs as the dominant class; 3) the human genome stands out among the eight genomes by having the largest number of recent highly active ME subfamilies, suggesting a greater impact of ME transposition on its recent evolution; and 4) at least 33% of the SS-MEs locate to genic regions, including protein coding regions, presenting significant potentials for impacting gene function. Our study, as the first of its kind, demonstrates that mobile elements evolve quite differently among these primates, suggesting differential ME transposition as an important mechanism in primate evolution.

SQuIRE reveals locus-specific regulation of interspersed repeat expression

Transposable elements (TEs) are interspersed repeat sequences that make up much of the human genome. Their expression has been implicated in development and disease. However, TE-derived RNA-seq reads are difficult to quantify. Past approaches have excluded these reads or aggregated RNA expression to subfamilies shared by similar TE copies, sacrificing quantitative accuracy or the genomic context necessary to understand the basis of TE transcription. As a result, the effects of TEs on gene expression and associated phenotypes are not well understood. Here, we present Software for Quantifying Interspersed Repeat Expression (SQuIRE), the first RNA-seq analysis pipeline that provides a quantitative and locus-specific picture of TE expression (https://github.com/wyang17/SQuIRE). SQuIRE is an accurate and user-friendly tool that can be used for a variety of species. We applied SQuIRE to RNA-seq from normal mouse tissues and a Drosophila model of amyotrophic lateral sclerosis. In both model organisms, we recapitulated previously reported TE subfamily expression levels and revealed locus-specific TE expression. We also identified differences in TE transcription patterns relating to transcript type, gene expression and RNA splicing that would be lost with other approaches using subfamily-level analyses. Altogether, our findings illustrate the importance of studying TE transcription with locus-level resolution.

Transmissible Cancers and Immune Downregulation in Tasmanian Devil (Sacrophilus harrisii) and Canine Populations

Known as devil facial tumor disease (DFTD) and canine transmissible venereal tumor (CTVT), transmissible cancer occurs in both Tasmanian devil and canine populations, respectively. Both malignancies show remarkable ability to be transmitted as allografts into subsequent hosts. How DFTD and CTVT avoid detection by immunocompetent hosts is of particular interest, given that these malignancies are rarely seen in other species in nature. Both of these transmissible cancers can downregulate the host immune system, enabling proliferation. DFTD is characterized by epigenetic modifications to the DNA promoter regions of β₂microglobulin, transporters associated with antigen processing 1 and 2, MHC I, and MHC II-crucial proteins required in the detection and surveillance of foreign material. Downregulation during DFTD may be achieved by altering the activity of histone deacetylases. DFTD has caused widespread destruction of devil populations, placing the species on the brink of extinction. CTVT demonstrates a proliferative phase, during which the tumor evades immune detection, allowing it to proliferate, and a regressive phase when hosts mount an effective immune response. Alteration of TGFβ signaling in CTVT likely impedes the antigen-processing capabilities of canine hosts in addition to hindering the ability of natural killer cells to detect immune system downregulation. Immunosuppressive cytokines such as CXCL7 may contribute to a favorable microenvironment that supports the proliferation of CTVT. When viewed from an evolutionary paradigm, both DFTD and CTVT may conform to a model of host-parasite coevolution. Furthermore, various genetic features, such as genetically active transposons in CTVT and chromosomal rearrangements in DFTD, play important roles in promoting the survival of these disease agents. Understanding the mode of transmission for these transmissible cancers may shed light on mechanisms for human malignancies and reveal opportunities for treatment in the future.

LINE-1 Evasion of Epigenetic Repression in Humans

Epigenetic silencing defends against LINE-1 (L1) retrotransposition in mammalian cells. However, the mechanisms that repress young L1 families and how L1 escapes to cause somatic genome mosaicism in the brain remain unclear. Here we report that a conserved Yin Yang 1 (YY1) transcription factor binding site mediates L1 promoter DNA methylation in pluripotent and differentiated cells. By analyzing 24 hippocampal neurons with three distinct single-cell genomic approaches, we characterized and validated a somatic L1 insertion bearing a 3' transduction. The source (donor) L1 for this insertion was slightly 5' truncated, lacked the YY1 binding site, and was highly mobile when tested in vitro. Locus-specific bisulfite sequencing revealed that the donor L1 and other young L1s with mutated YY1 binding sites were hypomethylated in embryonic stem cells, during neurodifferentiation, and in liver and brain tissue. These results explain how L1 can evade repression and retrotranspose in the human body.

Retrotransposon-induced mosaicism in the neural genome

Over the past decade, major discoveries in retrotransposon biology have depicted the neural genome as a dynamic structure during life. In particular, the retrotransposon LINE-1 (L1) has been shown to be transcribed and mobilized in the brain. Retrotransposition in the developing brain, as well as during adult neurogenesis, provides a milieu in which neural diversity can arise. Dysregulation of retrotransposon activity may also contribute to neurological disease. Here, we review recent reports of retrotransposon activity in the brain, and discuss the temporal nature of retrotransposition and its regulation in neural cells in response to stimuli. We also put forward hypotheses regarding the significance of retrotransposons for brain development and neurological function, and consider the potential implications of this phenomenon for neuropsychiatric and neurodegenerative conditions.

Dynamic Methylation of an L1 Transduction Family during Reprogramming and Neurodifferentiation

The retrotransposon LINE-1 (L1) is a significant source of endogenous mutagenesis in humans. In each individual genome, a few retrotransposition-competent L1s (RC-L1s) can generate new heritable L1 insertions in the early embryo, primordial germ line, and germ cells. L1 retrotransposition can also occur in the neuronal lineage and cause somatic mosaicism. Although DNA methylation mediates L1 promoter repression, the temporal pattern of methylation applied to individual RC-L1s during neurogenesis is unclear. Here, we identified a de novo L1 insertion in a human induced pluripotent stem cell (hiPSC) line via retrotransposon capture sequencing (RC-seq). The L1 insertion was full-length and carried 5' and 3' transductions. The corresponding donor RC-L1 was part of a large and recently active L1 transduction family and was highly mobile in a cultured-cell L1 retrotransposition reporter assay. Notably, we observed distinct and dynamic DNA methylation profiles for the de novo L1 and members of its extended transduction family during neuronal differentiation. These experiments reveal how a de novo L1 insertion in a pluripotent stem cell is rapidly recognized and repressed, albeit incompletely, by the host genome during neurodifferentiation, while retaining potential for further retrotransposition.

Transposon insertion profiling by sequencing (TIPseq) for mapping LINE-1 insertions in the human genome

Background

Transposable elements make up a significant portion of the human genome. Accurately locating these mobile DNAs is vital to understand their role as a source of structural variation and somatic mutation. To this end, laboratories have developed strategies to selectively amplify or otherwise enrich transposable element insertion sites in genomic DNA.

Results

Here we describe a technique, Transposon Insertion Profiling by sequencing (TIPseq), to map Long INterspersed Element 1 (LINE-1, L1) retrotransposon insertions in the human genome. This method uses vectorette PCR to amplify species-specific L1 (L1PA1) insertion sites followed by paired-end Illumina sequencing. In addition to providing a step-by-step molecular biology protocol, we offer users a guide to our pipeline for data analysis, TIPseqHunter. Our recent studies in pancreatic and ovarian cancer demonstrate the ability of TIPseq to identify invariant (fixed), polymorphic (inherited variants), as well as somatically-acquired L1 insertions that distinguish cancer genomes from a patient’s constitutional make-up.

Conclusions

TIPseq provides an approach for amplifying evolutionarily young, active transposable element insertion sites from genomic DNA. Our rationale and variations on this protocol may be useful to those mapping L1 and other mobile elements in complex genomes.

Electronic supplementary material

The online version of this article (10.1186/s13100-019-0148-5) contains supplementary material, which is available to authorized users.

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.

Transposable Elements

Transposable elements (TEs) are low-complexity elements (e.g., LINEs, SINEs, SVAs, and HERVs) that make up to two-thirds of the human genome. There is mounting evidence that TEs play an essential role in molecular functions that influence genomic plasticity and gene expression regulation. With the advent of next-generation sequencing approaches, our understanding of the relationship between TEs and psychiatric disorders will greatly improve. In this chapter, the Authors comprehensively summarize the state-of the-art of TE research in animal models and humans supporting a framework in which TEs play a functional role in mechanisms affecting a variety of behaviors, including neurodevelopmental, neuropsychiatric, and neurodegenerative disorders. Finally, the Authors discuss recent therapeutic applications raised from the increasing experimental evidence on TE functional mechanisms.

Reduced PRC2 function alters male germline epigenetic programming and paternal inheritance

BACKGROUND

Defining the mechanisms that establish and regulate the transmission of epigenetic information from parent to offspring is critical for understanding disease heredity. Currently, the molecular pathways that regulate epigenetic information in the germline and its transmission to offspring are poorly understood.

RESULTS

Here we provide evidence that Polycomb Repressive Complex 2 (PRC2) regulates paternal inheritance. Reduced PRC2 function in mice resulted in male sub-fertility and altered epigenetic and transcriptional control of retrotransposed elements in foetal male germ cells. Males with reduced PRC2 function produced offspring that over-expressed retrotransposed pseudogenes and had altered preimplantation embryo cleavage rates and cell cycle control.

CONCLUSION

This study reveals a novel role for the histone-modifying complex, PRC2, in paternal intergenerational transmission of epigenetic effects on offspring, with important implications for understanding disease inheritance.

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.

Population-wide sampling of retrotransposon insertion polymorphisms using deep sequencing and efficient detection

Active retrotransposons play important roles during evolution and continue to shape our genomes today, especially in genetic polymorphisms underlying a diverse set of diseases. However, studies of human retrotransposon insertion polymorphisms (RIPs) based on whole-genome deep sequencing at the population level have not been sufficiently undertaken, despite the obvious need for a thorough characterization of RIPs in the general population. Herein, we present a novel and efficient computational tool called Specific Insertions Detector (SID) for the detection of non-reference RIPs. We demonstrate that SID is suitable for high-depth whole-genome sequencing data using paired-end reads obtained from simulated and real datasets. We construct a comprehensive RIP database using a large population of 90 Han Chinese individuals with a mean ×68 depth per individual. In total, we identify 9342 recent RIPs, and 8433 of these RIPs are novel compared with dbRIP, including 5826 Alu, 2169 long interspersed nuclear element 1 (L1), 383 SVA, and 55 long terminal repeats. Among the 9342 RIPs, 4828 were located in gene regions and 5 were located in protein-coding regions. We demonstrate that RIPs can, in principle, be an informative resource to perform population evolution and phylogenetic analyses. Taking the demographic effects into account, we identify a weak negative selection on SVA and L1 but an approximately neutral selection for Alu elements based on the frequency spectrum of RIPs. SID is a powerful open-source program for the detection of non-reference RIPs. We built a non-reference RIP dataset that greatly enhanced the diversity of RIPs detected in the general population, and it should be invaluable to researchers interested in many aspects of human evolution, genetics, and disease. As a proof of concept, we demonstrate that the RIPs can be used as biomarkers in a similar way as single nucleotide polymorphisms.

Insights into the RNA binding mechanism of human L1-ORF1p: a molecular dynamics study

The recognition and binding of nucleic acids by ORF1p, an L1 retrotransposon protein, have not yet been clearly understood due to the lack of structural knowledge. The present study attempts to identify the probable single-stranded RNA binding pathway of trimeric ORF1p using computational methods like ligand mapping methodology combined with molecular dynamics simulations. Using the ligand mapping methodology, the possible RNA interacting sites on the surface of the trimeric ORF1p were identified. The crystal structure of the ORF1p timer and an RNA molecule of 29 nucleotide bases in length were used to generate the structure of the ORF1p complex based on information on predicted binding sites as well as the functional states of the CTD. The various complexes of ORF1p-RNA were generated using polyU, polyA and L1RNA sequences and were simulated for a period of 75 ns. The observed stable interaction pattern was used to propose the possible binding pathway. Based on the binding free energy for complex formation, both polyU and L1RNA complexes were identified as stable complexes, while the complex formed with polyA was the least stable one. Furthermore, the importance of the residues in the CC domain (Lys137 and Arg141), the RRM loop (Arg206, Arg210 and Arg211) and the CTD (Arg 261 and Arg262) of all three chains in stabilizing the wrapped RNA has been highlighted in this study. The presence of several electrostatic interactions including H-bond interactions increases the affinity towards RNA and hence plays a vital role in retaining the wrapped position of RNA around ORF1p. Altogether, this study presents one of the possible RNA binding pathways of ORF1p and clearly highlights the functional state of ORF1p visited during RNA binding.

A high throughput screen for active human transposable elements

Background

Transposable elements (TEs) are mobile genetic sequences that randomly propagate within their host’s genome. This mobility has the potential to affect gene transcription and cause disease. However, TEs are technically challenging to identify, which complicates efforts to assess the impact of TE insertions on disease. Here we present a targeted sequencing protocol and computational pipeline to identify polymorphic and novel TE insertions using next-generation sequencing: TE-NGS. The method simultaneously targets the three subfamilies that are responsible for the majority of recent TE activity (L1HS, AluYa5/8, and AluYb8/9) thereby obviating the need for multiple experiments and reducing the amount of input material required.

Results

Here we describe the laboratory protocol and detection algorithm, and a benchmark experiment for the reference genome NA12878. We demonstrate a substantial enrichment for on-target fragments, and high sensitivity and precision to both reference and NA12878-specific insertions. We report 17 previously unreported loci for this individual which are supported by orthogonal long-read evidence, and we identify 1470 polymorphic and novel TEs in 12 additional samples that were previously undocumented in databases of insertion polymorphisms.

Conclusions

We anticipate that future applications of TE-NGS alongside exome sequencing of patients with sporadic disease will reduce the number of unresolved cases, and improve estimates of the contribution of TEs to human genetic disease.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-4485-4) contains supplementary material, which is available to authorized users.

Analysis of LINE-1 Elements in DNA from Postmortem Brains of Individuals with Schizophrenia

Whereas some rare genetic variants convey high risk for schizophrenia (SZ), common alleles conveying even moderate risk remain elusive. Long interspersed element-1s (L1) are mobile retrotransposons comprising ~17% of the human genome. L1 retrotransposition can cause somatic mosaicism during neurodevelopment by insertional mutagenesis. We hypothesized that, compared to controls, patients diagnosed with schizophrenia (PDS) may have increased numbers of deleterious L1 insertions, perhaps occurring de novo, in brain-expressed genes of dorsolateral prefrontal cortex (DLPFC) neurons. Neuronal and non-neuronal nuclei were separated by fluorescence-activated cell sorting from postmortem DLPFC of 36 PDS and 26 age-matched controls. Genomic sequences flanking the 3'-side of L1s were amplified from neuronal DNA, and neuronal L1 libraries were sequenced. Aligned sequences were analyzed for L1 insertions using custom bioinformatics programs. Ontology and pathway analyses were done on lists of genes putatively disrupted by L1s in PDS and controls. Cellular or population allele frequencies of L1s were assessed by droplet digital PCR or Taqman genotyping. We observed a statistically significant increase in the proportion of intragenic novel L1s in DLPFC of PDS. We found over-representation of L1 insertions within the gene ontologies 'cell projection' and 'postsynaptic membrane' in the gene lists derived from PDS samples, but not from controls. Cellular allele frequencies of examined L1 insertions indicated heterozygosity in genomes of DLPFC cells. An L1 within ERI1 exoribonuclease family member 3 (ERI3) was found to associate with SZ. These results extend prior work documenting increased L1 genetic burden in the brains of PDS and also identify unique genes that may provide new insight into the pathophysiology of schizophrenia.

Truncated ORF1 proteins can suppress LINE-1 retrotransposition in trans

Long interspersed element 1 (L1) is an autonomous non-LTR retroelement that is active in mammalian genomes. Although retrotranspositionally incompetent and functional L1 loci are present in the same genomes, it remains unknown whether non-functional L1s have any trans effect on mobilization of active elements. Using bioinformatic analysis, we identified over a thousand of human L1 loci containing at least one stop codon in their ORF1 sequence. RNAseq analysis confirmed that many of these loci are expressed. We demonstrate that introduction of equivalent stop codons in the full-length human L1 sequence leads to the expression of truncated ORF1 proteins. When supplied in trans some truncated human ORF1 proteins suppress human L1 retrotransposition. This effect requires the N-terminus and coiled-coil domain (C-C) as mutations within the ORF1p C-C domain abolish the suppressive effect of truncated proteins on L1 retrotransposition. We demonstrate that the expression levels and length of truncated ORF1 proteins influence their ability to suppress L1 retrotransposition. Taken together these findings suggest that L1 retrotransposition may be influenced by coexpression of defective L1 loci and that these L1 loci may reduce accumulation of de novo L1 integration events.

Genome-wide characterization of non-reference transposable element insertion polymorphisms reveals genetic diversity in tropical and temperate maize

Background

Maize was originally domesticated in a tropical environment but is now widely cultivated at temperate latitudes. Temperate and tropical maize populations have diverged both genotypically and phenotypically. Tropical maize lines grown in temperate environments usually exhibit delayed flowering, pollination, and seed set, which reduces their grain yield relative to temperate adapted maize lines. One potential mechanism by which temperate maize may have adapted to a new environment is novel transposable element insertions, which can influence gene regulation. Recent advances in sequencing technology have made it possible to study variation in transposon content and insertion location in large sets of maize lines.

Results

In total, 274,408 non-redundant TEs (NRTEs) were identified using resequencing data generated from 83 maize inbred lines. The locations of DNA TEs and copia-superfamily retrotransposons showed significant positive correlations with gene density and genetic recombination rates, whereas gypsy-superfamily retrotransposons showed a negative correlation with these two parameters. Compared to tropical maize, temperate maize had fewer unique NRTEs but higher insertion frequency, lower background recombination rates, and higher linkage disequilibrium, with more NRTEs close to flowering and stress-related genes in the genome. Association mapping demonstrated that the presence/absence of 48 NRTEs was associated with flowering time and that expression of neighboring genes differed between haplotypes where a NRTE was present or absent.

Conclusions

This study suggests that NRTEs may have played an important role in creating the variation in gene regulation that enabled the rapid adaptation of maize to diverse environments.

Electronic supplementary material

The online version of this article (10.1186/s12864-017-4103-x) contains supplementary material, which is available to authorized users.

L1 Mosaicism in Mammals: Extent, Effects, and Evolution

The retrotransposon LINE-1 (long interspersed element 1, L1) is a transposable element that has extensively colonized the mammalian germline. L1 retrotransposition can also occur in somatic cells, causing genomic mosaicism, as well as in cancer. However, the extent of L1-driven mosaicism arising during ontogenesis is unclear. We discuss here recent experimental data which, at a minimum, fully substantiate L1 mosaicism in early embryonic development and neural cells, including post-mitotic neurons. We also consider the possible biological impact of somatic L1 insertions in neurons, the existence of donor L1s that are highly active ('hot') in specific spatiotemporal niches, and the evolutionary selection of donor L1s driving neuronal mosaicism.

McClintock: An Integrated Pipeline for Detecting Transposable Element Insertions in Whole-Genome Shotgun Sequencing Data

Transposable element (TE) insertions are among the most challenging types of variants to detect in genomic data because of their repetitive nature and complex mechanisms of replication . Nevertheless, the recent availability of large resequencing data sets has spurred the development of many new methods to detect TE insertions in whole-genome shotgun sequences. Here we report an integrated bioinformatics pipeline for the detection of TE insertions in whole-genome shotgun data, called McClintock (https://github.com/bergmanlab/mcclintock), which automatically runs and standardizes output for multiple TE detection methods. We demonstrate the utility of McClintock by evaluating six TE detection methods using simulated and real genome data from the model microbial eukaryote, Saccharomyces cerevisiae We find substantial variation among McClintock component methods in their ability to detect nonreference TEs in the yeast genome, but show that nonreference TEs at nearly all biologically realistic locations can be detected in simulated data by combining multiple methods that use split-read and read-pair evidence. In general, our results reveal that split-read methods detect fewer nonreference TE insertions than read-pair methods, but generally have much higher positional accuracy. Analysis of a large sample of real yeast genomes reveals that most McClintock component methods can recover known aspects of TE biology in yeast such as the transpositional activity status of families, target preferences, and target site duplication structure, albeit with varying levels of accuracy. Our work provides a general framework for integrating and analyzing results from multiple TE detection methods, as well as useful guidance for researchers studying TEs in yeast resequencing data.

Transposable elements in cancer

Transposable elements give rise to interspersed repeats, sequences that comprise most of our genomes. These mobile DNAs have been historically underappreciated - both because they have been presumed to be unimportant, and because their high copy number and variability pose unique technical challenges. Neither impediment now seems steadfast. Interest in the human mobilome has never been greater, and methods enabling its study are maturing at a fast pace. This Review describes the activity of transposable elements in human cancers, particularly long interspersed element-1 (LINE-1). LINE-1 sequences are self-propagating, protein-coding retrotransposons, and their activity results in somatically acquired insertions in cancer genomes. Altered expression of transposable elements and animation of genomic LINE-1 sequences appear to be hallmarks of cancer, and can be responsible for driving mutations in tumorigenesis.

Transposable elements in cancer

Transposable elements give rise to interspersed repeats, sequences that comprise most of our genomes. These mobile DNAs have been historically underappreciated - both because they have been presumed to be unimportant, and because their high copy number and variability pose unique technical challenges. Neither impediment now seems steadfast. Interest in the human mobilome has never been greater, and methods enabling its study are maturing at a fast pace. This Review describes the activity of transposable elements in human cancers, particularly long interspersed element-1 (LINE-1). LINE-1 sequences are self-propagating, protein-coding retrotransposons, and their activity results in somatically acquired insertions in cancer genomes. Altered expression of transposable elements and animation of genomic LINE-1 sequences appear to be hallmarks of cancer, and can be responsible for driving mutations in tumorigenesis.

Reading LINEs within the cocaine addicted brain

INTRODUCTION

Long interspersed element (LINE)-1 (L1) is a type of retrotransposon capable of mobilizing into new genomic locations. Often studied in Mendelian diseases or cancer, L1s may also cause somatic mutation in the developing central nervous system. Recent reports showed L1 transcription was activated in brains of cocaine-treated mice, and L1 retrotransposition was increased in cocaine-treated neuronal cell cultures. We hypothesized that the predisposition to cocaine addiction may result from inherited L1s or somatic L1 mobilization in the brain.

METHODS

Postmortem medial prefrontal cortex (mPFC) tissue from 30 CA and 30 control individuals was studied. An Alexafluor488-labeled NeuN antibody and fluorescence activated nuclei sorting were used to separate neuronal from non-neuronal cell nuclei. L1s and their 3' flanking sequences were amplified from neuronal and non-neuronal genomic DNA (gDNA) using L1-seq. L1 DNA libraries from the neuronal gDNA were sequenced on an Illumina HiSeq2000. Sequences aligned to the hg19 human genome build were analyzed for L1 insertions using custom "L1-seq" bioinformatics programs.

RESULTS

Previously uncataloged L1 insertions, some validated by PCR, were detected in neurons from both CA and control brain samples. Steady-state L1 mRNA levels in CA and control mPFC were also assessed. Gene ontology and pathway analyses were used to assess relationships between genes putatively disrupted by novel L1s in CA and control individuals. L1 insertions in CA samples were enriched in gene ontologies and pathways previously associated with CA.

CONCLUSIONS

We conclude that neurons in the mPFC harbor L1 insertions that have the potential to influence predisposition to CA.

Human transposon insertion profiling: Analysis, visualization and identification of somatic LINE-1 insertions in ovarian cancer

Mammalian genomes are replete with interspersed repeats reflecting the activity of transposable elements. These mobile DNAs are self-propagating, and their continued transposition is a source of both heritable structural variation as well as somatic mutation in human genomes. Tailored approaches to map these sequences are useful to identify insertion alleles. Here, we describe in detail a strategy to amplify and sequence long interspersed element-1 (LINE-1, L1) retrotransposon insertions selectively in the human genome, transposon insertion profiling by next-generation sequencing (TIPseq). We also report the development of a machine-learning-based computational pipeline, TIPseqHunter, to identify insertion sites with high precision and reliability. We demonstrate the utility of this approach to detect somatic retrotransposition events in high-grade ovarian serous carcinoma.

Study of Transposable Elements and Their Genomic Impact

Transposable elements (TEs) have been considered traditionally as junk DNA, i.e., DNA sequences that despite representing a high proportion of genomes had no evident cellular functions. However, over the last decades, it has become undeniable that not only TE-derived DNA sequences have (and had) a fundamental role during genome evolution, but also TEs have important implications in the origin and evolution of many genomic disorders. This concise review provides a brief overview of the different types of TEs that can be found in genomes, as well as a list of techniques and methods used to study their impact and mobilization. Some of these techniques will be covered in detail in this Method Book.

Evidence for L1-associated DNA rearrangements and negligible L1 retrotransposition in glioblastoma multiforme

Background

LINE-1 (L1) retrotransposons are a notable endogenous source of mutagenesis in mammals. Notably, cancer cells can support unusual L1 retrotransposition and L1-associated sequence rearrangement mechanisms following DNA damage. Recent reports suggest that L1 is mobile in epithelial tumours and neural cells but, paradoxically, not in brain cancers.

Results

Here, using retrotransposon capture sequencing (RC-seq), we surveyed L1 mutations in 14 tumours classified as glioblastoma multiforme (GBM) or as a lower grade glioma. In four GBM tumours, we characterised one probable endonuclease-independent L1 insertion, two L1-associated rearrangements and one likely Alu-Alu recombination event adjacent to an L1. These mutations included PCR validated intronic events in MeCP2 and EGFR. Despite sequencing L1 integration sites at up to 250× depth by RC-seq, we found no tumour-specific, endonuclease-dependent L1 insertions. Whole genome sequencing analysis of the tumours carrying the MeCP2 and EGFR L1 mutations also revealed no endonuclease-dependent L1 insertions. In a complementary in vitro assay, wild-type and endonuclease mutant L1 reporter constructs each mobilised very inefficiently in four cultured GBM cell lines.

Conclusions

These experiments altogether highlight the consistent absence of canonical L1 retrotransposition in GBM tumours and cultured cell lines, as well as atypical L1-associated sequence rearrangements following DNA damage in vivo.

Electronic supplementary material

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

Identifying wrong assemblies in de novo short read primary sequence assembly contigs

With the advent of short-reads-based genome sequencing approaches, large number of organisms are being sequenced all over the world. Most of these assemblies are done using some de novo short read assemblers and other related approaches. However, the contigs produced this way are prone to wrong assembly. So far, there is a conspicuous dearth of reliable tools to identify mis-assembled contigs. Mis-assemblies could result from incorrectly deleted or wrongly arranged genomic sequences. In the present work various factors related to sequence, sequencing and assembling have been assessed for their role in causing mis-assembly by using different genome sequencing data. Finally, some mis-assembly detecting tools have been evaluated for their ability to detect the wrongly assembled primary contigs, suggesting a lot of scope for improvement in this area. The present work also proposes a simple unsupervised learning-based novel approach to identify mis-assemblies in the contigs which was found performing reasonably well when compared to the already existing tools to report mis-assembled contigs. It was observed that the proposed methodology may work as a complementary system to the existing tools to enhance their accuracy.

Genome-wide characterization of non-reference transposons in crops suggests non-random insertion

BACKGROUND

Transposons (transposable elements or TEs) are DNA sequences that can change their position within the genome. A large number of TEs have been identified in reference genome of each crop(named accumulated TEs), which are the important part of genome. However, whether there existed TEs with different insert positions in resequenced crop accession genomes from those of reference genome (named non-reference transposable elements, non-ref TEs), and what the characteristics (such as the number, type and distribution) are. To identify and characterize crop non-ref TEs, we analyzed non-ref TEs in more than 125 accessions from rice (Oryza sativa), maize (Zea mays) and sorghum (Sorghum bicolor) using resequenced data with paired-end mapping methods.

RESULTS

We identified 13,066, 23,866 and 35,679 non-ref TEs in rice, maize and sorghum, respectively. Genome-wide characterization analysis shows that most of non-ref TEs were unique and non-ref TE classes shows different among rice, maize and sorghum. We found that non-ref TEs have a strong positive correlation with gene number and have a bias toward insertion near genes, but with a preference for avoiding coding regions in maize and sorghum. The genes affected by non-ref TE insertion were functionally enriched for stress response mechanisms in all three crops.

CONCLUSIONS

These observations suggest that transposon insertion is not a random event and it makes genomic diversity, which may affect the intraspecific adaption and evolution of crops.

Stress-induced transposon reactivation: a mediator or an estimator of allostatic load?

Transposons are playing an important role in the evolution of eukaryotic genomes. These endogenous virus-like elements often amplify within their host genomes in a species specific manner. Today we have limited understanding when and how these amplification events happens. What we do know is that cells have evolved multiple line of defenses to keep these potentially invasive elements under control, often involving epigenetic mechanisms such as DNA-methylation and histone modifications. Emerging evidence shows a strong link between transposon activity and human aging and diseases, as well as a role for transposons in normal brain development. Controlling transposon activity may therefore uphold the fine balance between health and disease. In this article we investigate this balance, and sets it in relation to allostatic load, which conceptualize the link between stress and the "wear and tear" of the organism that leads to aging and disease. We hypothesize that stress-induced retrotransposon reactivation in humans may be used to estimate allostatic load, and may be a possible mechanism in which transposons amplify within species genomes.

Epigenetic mechanisms and associated brain circuits in the regulation of positive emotions: A role for transposable elements

Epigenetic programming and reprogramming are at the heart of cellular differentiation and represent developmental and evolutionary mechanisms in both germline and somatic cell lines. Only about 2% of our genome is composed of protein-coding genes, while the remaining 98%, once considered "junk" DNA, codes for regulatory/epigenetic elements that control how genes are expressed in different tissues and across time from conception to death. While we already know that epigenetic mechanisms are at play in cancer development and in regulating metabolism (cellular and whole body), the role of epigenetics in the developing prenatal and postnatal brain, and in maintaining a proper brain activity throughout the various stages of life, in addition to having played a critical role in human evolution, is a relatively new domain of knowledge. Here we present the current state-of-the-art techniques and results of these studies within the domain of emotions, and then speculate on how genomic and epigenetic mechanisms can modify and potentially alter our emotional (limbic) brain and affect our social interactions. J. Comp. Neurol. 524:2944-2954, 2016. © 2016 Wiley Periodicals, Inc.

Activation of LINE-1 Retrotransposon Increases the Risk of Epithelial-Mesenchymal Transition and Metastasis in Epithelial Cancer

Epithelial cancers comprise 80-90% of human cancers. During the process of cancer progression, cells lose their epithelial characteristics and acquire stem-like mesenchymal features that are resistant to chemotherapy. This process, termed the epithelial-mesenchymal transition (EMT), plays a critical role in the development of metastases. Because of the unique migratory and invasive properties of cells undergoing the EMT, therapeutic control of the EMT offers great hope and new opportunities for treating cancer. In recent years, a plethora of genes and noncoding RNAs, including miRNAs, have been linked to the EMT and the acquisition of stem cell-like properties. Despite these advances, questions remain unanswered about the molecular processes underlying such a cellular transition. In this article, we discuss how expression of the normally repressed LINE-1 (or L1) retrotransposons activates the process of EMT and the development of metastases. L1 is rarely expressed in differentiated stem cells or adult somatic tissues. However, its expression is widespread in almost all epithelial cancers and in stem cells in their undifferentiated state, suggesting a link between L1 activity and the proliferative and metastatic behaviour of cancer cells. We present an overview of L1 activity in cancer cells including how genes involved in proliferation, invasive and metastasis are modulated by L1 expression. The role of L1 in the differential expression of the let-7 family of miRNAs (that regulate genes involved in the EMT and metastasis) is also discussed. We also summarize recent novel insights into the role of the L1-encoded reverse transcriptase enzyme in epithelial cell plasticity that suggest it might be a potential therapeutic target that could reverse the EMT and the metastasis-associated stem cell-like properties of cancer cells.

LINE1 insertions as a genomic risk factor for schizophrenia: Preliminary evidence from an affected family

Recent studies show that human-specific LINE1s (L1HS) play a key role in the development of the central nervous system (CNS) and its disorders, and that their transpositions within the human genome are more common than previously thought. Many polymorphic L1HS, that is, present or absent across individuals, are not annotated in the current release of the genome and are customarily termed "non-reference L1s." We developed an analytical workflow to identify L1 polymorphic insertions with next-generation sequencing (NGS) using data from a family in which SZ segregates. Our workflow exploits two independent algorithms to detect non-reference L1 insertions, performs local de novo alignment of the regions harboring predicted L1 insertions and resolves the L1 subfamily designation from the de novo assembled sequence. We found 110 non-reference L1 polymorphic loci exhibiting Mendelian inheritance, the vast majority of which are already reported in dbRIP and/or euL1db, thus, confirming their status as non-reference L1 polymorphic insertions. Four previously undetected L1 polymorphic loci were confirmed by PCR amplification and direct sequencing of the insert. A large fraction of our non-reference L1s is located within the open reading frame of protein-coding genes that belong to pathways already implicated in the pathogenesis of schizophrenia. The finding of these polymorphic variants among SZ offsprings is intriguing and suggestive of putative pathogenic role. Our data show the utility of NGS to uncover L1 polymorphic insertions, a neglected type of genetic variants with the potential to influence the risk to develop schizophrenia like SNVs and CNVs. © 2016 Wiley Periodicals, Inc.

Activation of individual L1 retrotransposon instances is restricted to cell-type dependent permissive loci

LINE-1 (L1) retrotransposons represent approximately one sixth of the human genome, but only the human-specific L1HS-Ta subfamily acts as an endogenous mutagen in modern humans, reshaping both somatic and germline genomes. Due to their high levels of sequence identity and the existence of many polymorphic insertions absent from the reference genome, the transcriptional activation of individual genomic L1HS-Ta copies remains poorly understood. Here we comprehensively mapped fixed and polymorphic L1HS-Ta copies in 12 commonly-used somatic cell lines, and identified transcriptional and epigenetic signatures allowing the unambiguous identification of active L1HS-Ta copies in their genomic context. Strikingly, only a very restricted subset of L1HS-Ta loci - some being polymorphic among individuals - significantly contributes to the bulk of L1 expression, and these loci are differentially regulated among distinct cell lines. Thus, our data support a local model of L1 transcriptional activation in somatic cells, governed by individual-, locus-, and cell-type-specific determinants.

DOI:http://dx.doi.org/10.7554/eLife.13926.001

Retrotransposons, also known as jumping genes, have invaded the genomes of most living organisms. These fragments of DNA have the ability to move or copy themselves from one location of a chromosome to another. Depending on where they insert themselves, retrotransposons can modify the sequence of nearby genes, which can alter or even abolish their activity. Although these genetic modifications have contributed significantly to the evolution of different species, retrotransposons can also have detrimental effects; for example, by causing new cases of genetic diseases.

Adult human cells have a number of mechanisms that work to keep the activity of retrotransposons at a very low level. However, in many types of cancers retrotransposons escape these defense mechanisms and ‘jump’ actively. This is thought to contribute to the development and spread of cancerous tumors.

To understand how jumping genes are mobilized, a fundamental question must be answered: is the high jumping gene activity observed in some cell types a result of activating many copies of the retrotransposons, or only a few of them? This question has been difficult to address because there are more than one hundred copies of retrotransposons that could potentially move in humans, many of which have not even been referenced in the human genome map. Furthermore, each copy is almost identical to another one, making it difficult to discriminate between them.

Philippe et al. have now developed an approach that can map where individual retrotransposons are located in the genome of normal and cancerous cells and measure how active these jumping genes are. This revealed that only a very restricted number of them are active in any given cell type. Moreover, different subsets of jumping genes are active in different cell types, and their locations in the genome often do not overlap.

Thus, whether jumping genes are activated depends on the cell type and their position in the genome. These results are in contrast to the prevalent view that retrotransposons are activated in a more widespread manner across the genome, at least in cancerous cells. Overall, Philippe et al.’s findings pave the way towards characterizing the chromosome regions in which retrotransposons are frequently activated and understanding how they contribute to cancer and other diseases.

DOI:http://dx.doi.org/10.7554/eLife.13926.002

Germline Chromothripsis Driven by L1-Mediated Retrotransposition and Alu/Alu Homologous Recombination

Chromothripsis (CTH) is a phenomenon where multiple localized double-stranded DNA breaks result in complex genomic rearrangements. Although the DNA-repair mechanisms involved in CTH have been described, the mechanisms driving the localized "shattering" process remain unclear. High-throughput sequence analysis of a familial germline CTH revealed an inserted SVAE retrotransposon associated with a 110-kb deletion displaying hallmarks of L1-mediated retrotransposition. Our analysis suggests that the SVAE insertion did not occur prior to or after, but concurrent with the CTH event. We also observed L1-endonuclease potential target sites in other breakpoints. In addition, we found four Alu elements flanking the 110-kb deletion and associated with an inversion. We suggest that chromatin looping mediated by homologous Alu elements may have brought distal DNA regions into close proximity facilitating DNA cleavage by catalytically active L1-endonuclease. Our data provide the first evidence that active and inactive human retrotransposons can serve as endogenous mutagens driving CTH in the germline.

Reprogramming triggers endogenous L1 and Alu retrotransposition in human induced pluripotent stem cells

Human induced pluripotent stem cells (hiPSCs) are capable of unlimited proliferation and can differentiate in vitro to generate derivatives of the three primary germ layers. Genetic and epigenetic abnormalities have been reported by Wissing and colleagues to occur during hiPSC derivation, including mobilization of engineered LINE-1 (L1) retrotransposons. However, incidence and functional impact of endogenous retrotransposition in hiPSCs are yet to be established. Here we apply retrotransposon capture sequencing to eight hiPSC lines and three human embryonic stem cell (hESC) lines, revealing endogenous L1, Alu and SINE-VNTR-Alu (SVA) mobilization during reprogramming and pluripotent stem cell cultivation. Surprisingly, 4/7 de novo L1 insertions are full length and 6/11 retrotransposition events occurred in protein-coding genes expressed in pluripotent stem cells. We further demonstrate that an intronic L1 insertion in the CADPS2 gene is acquired during hiPSC cultivation and disrupts CADPS2 expression. These experiments elucidate endogenous retrotransposition, and its potential consequences, in hiPSCs and hESCs.

Genetic and epigenetic abnormalities have been found to result from reprogramming of differentiated cells into human induced pluripotent stem cells (hiPSCs). Here, Klawitter et al. identify endogenous L1, Alu and SVA mobilization during reprogramming, highlighting the risk of insertional mutagens in hiPSCs.

Mobster: accurate detection of mobile element insertions in next generation sequencing data

Mobile elements are major drivers in changing genomic architecture and can cause disease. The detection of mobile elements is hindered due to the low mappability of their highly repetitive sequences. We have developed an algorithm, called Mobster, to detect non-reference mobile element insertions in next generation sequencing data from both whole genome and whole exome studies. Mobster uses discordant read pairs and clipped reads in combination with consensus sequences of known active mobile elements. Mobster has a low false discovery rate and high recall rate for both L1 and Alu elements. Mobster is available at http://sourceforge.net/projects/mobster.

The Fine LINE: Methylation Drawing the Cancer Landscape

LINE-1 (L1) is the most abundant mammalian transposable element that comprises nearly 20% of the genome, and nearly half of the mammalian genome has stemmed from L1-mediated mobilization. Expression and retrotransposition of L1 are suppressed by complex mechanisms, where the key role belongs to DNA methylation. Alterations in L1 methylation may lead to aberrant expression of L1 and have been described in numerous diseases. Accumulating evidence clearly indicates that loss of global DNA methylation observed in cancer development and progression is tightly associated with hypomethylation of L1 elements. Significant progress achieved in the last several years suggests that such parameters as L1 methylation status can be potentially utilized as clinical biomarkers for determination of the disease stage and in predicting the disease-free survival in cancer patients. In this paper, we summarize the current knowledge on L1 methylation, with specific emphasis given to success and challenges on the way of introduction of L1 into clinical practice.