Nanopore Sequencing Enables Comprehensive Transposable Element Epigenomic Profiling

Advances in understanding LINE-1 regulation and function in the human genome

LINE-1 (long interspersed nuclear element 1, L1) retrotransposons constitute ~17% of human DNA (~0.5 million genomic L1 copies) and exhibit context-dependent expression in different cell lines. Recent studies reveal that L1 is under multilayered control by diverse factors that either collaborate or compete with each other to ensure precise L1 activity. Remarkably, L1s have been co-opted as various transcription-dependent regulatory elements, such as promoters, enhancers, and topologically associating domain (TAD) boundaries, that regulate gene expression in zygotic genome activation, aging, cancer, and other disorders. This review highlights the regulation of L1 and its regulatory functions that influence disease and development.

Toward clinical long-read genome sequencing for rare diseases

Genetic diagnostics is driven by technological advances, forming a tight interface between research, clinic and industry, which enables rapid implementation of new technologies. Short-read genome and exome sequencing, the current state of the art in clinical genetics, can detect a broad spectrum of genetic variants across the genome. However, despite these advancements, more than half of individuals with rare diseases remain undiagnosed after genomic investigations. Long-read whole-genome sequencing (LR-WGS) is a promising technology that identifies previously difficult-to-detect variants while also enabling phasing and methylation analysis and has the potential of generating complete personal assemblies. To pave the way for clinical use of LR-WGS, the clinical genomic community must establish standardized protocols and quality parameters while also developing innovative tools for data analysis and interpretation. In this Perspective, we explore the key challenges and benefits in integrating LR-WGS into routine clinical diagnostics.

Transposable elements as genome regulators in normal and malignant haematopoiesis

Transposable elements (TEs) constitute over half of the human genome and have played a profound role in genome evolution. While most TEs have lost the ability to transpose, many retain functional elements that serve as drivers of genome innovation, including the emergence of novel genes and regulatory elements. Recent advances in experimental and bioinformatic methods have provided new insights into their roles in human biology, both in health and disease. In this review, we discuss the multifaceted roles of TEs in haematopoiesis, highlighting their contributions to both normal and pathological contexts. TEs influence gene regulation by reshaping gene-regulatory networks, modulating transcriptional activity, and creating novel regulatory elements. These activities play key roles in maintaining normal haematopoietic processes and supporting cellular regeneration. However, in haematological malignancies, TE reactivation can disrupt genomic integrity, induce structural variations, and dysregulate transcriptional programmes, thereby driving oncogenesis. By examining the impact of TE activity on genome regulation and variation, we highlight their pivotal roles in both normal haematopoietic processes and haematological cancers.

Structural features of somatic and germline retrotransposition events in humans

BACKGROUND

Transposons are DNA sequences able to move or copy themselves to other genomic locations leading to insertional mutagenesis. Although transposon-derived sequences account for half of the human genome, most elements are no longer transposition competent. Moreover, transposons are normally repressed through epigenetic silencing in healthy adult tissues but become derepressed in several human cancers, with high activity detected in colorectal cancer. Their impact on non-malignant and malignant tissue as well as the differences between somatic and germline retrotransposition remain poorly understood. With new sequencing technologies, including long read sequencing, we can access intricacies of retrotransposition, such as insertion sequence details and nested repeats, that have been previously challenging to characterize.

RESULTS

In this study, we investigate somatic and germline retrotransposition by analyzing long read sequencing from 56 colorectal cancers and 112 uterine leiomyomas. We identified 1495 somatic insertions in colorectal samples, while striking lack of insertions was detected in uterine leiomyomas. Our findings highlight differences between somatic and germline events, such as transposon type distribution, insertion length, and target site preference. Leveraging long-read sequencing, we provide an in-depth analysis of the twin-priming phenomenon, detecting it across transposable element types that remain active in humans, including Alus. Additionally, we detect an abundance of germline transposons in repetitive DNA, along with a relationship between replication timing and insertion target site.

CONCLUSIONS

Our study reveals a stark contrast in somatic transposon activity between colorectal cancers and uterine leiomyomas, and highlights differences between somatic and germline transposition. This suggests potentially different conditions in malignant and non-malignant tissues, as well as in germline and somatic tissues, which could be involved in the transposition process. Long-read sequencing provided important insights into transposon behavior, allowing detailed examination of structural features such as twin priming and nested elements.

Unraveling the hidden complexity of cancer through long-read sequencing

Cancer is fundamentally a disease of the genome, characterized by extensive genomic, transcriptomic, and epigenomic alterations. Most current studies predominantly use short-read sequencing, gene panels, or microarrays to explore these alterations; however, these technologies can systematically miss or misrepresent certain types of alterations, especially structural variants, complex rearrangements, and alterations within repetitive regions. Long-read sequencing is rapidly emerging as a transformative technology for cancer research by providing a comprehensive view across the genome, transcriptome, and epigenome, including the ability to detect alterations that previous technologies have overlooked. In this Perspective, we explore the current applications of long-read sequencing for both germline and somatic cancer analysis. We provide an overview of the computational methodologies tailored to long-read data and highlight key discoveries and resources within cancer genomics that were previously inaccessible with prior technologies. We also address future opportunities and persistent challenges, including the experimental and computational requirements needed to scale to larger sample sizes, the hurdles in sequencing and analyzing complex cancer genomes, and opportunities for leveraging machine learning and artificial intelligence technologies for cancer informatics. We further discuss how the telomere-to-telomere genome and the emerging human pangenome could enhance the resolution of cancer genome analysis, potentially revolutionizing early detection and disease monitoring in patients. Finally, we outline strategies for transitioning long-read sequencing from research applications to routine clinical practice.

OGT prevents DNA demethylation and suppresses the expression of transposable elements in heterochromatin by restraining TET activity genome-wide

