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

Comprehensive profiling of L1 retrotransposons in mouse

L1 elements are retrotransposons currently active in mammals. Although L1s are typically silenced in most normal tissues, elevated L1 expression is associated with a variety of conditions, including cancer, aging, infertility and neurological disease. These associations have raised interest in the mapping of human endogenous de novo L1 insertions, and a variety of methods have been developed for this purpose. Adapting these methods to mouse genomes would allow us to monitor endogenous in vivo L1 activity in controlled, experimental conditions using mouse disease models. Here, we use a modified version of transposon insertion profiling, called nanoTIPseq, to selectively enrich young mouse L1s. By linking this amplification step with nanopore sequencing, we identified >95% annotated L1s from C57BL/6 genomic DNA using only 200 000 sequencing reads. In the process, we discovered 82 unannotated L1 insertions from a single C57BL/6 genome. Most of these unannotated L1s were near repetitive sequence and were not found with short-read TIPseq. We used nanoTIPseq on individual mouse breast cancer cells and were able to identify the annotated and unannotated L1s, as well as new insertions specific to individual cells, providing proof of principle for using nanoTIPseq to interrogate retrotransposition activity at the single-cell level in vivo.

Regulation and function of transposable elements in cancer genomes

Over half of human genomic DNA is composed of repetitive sequences generated throughout evolution by prolific mobile genetic parasites called transposable elements (TEs). Long disregarded as "junk" or "selfish" DNA, TEs are increasingly recognized as formative elements in genome evolution, wired intimately into the structure and function of the human genome. Advances in sequencing technologies and computational methods have ushered in an era of unprecedented insight into how TE activity impacts human biology in health and disease. Here we discuss the current views on how TEs have shaped the regulatory landscape of the human genome, how TE activity is implicated in human cancers, and how recent findings motivate novel strategies to leverage TE activity for improved cancer therapy. Given the crucial role of methodological advances in TE biology, we pair our conceptual discussions with an in-depth review of the inherent technical challenges in studying repeats, specifically related to structural variation, expression analyses, and chromatin regulation. Lastly, we provide a catalog of existing and emerging assays and bioinformatic software that altogether are enabling the most sophisticated and comprehensive investigations yet into the regulation and function of interspersed repeats in cancer genomes.

Comprehensive profiling of L1 retrotransposons in mouse

L1 elements are retrotransposons currently active in mammals. Although L1s are typically silenced in most normal tissues, elevated L1 expression is associated with a variety of conditions, including cancer, aging, infertility, and neurological disease. These associations have raised interest in the mapping of human endogenous de novo L1 insertions, and a variety of methods have been developed for this purpose. Adapting these methods to mouse genomes would allow us to monitor endogenous in vivo L1 activity in controlled, experimental conditions using mouse disease models. Here we use a modified version of transposon insertion profiling, called nanoTIPseq, to selectively enrich young mouse L1s. By linking this amplification step with nanopore sequencing, we identified >95% annotated L1s from C57BL/6 genomic DNA using only 200,000 sequencing reads. In the process, we discovered 82 unannotated L1 insertions from a single C57BL/6 genome. Most of these unannotated L1s were near repetitive sequence and were not found with short-read TIPseq. We used nanoTIPseq on individual mouse breast cancer cells and were able to identify the annotated and unannotated L1s, as well as new insertions specific to individual cells, providing proof of principle for using nanoTIPseq to interrogate retrotransposition activity at the single cell level in vivo .

Transposon insertion profiling by sequencing (TIPseq) identifies novel LINE-1 insertions in human sperm

PURPOSE

Long interspersed nuclear element-1 (LINE-1 or L1) comprises 17% of the human genome. Retrotransposons may perturb gene integrity or alter gene expression by altering regulatory regions in the genome. The germline employs a number of mechanisms, including cytosine methylation, to repress retrotransposon transcription throughout most of life. Demethylation during germ cell and early embryo development de-represses retrotransposons. Intriguingly, de novo genetic variation appearing in sperm has been implicated in a number of disorders in offspring, including autism spectrum disorder, schizophrenia, and bipolar disorder. We hypothesize that human sperm exhibit de novo retrotransposition and employ a new sequencing method, single cell transposon insertion profiling by sequencing (scTIPseq) to map them in small amounts of human sperm.

