L1 retrotransposition is a common feature of mammalian hepatocarcinogenesis

LINE-1 retrotransposons contribute to mouse PV interneuron development

Retrotransposons are mobile DNA sequences duplicated via transcription and reverse transcription of an RNA intermediate. Cis-regulatory elements encoded by retrotransposons can also promote the transcription of adjacent genes. Somatic LINE-1 (L1) retrotransposon insertions have been detected in mammalian neurons. It is, however, unclear whether L1 sequences are mobile in only some neuronal lineages or therein promote neurodevelopmental gene expression. Here we report programmed L1 activation by SOX6, a transcription factor critical for parvalbumin (PV) interneuron development. Mouse PV interneurons permit L1 mobilization in vitro and in vivo, harbor unmethylated L1 promoters and express full-length L1 mRNAs and proteins. Using nanopore long-read sequencing, we identify unmethylated L1s proximal to PV interneuron genes, including a novel L1 promoter-driven Caps2 transcript isoform that enhances neuron morphological complexity in vitro. These data highlight the contribution made by L1 cis-regulatory elements to PV interneuron development and transcriptome diversity, uncovered due to L1 mobility in this milieu.

Regulation and function of transposable elements in cancer genomes

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

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.

Higher incidence of embryonic defects in mouse offspring conceived with assisted reproduction from fathers with sperm epimutations

Assisted reproductive technologies (ART) account for 1-6% of births in developed countries. While most children conceived are healthy, increases in birth and genomic imprinting defects have been reported; such abnormal outcomes have been attributed to underlying parental infertility and/or the ART used. Here, we assessed whether paternal genetic and lifestyle factors, that are associated with male infertility and affect the sperm epigenome, can influence ART outcomes. We examined how paternal factors, haploinsufficiency for Dnmt3L, an important co-factor for DNA methylation reactions, and/or diet-induced obesity, in combination with ART (superovulation, in vitro fertilization, embryo culture and embryo transfer), could adversely influence embryo development and DNA methylation patterning in mice. While male mice fed high-fat diets (HFD) gained weight and showed perturbed metabolic health, their sperm DNA methylation was minimally affected by the diet. In contrast, Dnmt3L haploinsufficiency induced a marked loss of DNA methylation in sperm; notably, regions affected were associated with neurodevelopmental pathways and enriched in young retrotransposons, sequences that can have functional consequences in the next generation. Following ART, placental imprinted gene methylation and growth parameters were impacted by one or both paternal factors. For embryos conceived by natural conception, abnormality rates were similar for WT and Dnmt3L+/- fathers. In contrast, paternal Dnmt3L+/- genotype, as compared to WT fathers, resulted in a 3-fold increase in the incidence of morphological abnormalities in embryos generated by ART. Together, the results indicate that embryonic morphological and epigenetic defects associated with ART may be exacerbated in offspring conceived by fathers with sperm epimutations.

Dysregulated Repeat Element Viral-like Immune Response in Hepatocellular Carcinoma

Purpose

Dysregulation of viral-like repeat RNAs are a common feature across many malignancies that are linked with immunological response, but the characterization of these in hepatocellular carcinoma (HCC) is understudied. In this study, we performed RNA in situ hybridization (RNA-ISH) of different repeat RNAs, immunohistochemistry (IHC) for immune cell subpopulations, and spatial transcriptomics to understand the relationship of HCC repeat expression, immune response, and clinical outcomes.

Experimental Design

RNA-ISH for LINE1, HERV-K, HERV-H, and HSATII repeats and IHC for T-cell, Treg, B-cell, macrophage, and immune checkpoint markers were performed on 43 resected HCC specimens. Spatial transcriptomics on tumor and vessel regions of interest was performed on 28 specimens from the same cohort.

Results

High HERV-K and high LINE1 expression were both associated with worse overall survival. There was a positive correlation between LINE1 expression and FOXP3 T-regulatory cells (r = 0.51 p < 0.001) as well as expression of the TIM3 immune checkpoint (r = 0.34, p = 0.03). Spatial transcriptomic profiling of HERV-K high and LINE-1 high tumors identified elevated expression of multiple genes previously associated with epithelial mesenchymal transition, cellular proliferation, and worse overall prognosis in HCC including SSX1, MAGEC2, and SPINK1.

Conclusion

Repeat RNAs may serve as useful prognostic biomarkers in HCC and may also serve as novel therapeutic targets. Additional study is needed to understand the mechanisms by which repeat RNAs impact HCC tumorigenesis.

DNA methylation markers in the detection of hepatocellular carcinoma

Hepatocellular carcinoma (HCC) is the most common primary liver malignancy and has a poor prognosis. Epigenetic modification has been shown to be deregulated during HCC development by dramatically impacting the differentiation, proliferation, and function of cells. One important epigenetic modification is DNA methylation during which methyl groups are added to cytosines without changing the DNA sequence itself. Studies found that methylated DNA markers can be specific for detection of HCC. On the basis of these findings, the utility of methylated DNA markers as novel biomarkers for early-stage HCC has been measured in blood, and indeed superior sensitivity and specificity have been found in several studies when compared to current surveillance methods. However, a variety of factors currently limit the immediate application of these exciting biomarkers. In this review, we provide a detailed rationalisation of the approach and basis for the use of methylation biomarkers for HCC detection and summarise recent studies on methylated DNA markers in HCC focusing on the importance of the aetiological cause of liver disease in the mechanisms leading to cancer.

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

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

The Regulation and Immune Signature of Retrotransposons in Cancer

Advances in sequencing technologies and the bioinformatic analysis of big data facilitate the study of jumping genes' activity in the human genome in cancer from a broad perspective. Retrotransposons, which move from one genomic site to another by a copy-and-paste mechanism, are regulated by various molecular pathways that may be disrupted during tumorigenesis. Active retrotransposons can stimulate type I IFN responses. Although accumulated evidence suggests that retrotransposons can induce inflammation, the research investigating the exact mechanism of triggering these responses is ongoing. Understanding these mechanisms could improve the therapeutic management of cancer through the use of retrotransposon-induced inflammation as a tool to instigate immune responses to tumors.

Comprehensive Analysis and Drug Modulation of Human Endogenous Retrovirus in Hepatocellular Carcinomas

BACKGROUND

Human endogenous retroviruses (HERVs) play an important role in the development of cancer and many diseases. Here, we comprehensively explored the impact of HERVs on hepatocellular carcinomas (HCCs).

METHODS

We employed Telescope to identify HERVs and quantify their expression in the total RNA sequencing data obtained from 254 HCC samples, comprising 254 tumor tissues and 34 matched normal tissues.

RESULTS

