Natural mutagenesis of human genomes by endogenous retrotransposons

Altering Genomic Integrity: Heavy Metal Exposure Promotes Transposable Element-Mediated Damage

Maintenance of genomic integrity is critical for cellular homeostasis and survival. The active transposable elements (TEs) composed primarily of three mobile element lineages LINE-1, Alu, and SVA comprise approximately 30% of the mass of the human genome. For the past 2 decades, studies have shown that TEs significantly contribute to genetic instability and that TE-caused damages are associated with genetic diseases and cancer. Different environmental exposures, including several heavy metals, influence how TEs interact with its host genome increasing their negative impact. This mini-review provides some basic knowledge on TEs, their contribution to disease, and an overview of the current knowledge on how heavy metals influence TE-mediated damage.

Response of transposable elements to environmental stressors

Transposable elements (TEs) comprise a group of repetitive sequences that bring positive, negative, as well as neutral effects to the host organism. Earlier considered as "junk DNA," TEs are now well-accepted driving forces of evolution and critical regulators of the expression of genetic information. Their activity is regulated by epigenetic mechanisms, including methylation of DNA and histone modifications. The loss of epigenetic control over TEs, exhibited as loss of DNA methylation and decondensation of the chromatin structure, may result in TEs reactivation, initiation of their insertional mutagenesis (retrotransposition) and has been reported in numerous human diseases, including cancer. Accumulating evidence suggests that these alterations are not the simple consequences of the disease, but often may drive the pathogenesis, as they can be detected early during disease development. Knowledge derived from the in vitro, in vivo, and epidemiological studies, clearly demonstrates that exposure to ubiquitous environmental stressors, many of which are carcinogens or suspected carcinogens, are capable of causing alterations in methylation and expression of TEs and initiate retrotransposition events. Evidence summarized in this review suggests that TEs are the sensitive endpoints for detection of effects caused by such environmental stressors, as ionizing radiation (terrestrial, space, and UV-radiation), air pollution (including particulate matter [PM]-derived and gaseous), persistent organic pollutants, and metals. Furthermore, the significance of these effects is characterized by their early appearance, persistence and presence in both, target organs and peripheral blood. Altogether, these findings suggest that TEs may potentially be introduced into safety and risk assessment and serve as biomarkers of exposure to environmental stressors. Furthermore, TEs also show significant potential to become invaluable surrogate biomarkers in clinic and possible targets for therapeutic modalities for disease treatment and prevention.

Regulatory roles of LINE-1-encoded reverse transcriptase in cancer onset and progression

LINE-1 retrotransposons encode the reverse transcriptase (RT) enzyme, required for their own mobility, the expression of which is inhibited in differentiated tissues while being active in tumors. Experimental evidence indicate that the inhibition of LINE-1-derived RT restores differentiation in cancer cells, inhibits tumor progression and yields globally reprogrammed transcription profiles. Newly emerging data suggest that LINE-1-encoded RT modulates the biogenesis of miRNAs, by governing the balance between the production of regulatory double-stranded RNAs and RNA:DNA hybrid molecules, with a direct impact on global gene expression. Abnormally high RT activity unbalances the transcriptome in cancer cells, while RT inhibition restores ‘normal’ miRNA profiles and their regulatory networks. This RT-dependent mechanism can target the myriad of transcripts - both coding and non-coding, sense and antisense - in eukaryotic transcriptomes, with a profound impact on cell fates. LINE-1-encoded RT emerges therefore as a key regulator of a previously unrecognized mechanism in tumorigenesis

The potential role of Alu Y in the development of resistance to SN38 (Irinotecan) or oxaliplatin in colorectal cancer

BACKGROUND

Irinotecan (SN38) and oxaliplatin are chemotherapeutic agents used in the treatment of colorectal cancer. However, the frequent development of resistance to these drugs represents a considerable challenge in the clinic. Alus as retrotransposons comprise 11% of the human genome. Genomic toxicity induced by carcinogens or drugs can reactivate Alus by altering DNA methylation. Whether or not reactivation of Alus occurs in SN38 and oxaliplatin resistance remains unknown.

RESULTS

We applied reduced representation bisulfite sequencing (RRBS) to investigate the DNA methylome in SN38 or oxaliplatin resistant colorectal cancer cell line models. Moreover, we extended the RRBS analysis to tumor tissue from 14 patients with colorectal cancer who either did or did not benefit from capecitabine + oxaliplatin treatment. For the clinical samples, we applied a concept of 'DNA methylation entropy' to estimate the diversity of DNA methylation states of the identified resistance phenotype-associated methylation loci observed in the cell line models. We identified different loci being characteristic for the different resistant cell lines. Interestingly, 53% of the identified loci were Alu sequences- especially the Alu Y subfamily. Furthermore, we identified an enrichment of Alu Y sequences that likely results from increased integration of new copies of Alu Y sequence in the drug-resistant cell lines. In the clinical samples, SOX1 and other SOX gene family members were shown to display variable DNA methylation states in their gene regions. The Alu Y sequences showed remarkable variation in DNA methylation states across the clinical samples.

CONCLUSION

Our findings imply a crucial role of Alu Y in colorectal cancer drug resistance. Our study underscores the complexity of colorectal cancer aggravated by mobility of Alu elements and stresses the importance of personalized strategies, using a systematic and dynamic view, for effective cancer therapy.

Short-term exposure to engineered nanomaterials affects cellular epigenome

Extensive incorporation of engineered nanomaterials (ENMs) into industrial and biomedical applications increases the risks of exposure to these potentially hazardous materials. While the geno- and cytotoxic effects of ENMs have been investigated, the potential of ENMs to target the cellular epigenome remains largely unknown. Our goal was to determine whether industry relevant ENMs can affect the epigenome at low cytotoxic doses. A panel of cells relevant to inhalation exposures such as human and murine macrophages (THP-1 and RAW264.7, respectively) and human small airway epithelial cells (SAEC) were exposed to printer-emitted engineered nanoparticles (PEPs), mild steel welding fumes (MS-WF), copper oxide (CuO) and titanium dioxide nanoparticles. Toxicological effects, including cytotoxicity, oxidative stress and inflammatory responses were assessed, taking into consideration in vitro dosimetry. The effects of ENMs on cellular epigenome were determined by addressing the global and transposable elements (TEs)-associated DNA methylation and expression of DNA methylation machinery and TEs. The percentage of ENMs-induced cytotoxicity for all cell lines was in the range of 0-15%. Oxidative stress was evident in SAEC after exposure to PEPs and in THP-1 when exposed to CuO. In addition, exposure to ENMs resulted in modest alterations in DNA methylation of two most abundant TEs in mammalian genomes, LINE-1 and Alu/SINE, their transcriptional reactivation, and decreased expression of DNA methylation machinery in a cell-, dose- and ENM-dependent manner. These results indicate that exposure to ENMs at environmentally relevant concentrations, aside from the geno- and cytotoxic effects, can also affect the epigenome of target cells.