METHODS

Cross-sectional case-control study of sperm samples (n=10 men; ages 32-55 years old) from consenting men undergoing IVF at NYU Langone Fertility Center. scTIPseq identified novel LINE-1 insertions in individual sperm and TIPseqHunter, a custom bioinformatics pipeline, compared the architecture of sperm LINE-1 to known LINE-1 insertions from the European database of Human specific LINE-1 (L1Hs) retrotransposon insertions (euL1db).

RESULTS

scTIPseq identified 17 novel insertions in sperm. New insertions were mainly intergenic or intronic. Only one sample did not exhibit new insertions. The location or number of novel insertions did not differ by paternal age.

CONCLUSION

This study for the first time reports novel LINE-1 insertions in human sperm, demonstrating the feasibility of scTIPseq, and identifies new contributors to genetic diversity in the human germ line.

TDP-43 safeguards the embryo genome from L1 retrotransposition

Transposable elements (TEs) are genomic parasites that propagate within the host genome and introduce mutations. Long interspersed nuclear element-1 (LINE-1 or L1) is the major TE class, which occupies nearly 20% of the mouse genome. L1 is highly active in mammalian preimplantation embryos, posing a major threat to genome integrity, but the mechanism of stage-specific protection against L1 retrotransposition is unknown. Here, we show that TAR DNA-binding protein 43 (TDP-43), mutations in which constitute a major risk factor for amyotrophic lateral sclerosis, inhibits L1 retrotransposition in mouse embryonic stem cells (mESCs) and preimplantation embryos. Knockdown of TDP-43 resulted in massive genomic L1 expansion and impaired cell growth in preimplantation embryos and ESCs. Functional analysis demonstrated that TDP-43 interacts with L1 open reading frame 1 protein (L1 ORF1p) to mediate genomic protection, and loss of this interaction led to derepression of L1 retrotransposition. Our results identify TDP-43 as a guardian of the embryonic genome.

Transposable element polymorphisms improve prediction of complex agronomic traits in rice

Transposon insertion polymorphisms can improve prediction of complex agronomic traits in rice compared to using SNPs only, especially when accessions to be predicted are less related to the training set. Transposon insertion polymorphisms (TIPs) are significant sources of genetic variation. Previous work has shown that TIPs can improve detection of causative loci on agronomic traits in rice. Here, we quantify the fraction of variance explained by single nucleotide polymorphisms (SNPs) compared to TIPs, and we explore whether TIPs can improve prediction of traits when compared to using only SNPs. We used eleven traits of agronomic relevance from by five different rice population groups (Aus, Indica, Aromatic, Japonica, and Admixed), 738 accessions in total. We assess prediction by applying data split validation in two scenarios. In the within-population scenario, we predicted performance of improved Indica varieties using the rest of Indica accessions. In the across population scenario, we predicted all Aromatic and Admixed accessions using the rest of populations. In each scenario, Bayes C and a Bayesian reproducible kernel Hilbert space regression were compared. We find that TIPs can explain an important fraction of total genetic variance and that they also improve genomic prediction. In the across population prediction scenario, TIPs outperformed SNPs in nine out of the eleven traits analyzed. In some traits like leaf senescence or grain width, using TIPs increased predictive correlation by 30-50%. Our results evidence, for the first time, that TIPs genotyping can improve prediction on complex agronomic traits in rice, especially when accessions to be predicted are less related to training accessions.

Somatic Mobilization: High Somatic Insertion Rate of mariner Transposable Element in Drosophila simulans