In total, 3357 locus-specific activations of HERVs were differentially expressed, and 180 were correlated with patient survival. Using these 180 HERVs for classification, we found four subgroups with survival correlation. Higher expression levels of the 180 HERVs were correlated with poorer survival, while age, AFP, some mutations, and copy and structural variants differed among subgroups. The differential expression of host genes in high expression of these 180 HERVs primarily involved the activation of pathways related to immunity and infection, lipid and atherosclerosis, MAPK and NF-kB signaling, and cytokine-cytokine receptor interactions. Conversely, there was a suppression of pathways associated with RNA processing, including nucleocytoplasmic transport, surveillance and ribosome biogenesis, and transcriptional misregulation in cancer pathways. Almost all genes involved in HERV activation restriction, KRAB zinc finger proteins, RNA nucleocytoplasmic transport, stemness, HLA and antigen processing and presentation, and immune checkpoints were overexpressed in cancerous tissues, and many over-expressed HERV-related nearby genes were correlated with high HERV activation and poor survival. Twenty-three immune and stromal cells showed higher expression in non-cancerous than cancerous tissues, and seven were correlated with HERV activation. Small-molecule modulation of alternative splicing (AS) altered the expression of survival-related HERVs and their activation-related genes, as well as nearby genes.

CONCLUSION

Comprehensive and integrated approaches for evaluating HERV expression and their correlation with specific pathways have the potential to provide new companion diagnostics and therapeutic strategies for HCC.

Reproducible evaluation of transposable element detectors with McClintock 2 guides accurate inference of Ty insertion patterns in yeast

BACKGROUND

Many computational methods have been developed to detect non-reference transposable element (TE) insertions using short-read whole genome sequencing data. The diversity and complexity of such methods often present challenges to new users seeking to reproducibly install, execute, or evaluate multiple TE insertion detectors.

RESULTS

We previously developed the McClintock meta-pipeline to facilitate the installation, execution, and evaluation of six first-generation short-read TE detectors. Here, we report a completely re-implemented version of McClintock written in Python using Snakemake and Conda that improves its installation, error handling, speed, stability, and extensibility. McClintock 2 now includes 12 short-read TE detectors, auxiliary pre-processing and analysis modules, interactive HTML reports, and a simulation framework to reproducibly evaluate the accuracy of component TE detectors. When applied to the model microbial eukaryote Saccharomyces cerevisiae, we find substantial variation in the ability of McClintock 2 components to identify the precise locations of non-reference TE insertions, with RelocaTE2 showing the highest recall and precision in simulated data. We find that RelocaTE2, TEMP, TEMP2 and TEBreak provide consistent estimates of [Formula: see text]50 non-reference TE insertions per strain and that Ty2 has the highest number of non-reference TE insertions in a species-wide panel of [Formula: see text]1000 yeast genomes. Finally, we show that best-in-class predictors for yeast applied to resequencing data have sufficient resolution to reveal a dyad pattern of integration in nucleosome-bound regions upstream of yeast tRNA genes for Ty1, Ty2, and Ty4, allowing us to extend knowledge about fine-scale target preferences revealed previously for experimentally-induced Ty1 insertions to spontaneous insertions for other copia-superfamily retrotransposons in yeast.

CONCLUSION

