Nanopore Sequencing Enables Comprehensive Transposable Element Epigenomic Profiling

Research progress of LINE-1 in the diagnosis, prognosis, and treatment of gynecologic tumors

The retrotransposon known as long interspersed nuclear element-1 (LINE-1), which is currently the sole autonomously mobile transposon in the human genome, can result in insertional mutations, chromosomal rearrangements, and genomic instability. In recent years, numerous studies have shown that LINE-1 is involved in the development of various diseases and also plays an important role in the immune regulation of the organism. The expression of LINE-1 in gynecologic tumors suggests that it is expected to be an independent indicator for early diagnosis and prognosis, and also, as a therapeutic target, LINE-1 is closely associated with gynecologic tumor prognosis. This article discusses the function of LINE-1 in the diagnosis, treatment, and prognosis of ovarian, cervical, and endometrial malignancies, as well as other gynecologic malignancies. It offers fresh perspectives on the early detection of tumors and the creation of novel anti-tumor medications.

LINE-1 repression in Epstein-Barr virus-associated gastric cancer through viral-host genome interaction

Long INterspersed Element 1 (LINE-1 or L1) acts as a major remodeling force in genome regulation and evolution. Accumulating evidence shows that virus infection impacts L1 expression, potentially impacting host antiviral response and diseases. The underlying regulation mechanism is unclear. Epstein-Barr virus (EBV), a double-stranded DNA virus linked to B-cell and epithelial malignancies, is known to have viral-host genome interaction, resulting in transcriptional rewiring in EBV-associated gastric cancer (EBVaGC). By analyzing publicly available datasets from the Gene Expression Omnibus (GEO), we found that EBVaGC has L1 transcriptional repression compared with EBV-negative gastric cancer (EBVnGC). More specifically, retrotransposition-associated young and full-length L1s (FL-L1s) were among the most repressed L1s. Epigenetic alterations, especially increased H3K9me3, were observed on FL-L1s. H3K9me3 deposition was potentially attributed to increased TASOR expression, a key component of the human silencing hub (HUSH) complex for H3K9 trimethylation. The 4C- and HiC-seq data indicated that the viral DNA interacted in the proximity of the TASOR enhancer, strengthening the loop formation between the TASOR enhancer and its promoter. These results indicated that EBV infection is associated with increased H3K9me3 deposition, leading to L1 repression. This study uncovers a regulation mechanism of L1 expression by chromatin topology remodeling associated with viral-host genome interaction in EBVaGC.

Chromosome-level, nanopore-only genome and allele-specific DNA methylation of Pallas's cat, Otocolobus manul

Pallas's cat, or the manul cat (Otocolobus manul), is a small felid native to the grasslands and steppes of central Asia. Population strongholds in Mongolia and China face growing challenges from climate change, habitat fragmentation, poaching, and other sources. These threats, combined with O. manul's zoo collection popularity and value in evolutionary biology, necessitate improvement of species genomic resources. We used standalone nanopore sequencing to assemble a 2.5 Gb, 61-contig nuclear assembly and 17097 bp mitogenome for O. manul. The primary nuclear assembly had 56× sequencing coverage, a contig N50 of 118 Mb, and a 94.7% BUSCO completeness score for Carnivora-specific genes. High genome collinearity within Felidae permitted alignment-based scaffolding onto the fishing cat (Prionailurus viverrinus) reference genome. Manul contigs spanned all 19 felid chromosomes with an inferred total gap length of less than 400 kilobases. Modified basecalling and variant phasing produced an alternate pseudohaplotype assembly and allele-specific DNA methylation calls; 61 differentially methylated regions were identified between haplotypes. Nearest features included classical imprinted genes, non-coding RNAs, and putative novel imprinted loci. The assembled mitogenome successfully resolved existing discordance between Felinae nuclear and mtDNA phylogenies. All assembly drafts were generated from 158 Gb of sequence using seven minION flow cells.

Rapid evolution of piRNA clusters in the Drosophila melanogaster ovary

Animal genomes are parasitized by a horde of transposable elements (TEs) whose mutagenic activity can have catastrophic consequences. The piRNA pathway is a conserved mechanism to repress TE activity in the germline via a specialized class of small RNAs associated with effector Piwi proteins called piwi-associated RNAs (piRNAs). piRNAs are produced from discrete genomic regions called piRNA clusters (piCs). While piCs are generally enriched for TE sequences and the molecular processes by which they are transcribed and regulated are relatively well understood in Drosophila melanogaster, much less is known about the origin and evolution of piCs in this or any other species. To investigate piC evolution, we use a population genomics approach to compare piC activity and sequence composition across 8 geographically distant strains of D. melanogaster with high quality long-read genome assemblies. We perform extensive annotations of ovary piCs and TE content in each strain and test predictions of two proposed models of piC evolution. The 'de novo' model posits that individual TE insertions can spontaneously attain the status of a small piC to generate piRNAs silencing the entire TE family. The 'trap' model envisions large and evolutionary stable genomic clusters where TEs tend to accumulate and serves as a long-term "memory" of ancient TE invasions and produce a great variety of piRNAs protecting against related TEs entering the genome. It remains unclear which model best describes the evolution of piCs. Our analysis uncovers extensive variation in piC activity across strains and signatures of rapid birth and death of piCs in natural populations. Most TE families inferred to be recently or currently active show an enrichment of strain-specific insertions into large piCs, consistent with the trap model. By contrast, only a small subset of active LTR retrotransposon families is enriched for the formation of strain-specific piCs, suggesting that these families have an inherent proclivity to form de novo piCs. Thus, our findings support aspects of both 'de novo' and 'trap' models of piC evolution. We propose that these two models represent two extreme stages along an evolutionary continuum, which begins with the emergence of piCs de novo from a few specific LTR retrotransposon insertions that subsequently expand by accretion of other TE insertions during evolution to form larger 'trap' clusters. Our study shows that piCs are evolutionarily labile and that TEs themselves are the major force driving the formation and evolution of piCs.

Resolution of structural variation in diverse mouse genomes reveals chromatin remodeling due to transposable elements