O-GlcNAc transferase (OGT) interacts robustly with all three mammalian TET methylcytosine dioxygenases. Here we show that deletion of the Ogt gene in mouse embryonic stem (mES) cells results in a widespread increase in the TET product 5-hydroxymethylcytosine in both euchromatic and heterochromatic compartments, with a concomitant reduction in the TET substrate 5-methylcytosine at the same genomic regions. mES cells treated with an OGT inhibitor also displayed increased 5-hydroxymethylcytosine, and attenuating the TET1-OGT interaction in mES cells resulted in a genome-wide decrease of 5-methylcytosine, indicating that OGT restrains TET activity and limits inappropriate DNA demethylation in a manner that requires the TET-OGT interaction and the catalytic activity of OGT. DNA hypomethylation in OGT-deficient cells was accompanied by derepression of transposable elements predominantly located in heterochromatin. We suggest that OGT protects the genome against TET-mediated DNA demethylation and loss of heterochromatin integrity, preventing the aberrant increase in transposable element expression noted in cancer, autoimmune-inflammatory diseases, cellular senescence and aging.

Computational analysis of DNA methylation from long-read sequencing

DNA methylation is a critical epigenetic mechanism in numerous biological processes, including gene regulation, development, ageing and the onset of various diseases such as cancer. Studies of methylation are increasingly using single-molecule long-read sequencing technologies to simultaneously measure epigenetic states such as DNA methylation with genomic variation. These long-read data sets have spurred the continuous development of advanced computational methods to gain insights into the roles of methylation in regulating chromatin structure and gene regulation. In this Review, we discuss the computational methods for calling methylation signals, contrasting methylation between samples, analysing cell-type diversity and gaining additional genomic insights, and then further discuss the challenges and future perspectives of tool development for DNA methylation research.

Third generation sequencing transforming plant genome research: Current trends and challenges

In recent years, third-generation sequencing (TGS) technologies have transformed genomics and transcriptomics research, providing novel opportunities for significant discoveries. The long-read sequencing platforms, with their unique advantages over next-generation sequencing (NGS), including a definitive protocol, reduced operational time, and real-time sequencing, possess the potential to transform plant genomics. TGS optimizes and enhances the efficiency of data analysis by removing the necessity for time-consuming assembly tools. The current review examines the development and application of bioinformatics tools for data analysis and annotation, driven by the rapid advancement of TGS platforms like Oxford Nanopore Technologies and Pacific Biosciences. Transcriptome analysis utilizing TGS has been extensively employed to elucidate complex plant transcriptomes and genomes, particularly those characterized by high frequencies of duplicated genomes and repetitive sequences. As a result, current methodologies that allow for generating transcriptomes and comprehensive whole-genome sequences of complex plant genomes employing tailored hybrid sequencing techniques that integrate NGS and TGS technologies have been emphasized herein. This paper, thus, articulates a vision for a future in which TGS effectively addresses the challenges faced in plant research, offering a comprehensive understanding of its advantages, applications, limitations, and promising prospects.

Nanopore Long-Read Sequencing Unveils Genomic Disruptions in Alzheimer's Disease

Studies in laboratory models and postmortem human brain tissue from patients with Alzheimer's disease have revealed disruption of basic cellular processes such as DNA repair and epigenetic control as drivers of neurodegeneration. While genomic alterations in regions of the genome that are rich in repetitive sequences, often termed "dark regions," are difficult to resolve using traditional sequencing approaches, long-read technologies offer promising new avenues to explore previously inaccessible regions of the genome. In the current study, we leverage nanopore-based long-read whole-genome sequencing of DNA extracted from postmortem human frontal cortex at early and late stages of Alzheimer's disease, as well as age-matched controls, to analyze retrotransposon insertion events, non-allelic homologous recombination (NAHR), structural variants and DNA methylation within retrotransposon loci and other repetitive/dark regions of the human genome. Interestingly, we find that retrotransposon insertion events and repetitive element-associated NAHR are particularly enriched within centromeric and pericentromeric regions of DNA in the aged human brain, and that ribosomal DNA (rDNA) is subject to a high degree of NAHR compared to other regions of the genome. We detect a trending increase in potential somatic retrotransposition events of the small interfering nuclear element (SINE) AluY in late-stage Alzheimer's disease, and differential changes in methylation within repetitive elements and retrotransposons according to disease stage. Taken together, our analysis provides the first long-read DNA sequencing-based analysis of retrotransposon sequences, NAHR, structural variants, and DNA methylation in the aged brain, and points toward transposable elements, centromeric/pericentromeric regions and rDNA as hotspots for genomic variation.

TE-Seq: A Transposable Element Annotation and RNA-Seq Pipeline

The recognition that transposable elements (TEs) play important roles in many biological processes has elicited growing interest in analyzing sequencing data derived from this mobile 'dark genome'. The TE-Seq pipeline conducts an end-to-end analysis of RNA-sequencing data, examining both genes and TEs. It implements the most current computational methods tailor-made for TEs, enabling a comprehensive analysis of TE expression at both the individual element level and at the TE clade level. If supplied with long-read DNA sequencing data, it creates a TE-complete genome incorporating non-reference (polymorphic) TE-loci, enabling the functional characterization of the evolutionarily youngest mobile elements in the genome.

Comparative Genomics Reveals LINE-1 Recombination with Diverse RNAs

Long interspersed element-1 (LINE-1, L1) retrotransposons are the most abundant protein-coding transposable elements (TE) in mammalian genomes, and have shaped genome content over 170 million years of evolution. LINE-1 is self-propagating and mobilizes other sequences, including Alu elements. Occasionally, LINE-1 forms chimeric insertions with non-coding RNAs and mRNAs. U6 spliceosomal small nuclear RNA/LINE-1 chimeras are best known, though there are no comprehensive catalogs of LINE-1 chimeras. To address this, we developed TiMEstamp, a computational pipeline that leverages multiple sequence alignments (MSA) to estimate the age of LINE-1 insertions and identify candidate chimeric insertions where an adjacent sequence arrives contemporaneously. Candidates were refined by detecting hallmark features of L1 retrotransposition, such as target site duplication (TSD). Applying this pipeline to the human genome, we recovered all known species of LINE-1 chimeras and discovered new chimeric insertions involving small RNAs, Alu elements, and mRNA fragments. Some insertions are compatible with known mechanisms, such as RNA ligation. Other structures nominate novel mechanisms, such as trans-splicing. We also see evidence that LINE-1 loci with defunct promoters can acquire regulatory elements from nearby genes to restore retrotransposition activity. These discoveries highlight the recombinatory potential of LINE-1 RNA with implications for genome evolution and TE domestication.