Genome-wide amplification of proviral sequences reveals new polymorphic HERV-K(HML-2) proviruses in humans and chimpanzees that are absent from genome assemblies

BACKGROUND

To date, the human population census of proviruses of the Betaretrovirus-like human endogenous retroviral (HERV-K) (HML-2) family has been compiled from a limited number of complete genomes, making it certain that rare polymorphic loci are under-represented and are yet to be described.

RESULTS

Here we describe a suppression PCR-based method called genome-wide amplification of proviral sequences (GAPS) that selectively amplifies DNA fragments containing the termini of HERV-K(HML-2) proviral sequences and their flanking genomic sequences. We analysed the HERV-K(HML-2) proviral content of 101 unrelated humans, 4 common chimpanzees and three centre d'etude du polymorphisme humain (CEPH) pedigrees (44 individuals). The technique isolated HERV-K(HML-2) proviruses that had integrated in the genomes of the great apes throughout their divergence and included evolutionarily young elements still unfixed for presence/absence.

CONCLUSIONS

By examining the HERV-K(HML-2) proviral content of 145 humans we detected a new insertionally polymorphic Type I HERV-K(HML-2) provirus. We also observed provirus versus solo long terminal repeat (LTR) polymorphism within humans at a previously unreported, but ancient, locus. Finally, we report two novel chimpanzee specific proviruses, one of which is dimorphic for a provirus versus solo LTR. Thus GAPS enables the isolation of uncharacterised HERV-K(HML-2) proviral sequences and provides a direct means to assess inter-individual genetic variation associated with HERV-K(HML-2) proviruses.

Silencing of LINE-1 retrotransposons contributes to variation in small noncoding RNA expression in human cancer cells

Noncoding RNAs are key players in the maintenance of genomic integrity, particularly in silencing the expression of repetitive elements, some of which are retrotransposable and capable of causing genomic instability. Recent computational studies suggest an association between L1 expression and the generation of small RNAs. However, whether L1 expression has a role in the activation of small RNA expression has yet to be determined experimentally.; Here we report a global analysis of small RNAs in deep sequencing from L1-active and L1-silenced breast cancer cells. We found that cells in which L1 expression was silenced exhibited greatly increased expression of a number of miRNAs and in particular, members of the let-7 family. In addition, we found differential expression of a few piRNAs that might potentially regulate gene expression. We also report the identification of several repeat RNAs against LTRs, LINEs and SINE elements. Although most of the repeat RNAs mapped to L1 elements, in general we found no significant differences in the expression levels of repeat RNAs in the presence or absence of L1 expression except for a few RNAs targeting subclasses of L1 elements. These differentially expressed small RNAs may function in human genome defence responses.

piRNA involvement in genome stability and human cancer

PIWI-interacting RNAs (piRNAs) are a large family of small, single-stranded, non-coding RNAs present throughout the animal kingdom. They form complexes with several members of the PIWI clade of Argonaute proteins and carry out regulatory functions. Their best established biological role is the inhibition of transposon mobilization, which they enforce both at the transcriptional level, through regulation of heterochromatin formation, and by promoting transcript degradation. In this capacity, piRNAs and PIWI proteins are at the heart of the germline cells’ efforts to preserve genome integrity. Additional regulatory roles of piRNAs and PIWI proteins in gene expression are becoming increasingly apparent.

PIWI proteins and piRNAs are often detected in human cancers deriving from germline cells as well as somatic tissues. Their detection in cancer correlates with poorer clinical outcomes, suggesting that they play a functional role in the biology of cancer. Nonetheless, the currently available information, while highly suggestive, is still not sufficient to entirely discriminate between a ‘passenger’ role for the ectopic expression of piRNAs and PIWI proteins in cancer from a ‘driver’ role in the pathogenesis of these diseases. In this article, we review some of the key available evidence for the role of piRNAs and PIWI in human cancer and discuss ways in which our understanding of their functions may be improved.

Phosphorylation of ORF1p is required for L1 retrotransposition

Although members of the L1 (LINE-1) clade of non-LTR retrotransposons can be deleterious, the L1 clade has remained active in most mammals for ∼100 million years and generated almost 40% of the human genome. The details of L1-host interaction are largely unknown, however. Here we report that L1 activity requires phosphorylation of the protein encoded by the L1 ORF1 (ORF1p). Critical phospho-acceptor residues (two serines and two threonines) reside in four conserved proline-directed protein kinase (PDPK) target sites. The PDPK family includes mitogen-activated protein kinases and cyclin-dependent kinases. Mutation of any PDPK phospho-acceptor inhibits L1 retrotransposition. The phosphomimetic aspartic acid can restore activity at the two serine sites, but not at either threonine site, where it is strongly inhibitory. ORF1p also contains conserved PDPK docking sites, which promote specific interaction of PDPKs with their targets. As expected, mutations in these sites also inhibit L1 activity. PDPK mutations in ORF1p that inactivate L1 have no significant effect on the ability of ORF1p to anneal RNA in vitro, an important biochemical property of the protein. We show that phosphorylated PDPK sites in ORF1p are required for an interaction with the peptidyl prolyl isomerase 1 (Pin1), a critical component of PDPK-mediated regulation. Pin1 acts via isomerization of proline side chains at phosphorylated PDPK motifs, thereby affecting substrate conformation and activity. Our demonstration that L1 activity is dependent on and integrated with cellular phosphorylation regulatory cascades significantly increases our understanding of interactions between L1 and its host.

Sequencing, identification and mapping of primed L1 elements (SIMPLE) reveals significant variation in full length L1 elements between individuals

BACKGROUND

There are over a half a million copies of L1 retroelements in the human genome which are responsible for as much as 0.5% of new human genetic diseases. Most new L1 inserts arise from young source elements that are polymorphic in the human genome. Highly active polymorphic "hot" L1 source elements have been shown to be capable of extremely high levels of mobilization and result in numerous instances of disease. Additionally, hot polymorphic L1s have been described to be highly active within numerous cancer genomes. These hot L1s result in mutagenesis by insertion of new L1 copies elsewhere in the genome, but also have been shown to generate additional full length L1 insertions which are also hot and able to further retrotranspose. Through this mechanism, hot L1s may amplify within a tumor and result in a continued cycle of mutagenesis.

RESULTS AND CONCLUSIONS