Diverse inbred mouse strains are important biomedical research models, yet genome characterization of many strains is fundamentally lacking in comparison with humans. In particular, catalogs of structural variants (SVs) (variants ≥ 50 bp) are incomplete, limiting the discovery of causative alleles for phenotypic variation. Here, we resolve genome-wide SVs in 20 genetically distinct inbred mice with long-read sequencing. We report 413,758 site-specific SVs affecting 13% (356 Mbp) of the mouse reference assembly, including 510 previously unannotated coding variants. We substantially improve the Mus musculus transposable element (TE) callset, and we find that TEs comprise 39% of SVs and account for 75% of altered bases. We further utilize this callset to investigate how TE heterogeneity affects mouse embryonic stem cells and find multiple TE classes that influence chromatin accessibility. Our work provides a comprehensive analysis of SVs found in diverse mouse genomes and illustrates the role of TEs in epigenetic differences.

Widespread somatic L1 retrotransposition in normal colorectal epithelium

Throughout an individual's lifetime, genomic alterations accumulate in somatic cells1-11. However, the mutational landscape induced by retrotransposition of long interspersed nuclear element-1 (L1), a widespread mobile element in the human genome12-14, is poorly understood in normal cells. Here we explored the whole-genome sequences of 899 single-cell clones established from three different cell types collected from 28 individuals. We identified 1,708 somatic L1 retrotransposition events that were enriched in colorectal epithelium and showed a positive relationship with age. Fingerprinting of source elements showed 34 retrotransposition-competent L1s. Multidimensional analysis demonstrated that (1) somatic L1 retrotranspositions occur from early embryogenesis at a substantial rate, (2) epigenetic on/off of a source element is preferentially determined in the early organogenesis stage, (3) retrotransposition-competent L1s with a lower population allele frequency have higher retrotransposition activity and (4) only a small fraction of L1 transcripts in the cytoplasm are finally retrotransposed in somatic cells. Analysis of matched cancers further suggested that somatic L1 retrotransposition rate is substantially increased during colorectal tumourigenesis. In summary, this study illustrates L1 retrotransposition-induced somatic mosaicism in normal cells and provides insights into the genomic and epigenomic regulation of transposable elements over the human lifetime.

Repetitive elements in aging and neurodegeneration

Repetitive elements (REs), such as transposable elements (TEs) and satellites, comprise much of the genome. Here, we review how TEs and (peri)centromeric satellite DNA may contribute to aging and neurodegenerative disorders, including amyotrophic lateral sclerosis (ALS). Alterations in RE expression, retrotransposition, and chromatin microenvironment may shorten lifespan, elicit neurodegeneration, and impair memory and movement. REs may cause these phenotypes via DNA damage, protein sequestration, insertional mutagenesis, and inflammation. We discuss several TE families, including gypsy, HERV-K, and HERV-W, and how TEs interact with various factors, including transactive response (TAR) DNA-binding protein 43 kDa (TDP-43) and the siRNA and piwi-interacting (pi)RNA systems. Studies of TEs in neurodegeneration have focused on Drosophila and, thus, further examination in mammals is needed. We suggest that therapeutic silencing of REs could help mitigate neurodegenerative disorders.

High-resolution Nanopore methylome-maps reveal random hyper-methylation at CpG-poor regions as driver of chemoresistance in leukemias

Aberrant DNA methylation at CpG dinucleotides is a cancer hallmark that is associated with the emergence of resistance to anti cancer treatment, though molecular mechanisms and biological significance remain elusive. Genome scale methylation maps by currently used methods are based on chemical modification of DNA and are best suited for analyses of methylation at CpG rich regions (CpG islands). We report the first high coverage whole-genome map in cancer using the long read nanopore technology, which allows simultaneous DNA-sequence and -methylation analyses on native DNA. We analyzed clonal epigenomic/genomic evolution in Acute Myeloid Leukemias (AMLs) at diagnosis and relapse, after chemotherapy. Long read sequencing coupled to a novel computational method allowed definition of differential methylation at unprecedented resolution, and showed that the relapse methylome is characterized by hypermethylation at both CpG islands and sparse CpGs regions. Most differentially methylated genes, however, were not differentially expressed nor enriched for chemoresistance genes. A small fraction of under-expressed and hyper-methylated genes at sparse CpGs, in the gene body, was significantly enriched in transcription factors (TFs). Remarkably, these few TFs supported large gene-regulatory networks including 50% of all differentially expressed genes in the relapsed AMLs and highly-enriched in chemoresistance genes. Notably, hypermethylated regions at sparse CpGs were poorly conserved in the relapsed AMLs, under-represented at their genomic positions and showed higher methylation entropy, as compared to CpG islands. Analyses of available datasets confirmed TF binding to their target genes and conservation of the same gene-regulatory networks in large patient cohorts. Relapsed AMLs carried few patient specific structural variants and DNA mutations, apparently not involved in drug resistance. Thus, drug resistance in AMLs can be mainly ascribed to the selection of random epigenetic alterations at sparse CpGs of a few transcription factors, which then induce reprogramming of the relapsing phenotype, independently of clonal genomic evolution.

TrEMOLO: accurate transposable element allele frequency estimation using long-read sequencing data combining assembly and mapping-based approaches

Transposable Element MOnitoring with LOng-reads (TrEMOLO) is a new software that combines assembly- and mapping-based approaches to robustly detect genetic elements called transposable elements (TEs). Using high- or low-quality genome assemblies, TrEMOLO can detect most TE insertions and deletions and estimate their allele frequency in populations. Benchmarking with simulated data revealed that TrEMOLO outperforms other state-of-the-art computational tools. TE detection and frequency estimation by TrEMOLO were validated using simulated and experimental datasets. Therefore, TrEMOLO is a comprehensive and suitable tool to accurately study TE dynamics. TrEMOLO is available under GNU GPL3.0 at https://github.com/DrosophilaGenomeEvolution/TrEMOLO .

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.

The interferon stimulated gene-encoded protein HELZ2 inhibits human LINE-1 retrotransposition and LINE-1 RNA-mediated type I interferon induction