Stratification system with dual human endogenous retroviruses for predicting immunotherapy efficacy in metastatic clear-cell renal cell carcinoma

BACKGROUND

Endogenous retrovirus (ERV) elements are genomic footprints of ancestral retroviral infections within the human genome. While the dysregulation of ERV transcription has been linked to immune cell infiltration in various cancers, its relationship with immune checkpoint inhibitor (ICI) response in solid tumors, particularly metastatic clear-cell renal cell carcinoma (ccRCC), remains inadequately explored.

METHODS

This study analyzed patients with metastatic ccRCC from two prospective clinical trials, encompassing 181 patients receiving nivolumab in the CheckMate trials (-009 to -010 and -025) and 48 patients treated with the ipilimumab-nivolumab combination in the BIONIKK trial. ERV expression was quantified using the ERVmap algorithm from RNA sequencing data. Our primary objective was to correlate ERV expression with progression-free survival, with overall survival and time-to-second-treatment survival as secondary endpoints. We used bootstrap methods with univariate Cox regression on 666 substantially expressed ERVs to evaluate their prognostic significance and stability.

RESULTS

Our analysis centered on two ERVs, E4421_chr17 and E1659_chr4, which consistently exhibited opposing prognostic impacts across both cohorts. We developed a stratification system based on their median expression levels, categorizing patients into four ERV subgroups. These subgroups were further consolidated into a three-tier risk model that significantly correlated with ICI treatment outcomes. The most responsive ERV risk category showed enhanced endothelial cell infiltration, whereas the resistant category was characterized by higher levels of myeloid dendritic cells, regulatory T cells, myeloid-derived suppressor cells, and markers of T-cell exhaustion. Notably, this ERV-based classification outperformed traditional transcriptomic signatures in predicting ICI efficacy and showed further improvement when combined with epigenetic DNA methylation markers.

CONCLUSIONS

Our findings introduce a dual ERV-based stratification system that effectively categorizes patient risk and predicts clinical outcomes for ccRCC patients undergoing ICI therapy. Beyond enhancing the predictive precision of existing transcriptomic models, this system paves the way for more targeted and individualized approaches in the realm of precision oncology.

SUMMER: an integrated nanopore sequencing pipeline for variants detection and clinical annotation on the human genome

Long-read sequencing has emerged as a transformative technology in recent years, offering significant potential for the molecular diagnosis of unresolved genetic disorders. Despite its promise, the comprehensive detection and clinical annotation of genomic variants remain intricate and technically demanding. We present SUMMER, an integrated and structured workflow specifically designed to process raw Nanopore sequencing reads. SUMMER facilitates an in-depth analysis of multiple variant types, including SNV, SV, short tandem repeat and mobile element insertion. For clinical applications, SUMMER employs SvAnna to prioritize SV candidates based on phenotype relevance and utilizes Straglr to provide reference distributions of non-pathogenic unit counts for 55 known pathogenic short tandem repeats. By addressing critical challenges in variant detection and annotation, SUMMER seeks to advance the clinical utility of long-read sequencing in diagnostic genomics. SUMMER is available on the web at https://github.com/carolhuaxia/summer .

Shedding light on DNA methylation and its clinical implications: the impact of long-read-based nanopore technology

DNA methylation is an essential epigenetic mechanism for regulation of gene expression, through which many physiological (X-chromosome inactivation, genetic imprinting, chromatin structure and miRNA regulation, genome defense, silencing of transposable elements) and pathological processes (cancer and repetitive sequences-associated diseases) are regulated. Nanopore sequencing has emerged as a novel technique that can analyze long strands of DNA (long-read sequencing) without chemically treating the DNA. Interestingly, nanopore sequencing can also extract epigenetic status of the nucleotides (including both 5-Methylcytosine and 5-hydroxyMethylcytosine), and a large variety of bioinformatic tools have been developed for improving its detection properties. Out of all genomic regions, long read sequencing provides advantages in studying repetitive elements, which are difficult to characterize through other sequencing methods. Transposable elements are repetitive regions of the genome that are silenced and usually display high levels of DNA methylation. Their demethylation and activation have been observed in many cancers. Due to their repetitive nature, it is challenging to accurately estimate DNA methylation levels within transposable elements using short sequencing technologies. The advantage to sequence native DNA (without PCR amplification biases or harsh bisulfite treatment) and long and ultra long reads coupled with epigenetic states of the DNA allows to accurately estimate DNA methylation levels in transposable elements. This is a big step forward for epigenomic studies, and unsolved questions regarding gene expression and transposable elements silencing through DNA methylation can now be answered.

Critical assessment of nanopore sequencing for the detection of multiple forms of DNA modifications

While nanopore sequencing is increasingly used for mapping DNA modifications, it is important to recognize false positive calls as they can mislead biological interpretations. To assist biologists and methods developers, we describe a framework for rigorous evaluation that highlights the use of false discovery rate with rationally designed negative controls capturing both general background and confounding modifications. Our critical assessment across multiple forms of DNA modifications highlights that while nanopore sequencing performs reliably for high-abundance modifications, including 5-methylcytosine (5mC) at CpG sites in mammalian cells and 5-hydroxymethylcytosine (5hmC) in mammalian brain cells, it makes a significant proportion of false positive detections for low-abundance modifications, such as 5mC at CpH sites, 5hmC and N6-methyldeoxyadenine (6mA) in most mammal cell types. This study highlights the urgent need to incorporate this framework in future methods development and biological studies, and advocates prioritizing nanopore sequencing for mapping abundant over rare modifications in biomedical applications.

Structural variations in livestock genomes and their associations with phenotypic traits: a review

Genomic structural variation (SV) refers to differences in gene sequences between individuals on a genomic scale. It is widely distributed in the genome, primarily in the form of insertions, deletions, duplications, inversions, and translocations. Due to its characterization by long segments and large coverage, SVs significantly impact the genetic characteristics and production performance of livestock, playing a crucial role in studying breed diversity, biological evolution, and disease correlation. Research on SVs contributes to an enhanced understanding of chromosome function and genetic characteristics and is important for understanding hereditary diseases mechanisms. In this article, we review the concept, classification, main formation mechanisms, detection methods, and advancement of research on SVs in the genomes of cattle, buffalo, equine, sheep, and goats, aiming to reveal the genetic basis of differences in phenotypic traits and adaptive genetic mechanisms through genomic research, which will provide a theoretical basis for better understanding and utilizing the genetic resources of herbivorous livestock.