We have developed a method to detect full-length, polymorphic L1 elements using a targeted next generation sequencing approach, Sequencing Identification and Mapping of Primed L1 Elements (SIMPLE). SIMPLE has 94% sensitivity and detects nearly all full-length L1 elements in a genome. SIMPLE will allow researchers to identify hot mutagenic full-length L1s as potential drivers of genome instability. Using SIMPLE we find that the typical individual has approximately 100 non-reference, polymorphic L1 elements in their genome. These elements are at relatively low population frequencies relative to previously identified polymorphic L1 elements and demonstrate the tremendous diversity in potentially active L1 elements in the human population.

U6 snRNA Pseudogenes: Markers of Retrotransposition Dynamics in Mammals

Transposable elements comprise more than 45% of the human genome and long interspersed nuclear element 1 (LINE-1 or L1) is the only autonomous mobile element remaining active. Since its identification, it has been proposed that L1 contributes to the mobilization and amplification of other cellular RNAs and more recently, experimental demonstrations of this function has been described for many transcripts such as Alu, a nonautonomous mobile element, cellular mRNAs, or small noncoding RNAs. Detailed examination of the mobilization of various cellular RNAs revealed distinct pathways by which they could be recruited during retrotransposition; template choice or template switching. Here, by analyzing genomic structures and retrotransposition signatures associated with small nuclear RNA (snRNA) sequences, we identified distinct recruiting steps during the L1 retrotransposition cycle for the formation of snRNA-processed pseudogenes. Interestingly, some of the identified recruiting steps take place in the nucleus. Moreover, after comparison to other vertebrate genomes, we established that snRNA amplification by template switching is common to many LINE families from several LINE clades. Finally, we suggest that U6 snRNA copies can serve as markers of L1 retrotransposition dynamics in mammalian genomes.

Landscape of insertion polymorphisms in the human genome

Nucleotide substitutions, small (<50 bp) insertions or deletions (indels), and large (>50 bp) deletions are well-known causes of genetic variation within the human genome. We recently reported a previously unrecognized form of polymorphic insertions, termed templated sequence insertion polymorphism (TSIP), in which the inserted sequence was templated from a distant genomic region, and was inserted in the genome through reverse transcription of an RNA intermediate. TSIPs can be grouped into two classes based on nucleotide sequence features at the insertion junctions; class 1 TSIPs show target site duplication, polyadenylation, and preference for insertion at a 5′-TTTT/A-3′ sequence, suggesting a LINE-1 based insertion mechanism, whereas class 2 TSIPs show features consistent with repair of a DNA double strand break by nonhomologous end joining. To gain a more complete picture of TSIPs throughout the human population, we evaluated whole-genome sequence from 52 individuals, and identified 171 TSIPs. Most individuals had 25–30 TSIPs, and common (present in >20% of individuals) TSIPs were found in individuals throughout the world, whereas rare TSIPs tended to cluster in specific geographic regions. The number of rare TSIPs was greater than the number of common TSIPs, suggesting that TSIP generation is an ongoing process. Intriguingly, mitochondrial sequences were a frequent template for class 2 insertions, used more commonly than any nuclear chromosome. Similar to single nucleotide polymorphisms and indels, we suspect that these TSIPs may be important for the generation of human diversity and genetic diseases, and can be useful in tracking historical migration of populations.

Association of new deletion/duplication region at chromosome 1p21 with intellectual disability, severe speech deficit and autism spectrum disorder-like behavior: an all-in approach to solving the DPYD enigma

We describe an as yet unreported neocentric small supernumerary marker chromosome (sSMC) derived from chromosome 1p21.3p21.2. It was present in 80% of the lymphocytes in a male patient with intellectual disability, severe speech deficit, mild dysmorphic features, and hyperactivity with elements of autism spectrum disorder (ASD).

Several important neurodevelopmental genes are affected by the 3.56 Mb copy number gain of 1p21.3p21.2, which may be considered reciprocal in gene content to the recently recognized 1p21.3 microdeletion syndrome. Both 1p21.3 deletions and the presented duplication display overlapping symptoms, fitting the same disorder category. Contribution of coding and non-coding genes to the phenotype is discussed in the light of cellular and intercellular homeostasis disequilibrium. In line with this the presented 1p21.3p21.2 copy number gain correlated to 1p21.3 microdeletion syndrome verifies the hypothesis of a cumulative effect of the number of deregulated genes - homeostasis disequilibrium leading to overlapping phenotypes between microdeletion and microduplication syndromes.

Although miR-137 appears to be the major player in the 1p21.3p21.2 region, deregulation of the DPYD (dihydropyrimidine dehydrogenase) gene may potentially affect neighboring genes underlying the overlapping symptoms present in both the copy number loss and copy number gain of 1p21. Namely, the all-in approach revealed that DPYD is a complex gene whose expression is epigenetically regulated by long non-coding RNAs (lncRNAs) within the locus. Furthermore, the long interspersed nuclear element-1 (LINE-1) L1MC1 transposon inserted in DPYD intronic transcript 1 (DPYD-IT1) lncRNA with its parasites, TcMAR-Tigger5b and pair of Alu repeats appears to be the “weakest link” within the DPYD gene liable to break. Identification of the precise mechanism through which DPYD is epigenetically regulated, and underlying reasons why exactly the break (FRA1E) happens, will consequently pave the way toward preventing severe toxicity to the antineoplastic drug 5-fluorouracil (5-FU) and development of the causative therapy for the dihydropyrimidine dehydrogenase deficiency.

Retrotransposition of long interspersed nucleotide element-1 is associated with colitis but not tumors in a murine colitic cancer model

Long interspersed element-1 (L1) is a transposable element that can move within the genome, potentially leading to genome diversity and modified gene function. Although L1 activity in somatic cells is normally suppressed through DNA methylation, some L1s are activated in tumors including colorectal carcinoma. However, how L1-retrotransposition (L1-RTP) is involved in gastrointestinal disorders remains to be elucidated. We hypothesized that L1-RTP in somatic cells might contribute to colitis-associated cancer (CAC). To address this, we employed an experimental model of CAC using transgenic L1-reporter mice carrying a human L1-EGFP reporter gene. Mice were subjected to repeated cycles of colitis induced by administration of dextran sodium sulfate (DSS) in drinking water with injection of carcinogen azoxymethane (AOM). L1-RTP levels were measured by a quantitative polymerase chain reaction targeting the newly inserted reporter EGFP in various tissues and cell types, including samples obtained by laser microdissection and cell sorting with flow cytometry. DNA methylation levels of the human L1 promoter were analyzed by bisulfite pyrosequencing. AOM+DSS-treated mice exhibited significantly higher levels of L1-RTP in whole colon tissue during the acute phase of colitis when compared with control naïve mice. L1-RTP levels in whole colon tissue were positively correlated with the histological severity of colitis and degree of neutrophil infiltration into the lamina propria (LP), but not with tumor development in the colon. L1-RTP was enriched in LP mesenchymal cells rather than epithelial cells (ECs), myeloid, or lymphoid cells. DNA methylation levels of the human L1 promoter region showed a negative correlation with L1-RTP levels. L1-RTP was absent from most tumors found in 22-week-old mice. In conclusion, we demonstrated that L1-RTP was induced in the mouse CAC mucosa in accordance with the acute inflammatory response; however, retrotransposition appears not to have direct relevance to colitis-induced cancer initiation.

