Antisense promoter of human L1 retrotransposon drives transcription of adjacent cellular genes

Abnormal activation of genomic LINE1 elements caused by DNA demethylation contributes to lncRNA CASC9 overexpression in esophageal squamous cell carcinoma

Long noncoding RNA (lncRNA) cancer susceptibility 9 (CASC9) has been found to be overexpressed and functions as an oncogene in many cancer types. We investigated the molecular mechanism underlying CASC9 overexpression in esophageal squamous cell carcinoma (ESCC). Transcripts containing exons 2 and 6 and exons 4 and 6 showed the highest CASC9 expression levels in ESCC, no transcripts were detected in the normal esophageal epithelial Het1A cell line. The Long Interspersed Nuclear Element-1 (LINE1 or L1) element in the genome was found to participate in the evolution of lncRNA CASC9, the antisense promoter (ASP) of L1 provides the cis-regulatory elements necessary for CASC9 activation, and the antisense chain of L1 participates in the formation of exons of CASC9. The activation of the antisense promoter was due to the aberrant hypomethylation of L1 elements. An active enhancer element was identified in the downstream region of CASC9 gene by ChIP-seq and ChIP-qPCR. The interaction between ASP and the enhancer elements was confirmed by chromosome conformation capture (3C). Thus, our results suggest that the L1 ASP activation due to aberrant hypomethylation and downstream enhancer interaction plays a key role in the overexpression of lncRNA CASC9 in ESCC.

Targeting transposable elements in cancer: developments and opportunities

Transposable elements (TEs), comprising nearly 50% of the human genome, have transitioned from being perceived as "genomic junk" to key players in cancer progression. Contemporary research links TE regulatory disruptions with cancer development, underscoring their therapeutic potential. Advances in long-read sequencing, computational analytics, single-cell sequencing, proteomics, and CRISPR-Cas9 technologies have enriched our understanding of TEs' clinical implications, notably their impact on genome architecture, gene regulation, and evolutionary processes. In cancer, TEs, including long interspersed element-1 (LINE-1), Alus, and long terminal repeat (LTR) elements, demonstrate altered patterns, influencing both tumorigenic and tumor-suppressive mechanisms. TE-derived nucleic acids and tumor antigens play critical roles in tumor immunity, bridging innate and adaptive responses. Given their central role in oncology, TE-targeted therapies, particularly through reverse transcriptase inhibitors and epigenetic modulators, represent a novel avenue in cancer treatment. Combining these TE-focused strategies with existing chemotherapy or immunotherapy regimens could enhance efficacy and offer a new dimension in cancer treatment. This review delves into recent TE detection advancements, explores their multifaceted roles in tumorigenesis and immune regulation, discusses emerging diagnostic and therapeutic approaches centered on TEs, and anticipates future directions in cancer research.

RepEnTools: an automated repeat enrichment analysis package for ChIP-seq data reveals hUHRF1 Tandem-Tudor domain enrichment in young repeats

BACKGROUND

Repeat elements (REs) play important roles for cell function in health and disease. However, RE enrichment analysis in short-read high-throughput sequencing (HTS) data, such as ChIP-seq, is a challenging task.

RESULTS

Here, we present RepEnTools, a software package for genome-wide RE enrichment analysis of ChIP-seq and similar chromatin pulldown experiments. Our analysis package bundles together various software with carefully chosen and validated settings to provide a complete solution for RE analysis, starting from raw input files to tabular and graphical outputs. RepEnTools implementations are easily accessible even with minimal IT skills (Galaxy/UNIX). To demonstrate the performance of RepEnTools, we analysed chromatin pulldown data by the human UHRF1 TTD protein domain and discovered enrichment of TTD binding on young primate and hominid specific polymorphic repeats (SVA, L1PA1/L1HS) overlapping known enhancers and decorated with H3K4me1-K9me2/3 modifications. We corroborated these new bioinformatic findings with experimental data by qPCR assays using newly developed primate and hominid specific qPCR assays which complement similar research tools. Finally, we analysed mouse UHRF1 ChIP-seq data with RepEnTools and showed that the endogenous mUHRF1 protein colocalizes with H3K4me1-H3K9me3 on promoters of REs which were silenced by UHRF1. These new data suggest a functional role for UHRF1 in silencing of REs that is mediated by TTD binding to the H3K4me1-K9me3 double mark and conserved in two mammalian species.

CONCLUSIONS

RepEnTools improves the previously available programmes for RE enrichment analysis in chromatin pulldown studies by leveraging new tools, enhancing accessibility and adding some key functions. RepEnTools can analyse RE enrichment rapidly, efficiently, and accurately, providing the community with an up-to-date, reliable and accessible tool for this important type of analysis.

Locus-level L1 DNA methylation profiling reveals the epigenetic and transcriptional interplay between L1s and their integration sites

Long interspersed element 1 (L1) retrotransposons are implicated in human disease and evolution. Their global activity is repressed by DNA methylation, but deciphering the regulation of individual copies has been challenging. Here, we combine short- and long-read sequencing to unveil L1 methylation heterogeneity across cell types, families, and individual loci and elucidate key principles involved. We find that the youngest primate L1 families are specifically hypomethylated in pluripotent stem cells and the placenta but not in most tumors. Locally, intronic L1 methylation is intimately associated with gene transcription. Conversely, the L1 methylation state can propagate to the proximal region up to 300 bp. This phenomenon is accompanied by the binding of specific transcription factors, which drive the expression of L1 and chimeric transcripts. Finally, L1 hypomethylation alone is typically insufficient to trigger L1 expression due to redundant silencing pathways. Our results illuminate the epigenetic and transcriptional interplay between retrotransposons and their host genome.

LINE-1 retrotransposition and its deregulation in cancers: implications for therapeutic opportunities

Long interspersed element 1 (LINE-1) is the only protein-coding transposon that is active in humans. LINE-1 propagates in the genome using RNA intermediates via retrotransposition. This activity has resulted in LINE-1 sequences occupying approximately one-fifth of our genome. Although most copies of LINE-1 are immobile, ∼100 copies are retrotransposition-competent. Retrotransposition is normally limited via epigenetic silencing, DNA repair, and other host defense mechanisms. In contrast, LINE-1 overexpression and retrotransposition are hallmarks of cancers. Here, we review mechanisms of LINE-1 regulation and how LINE-1 may promote genetic heterogeneity in tumors. Finally, we discuss therapeutic strategies to exploit LINE-1 biology in cancers.

An intronic LINE-1 regulates IFNAR1 expression in human immune cells

BACKGROUND

Despite their origins as selfish parasitic sequences, some transposons in the human genome have been co-opted to serve as regulatory elements, contributing to the evolution of transcriptional networks. Most well-characterized examples of transposon-derived regulatory elements derive from endogenous retroviruses (ERVs), due to the intrinsic regulatory activity of proviral long terminal repeat regions. However, one subclass of transposable elements, the Long Interspersed Nuclear Elements (LINEs), have been largely overlooked in the search for functional regulatory transposons, and considered to be broadly epigenetically repressed.

RESULTS

We examined the chromatin state of LINEs by analyzing epigenomic data from human immune cells. Many LINEs are marked by the repressive H3K9me3 modification, but a subset exhibits evidence of enhancer activity in human immune cells despite also showing evidence of epigenetic repression. We hypothesized that these competing forces of repressive and activating epigenetic marks might lead to inducible enhancer activity. We investigated a specific L1M2a element located within the first intron of Interferon Alpha/Beta Receptor 1 (IFNAR1). This element shows epigenetic signatures of B cell-specific enhancer activity, despite being repressed by the Human Silencing Hub (HUSH) complex. CRISPR deletion of the element in B lymphoblastoid cells revealed that the element acts as an enhancer that regulates both steady state and interferon-inducible expression of IFNAR1.

CONCLUSIONS

Our study experimentally demonstrates that an L1M2a element was co-opted to function as an interferon-inducible enhancer of IFNAR1, creating a feedback loop wherein IFNAR1 is transcriptionally upregulated by interferon signaling. This finding suggests that other LINEs may exhibit cryptic cell type-specific or context-dependent enhancer activity. LINEs have received less attention than ERVs in the effort to understand the contribution of transposons to the regulatory landscape of cellular genomes, but these are likely important, lineage-specific players in the rapid evolution of immune system regulatory networks and deserve further study.

Roles of Human Endogenous Retroviruses and Endogenous Virus-Like Elements in Cancer Development and Innate Immunity