Single cell long read whole genome sequencing reveals somatic transposon activity in human brain

The advent of single cell DNA sequencing revealed astonishing dynamics of genomic variability, but failed at characterizing smaller to mid size variants that on the germline level have a profound impact. In this work we discover novel dynamics in three brains utilizing single cell long-read sequencing. This provides key insights into the dynamic of the genomes of individual cells and further highlights brain specific activity of transposable elements.

Detecting transposable elements in long-read genomes using sTELLeR

MOTIVATION

Repeat elements, such as transposable elements (TE), are highly repetitive DNA sequences that compose around 50% of the genome. TEs such as Alu, SVA, HERV, and L1 elements can cause disease through disrupting genes, causing frameshift mutations or altering splicing patters. These are elements challenging to characterize using short-read genome sequencing, due to its read length and TEs repetitive nature. Long-read genome sequencing (lrGS) enables bridging of TEs, allowing increased resolution across repetitive DNA sequences. lrGS therefore present an opportunity for improved TE detection and analysis not only from a research perspective but also for future clinical detection. When choosing an lrGS TE caller, parameters such as runtime, CPU hours, sensitivity, precision, and compatibility with inclusion into pipelines are crucial for efficient detection.

RESULTS

We therefore developed sTELLeR, (s) Transposable ELement in Long (e) Read, for accurate, fast, and effective TE detection. Particularly, sTELLeR exhibit higher precision and sensitivity for calling of Alu elements than similar tools. The caller is 5-48× as fast and uses <2% of the CPU hours compared to competitive callers. The caller is haplotype aware and output results in a variant call format (VCF) file, enabling compatibility with other variant callers and downstream analysis.

AVAILABILITY AND IMPLEMENTATION

sTELLeR is a python-based tool and is available at https://github.com/kristinebilgrav/sTELLeR. Altogether, we show that sTELLeR is a fast, sensitive, and precise caller for detection of TE elements, and can easily be implemented into variant calling workflows.

A unified framework to analyze transposable element insertion polymorphisms using graph genomes

Transposable elements are ubiquitous mobile DNA sequences generating insertion polymorphisms, contributing to genomic diversity. We present GraffiTE, a flexible pipeline to analyze polymorphic mobile elements insertions. By integrating state-of-the-art structural variant detection algorithms and graph genomes, GraffiTE identifies polymorphic mobile elements from genomic assemblies or long-read sequencing data, and genotypes these variants using short or long read sets. Benchmarking on simulated and real datasets reports high precision and recall rates. GraffiTE is designed to allow non-expert users to perform comprehensive analyses, including in models with limited transposable element knowledge and is compatible with various sequencing technologies. Here, we demonstrate the versatility of GraffiTE by analyzing human, Drosophila melanogaster, maize, and Cannabis sativa pangenome data. These analyses reveal the landscapes of polymorphic mobile elements and their frequency variations across individuals, strains, and cultivars.

Mini-heterochromatin domains constrain the cis-regulatory impact of SVA transposons in human brain development and disease

SVA (SINE (short interspersed nuclear element)-VNTR (variable number of tandem repeats)-Alu) retrotransposons remain active in humans and contribute to individual genetic variation. Polymorphic SVA alleles harbor gene regulatory potential and can cause genetic disease. However, how SVA insertions are controlled and functionally impact human disease is unknown. Here we dissect the epigenetic regulation and influence of SVAs in cellular models of X-linked dystonia parkinsonism (XDP), a neurodegenerative disorder caused by an SVA insertion at the TAF1 locus. We demonstrate that the KRAB zinc finger protein ZNF91 establishes H3K9me3 and DNA methylation over SVAs, including polymorphic alleles, in human neural progenitor cells. The resulting mini-heterochromatin domains attenuate the cis-regulatory impact of SVAs. This is critical for XDP pathology; removal of local heterochromatin severely aggravates the XDP molecular phenotype, resulting in increased TAF1 intron retention and reduced expression. Our results provide unique mechanistic insights into how human polymorphic transposon insertions are recognized and how their regulatory impact is constrained by an innate epigenetic defense system.

Targeting transposable elements in cancer: developments and opportunities

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

DNA methylation governs the sensitivity of repeats to restriction by the HUSH-MORC2 corepressor

The human silencing hub (HUSH) complex binds to transcripts of LINE-1 retrotransposons (L1s) and other genomic repeats, recruiting MORC2 and other effectors to remodel chromatin. How HUSH and MORC2 operate alongside DNA methylation, a central epigenetic regulator of repeat transcription, remains largely unknown. Here we interrogate this relationship in human neural progenitor cells (hNPCs), a somatic model of brain development that tolerates removal of DNA methyltransferase DNMT1. Upon loss of MORC2 or HUSH subunit TASOR in hNPCs, L1s remain silenced by robust promoter methylation. However, genome demethylation and activation of evolutionarily-young L1s attracts MORC2 binding, and simultaneous depletion of DNMT1 and MORC2 causes massive accumulation of L1 transcripts. We identify the same mechanistic hierarchy at pericentromeric α-satellites and clustered protocadherin genes, repetitive elements important for chromosome structure and neurodevelopment respectively. Our data delineate the epigenetic control of repeats in somatic cells, with implications for understanding the vital functions of HUSH-MORC2 in hypomethylated contexts throughout human development.

DeepBAM: a high-accuracy single-molecule CpG methylation detection tool for Oxford nanopore sequencing

Recent nanopore sequencing system (R10.4) has enhanced base calling accuracy and is being increasingly utilized for detecting CpG methylation state. However, the robustness and universality of the methylation calling model in officially supplied Dorado remains poorly tested. In this study, we obtained heterogeneous datasets from human and plant sources to carry out comprehensive evaluations, which showed that Dorado performed significantly different across datasets. We therefore developed deep neural networks and implemented several optimizations in training a new model called DeepBAM. DeepBAM achieved superior and more stable performances compared with Dorado, including higher area under the ROC curves (98.47% on average and up to 7.36% improvement) and F1 scores (94.97% on average and up to 16.24% improvement) across the datasets. DeepBAM-based whole genome methylation frequencies have achieved >0.95 correlations with BS-seq on four of five datasets, outperforming Dorado in all instances. It enables unraveling allele-specific methylation patterns, including regions of transposable elements. The enhanced performance of DeepBAM paves the way for broader applications of nanopore sequencing in CpG methylation studies.

