Nanopore Sequencing Enables Comprehensive Transposable Element Epigenomic Profiling

Whole-genome analysis reveals the contribution of non-coding de novo transposon insertions to autism spectrum disorder

BACKGROUND

Retrotransposons have been implicated as causes of Mendelian disease, but their role in autism spectrum disorder (ASD) has not been systematically defined, because they are only called with adequate sensitivity from whole genome sequencing (WGS) data and a large enough cohort for this analysis has only recently become available.

RESULTS

We analyzed WGS data from a cohort of 2288 ASD families from the Simons Simplex Collection by establishing a scalable computational pipeline for retrotransposon insertion detection. We report 86,154 polymorphic retrotransposon insertions-including > 60% not previously reported-and 158 de novo retrotransposition events. The overall burden of de novo events was similar between ASD individuals and unaffected siblings, with 1 de novo insertion per 29, 117, and 206 births for Alu, L1, and SVA respectively, and 1 de novo insertion per 21 births total. However, ASD cases showed more de novo L1 insertions than expected in ASD genes. Additionally, we observed exonic insertions in loss-of-function intolerant genes, including a likely pathogenic exonic insertion in CSDE1, only in ASD individuals.

CONCLUSIONS

These findings suggest a modest, but important, impact of intronic and exonic retrotransposon insertions in ASD, show the importance of WGS for their analysis, and highlight the utility of specific bioinformatic tools for high-throughput detection of retrotransposon insertions.

Nanopore sequencing technology, bioinformatics and applications

Rapid advances in nanopore technologies for sequencing single long DNA and RNA molecules have led to substantial improvements in accuracy, read length and throughput. These breakthroughs have required extensive development of experimental and bioinformatics methods to fully exploit nanopore long reads for investigations of genomes, transcriptomes, epigenomes and epitranscriptomes. Nanopore sequencing is being applied in genome assembly, full-length transcript detection and base modification detection and in more specialized areas, such as rapid clinical diagnoses and outbreak surveillance. Many opportunities remain for improving data quality and analytical approaches through the development of new nanopores, base-calling methods and experimental protocols tailored to particular applications.

Beyond sequencing: machine learning algorithms extract biology hidden in Nanopore signal data

Nanopore sequencing provides signal data corresponding to the nucleotide motifs sequenced. Through machine learning-based methods, these signals are translated into long-read sequences that overcome the read size limit of short-read sequencing. However, analyzing the raw nanopore signal data provides many more opportunities beyond just sequencing genomes and transcriptomes: algorithms that use machine learning approaches to extract biological information from these signals allow the detection of DNA and RNA modifications, the estimation of poly(A) tail length, and the prediction of RNA secondary structures. In this review, we discuss how developments in machine learning methodologies contributed to more accurate basecalling and lower error rates, and how these methods enable new biological discoveries. We argue that direct nanopore sequencing of DNA and RNA provides a new dimensionality for genomics experiments and highlight challenges and future directions for computational approaches to extract the additional information provided by nanopore signal data.

DNA methylation-calling tools for Oxford Nanopore sequencing: a survey and human epigenome-wide evaluation

BACKGROUND

Nanopore long-read sequencing technology greatly expands the capacity of long-range, single-molecule DNA-modification detection. A growing number of analytical tools have been developed to detect DNA methylation from nanopore sequencing reads. Here, we assess the performance of different methylation-calling tools to provide a systematic evaluation to guide researchers performing human epigenome-wide studies.

RESULTS