Although transposable elements (TEs) are usually silent in somatic tissues, they are sometimes mobilized in the soma and can potentially have biological consequences. The mariner element is one of the TEs involved in somatic mobilization (SM) in Drosophila and has a high rate of somatic excision. It is also known that temperature is an important factor in the increase of the mariner element SM in the fly. However, it is important to emphasize that excision is only one step of TE transposition, and the final step in this process is insertion. In the present study, we used an assay based on sequencing of the mariner flanking region and developed a pipeline to identify novel mariner insertions in Drosophila simulans at 20 and 28 °C. We found that flies carrying two mariner copies (one autonomous and one non-autonomous) had an average of 236.4 (±99.3) to 279 (±107.7) new somatic insertions at 20 °C and an average of 172.7 (±95.3) to 252.6 (±67.3) at 28 °C. In addition, we detected fragments containing mariner and others without mariner in the same regions with low-coverage long-read sequencing, indicating the process of excision and insertion. In conclusion, a low number of autonomous copies of the mariner transposon can promote a high rate of new somatic insertions during the developmental stages of Drosophila. Additionally, the developed method seems to be sensitive and adequate for the verification and estimation of somatic insertion.

Transposable elements that have recently been mobile in the human genome

Background

Transposable elements (TE) comprise nearly half of the human genome and their insertions have profound effects to human genetic diversification and as well as disease. Despite their abovementioned significance, there is no consensus on the TE subfamilies that remain active in the human genome. In this study, we therefore developed a novel statistical test for recently mobile subfamilies (RMSs), based on patterns of overlap with > 100,000 polymorphic indels.

Results

Our analysis produced a catalogue of 20 high-confidence RMSs, which excludes many false positives in public databases. Intriguingly though, it includes HERV-K, an LTR subfamily previously thought to be extinct. The RMS catalogue is strongly enriched for contributions to germline genetic disorders (P = 1.1e-10), and thus constitutes a valuable resource for diagnosing disorders of unknown aetiology using targeted TE-insertion screens. Remarkably, RMSs are also highly enriched for somatic insertions in diverse cancers (P = 2.8e-17), thus indicating strong correlations between germline and somatic TE mobility. Using CRISPR/Cas9 deletion, we show that an RMS-derived polymorphic TE insertion increased the expression of RPL17, a gene associated with lower survival in liver cancer. More broadly, polymorphic TE insertions from RMSs were enriched near genes with allele-specific expression, suggesting widespread effects on gene regulation.

Conclusions

By using a novel statistical test we have defined a catalogue of 20 recently mobile transposable element subfamilies. We illustrate the gene regulatory potential of RMS-derived polymorphic TE insertions, using CRISPR/Cas9 deletion in vitro on a specific candidate, as well as by genome wide analysis of allele-specific expression. Our study presents novel insights into TE mobility and regulatory potential and provides a key resource for human disease genetics and population history studies.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-08085-0.

Genomic approaches to trace the history of human brain evolution with an emerging opportunity for transposon profiling of ancient humans

Transposable elements (TEs) significantly contribute to shaping the diversity of the human genome, and lines of evidence suggest TEs as one of driving forces of human brain evolution. Existing computational approaches, including cross-species comparative genomics and population genetic modeling, can be adapted for the study of the role of TEs in evolution. In particular, diverse ancient and archaic human genome sequences are increasingly available, allowing reconstruction of past human migration events and holding the promise of identifying and tracking TEs among other evolutionarily important genetic variants at an unprecedented spatiotemporal resolution. However, highly degraded short DNA templates and other unique challenges presented by ancient human DNA call for major changes in current experimental and computational procedures to enable the identification of evolutionarily important TEs. Ancient human genomes are valuable resources for investigating TEs in the evolutionary context, and efforts to explore ancient human genomes will potentially provide a novel perspective on the genetic mechanism of human brain evolution and inspire a variety of technological and methodological advances. In this review, we summarize computational and experimental approaches that can be adapted to identify and validate evolutionarily important TEs, especially for human brain evolution. We also highlight strategies that leverage ancient genomic data and discuss unique challenges in ancient transposon genomics.

Intronic LINE-1 insertion in SLCO1B3 as a highly prevalent cause of rotor syndrome in East Asian population