Local read haplotagging enables accurate long-read small variant calling

Long-read sequencing technology has enabled variant detection in difficult-to-map regions of the genome and enabled rapid genetic diagnosis in clinical settings. Rapidly evolving third-generation sequencing platforms like Pacific Biosciences (PacBio) and Oxford Nanopore Technologies (ONT) are introducing newer platforms and data types. It has been demonstrated that variant calling methods based on deep neural networks can use local haplotyping information with long-reads to improve the genotyping accuracy. However, using local haplotype information creates an overhead as variant calling needs to be performed multiple times which ultimately makes it difficult to extend to new data types and platforms as they get introduced. In this work, we have developed a local haplotype approximate method that enables state-of-the-art variant calling performance with multiple sequencing platforms including PacBio Revio system, ONT R10.4 simplex and duplex data. This addition of local haplotype approximation simplifies long-read variant calling with DeepVariant.

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.

Viral genome sequencing methods: benefits and pitfalls of current approaches

Whole genome sequencing of viruses provides high-resolution molecular insights, enhancing our understanding of viral genome function and phylogeny. Beyond fundamental research, viral sequencing is increasingly vital for pathogen surveillance, epidemiology, and clinical applications. As sequencing methods rapidly evolve, the diversity of viral genomics applications and catalogued genomes continues to expand. Advances in long-read, single molecule, real-time sequencing methodologies present opportunities to sequence contiguous, haplotype resolved viral genomes in a range of research and applied settings. Here we present an overview of nucleic acid sequencing methods and their applications in studying viral genomes. We emphasise the advantages of different viral sequencing approaches, with a particular focus on the benefits of third-generation sequencing technologies in elucidating viral evolution, transmission networks, and pathogenesis.

Rapid evolution of piRNA clusters in the Drosophila melanogaster ovary

The piRNA pathway is a highly conserved mechanism to repress transposable element (TE) activity in the animal germline via a specialized class of small RNAs called piwi-interacting RNAs (piRNAs). piRNAs are produced from discrete genomic regions called piRNA clusters (piCs). Although the molecular processes by which piCs function 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 origin and evolution, we use a population genomic approach to compare piC activity and sequence composition across eight geographically distant strains of D. melanogaster with high-quality long-read genome assemblies. We perform annotations of ovary piCs and genome-wide TE content in each strain. Our analysis uncovers extensive variation in piC activity across strains and signatures of rapid birth and death of piCs. Most TEs inferred to be recently active show an enrichment of insertions into old and large piCs, consistent with the previously proposed "trap" model of piC evolution. In contrast, a small subset of active LTR families is enriched for the formation of new piCs, suggesting that these TEs have higher proclivity to form piCs. Thus, our findings uncover processes leading to the origin of piCs. We propose that piC evolution begins with the emergence of piRNAs from individual insertions of a few select TE families prone to seed new piCs that subsequently expand by accretion of insertions from most other TE families during evolution to form larger "trap" clusters. Our study shows that TEs themselves are the major force driving the rapid evolution of piCs.

Structure and transcription of integrated HPV DNA in vulvar carcinomas

HPV infections are associated with a fraction of vulvar cancers. Through hybridization capture and DNA sequencing, HPV DNA was detected in five of thirteen vulvar cancers. HPV16 DNA was integrated into human DNA in three of the five. The insertions were in introns of human NCKAP1, C5orf67, and LRP1B. Integrations in NCKAP1 and C5orf67 were flanked by short direct repeats in the human DNA, consistent with HPV DNA insertions at sites of abortive, staggered, endonucleolytic incisions. The insertion in C5orf67 was present as a 36 kbp, human-HPV-hetero-catemeric DNA as either an extrachromosomal circle or a tandem repeat within the human genome. The human circularization/repeat junction was defined at single nucleotide resolution. The integrated viral DNA segments all retained an intact upstream regulatory region and the adjacent viral E6 and E7 oncogenes. RNA sequencing revealed that the only HPV genes consistently transcribed from the integrated viral DNAs were E7 and E6*I. The other two HPV DNA+ tumors had coinfections, but no evidence for integration. HPV-positive and HPV-negative vulvar cancers exhibited contrasting human, global gene expression patterns partially overlapping with previously observed differences between HPV-positive and HPV-negative cervical and oropharyngeal cancers. A substantial fraction of the differentially expressed genes involved immune system function. Thus, transcription and HPV DNA integration in vulvar cancers resemble those in other HPV-positive cancers. This study emphasizes the power of hybridization capture coupled with DNA and RNA sequencing to identify a broad spectrum of HPV types, determine human genome integration status of viral DNAs, and elucidate their structures.

Neural-net-based cell deconvolution from DNA methylation reveals tumor microenvironment associated with cancer prognosis

DNA methylation is a pivotal epigenetic modification that defines cellular identity. While cell deconvolution utilizing this information is considered useful for clinical practice, current methods for deconvolution are limited in their accuracy and resolution. In this study, we collected DNA methylation data from 945 human samples derived from various tissues and tumor-infiltrating immune cells and trained a neural network model with them. The model, termed MEnet, predicted abundance of cell population together with the detailed immune cell status from bulk DNA methylation data, and showed consistency to those of flow cytometry and histochemistry. MEnet was superior to the existing methods in the accuracy, speed, and detectable cell diversity, and could be applicable for peripheral blood, tumors, cell-free DNA, and formalin-fixed paraffin-embedded sections. Furthermore, by applying MEnet to 72 intrahepatic cholangiocarcinoma samples, we identified immune cell profiles associated with cancer prognosis. We believe that cell deconvolution by MEnet has the potential for use in clinical settings.

Comprehensive profiling of L1 retrotransposons in mouse

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

Regulation and function of transposable elements in cancer genomes

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