Some interferon stimulated genes (ISGs) encode proteins that inhibit LINE-1 (L1) retrotransposition. Here, we use immunoprecipitation followed by liquid chromatography-tandem mass spectrometry to identify proteins that associate with the L1 ORF1-encoded protein (ORF1p) in ribonucleoprotein particles. Three ISG proteins that interact with ORF1p inhibit retrotransposition: HECT and RLD domain containing E3 ubiquitin-protein ligase 5 (HERC5); 2'-5'-oligoadenylate synthetase-like (OASL); and helicase with zinc finger 2 (HELZ2). HERC5 destabilizes ORF1p, but does not affect its cellular localization. OASL impairs ORF1p cytoplasmic foci formation. HELZ2 recognizes sequences and/or structures within the L1 5'UTR to reduce L1 RNA, ORF1p, and ORF1p cytoplasmic foci levels. Overexpression of WT or reverse transcriptase-deficient L1s lead to a modest induction of IFN-α expression, which is abrogated upon HELZ2 overexpression. Notably, IFN-α expression is enhanced upon overexpression of an ORF1p RNA binding mutant, suggesting ORF1p binding might protect L1 RNA from "triggering" IFN-α induction. Thus, ISG proteins can inhibit retrotransposition by different mechanisms.

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.

Experimental Approaches to Study Somatic Transposition in Drosophila Using Whole-Genome DNA Sequencing

The extent of transposable element (TE) mobilization in different somatic tissues and throughout diverse species is not well understood. Somatic transposition is often challenging to study as it generates de novo TE insertions that represent rare genetic variants present in heterogenous tissues. Here, we describe experimental approaches that can be applied to address TE mobility in somatic tissues with the use of short- and long-read whole-genome DNA sequencing. Focusing on the analysis of the Drosophila melanogaster intestinal and head tissues, we provide instructions on how to design, perform, and validate experiments that aim at detecting somatic transposition. In addition to providing examples of protocols, this chapter intends to deliver general experimental guidelines that may be adapted to other fly tissues or to other species.

Targeted Nanopore Resequencing and Methylation Analysis of LINE-1 Retrotransposons

Retrotransposition of LINE-1 (L1) elements represents a major source of insertional polymorphisms in mammals, and their mutagenic activity is restricted by silencing mechanisms, such as DNA methylation. Despite a very high level of sequence identity between copies, their internal sequence contains small nucleotide polymorphisms (SNPs) that can alter their activity. Such internal SNPs can also appear in different alleles of a given L1 locus. Given their repetitive nature and relatively long size, short-read sequencing approaches have limited access to L1 internal sequence or DNA methylation state. Here, we describe a targeted method to specifically sequence more than a hundred L1-containing loci in parallel and measure their DNA methylation levels using nanopore long-read sequencing. Each targeted locus is sequenced at high coverage (~45X) with unambiguously mapped reads spanning the entire L1 element, as well as its flanking sequences over several kilobases. Our protocol, modified from the nanopore Cas9 targeted sequencing (nCATS) strategy, provides a full and haplotype-resolved L1 sequence and DNA methylation levels. It introduces a streamlined and multiplex approach to synthesize guide RNAs and a quantitative PCR (qPCR)-based quality check during library preparation for cost-effective L1 sequencing. More generally, this method can be applied to any type of transposable elements and organisms.

Nanopore Sequencing to Identify Transposable Element Insertions and Their Epigenetic Modifications

Over the past 20 years, high-throughput genomic assays have fundamentally changed how transposable elements (TEs) are studied. While short-read DNA sequencing has been at the heart of these efforts, novel technologies that generate longer reads are driving a shift in the field. Long-read sequencing now permits locus-specific approaches to locate individual TE insertions and understand their epigenetic and transcriptional regulation, while still profiling TE activity genome-wide. Here we provide detailed guidelines to implement Oxford Nanopore Technologies (ONT) sequencing to identify polymorphic TE insertions and profile TE epigenetic landscapes. Using human long interspersed element-1 (LINE-1, L1) as an example, we explain the procedures involved, including final visualization, and potential bottlenecks and pitfalls. ONT sequencing will be, in our view, a workhorse technology for the foreseeable future in the TE field.

Single-Cell DNA Methylation Analysis in Cancer

Morphological, transcriptomic, and genomic defects are well-explored parameters of cancer biology. In more recent years, the impact of epigenetic influences, such as DNA methylation, is becoming more appreciated. Aberrant DNA methylation has been implicated in many types of cancers, influencing cell type, state, transcriptional regulation, and genomic stability to name a few. Traditionally, large populations of cells from the tissue of interest are coalesced for analysis, producing averaged methylome data. Considering the inherent heterogeneity of cancer, analysing populations of cells as a whole denies the ability to discover novel aberrant methylation patterns, identify subpopulations, and trace cell lineages. Due to recent advancements in technology, it is now possible to obtain methylome data from single cells. This has both research and clinical implications, ranging from the identification of biomarkers to improved diagnostic tools. As with all emerging technologies, distinct experimental, bioinformatic, and practical challenges present themselves. This review begins with exploring the potential impact of single-cell sequencing on understanding cancer biology and how it could eventually benefit a clinical setting. Following this, the techniques and experimental approaches which made this technology possible are explored. Finally, the present challenges currently associated with single-cell DNA methylation sequencing are described.

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.

Retrotransposon activation during Drosophila metamorphosis conditions adult antiviral responses

Retrotransposons are one type of mobile genetic element that abundantly reside in the genomes of nearly all animals. Their uncontrolled activation is linked to sterility, cancer and other pathologies, thereby being largely considered detrimental. Here we report that, within a specific time window of development, retrotransposon activation can license the host's immune system for future antiviral responses. We found that the mdg4 (also known as Gypsy) retrotransposon selectively becomes active during metamorphosis at the Drosophila pupal stage. At this stage, mdg4 activation educates the host's innate immune system by inducing the systemic antiviral function of the nuclear factor-κB protein Relish in a dSTING-dependent manner. Consequently, adult flies with mdg4, Relish or dSTING silenced at the pupal stage are unable to clear exogenous viruses and succumb to viral infection. Altogether, our data reveal that hosts can establish a protective antiviral response that endows a long-term benefit in pathogen warfare due to the developmental activation of mobile genetic elements.

A review of strategies used to identify transposition events in plant genomes