Human endogenous retroviruses (HERVs) are remnants of ancient retroviral infections in the host genome. Although mutations and silencing mechanisms impair their original role in viral replication, HERVs are believed to play roles in various biological processes. Long interspersed nuclear elements (LINEs) are non-LTR retrotransposons that have a lifecycle resembling that of retroviruses. Although LINE expression is typically silenced in somatic cells, it also contributes to various biological processes. The aberrant expression of HERVs and LINEs is closely associated with the development of cancer and/or immunological diseases, suggesting that they are integrated into various pathways related to the diseases. HERVs/LINEs control gene expression depending on the context as promoter/enhancer elements. Some RNAs and proteins derived from HERVs/LINEs have oncogenic potential, whereas others stimulate innate immunity. Non-retroviral endogenous viral elements (nrEVEs) are a novel type of virus-like element in the genome. nrEVEs may also be involved in host immunity. This article provides a current understanding of how these elements impact cellular physiology in cancer development and innate immunity, and provides perspectives for future studies.

LINE-1 retrotransposons drive human neuronal transcriptome complexity and functional diversification

The genetic mechanisms underlying the expansion in size and complexity of the human brain remain poorly understood. Long interspersed nuclear element-1 (L1) retrotransposons are a source of divergent genetic information in hominoid genomes, but their importance in physiological functions and their contribution to human brain evolution are largely unknown. Using multiomics profiling, we here demonstrate that L1 promoters are dynamically active in the developing and the adult human brain. L1s generate hundreds of developmentally regulated and cell type-specific transcripts, many that are co-opted as chimeric transcripts or regulatory RNAs. One L1-derived long noncoding RNA, LINC01876, is a human-specific transcript expressed exclusively during brain development. CRISPR interference silencing of LINC01876 results in reduced size of cerebral organoids and premature differentiation of neural progenitors, implicating L1s in human-specific developmental processes. In summary, our results demonstrate that L1-derived transcripts provide a previously undescribed layer of primate- and human-specific transcriptome complexity that contributes to the functional diversification of the human brain.

Molecular Mechanisms of Alu and LINE-1 Interspersed Repetitive Sequences Reveal Diseases of Visual System Dysfunction

BACKGROUND

Short interspersed nuclear elements (SINEs) and long interspersed nuclear elements (LINE-1s) are the abundant and well-characterized repetitive elements in the human genome.

METHODS

For this review, all relevant original research studies were assessed by searching electronic databases, including PubMed, Google Scholar, and Web of Science, by using relevant keywords. Accumulating evidence indicates that the disorder of gene expression regulated by these repetitive sequences is one of the causes of the diseases of visual system dysfunction, including retinal degenerations, glaucoma, retinitis punctata albescens, retinitis pigmentosa, geographic atrophy, and age-related macular degeneration, suggesting that SINEs and LINE-1s may have great potential implications in ophthalmology.

RESULTS

Alu elements belonging to the SINEs are present in more than one million copies, comprising 10% of the human genome.

CONCLUSION

This study offers recent advances in Alu and LINE-1 mechanisms in the development of eye diseases. The current study could advance our knowledge of the roles of SINEs and LINE-1s in the developing process of eye diseases, suggesting new diagnostic biomarkers, therapeutic strategies, and significant points for future studies.

Retrotransposons and Telomeres

Transposable elements (TEs) comprise a significant part of eukaryotic genomes being a major source of genome instability and mutagenesis. Cellular defense systems suppress the TE expansion at all stages of their life cycle. Piwi proteins and Piwi-interacting RNAs (piRNAs) are key elements of the anti-transposon defense system, which control TE activity in metazoan gonads preventing inheritable transpositions and developmental defects. In this review, we discuss various regulatory mechanisms by which small RNAs combat TE activity. However, active transposons persist, suggesting these powerful anti-transposon defense mechanisms have a limited capacity. A growing body of evidence suggests that increased TE activity coincides with genome reprogramming and telomere lengthening in different species. In the Drosophila fruit fly, whose telomeres consist only of retrotransposons, a piRNA-mediated mechanism is required for telomere maintenance and their length control. Therefore, the efficacy of protective mechanisms must be finely balanced in order not only to suppress the activity of transposons, but also to maintain the proper length and stability of telomeres. Structural and functional relationship between the telomere homeostasis and LINE1 retrotransposon in human cells indicates a close link between selfish TEs and the vital structure of the genome, telomere. This relationship, which permits the retention of active TEs in the genome, is reportedly a legacy of the retrotransposon origin of telomeres. The maintenance of telomeres and the execution of other crucial roles that TEs acquired during the process of their domestication in the genome serve as a type of payment for such a "service."

Identification of epigenetically active L1 promoters in the human brain and their relationship with psychiatric disorders

Long interspersed nuclear element-1 (LINE-1, L1) affects the transcriptome landscape in multiple ways. Promoter activity within its 5'UTR plays a critical role in regulating diverse L1 activities. However, the epigenetic status of L1 promoters in adult brain cells and their relationship with psychiatric disorders remain poorly understood. Here, we examined DNA methylation and hydroxymethylation of the full-length L1s in neurons and nonneurons and identified "epigenetically active" L1s. Notably, some of epigenetically active L1s were retrotransposition competent, which even had chimeric transcripts from the antisense promoters at their 5'UTRs. We also identified differentially methylated L1s in the prefrontal cortices of patients with psychiatric disorders. In nonneurons of bipolar disorder patients, one L1 was significantly hypomethylated and showed an inverse correlation with the expression level of the overlapping gene NREP. Finally, we observed that altered DNA methylation levels of L1 in patients with psychiatric disorders were not affected by the surrounding genomic regions but originated from the L1 sequences. These results suggested that altered epigenetic regulation of the L1 5'UTR in the brain was involved in the pathophysiology of psychiatric disorders.

Human LINE-1 retrotransposons: impacts on the genome and regulation by host factors

Genome sequencing revealed that nearly half of the human genome is comprised of transposable elements. Although most of these elements have been rendered inactive due to mutations, full-length intact long interspersed element-1 (LINE-1 or L1) copies retain the ability to mobilize through RNA intermediates by a so-called "copy-and-paste" mechanism, termed retrotransposition. L1 is the only known autonomous mobile genetic element in the genome, and its retrotransposition contributes to inter- or intra-individual genetic variation within the human population. However, L1 retrotransposition also poses a threat to genome integrity due to gene disruption and chromosomal instability. Moreover, recent studies suggest that aberrant L1 expression can impact human health by causing diseases such as cancer and chronic inflammation that might lead to autoimmune disorders. To counteract these adverse effects, the host cells have evolved multiple layers of defense mechanisms at the epigenetic, RNA and protein levels. Intriguingly, several host factors have also been reported to facilitate L1 retrotransposition, suggesting that there is competition between negative and positive regulation of L1 by host factors. Here, we summarize the known host proteins that regulate L1 activity at different stages of the replication cycle and discuss how these factors modulate disease-associated phenotypes caused by L1.

Pan cancer characterization of genes whose expression has been associated with LINE-1 antisense promoter activity

BACKGROUND

Long interspersed nuclear element-1 (LINE-1 or L1) comprises 17% of the human genome. As the only autonomous and active retrotransposons, L1 may take part in cancer initiation and progression in some ways. The studies of L1 in cancer mainly focus on the impact of L1 insertion into the new genome locus. The L1 5´ untranslated region (UTR) also contains antisense promoter (ASP) activity, generating L1-gene chimeric transcripts to a neighbor exon. Some of these ASP-associated genes have been reported to be overexpressed in cancer and promote cancer cell growth. However, little is known about overall expression patterns and the roles of L1 ASP-associated genes in human cancers.

RESULTS

L1 ASP-associated genes were frequently dysregulated in cancer and associated with the cell cycle, the PI3K/AKT pathway, and the GTPase signaling pathway. The expression of L1 ASP-associated genes was correlated with tumor patient prognosis. Hub L1 ASP-associated genes CENPU and MCM2 showed a correlation with immune infiltration, clinical T stage, and cancer stemness in pan-cancer. Knockdown of L1 ASP-associated gene LINC00491 resulted in a significant decrease in tumor growth and migration ability.

CONCLUSIONS

The expression of L1 ASP-associated genes is significantly dysregulated at the pan-cancer level, which is closely related to the tumor microenvironment, progression, and patient prognosis. Hub genes CENPU and MCM2 are expected to be new tumor diagnostic markers and therapeutic targets.