The Application of Long-Read Sequencing to Cancer

Cancer is a multifaceted disease arising from numerous genomic aberrations that have been identified as a result of advancements in sequencing technologies. While next-generation sequencing (NGS), which uses short reads, has transformed cancer research and diagnostics, it is limited by read length. Third-generation sequencing (TGS), led by the Pacific Biosciences and Oxford Nanopore Technologies platforms, employs long-read sequences, which have marked a paradigm shift in cancer research. Cancer genomes often harbour complex events, and TGS, with its ability to span large genomic regions, has facilitated their characterisation, providing a better understanding of how complex rearrangements affect cancer initiation and progression. TGS has also characterised the entire transcriptome of various cancers, revealing cancer-associated isoforms that could serve as biomarkers or therapeutic targets. Furthermore, TGS has advanced cancer research by improving genome assemblies, detecting complex variants, and providing a more complete picture of transcriptomes and epigenomes. This review focuses on TGS and its growing role in cancer research. We investigate its advantages and limitations, providing a rigorous scientific analysis of its use in detecting previously hidden aberrations missed by NGS. This promising technology holds immense potential for both research and clinical applications, with far-reaching implications for cancer diagnosis and treatment.

DNA Origami-Enabled Gene Localization of Repetitive Sequences

Repetitive sequences, which make up over 50% of human DNA, have diverse applications in disease diagnosis, forensic identification, paternity testing, and population genetic analysis due to their crucial functions for gene regulation. However, representative detection technologies such as sequencing and fluorescence imaging suffer from time-consuming protocols, high cost, and inaccuracy of the position and order of repetitive sequences. Here, we develop a precise and cost-effective strategy that combines the high resolution of atomic force microscopy with the shape customizability of DNA origami for repetitive sequence-specific gene localization. "Tri-block" DNA structures were specifically designed to connect repetitive sequences to DNA origami tags, thereby revealing precise genetic information in terms of position and sequence for high-resolution and high-precision visualization of repetitive sequences. More importantly, we achieved the results of simultaneous detection of different DNA repetitive sequences on the gene template with a resolution of ∼6.5 nm (19 nt). This strategy is characterized by high efficiency, high precision, low operational complexity, and low labor/time costs, providing a powerful complement to sequencing technologies for gene localization of repetitive sequences.

Nanopore DNA sequencing technologies and their applications towards single-molecule proteomics

Sequencing of nucleic acids with nanopores has emerged as a powerful tool offering rapid readout, high accuracy, low cost and portability. This label-free method for sequencing at the single-molecule level is an achievement on its own. However, nanopores also show promise for the technologically even more challenging sequencing of polypeptides, something that could considerably benefit biological discovery, clinical diagnostics and homeland security, as current techniques lack portability and speed. Here we survey the biochemical innovations underpinning commercial and academic nanopore DNA/RNA sequencing techniques, and explore how these advances can fuel developments in future protein sequencing with nanopores.

Pharmacological inhibition of neddylation impairs long interspersed element 1 retrotransposition

Aberrant long interspersed element 1 (LINE-1 or L1) activity can cause insertional mutagenesis and chromosomal rearrangements and has been detected in several types of cancers. Here, we show that neddylation, a post-translational modification process, is essential for L1 transposition. The antineoplastic drug MLN4924 is an L1 inhibitor that suppresses NEDD8-activating enzyme activity. Neddylation inhibition by MLN4924 selectively impairs ORF2p-mediated L1 reverse transcription and blocks the generation of L1 cDNA. Consistent with these results, MLN4924 treatment suppresses the retrotransposition activity of the non-autonomous retrotransposons short interspersed nuclear element R/variable number of tandem repeat/Alu and Alu, which rely on the reverse transcription activity of L1 ORF2p. The E2 enzyme UBE2M in the neddylation pathway, rather than UBE2F, is required for L1 ORF2p and retrotransposition. Interference with the functions of certain neddylation-dependent Cullin-really interesting new gene E3 ligases disrupts L1 reverse transcription and transposition activity. Our findings provide insights into the regulation of L1 retrotransposition and the identification of therapeutic targets for L1 dysfunctions.

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.

Towards targeting transposable elements for cancer therapy

Transposable elements (TEs) represent almost half of the human genome. Historically deemed 'junk DNA', recent technological advancements have stimulated a wave of research into the functional impact of TEs on gene-regulatory networks in evolution and development, as well as in diseases including cancer. The genetic and epigenetic evolution of cancer involves the exploitation of TEs, whereby TEs contribute directly to cancer-specific gene activities. This Review provides a perspective on the role of TEs in cancer as being a 'double-edged sword', both promoting cancer evolution and representing a vulnerability that could be exploited in cancer therapy. We discuss how TEs affect transcriptome regulation and other cellular processes in cancer. We highlight the potential of TEs as therapeutic targets for cancer. We also summarize technical hurdles in the characterization of TEs with genomic assays. Last, we outline open questions and exciting future research avenues.

OGT prevents DNA demethylation and suppresses the expression of transposable elements in heterochromatin by restraining TET activity genome-wide

The O- GlcNAc transferase OGT interacts robustly with all three mammalian TET methylcytosine dioxygenases. We show here that deletion of the Ogt gene in mouse embryonic stem cells (mESC) results in a widespread increase in the TET product 5-hydroxymethylcytosine (5hmC) in both euchromatic and heterochromatic compartments, with concomitant reduction of the TET substrate 5-methylcytosine (5mC) at the same genomic regions. mESC engineered to abolish the TET1-OGT interaction likewise displayed a genome-wide decrease of 5mC. DNA hypomethylation in OGT-deficient cells was accompanied by de-repression of transposable elements (TEs) predominantly located in heterochromatin, and this increase in TE expression was sometimes accompanied by increased cis -expression of genes and exons located 3' of the expressed TE. Thus, the TET-OGT interaction prevents DNA demethylation and TE expression in heterochromatin by restraining TET activity genome-wide. We suggest that OGT protects the genome against DNA hypomethylation and impaired heterochromatin integrity, preventing the aberrant increase in TE expression observed in cancer, autoimmune-inflammatory diseases, cellular senescence and ageing.

Evolutionary insights from profiling LINE-1 activity at allelic resolution in a single human genome