Transposable elements (TEs) were initially considered redundant and dubbed 'junk DNA'. However, more recently they were recognized as an essential element of genome plasticity. In nature, they frequently become active upon exposition of the host to stress conditions. Even though most transposition events are neutral or even deleterious, occasionally they may happen to be beneficial, resulting in genetic novelty providing better fitness to the host. Hence, TE mobilization may promote adaptability and, in the long run, act as a significant evolutionary force. There are many examples of TE insertions resulting in increased tolerance to stresses or in novel features of crops which are appealing to the consumer. Possibly, TE-driven de novo variability could be utilized for crop improvement. However, in order to systematically study the mechanisms of TE/host interactions, it is necessary to have suitable tools to globally monitor any ongoing TE mobilization. With the development of novel potent technologies, new high-throughput strategies for studying TE dynamics are emerging. Here, we present currently available methods applied to monitor the activity of TEs in plants. We divide them on the basis of their operational principles, the position of target molecules in the process of transposition and their ability to capture real cases of actively transposing elements. Their possible theoretical and practical drawbacks are also discussed. Finally, conceivable strategies and combinations of methods resulting in an improved performance are proposed.

The critical role of epigenetic mechanisms involved in nanotoxicology

Over the past decades, nanomaterials (NMs) have been widely applied in the cosmetic, food, engineering, and medical fields. Along with the prevalence of NMs, the toxicological characteristics exhibited by these materials on health and the environment have gradually attracted attentions. A growing number of evidences have indicated that epigenetics holds an essential role in the onset and development of various diseases. NMs could cause epigenetic alterations such as DNA methylation, noncoding RNA (ncRNA) expression, and histone modifications. NMs might alternate either global DNA methylation or the methylation of specific genes to affect the biological function. Abnormal upregulation or downregulation of ncRNAs might also be a potential mechanism for the toxic effects caused by NMs. In parallel, the phosphorylation, acetylation, and methylation of histones also take an important part in the process of NMs-induced toxicity. As the adverse effects of NMs continue to be explored, mechanisms such as chromosomal remodeling, genomic imprinting, and m⁶ A modification are also gradually coming into the limelight. Since the epigenetic alterations often occur in the early development of diseases, thus the relevant studies not only provide insight into the pathogenesis of diseases, but also screen for the prospective biomarkers for early diagnosis and prevention. This review summarizes the epigenetic alterations elicited by NMs, hoping to provide a clue for nanotoxicity studies and security evaluation of NMs. This article is categorized under: Toxicology and Regulatory Issues in Nanomedicine > Toxicology of Nanomaterials.

Transcription start site signal profiling improves transposable element RNA expression analysis at locus-level

The transcriptional activity of Transposable Elements (TEs) has been involved in numerous pathological processes, including neurodegenerative diseases such as amyotrophic lateral sclerosis and frontotemporal lobar degeneration. The TE expression analysis from short-read sequencing technologies is, however, challenging due to the multitude of similar sequences derived from singular TEs subfamilies and the exaptation of TEs within longer coding or non-coding RNAs. Specialised tools have been developed to quantify the expression of TEs that either relies on probabilistic re-distribution of multimapper count fractions or allow for discarding multimappers altogether. Until now, the benchmarking across those tools was largely limited to aggregated expression estimates over whole TEs subfamilies. Here, we compared the performance of recently published tools (SQuIRE, TElocal, SalmonTE) with simplistic quantification strategies (featureCounts in unique, fraction and random modes) at the individual loci level. Using simulated datasets, we examined the false discovery rate and the primary driver of those false positive hits in the optimal quantification strategy. Our findings suggest a high false discovery number that exceeds the total number of correctly recovered active loci for all the quantification strategies, including the best performing tool TElocal. As a remedy, filtering based on the minimum number of read counts or baseMean expression improves the F1 score and decreases the number of false positives. Finally, we demonstrate that additional profiling of Transcription Start Site mapping statistics (using a k-means clustering approach) significantly improves the performance of TElocal while reporting a reliable set of detected and differentially expressed TEs in human simulated RNA-seq data.

Local assembly of long reads enables phylogenomics of transposable elements in a polyploid cell line