IS21 family transposase cleaved donor complex traps two right-handed superhelical crossings

Transposases are ubiquitous enzymes that catalyze DNA rearrangement events with broad impacts on gene expression, genome evolution, and the spread of drug-resistance in bacteria. Here, we use biochemical and structural approaches to define the molecular determinants by which IstA, a transposase present in the widespread IS21 family of mobile elements, catalyzes efficient DNA transposition. Solution studies show that IstA engages the transposon terminal sequences to form a high-molecular weight complex and promote DNA integration. A 3.4 Å resolution structure of the transposase bound to transposon ends corroborates our biochemical findings and reveals that IstA self-assembles into a highly intertwined tetramer that synapses two supercoiled terminal inverted repeats. The three-dimensional organization of the IstA•DNA cleaved donor complex reveals remarkable similarities with retroviral integrases and classic transposase systems, such as Tn7 and bacteriophage Mu, and provides insights into IS21 transposition.

Strand Asymmetries Across Genomic Processes

Across biological systems, a number of genomic processes, including transcription, replication, DNA repair, and transcription factor binding, display intrinsic directionalities. These directionalities are reflected in the asymmetric distribution of nucleotides, motifs, genes, transposon integration sites, and other functional elements across the two complementary strands. Strand asymmetries, including GC skews and mutational biases, have shaped the nucleotide composition of diverse organisms. The investigation of strand asymmetries often serves as a method to understand underlying biological mechanisms, including protein binding preferences, transcription factor interactions, retrotransposition, DNA damage and repair preferences, transcription-replication collisions, and mutagenesis mechanisms. Research into this subject also enables the identification of functional genomic sites, such as replication origins and transcription start sites. Improvements in our ability to detect and quantify DNA strand asymmetries will provide insights into diverse functionalities of the genome, the contribution of different mutational mechanisms in germline and somatic mutagenesis, and our knowledge of genome instability and evolution, which all have significant clinical implications in human disease, including cancer. In this review, we describe key developments that have been made across the field of genomic strand asymmetries, as well as the discovery of associated mechanisms.

An update on post-transcriptional regulation of retrotransposons

Retrotransposons, including LINE-1, Alu, SVA, and endogenous retroviruses, are one of the major constituents of human genomic repetitive sequences. Through the process of retrotransposition, some of them occasionally insert into new genomic locations by a copy-paste mechanism involving RNA intermediates. Irrespective of de novo genomic insertions, retrotransposon expression can lead to DNA double-strand breaks and stimulate cellular innate immunity through endogenous patterns. As a result, retrotransposons are tightly regulated by multi-layered regulatory processes to prevent the dangerous effects of their expression. In recent years, significant progress was made in revealing how retrotransposon biology intertwines with general post-transcriptional RNA metabolism. Here, I summarize current knowledge on the involvement of post-transcriptional factors in the biology of retrotransposons, focusing on LINE-1. I emphasize general RNA metabolisms such as methylation of adenine (m⁶ A), RNA 3'-end polyadenylation and uridylation, RNA decay and translation regulation. I discuss the effects of retrotransposon RNP sequestration in cytoplasmic bodies and autophagy. Finally, I summarize how innate immunity restricts retrotransposons and how retrotransposons make use of cellular enzymes, including the DNA repair machinery, to complete their replication cycles.

Transposable Elements Co-Option in Genome Evolution and Gene Regulation

The genome is no longer deemed as a fixed and inert item but rather as a moldable matter that is continuously evolving and adapting. Within this frame, Transposable Elements (TEs), ubiquitous, mobile, repetitive elements, are considered an alive portion of the genomes to date, whose functions, although long considered "dark", are now coming to light. Here we will review that, besides the detrimental effects that TE mobilization can induce, TEs have shaped genomes in their current form, promoting genome sizing, genomic rearrangements and shuffling of DNA sequences. Although TEs are mostly represented in the genomes by evolutionarily old, short, degenerated, and sedentary fossils, they have been thoroughly co-opted by the hosts as a prolific and original source of regulatory instruments for the control of gene transcription and genome organization in the nuclear space. For these reasons, the deregulation of TE expression and/or activity is implicated in the onset and progression of several diseases. It is likely that we have just revealed the outermost layers of TE functions. Further studies on this portion of the genome are required to unlock novel regulatory functions that could also be exploited for diagnostic and therapeutic approaches.

LINE-1 Retrotransposition Assays in Embryonic Stem Cells

The ongoing mobilization of active non-long terminal repeat (LTR) retrotransposons continues to impact the genomes of most mammals, including humans and rodents. Non-LTR retrotransposons mobilize using an intermediary RNA and a copy-and-paste mechanism termed retrotransposition. Non-LTR retrotransposons are subdivided into long and short interspersed elements (LINEs and SINEs, respectively), depending on their size and autonomy; while active class 1 LINEs (LINE-1s or L1s) encode the enzymatic machinery required to mobilize in cis, active SINEs use the enzymatic machinery of active LINE-1s to mobilize in trans. The mobilization mechanism used by LINE-1s/SINEs was exploited to develop ingenious plasmid-based retrotransposition assays in cultured cells, which typically exploit a reporter gene that can only be activated after a round of retrotransposition. Retrotransposition assays, in cis or in trans, are instrumental tools to study the biology of mammalian LINE-1s and SINEs. In fact, these and other biochemical/genetic assays were used to uncover that endogenous mammalian LINE-1s/SINEs naturally retrotranspose during early embryonic development. However, embryonic stem cells (ESCs) are typically used as a cellular model in these and other studies interrogating LINE-1/SINE expression/regulation during early embryogenesis. Thus, human and mouse ESCs represent an excellent model to understand how active retrotransposons are regulated and how their activity impacts the germline. Here, we describe robust and quantitative protocols to study human/mouse LINE-1 (in cis) and SINE (in trans) retrotransposition using (human and mice) ESCs. These protocols are designed to study the mobilization of active non-LTR retrotransposons in a cellular physiologically relevant context.

Transcriptomic study reveals lncRNA-mediated downregulation of innate immune and inflammatory response in the SARS-CoV-2 vaccination breakthrough infections

INTRODUCTION

The emergence of multiple variants of concerns (VOCs) with higher number of Spike mutations have led to enhanced immune escape by the SARS-CoV-2. With the increasing number of vaccination breakthrough (VBT) infections, it is important to understand the possible reason/s of the breakthrough infections.

METHODS

We performed transcriptome sequencing of 57 VBT and unvaccinated COVID-19 patients, followed by differential expression and co-expression analysis of the lncRNAs and the mRNAs. The regulatory mechanism was highlighted by analysis towards repeat element distribution within the co-expressed lncRNAs, followed by repeats driven homologous interaction between the lncRNAs and the promoter regions of genes from the same topologically associated domains (TAD).

RESULTS

We identified 727 differentially expressed lncRNAs (153 upregulated and 574 downregulated) and 338 mRNAs (34 up- and 334 downregulated) in the VBT patients. This includes LUCAT1, MALAT1, ROR1-AS1, UGDH-AS1 and LINC00273 mediated modulation of immune response, whereas MALAT1, NEAT1 and GAS5 regulated inflammatory response in the VBT. LncRNA-mRNA co-expression analysis highlighted 34 lncRNAs interacting with 267 mRNAs. We also observed a higher abundance of Alu, LINE1 and LTRs within the interacting lncRNAs of the VBT patients. These interacting lncRNAs have higher interaction with the promoter region of the genes from the same TAD, compared to the non-interacting lncRNAs with the enrichment of Alu and LINE1 in the gene promoter.

DISCUSSION

Significant downregulation and GSEA of the TAD gene suggest Alu and LINE1 driven homologous interaction between the lncRNAs and the TAD genes as a possible mechanism of lncRNA-mediated suppression of innate immune/inflammatory responses and activation of adaptive immune response. The lncRNA-mediated suppression of innate immune/inflammatory responses and activation of adaptive immune response might explain the SARS-CoV-2 breakthrough infections with milder symptoms in the VBT. Besides, the study also highlights repeat element mediated regulation of genes in 3D as another possible way of lncRNA-mediated immune-regulation modulating vaccination breakthroughs milder disease phenotype and shorter hospital stay.

SCIFER: approach for analysis of LINE-1 mRNA expression in single cells at a single locus resolution

BACKGROUND