Rotor syndrome is caused by digenic loss-of-function variants in SLCO1B1 and SLCO1B3 but only a few studies have reported co-occurring inactivating variants from both genes. A rotor syndrome-causing long interspersed element-1 (LINE-1) insertion in SLCO1B3 had been reported to be highly prevalent in the Japanese population but there has been no additional report. In spite of its known association with various human diseases, LINE-1 is hard to detect with current sequencing technologies. In this study, we aimed to devise a method to screen the LINE-1 insertion variant and investigate the frequency of this variant in various populations. A chimeric sequence, that was generated by concatenating the reference sequence at the junction and a part of inserted LINE-1 sequence, was searched from 725 raw sequencing data files. In cases containing the chimeric sequence, confirmatory long-range PCR and gap-PCR were performed. In total, 95 (13.1%) of 725 patients were positive for the chimeric sequence, and all were confirmed to have the SLCO1B3 LINE-1 insertion by PCR-based tests. The same chimeric sequence was searched from the 1000 Genomes Project data repository and the carrier frequency was remarkably high in the East Asian populations (10.1%), especially in Southern Han Chinese (18.5%), but almost absent in other populations. This SLCO1B3 LINE-1 insertion should be screened in a population-specific manner under suspicion of Rotor syndrome and the methods proposed in this study would enable this in a simple way.

Cas9 targeted enrichment of mobile elements using nanopore sequencing

Mobile element insertions (MEIs) are repetitive genomic sequences that contribute to genetic variation and can lead to genetic disorders. Targeted and whole-genome approaches using short-read sequencing have been developed to identify reference and non-reference MEIs; however, the read length hampers detection of these elements in complex genomic regions. Here, we pair Cas9-targeted nanopore sequencing with computational methodologies to capture active MEIs in human genomes. We demonstrate parallel enrichment for distinct classes of MEIs, averaging 44% of reads on-targeted signals and exhibiting a 13.4-54x enrichment over whole-genome approaches. We show an individual flow cell can recover most MEIs (97% L1Hs, 93% AluYb, 51% AluYa, 99% SVA_F, and 65% SVA_E). We identify seventeen non-reference MEIs in GM12878 overlooked by modern, long-read analysis pipelines, primarily in repetitive genomic regions. This work introduces the utility of nanopore sequencing for MEI enrichment and lays the foundation for rapid discovery of elusive, repetitive genetic elements.

RIP-seq reveals LINE-1 ORF1p association with p-body enriched mRNAs

BACKGROUND

Long INterspersed Element-1 (LINE-1) is an autonomous retroelement able to "copy-and-paste" itself into new loci of the host genome through a process called retrotransposition. The LINE-1 bicistronic mRNA codes for two proteins, ORF1p, a nucleic acid chaperone, and ORF2p, a protein with endonuclease and reverse transcriptase activity. Both proteins bind LINE-1 mRNA in cis and are necessary for retrotransposition. While LINE-1 transcription is usually repressed in most healthy somatic cells through a plethora of mechanisms, ORF1p expression has been observed in nearly 50% of tumors, and new LINE-1 insertions have been documented in a similar fraction of tumors, including prostate cancer.

RESULTS

Here, we utilized RNA ImmunoPrecipitation (RIP) and the L1EM analysis software to identify ORF1p bound RNA in prostate cancer cells. We identified LINE-1 loci that were expressed in parental androgen sensitive and androgen independent clonal derivatives. In all androgen independent cells, we found higher levels of LINE-1 RNA, as well as unique expression patterns of LINE-1 loci. Interestingly, we observed that ORF1p bound many non-LINE-1 mRNA in all prostate cancer cell lines evaluated, and polyA RNA, and RNA localized in p-bodies were especially enriched. Furthermore, the expression levels of RNAs identified in our ORF1p RIP correlated with RNAs expressed in LINE-1 positive tumors from The Cancer Genome Atlas (TCGA).

CONCLUSION