Animal cell lines often undergo extreme genome restructuring events, including polyploidy and segmental aneuploidy that can impede de novo whole-genome assembly (WGA). In some species like Drosophila, cell lines also exhibit massive proliferation of transposable elements (TEs). To better understand the role of transposition during animal cell culture, we sequenced the genome of the tetraploid Drosophila S2R+ cell line using long-read and linked-read technologies. WGAs for S2R+ were highly fragmented and generated variable estimates of TE content across sequencing and assembly technologies. We therefore developed a novel WGA-independent bioinformatics method called TELR that identifies, locally assembles, and estimates allele frequency of TEs from long-read sequence data (https://github.com/bergmanlab/telr). Application of TELR to a ∼130x PacBio dataset for S2R+ revealed many haplotype-specific TE insertions that arose by transposition after initial cell line establishment and subsequent tetraploidization. Local assemblies from TELR also allowed phylogenetic analysis of paralogous TEs, which revealed that proliferation of TE families in vitro can be driven by single or multiple source lineages. Our work provides a model for the analysis of TEs in complex heterozygous or polyploid genomes that are recalcitrant to WGA and yields new insights into the mechanisms of genome evolution in animal cell culture.

Epigenetic tumor heterogeneity in the era of single-cell profiling with nanopore sequencing

Nanopore sequencing has brought the technology to the next generation in the science of sequencing. This is achieved through research advancing on: pore efficiency, creating mechanisms to control DNA translocation, enhancing signal-to-noise ratio, and expanding to long-read ranges. Heterogeneity regarding epigenetics would be broad as mutations in the epigenome are sensitive to cause new challenges in cancer research. Epigenetic enzymes which catalyze DNA methylation and histone modification are dysregulated in cancer cells and cause numerous heterogeneous clones to evolve. Detection of this heterogeneity in these clones plays an indispensable role in the treatment of various cancer types. With single-cell profiling, the nanopore sequencing technology could provide a simple sequence at long reads and is expected to be used soon at the bedside or doctor's office. Here, we review the advancements of nanopore sequencing and its use in the detection of epigenetic heterogeneity in cancer.

Ongoing endeavors to detect mobilization of transposable elements

Transposable elements (TEs) are DNA sequences capable of mobilization from one location to another in the genome. Since the discovery of ‘Dissociation (Dc) locus’ by Barbara McClintock in maize (1), mounting evidence in the era of genomics indicates that a significant fraction of most eukaryotic genomes is composed of TE sequences, involving in various aspects of biological processes such as development, physiology, diseases and evolution. Although technical advances in genomics have discovered numerous functional impacts of TE across species, our understanding of TEs is still ongoing process due to challenges resulted from complexity and abundance of TEs in the genome. In this mini-review, we briefly summarize biology of TEs and their impacts on the host genome, emphasizing importance of understanding TE landscape in the genome. Then, we introduce recent endeavors especially in vivo retrotransposition assays and long read sequencing technology for identifying de novo insertions/TE polymorphism, which will broaden our knowledge of extraordinary relationship between genomic cohabitants and their host.

The Methylation Inhibitor 5-Aza-2'-Deoxycytidine Induces Genome-Wide Hypomethylation in Rice

DNA methylation is a conserved epigenetic modification which is vital for regulating gene expression and maintaining genome stability in both mammals and plants. Homozygous mutation of rice methyltransferase 1 (met1) gene can cause host death in rice, making it difficult to obtain plant material needed for hypomethylation research. To circumvent this challenge, the methylation inhibitor, 5-Aza-2'-deoxycytidine (AzaD), is used as a cytosine nucleoside analogue to reduce genome wide hypomethylation and is widely used in hypomethylation research. However, how AzaD affects plant methylation profiles at the genome scale is largely unknown. Here, we treated rice seedlings with AzaD and compared the AzaD treatment with osmet1-2 mutants, illustrating that there are similar CG hypomethylation and distribution throughout the whole genome. Along with global methylation loss class I transposable elements (TEs) which are farther from genes compared with class II TEs, were more significantly activated, and the RNA-directed DNA Methylation (RdDM) pathway was activated in specific genomic regions to compensate for severe CG loss. Overall, our results suggest that AzaD is an effective DNA methylation inhibitor that can influence genome wide methylation and cause a series of epigenetic variations.

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.

Ongoing endeavors to detect mobilization of transposable elements

Transposable elements (TEs) are DNA sequences capable of mobilization from one location to another in the genome. Since the discovery of 'Dissociation (Dc) locus' by Barbara McClintock in maize (1), mounting evidence in the era of genomics indicates that a significant fraction of most eukaryotic genomes is composed of TE sequences, involving in various aspects of biological processes such as development, physiology, diseases and evolution. Although technical advances in genomics have discovered numerous functional impacts of TE across species, our understanding of TEs is still ongoing process due to challenges resulted from complexity and abundance of TEs in the genome. In this mini-review, we briefly summarize biology of TEs and their impacts on the host genome, emphasizing importance of understanding TE landscape in the genome. Then, we introduce recent endeavors especially in vivo retrotransposition assays and long read sequencing technology for identifying de novo insertions/TE polymorphism, which will broaden our knowledge of extraordinary relationship between genomic cohabitants and their host. [BMB Reports 2022; 55(7): 305-315].

The Tumor-Specific Expression of L1 Retrotransposons Independently Correlates with Time to Relapse in Hormone-Negative Breast Cancer Patients

BACKGROUND

Long-Interspersed Nuclear Element (L1) retrotransposons are silenced in healthy tissues but unrepressed in cancer. Even if L1 reactivation has been associated with reduced overall survival in breast cancer (BC) patients, a comprehensive correlation with clinicopathological features is still missing.

METHODS

Using quantitative, reverse-transcription PCR, we assessed L1 mRNA expression in 12 BC cells, 210 BC patients and in 47 normal mammary tissues. L1 expression was then correlated with molecular and clinicopathological data.

RESULTS

We identified a tumor-exclusive expression of L1s, absent in normal mammary cells and tissues. A positive correlation between L1 expression and tumor dedifferentiation, lymph-node involvement and increased immune infiltration was detected. Molecular subtyping highlighted an enrichment of L1s in basal-like cells and cancers. By exploring disease-free survival, we identified L1 overexpression as an independent biomarker for patients with a high risk of recurrence in hormone-receptor-negative BCs.

CONCLUSIONS

Overall, L1 reactivation identified BCs with aggressive features and patients with a worse clinical fate.

Genomic features underlie the co-option of SVA transposons as cis-regulatory elements in human pluripotent stem cells

Domestication of transposable elements (TEs) into functional cis-regulatory elements is a widespread phenomenon. However, the mechanisms behind why some TEs are co-opted as functional enhancers while others are not are underappreciated. SINE-VNTR-Alus (SVAs) are the youngest group of transposons in the human genome, where ~3,700 copies are annotated, nearly half of which are human-specific. Many studies indicate that SVAs are among the most frequently co-opted TEs in human gene regulation, but the mechanisms underlying such processes have not yet been thoroughly investigated. Here, we leveraged CRISPR-interference (CRISPRi), computational and functional genomics to elucidate the genomic features that underlie SVA domestication into human stem-cell gene regulation. We found that ~750 SVAs are co-opted as functional cis-regulatory elements in human induced pluripotent stem cells. These SVAs are significantly closer to genes and harbor more transcription factor binding sites than non-co-opted SVAs. We show that a long DNA motif composed of flanking YY1/2 and OCT4 binding sites is enriched in the co-opted SVAs and that these two transcription factors bind consecutively on the TE sequence. We used CRISPRi to epigenetically repress active SVAs in stem cell-like NCCIT cells. Epigenetic perturbation of active SVAs strongly attenuated YY1/OCT4 binding and influenced neighboring gene expression. Ultimately, SVA repression resulted in ~3,000 differentially expressed genes, 131 of which were the nearest gene to an annotated SVA. In summary, we demonstrated that SVAs modulate human gene expression, and uncovered that location and sequence composition contribute to SVA domestication into gene regulatory networks.

Methylartist: Tools for Visualising Modified Bases from Nanopore Sequence Data

Summary

Methylartist is a consolidated suite of tools for processing, visualizing and analysing nanopore-derived modified base calls. All detectable methylation types (e.g. 5mCpG, 5hmC, 6mA) are supported, enabling integrated study of base pairs when modified naturally or as part of an experimental protocol.

Availability and implementation

Methylartist is implemented in Python and is installable via PyPI and bioconda. Source code and test data are available at https://github.com/adamewing/methylartist.

Supplementary information

Supplementary data are available at Bioinformatics online.

Single-Cell Epigenomics Reveals Mechanisms of Cancer Progression

Cancer initiation is driven by the cooperation between genetic and epigenetic aberrations that disrupt gene regulatory programs critical to maintaining specialized cellular functions. After initiation, cells acquire additional genetic and epigenetic alterations influenced by tumor-intrinsic and -extrinsic mechanisms, which increase intratumoral heterogeneity, reshape the cell's underlying gene regulatory networks and promote cancer evolution. Furthermore, environmental or therapeutic insults drive the selection of heterogeneous cell states, with implications for cancer initiation, maintenance, and drug resistance. The advancement of single-cell genomics has begun to uncover the full repertoire of chromatin and gene expression states (cell states) that exist within individual tumors. These single-cell analyses suggest that cells diversify in their regulatory states upon transformation by co-opting damage-induced and nonlineage regulatory programs that can lead to epigenomic plasticity. Here, we review these recent studies related to regulatory state changes in cancer progression and highlight the growing single-cell epigenomics toolkit poised to address unresolved questions in the field.

Locus-specific chromatin profiling of evolutionarily young transposable elements

Despite a vast expansion in the availability of epigenomic data, our knowledge of the chromatin landscape at interspersed repeats remains highly limited by difficulties in mapping short-read sequencing data to these regions. In particular, little is known about the locus-specific regulation of evolutionarily young transposable elements (TEs), which have been implicated in genome stability, gene regulation and innate immunity in a variety of developmental and disease contexts. Here we propose an approach for generating locus-specific protein-DNA binding profiles at interspersed repeats, which leverages information on the spatial proximity between repetitive and non-repetitive genomic regions. We demonstrate that the combination of HiChIP and a newly developed mapping tool (PAtChER) yields accurate protein enrichment profiles at individual repetitive loci. Using this approach, we reveal previously unappreciated variation in the epigenetic profiles of young TE loci in mouse and human cells. Insights gained using our method will be invaluable for dissecting the molecular determinants of TE regulation and their impact on the genome.

Locus-specific expression of transposable elements in single cells with CELLO-seq

Transposable elements (TEs) regulate diverse biological processes, from early development to cancer. Expression of young TEs is difficult to measure with next-generation, single-cell sequencing technologies because their highly repetitive nature means that short complementary DNA reads cannot be unambiguously mapped to a specific locus. Single CELl LOng-read RNA-sequencing (CELLO-seq) combines long-read single cell RNA-sequencing with computational analyses to measure TE expression at unique loci. We used CELLO-seq to assess the widespread expression of TEs in two-cell mouse blastomeres as well as in human induced pluripotent stem cells. Across both species, old and young TEs showed evidence of locus-specific expression with simulations demonstrating that only a small number of very young elements in the mouse could not be mapped back to the reference with high confidence. Exploring the relationship between the expression of individual elements and putative regulators revealed large heterogeneity, with TEs within a class showing different patterns of correlation and suggesting distinct regulatory mechanisms.

Locus specific reduction of L1 expression in the cortices of individuals with amyotrophic lateral sclerosis

The activation and dysregulation of retrotransposons has been identified in the CNS of individuals with the fatal neurodegenerative disorder Amyotrophic lateral sclerosis (ALS). This includes elements from multiple different families and subfamilies of retrotransposons, however there is limited knowledge of the specific loci from which this expression occurs in ALS. The long interspersed element-1 (L1) is the only autonomous retrotransposon in the human genome and members of this family of elements maintain the ability to mobilise. Despite L1s contributing to 17% of the human genome only 80-100 L1s encode the required proteins for mobilisation and are retrotransposition competent. Identifying the specific loci from which L1 expression occurs will inform on the potential functional consequences of their expression, such as the potential for somatic retrotransposition or DNA damage caused by the endonuclease activity of the ORF2 protein of the L1. Here we characterised L1 loci expression using the L1EM tool ( https://github.com/FenyoLab/L1EM ) in RNA sequencing data from 518 samples across four tissues (motor cortex, frontal cortex, cerebellum and cervical spinal cord) in the Target ALS cohort obtained from the New York Genome Center. There was a significant reduction in total intact L1 expression (those that encode functional proteins) in two brain regions of individuals with ALS compared to controls and clustering of the ALS brain regions occurred based on their intact L1 expression profile. Although overall the levels of L1 expression were reduced in ALS/ALS with other neurological disorder (ND) there were individuals in which L1s were expressed at much higher levels than the rest of the ALS/ALSND cohort. Expressed L1 loci were more frequently located in introns compared to those not expressed and the level of L1 expression positively correlated with the expression of the gene in which it was located. Significant differences were observed in the expression profiles of L1s in ALS and specific features of these elements, such as location in the genome and whether or not they are intact, were significantly associated with those that were expressed in the cohort.

A hidden layer of structural variation in transposable elements reveals potential genetic modifiers in human disease-risk loci

Genome-wide association studies (GWAS) have been highly informative in discovering disease-associated loci but are not designed to capture all structural variations in the human genome. Using long-read sequencing data, we discovered widespread structural variation within SINE-VNTR-Alu (SVA) elements, a class of great ape-specific transposable elements with gene-regulatory roles, which represents a major source of structural variability in the human population. We highlight the presence of structurally variable SVAs (SV-SVAs) in neurological disease-associated loci, and we further associate SV-SVAs to disease-associated SNPs and differential gene expression using luciferase assays and expression quantitative trait loci data. Finally, we genetically deleted SV-SVAs in the BIN1 and CD2AP Alzheimer's disease-associated risk loci and in the BCKDK Parkinson's disease-associated risk locus and assessed multiple aspects of their gene-regulatory influence in a human neuronal context. Together, this study reveals a novel layer of genetic variation in transposable elements that may contribute to identification of the structural variants that are the actual drivers of disease associations of GWAS loci.

Transcriptional states of retroelement-inserted regions and specific KRAB zinc finger protein association are correlated with DNA methylation of retroelements in human male germ cells

DNA methylation, repressive histone modifications, and PIWI-interacting RNAs are essential for controlling retroelement silencing in mammalian germ lines. Dysregulation of retroelement silencing is associated with male sterility. Although retroelement silencing mechanisms have been extensively studied in mouse germ cells, little progress has been made in humans. Here, we show that the Krüppel-associated box domain zinc finger proteins are associated with DNA methylation of retroelements in human primordial germ cells. Further, we show that the hominoid-specific retroelement SINE-VNTR-Alus (SVA) is subjected to transcription-directed de novo DNA methylation during human spermatogenesis. The degree of de novo DNA methylation in SVAs varies among human individuals, which confers significant inter-individual epigenetic variation in sperm. Collectively, our results highlight potential molecular mechanisms for the regulation of retroelements in human male germ cells.

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.

Comprehensive genetic diagnosis of tandem repeat expansion disorders with programmable targeted nanopore sequencing

More than 50 neurological and neuromuscular diseases are caused by short tandem repeat (STR) expansions, with 37 different genes implicated to date. We describe the use of programmable targeted long-read sequencing with Oxford Nanopore's ReadUntil function for parallel genotyping of all known neuropathogenic STRs in a single assay. Our approach enables accurate, haplotype-resolved assembly and DNA methylation profiling of STR sites, from a list of predetermined candidates. This correctly diagnoses all individuals in a small cohort (n = 37) including patients with various neurogenetic diseases (n = 25). Targeted long-read sequencing solves large and complex STR expansions that confound established molecular tests and short-read sequencing and identifies noncanonical STR motif conformations and internal sequence interruptions. We observe a diversity of STR alleles of known and unknown pathogenicity, suggesting that long-read sequencing will redefine the genetic landscape of repeat disorders. Last, we show how the inclusion of pharmacogenomic genes as secondary ReadUntil targets can further inform patient care.

Making sense of IL-6 signalling cues in pathophysiology

Unravelling the molecular mechanisms that account for functional pleiotropy is a major challenge for researchers in cytokine biology. Cytokine-receptor cross-reactivity and shared signalling pathways are considered primary drivers of cytokine pleiotropy. However, reports epitomized by studies of Jak-STAT cytokine signalling identify interesting biochemical and epigenetic determinants of transcription factor regulation that affect the delivery of signal-dependent cytokine responses. Here, a regulatory interplay between STAT transcription factors and their convergence to specific genomic enhancers support the fine-tuning of cytokine responses controlling host immunity, functional identity, and tissue homeostasis and repair. In this review, we provide an overview of the signalling networks that shape the way cells sense and interpret cytokine cues. With an emphasis on the biology of interleukin-6, we highlight the importance of these mechanisms to both physiological processes and pathophysiological outcomes.

Long-read sequencing of the zebrafish genome reorganizes genomic architecture

BACKGROUND

Nanopore sequencing technology has revolutionized the field of genome biology with its ability to generate extra-long reads that can resolve regions of the genome that were previously inaccessible to short-read sequencing platforms. Over 50% of the zebrafish genome consists of difficult to map, highly repetitive, low complexity elements that pose inherent problems for short-read sequencers and assemblers.

RESULTS

We used long-read nanopore sequencing to generate a de novo assembly of the zebrafish genome and compared our assembly to the current reference genome, GRCz11. The new assembly identified 1697 novel insertions and deletions over one kilobase in length and placed 106 previously unlocalized scaffolds. We also discovered additional sites of retrotransposon integration previously unreported in GRCz11 and observed the expression of these transposable elements in adult zebrafish under physiologic conditions, implying they have active mobility in the zebrafish genome and contribute to the ever-changing genomic landscape.

CONCLUSIONS

We used nanopore sequencing to improve upon and resolve the issues plaguing the current zebrafish reference assembly, GRCz11. Zebrafish is a prominent model of human disease, and our corrected assembly will be useful for studies relying on interspecies comparisons and precise linkage of genetic events to disease phenotypes.

Elucidating Hexanucleotide Repeat Number and Methylation within the X-Linked Dystonia-Parkinsonism (XDP)-Related SVA Retrotransposon in TAF1 with Nanopore Sequencing

Background: X-linked dystonia-parkinsonism (XDP) is an adult-onset neurodegenerative disorder characterized by progressive dystonia and parkinsonism. It is caused by a SINE-VNTR-Alu (SVA) retrotransposon insertion in the TAF1 gene with a polymorphic (CCCTCT)n domain that acts as a genetic modifier of disease onset and expressivity. Methods: Herein, we used Nanopore sequencing to investigate SVA genetic variability and methylation. We used blood-derived DNA from 96 XDP patients for amplicon-based deep Nanopore sequencing and validated it with fragment analysis which was performed using fluorescence-based PCR. To detect methylation from blood- and brain-derived DNA, we used a Cas9-targeted approach. Results: High concordance was observed for hexanucleotide repeat numbers detected with Nanopore sequencing and fragment analysis. Within the SVA locus, there was no difference in genetic variability other than variations of the repeat motif between patients. We detected high CpG methylation frequency (MF) of the SVA and flanking regions (mean MF = 0.94, SD = ±0.12). Our preliminary results suggest only subtle differences between the XDP patient and the control in predicted enhancer sites directly flanking the SVA locus. Conclusions: Nanopore sequencing can reliably detect SVA hexanucleotide repeat numbers, methylation and, lastly, variation in the repeat motif.

Fast nanopore sequencing data analysis with SLOW5

Nanopore sequencing depends on the FAST5 file format, which does not allow efficient parallel analysis. Here we introduce SLOW5, an alternative format engineered for efficient parallelization and acceleration of nanopore data analysis. Using the example of DNA methylation profiling of a human genome, analysis runtime is reduced from more than two weeks to approximately 10.5 h on a typical high-performance computer. SLOW5 is approximately 25% smaller than FAST5 and delivers consistent improvements on different computer architectures.

Nanopore sequencing data are rapidly analyzed with parallel data access.

Emerging roles for endogenous retroviruses in immune epigenetic regulation

In recent years, there has been significant progress toward understanding the transcriptional networks underlying mammalian immune responses, fueled by advances in regulatory genomic technologies. Epigenomic studies profiling immune cells have generated detailed genome-wide maps of regulatory elements that will be key to deciphering the regulatory networks underlying cellular immune responses and autoimmune disorders. Unbiased analyses of these genomic maps have uncovered endogenous retroviruses as an unexpected ally in the regulation of human immune systems. Despite their parasitic origins, studies are finding an increasing number of examples of retroviral sequences having been co-opted for beneficial immune function and regulation by the host cell. Here, we review how endogenous retroviruses have given rise to numerous regulatory elements that shape the epigenetic landscape of host immune responses. We will discuss the implications of these elements on the function, dysfunction, and evolution of innate immunity.

Single-Run Catalysis and Kinetic Control of Human Telomerase Holoenzyme

Genome stability in eukaryotic cells relies on proper maintenance of telomeres at the termini of linear chromosomes. Human telomerase holoenzyme is required for maintaining telomere stability in a majority of proliferative human cells, making it essential for control of cell division and aging, stem cell maintenance, and development and survival of tumor or cancer. A dividing human cell usually contains a limited number of active telomerase holoenzymes. Recently, we discovered that a human telomerase catalytic site undergoes catalysis-dependent shut-off and an inactive site can be reactivated by cellular fractions containing human intracellular telomerase-activating factors (hiTAFs). Such ON-OFF control of human telomerase activity suggests a dynamic switch between inactive and active pools of the holoenzymes. In this review, we will link the ON-OFF control to the thermodynamic and kinetic properties of human telomerase holoenzymes, and discuss its potential contributions to the maintenance of telomere length equilibrium. This treatment suggests probabilistic fluctuations in the number of active telomerase holoenzymes as well as the number of telomeres that are extended in a limited number of cell cycles, and may be an important component of a fully quantitative model for the dynamic control of telomerase activities and telomere lengths in different types of eukaryotic cells.

Third-Generation Sequencing: The Spearhead towards the Radical Transformation of Modern Genomics

Although next-generation sequencing (NGS) technology revolutionized sequencing, offering a tremendous sequencing capacity with groundbreaking depth and accuracy, it continues to demonstrate serious limitations. In the early 2010s, the introduction of a novel set of sequencing methodologies, presented by two platforms, Pacific Biosciences (PacBio) and Oxford Nanopore Sequencing (ONT), gave birth to third-generation sequencing (TGS). The innovative long-read technologies turn genome sequencing into an ease-of-handle procedure by greatly reducing the average time of library construction workflows and simplifying the process of de novo genome assembly due to the generation of long reads. Long sequencing reads produced by both TGS methodologies have already facilitated the decipherment of transcriptional profiling since they enable the identification of full-length transcripts without the need for assembly or the use of sophisticated bioinformatics tools. Long-read technologies have also provided new insights into the field of epitranscriptomics, by allowing the direct detection of RNA modifications on native RNA molecules. This review highlights the advantageous features of the newly introduced TGS technologies, discusses their limitations and provides an in-depth comparison regarding their scientific background and available protocols as well as their potential utility in research and clinical applications.

Transposons Hidden in Arabidopsis thaliana Genome Assembly Gaps and Mobilization of Non-Autonomous LTR Retrotransposons Unravelled by Nanotei Pipeline

Long-read data is a great tool to discover new active transposable elements (TEs). However, no ready-to-use tools were available to gather this information from low coverage ONT datasets. Here, we developed a novel pipeline, nanotei, that allows detection of TE-contained structural variants, including individual TE transpositions. We exploited this pipeline to identify TE insertion in the Arabidopsis thaliana genome. Using nanotei, we identified tens of TE copies, including ones for the well-characterized ONSEN retrotransposon family that were hidden in genome assembly gaps. The results demonstrate that some TEs are inaccessible for analysis with the current A. thaliana (TAIR10.1) genome assembly. We further explored the mobilome of the ddm1 mutant with elevated TE activity. Nanotei captured all TEs previously known to be active in ddm1 and also identified transposition of non-autonomous TEs. Of them, one non-autonomous TE derived from (AT5TE33540) belongs to TR-GAG retrotransposons with a single open reading frame (ORF) encoding the GAG protein. These results provide the first direct evidence that TR-GAGs and other non-autonomous LTR retrotransposons can transpose in the plant genome, albeit in the absence of most of the encoded proteins. In summary, nanotei is a useful tool to detect active TEs and their insertions in plant genomes using low-coverage data from Nanopore genome sequencing.

Application of third-generation sequencing in cancer research

In the past several years, nanopore sequencing technology from Oxford Nanopore Technologies (ONT) and single-molecule real-time (SMRT) sequencing technology from Pacific BioSciences (PacBio) have become available to researchers and are currently being tested for cancer research. These methods offer many advantages over most widely used high-throughput short-read sequencing approaches and allow the comprehensive analysis of transcriptomes by identifying full-length splice isoforms and several other posttranscriptional events. In addition, these platforms enable structural variation characterization at a previously unparalleled resolution and direct detection of epigenetic marks in native DNA and RNA. Here, we present a comprehensive summary of important applications of these technologies in cancer research, including the identification of complex structure variants, alternatively spliced isoforms, fusion transcript events, and exogenous RNA. Furthermore, we discuss the impact of the newly developed nanopore direct RNA sequencing (RNA-Seq) approach in advancing epitranscriptome research in cancer. Although the unique challenges still present for these new single-molecule long-read methods, they will unravel many aspects of cancer genome complexity in unprecedented ways and present an encouraging outlook for continued application in an increasing number of different cancer research settings.

Mutagenesis of human genomes by endogenous mobile elements on a population scale

Several large-scale Illumina whole-genome sequencing (WGS) and whole-exome sequencing (WES) projects have emerged recently that have provided exceptional opportunities to discover mobile element insertions (MEIs) and study the impact of these MEIs on human genomes. However, these projects also have presented major challenges with respect to the scalability and computational costs associated with performing MEI discovery on tens or even hundreds of thousands of samples. To meet these challenges, we have developed a more efficient and scalable version of our mobile element locator tool (MELT) called CloudMELT. We then used MELT and CloudMELT to perform MEI discovery in 57,919 human genomes and exomes, leading to the discovery of 104,350 nonredundant MEIs. We leveraged this collection (1) to examine potentially active L1 source elements that drive the mobilization of new Alu, L1, and SVA MEIs in humans; (2) to examine the population distributions and subfamilies of these MEIs; and (3) to examine the mutagenesis of GENCODE genes, ENCODE-annotated features, and disease genes by these MEIs. Our study provides new insights on the L1 source elements that drive MEI mutagenesis and brings forth a better understanding of how this mutagenesis impacts human genomes.