Endogenous expression of L1 mRNA is the first step in an L1-initiated mutagenesis event. However, the contribution of individual cell types to patterns of organ-specific L1 mRNA expression remains poorly understood, especially at single-locus resolution. We introduce a method to quantify expression of mobile elements at the single-locus resolution in scRNA-Seq datasets called Single Cell Implementation to Find Expressed Retrotransposons (SCIFER). SCIFER aligns scRNA-Seq reads uniquely to the genome and extracts alignments from single cells by cell-specific barcodes. In contrast to the alignment performed using default parameters, this alignment strategy increases accuracy of L1 locus identification by retaining only reads that are uniquely mapped to individual L1 loci. L1 loci expressed in single cells are unambiguously identified using a list of L1 loci manually validated to be expressed in bulk RNA-Seq datasets generated from the same cell line or organ.

RESULTS

Validation of SCIFER using MCF7 cells determined technical parameters needed for optimal detection of L1 expression in single cells. We show that unsupervised analysis of L1 expression in single cells exponentially inflates both the levels of L1 expression and the number of expressed L1 loci. Application of SCIFER to analysis of scRNA-Seq datasets generated from mouse and human testes identified that mouse Round Spermatids and human Spermatogonia, Spermatocytes, and Round Spermatids express the highest levels of L1 mRNA. Our analysis also determined that similar to mice, human testes from unrelated individuals share as much as 80% of expressed L1 loci. Additionally, SCIFER determined that individual mouse cells co-express different L1 sub-families and different families of transposable elements, experimentally validating their co-existence in the same cell.

CONCLUSIONS

SCIFER detects mRNA expression of individual L1 loci in single cells. It is compatible with scRNA-Seq datasets prepared using traditional sequencing methods. Validated using a human cancer cell line, SCIFER analysis of mouse and human testes identified key cell types supporting L1 expression in these species. This will further our understanding of differences and similarities in endogenous L1 mRNA expression patterns in mice and humans.

Comprehensive analysis of full-length transcripts reveals novel splicing abnormalities and oncogenic transcripts in liver cancer

Genes generate transcripts of various functions by alternative splicing. However, in most transcriptome studies, short-reads sequencing technologies (next-generation sequencers) have been used, leaving full-length transcripts unobserved directly. Although long-reads sequencing technologies would enable the sequencing of full-length transcripts, the data analysis is difficult. In this study, we developed an analysis pipeline named SPLICE and analyzed cDNA sequences from 42 pairs of hepatocellular carcinoma (HCC) and matched non-cancerous livers with an Oxford Nanopore sequencer. Our analysis detected 46,663 transcripts from the protein-coding genes in the HCCs and the matched non-cancerous livers, of which 5,366 (11.5%) were novel. A comparison of expression levels identified 9,933 differentially expressed transcripts (DETs) in 4,744 genes. Interestingly, 746 genes with DETs, including the LINE1-MET transcript, were not found by a gene-level analysis. We also found that fusion transcripts of transposable elements and hepatitis B virus (HBV) were overexpressed in HCCs. In vitro experiments on DETs showed that LINE1-MET and HBV-human transposable elements promoted cell growth. Furthermore, fusion gene detection showed novel recurrent fusion events that were not detected in the short-reads. These results suggest the efficiency of full-length transcriptome studies and the importance of splicing variants in carcinogenesis.

Genes generate transcripts of various functions by alternative splicing. However, in most transcriptome studies, short-reads sequencing technologies (next-generation sequencers) have been used, leaving full-length transcripts unobserved directly. In this study, we developed an analysis pipeline named SPLICE for long-read transcriptome sequencing and analyzed cDNA sequences from 42 pairs of hepatocellular carcinoma (HCC), and matched non-cancerous livers with an Oxford Nanopore sequencer. Our analysis detected 5,366 novel transcripts and 9,933 differentially expressed transcripts in 4,744 genes between HCCs and non-cancerous livers. An analysis of hepatitis B virus (HBV) transcripts showed that fusion transcripts of the HBV gene and human transposable elements were overexpressed in HBV-infected HCCs. We also identified fusion genes that were not found in the short-reads. These results suggest that long-reads sequencing technologies provide a fuller understanding of cancer transcripts and that our method contributes to the analysis of transcriptome sequences by such technologies.

LINE-1 promotes tumorigenicity and exacerbates tumor progression via stimulating metabolism reprogramming in non-small cell lung cancer

Background

Long Interspersed Nuclear Element-1 (LINE-1, L1) is increasingly regarded as a genetic risk for lung cancer. Transcriptionally active LINE-1 forms a L1-gene chimeric transcript (LCTs), through somatic L1 retrotransposition (LRT) or L1 antisense promoter (L1-ASP) activation, to play an oncogenic role in cancer progression.

Methods

Here, we developed Retrotransposon-gene fusion estimation program (ReFuse), to identify and quantify LCTs in RNA sequencing data from TCGA lung cancer cohort (n = 1146) and a single cell RNA sequencing dataset then further validated those LCTs in an independent cohort (n = 134). We next examined the functional roles of a cancer specific LCT (L1-FGGY) in cell proliferation and tumor progression in LUSC cell lines and mice.

Results

The LCT events correspond with specific metabolic processes and mitochondrial functions and was associated with genomic instability, hypomethylation, tumor stage and tumor immune microenvironment (TIME). Functional analysis of a tumor specific and frequent LCT involving FGGY (L1-FGGY) reveal that the arachidonic acid (AA) metabolic pathway was activated by the loss of FGGY through the L1-FGGY chimeric transcript to promote tumor growth, which was effectively targeted by a combined use of an anti-HIV drug (NVR) and a metabolic inhibitor (ML355). Lastly, we identified a set of transcriptomic signatures to stratify the LUSC patients with a higher risk for poor outcomes who may benefit from treatments using NVR alone or combined with an anti-metabolism drug.

Conclusions

This study is the first to characterize the role of L1 in metabolic reprogramming of lung cancer and provide rationale for L1-specifc prognosis and potential for a therapeutic strategy for treating lung cancer.

Trial registration

Study on the mechanisms of the mobile element L1-FGGY promoting the proliferation, invasion and immune escape of lung squamous cell carcinoma through the 12-LOX/Wnt pathway, Ek2020111. Registered 27 March 2020 ‐ Retrospectively registered.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12943-022-01618-5.

The Role of Transposable Elements of the Human Genome in Neuronal Function and Pathology

Transposable elements (TEs) have been extensively studied for decades. In recent years, the introduction of whole-genome and whole-transcriptome approaches, as well as single-cell resolution techniques, provided a breakthrough that uncovered TE involvement in host gene expression regulation underlying multiple normal and pathological processes. Of particular interest is increased TE activity in neuronal tissue, and specifically in the hippocampus, that was repeatedly demonstrated in multiple experiments. On the other hand, numerous neuropathologies are associated with TE dysregulation. Here, we provide a comprehensive review of literature about the role of TEs in neurons published over the last three decades. The first chapter of the present review describes known mechanisms of TE interaction with host genomes in general, with the focus on mammalian and human TEs; the second chapter provides examples of TE exaptation in normal neuronal tissue, including TE involvement in neuronal differentiation and plasticity; and the last chapter lists TE-related neuropathologies. We sought to provide specific molecular mechanisms of TE involvement in neuron-specific processes whenever possible; however, in many cases, only phenomenological reports were available. This underscores the importance of further studies in this area.

Subfamily-specific differential contribution of individual monomers and the tether sequence to mouse L1 promoter activity

BACKGROUND

The internal promoter in L1 5'UTR is critical for autonomous L1 transcription and initiating retrotransposition. Unlike the human genome, which features one contemporarily active subfamily, four subfamilies (A_I, Gf_I and Tf_I/II) have been amplifying in the mouse genome in the last one million years. Moreover, mouse L1 5'UTRs are organized into tandem repeats called monomers, which are separated from ORF1 by a tether domain. In this study, we aim to compare promoter activities across young mouse L1 subfamilies and investigate the contribution of individual monomers and the tether sequence.

RESULTS