McClintock ( https://github.com/bergmanlab/mcclintock/ ) provides a user-friendly pipeline for the identification of TEs in short-read WGS data using multiple TE detectors, which should benefit researchers studying TE insertion variation in a wide range of different organisms. Application of the improved McClintock system to simulated and empirical yeast genome data reveals best-in-class methods and novel biological insights for one of the most widely-studied model eukaryotes and provides a paradigm for evaluating and selecting non-reference TE detectors in other species.

Impact of retrotransposon protein L1 ORF1p expression on oncogenic pathways in hepatocellular carcinoma: the role of cytoplasmic PIN1 upregulation

BACKGROUND

Molecular characterisation of hepatocellular carcinoma (HCC) is central to the development of novel therapeutic strategies for the disease. We have previously demonstrated mutagenic consequences of Long-Interspersed Nuclear Element-1 (LINE1s/L1) retrotransposition. However, the role of L1 in HCC, besides somatic mutagenesis, is not well understood.

METHODS

We analysed L1 expression in the TCGA-HCC RNAseq dataset (n = 372) and explored potential relationships between L1 expression and clinical features. The findings were confirmed by immunohistochemical (IHC) analysis of an independent human HCC cohort (n = 48) and functional mechanisms explored using in vitro and in vivo model systems.

RESULTS

We observed positive associations between L1 and activated TGFβ-signalling, TP53 mutation, alpha-fetoprotein and tumour invasion. IHC confirmed a positive association between pSMAD3, a surrogate for TGFβ-signalling status, and L1 ORF1p (P < 0.0001, n = 32). Experimental modulation of L1 ORF1p levels revealed an influence of L1 ORF1p on key hepatocarcinogenesis-related pathways. Reduction in cell migration and invasive capacity was observed upon L1 ORF1 knockdown, both in vitro and in vivo. In particular, L1 ORF1p increased PIN1 cytoplasmic localisation. Blocking PIN1 activity abrogated L1 ORF1p-induced NF-κB-mediated inflammatory response genes while further activated TGFβ-signalling confirming differential alteration of PIN1 activity in cellular compartments by L1 ORF1p.

DISCUSSION

Our data demonstrate a causal link between L1 ORF1p and key oncogenic pathways mediated by PIN1, presenting a novel therapeutic avenue.

Reproducible evaluation of short-read transposable element detectors and species-wide data mining of insertion patterns in yeast

Background

Many computational methods have been developed to detect non-reference transposable element (TE) insertions using short-read whole genome sequencing data. The diversity and complexity of such methods often present challenges to new users seeking to reproducibly install, execute or evaluate multiple TE insertion detectors.

Results

We previously developed the McClintock meta-pipeline to facilitate the installation, execution, and evaluation of six first-generation short-read TE detectors. Here, we report a completely re-implemented version of McClintock written in Python using Snakemake and Conda that improves its installation, error handling, speed, stability, and extensibility. McClintock 2 now includes 12 short-read TE detectors, auxiliary pre-processing and analysis modules, interactive HTML reports, and a simulation framework to reproducibly evaluate the accuracy of component TE detectors. When applied to the model microbial eukaryote Saccharomyces cerevisiae , we find substantial variation in the ability of McClintock 2 components to identify the precise locations of non-reference TE insertions, with RelocaTE2 showing the highest recall and precision in simulated data. We find that RelocaTE2, TEMP, TEMP2 and TEBreak provide a consistent and biologically meaningful view of non-reference TE insertions in a species-wide panel of ∼ 1000 yeast genomes, as evaluated by coverage-based abundance estimates and expected patterns of tRNA promoter targeting. Finally, we show that best-in-class predictors for yeast have sufficient resolution to reveal a dyad pattern of integration in nucleosome-bound regions upstream of yeast tRNA genes for Ty1, Ty2, and Ty4, allowing us to extend knowledge aboutfine-scale target preferences first revealed experimentally for Ty1 to natural insertions and related copia -superfamily retrotransposons in yeast.

Conclusion

McClintock ( https://github.com/bergmanlab/mcclintock/ ) provides a user-friendly pipeline for the identification of TEs in short-read WGS data using multiple TE detectors, which should benefit researchers studying TE insertion variation in a wide range of different organisms. Application of the improved McClintock system to simulated and empirical yeast genome data reveals best-in-class methods and novel biological insights for one of the most widely-studied model eukaryotes and provides a paradigm for evaluating and selecting non-reference TE detectors for other species.

Retrotransposon insertions associated with risk of neurologic and psychiatric diseases

Transposable elements (TEs) are active in neuronal cells raising the question whether TE insertions contribute to risk of neuropsychiatric disease. While genome-wide association studies (GWAS) serve as a tool to discover genetic loci associated with neuropsychiatric diseases, unfortunately GWAS do not directly detect structural variants such as TEs. To examine the role of TEs in psychiatric and neurologic disease, we evaluated 17,000 polymorphic TEs and find 76 are in linkage disequilibrium with disease haplotypes (P < 10-6 ) defined by GWAS. From these 76 polymorphic TEs, we identify potentially causal candidates based on having insertions in genomic regions of regulatory chromatin and on having associations with altered gene expression in brain tissues. We show that lead candidate insertions have regulatory effects on gene expression in human neural stem cells altering the activity of a minimal promoter. Taken together, we identify 10 polymorphic TE insertions that are potential candidates on par with other variants for having a causal role in neurologic and psychiatric disorders.

Evolvability and emergence of tumor heterogeneity as a space-time function

The loss of control of cell proliferation, apoptosis regulation and contact inhibition leads to tumor development. While benign tumors are restricted to their primary space, i.e. where these tumors first originate, the metastatic tumors not only disseminate- facilitated by hypoxia-driven neovascularization- to distant secondary sites but also show substantial changes in metabolism, tissue architectures, gene expression profiles and immune phenotypes. All these alterations result in radio-, chemo- and immune-resistance rendering these metastatic tumor cells refractory to therapy. Since the beginning of the transformation, these factors- which influence each other- are incorporated to the developing and metastasizing tumor. As a result, the complexities in the heterogeneity of tumor progressively increase. This space-time function in the heterogeneity of tumors is generated by various conditions and factors at the genetic as well as microenvironmental levels, for example, endogenous retroviruses, methylation and epigenetic dysregulation that may be etiology-specific, cancer associated inflammation, remodeling of the extracellular matrix and mesenchymal cell shifted functions. On the one hand, these factors may cause de-differentiation of the tumor cells leading to cancer stem cells that contribute to radio-, chemo- and immune-resistance and recurrence of tumors. On the other hand, they may also enhance the heterogeneity under specific microenvironment-driven proliferation. In this editorial, we intend to underline the importance of heterogeneity in cancer progress, its evaluation and its use in correlation with the tumor evolution in a specific patient as a field of research for achieving precise patient-tailored treatments and amelioration of diagnostic (monitoring) tools and prognostic capacity.

Retrotransposon instability dominates the acquired mutation landscape of mouse induced pluripotent stem cells

Induced pluripotent stem cells (iPSCs) can in principle differentiate into any cell of the body, and have revolutionized biomedical research and regenerative medicine. Unlike their human counterparts, mouse iPSCs (miPSCs) are reported to silence transposable elements and prevent transposable element-mediated mutagenesis. Here we apply short-read or Oxford Nanopore Technologies long-read genome sequencing to 38 bulk miPSC lines reprogrammed from 10 parental cell types, and 18 single-cell miPSC clones. While single nucleotide variants and structural variants restricted to miPSCs are rare, we find 83 de novo transposable element insertions, including examples intronic to Brca1 and Dmd. LINE-1 retrotransposons are profoundly hypomethylated in miPSCs, beyond other transposable elements and the genome overall, and harbor alternative protein-coding gene promoters. We show that treatment with the LINE-1 inhibitor lamivudine does not hinder reprogramming and efficiently blocks endogenous retrotransposition, as detected by long-read genome sequencing. These experiments reveal the complete spectrum and potential significance of mutations acquired by miPSCs.

Characterisation of retrotransposon insertion polymorphisms in whole genome sequencing data from individuals with amyotrophic lateral sclerosis

The genetics of an individual is a crucial factor in understanding the risk of developing the neurodegenerative disease amyotrophic lateral sclerosis (ALS). There is still a large proportion of the heritability of ALS, particularly in sporadic cases, to be understood. Among others, active transposable elements drive inter-individual variability, and in humans long interspersed element 1 (LINE1, L1), Alu and SINE-VNTR-Alu (SVA) retrotransposons are a source of polymorphic insertions in the population. We undertook a pilot study to characterise the landscape of non-reference retrotransposon insertion polymorphisms (non-ref RIPs) in 15 control and 15 ALS individuals' whole genomes from Project MinE, an international project to identify potential genetic causes of ALS. The combination of two bioinformatics tools (mobile element locator tool (MELT) and TEBreak) identified on average 1250 Alu, 232 L1 and 77 SVA non-ref RIPs per genome across the 30 analysed. Further PCR validation of individual polymorphic retrotransposon insertions showed a similar level of accuracy for MELT and TEBreak. Our preliminary study did not identify a specific RIP or a significant difference in the total number of non-ref RIPs in ALS compared to control genomes. The use of multiple bioinformatic tools improved the accuracy of non-ref RIP detection and our study highlights the potential importance of studying these elements further in ALS.

Population Epigenetics: The Extent of DNA Methylation Variation in Wild Animal Populations

Population epigenetics explores the extent of epigenetic variation and its dynamics in natural populations encountering changing environmental conditions. In contrast to population genetics, the basic concepts of this field are still in their early stages, especially in animal populations. Epigenetic variation may play a crucial role in phenotypic plasticity and local adaptation as it can be affected by the environment, it is likely to have higher spontaneous mutation rate than nucleotide sequences do, and it may be inherited via non-mendelian processes. In this review, we aim to bring together natural animal population epigenetic studies to generate new insights into ecological epigenetics and its evolutionary implications. We first provide an overview of the extent of DNA methylation variation and its autonomy from genetic variation in wild animal population. Second, we discuss DNA methylation dynamics which create observed epigenetic population structures by including basic population genetics processes. Then, we highlight the relevance of DNA methylation variation as an evolutionary mechanism in the extended evolutionary synthesis. Finally, we suggest new research directions by highlighting gaps in the knowledge of the population epigenetics field. As for our results, DNA methylation diversity was found to reveal parameters that can be used to characterize natural animal populations. Some concepts of population genetics dynamics can be applied to explain the observed epigenetic structure in natural animal populations. The set of recent advancements in ecological epigenetics, especially in transgenerational epigenetic inheritance in wild animal population, might reshape the way ecologists generate predictive models of the capacity of organisms to adapt to changing environments.

Somatic retrotransposition in the developing rhesus macaque brain

The retrotransposon LINE-1 (L1) is central to the recent evolutionary history of the human genome and continues to drive genetic diversity and germline pathogenesis. However, the spatiotemporal extent and biological significance of somatic L1 activity are poorly defined and are virtually unexplored in other primates. From a single L1 lineage active at the divergence of apes and Old World monkeys, successive L1 subfamilies have emerged in each descendant primate germline. As revealed by case studies, the presently active human L1 subfamily can also mobilize during embryonic and brain development in vivo. It is unknown whether nonhuman primate L1s can similarly generate somatic insertions in the brain. Here we applied approximately 40× single-cell whole-genome sequencing (scWGS), as well as retrotransposon capture sequencing (RC-seq), to 20 hippocampal neurons from two rhesus macaques (Macaca mulatta). In one animal, we detected and PCR-validated a somatic L1 insertion that generated target site duplications, carried a short 5′ transduction, and was present in ∼7% of hippocampal neurons but absent from cerebellum and nonbrain tissues. The corresponding donor L1 allele was exceptionally mobile in vitro and was embedded in PRDM4, a gene expressed throughout development and in neural stem cells. Nanopore long-read methylome and RNA-seq transcriptome analyses indicated young retrotransposon subfamily activation in the early embryo, followed by repression in adult tissues. These data highlight endogenous macaque L1 retrotransposition potential, provide prototypical evidence of L1-mediated somatic mosaicism in a nonhuman primate, and allude to L1 mobility in the brain over the past 30 million years of human evolution.

Mutational signatures and processes in hepatobiliary cancers

The evolutionary history of hepatobiliary cancers is embedded in their genomes. By analysing their catalogue of somatic mutations and the DNA sequence context in which they occur, it is possible to infer the mechanisms underpinning tumorigenesis. These mutational signatures reflect the exogenous and endogenous origins of genetic damage as well as the capacity of hepatobiliary cells to repair and replicate DNA. Genomic analysis of thousands of patients with hepatobiliary cancers has highlighted the diversity of mutagenic processes active in these malignancies, highlighting a prominent source of the inter-cancer-type, inter-patient, intertumour and intratumoural heterogeneity that is observed clinically. However, a substantial proportion of mutational signatures detected in hepatocellular carcinoma and biliary tract cancer remain of unknown cause, emphasizing the important contribution of processes yet to be identified. Exploiting mutational signatures to retrospectively understand hepatobiliary carcinogenesis could advance preventative management of these aggressive tumours as well as potentially predict treatment response and guide the development of therapies targeting tumour evolution.

MER4 endogenous retrovirus correlated with better efficacy of anti-PD1/PD-L1 therapy in non-small cell lung cancer

Background

Endogenous retroviruses (ERVs) are highly expressed in various cancer types and are associated with increased innate immune response and better efficacy of antiprogrammed death-1/ligand-1 (anti-PD1/PD-L1)-directed immune checkpoint inhibitors (ICI) in preclinical models. However, their role in human non-small cell lung cancer (NSCLC) remains unknown.

Methods

We conducted a retrospective study of patients receiving ICI for advanced NSCLC in two independent cohorts. ERV expression was determined by RNA sequencing. The primary endpoint was progression-free survival (PFS) under ICI. The secondary endpoint was overall survival (OS) from ICI initiation. We studied expression of 6205 ERVs. Multivariate Cox regression model with lasso penalty was estimated on the training set to select ERVs significantly associated with survival. The predictive power of these ERVs was compared with that of previously described transcriptomic signatures.

Results

We studied two independent cohorts of 89 and 70 patients, used as training and validation sets. Clinicopathological characteristics included 75% of patients with non-squamous NSCLC. We selected four ERVs significantly associated with PFS. Only high MER4 ERV was associated with better PFS and OS in both cohorts. From a biological point of view, high MER4 expression is associated with higher infiltration of eosinophils and inflammatory gene signatures, while low MER4 expression is associated with enrichment in metabolism and proliferation signatures. Adding MER4 to previously described transcriptomic signatures of response to ICI improved their predictive power.

Conclusions

MER4 ERV expression is useful to stratify risk and predict PFS and OS in patients treated with ICI for NSCLC. It also improves the predictive power of other known transcriptomic signatures.

Identification of differentially expressed genomic repeats in primary hepatocellular carcinoma and their potential links to biological processes and survival

Hepatocellular carcinoma (HCC) is one of the deadliest cancers. Research on HCC so far primarily focused on genes and provided limited information on genomic repeats, which constitute more than half of the human genome and contribute to genomic stability. In line with this, repeat dysregulation was significantly shown to be pathological in various cancers and other diseases. In this study, we aimed to determine the full repeat expression profile of HCC for the first time. We utilised two independent RNA-seq datasets obtained from primary HCC tumours with matched normal tissues of 20 and 17 HCC patients, respectively. We quantified repeat expressions and analysed their differential expression. We also identified repeats that are cooperatively expressed with genes by constructing a gene coexpression network. Our results indicated that HCC tumours in both datasets harbour 24 differentially expressed repeats and even more elements were coexpressed with genes involved in various metabolic pathways. We discovered that two L1 elements (L1M3b, L1M3de) were downregulated and a handful of HERV subfamily repeats (HERV-Fc1-int, HERV3-int, HERVE_a-int, HERVK11D-int, HERVK14C-int, HERVL18-int) were upregulated with the exception of HERV1_LTRc, which was downregulated. Various LTR elements (LTR32, LTR9, LTR4, LTR52-int, LTR70) and MER elements (MER11C, MER11D, MER57C1, MER9a1, MER74C) were implicated along with few other subtypes including Charlie12, MLT2A2, Tigger15a, Tigger 17b. The only satellite repeat differentially expressed in both datasets was GSATII, whose expression was upregulated in 33 (>90%) out of 37 patients. Notably, GSATII expression correlated with HCC survival genes. Elements discovered here promise future studies to be considered for biomarker and HCC therapy research. The coexpression pattern of the GSATII satellite with HCC survival genes and the fact that it has been upregulated in the vast majority of patients make this repeat particularly stand out for HCC.

HCV Activates Somatic L1 Retrotransposition-A Potential Hepatocarcinogenesis Pathway

Hepatitis C virus (HCV) is a common cause of hepatocellular carcinoma (HCC). The activation and mutagenic consequences of L1 retrotransposons in virus-associated-HCC have been documented. However, the direct influence of HCV upon L1 elements is unclear, and is the focus of the present study. L1 transcript expression was evaluated in a publicly available liver tissue RNA-seq dataset from patients with chronic HCV hepatitis (CHC), as well as healthy controls. L1 transcript expression was significantly higher in CHC than in controls. L1orf1p (a L1 encoded protein) expression was observed in six out of 11 CHC livers by immunohistochemistry. To evaluate the influence of HCV on retrotransposition efficiency, in vitro engineered-L1 retrotransposition assays were employed in Huh7 cells in the presence and absence of an HCV replicon. An increased retrotransposition rate was observed in the presence of replicating HCV RNA, and persisted in cells after viral clearance due to sofosbuvir (PSI7977) treatment. Increased retrotransposition could be due to dysregulation of the DNA-damage repair response, including homologous recombination, due to HCV infection. Altogether these data suggest that L1 expression can be activated before oncogenic transformation in CHC patients, with HCV-upregulated retrotransposition potentially contributing to HCC genomic instability and a risk of transformation that persists post-viral clearance.

No evidence of human genome integration of SARS-CoV-2 found by long-read DNA sequencing

A recent study proposed that severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) hijacks the LINE-1 (L1) retrotransposition machinery to integrate into the DNA of infected cells. If confirmed, this finding could have significant clinical implications. Here, we apply deep (>50×) long-read Oxford Nanopore Technologies (ONT) sequencing to HEK293T cells infected with SARS-CoV-2 and do not find the virus integrated into the genome. By examining ONT data from separate HEK293T cultivars, we completely resolve 78 L1 insertions arising in vitro in the absence of L1 overexpression systems. ONT sequencing applied to hepatitis B virus (HBV)-positive liver cancer tissues located a single HBV insertion. These experiments demonstrate reliable resolution of retrotransposon and exogenous virus insertions by ONT sequencing. That we find no evidence of SARS-CoV-2 integration suggests that such events are, at most, extremely rare in vivo and therefore are unlikely to drive oncogenesis or explain post-recovery detection of the virus.

Value of methylation markers in colorectal cancer (Review)

Colorectal cancer (CRC) is a multifactorial and multistage process that occurs due to both genetic and epigenetic variations in normal epithelial cells. Analysis of the CRC epigenome has revealed that almost all CRC types have a large number of abnormally methylated genes. Hypermethylation of cell‑free DNA from CRC in the blood or stool is considered as a potential non‑invasive cancer biomarker, and various methylation markers have shown high sensitivity and specificity. The aim of the present review was to examine potential methylation markers in CRC that have been used or are expected to be used in the clinical setting, focusing on their screening, predictive, prognostic and therapeutic roles in CRC.

Paternal MTHFR deficiency leads to hypomethylation of young retrotransposons and reproductive decline across two successive generations

5,10-Methylenetetrahydrofolate reductase (MTHFR) is a crucial enzyme in the folate metabolic pathway with a key role in generating methyl groups. As MTHFR deficiency impacts male fertility and sperm DNA methylation, there is the potential for epimutations to be passed to the next generation. Here, we assessed whether the impact of MTHFR deficiency on testis morphology and sperm DNA methylation is exacerbated across generations in mouse. Although MTHFR deficiency in F1 fathers has only minor effects on sperm counts and testis weights and histology, F2 generation sons show further deterioration in reproductive parameters. Extensive loss of DNA methylation is observed in both F1 and F2 sperm, with >80% of sites shared between generations, suggestive of regions consistently susceptible to MTHFR deficiency. These regions are generally methylated during late embryonic germ cell development and are enriched in young retrotransposons. As retrotransposons are resistant to reprogramming of DNA methylation in embryonic germ cells, their hypomethylated state in the sperm of F1 males could contribute to the worsening reproductive phenotype observed in F2 MTHFR-deficient males, compatible with the intergenerational passage of epimutations.

HBV Integration Induces Complex Interactions between Host and Viral Genomic Functions at the Insertion Site

Hepatitis B virus (HBV), one of the well-known DNA oncogenic viruses, is the leading cause of hepatocellular carcinoma (HCC). In infected hepatocytes, HBV DNA can be integrated into the host genome through an insertional mutagenesis process inducing tumorigenesis. Dissection of the genomic features surrounding integration sites will deepen our understanding of mechanisms underlying integration. Moreover, the quantity and biological activity of integration sites may reflect the DNA damage within affected cells or the potential survival benefits they may confer. The well-known human genomic features include repeat elements, particular regions (such as telomeres), and frequently interrupted genes (e.g., telomerase reverse transcriptase [i.e. TERT], lysine methyltransferase 2B [i.e. KMT2B], cyclin E1 [CCNE1], and cyclin A2 [CCNA2]). Consequently, distinct genomic features within diverse integrations differentiate their biological functions. Meanwhile, accumulating evidence has shown that viral proteins produced by integrants may cause cell damage even after the suppression of HBV replication. The integration-derived gene products can also serve as tumor markers, promoting the development of novel therapeutic strategies for HCC. Viral integrants can be single copy or multiple copies of different fragments with complicated rearrangement, which warrants elucidation of the whole viral integrant arrangement in future studies. All of these considerations underlie an urgent need to develop novel methodology and technology for sequence characterization and function evaluation of integration events in chronic hepatitis B-associated disease progression by monitoring both host genomic features and viral integrants. This endeavor may also serve as a promising solution for evaluating the risk of tumorigenesis and as a companion diagnostic for designing therapeutic strategies targeting integration-related disease complications.

DNA methylation profile of liver of mice conceived by in vitro fertilization

Offspring generated by in vitro fertilization (IVF) are believed to be healthy but display a possible predisposition to chronic diseases, like hypertension and glucose intolerance. Since epigenetic changes are believed to underlie such phenotype, this study aimed at describing global DNA methylation changes in the liver of adult mice generated by natural mating (FB group) or by IVF. Embryos were generated by IVF or natural mating. At 30 weeks of age, mice were sacrificed. The liver was removed, and global DNA methylation was assessed using whole-genome bisulfite sequencing (WGBS). Genomic Regions for Enrichment Analysis Tool (GREAT) and G:Profilerβ were used to identify differentially methylated regions (DMRs) and for functional enrichment analysis. Overrepresented gene ontology terms were summarized with REVIGO, while canonical pathways (CPs) were identified with Ingenuity® Pathway Analysis. Overall, 2692 DMRs (4.91%) were different between the groups. The majority of DMRs (84.92%) were hypomethylated in the IVF group. Surprisingly, only 0.16% of CpG islands were differentially methylated and only a few DMRs were located on known gene promoters (n = 283) or enhancers (n = 190). Notably, the long-interspersed element (LINE), short-interspersed element (SINE), and long terminal repeat (LTR1) transposable elements showed reduced methylation (P < 0.05) in IVF livers. Cellular metabolic process, hepatic fibrosis, and insulin receptor signaling were some of the principal biological processes and CPs modified by IVF. In summary, IVF modifies the DNA methylation signature in the adult liver, resulting in hypomethylation of genes involved in metabolism and gene transcription regulation. These findings may shed light on the mechanisms underlying the developmental origin of health and disease.

LINE-1 transcription in round spermatids is associated with accretion of 5-carboxylcytosine in their open reading frames

Chromatin of male and female gametes undergoes a number of reprogramming events during the transition from germ cell to embryonic developmental programs. Although the rearrangement of DNA methylation patterns occurring in the zygote has been extensively characterized, little is known about the dynamics of DNA modifications during spermatid maturation. Here, we demonstrate that the dynamics of 5-carboxylcytosine (5caC) correlate with active transcription of LINE-1 retroelements during murine spermiogenesis. We show that the open reading frames of active and evolutionary young LINE-1s are 5caC-enriched in round spermatids and 5caC is eliminated from LINE-1s and spermiogenesis-specific genes during spermatid maturation, being simultaneously retained at promoters and introns of developmental genes. Our results reveal an association of 5caC with activity of LINE-1 retrotransposons suggesting a potential direct role for this DNA modification in fine regulation of their transcription.

Robust expression of LINE-1 retrotransposon encoded proteins in oral squamous cell carcinoma

BACKGROUND

Oral Squamous Cell Carcinoma (OSCC) results from a series of genetic alteration in squamous cells. This particular type of cancer considers one of the most aggressive malignancies to control because of its frequent local invasions to the regional lymph node. Although several biomarkers have been reported, the key marker used to predict the behavior of the disease is largely unknown. Here we report Long INterpersed Element-1 (LINE1 or L1) retrotransposon activity in post-operative oral cancer samples. L1 is the only active retrotransposon occupying around 17% of the human genome with an estimated 500,000 copies. An active L1 encodes two proteins (L1ORF1p and L1ORF2p); both of which are critical in the process of retrotransposition. Several studies report that the L1 retrotransposon is highly active in many cancers. L1 activity is generally determined by assaying L1ORF1p because of its high expression and availability of the antibody. However, due to its lower expression and unavailability of a robust antibody, detection of L1ORF2p has been limited. L1ORF2p is the crucial protein in the process of retrotransposition as it provides endonuclease and reverse transcriptase (RT) activity.

METHODS

Immunohistochemistry and Western blotting were performed on the post-operative oral cancer samples and murine tissues.

RESULTS

Using in house novel antibodies against both the L1 proteins (L1ORF1p and L1ORF2p), we found L1 retrotransposon is extremely active in post-operative oral cancer tissues. Here, we report a novel human L1ORF2p antibody generated using an 80-amino-acid stretch from the RT domain, which is highly conserved among different species. The antibody detects significant L1ORF2p expression in human oral squamous cell carcinoma (OSCC) samples and murine germ tissues.

CONCLUSIONS

We report exceptionally high L1ORF1p and L1ORF2p expression in post-operative oral cancer samples. The novel L1ORF2p antibody reported in this study will serve as a useful tool to understand why L1 activity is deregulated in OSCC and how it contributes to the progression of this particular cancer. Cross-species reactivity of L1ORF2p antibody due to the conserved epitope will be useful to study the retrotransposon biology in mice and rat germ tissues.

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.

Unraveling the features of somatic transposition in the Drosophila intestine

Transposable elements (TEs) play a significant role in evolution, contributing to genetic variation. However, TE mobilization in somatic cells is not well understood. Here, we address the prevalence of transposition in a somatic tissue, exploiting the Drosophila midgut as a model. Using whole-genome sequencing of in vivo clonally expanded gut tissue, we have mapped hundreds of high-confidence somatic TE integration sites genome-wide. We show that somatic retrotransposon insertions are associated with inactivation of the tumor suppressor Notch, likely contributing to neoplasia formation. Moreover, applying Oxford Nanopore long-read sequencing technology we provide evidence for tissue-specific differences in retrotransposition. Comparing somatic TE insertional activity with transcriptomic and small RNA sequencing data, we demonstrate that transposon mobility cannot be simply predicted by whole tissue TE expression levels or by small RNA pathway activity. Finally, we reveal that somatic TE insertions in the adult fly intestine are enriched in genic regions and in transcriptionally active chromatin. Together, our findings provide clear evidence of ongoing somatic transposition in Drosophila and delineate previously unknown features underlying somatic TE mobility in vivo.

Intercellular viral spread and intracellular transposition of Drosophila gypsy

It has become increasingly clear that retrotransposons (RTEs) are more widely expressed in somatic tissues than previously appreciated. RTE expression has been implicated in a myriad of biological processes ranging from normal development and aging, to age related diseases such as cancer and neurodegeneration. Long Terminal Repeat (LTR)-RTEs are evolutionary ancestors to, and share many features with, exogenous retroviruses. In fact, many organisms contain endogenous retroviruses (ERVs) derived from exogenous retroviruses that integrated into the germ line. These ERVs are inherited in Mendelian fashion like RTEs, and some retain the ability to transmit between cells like viruses, while others develop the ability to act as RTEs. The process of evolutionary transition between LTR-RTE and retroviruses is thought to involve multiple steps by which the element loses or gains the ability to transmit copies between cells versus the ability to replicate intracellularly. But, typically, these two modes of transmission are incompatible because they require assembly in different sub-cellular compartments. Like murine IAP/IAP-E elements, the gypsy family of retroelements in arthropods appear to sit along this evolutionary transition. Indeed, there is some evidence that gypsy may exhibit retroviral properties. Given that gypsy elements have been found to actively mobilize in neurons and glial cells during normal aging and in models of neurodegeneration, this raises the question of whether gypsy replication in somatic cells occurs via intracellular retrotransposition, intercellular viral spread, or some combination of the two. These modes of replication in somatic tissues would have quite different biological implications. Here, we demonstrate that Drosophila gypsy is capable of both cell-associated and cell-free viral transmission between cultured S2 cells of somatic origin. Further, we demonstrate that the ability of gypsy to move between cells is dependent upon a functional copy of its viral envelope protein. This argues that the gypsy element has transitioned from an RTE into a functional endogenous retrovirus with the acquisition of its envelope gene. On the other hand, we also find that intracellular retrotransposition of the same genomic copy of gypsy can occur in the absence of the Env protein. Thus, gypsy exhibits both intracellular retrotransposition and intercellular viral transmission as modes of replicating its genome.

uhrf1 and dnmt1 Loss Induces an Immune Response in Zebrafish Livers Due to Viral Mimicry by Transposable Elements

Activation of transposable elements (TEs) can cause cellular damage. Cytoplasmic nucleic acid sensing pathways evolved to detect pathogens, but can also serve to cull cells with inappropriate TE activation as TEs can be viral mimetics. Epigenetic silencing of TEs is mediated in part by DNA methylation, but it is not clear if TE activation or the immune system contribute to the cellular damage caused by loss of DNA methylation. Here, we provide mechanistic insight into the observation of an activated interferon response in the liver of zebrafish larvae with deletion in critical components of the DNA methylation machinery, uhrf1 and dnmt1. We focus on dissecting the relationship between DNA methylation, TE activation and induction of an immune response through cytoplasmic DNA and double stranded RNA sensing pathways and identify tnfa as a mediator of cell death in the liver of these mutants. Integrated RNAseq and methylome analysis identified LTR transposons as the most upregulated in these mutants and also the most methylated in control larvae, indicating a direct role of DNA methylation in suppressing this TE subclass. RNAseq analysis from these same samples revealed expression signatures of a type-I interferon response and of tnfa activation, mimicking the pattern of gene expression in virally infected cells. CRISPR/Cas9 mediated depletion of the cellular antiviral sensors sting and mavs reduced expression of interferon response genes and tnfa depletion dramatically reduced cell death in uhrf1 mutant livers. This suggests that the antiviral response induced by DNA hypomethylation and TE activation in the liver is mediated by the signaling pathways activated by both cytoplasmic double stranded RNA and DNA and that tnfa mediates cell death as a potential mechanism to eliminate these damaged cells.

SINE-B1 Distribution and Chromosome Rearrangements in the South American Proechimys gr. goeldii (Echimyidae, Rodentia)

Proechimys species are remarkable for their extensive chromosome rearrangements, representing a good model to understand genome evolution. Herein, we cytogenetically analyzed 3 different cytotypes of Proechimys gr. goeldii to assess their evolutionary relationship. We also mapped the transposable element SINE-B1 on the chromosomes of P. gr. goeldii in order to investigate its distribution among individuals and evaluate its possible contribution to karyotype remodeling in this species. SINE-B1 showed a dispersed distribution along chromosome arms and was also detected at the pericentromeric regions of some chromosomes, including pair 1 and the sex chromosomes, which are involved in chromosome rearrangements. In addition, we describe a new cytotype for P. gr. goeldii, reinforcing the significant role of gross chromosomal rearrangements during the evolution of the genus. The results of FISH with SINE-B1 suggest that this issue should be more deeply investigated for a better understanding of its role in the mechanisms involved in the wide variety of Proechimys karyotypes.

SeqURE - a new copy-capture based method for sequencing of unknown Retroposition events

Background

Retroelements (REs) occupy a significant part of all eukaryotic genomes including humans. The majority of retroelements in the human genome are inactive and unable to retrotranspose. Dozens of active copies are repressed in most normal tissues by various cellular mechanisms. These copies can become active in normal germline and brain tissues or in cancer, leading to new retroposition events. The consequences of such events and their role in normal cell functioning and carcinogenesis are not yet fully understood. If new insertions occur in a small portion of cells they can be found only with the use of specific methods based on RE enrichment and high-throughput sequencing. The downside of the high sensitivity of such methods is the presence of various artifacts imitating real insertions, which in many cases cannot be validated due to lack of the initial template DNA. For this reason, adequate assessment of rare (< 1%) subclonal cancer specific RE insertions is complicated.

Results

Here we describe a new copy-capture technique which we implemented in a method called SeqURE for Sequencing Unknown of Retroposition Events that allows for efficient and reliable identification of new genomic RE insertions. The method is based on the capture of copies of target molecules (copy-capture), selective amplification and sequencing of genomic regions adjacent to active RE insertions from both sides. Importantly, the template genomic DNA remains intact and can be used for validation experiments. In addition, we applied a novel system for testing method sensitivity and precisely showed the ability of the developed method to reliably detect insertions present in 1 out of 100 cells and a substantial portion of insertions present in 1 out of 1000 cells. Using advantages of the method we showed the absence of somatic Alu insertions in colorectal cancer samples bearing tumor-specific L1HS insertions.

Conclusions

This study presents the first description and implementation of the copy-capture technique and provides the first methodological basis for the quantitative assessment of RE insertions present in a small portion of cells.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13100-020-00228-6.

Gestational arsenic exposure induces site-specific DNA hypomethylation in active retrotransposon subfamilies in offspring sperm in mice

Background

Environmental impacts on a fetus can disrupt germ cell development leading to epimutations in mature germ cells. Paternal inheritance of adverse health effects through sperm epigenomes, including DNA methylomes, has been recognized in human and animal studies. However, the impacts of gestational exposure to a variety of environmental factors on the germ cell epigenomes are not fully investigated. Arsenic, a naturally occurring contaminant, is one of the most concerning environmental chemicals, that is causing serious health problems, including an increase in cancer, in highly contaminated areas worldwide. We previously showed that gestational arsenic exposure of pregnant C3H mice paternally induces hepatic tumor increase in the second generation (F2). In the present study, we have investigated the F1 sperm DNA methylomes genome-widely by one-base resolution analysis using a reduced representation bisulfite sequencing (RRBS) method.

Results

We have clarified that gestational arsenic exposure increases hypomethylated cytosines in all the chromosomes and they are significantly overrepresented in the retrotransposon LINEs and LTRs, predominantly in the intergenic regions. Closer analyses of detailed annotated DNA sequences showed that hypomethylated cytosines are especially accumulated in the promoter regions of the active full-length L1MdA subfamily in LINEs, and 5′LTRs of the active IAPE subfamily in LTRs. This is the first report that has identified the specific positions of methylomes altered in the retrotransposon elements by environmental exposure, by genome-wide methylome analysis.

Conclusion

Lowered DNA methylation potentially enhances L1MdA retrotransposition and cryptic promoter activity of 5′LTR for coding genes and non-coding RNAs. The present study has illuminated the environmental impacts on sperm DNA methylome establishment that can lead to augmented retrotransposon activities in germ cells and can cause harmful effects in the following generation.

Nanopore Sequencing Enables Comprehensive Transposable Element Epigenomic Profiling

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

The tumor suppressor microRNA let-7 inhibits human LINE-1 retrotransposition

Nearly half of the human genome is made of transposable elements (TEs) whose activity continues to impact its structure and function. Among them, Long INterspersed Element class 1 (LINE-1 or L1) elements are the only autonomously active TEs in humans. L1s are expressed and mobilized in different cancers, generating mutagenic insertions that could affect tumor malignancy. Tumor suppressor microRNAs are ∼22nt RNAs that post-transcriptionally regulate oncogene expression and are frequently downregulated in cancer. Here we explore whether they also influence L1 mobilization. We show that downregulation of let-7 correlates with accumulation of L1 insertions in human lung cancer. Furthermore, we demonstrate that let-7 binds to the L1 mRNA and impairs the translation of the second L1-encoded protein, ORF2p, reducing its mobilization. Overall, our data reveals that let-7, one of the most relevant microRNAs, maintains somatic genome integrity by restricting L1 retrotransposition.

Human Long INterspersed Element class 1 (LINE-1) elements are expressed and mobilized in many types of cancer, contributing to malignancy. Here the authors show that the tumor suppressor microRNA let-7 targets the LINE-1 mRNA and reduces LINE-1 mobilization.

Gestational arsenic exposure and paternal intergenerational epigenetic inheritance

A growing body of evidence has shown that gestational exposure to environmental factors such as imbalanced diet, environmental chemicals, and stress can lead to late-onset health effects in offspring and that some of these effects are heritable by the next generation and subsequent generations. Furthermore, altered epigenetic modifications in DNA methylation, histone modifications and small RNAs in a single sperm genome have been shown to transmit disease phenotypes acquired from the environment to later generations. Recently, our group found that gestational exposure of F0 pregnant dams to an inorganic arsenic, sodium arsenite, increases the incidence of hepatic tumors in male F2 mice, and the effects are paternally transmitted to the F2. Here, we first overview the epigenetic changes involved in paternal intergenerational and transgenerational inheritance caused by exposure to environmental factors. Then, we discuss our recent studies regarding paternal inheritance of the tumor-augmenting effects in F2 mice by gestational arsenite exposure, in which we investigated alterations of DNA methylation status in F2 tumors and causative F1 sperm. We also discuss the possible targets of the F2 effects. Finally, we discuss future perspectives on the studies that are needed to fully understand the health effects of arsenic exposure.

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

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

Transposon expression in the Drosophila brain is driven by neighboring genes and diversifies the neural transcriptome

Somatic transposon expression in neural tissue is commonly considered as a measure of mobilization and has therefore been linked to neuropathology and organismal individuality. We combined genome sequencing data with single-cell mRNA sequencing of the same inbred fly strain to map transposon expression in the Drosophila midbrain and found that transposon expression patterns are highly stereotyped. Every detected transposon is resident in at least one cellular gene with a matching expression pattern. Bulk RNA sequencing from fly heads of the same strain revealed that coexpression is a physical link in the form of abundant chimeric transposon-gene mRNAs. We identified 264 genes where transposons introduce cryptic splice sites into the nascent transcript and thereby significantly expand the neural transcript repertoire. Some genes exclusively produce chimeric mRNAs with transposon sequence; on average, 11.6% of the mRNAs produced from a given gene are chimeric. Conversely, most transposon-containing transcripts are chimeric, which suggests that somatic expression of these transposons is largely driven by cellular genes. We propose that chimeric mRNAs produced by alternative splicing into polymorphic transposons, rather than transposon mobilization, may contribute to functional differences between individual cells and animals.

Promoter methylation, transcription, and retrotransposition of LINE-1 in colorectal adenomas and adenocarcinomas

BACKGROUND

The methylation of the CpG islands of the LINE-1 promoter is a tight control mechanism on the function of mobile elements. However, simultaneous quantification of promoter methylation and transcription of LINE-1 has not been performed in progressive stages of colorectal cancer. In addition, the insertion of mobile elements in the genome of advanced adenoma stage, a precancerous stage before colorectal carcinoma has not been emphasized. In this study, we quantify promoter methylation and transcripts of LINE-1 in three stages of colorectal non-advanced adenoma, advanced adenoma, and adenocarcinoma. In addition, we analyze the insertion of LINE-1, Alu, and SVA elements in the genome of patient tumors with colorectal advanced adenomas.

METHODS

LINE-1 hypomethylation status was evaluated by absolute quantitative analysis of methylated alleles (AQAMA) assay. To quantify the level of transcripts for LINE-1, quantitative RT-PCR was performed. To find mobile element insertions, the advanced adenoma tissue samples were subjected to whole genome sequencing and MELT analysis.

RESULTS

We found that the LINE-1 promoter methylation in advanced adenoma and adenocarcinoma was significantly lower than that in non-advanced adenomas. Accordingly, the copy number of LINE-1 transcripts in advanced adenoma was significantly higher than that in non-advanced adenomas, and in adenocarcinomas was significantly higher than that in the advanced adenomas. Whole-genome sequencing analysis of colorectal advanced adenomas revealed that at this stage polymorphic insertions of LINE-1, Alu, and SVA comprise approximately 16%, 51%, and 74% of total insertions, respectively.

CONCLUSIONS

Our correlative analysis showing a decreased methylation of LINE-1 promoter accompanied by the higher level of LINE-1 transcription, and polymorphic genomic insertions in advanced adenoma, suggests that the early and advanced polyp stages may host very important pathogenic processes concluding to cancer.

From genome sequencing to the discovery of potential biomarkers in liver disease

Chronic liver disease progresses through several stages, fatty liver, steatohepatitis, cirrhosis, and eventually, it leads to hepatocellular carcinoma (HCC) over a long period of time. Since a large proportion of patients with HCC are accompanied by cirrhosis, it is considered to be an important factor in the diagnosis of liver cancer. This is because cirrhosis leads to an irreversible harmful effect, but the early stages of chronic liver disease could be reversed to a healthy state. Therefore, the discovery of biomarkers that could identify the early stages of chronic liver disease is important to prevent serious liver damage. Biomarker discovery at liver cancer and cirrhosis has enhanced the development of sequencing technology. Next generation sequencing (NGS) is one of the representative technical innovations in the biological field in the recent decades and it is the most important thing to design for research on what type of sequencing methods are suitable and how to handle the analysis steps for data integration. In this review, we comprehensively summarized NGS techniques for identifying genome, transcriptome, DNA methylome and 3D/4D chromatin structure, and introduced framework of processing data set and integrating multi-omics data for uncovering biomarkers.

KRAB-zinc finger protein gene expansion in response to active retrotransposons in the murine lineage

The Krüppel-associated box zinc finger protein (KRAB-ZFP) family diversified in mammals. The majority of human KRAB-ZFPs bind transposable elements (TEs), however, since most TEs are inactive in humans it is unclear whether KRAB-ZFPs emerged to suppress TEs. We demonstrate that many recently emerged murine KRAB-ZFPs also bind to TEs, including the active ETn, IAP, and L1 families. Using a CRISPR/Cas9-based engineering approach, we genetically deleted five large clusters of KRAB-ZFPs and demonstrate that target TEs are de-repressed, unleashing TE-encoded enhancers. Homozygous knockout mice lacking one of two KRAB-ZFP gene clusters on chromosome 2 and chromosome 4 were nonetheless viable. In pedigrees of chromosome 4 cluster KRAB-ZFP mutants, we identified numerous novel ETn insertions with a modest increase in mutants. Our data strongly support the current model that recent waves of retrotransposon activity drove the expansion of KRAB-ZFP genes in mice and that many KRAB-ZFPs play a redundant role restricting TE activity.

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

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

LINE-1 ORF2p expression is nearly imperceptible in human cancers

BACKGROUND

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

RESULTS

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

CONCLUSIONS

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

L1 Retrotransposon Heterogeneity in Ovarian Tumor Cell Evolution

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

Defective modulation of LINE-1 retrotransposition by cancer-associated SAMHD1 mutants

Long interspersed nuclear elements (LINE-1) is now considered as the only active autonomous mobile DNA in humans, LINE-1 retrotransposition activities are associated with and fluctuate during cancer initiation and progression; however, the mechanism underlying the increased LINE-1 activity in cancer is poorly understood. SAMHD1 has been reported to be a potent inhibitor of LINE-1 retrotransposition, and SAMHD1 mutations are frequently associated with cancer development. To gain insights on whether cancer-related SAMHD1 mutants affect LINE-1 activity, we explored the biochemical and cellular properties of some human mutants known correlate with the development of cancer. Most of the tested SAMHD1 cancer-related mutations were defective in LINE-1 inhibition. Interestingly we also found that SAMHD1 mutant K288T was defective for dNTPase activity but showed potent activity against LINE-1 retrotransposition. These findings suggest that LINE-1 inhibition does not depend solely on the dNTPase activity of SAMHD1. In contrast, SAMHD1's ability to inhibit ORF2p-mediated LINE-1 RNP reverse transcription was correlated with SAMHD1-mediated LINE-1 inhibition. Together, our data could also facilitate the deeper understanding for the inhibition of endogenous LINE-1 elements by SAMHD1.

Host-transposon interactions: conflict, cooperation, and cooption

Transposable elements (TEs) are mobile DNA sequences that colonize genomes and threaten genome integrity. As a result, several mechanisms appear to have emerged during eukaryotic evolution to suppress TE activity. However, TEs are ubiquitous and account for a prominent fraction of most eukaryotic genomes. We argue that the evolutionary success of TEs cannot be explained solely by evasion from host control mechanisms. Rather, some TEs have evolved commensal and even mutualistic strategies that mitigate the cost of their propagation. These coevolutionary processes promote the emergence of complex cellular activities, which in turn pave the way for cooption of TE sequences for organismal function.