Transposable elements have created the majority of the sequence in many genomes. In mammals, LINE-1 retrotransposons have been expanding for more than 100 million years as distinct, consecutive lineages; however, the drivers of this recurrent lineage emergence and disappearance are unknown. Most human genome assemblies provide a record of this ancient evolution, but fail to resolve ongoing LINE-1 retrotranspositions. Utilizing the human CHM1 long-read-based haploid assembly, we identified and cloned all full-length, intact LINE-1s, and found 29 LINE-1s with measurable in vitro retrotransposition activity. Among individuals, these LINE-1s varied in their presence, their allelic sequences, and their activity. We found that recently retrotransposed LINE-1s tend to be active in vitro and polymorphic in the population relative to more ancient LINE-1s. However, some rare allelic forms of old LINE-1s retain activity, suggesting older lineages can persist longer than expected. Finally, in LINE-1s with in vitro activity and in vivo fitness, we identified mutations that may have increased replication in ancient genomes and may prove promising candidates for mechanistic investigations of the drivers of LINE-1 evolution and which LINE-1 sequences contribute to human disease.

Generation of somatic de novo structural variation as a hallmark of cellular senescence in human lung fibroblasts

Cellular senescence is characterized by replication arrest in response to stress stimuli. Senescent cells accumulate in aging tissues and can trigger organ-specific and possibly systemic dysfunction. Although senescent cell populations are heterogeneous, a key feature is that they exhibit epigenetic changes. Epigenetic changes such as loss of repressive constitutive heterochromatin could lead to subsequent LINE-1 derepression, a phenomenon often described in the context of senescence or somatic evolution. LINE-1 elements decode the retroposition machinery and reverse transcription generates cDNA from autonomous and non-autonomous TEs that can potentially reintegrate into genomes and cause structural variants. Another feature of cellular senescence is mitochondrial dysfunction caused by mitochondrial damage. In combination with impaired mitophagy, which is characteristic of senescent cells, this could lead to cytosolic mtDNA accumulation and, as a genomic consequence, integrations of mtDNA into nuclear DNA (nDNA), resulting in mitochondrial pseudogenes called numts. Thus, both phenomena could cause structural variants in aging genomes that go beyond epigenetic changes. We therefore compared proliferating and senescent IMR-90 cells in terms of somatic de novo numts and integrations of a non-autonomous composite retrotransposons - the so-called SVA elements-that hijack the retropositional machinery of LINE-1. We applied a subtractive and kinetic enrichment technique using proliferating cell DNA as a driver and senescent genomes as a tester for the detection of nuclear flanks of de novo SVA integrations. Coupled with deep sequencing we obtained a genomic readout for SVA retrotransposition possibly linked to cellular senescence in the IMR-90 model. Furthermore, we compared the genomes of proliferative and senescent IMR-90 cells by deep sequencing or after enrichment of nuclear DNA using AluScan technology. A total of 1,695 de novo SVA integrations were detected in senescent IMR-90 cells, of which 333 were unique. Moreover, we identified a total of 81 de novo numts with perfect identity to both mtDNA and nuclear hg38 flanks. In summary, we present evidence for possible age-dependent structural genomic changes by paralogization that go beyond epigenetic modifications. We hypothesize, that the structural variants we observe potentially impact processes associated with replicative aging of IMR-90 cells.

The epitranscriptome of high-grade gliomas: a promising therapeutic target with implications from the tumor microenvironment to endogenous retroviruses

Glioblastoma (GBM) comprises 45.6% of all primary malignant brain cancers and is one of the most common and aggressive intracranial tumors in adults. Intratumoral heterogeneity with a wide range of proteomic, genetic, and epigenetic dysregulation contributes to treatment resistance and poor prognosis, thus demanding novel therapeutic approaches. To date, numerous clinical trials have been developed to target the proteome and epigenome of high-grade gliomas with promising results. However, studying RNA modifications, or RNA epitranscriptomics, is a new frontier within neuro-oncology. RNA epitranscriptomics was discovered in the 1970s, but in the last decade, the extent of modification of mRNA and various non-coding RNAs has emerged and been implicated in transposable element activation and many other oncogenic processes within the tumor microenvironment. This review provides background information and discusses the therapeutic potential of agents modulating epitranscriptomics in high-grade gliomas. A particular emphasis will be placed on how combination therapies that include immune agents targeting hERV-mediated viral mimicry could improve the treatment of GBM.

Stress Induced Activation of LTR Retrotransposons in the Drosophila melanogaster Genome

Background: Retrotransposons with long terminal repeats (LTR retrotransposons) are widespread in all groups of eukaryotes and are often both the cause of new mutations and the source of new sequences. Apart from their high activity in generative and differentiation-stage tissues, LTR retrotransposons also become more active in response to different stressors. The precise causes of LTR retrotransposons' activation in response to stress, however, have not yet been thoroughly investigated. Methods: We used RT-PCR to investigate the transcriptional profile of LTR retrotransposons and piRNA clusters in response to oxidative and chronic heat stresses. We used Oxford Nanopore sequencing to investigate the genomic environment of new insertions of the retrotransposons. We used bioinformatics methods to find the stress-induced transcription factor binding sites in LTR retrotransposons. Results: We studied the transposition activity and transcription level of LTR retrotransposons in response to oxidative and chronic heat stress and assessed the contribution of various factors that can affect the increase in their expression under stress conditions: the state of the piRNA-interference system, the influence of the genomic environment on individual copies, and the presence of the stress-induced transcription factor binding sites in retrotransposon sequences. Conclusions: The main reason for the activation of LTR retrotransposons under stress conditions is the presence of transcription factor binding sites in their regulatory sequences, which are triggered in response to stress and are necessary for tissue regeneration processes. Stress-induced transposable element activation can function as a trigger mechanism, triggering multiple signal pathways and resulting in a polyvariant cell response.

Detection and annotation of transposable element insertions and deletions on the human genome using nanopore sequencing