We observed an inverse relationship between subfamily age and the average number of monomers among evolutionarily young mouse L1 subfamilies. The youngest subgroup (A_I and Tf_I/II) on average carry 3-4 monomers in the 5'UTR. Using a single-vector dual-luciferase reporter assay, we compared promoter activities across six L1 subfamilies (A_I/II, Gf_I and Tf_I/II/III) and established their antisense promoter activities in a mouse embryonic fibroblast cell line and a mouse embryonal carcinoma cell line. Using consensus promoter sequences for three subfamilies (A_I, Gf_I and Tf_I), we dissected the differential roles of individual monomers and the tether domain in L1 promoter activity. We validated that, across multiple subfamilies, the second monomer consistently enhances the overall promoter activity. For individual promoter components, monomer 2 is consistently more active than the corresponding monomer 1 and/or the tether for each subfamily. Importantly, we revealed intricate interactions between monomer 2, monomer 1 and tether domains in a subfamily-specific manner. Furthermore, using three-monomer 5'UTRs, we established a complex nonlinear relationship between the length of the outmost monomer and the overall promoter activity.

CONCLUSIONS

The laboratory mouse is an important mammalian model system for human diseases as well as L1 biology. Our study extends previous findings and represents an important step toward a better understanding of the molecular mechanism controlling mouse L1 transcription as well as L1's impact on development and disease.

A hypothesis: Retrotransposons as a relay of epigenetic marks in intergenerational epigenetic inheritance

Epigenetic marks in gametes, which both respond to the parental environmental factors and shape offspring phenotypes, are usually positioned to mediate intergenerational or transgenerational epigenetic inheritance. Nonetheless, the mechanisms through which gametic epigenetic signatures encode parental acquired phenotypes, and further initiate a cascade of molecular events to affect offspring phenotypes during early embryonic development, remain unclear. Retrotransposons are mobile DNA elements that could resist to genomic epigenetic reprogramming at specific loci and rewire the core regulatory networks of embryogenesis. Increasing evidences show that retrotransposons in the embryonic genome could interact with gametic epigenetic marks, which provides a tentative possibility that retrotransposons may serve as a relay of gametic epigenetic marks to transmit parental acquired traits. Here, we summarize the recent progress in exploring the crosstalk between gametic epigenetic marks and retrotransposons, and the regulation of gene expression and early embryonic development by retrotransposons. Accordingly, deciphering the mystery of interactions between gametic epigenetic marks and retrotransposons during early embryonic development will provide valuable insights into the intergenerational or transgenerational transmission of acquired traits.

L1 chimeric transcripts are expressed in healthy brain and their deregulation in glioma follows that of their host locus

Besides the consequences of retrotransposition, long interspersed element 1 (L1) retrotransposons can affect the host genome through their antisense promoter. In addition to the sense promoter, the evolutionarily recent L1 retrotransposons, which are present in several thousand copies, also possess an anti-sense promoter that can produce L1 chimeric transcripts (LCT) composed of the L1 5′ UTR followed by the adjacent genomic sequence. The full extent to which LCT expression occurs in a given tissue and whether disruption of the defense mechanisms that normally control L1 retrotransposons affects their expression and function in cancer cells, remain to be established. By using CLIFinder, a dedicated bioinformatics tool, we found that LCT expression was widespread in normal brain and aggressive glioma samples, and that approximately 17% of recent L1 retrotransposons, from the L1PA1 to L1PA7 subfamilies, were involved in their production. Importantly, the transcriptional activities of the L1 antisense promoters and of their host loci were coupled. Accordingly, we detected LCT-producing L1 retrotransposons mainly in transcriptionally active genes and genomic loci. Moreover, changes in the host genomic locus expression level in glioma were associated with a similar change in LCT expression level, regardless of the L1 promoter methylation status. Our findings support a model in which the host genomic locus transcriptional activity is the main driving force of LCT expression. We hypothesize that this model is more applicable when host gene and LCT are transcribed from the same strand.

Analysis of epigenetic features characteristic of L1 loci expressed in human cells

Only a select few L1 loci in the human genome are expressed in any given cell line or organ, likely to minimize damage done to the genome. The epigenetic features and requirements of expressed L1 loci are currently unknown. Using human cells and comprehensive epigenetic analysis of individual expressed and unexpressed L1 loci, we determined that endogenous L1 transcription depends on a combination of epigenetic factors, including open chromatin, activating histone modifications, and hypomethylation at the L1 promoter. We demonstrate that the L1 promoter seems to require interaction with enhancer elements for optimal function. We utilize epigenetic context to predict the expression status of L1Hs loci that are poorly mappable with RNA-Seq. Our analysis identified a population of ‘transitional’ L1 loci that likely have greater potential to be activated during the epigenetic dysregulation seen in tumors and during aging because they are the most responsive to targeted CRISPR-mediated delivery of trans-activating domains. We demonstrate that an engineered increase in endogenous L1 mRNA expression increases Alu mobilization. Overall, our findings present the first global and comprehensive analysis of epigenetic status of individual L1 loci based on their expression status and demonstrate the importance of epigenetic context for L1 expression heterogeneity.

The non-LTR retrotransposons of Entamoeba histolytica: genomic organization and biology

Genome sequence analysis of Entamoeba species revealed various classes of transposable elements. While E. histolytica and E. dispar are rich in non-long terminal repeat (LTR) retrotransposons, E. invadens contains predominantly DNA transposons. Non-LTR retrotransposons of E. histolytica constitute three families of long interspersed nuclear elements (LINEs), and their short, nonautonomous partners, SINEs. They occupy ~ 11% of the genome. The EhLINE1/EhSINE1 family is the most abundant and best studied. EhLINE1 is 4.8 kb, with two ORFs that encode functions needed for retrotransposition. ORF1 codes for the nucleic acid-binding protein, and ORF2 has domains for reverse transcriptase (RT) and endonuclease (EN). Most copies of EhLINEs lack complete ORFs. ORF1p is expressed constitutively, but ORF2p is not detected. Retrotransposition could be demonstrated upon ectopic over expression of ORF2p, showing that retrotransposition machinery is functional. The newly retrotransposed sequences showed a high degree of recombination. In transcriptomic analysis, RNA-Seq reads were mapped to individual EhLINE1 copies. Although full-length copies were transcribed, no full-length 4.8 kb transcripts were seen. Rather, sense transcripts mapped to ORF1, RT and EN domains. Intriguingly, there was strong antisense transcription almost exclusively from the RT domain. These unique features of EhLINE1 could serve to attenuate retrotransposition in E. histolytica.

Unbiased proteomic mapping of the LINE-1 promoter using CRISPR Cas9

BACKGROUND

The autonomous retroelement Long Interspersed Element-1 (LINE-1) mobilizes though a copy and paste mechanism using an RNA intermediate (retrotransposition). Throughout human evolution, around 500,000 LINE-1 sequences have accumulated in the genome. Most of these sequences belong to ancestral LINE-1 subfamilies, including L1PA2-L1PA7, and can no longer mobilize. Only a small fraction of LINE-1 sequences, approximately 80 to 100 copies belonging to the L1Hs subfamily, are complete and still capable of retrotransposition. While silenced in most cells, many questions remain regarding LINE-1 dysregulation in cancer cells.

RESULTS

Here, we optimized CRISPR Cas9 gRNAs to specifically target the regulatory sequence of the L1Hs 5'UTR promoter. We identified three gRNAs that were more specific to L1Hs, with limited binding to older LINE-1 sequences (L1PA2-L1PA7). We also adapted the C-BERST method (dCas9-APEX2 Biotinylation at genomic Elements by Restricted Spatial Tagging) to identify LINE-1 transcriptional regulators in cancer cells. Our LINE-1 C-BERST screen revealed both known and novel LINE-1 transcriptional regulators, including CTCF, YY1 and DUSP1.

CONCLUSION

Our optimization and evaluation of gRNA specificity and application of the C-BERST method creates a tool for studying the regulatory mechanisms of LINE-1 in cancer. Further, we identified the dual specificity protein phosphatase, DUSP1, as a novel regulator of LINE-1 transcription.

TFs for TEs: the transcription factor repertoire of mammalian transposable elements

In this review, Hermant and Torres-Padilla summarize and discuss the transcription factors known to be involved in the sequence-specific recognition and transcriptional activation of specific transposable element families or subfamilies.