We compare seven analytic tools for detecting DNA methylation from nanopore long-read sequencing data generated from human natural DNA at a whole-genome scale. We evaluate the per-read and per-site performance of CpG methylation prediction across different genomic contexts, CpG site coverage, and computational resources consumed by each tool. The seven tools exhibit different performances across the evaluation criteria. We show that the methylation prediction at regions with discordant DNA methylation patterns, intergenic regions, low CG density regions, and repetitive regions show room for improvement across all tools. Furthermore, we demonstrate that 5hmC levels at least partly contribute to the discrepancy between bisulfite and nanopore sequencing. Lastly, we provide an online DNA methylation database ( https://nanome.jax.org ) to display the DNA methylation levels detected by nanopore sequencing and bisulfite sequencing data across different genomic contexts.

CONCLUSIONS

Our study is the first systematic benchmark of computational methods for detection of mammalian whole-genome DNA modifications in nanopore sequencing. We provide a broad foundation for cross-platform standardization and an evaluation of analytical tools designed for genome-scale modified base detection using nanopore sequencing.

Primate-specific retrotransposons and the evolution of circadian networks in the human brain

The circadian rhythm of the human brain is attuned to sleep-wake cycles that entail global alterations in neuronal excitability. This periodicity involves a highly coordinated regulation of gene expression. A growing number of studies are documenting a fascinating connection between primate-specific retrotransposons (Alu elements) and key epigenetic regulatory processes in the primate brain. Collectively, these studies indicate that Alu elements embedded in the human neuronal genome mediate post-transcriptional processes that unite human-specific neuroepigenetic landscapes and circadian rhythm. Here, we review evidence linking Alu retrotransposon-mediated posttranscriptional pathways to circadian gene expression. We hypothesize that Alu retrotransposons participate in the organization of circadian brain function through multidimensional neuroepigenetic pathways. We anticipate that these pathways are closely tied to the evolution of human cognition and their perturbation contributes to the manifestation of human-specific neurological diseases. Finally, we address current challenges and accompanying opportunities in studying primate- and human-specific transposable elements.

Characterization of full-length LINE-1 insertions in 154 genomes

Long interspersed nuclear elements (LINEs) are retrotransposons that contribute to genetic variation in the human genome. LINE-1 elements in larger-scale studies are challenging to identify using sequencing technologies due to cost and scalability. We developed an approach using optical mapping for detection of full-length LINE-1 insertions and 10× sequencing for confirmation. We found 51 true positive full-length LINE-1 insertions, of which 4 are novel insertions, in NA12878. Repeating our analysis on a larger sample set representing 26 populations, we identified 329 full-length LINE-1 elements, of which 123 are novel. 24.8% of these 329 LINE-1 insertions were shared amongst all 5 superpopulations (AFR, AMR, EUR, EAS, SAS). The African superpopulation has a higher percentage of population-specific LINE-1 insertions than any other superpopulation. These data indicate that our approach can provide high-speed, cost-effective, and increased accuracy for LINE-1 detection. These data also provide an insight into variations of LINE-1 elements between different populations.

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

BACKGROUND

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

RESULTS

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

CONCLUSION

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

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

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

DNA methylation detection technology and plasma-based methylation biomarkers in screening of gastrointestinal carcinoma

DNA methylation is of paramount importance for the evolution of human cancers. Its high sensitivity and specificity make it a potential biomarker for early cancer screening in the context of an increasing global burden of gastrointestinal (GI) carcinoma. More DNA methylation biomarkers are emerging with the development of liquid biopsy and sensitive DNA methylation detection technology. This review provides an overview of DNA methylation, focusing on the presentation and comparison of 5-methylcytosine detection technologies, and introduces the promising plasma-based cell-free DNA (cfDNA) methylation biomarkers published in recent years for early screening of GI carcinoma. Finally, we summarize and discuss the future of plasma cfDNA methylation markers detection as a clinical tool for early screening of GI carcinoma.

Detection of Cell Types Contributing to Cancer From Circulating, Cell-Free Methylated DNA

Detection of cellular changes in tissue biopsies has been the basis for cancer diagnostics. However, tissue biopsies are invasive and limited by inaccuracies due to sampling locations, restricted sampling frequency, and poor representation of tissue heterogeneity. Liquid biopsies are emerging as a complementary approach to traditional tissue biopsies to detect dynamic changes in specific cell populations. Cell-free DNA (cfDNA) fragments released into the circulation from dying cells can be traced back to the tissues and cell types they originated from using DNA methylation, an epigenetic regulatory mechanism that is highly cell-type specific. Decoding changes in the cellular origins of cfDNA over time can reveal altered host tissue homeostasis due to local cancer invasion and metastatic spread to distant organs as well as treatment responses. In addition to host-derived cfDNA, changes in cancer cells can be detected from cell-free, circulating tumor DNA (ctDNA) by monitoring DNA mutations carried by cancer cells. Here, we will discuss computational approaches to identify and validate robust biomarkers of changed tissue homeostasis using cell-free, methylated DNA in the circulation. We highlight studies performing genome-wide profiling of cfDNA methylation and those that combine genetic and epigenetic markers to further identify cell-type specific signatures. Finally, we discuss opportunities and current limitations of these approaches for implementation in clinical oncology.

Roles of TET and TDG in DNA demethylation in proliferating and non-proliferating immune cells

BACKGROUND

TET enzymes mediate DNA demethylation by oxidizing 5-methylcytosine (5mC) in DNA to 5-hydroxymethylcytosine (5hmC), 5-formylcytosine (5fC), and 5-carboxylcytosine (5caC). Since these oxidized methylcytosines (oxi-mCs) are not recognized by the maintenance methyltransferase DNMT1, DNA demethylation can occur through "passive," replication-dependent dilution when cells divide. A distinct, replication-independent ("active") mechanism of DNA demethylation involves excision of 5fC and 5caC by the DNA repair enzyme thymine DNA glycosylase (TDG), followed by base excision repair.

RESULTS

Here by analyzing inducible gene-disrupted mice, we show that DNA demethylation during primary T cell differentiation occurs mainly through passive replication-dependent dilution of all three oxi-mCs, with only a negligible contribution from TDG. In addition, by pyridine borane sequencing (PB-seq), a simple recently developed method that directly maps 5fC/5caC at single-base resolution, we detect the accumulation of 5fC/5caC in TDG-deleted T cells. We also quantify the occurrence of concordant demethylation within and near enhancer regions in the Il4 locus. In an independent system that does not involve cell division, macrophages treated with liposaccharide accumulate 5hmC at enhancers and show altered gene expression without DNA demethylation; loss of TET enzymes disrupts gene expression, but loss of TDG has no effect. We also observe that mice with long-term (1 year) deletion of Tdg are healthy and show normal survival and hematopoiesis.

CONCLUSIONS

We have quantified the relative contributions of TET and TDG to cell differentiation and DNA demethylation at representative loci in proliferating T cells. We find that TET enzymes regulate T cell differentiation and DNA demethylation primarily through passive dilution of oxi-mCs. In contrast, while we observe a low level of active, replication-independent DNA demethylation mediated by TDG, this process does not appear to be essential for immune cell activation or differentiation.

Cas9 targeted enrichment of mobile elements using nanopore sequencing

Mobile element insertions (MEIs) are repetitive genomic sequences that contribute to genetic variation and can lead to genetic disorders. Targeted and whole-genome approaches using short-read sequencing have been developed to identify reference and non-reference MEIs; however, the read length hampers detection of these elements in complex genomic regions. Here, we pair Cas9-targeted nanopore sequencing with computational methodologies to capture active MEIs in human genomes. We demonstrate parallel enrichment for distinct classes of MEIs, averaging 44% of reads on-targeted signals and exhibiting a 13.4-54x enrichment over whole-genome approaches. We show an individual flow cell can recover most MEIs (97% L1Hs, 93% AluYb, 51% AluYa, 99% SVA_F, and 65% SVA_E). We identify seventeen non-reference MEIs in GM12878 overlooked by modern, long-read analysis pipelines, primarily in repetitive genomic regions. This work introduces the utility of nanopore sequencing for MEI enrichment and lays the foundation for rapid discovery of elusive, repetitive genetic elements.

Cancer Biomarkers Discovery of Methylation Modification With Direct High-Throughput Nanopore Sequencing

Cancer is a complex disease, driven by a combination of genetic and epigenetic alterations. DNA and RNA methylation modifications are the most common epigenetic events that play critical roles in cancer development and progression. Bisulfite converted sequencing is a widely used technique to detect base modifications in DNA methylation, but its main drawbacks lie in DNA degradation, lack of specificity, or short reads with low sequence diversity. The nanopore sequencing technology can directly detect base modifications in native DNA as well as RNA without harsh chemical treatment, compared to bisulfite sequencing. Furthermore, CRISPR/Cas9-targeted enrichment nanopore sequencing techniques are straightforward and cost-effective when targeting genomic regions are of interest. In this review, we mainly focus on DNA and RNA methylation modification detection in cancer with the current nanopore sequencing approaches. We also present the respective strengths, weaknesses of nanopore sequencing techniques, and their future translational applications in identification of epigenetic biomarkers for cancer detection and prognosis.

Endogenous retroviruses in the origins and treatment of cancer

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

Unraveling the features of somatic transposition in the Drosophila intestine

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

Virus-derived variation in diverse human genomes

Acquisition of genetic material from viruses by their hosts can generate inter-host structural genome variation. We developed computational tools enabling us to study virus-derived structural variants (SVs) in population-scale whole genome sequencing (WGS) datasets and applied them to 3,332 humans. Although SVs had already been cataloged in these subjects, we found previously-overlooked virus-derived SVs. We detected non-germline SVs derived from squirrel monkey retrovirus (SMRV), human immunodeficiency virus 1 (HIV-1), and human T lymphotropic virus (HTLV-1); these variants are attributable to infection of the sequenced lymphoblastoid cell lines (LCLs) or their progenitor cells and may impact gene expression results and the biosafety of experiments using these cells. In addition, we detected new heritable SVs derived from human herpesvirus 6 (HHV-6) and human endogenous retrovirus-K (HERV-K). We report the first solo-direct repeat (DR) HHV-6 likely to reflect DR rearrangement of a known full-length endogenous HHV-6. We used linkage disequilibrium between single nucleotide variants (SNVs) and variants in reads that align to HERV-K, which often cannot be mapped uniquely using conventional short-read sequencing analysis methods, to locate previously-unknown polymorphic HERV-K loci. Some of these loci are tightly linked to trait-associated SNVs, some are in complex genome regions inaccessible by prior methods, and some contain novel HERV-K haplotypes likely derived from gene conversion from an unknown source or introgression. These tools and results broaden our perspective on the coevolution between viruses and humans, including ongoing virus-to-human gene transfer contributing to genetic variation between humans.

Variable number tandem repeats - Their emerging role in sickness and health

Understanding the mechanisms regulating tissue specific and stimulus inducible regulation is at the heart of understanding human biology and how this translates to wellbeing, the ageing process, and disease progression. Polymorphic DNA variation is superimposed as an extra layer of complexity in such processes which underpin our individuality and are the focus of personalized medicine. This review focuses on the role and action of repetitive DNA, specifically variable number tandem repeats and SINE-VNTR-Alu domains, highlighting their role in modification of gene structure and gene expression in addition to their polymorphic nature being a genetic modifier of disease risk and progression. Although the literature focuses on their role in disease, it illustrates their potential to be major contributors to normal physiological function. To date, these elements have been under-reported in genomic analysis due to the difficulties in their characterization with short read DNA sequencing methods. However, recent advances in long read sequencing methods should resolve these problems allowing for a greater understanding of their contribution to a host of genomic and functional mechanisms underlying physiology and disease.