The evidence for increased L1 activity in the site of human adult brain neurogenesis

Retroelement activity is a common source of polymorphisms in human genome. The mechanism whereby retroelements contribute to the intraindividual genetic heterogeneity by inserting into the DNA of somatic cells is gaining increasing attention. Brain tissues are suspected to accumulate genetic heterogeneity as a result of the retroelements somatic activity. This study aims to expand our understanding of the role retroelements play in generating somatic mosaicism of neural tissues. Whole-genome Alu and L1 profiling of genomic DNA extracted from the cerebellum, frontal cortex, subventricular zone, dentate gyrus, and the myocardium revealed hundreds of somatic insertions in each of the analyzed tissues. Interestingly, the highest concentration of such insertions was detected in the dentate gyrus—the hotspot of adult neurogenesis. Insertions of retroelements and their activity could produce genetically diverse neuronal subsets, which can be involved in hippocampal-dependent learning and memory.

Library Construction for High-Throughput Mobile Element Identification and Genotyping

Mobile genetic elements are discrete DNA elements that can move around and copy themselves in a genome. As a ubiquitous component of the genome, mobile elements contribute to both genetic and epigenetic variation. Therefore, it is important to determine the genome-wide distribution of mobile elements. Here we present a targeted high-throughput sequencing protocol called Mobile Element Scanning (ME-Scan) for genome-wide mobile element detection. We will describe oligonucleotides design, sequencing library construction, and computational analysis for the ME-Scan protocol.

A droplet digital PCR detection method for rare L1 insertions in tumors

Background

The active human mobile element, long interspersed element 1 (L1) currently populates human genomes in excess of 500,000 copies per haploid genome. Through its mobility via a process called target primed reverse transcription (TPRT), L1 mobilization has resulted in over 100 de novo cases of human disease and has recently been associated with various cancer types. Large advances in high-throughput sequencing (HTS) technology have allowed for an increased understanding of the role of L1 in human cancer; however, researchers are still limited by the ability to validate potentially rare L1 insertion events detected by HTS that may occur in only a small fraction of tumor cells. Additionally, HTS detection of rare events varies greatly as a function of read depth, and new tools for de novo element discovery are needed to fill in gaps created by HTS.

Results

We have employed droplet digital PCR (ddPCR) to detect rare L1 loci in mosaic human genomes. Our assay allows for the detection of L1 insertions as rare as one cell in every 10,000.

Conclusions

ddPCR represents a robust method to be used alongside HTS techniques for detecting, validating and quantitating rare L1 insertion events in tumors and other tissues.

Electronic supplementary material

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

T-lex2: genotyping, frequency estimation and re-annotation of transposable elements using single or pooled next-generation sequencing data

Transposable elements (TEs) constitute the most active, diverse and ancient component in a broad range of genomes. Complete understanding of genome function and evolution cannot be achieved without a thorough understanding of TE impact and biology. However, in-depth analysis of TEs still represents a challenge due to the repetitive nature of these genomic entities. In this work, we present a broadly applicable and flexible tool: T-lex2. T-lex2 is the only available software that allows routine, automatic and accurate genotyping of individual TE insertions and estimation of their population frequencies both using individual strain and pooled next-generation sequencing data. Furthermore, T-lex2 also assesses the quality of the calls allowing the identification of miss-annotated TEs and providing the necessary information to re-annotate them. The flexible and customizable design of T-lex2 allows running it in any genome and for any type of TE insertion. Here, we tested the fidelity of T-lex2 using the fly and human genomes. Overall, T-lex2 represents a significant improvement in our ability to analyze the contribution of TEs to genome function and evolution as well as learning about the biology of TEs. T-lex2 is freely available online at http://sourceforge.net/projects/tlex.

Cohesin Rad21 mediates loss of heterozygosity and is upregulated via Wnt promoting transcriptional dysregulation in gastrointestinal tumors

Loss of heterozygosity (LOH) of the adenomatous polyposis coli (APC) gene triggers a series of molecular events leading to intestinal adenomagenesis. Haploinsufficiency of the cohesin Rad21 influences multiple initiating events in colorectal cancer (CRC). We identify Rad21 as a gatekeeper of LOH and a β-catenin target gene and provide evidence that Wnt pathway activation drives RAD21 expression in human CRC. Genome-wide analyses identified Rad21 as a key transcriptional regulator of critical CRC genes and long interspersed element (LINE-1 or L1) retrotransposons. Elevated RAD21 expression tracks with reactivation of L1 expression in human sporadic CRC, implicating cohesin-mediated L1 expression in global genomic instability and gene dysregulation in cancer.

euL1db: the European database of L1HS retrotransposon insertions in humans