Transposable elements (TEs) are genetic elements capable of changing position within the genome. Although their mobilization can constitute a threat to genome integrity, nearly half of modern mammalian genomes are composed of remnants of TE insertions. The first critical step for a successful transposition cycle is the generation of a full-length transcript. TEs have evolved cis-regulatory elements enabling them to recruit host-encoded factors driving their own, selfish transcription. TEs are generally transcriptionally silenced in somatic cells, and the mechanisms underlying their repression have been extensively studied. However, during germline formation, preimplantation development, and tumorigenesis, specific TE families are highly expressed. Understanding the molecular players at stake in these contexts is of utmost importance to establish the mechanisms regulating TEs, as well as the importance of their transcription to the biology of the host. Here, we review the transcription factors known to be involved in the sequence-specific recognition and transcriptional activation of specific TE families or subfamilies. We discuss the diversity of TE regulatory elements within mammalian genomes and highlight the importance of TE mobilization in the dispersal of transcription factor-binding sites over the course of evolution.

The RNA editing enzyme ADAR2 restricts L1 mobility

Adenosine deaminases acting on RNA (ADARs) are enzymes that convert adenosines to inosines in double-stranded RNAs (RNA editing A-to-I). ADAR1 and ADAR2 were previously reported as HIV-1 proviral factors. The aim of this study was to investigate the composition of the ADAR2 ribonucleoprotein complex during HIV-1 expression. By using a dual-tag affinity purification procedure in cells expressing HIV-1 followed by mass spectrometry analysis, we identified 10 non-ribosomal ADAR2-interacting factors. A significant fraction of these proteins was previously found associated to the Long INterspersed Element 1 (LINE1 or L1) ribonucleoparticles and to regulate the life cycle of L1 retrotransposons. Considering that we previously demonstrated that ADAR1 is an inhibitor of LINE-1 retrotransposon activity, we investigated whether also ADAR2 played a similar function. To reach this goal, we performed specific cell culture retrotransposition assays in cells overexpressing or ablated for ADAR2. These experiments unveil a novel function of ADAR2 as suppressor of L1 retrotransposition. Furthermore, we showed that ADAR2 binds the basal L1 RNP complex.

Overall, these data support the role of ADAR2 as regulator of L1 life cycle.

Intracellular DNA transfer events restricted to the genus Convallaria within the Asparagaceae family: Possible mechanisms and potential as genetic markers for biographical studies

Intracellular gene transfer among plant genomes is a common phenomenon. Due to their high conservation and high plastid membrane integrity, chloroplast (cp) genomes incorporate foreign genetic material very rarely. Convallaria is a small monocotyledonous genus consisting of C. keiskei, C. majalis and C. montana. Here, we characterized, analyzed and identified 3.3 and 3.7 kb of mitochondrial DNA sequences in the plastome (MCP) of C. majalis and C. montana, respectively. We identified 6 bp and 23 bp direct repeats and mitochondrial pseudogenes, with rps3, rps19 and rpl10 identified in the MCP region. Additionally, we developed novel plastid molecular genetic markers to differentiate Convallaria spp. based on 21 populations. BEAST and biogeographical analyses suggested that Convallaria separated into Eurasian and North American lineages during the middle Pliocene and originated in East Asia. Vicariance in the genus was followed by dispersal into Europe and southeastern North America. These analyses indicate that the MCP event was restricted to the genus Convallaria of Asparagaceae, in contrast to similar events that occurred in its common ancestors with other families of land plants. However, further mitochondrial and population studies are necessary to understand the integration of the MCP region and gene flow in the genus Convallaria.

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

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

Endogenous retroviruses in the origins and treatment of cancer

Endogenous retroviruses (ERVs) are emerging as promising therapeutic targets in cancer. As remnants of ancient retroviral infections, ERV-derived regulatory elements coordinate expression from gene networks, including those underpinning embryogenesis and immune cell function. ERV activation can promote an interferon response, a phenomenon termed viral mimicry. Although ERV expression is associated with cancer, and provisionally with autoimmune and neurodegenerative diseases, ERV-mediated inflammation is being explored as a way to sensitize tumors to immunotherapy. Here we review ERV co-option in development and innate immunity, the aberrant contribution of ERVs to tumorigenesis, and the wider biomedical potential of therapies directed at ERVs.

Sensing of transposable elements by the antiviral innate immune system

Around half of the genome in mammals is composed of transposable elements (TEs) such as DNA transposons and retrotransposons. Several mechanisms have evolved to prevent their activity and the detrimental impact of their insertional mutagenesis. Despite these potentially negative effects, TEs are essential drivers of evolution, and in certain settings, beneficial to their hosts. For instance, TEs have rewired the antiviral gene regulatory network and are required for early embryonic development. However, due to structural similarities between TE-derived and viral nucleic acids, cells can misidentify TEs as invading viruses and trigger the major antiviral innate immune pathway, the type I interferon (IFN) response. This review will focus on the different settings in which the role of TE-mediated IFN activation has been documented, including cancer and senescence. Importantly, TEs may also play a causative role in the development of complex autoimmune diseases characterised by constitutive type I IFN activation. All these observations suggest the presence of strong but opposing forces driving the coevolution of TEs and antiviral defence. A better biological understanding of the TE replicative cycle as well as of the antiviral nucleic acid sensing mechanisms will provide insights into how these two biological processes interact and will help to design better strategies to treat human diseases characterised by aberrant TE expression and/or type I IFN activation.

Integrated transcription factor profiling with transcriptome analysis identifies L1PA2 transposons as global regulatory modulators in a breast cancer model

While transposons are generally silenced in somatic tissues, many transposons escape epigenetic repression in epithelial cancers, become transcriptionally active and contribute to the regulation of human gene expression. We have developed a bioinformatic pipeline for the integrated analysis of transcription factor binding and transcriptomic data to identify transposon-derived promoters that are activated in specific diseases and developmental states. We applied this pipeline to a breast cancer model, and found that the L1PA2 transposon subfamily contributes abundant regulatory sequences to co-ordinated transcriptional regulation in breast cancer. Transcription factor profiling demonstrates that over 27% of L1PA2 transposons harbour co-localised binding sites of functionally interacting, cancer-associated transcription factors in MCF7 cells, a cell line used to model breast cancer. Transcriptomic analysis reveals that L1PA2 transposons also contribute transcription start sites to up-regulated transcripts in MCF7 cells, including some transcripts with established oncogenic properties. In addition, we verified the utility of our pipeline on other transposon subfamilies, as well as on leukemia and lung carcinoma cell lines. We demonstrate that the normally quiescent regulatory activities of transposons can be activated and alter the cancer transcriptome. In particular, the L1PA2 subfamily contributes abundant regulatory sequences, and likely plays a global role in modulating breast cancer transcriptional regulation. Understanding the regulatory impact of L1PA2 on breast cancer genomes provides additional insights into cancer genome regulation, and may provide novel biomarkers for disease diagnosis, prognosis and therapy.

Genomic Analysis Revealed a Convergent Evolution of LINE-1 in Coat Color: A Case Study in Water Buffaloes (Bubalus bubalis)

Visible pigmentation phenotypes can be used to explore the regulation of gene expression and the evolution of coat color patterns in animals. Here, we performed whole-genome and RNA sequencing and applied genome-wide association study, comparative population genomics and biological experiments to show that the 2,809-bp-long LINE-1 insertion in the ASIP (agouti signaling protein) gene is the causative mutation for the white coat phenotype in swamp buffalo (Bubalus bubalis). This LINE-1 insertion (3' truncated and containing only 5' UTR) functions as a strong proximal promoter that leads to a 10-fold increase in the transcription of ASIP in white buffalo skin. The 165 bp of 5' UTR transcribed from the LINE-1 is spliced into the first coding exon of ASIP, resulting in a chimeric transcript. The increased expression of ASIP prevents melanocyte maturation, leading to the absence of pigment in white buffalo skin and hairs. Phylogenetic analyses indicate that the white buffalo-specific ASIP allele originated from a recent genetic transposition event in swamp buffalo. Interestingly, as a similar LINE-1 insertion has been identified in the cattle ASIP gene, we discuss the convergent mechanism of coat color evolution in the Bovini tribe.

Role of Transposable Elements in Gene Regulation in the Human Genome