Our results show a significant remodeling of LINE-1 loci expression in androgen independent cell lines when compared to parental androgen dependent cells. Additionally, we found that ORF1p bound a significant amount of non-LINE-1 mRNA, and that the enriched ORF1p bound mRNAs are also amplified in LINE-1 expressing TCGA prostate tumors, indicating the biological relevance of our findings to prostate cancer.

The genetics of ductal adenocarcinoma of the pancreas in the year 2020: dramatic progress, but far to go

The publication of the "Pan-Cancer Atlas" by the Pan-Cancer Analysis of Whole Genomes Consortium, a partnership formed by The Cancer Genome Atlas (TCGA) and International Cancer Genome Consortium (ICGC), provides a wonderful opportunity to reflect on where we stand in our understanding of the genetics of pancreatic cancer, as well as on the opportunities to translate this understanding to patient care. From germline variants that predispose to the development of pancreatic cancer, to somatic mutations that are therapeutically targetable, genetics is now providing hope, where there once was no hope, for those diagnosed with pancreatic cancer.

TranspoScope: interactive visualization of retrotransposon insertions

MOTIVATION

Retrotransposition is an important force in shaping the human genome and is involved in prenatal development, disease and aging. Current genome browsers are not optimized for visualizing the experimental evidence for retrotransposon insertions.

RESULTS

We have developed a specialized browser to visualize the evidence for retrotransposon insertions for both targeted and whole-genome sequencing data.

AVAILABILITY AND IMPLEMENTATION

TranspoScope's source code, as well as installation instructions, are available at https://github.com/FenyoLab/transposcope.

Genomic Indexing by Somatic Gene Recombination of mRNA/ncRNA - Does It Play a Role in Genomic Mosaicism, Memory Formation, and Alzheimer's Disease?

Recent evidence indicates that genomic individuality of neurons, characterized by DNA-content variation, is a common if not universal phenomenon in the human brain that occurs naturally but can also show aberrancies that have been linked to the pathomechanism of Alzheimer’s disease and related neurodegenerative disorders. Etiologically, this genomic mosaic has been suggested to arise from defects of cell cycle regulation that may occur either during brain development or in the mature brain after terminal differentiation of neurons. Here, we aim to draw attention towards another mechanism that can give rise to genomic individuality of neurons, with far-reaching consequences. This mechanism has its origin in the transcriptome rather than in replication defects of the genome, i.e., somatic gene recombination of RNA. We continue to develop the concept that somatic gene recombination of RNA provides a physiological process that, through integration of intronless mRNA/ncRNA into the genome, allows a particular functional state at the level of the individual neuron to be indexed. By insertion of defined RNAs in a somatic recombination process, the presence of specific mRNA transcripts within a definite temporal context can be “frozen” and can serve as an index that can be recalled at any later point in time. This allows information related to a specific neuronal state of differentiation and/or activity relevant to a memory trace to be fixed. We suggest that this process is used throughout the lifetime of each neuron and might have both advantageous and deleterious consequences.

Integrated Mobile Element Scanning (ME-Scan) method for identifying multiple types of polymorphic mobile element insertions

Background

Mobile elements are ubiquitous components of mammalian genomes and constitute more than half of the human genome. Polymorphic mobile element insertions (pMEIs) are a major source of human genomic variation and are gaining research interest because of their involvement in gene expression regulation, genome integrity, and disease.

Results

Building on our previous Mobile Element Scanning (ME-Scan) protocols, we developed an integrated ME-Scan protocol to identify three major active families of human mobile elements, AluYb, L1HS, and SVA. This approach selectively amplifies insertion sites of currently active retrotransposons for Illumina sequencing. By pooling the libraries together, we can identify pMEIs from all three mobile element families in one sequencing run. To demonstrate the utility of the new ME-Scan protocol, we sequenced 12 human parent-offspring trios. Our results showed high sensitivity (> 90%) and accuracy (> 95%) of the protocol for identifying pMEIs in the human genome. In addition, we also tested the feasibility of identifying somatic insertions using the protocol.