Retrotransposons account for almost half of our genome. They are mobile genetics elements—also known as jumping genes—but only the L1HS subfamily of Long Interspersed Nuclear Elements (LINEs) has retained the ability to jump autonomously in modern humans. Their mobilization in germline—but also some somatic tissues—contributes to human genetic diversity and to diseases, such as cancer. Here, we present euL1db, the European database of L1HS retrotransposon insertions in humans (available at http://euL1db.unice.fr). euL1db provides a curated and comprehensive summary of L1HS insertion polymorphisms identified in healthy or pathological human samples and published in peer-reviewed journals. A key feature of euL1db is its sample-wise organization. Hence L1HS insertion polymorphisms are connected to samples, individuals, families and clinical conditions. The current version of euL1db centralizes results obtained in 32 studies. It contains >900 samples, >140 000 sample-wise insertions and almost 9000 distinct merged insertions. euL1db will help understanding the link between L1 retrotransposon insertion polymorphisms and phenotype or disease.

A proposal for the reference-based annotation of de novo transposable element insertions

Understanding the causes and consequences of transposable element (TE) activity in the genomic era requires sophisticated bioinformatics approaches to accurately identify individual insertion sites. Next-generation sequencing technology now makes it possible to rapidly identify new TE insertions using resequencing data, opening up new possibilities to study the nature of TE-induced mutation and the target site preferences of different TE families. While the identification of new TE insertion sites is seemingly a simple task, the mechanisms of transposition present unique challenges for the annotation of de novo transposable element insertions mapped to a reference genome. Here I discuss these challenges and propose a framework for the annotation of de novo TE insertions that accommodates known mechanisms of TE insertion and established coordinate systems for genome annotation.

LINE-1 retrotransposons and let-7 miRNA: partners in the pathogenesis of cancer?

Long interspersed nuclear element-1 (LINE-1 or L1) retrotransposons are insertional mutagens capable of altering the genomic landscape in many ways. Activation of the normally silent LINE-1 retrotransposon is associated with a high level of cancer-associated DNA damage and genomic instability. Studies of LINE-1 have so far focused mainly on changes in gene expression, and our knowledge of its impact on functional non-coding RNAs is in its infancy. However, current evidence suggests that a significant number of human miRNAs originate from retrotransposon sequences. Furthermore, LINE-1 is generally not expressed in normal tissues while its expression is widespread in epithelial cancers. Based on our recent studies, we demonstrate a functional link between aberrant LINE-1 expression and deregulation of let-7 miRNA expression. Since the expression of let-7 is modulated by LINE-1 activity, we discuss possible mechanisms for this effect and how the silencing of LINE-1 activation could provide new therapeutic options for cancer treatment. Based on the deep sequencing of small RNAs in parallel with gene expression profiling in breast cancer cells, we have identified potential pathways linking L1 activity to let-7 processing and maturation and ultimately to the control of stemness in human cancer cells.

Deep sequencing reveals low incidence of endogenous LINE-1 retrotransposition in human induced pluripotent stem cells

Long interspersed element-1 (LINE-1 or L1) retrotransposition induces insertional mutations that can result in diseases. It was recently shown that the copy number of L1 and other retroelements is stable in induced pluripotent stem cells (iPSCs). However, by using an engineered reporter construct over-expressing L1, another study suggests that reprogramming activates L1 mobility in iPSCs. Given the potential of human iPSCs in therapeutic applications, it is important to clarify whether these cells harbor somatic insertions resulting from endogenous L1 retrotransposition. Here, we verified L1 expression during and after reprogramming as well as potential somatic insertions driven by the most active human endogenous L1 subfamily (L1Hs). Our results indicate that L1 over-expression is initiated during the reprogramming process and is subsequently sustained in isolated clones. To detect potential somatic insertions in iPSCs caused by L1Hs retotransposition, we used a novel sequencing strategy. As opposed to conventional sequencing direction, we sequenced from the 3′ end of L1Hs to the genomic DNA, thus enabling the direct detection of the polyA tail signature of retrotransposition for verification of true insertions. Deep coverage sequencing thus allowed us to detect seven potential somatic insertions with low read counts from two iPSC clones. Negative PCR amplification in parental cells, presence of a polyA tail and absence from seven L1 germline insertion databases highly suggested true somatic insertions in iPSCs. Furthermore, these insertions could not be detected in iPSCs by PCR, likely due to low abundance. We conclude that L1Hs retrotransposes at low levels in iPSCs and therefore warrants careful analyses for genotoxic effects.

SIRT6 represses LINE1 retrotransposons by ribosylating KAP1 but this repression fails with stress and age

L1 retrotransposons are an abundant class of transposable elements which pose a threat to genome stability and may play a role in age-related pathologies such as cancer. Recent evidence indicates that L1s become more active in somatic tissues during the course of aging; the mechanisms underlying this phenomenon remain unknown, however. Here we report that the longevity regulating protein, SIRT6, is a powerful repressor of L1 activity. Specifically, SIRT6 binds to the 5′UTR of L1 loci, where it mono-ADP ribosylates the nuclear corepressor protein, KAP1, and facilitates KAP1 interaction with the heterochromatin factor, HP1α, thereby contributing to the packaging of L1 elements into transcriptionally repressive heterochromatin. During the course of aging, and also in response to DNA damage, however, we find that SIRT6 is depleted from L1 loci, allowing for the activation of these previously silenced retroelements.

Repetitive elements dynamics in cell identity programming, maintenance and disease

The days of 'junk DNA' seem to be over. The rapid progress of genomics technologies has been unveiling unexpected mechanisms by which repetitive DNA and in particular transposable elements (TEs) have evolved, becoming key issues in understanding genome structure and function. Indeed, rather than 'parasites', recent findings strongly suggest that TEs may have a positive function by contributing to tissue specific transcriptional programs, in particular as enhancer-like elements and/or modules for regulation of higher order chromatin structure. Further, it appears that during development and aging genomes experience several waves of TEs activation, and this contributes to individual genome shaping during lifetime. Interestingly, TEs activity is major target of epigenomic regulation. These findings are shedding new light on the genome-phenotype relationship and set the premises to help to explain complex disease manifestation, as consequence of TEs activity deregulation.

Epigenetic diversity in hematopoietic neoplasms

Tumor cell populations display a remarkable extent of variability in non-genetic characteristics such as DNA methylation, histone modification patterns, and differentiation levels of individual cells. It remains to be elucidated whether non-genetic heterogeneity is simply a byproduct of tumor evolution or instead a manifestation of a higher-order tissue organization that is maintained within the neoplasm to establish a differentiation hierarchy, a favorable microenvironment, or a buffer against changing selection pressures during tumorigenesis. Here, we review recent findings on epigenetic diversity, particularly heterogeneity in DNA methylation patterns in hematologic malignancies. We also address the implications of epigenetic heterogeneity for the clonal evolution of tumors and discuss its effects on gene expression and other genome functions in cancer.

Tangram: a comprehensive toolbox for mobile element insertion detection

Background

Mobile elements (MEs) constitute greater than 50% of the human genome as a result of repeated insertion events during human genome evolution. Although most of these elements are now fixed in the population, some MEs, including ALU, L1, SVA and HERV-K elements, are still actively duplicating. Mobile element insertions (MEIs) have been associated with human genetic disorders, including Crohn’s disease, hemophilia, and various types of cancer, motivating the need for accurate MEI detection methods. To comprehensively identify and accurately characterize these variants in whole genome next-generation sequencing (NGS) data, a computationally efficient detection and genotyping method is required. Current computational tools are unable to call MEI polymorphisms with sufficiently high sensitivity and specificity, or call individual genotypes with sufficiently high accuracy.

Results

Here we report Tangram, a computationally efficient MEI detection program that integrates read-pair (RP) and split-read (SR) mapping signals to detect MEI events. By utilizing SR mapping in its primary detection module, a feature unique to this software, Tangram is able to pinpoint MEI breakpoints with single-nucleotide precision. To understand the role of MEI events in disease, it is essential to produce accurate individual genotypes in clinical samples. Tangram is able to determine sample genotypes with very high accuracy. Using simulations and experimental datasets, we demonstrate that Tangram has superior sensitivity, specificity, breakpoint resolution and genotyping accuracy, when compared to other, recently developed MEI detection methods.

Conclusions

Tangram serves as the primary MEI detection tool in the 1000 Genomes Project, and is implemented as a highly portable, memory-efficient, easy-to-use C++ computer program, built under an open-source development model.

Increased expression and copy number amplification of LINE-1 and SINE B1 retrotransposable elements in murine mammary carcinoma progression

In higher eukaryotic genomes, Long Interspersed Nuclear Element 1 (LINE-1) retrotransposons and endogenous retroviruses represent large families of repeated elements encoding reverse transcriptase (RT) proteins. Short Interspersed Nuclear Element B1 (SINE B1) retrotrasposons do not encode RT, but use LINE-1-derived RT for their retrotransposition. We previously showed that many cancer types have an abundant endogenous RT activity. Inhibition of that activity, by either RNA interference-dependent silencing of active LINE-1 elements or by RT inhibitory drugs, reduced proliferation and promoted differentiation in cancer cells, indicating that LINE-1-encoded RT is required for tumor progression. Using MMTV-PyVT transgenic mice as a well-defined model of breast cancer progression, we now report that both LINE-1 and SINE B1 retrotransposons are up-regulated at a very early stage of tumorigenesis; LINE-1-encoded RT product and enzymatic activity were detected in tumor tissues as early as stage 1, preceding the widespread appearance of histological alterations and specific cancer markers, and further increased in later progression stages, while neither was present in non-pathological breast tissues. Importantly, both LINE-1 and SINE B1 retrotransposon families undergo copy number amplification during tumor progression. These findings therefore indicate that RT activity is distinctive of breast cancer cells and that, furthermore, LINE-1 and SINE B1 undergo copy number amplification during cancer progression.

Dynamic Alu methylation during normal development, aging, and tumorigenesis

DNA methylation primarily occurs on CpG dinucleotides and plays an important role in transcriptional regulations during tissue development and cell differentiation. Over 25% of CpG dinucleotides in the human genome reside within Alu elements, the most abundant human repeats. The methylation of Alu elements is an important mechanism to suppress Alu transcription and subsequent retrotransposition. Decades of studies revealed that Alu methylation is highly dynamic during early development and aging. Recently, many environmental factors were shown to have a great impact on Alu methylation. In addition, aberrant Alu methylation has been documented to be an early event in many tumors and Alu methylation levels have been associated with tumor aggressiveness. The assessment of the Alu methylation has become an important approach for early diagnosis and/or prognosis of cancer. This review focuses on the dynamic Alu methylation during development, aging, and tumor genesis. The cause and consequence of Alu methylation changes will be discussed.

Potential for genomic instability associated with retrotranspositionally-incompetent L1 loci

Expression of the L1 retrotransposon can damage the genome through insertional mutagenesis and the generation of DNA double-strand breaks (DSBs). The majority of L1 loci in the human genome are 5'-truncated and therefore incapable of retrotransposition. While thousands of full-length L1 loci remain, most are retrotranspositionally-incompetent due to inactivating mutations. However, mutations leading to premature stop codons within the L1 ORF2 sequence may yield truncated proteins that retain a functional endonuclease domain. We demonstrate that some truncated ORF2 proteins cause varying levels of toxicity and DNA damage when chronically overexpressed in mammalian cells. Furthermore, transfection of some ORF2 constructs containing premature stop codons supported low levels of Alu retrotransposition, demonstrating the potential for select retrotranspositionally-incompetent L1 loci to generate genomic instability. This result suggests yet another plausible explanation for the relative success of Alu elements in populating the human genome. Our data suggest that a subset of retrotranspositionally-incompetent L1s, previously considered to be harmless to genomic integrity, may have the potential to cause chronic DNA damage by introducing DSBs and mobilizing Alu. These results imply that the number of known L1 loci in the human genome that potentially threaten its stability may not be limited to the retrotranspositionally active loci.

Distinct isoform of FABP7 revealed by screening for retroelement-activated genes in diffuse large B-cell lymphoma

Remnants of ancient transposable elements (TEs) are abundant in mammalian genomes. These sequences harbor multiple regulatory motifs and hence are capable of influencing expression of host genes. In response to environmental changes, TEs are known to be released from epigenetic repression and to become transcriptionally active. Such activation could also lead to lineage-inappropriate activation of oncogenes, as one study described in Hodgkin lymphoma. However, little further evidence for this mechanism in other cancers has been reported. Here, we reanalyzed whole transcriptome data from a large cohort of patients with diffuse large B-cell lymphoma (DLBCL) compared with normal B-cell centroblasts to detect genes ectopically expressed through activation of TE promoters. We have identified 98 such TE-gene chimeric transcripts that were exclusively expressed in primary DLBCL cases and confirmed several in DLBCL-derived cell lines. We further characterized a TE-gene chimeric transcript involving a fatty acid-binding protein gene (LTR2-FABP7), normally expressed in brain, that was ectopically expressed in a subset of DLBCL patients through the use of an endogenous retroviral LTR promoter of the LTR2 family. The LTR2-FABP7 chimeric transcript encodes a novel chimeric isoform of the protein with characteristics distinct from native FABP7. In vitro studies reveal a dependency for DLBCL cell line proliferation and growth on LTR2-FABP7 chimeric protein expression. Taken together, these data demonstrate the significance of TEs as regulators of aberrant gene expression in cancer and suggest that LTR2-FABP7 may contribute to the pathogenesis of DLBCL in a subgroup of patients.

Mobile DNA in cancer. Extensive transduction of nonrepetitive DNA mediated by L1 retrotransposition in cancer genomes

Long interspersed nuclear element-1 (L1) retrotransposons are mobile repetitive elements that are abundant in the human genome. L1 elements propagate through RNA intermediates. In the germ line, neighboring, nonrepetitive sequences are occasionally mobilized by the L1 machinery, a process called 3' transduction. Because 3' transductions are potentially mutagenic, we explored the extent to which they occur somatically during tumorigenesis. Studying cancer genomes from 244 patients, we found that tumors from 53% of the patients had somatic retrotranspositions, of which 24% were 3' transductions. Fingerprinting of donor L1s revealed that a handful of source L1 elements in a tumor can spawn from tens to hundreds of 3' transductions, which can themselves seed further retrotranspositions. The activity of individual L1 elements fluctuated during tumor evolution and correlated with L1 promoter hypomethylation. The 3' transductions disseminated genes, exons, and regulatory elements to new locations, most often to heterochromatic regions of the genome.

Retrotransposition of long interspersed element 1 induced by methamphetamine or cocaine

Long interspersed element 1 (L1) is a retroelement constituting ∼17% of the human genome. A single human cell has 80–100 copies of L1 capable of retrotransposition (L1-RTP), ∼10% of which are “hot L1” copies, meaning they are primed for “jumping” within the genome. Recent studies demonstrated induction of L1 activity by drugs of abuse or low molecular weight compounds, but little is known about the underlying mechanism. The aim of this study was to identify the mechanism and effects of methamphetamine (METH) and cocaine on L1-RTP. Our results revealed that METH and cocaine induced L1-RTP in neuronal cell lines. This effect was found to be reverse transcriptase-dependent. However, METH and cocaine did not induce double-strand breaks. RNA interference experiments combined with add-back of siRNA-resistant cDNAs revealed that the induction of L1-RTP by METH or cocaine depends on the activation of cAMP response element-binding protein (CREB). METH or cocaine recruited the L1-encoded open reading frame 1 (ORF1) to chromatin in a CREB-dependent manner. These data suggest that the cellular cascades underlying METH- and cocaine-induced L1-RTP are different from those behind L1-RTP triggered by DNA damage; CREB is involved in drug-induced L1-RTP. L1-RTP caused by drugs of abuse is a novel type of genomic instability, and analysis of this phenomenon might be a novel approach to studying substance-use disorders.

Retrotransposons in pluripotent cells: Impact and new roles in cellular plasticity

Transposable Elements are pieces of DNA able to mobilize from one location to another within genomes. Although they constitute more than 50% of the human genome, they have been classified as selfish DNA, with the only mission to spread within genomes and generate more copies of themselves that will ensure their presence over generations. Despite their remarkable prevalence, only a minor group of transposable elements remain active in the human genome and can sporadically be associated with the generation of a genetic disorder due to their ongoing mobility. Most of the transposable elements identified in the human genome corresponded to fixed insertions that no longer move in genomes. As selfish DNA, transposable element insertions accumulate in cell types where genetic information can be passed to the next generation. Indeed, work from different laboratories has demonstrated that the main heritable load of TE accumulation in humans occurs during early embryogenesis. Thus, active transposable elements have a clear impact on our pluripotent genome. However, recent findings suggest that the main proportion of fixed non-mobile transposable elements might also have emerging roles in cellular plasticity. In this concise review, we provide an overview of the impact of currently active transposable elements in our pluripotent genome and further discuss new roles of transposable elements (active or not) in regulating pluripotency. This article is part of a Special Issue entitled: Stress as a fundamental theme in cell plasticity.

[Retrotransposons: selfish DNA or active epigenetic players in somatic cells?]

Transposable elements (TE) represent around 40% of the human genome. They are endogenous mobile DNA sequences able to jump and duplicate in the host genome. TE have long been considered as "junk" DNA but are now believed to be important regulators of gene expression by participating to the establishment of the DNA methylation profile. Recent advances in genome sequencing reveals a higher transposition frequency and TE driven gene expression in somatic cells than previously thought. As TE propagation is deleterious and may be involved in oncogenic mechanisms, host cells have developed silencing mechanisms mainly described in germinal and embryonic cells. However, somatic cells are also proned to TE transposition and use specific mechanisms involving tumor suppressor proteins including p53, Rb and PLZF. These transcription factors specifically target genomic retrotransposon sequences, histone deacetylase and DNA methylase activities, inducing epigenetic modifications related to gene silencing. Thus, these transcription factors negatively regulate TE expression by the formation of DNA methylation profil in somatic cells possibly associated with oncogenic mechanisms.

Mobile DNA elements in the generation of diversity and complexity in the brain

Mobile elements are DNA sequences that can change their position (retrotranspose) within the genome. Although its biological function is largely unappreciated, DNA derived from mobile elements comprises nearly half of the human genome. It has long been thought that neuronal genomes are invariable; however, recent studies have demonstrated that mobile elements actively retrotranspose during neurogenesis, thereby creating genomic diversity between neurons. In addition, mounting data demonstrate that mobile elements are misregulated in certain neurological disorders, including Rett syndrome and schizophrenia.

The enemy within: endogenous retroelements and autoimmune disease

Inappropriate or chronic detection of self nucleic acids by the innate immune system underlies many human autoimmune diseases. We discuss here an unexpected source of endogenous immunostimulatory nucleic acids: the reverse-transcribed cDNA of endogenous retroelements. The interplay between innate immune sensing and clearance of retroelement cDNA has important implications for the understanding of immune responses to infectious retroviruses such as human immunodeficiency virus (HIV). Furthermore, the detection of cDNA by the innate immune system reveals an evolutionary tradeoff: selection for a vigorous, sensitive response to infectious retroviruses may predispose the inappropriate detection of endogenous retroelements. We propose that this tradeoff has placed unique constraints on the sensitivity of the DNA-activated antiviral response, with implications for the interactions of DNA viruses and retroviruses with their hosts. Finally, we discuss how better understanding of the intersection of retroelement biology and innate immunity can guide the way to novel therapies for specific autoimmune diseases.

Regulation of L1 expression and retrotransposition by melatonin and its receptor: implications for cancer risk associated with light exposure at night

Expression of long interspersed element-1 (L1) is upregulated in many human malignancies. L1 can introduce genomic instability via insertional mutagenesis and DNA double-strand breaks, both of which may promote cancer. Light exposure at night, a recently recognized carcinogen, is associated with an increased risk of cancer in shift workers. We report that melatonin receptor 1 inhibits mobilization of L1 in cultured cells through downregulation of L1 mRNA and ORF1 protein. The addition of melatonin receptor antagonists abolishes the MT1 effect on retrotransposition in a dose-dependent manner. Furthermore, melatonin-rich, but not melatonin-poor, human blood collected at different times during the circadian cycle suppresses endogenous L1 mRNA during in situ perfusion of tissue-isolated xenografts of human cancer. Supplementation of human blood with exogenous melatonin or melatonin receptor antagonist during the in situ perfusion establishes a receptor-mediated action of melatonin on L1 expression. Combined tissue culture and in vivo data support that environmental light exposure of the host regulates expression of L1 elements in tumors. Our data imply that light-induced suppression of melatonin production in shift workers may increase L1-induced genomic instability in their genomes and suggest a possible connection between L1 activity and increased incidence of cancer associated with circadian disruption.

Somatic retrotransposition in human cancer revealed by whole-genome and exome sequencing

Retrotransposons constitute a major source of genetic variation, and somatic retrotransposon insertions have been reported in cancer. Here, we applied TranspoSeq, a computational framework that identifies retrotransposon insertions from sequencing data, to whole genomes from 200 tumor/normal pairs across 11 tumor types as part of The Cancer Genome Atlas (TCGA) Pan-Cancer Project. In addition to novel germline polymorphisms, we find 810 somatic retrotransposon insertions primarily in lung squamous, head and neck, colorectal, and endometrial carcinomas. Many somatic retrotransposon insertions occur in known cancer genes. We find that high somatic retrotransposition rates in tumors are associated with high rates of genomic rearrangement and somatic mutation. Finally, we developed TranspoSeq-Exome to interrogate an additional 767 tumor samples with hybrid-capture exome data and discovered 35 novel somatic retrotransposon insertions into exonic regions, including an insertion into an exon of the PTEN tumor suppressor gene. The results of this large-scale, comprehensive analysis of retrotransposon movement across tumor types suggest that somatic retrotransposon insertions may represent an important class of structural variation in cancer.

Efficient screening of long terminal repeat retrotransposons that show high insertion polymorphism via high-throughput sequencing of the primer binding site

Retrotransposons have been used frequently for the development of molecular markers by using their insertion polymorphisms among cultivars, because multiple copies of these elements are dispersed throughout the genome and inserted copies are inherited genetically. Although a large number of long terminal repeat (LTR) retrotransposon families exist in the higher eukaryotic genomes, the identification of families that show high insertion polymorphism has been challenging. Here, we performed an efficient screening of these retrotransposon families using an Illumina HiSeq2000 sequencing platform with comprehensive LTR library construction based on the primer binding site (PBS), which is located adjacent to the 5' LTR and has a motif that is universal and conserved among LTR retrotransposon families. The paired-end sequencing library of the fragments containing a large number of LTR sequences and their insertion sites was sequenced for seven strawberry (Fragaria × ananassa Duchesne) cultivars and one diploid wild species (Fragaria vesca L.). Among them, we screened 24 families with a "unique" insertion site that appeared only in one cultivar and not in any others, assuming that this type of insertion should have occurred quite recently. Finally, we confirmed experimentally the selected LTR families showed high insertion polymorphisms among closely related cultivars.

Exposure of hepatocellular carcinoma cells to low-level As₂O₃‎ causes an extra toxicity pathway via L1 retrotransposition ‎induction

Various mechanisms have been proposed for toxicity and carcinogenesis pattern of arsenic, a naturally occurring metalloid. The extent to which the long interspersed element-1 (LINE-1) retrotransposon, an ubiquitous retroelement with autonomous mobility, can be influenced upon exposure to low-level arsenic remains to be elucidated. The aim of this study was to evaluate the possible effect of low-level As2O3 on L1 retrotransposition alteration in human hepatocellular carcinoma cells (HepG2). L1 retrotransposition in HepG2 cells was performed by the in vitro retrotransposition assay using an EGFP-tagged L1RP. Following determination of non-cytotoxic concentrations of arsenic by a MTT assay, the cells were transfected with pL1RP-EGFP and then exposed to 0.25, 0.50 and 0.75 μM of As2O3. The amount of EGFP and its copy number in retrotransposed cells were evaluated by FACS and qPCR analysis in treated vs. control cells, respectively. Significant increase in retrotransposition frequency was found after 12 days exposure to 0.50 and 0.75 μM of As2O3 by FACS analysis (P<0.05). Obtained results were further confirmed by real time PCR, which showed significant induction of retrotransposition in all mentioned concentrations. Our findings indicate that low-level long-term As2O3 exposure may pave activation of L1 retrotransposon.

Mobile elements and viral integrations prompt considerations for bacterial DNA integration as a novel carcinogen

Insertional mutagenesis has been repeatedly demonstrated in cancer genomes and has a role in oncogenesis. Mobile genetic elements can induce cancer development by random insertion into cancer related genes or by inducing translocations. L1s are typically implicated in cancers of an epithelial cell origin, while Alu elements have been implicated in leukemia as well as epithelial cell cancers. Likewise, viral infections have a significant role in cancer development predominantly through integration into the human genome and mutating or deregulating cancer related genes. Human papilloma virus is the best-known example of viral integrations contributing to carcinogenesis. However, hepatitis B virus, Epstein-Barr virus, and Merkel cell polyomavirus also integrate into the human genome and disrupt cancer related genes. Thus far, the role of microbes in cancer has primarily been attributed to mutations induced through chronic inflammation or toxins, as is the case with Helicobacter pylori and enterotoxigenic Bacteroides fragilis. We hypothesize that like mobile elements and viral DNA, bacterial and parasitic DNA may also integrate into the human somatic genome and be oncogenic. Until recently it was believed that bacterial DNA could not integrate into the human genome, but new evidence demonstrates that bacterial insertional mutagenesis may occur in cancer cells. Although this work does not show causation between bacterial insertions and cancer, it prompts more research in this area. Promising new sequencing technologies may reduce the risk of artifactual chimeric sequences, thus diminishing some of the challenges of identifying novel insertions in the somatic human genome.

Unfixed endogenous retroviral insertions in the human population

One lineage of human endogenous retroviruses (HERVs), HERV-K(HML2), is upregulated in many cancers, some autoimmune/inflammatory diseases, and HIV-infected cells. Despite 3 decades of research, it is not known if these viruses play a causal role in disease, and there has been recent interest in whether they can be used as immunotherapy targets. Resolution of both these questions will be helped by an ability to distinguish between the effects of different integrated copies of the virus (loci). Research so far has concentrated on the 20 or so recently integrated loci that, with one exception, are in the human reference genome sequence. However, this viral lineage has been copying in the human population within the last million years, so some loci will inevitably be present in the human population but absent from the reference sequence. We therefore performed the first detailed search for such loci by mining whole-genome sequences generated by next-generation sequencing. We found a total of 17 loci, and the frequency of their presence ranged from only 2 of the 358 individuals examined to over 95% of them. On average, each individual had six loci that are not in the human reference genome sequence. Comparing the number of loci that we found to an expectation derived from a neutral population genetic model suggests that the lineage was copying until at least ∼250,000 years ago.

IMPORTANCE About 5% of the human genome sequence is composed of the remains of retroviruses that over millions of years have integrated into the chromosomes of egg and/or sperm precursor cells. There are indications that protein expression of these viruses is higher in some diseases, and we need to know (i) whether these viruses have a role in causing disease and (ii) whether they can be used as immunotherapy targets in some of them. Answering both questions requires a better understanding of how individuals differ in the viruses that they carry. We carried out the first careful search for new viruses in some of the many human genome sequences that are now available thanks to advances in sequencing technology. We also compared the number that we found to a theoretical expectation to see if it is likely that these viruses are still replicating in the human population today.