Transposable elements (TEs), also known as mobile elements (MEs), are interspersed repeats that constitute a major fraction of the genomes of higher organisms. As one of their important functional impacts on gene function and genome evolution, TEs participate in regulating the expression of genes nearby and even far away at transcriptional and post-transcriptional levels. There are two known principal ways by which TEs regulate the expression of genes. First, TEs provide cis-regulatory sequences in the genome with their intrinsic regulatory properties for their own expression, making them potential factors for regulating the expression of the host genes. TE-derived cis-regulatory sites are found in promoter and enhancer elements, providing binding sites for a wide range of trans-acting factors. Second, TEs encode for regulatory RNAs with their sequences showed to be present in a substantial fraction of miRNAs and long non-coding RNAs (lncRNAs), indicating the TE origin of these RNAs. Furthermore, TEs sequences were found to be critical for regulatory functions of these RNAs, including binding to the target mRNA. TEs thus provide crucial regulatory roles by being part of cis-regulatory and regulatory RNA sequences. Moreover, both TE-derived cis-regulatory sequences and TE-derived regulatory RNAs have been implicated in providing evolutionary novelty to gene regulation. These TE-derived regulatory mechanisms also tend to function in a tissue-specific fashion. In this review, we aim to comprehensively cover the studies regarding these two aspects of TE-mediated gene regulation, mainly focusing on the mechanisms, contribution of different types of TEs, differential roles among tissue types, and lineage-specificity, based on data mostly in humans.

Transcriptomic analysis of Entamoeba histolytica reveals domain-specific sense strand expression of LINE-encoded ORFs with massive antisense expression of RT domain

LINEs are retrotransposable elements found in diverse organisms. Their activity is kept in check by several mechanisms, including transcriptional silencing. Here we have analyzed the transcription status of LINE1 copies in the early-branching parasitic protist Entamoeba histolytica. Full-length EhLINE1 encodes ORF1, and ORF2 with reverse transcriptase (RT) and endonuclease (EN) domains. RNA-Seq analysis of EhLINE1 copies (both truncated and full-length) showed unique features. Firstly, although 20/41 transcribed copies were full-length, we failed to detect any full-length transcripts. Rather, sense-strand transcripts mapped to the functional domains- ORF1, RT and EN. Secondly, there was strong antisense transcription specifically from RT domain. No antisense transcripts were seen from ORF1. Antisense RT transcripts did not encode known functional peptides. They could possibly be involved in attenuating translation of RT domain, as we failed to detect ORF2p, whereas ORF1p was detectable. Lack of full-length transcripts and strong antisense RT expression may serve to limit EhLINE1 retrotransposition.

Effects of activation of the LINE-1 antisense promoter on the growth of cultured cells

Long interspersed element 1 (LINE-1, or L1) is a retrotransposon that constitutes ~ 17% of the human genome. Although ~ 6000 full-length L1s spread throughout the human genome, their biological significance remains undetermined. The L1 5′ untranslated region has bidirectional promoter activity with a sense promoter driving L1 mRNA production and an antisense promoter (ASP) driving the production of L1-gene chimeric RNAs. Here, we stimulated L1 ASP activity using CRISPR-Cas9 technology to evaluate its biological impacts. Activation of the L1 ASP upregulated the expression of L1 ASP-driven ORF0 and enhanced cell growth. Furthermore, the exogenous expression of ORF0 also enhanced cell growth. These results indicate that activation of L1 ASP activity fuels cell growth at least through ORF0 expression. To our knowledge, this is the first report demonstrating the role of the L1 ASP in a biological context. Considering that L1 sequences are desilenced in various tumor cells, our results indicate that activation of the L1 ASP may be a cause of tumor growth; therefore, interfering with L1 ASP activity may be a potential strategy to suppress the growth.

Retroelement-derived RNA and its role in the brain

Comprising ~40% of the human genome, retroelements are mobile genetic elements which are transcribed into RNA, then reverse-transcribed into DNA and inserted into a new site in the genome. Retroelements are referred to as "genetic parasites", residing among host genes and relying on host machinery for transcription and evolutionary propagation. The healthy brain has the highest expression of retroelement-derived sequences compared to other somatic tissue, which leads to the question: how does retroelement-derived RNA influence human traits and cellular states? While the functional importance of upregulating retroelement expression in the brain is an active area of research, RNA species derived from retroelements influence both self- and host gene expression by contributing to chromatin remodeling, alternative splicing, somatic mosaicism and translational repression. Here, we review the emerging evidence that the functional importance of RNA derived from retroelements is multifaceted. Retroelements can influence organismal states through the seeding of epigenetic states in chromatin, the production of structured RNA and even catalytically active ribozymes, the generation of cytoplasmic ssDNA and RNA/DNA hybrids, the production of viral-like proteins, and the generation of somatic mutations. Comparative sequencing suggests that retroelements can contribute to intraspecies variation through these mechanisms to alter transcript identity and abundance. In humans, an increasing number of neurodevelopmental and neurodegenerative conditions are associated with dysregulated retroelements, including Aicardi-Goutieres syndrome (AGS), Rett syndrome (RTT), Amyotrophic Lateral Sclerosis (ALS), Alzheimer's disease (AD), multiple sclerosis (MS), schizophrenia (SZ), and aging. Taken together, these concepts suggest a larger functional role for RNA derived from retroelements. This review aims to define retroelement-derived RNA, discuss how it impacts the mammalian genome, as well as summarize data supporting phenotypic consequences of this unique RNA subset in the brain.

Impact of Tumor LINE-1 Methylation Level and Neoadjuvant Treatment and Its Association with Colorectal Cancer Survival

Recent studies suggest that long-interspersed nucleotide element-1 (LINE-1) hypomethylation is commonly found in colorectal cancer (CRC), and is associated with worse prognosis. However, the utility of LINE-1 methylation on the prognosis of CRC is still controversial, and may be due to the fact that some clinical and pathological features may affect LINE-1 methylation. Thus, the aim of this study was to assess the prognostic value of tumor LINE-1 methylation in CRC, through their association with the CRC clinical and pathological characteristics. Survival of sixty-seven CRC patients was evaluated according to the median of tumor LINE-1 methylation, as well as pathological and oncological variables. We also studied the association between LINE-1 methylation and pathological features, and finally, we assessed the overall and disease-free survival of LINE1 methylation, stratified by neoadjuvant treatment and further checked by multivariate Cox regression to assess the statistical interactions. LINE-1 was hypomethylated in the CRC tumor with respect to the tumor adjacent-free area (p < 0.05), without association with any other clinical and oncological features, nor with overall and disease-free survival rates for CRC. Relevantly, in neoadjuvant treatment, LINE-1 methylation was associated with survival rates. Thus, disease-free and overall survival rates of treated CRC patients were worse in the hypomethylated LINE-1 tumors than those with normal LINE-1 methylation (p = 0.004 and 0.0049, respectively). Indeed, LINE-1 was hypermethylated more in the treated patients than in the non-treated patients (p < 0.05). The present study showed that tumor LINE-1 hypomethylation was associated with worse survival rates in only treated patients. Our data suggest an interactive effect of neoadjuvant treatment and tumor LINE-1 methylation, which could be a specific-tissue biomarker to predict survival of the treated patients, and help to personalize treatment in CRC.

Measuring and interpreting transposable element expression

Transposable elements (TEs) are insertional mutagens that contribute greatly to the plasticity of eukaryotic genomes, influencing the evolution and adaptation of species as well as physiology or disease in individuals. Measuring TE expression helps to understand not only when and where TE mobilization can occur but also how this process alters gene expression, chromatin accessibility or cellular signalling pathways. Although genome-wide gene expression assays such as RNA sequencing include transposon-derived transcripts, most computational analytical tools discard or misinterpret TE-derived reads. Emerging approaches are improving the identification of expressed TE loci and helping to discriminate TE transcripts that permit TE mobilization from chimeric gene-TE transcripts or pervasive transcription. Here we review the main challenges associated with the detection of TE expression, including mappability, insertional and internal sequence polymorphisms, and the diversity of the TE transcriptional landscape, as well as the different experimental and computational strategies to solve them.

Identification and characterization of microRNAs (miRNAs) and their transposable element origins in the saltwater crocodile, Crocodylus porosus

MicroRNAs (miRNAs) are 18-24 nucleotide regulatory RNAs. They are involved in the regulation of genetic and biological pathways through post transcriptional gene silencing and/or translational repression. Data suggests a slow evolutionary rate for the saltwater crocodile (Crocodylus porosus) over the past several million years when compared to birds, the closest extant relatives of crocodilians. Understanding gene regulation in the saltwater crocodile in the context of relatively slow genomic change thus holds potential for the investigation of genomics, evolution, and adaptation. Utilizing eleven tissue types and sixteen small RNA libraries, we report 644 miRNAs in the saltwater crocodile with >78% of miRNAs being novel to crocodilians. We also identified potential targets for the miRNAs and analyzed the relationship of the miRNA repertoire to transposable elements (TEs). Results suggest an increased association of DNA transposons with miRNAs when compared to retrotransposons. This work reports the first comprehensive analysis of miRNAs in Crocodylus porosus and addresses the potential impacts of miRNAs in regulating the genome in the saltwater crocodile. In addition, the data suggests a supporting role of TEs as a source for miRNAs, adding to the increasing evidence that TEs play a significant role in the evolution of gene regulation.