Conclusions

The integrated ME-Scan protocol is a cost-effective way to identify novel pMEIs in the human genome. In addition, by developing the protocol to detect three mobile element families, we demonstrate the flexibility of the ME-Scan protocol. We present instructions for the library design, a sequencing protocol, and a computational pipeline for downstream analyses as a complete framework that will allow researchers to easily adapt the ME-Scan protocol to their own projects in other genomes.

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'.

Cell fitness screens reveal a conflict between LINE-1 retrotransposition and DNA replication

LINE-1 retrotransposon overexpression is a hallmark of human cancers. We identified a colorectal cancer wherein a fast-growing tumor subclone downregulated LINE-1, prompting us to examine how LINE-1 expression affects cell growth. We find that nontransformed cells undergo a TP53-dependent growth arrest and activate interferon signaling in response to LINE-1. TP53 inhibition allows LINE-1+ cells to grow, and genome-wide-knockout screens show that these cells require replication-coupled DNA-repair pathways, replication-stress signaling and replication-fork restart factors. Our findings demonstrate that LINE-1 expression creates specific molecular vulnerabilities and reveal a retrotransposition-replication conflict that may be an important determinant of cancer growth.

LINE-1 ORF2p expression is nearly imperceptible in human cancers

BACKGROUND

Long interspersed element-1 (LINE-1, L1) is the major driver of mobile DNA activity in modern humans. When expressed, LINE-1 loci produce bicistronic transcripts encoding two proteins essential for retrotransposition, ORF1p and ORF2p. Many types of human cancers are characterized by L1 promoter hypomethylation, L1 transcription, L1 ORF1p protein expression, and somatic L1 retrotransposition. ORF2p encodes the endonuclease and reverse transcriptase activities required for L1 retrotransposition. Its expression is poorly characterized in human tissues and cell lines.

RESULTS

We report mass spectrometry-based tumor proteome profiling studies wherein ORF2p eludes detection. To test whether ORF2p could be detected with specific reagents, we developed and validated five rabbit monoclonal antibodies with immunoreactivity for specific epitopes on the protein. These reagents readily detect ectopic ORF2p expressed from bicistronic L1 constructs. However, endogenous ORF2p is not detected in human tumor samples or cell lines by western blot, immunoprecipitation, or immunohistochemistry despite high levels of ORF1p expression. Moreover, we report endogenous ORF1p-associated interactomes, affinity isolated from colorectal cancers, wherein we similarly fail to detect ORF2p. These samples include primary tumors harboring hundreds of somatically acquired L1 insertions. The new data are available via ProteomeXchange with identifier PXD013743.

CONCLUSIONS

Although somatic retrotransposition provides unequivocal genetic evidence for the expression of ORF2p in human cancers, we are unable to directly measure its presence using several standard methods. Experimental systems have previously indicated an unequal stoichiometry between ORF1p and ORF2p, but in vivo, the expression of these two proteins may be more strikingly uncoupled. These findings are consistent with observations that ORF2p is not tolerable for cell growth.

Diseases of the nERVous system: retrotransposon activity in neurodegenerative disease

Transposable Elements (TEs) are mobile genetic elements whose sequences constitute nearly half of the human genome. Each TE copy can be present in hundreds to thousands of locations within the genome, complicating the genetic and genomic studies of these highly repetitive sequences. The recent development of better tools for evaluating TE derived sequences in genomic studies has enabled an increasing appreciation for the contribution of TEs to human development and disease. While some TEs have contributed novel and beneficial host functions, this review will summarize the evidence for detrimental TE activity in neurodegenerative disorders. Much of the evidence for pathogenicity implicates endogenous retroviruses (ERVs), a subset of TEs that entered the genome by retroviral infections of germline cells in our evolutionary ancestors and have since been passed down as a substantial fraction of the human genome. Human specific ERVs (HERVs) represent some of the youngest ERVs in the genome, and thus are presumed to retain greater function and resultant pathogenic potential.