Repetitive sequences represent about 45% of the human genome. Some are transposable elements (TEs) with the ability to change their position in the genome, creating genetic variability both as insertions or deletions, with potential pathogenic consequences. We used long-read nanopore sequencing to identify TE variants in the genomes of 24 patients with antithrombin deficiency. We identified 7 344 TE insertions and 3 056 TE deletions, 2 926 were not previously described in publicly available databases. The insertions affected 3 955 genes, with 6 insertions located in exons, 3 929 in introns, and 147 in promoters. Potential functional impact was evaluated with gene annotation and enrichment analysis, which suggested a strong relationship with neuron-related functions and autism. We conclude that this study encourages the generation of a complete map of TEs in the human genome, which will be useful for identifying new TEs involved in genetic disorders.

Comprehensive profiling of L1 retrotransposons in mouse

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

capTEs enables locus-specific dissection of transcriptional outputs from reference and nonreference transposable elements

Transposable elements (TEs) serve as both insertional mutagens and regulatory elements in cells, and their aberrant activity is increasingly being revealed to contribute to diseases and cancers. However, measuring the transcriptional consequences of nonreference and young TEs at individual loci remains challenging with current methods, primarily due to technical limitations, including short read lengths generated and insufficient coverage in target regions. Here, we introduce a long-read targeted RNA sequencing method, Cas9-assisted profiling TE expression sequencing (capTEs), for quantitative analysis of transcriptional outputs for individual TEs, including transcribed nonreference insertions, noncanonical transcripts from various transcription patterns and their correlations with expression changes in related genes. This method selectively identified TE-containing transcripts and outputted data with up to 90% TE reads, maintaining a comparable data yield to whole-transcriptome sequencing. We applied capTEs to human cancer cells and found that internal and inserted Alu elements may employ distinct regulatory mechanisms to upregulate gene expression. We expect that capTEs will be a critical tool for advancing our understanding of the biological functions of individual TEs at the locus level, revealing their roles as both mutagens and regulators in biological and pathogenic processes.

Local read haplotagging enables accurate long-read small variant calling

Long-read sequencing technology has enabled variant detection in difficult-to-map regions of the genome and enabled rapid genetic diagnosis in clinical settings. Rapidly evolving third-generation sequencing platforms like Pacific Biosciences (PacBio) and Oxford nanopore technologies (ONT) are introducing newer platforms and data types. It has been demonstrated that variant calling methods based on deep neural networks can use local haplotyping information with long-reads to improve the genotyping accuracy. However, using local haplotype information creates an overhead as variant calling needs to be performed multiple times which ultimately makes it difficult to extend to new data types and platforms as they get introduced. In this work, we have developed a local haplotype approximate method that enables state-of-the-art variant calling performance with multiple sequencing platforms including PacBio Revio system, ONT R10.4 simplex and duplex data. This addition of local haplotype approximation makes DeepVariant a universal variant calling solution for long-read sequencing platforms.

Genomics in the long-read sequencing era

Long-read sequencing (LRS) technologies have provided extremely powerful tools to explore genomes. While in the early years these methods suffered technical limitations, they have recently made significant progress in terms of read length, throughput, and accuracy and bioinformatics tools have strongly improved. Here, we aim to review the current status of LRS technologies, the development of novel methods, and the impact on genomics research. We will explore the most impactful recent findings made possible by these technologies focusing on high-resolution sequencing of genomes and transcriptomes and the direct detection of DNA and RNA modifications. We will also discuss how LRS methods promise a more comprehensive understanding of human genetic variation, transcriptomics, and epigenetics for the coming years.

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.

Transient naive reprogramming corrects hiPS cells functionally and epigenetically

Cells undergo a major epigenome reconfiguration when reprogrammed to human induced pluripotent stem cells (hiPS cells). However, the epigenomes of hiPS cells and human embryonic stem (hES) cells differ significantly, which affects hiPS cell function1-8. These differences include epigenetic memory and aberrations that emerge during reprogramming, for which the mechanisms remain unknown. Here we characterized the persistence and emergence of these epigenetic differences by performing genome-wide DNA methylation profiling throughout primed and naive reprogramming of human somatic cells to hiPS cells. We found that reprogramming-induced epigenetic aberrations emerge midway through primed reprogramming, whereas DNA demethylation begins early in naive reprogramming. Using this knowledge, we developed a transient-naive-treatment (TNT) reprogramming strategy that emulates the embryonic epigenetic reset. We show that the epigenetic memory in hiPS cells is concentrated in cell of origin-dependent repressive chromatin marked by H3K9me3, lamin-B1 and aberrant CpH methylation. TNT reprogramming reconfigures these domains to a hES cell-like state and does not disrupt genomic imprinting. Using an isogenic system, we demonstrate that TNT reprogramming can correct the transposable element overexpression and differential gene expression seen in conventional hiPS cells, and that TNT-reprogrammed hiPS and hES cells show similar differentiation efficiencies. Moreover, TNT reprogramming enhances the differentiation of hiPS cells derived from multiple cell types. Thus, TNT reprogramming corrects epigenetic memory and aberrations, producing hiPS cells that are molecularly and functionally more similar to hES cells than conventional hiPS cells. We foresee TNT reprogramming becoming a new standard for biomedical and therapeutic applications and providing a novel system for studying epigenetic memory.

Nanopore sequencing technology and its applications

Since the development of Sanger sequencing in 1977, sequencing technology has played a pivotal role in molecular biology research by enabling the interpretation of biological genetic codes. Today, nanopore sequencing is one of the leading third-generation sequencing technologies. With its long reads, portability, and low cost, nanopore sequencing is widely used in various scientific fields including epidemic prevention and control, disease diagnosis, and animal and plant breeding. Despite initial concerns about high error rates, continuous innovation in sequencing platforms and algorithm analysis technology has effectively addressed its accuracy. During the coronavirus disease (COVID-19) pandemic, nanopore sequencing played a critical role in detecting the severe acute respiratory syndrome coronavirus-2 virus genome and containing the pandemic. However, a lack of understanding of this technology may limit its popularization and application. Nanopore sequencing is poised to become the mainstream choice for preventing and controlling COVID-19 and future epidemics while creating value in other fields such as oncology and botany. This work introduces the contributions of nanopore sequencing during the COVID-19 pandemic to promote public understanding and its use in emerging outbreaks worldwide. We discuss its application in microbial detection, cancer genomes, and plant genomes and summarize strategies to improve its accuracy.