Differential expressions of L1-chimeric transcripts in normal and matched-cancer tissues

L1s are a cis-regulatory elements and contain bidirectional internal promoters within the 5' untranslated region (UTR). L1s provide bidirectional promoters that generate alternative transcripts and affect differential expressions in the human genome. In particular, L1 antisense promoters (L1ASPs) could produce aberrant transcripts in cancer tissues compared to normal tissues. In this study, we identified the L1-chimeric transcripts derived from L1ASPs and analyzed relative expression of L1-chimeric transcripts between normal and matched-cancer tissues. First, we collected 425 L1-chimeric transcripts by referring to previous studies. Through the manual inspection, we identified 144 L1-chimeric transcripts derived from 44 L1 antisense promoters, suggesting that the antisense promoter acted as an alternative promoter. We analyzed relative gene expression levels of 16 L1-chimeric transcripts between matched cancer-normal tissue pair (lung, liver, gastric, kidney, thyroid, breast, ovary, uterus, and prostate) using real-time quantitative PCR (RT-qPCR) and investigated putative transcription factor binding motifs to determine activity of L1ASPs. Taken together, we propose that L1ASPs could contribute to the differential gene expression between normal and cancer tissues.

Long INterspersed element-1 mobility as a sensor of environmental stresses

Long INterspersed element (LINE-1, L1) retrotransposons are the most abundant transposable elements in the human genome, constituting approximately 17%. They move by a "copy-paste" mechanism, involving reverse transcription of an RNA intermediate and insertion of its cDNA copy at a new site in the genome. L1 retrotransposition (L1-RTP) can cause insertional mutations, alter gene expression, transduce exons, and induce epigenetic dysregulation. L1-RTP is generally repressed; however, a number of observations collected over about 15 years revealed that it can occur in response to environmental stresses. Moreover, emerging evidence indicates that L1-RTP can play a role in the onset of several neurological and oncological diseases in humans. In recent years, great attention has been paid to the exposome paradigm, which proposes that health effects of an environmental factor should be evaluated considering both cumulative environmental exposures and the endogenous processes resulting from the biological response. L1-RTP could be an endogenous process considered for this application. Here, we summarize the current understanding of environmental factors that can affect the retrotransposition of human L1 elements. Evidence indicates that L1-RTP alteration is triggered by numerous and various environmental stressors, such as chemical agents (heavy metals, carcinogens, oxidants, and drugs), physical agents (ionizing and non-ionizing radiations), and experiential factors (voluntary exercise, social isolation, maternal care, and environmental light/dark cycles). These data come from in vitro studies on cell lines and in vivo studies on transgenic animals: future investigations should be focused on physiologically relevant models to gain a better understanding of this topic.

sRNA/L1 retrotransposition: using siRNAs and miRNAs to expand the applications of the cell culture-based LINE-1 retrotransposition assay

The cell culture-based retrotransposition reporter assay has been (and is) an essential tool for the study of vertebrate Long INterspersed Elements (LINEs). Developed more than 20 years ago, this assay has been instrumental in characterizing the role of LINE-encoded proteins in retrotransposition, understanding how ribonucleoprotein particles are formed, how host factors regulate LINE mobilization, etc. Moreover, variations of the conventional assay have been developed to investigate the biology of other currently active human retrotransposons, such as Alu and SVA. Here, we describe a protocol that allows combination of the conventional cell culture-based LINE-1 retrotransposition reporter assay with short interfering RNAs (siRNAs) and microRNA (miRNAs) mimics or inhibitors, which has allowed us to uncover specific miRNAs and host factors that regulate retrotransposition. The protocol described here is highly reproducible, quantitative, robust and flexible, and allows the study of several small RNA classes and various retrotransposons. To illustrate its utility, here we show that siRNAs to Fanconi anaemia proteins (FANC-A and FANC-C) and an inhibitor of miRNA-20 upregulate and downregulate human L1 retrotransposition, respectively. This article is part of a discussion meeting issue 'Crossroads between transposons and gene regulation'.

Long Interspersed Nuclear Elements 1 (LINE1): The chimeric transcript L1-MET and its involvement in cancer

Long interspersed nuclear elements 1 (LINE1) are non-LTR retrotransposons that represent the greatest remodeling force of the human genome during evolution. Genomically, LINE1 are constituted by a 5´ untranslated region (UTR), where the promoter regions are located, three open reading frames (ORF0, ORF1, and ORF2) and one 3´UTR, which has a poly(A) tail that harbors the short interspersed nuclear elements (SINEs) Alu and SVA. Although the intrinsic nature of LINE1 is to be copied and inserted into the genome, an increase in their mobility produces genomic instability. In response to this, the cell has "designed" many mechanisms controlling the retrotransposition levels of LINE1; however, alterations in these regulation systems can increase LINE1 mobility and the formation of chimeric genes. Evidence indicates that 988 human genes have LINE1 inserted in their sequence, resulting in the transcriptional control of genes by their own promoters, as well as by the LINE1 antisense promoter (ASP). To date, very little is known about the biologic impact of this and the L1-MET chimera is a more or less studied case. ASP hypomethylation has been observed in all studied cancer types, leading to increased L1-MET expression. In specific types of cancer, this L1-MET increase controls both low and high MET protein levels. It remains to be clarified if this protein product is a chimeric protein.

Comparative analysis on the expression of L1 loci using various RNA-Seq preparations

Background

Retrotransposons are one of the oldest evolutionary forces shaping mammalian genomes, with the ability to mobilize from one genomic location to another. This mobilization is also a significant factor in human disease. The only autonomous human retroelement, L1, has propagated to make up 17% of the human genome, accumulating over 500,000 copies. The majority of these loci are truncated or defective with only a few reported to remain capable of retrotransposition. We have previously published a strand-specific RNA-Seq bioinformatics approach to stringently identify at the locus-specific level the few expressed full-length L1s using cytoplasmic RNA. With growing repositories of RNA-Seq data, there is potential to mine these datasets to identify and study expressed L1s at single-locus resolution, although many datasets are not strand-specific or not generated from cytoplasmic RNA.

Results

We developed whole-cell, cytoplasmic and nuclear RNA-Seq datasets from 22Rv1 prostate cancer cells to test the influence of different preparations on the quality and effort needed to measure L1 expression. We found that there was minimal data loss in the identification of full-length expressed L1 s using whole cell, strand-specific RNA-Seq data compared to cytoplasmic, strand-specific RNA-Seq data. However, this was only possible with an increased amount of manual curation of the bioinformatics output to eliminate increased background. About half of the data was lost when the sequenced datasets were non-strand specific.

Conclusions

The results of these studies demonstrate that with rigorous manual curation the utilization of stranded RNA-Seq datasets allow identification of expressed L1 loci from either cytoplasmic or whole-cell RNA-Seq datasets.

Restricted and non-essential redundancy of RNAi and piRNA pathways in mouse oocytes

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

Retrotransposons are mobile genomic parasites causing mutations. Germ cells need protection against retrotransposons to prevent heritable transmission of their new insertions. The piRNA pathway is an ancient germline defense system analogous to acquired immunity: once a retrotransposon jumps into a piRNA-producing locus, which provides a kind of a “genomic sensor” for actively transposing elements, it is recognized and suppressed. Remarkably, the murine piRNA pathway is essential for spermatogenesis but not oocyte development. In contrast, zebrafish lacking the piRNA pathway do not develop any germ cells. It was hypothesized that RNA interference pathway could rescue oocyte development in mice lacking the piRNA pathway. RNA interference also targets retrotransposons and is particularly enhanced in mouse oocytes. To test this hypothesis, we engineered mice lacking both pathways and observed that oocytes in these mice develop normally, which argues against the hypothesis. Furthermore, analysis of individual retrotransposon groups revealed that in specific cases the two pathways mutually compensate each other. However, this redundancy apparently evolved stochastically and is restricted to specific retrotransposon groups. Finally, our results indicate that there must be yet another layer of retrotransposon silencing in mouse oocytes, which prevents high retrotransposon activity in the absence of piRNA and RNA interference pathways.