Mapping DNA methylation with high-throughput nanopore sequencing

Fc gamma receptors: Their evolution, genomic architecture, genetic variation, and impact on human disease

Fc gamma receptors (FcγRs) are a family of receptors that bind IgG antibodies and interface at the junction of humoral and innate immunity. Precise regulation of receptor expression provides the necessary balance to achieve healthy immune homeostasis by establishing an appropriate immune threshold to limit autoimmunity but respond effectively to infection. The underlying genetics of the FCGR gene family are central to achieving this immune threshold by regulating affinity for IgG, signaling efficacy, and receptor expression. The FCGR gene locus was duplicated during evolution, retaining very high homology and resulting in a genomic region that is technically difficult to study. Here, we review the recent evolution of the gene family in mammals, its complexity and variation through copy number variation and single-nucleotide polymorphism, and impact of these on disease incidence, resolution, and therapeutic antibody efficacy. We also discuss the progress and limitations of current approaches to study the region and emphasize how new genomics technologies will likely resolve much of the current confusion in the field. This will lead to definitive conclusions on the impact of genetic variation within the FCGR gene locus on immune function and disease.

A strategy for studying epigenetic diversity in natural populations: proof of concept in poplar and oak

In the last 20 years, several techniques have been developed for quantifying DNA methylation, the most studied epigenetic marks in eukaryotes, including the gold standard method, whole-genome bisulfite sequencing (WGBS). WGBS quantifies genome-wide DNA methylation but has several inconveniences rendering it less suitable for population-scale epigenetic studies. The high cost of deep sequencing and the large amounts of data generated prompted us to seek an alternative approach. Restricting studies to parts of the genome would be a satisfactory alternative had there not been a major limitation: the need to select upstream targets corresponding to differentially methylated regions as targets. Given the need to study large numbers of samples, we propose a strategy for investigating DNA methylation variation in natural populations, taking into account the structural complexity of genomes, their size, and their content in unique coding regions versus repeated regions as transposable elements. We first identified regions of highly variable DNA methylation in a subset of genotypes representative of the biological diversity in the population by WGBS. We then analysed the variations of DNA methylation in these targeted regions at the population level by sequencing capture bisulfite (SeqCapBis). The entire strategy was then validated by applying it to another species. Our strategy was developed as a proof of concept on natural populations of two forest species: Populus nigra and Quercus petraea.

Nanopore sequencing demonstrates the roles of spermatozoal DNA N6-methyladenine in mediating transgenerational lipid metabolism disorder induced by excessive folate consumpton

Increased consumption of folic acid is prevalent due to its beneficial effects, but growing evidence emphasizes the side effects pointing to excessive dietary folate intake. The effects of excessive paternal folic acid consumption on offspring and its transgenerational inheritance mechanism have not been elucidated. We hypothesize that excessive folic acid consumption will alter sperm DNA N6-methyladenine (6mA) and 5-methylcytosine (5mC) methylation and heritably influence offspring metabolic homeostasis. Here, we fed roosters either folic acid-control or folic acid-excess diet throughout life. Paternal chronic folic acid excessive supplementation increased hepatic lipogenesis and lipid accumulation but reduced lipolysis both in the roosters and their offspring, which was further confirmed to be induced by one-carbon metabolism inhibition and gene expression alteration associated with the Peroxisome proliferator-activated receptor pathway. Based on the spermatozoal genome-wide DNA methylome identified by Nanopore sequencing, multi-omics association analysis of spermatozoal and hepatic DNA methylome, transcriptome, and metabolome suggested that differential spermatozoal DNA 6mA and 5mC methylation could be involved in regulating lipid metabolism-related gene expression in offspring chickens. This model suggests that sperm DNA N6-methyladenine and 5-methylcytosine methylation were involved in epigenetic transmission and that paternal dietary excess folic acid leads to hepatic lipid accumulation in offspring.

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.

Adapting nanopore sequencing basecalling models for modification detection via incremental learning and anomaly detection

We leverage machine learning approaches to adapt nanopore sequencing basecallers for nucleotide modification detection. We first apply the incremental learning (IL) technique to improve the basecalling of modification-rich sequences, which are usually of high biological interest. With sequence backbones resolved, we further run anomaly detection (AD) on individual nucleotides to determine their modification status. By this means, our pipeline promises the single-molecule, single-nucleotide, and sequence context-free detection of modifications. We benchmark the pipeline using control oligos, further apply it in the basecalling of densely-modified yeast tRNAs and E.coli genomic DNAs, the cross-species detection of N6-methyladenosine (m6A) in mammalian mRNAs, and the simultaneous detection of N1-methyladenosine (m1A) and m6A in human mRNAs. Our IL-AD workflow is available at: https://github.com/wangziyuan66/IL-AD .

The Ruminant Telomere-to-Telomere (RT2T) Consortium

Telomere-to-telomere (T2T) assemblies reveal new insights into the structure and function of the previously 'invisible' parts of the genome and allow comparative analyses of complete genomes across entire clades. We present here an open collaborative effort, termed the 'Ruminant T2T Consortium' (RT2T), that aims to generate complete diploid assemblies for numerous species of the Artiodactyla suborder Ruminantia to examine chromosomal evolution in the context of natural selection and domestication of species used as livestock.

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.

Simultaneous assessment of human genome and methylome data in a single experiment using limited deamination of methylated cytosine

Multiomics require concerted recording of independent information, ideally from a single experiment. In this study, we introduce RIMS-seq2, a high-throughput technique to simultaneously sequence genomes and overlay methylation information while requiring only a small modification of the experimental protocol for high-throughput DNA sequencing to include a controlled deamination step. Importantly, the rate of deamination of 5-methylcytosine is negligible and thus does not interfere with standard DNA sequencing and data processing. Thus, RIMS-seq2 libraries from whole- or targeted-genome sequencing show the same germline variation calling accuracy and sensitivity compared with standard DNA-seq. Additionally, regional methylation levels provide an accurate map of the human methylome.

Emerging methylation-based approaches in microbiome engineering

Bacterial epigenetics, particularly through DNA methylation, exerts significant influence over various biological processes such as DNA replication, uptake, and gene regulation in bacteria. In this review, we explore recent advances in characterizing bacterial epigenomes, accompanied by emerging strategies that harness bacterial epigenetics to elucidate and engineer diverse bacterial species with precision and effectiveness. Furthermore, we delve into the potential of epigenetic modifications to steer microbial functions and influence community dynamics, offering promising opportunities for understanding and modulating microbiomes. Additionally, we investigate the extensive diversity of DNA methyltransferases and emphasize their potential utility in the context of the human microbiome. In summary, this review highlights the potential of DNA methylation as a powerful toolkit for engineering microbiomes.

Rockfish: A transformer-based model for accurate 5-methylcytosine prediction from nanopore sequencing

DNA methylation plays an important role in various biological processes, including cell differentiation, ageing, and cancer development. The most important methylation in mammals is 5-methylcytosine mostly occurring in the context of CpG dinucleotides. Sequencing methods such as whole-genome bisulfite sequencing successfully detect 5-methylcytosine DNA modifications. However, they suffer from the serious drawbacks of short read lengths and might introduce an amplification bias. Here we present Rockfish, a deep learning algorithm that significantly improves read-level 5-methylcytosine detection by using Nanopore sequencing. Rockfish is compared with other methods based on Nanopore sequencing on R9.4.1 and R10.4.1 datasets. There is an increase in the single-base accuracy and the F1 measure of up to 5 percentage points on R.9.4.1 datasets, and up to 0.82 percentage points on R10.4.1 datasets. Moreover, Rockfish shows a high correlation with whole-genome bisulfite sequencing, requires lower read depth, and achieves higher confidence in biologically important regions such as CpG-rich promoters while being computationally efficient. Its superior performance in human and mouse samples highlights its versatility for studying 5-methylcytosine methylation across varied organisms and diseases. Finally, its adaptable architecture ensures compatibility with new versions of pores and chemistry as well as modification types.

Single-stranded pre-methylated 5mC adapters uncover the methylation profile of plasma ultrashort Single-stranded cell-free DNA

Whole-genome bisulfite sequencing (BS-Seq) measures cytosine methylation changes at single-base resolution and can be used to profile cell-free DNA (cfDNA). In plasma, ultrashort single-stranded cfDNA (uscfDNA, ∼50 nt) has been identified together with 167 bp double-stranded mononucleosomal cell-free DNA (mncfDNA). However, the methylation profile of uscfDNA has not been described. Conventional BS-Seq workflows may not be helpful because bisulfite conversion degrades larger DNA into smaller fragments, leading to erroneous categorization as uscfDNA. We describe the '5mCAdpBS-Seq' workflow in which pre-methylated 5mC (5-methylcytosine) single-stranded adapters are ligated to heat-denatured cfDNA before bisulfite conversion. This method retains only DNA fragments that are unaltered by bisulfite treatment, resulting in less biased uscfDNA methylation analysis. Using 5mCAdpBS-Seq, uscfDNA had lower levels of DNA methylation (∼15%) compared to mncfDNA and was enriched in promoters and CpG islands. Hypomethylated uscfDNA fragments were enriched in upstream transcription start sites (TSSs), and the intensity of enrichment was correlated with expressed genes of hemopoietic cells. Using tissue-of-origin deconvolution, we inferred that uscfDNA is derived primarily from eosinophils, neutrophils, and monocytes. As proof-of-principle, we show that characteristics of the methylation profile of uscfDNA can distinguish non-small cell lung carcinoma from non-cancer samples. The 5mCAdpBS-Seq workflow is recommended for any cfDNA methylation-based investigations.

Effective training of nanopore callers for epigenetic marks with limited labelled data

Nanopore sequencing platforms combined with supervised machine learning (ML) have been effective at detecting base modifications in DNA such as 5-methylcytosine (5mC) and N6-methyladenine (6mA). These ML-based nanopore callers have typically been trained on data that span all modifications on all possible DNA [Formula: see text]-mer backgrounds-a complete training dataset. However, as nanopore technology is pushed to more and more epigenetic modifications, such complete training data will not be feasible to obtain. Nanopore calling has historically been performed with hidden Markov models (HMMs) that cannot make successful calls for [Formula: see text]-mer contexts not seen during training because of their independent emission distributions. However, deep neural networks (DNNs), which share parameters across contexts, are increasingly being used as callers, often outperforming their HMM cousins. It stands to reason that a DNN approach should be able to better generalize to unseen [Formula: see text]-mer contexts. Indeed, herein we demonstrate that a common DNN approach (DeepSignal) outperforms a common HMM approach (Nanopolish) in the incomplete data setting. Furthermore, we propose a novel hybrid HMM-DNN approach, amortized-HMM, that outperforms both the pure HMM and DNN approaches on 5mC calling when the training data are incomplete. This type of approach is expected to be useful for calling other base modifications such as 5-hydroxymethylcytosine and for the simultaneous calling of different modifications, settings in which complete training data are not likely to be available.

Identification of transcriptionally-active human papillomavirus integrants through nanopore sequencing reveals viable targets for gene therapy against cervical cancer

Integration of the human papillomavirus (HPV) genome into the cellular genome is a key event that leads to constitutive expression of viral oncoprotein E6/E7 and drives the progression of cervical cancer. However, HPV integration patterns differ on a case-by-case basis among related malignancies. Next-generation sequencing technologies still face challenges for interrogating HPV integration sites. In this study, utilizing Nanopore long-read sequencing, we identified 452 and 108 potential integration sites from the cervical cancer cell lines (CaSki and HeLa) and five tissue samples, respectively. Based on long Nanopore chimeric reads, we were able to analyze the methylation status of the HPV long control region (LCR), which controls oncogene E6/E7 expression, and to identify transcriptionally-active integrants among the numerous integrants. As a proof of concept, we identified an active HPV integrant in between RUNX2 and CLIC5 on chromosome 6 in the CaSki cell line, which was supported by ATAC-seq, H3K27Ac ChIP-seq, and RNA-seq analysis. Knockout of the active HPV integrant, by the CRISPR/Cas9 system, dramatically crippled cell proliferation and induced cell senescence. In conclusion, identifying transcriptionally-active HPV integrants with Nanopore sequencing can provide viable targets for gene therapy against HPV-associated cancers.

DNA modifications: Biomarkers for the exposome?

The concept of the exposome is the encompassing of all the environmental exposures, both exogenous and endogenous, across the life course. Many, if not all, of these exposures can result in the generation of reactive species, and/or the modulation of cellular processes, that can lead to a breadth of modifications of DNA, the nature of which may be used to infer their origin. Because of their role in cell function, such modifications have been associated with various major human diseases, including cancer, and so their assessment is crucial. Historically, most methods have been able to only measure one or a few DNA modifications at a time, limiting the information available. With the development of DNA adductomics, which aims to determine the totality of DNA modifications, a far more comprehensive picture of the DNA adduct burden can be gained. Importantly, DNA adductomics can facilitate a "top-down" investigative approach whereby patterns of adducts may be used to trace and identify the originating exposure source. This, together with other 'omic approaches, represents a major tool for unraveling the complexities of the exposome and hence allow a better a understanding of the environmental origins of disease.

Comprehensive genetic analysis of facioscapulohumeral muscular dystrophy by Nanopore long-read whole-genome sequencing

BACKGROUND

Facioscapulohumeral muscular dystrophy (FSHD) is a high-prevalence autosomal dominant neuromuscular disease characterized by significant clinical and genetic heterogeneity. Genetic diagnosis of FSHD remains a challenge because it cannot be detected by standard sequencing methods and requires a complex diagnosis workflow.

METHODS

We developed a comprehensive genetic FSHD detection method based on Oxford Nanopore Technologies (ONT) whole-genome sequencing. Using a case-control design, we applied this procedure to 29 samples and compared the results with those from optical genome mapping (OGM), bisulfite sequencing (BSS), and whole-exome sequencing (WES).

RESULTS

Using our ONT-based method, we identified 59 haplotypes (35 4qA and 24 4qB) among the 29 samples (including a mosaic sample), as well as the number of D4Z4 repeat units (RUs). The pathogenetic D4Z4 RU contraction identified by our ONT-based method showed 100% concordance with OGM results. The methylation levels of the most distal D4Z4 RU and the double homeobox 4 gene (DUX4) detected by ONT sequencing are highly consistent with the BSS results and showed excellent diagnostic efficiency. Additionally, our ONT-based method provided an independent methylation profile analysis of two permissive 4qA alleles, reflecting a more accurate scenario than traditional BSS. The ONT-based method detected 17 variations in three FSHD2-related genes from nine samples, showing 100% concordance with WES.

CONCLUSIONS

Our ONT-based FSHD detection method is a comprehensive method for identifying pathogenetic D4Z4 RU contractions, methylation level alterations, allele-specific methylation of two 4qA haplotypes, and variations in FSHD2-related genes, which will all greatly improve genetic testing for FSHD.

The fungicide pyraclostrobin affects gene expression by altering the DNA methylation pattern in Magnaporthe oryzae

Introduction

Rice blast disease caused by Magnaporthe oryzae has long been the main cause of rice (Oryza sativa L.) yield reduction worldwide. The quinone external inhibitor pyraclostrobin is widely used as a fungicide to effectively control the spread of pathogenic fungi, including M. oryzae. However, M. oryzae can develop resistance through multiple levels of mutation, such as target protein cytb mutation G143A/S, leading to a decrease in the effectiveness of the biocide after a period of application. Therefore, uncovering the possible mutational mechanisms from multiple perspectives will further provide feasible targets for drug development.

Methods

In this work, we determined the gene expression changes in M. oryzae in response to pyraclostrobin stress and their relationship with DNA methylation by transcriptome and methylome.

Results

The results showed that under pyraclostrobin treatment, endoplasmic reticulum (ER)-associated and ubiquitin-mediated proteolysis were enhanced, suggesting that more aberrant proteins may be generated that need to be cleared. DNA replication and repair processes were inhibited. Glutathione metabolism was enhanced, while lipid metabolism was impaired. The number of alternative splicing events increased. These changes may be related to the elevated methylation levels of cytosine and adenine in gene bodies. Both hypermethylation and hypomethylation of differentially methylated genes (DMGs) mainly occurred in exons and promoters. Some DMGs and differentially expressed genes (DEGs) were annotated to the same pathways by GO and KEGG, including protein processing in the ER, ubiquitin-mediated proteolysis, RNA transport and glutathione metabolism, suggesting that pyraclostrobin may affect gene expression by altering the methylation patterns of cytosine and adenine.

Discussion

Our results revealed that 5mC and 6mA in the gene body are associated with gene expression and contribute to adversity adaptation in M. oryzae. This enriched the understanding for potential mechanism of quinone inhibitor resistance, which will facilitate the development of feasible strategies for maintaining the high efficacy of this kind of fungicide.

Population epigenetics: DNA methylation in the plant omics era

DNA methylation plays an important role in many biological processes. The mechanisms underlying the establishment and maintenance of DNA methylation are well understood thanks to decades of research using DNA methylation mutants, primarily in Arabidopsis (Arabidopsis thaliana) accession Col-0. Recent genome-wide association studies (GWASs) using the methylomes of natural accessions have uncovered a complex and distinct genetic basis of variation in DNA methylation at the population level. Sequencing following bisulfite treatment has served as an excellent method for quantifying DNA methylation. Unlike studies focusing on specific accessions with reference genomes, population-scale methylome research often requires an additional round of sequencing beyond obtaining genome assemblies or genetic variations from whole-genome sequencing data, which can be cost prohibitive. Here, we provide an overview of recently developed bisulfite-free methods for quantifying methylation and cost-effective approaches for the simultaneous detection of genetic and epigenetic information. We also discuss the plasticity of DNA methylation in a specific Arabidopsis accession, the contribution of DNA methylation to plant adaptation, and the genetic determinants of variation in DNA methylation in natural populations. The recently developed technology and knowledge will greatly benefit future studies in population epigenomes.

Mini-review of DNA Methylation Detection Techniques and Their Potential Applications in Disease Diagnosis, Prognosis, and Treatment

DNA methylation is the dominant epigenetic mechanism for regulating gene expression in mammals, playing crucial roles in development, differentiation, and tissue homeostasis. Aberrations in DNA methylation are closely associated with the potential onset of various diseases. Consequently, numerous DNA methylation detection techniques have been successively developed. These methods not only facilitate the exploration of disease mechanisms but also hold significant promise for the development of diagnostic and prognostic strategies. In this Perspective, we present a comprehensive overview of commonly employed DNA methylation detection techniques as well as biosensing based on their underlying analytical techniques. For its medical applications, we begin by examining the pathogenesis of different diseases and then proceed to discuss how relevant technologies are applied in the context of these specific medical conditions. Additionally, we briefly discuss the current limitations of these techniques and highlight future challenges in advancing methylation detection and analysis methodologies.

Comparison of CcrM-dependent methylation in Caulobacter crescentus and Brucella abortus by nanopore sequencing

Bacteria rely on DNA methylation for restriction-modification systems and epigenetic control of gene expression. Here, we use direct detection of methylated bases by nanopore sequencing to monitor global DNA methylation in Alphaproteobacteria, where use of this technique has not yet been reported. One representative of this order, Caulobacter crescentus, relies on DNA methylation to control cell cycle progression, but it is unclear whether other members of this order, such as Brucella abortus, depend on the same systems. We addressed these questions by first measuring CcrM-dependent DNA methylation in Caulobacter and show excellent correlation between nanopore-based detection and previously published results. We then directly measure the impact of Lon-mediated CcrM degradation on the epigenome, verifying that loss of Lon results in pervasive methylation. We also show that the AlkB demethylase has no global impact on DNA methylation during normal growth. Next, we report on the global DNA methylation in Brucella abortus for the first time and find that CcrM-dependent methylation is reliant on Lon but impacts the two chromosomes differently. Finally, we explore the impact of the MucR transcription factor, known to compete with CcrM methylation, on the Brucella methylome and share the results with a publicly available visualization package. Our work demonstrates the utility of nanopore-based sequencing for epigenome measurements in Alphaproteobacteria and reveals new features of CcrM-dependent methylation in a zoonotic pathogen.

When Livestock Genomes Meet Third-Generation Sequencing Technology: From Opportunities to Applications

Third-generation sequencing technology has found widespread application in the genomic, transcriptomic, and epigenetic research of both human and livestock genetics. This technology offers significant advantages in the sequencing of complex genomic regions, the identification of intricate structural variations, and the production of high-quality genomes. Its attributes, including long sequencing reads, obviation of PCR amplification, and direct determination of DNA/RNA, contribute to its efficacy. This review presents a comprehensive overview of third-generation sequencing technologies, exemplified by single-molecule real-time sequencing (SMRT) and Oxford Nanopore Technology (ONT). Emphasizing the research advancements in livestock genomics, the review delves into genome assembly, structural variation detection, transcriptome sequencing, and epigenetic investigations enabled by third-generation sequencing. A comprehensive analysis is conducted on the application and potential challenges of third-generation sequencing technology for genome detection in livestock. Beyond providing valuable insights into genome structure analysis and the identification of rare genes in livestock, the review ventures into an exploration of the genetic mechanisms underpinning exemplary traits. This review not only contributes to our understanding of the genomic landscape in livestock but also provides fresh perspectives for the advancement of research in this domain.

Extraction and selection of high-molecular-weight DNA for long-read sequencing from Chlamydomonas reinhardtii

Recent advances in long-read sequencing technologies have enabled the complete assembly of eukaryotic genomes from telomere to telomere by allowing repeated regions to be fully sequenced and assembled, thus filling the gaps left by previous short-read sequencing methods. Furthermore, long-read sequencing can also help characterizing structural variants, with applications in the fields of genome evolution or cancer genomics. For many organisms, the main bottleneck to sequence long reads remains the lack of robust methods to obtain high-molecular-weight (HMW) DNA. For this purpose, we developed an optimized protocol to extract DNA suitable for long-read sequencing from the unicellular green alga Chlamydomonas reinhardtii, based on CTAB/phenol extraction followed by a size selection step for long DNA molecules. We provide validation results for the extraction protocol, as well as statistics obtained with Oxford Nanopore Technologies sequencing.

Characterization of Limnospira platensis PCC 9108 R-M and CRISPR-Cas systems

The filamentous cyanobacterium Limnospira platensis, formerly known as Arthrospira platensis or spirulina, is one of the most commercially important species of microalgae. Due to its high nutritional value, pharmacological and industrial applications it is extensively cultivated on a large commercial scale. Despite its widespread use, its precise manipulation is still under development due to the lack of effective genetic protocols. Genetic transformation of Limnospira has been attempted but the methods reported have not been generally reproducible in other laboratories. Knowledge of the transformation defense mechanisms is essential for understanding its physiology and for broadening their applications. With the aim to understand more about the genetic defenses of L. platensis, in this work we have identified the restriction-modification and CRISPR-Cas systems and we have cloned and characterized thirteen methylases. In parallel, we have also characterized the methylome and orphan methyltransferases using genome-wide analysis of DNA methylation patterns and RNA-seq. The identification and characterization of these enzymes will be a valuable resource to know how this strain avoids being genetically manipulated and for further genomics studies.

Tissue-specific RNA methylation prediction from gene expression data using sparse regression models

N6-methyladenosine (m6A) is a highly prevalent and conserved post-transcriptional modification observed in mRNA and long non-coding RNA (lncRNA). Identifying potential m6A sites within RNA sequences is crucial for unraveling the potential influence of the epitranscriptome on biological processes. In this study, we introduce Exp2RM, a novel approach that formulates single-site-based tissue-specific elastic net models for predicting tissue-specific methylation levels utilizing gene expression data. The resulting ensemble model demonstrates robust predictive performance for tissue-specific methylation levels, with an average R-squared value of 0.496 and a median R-squared value of 0.482 across all 22 human tissues. Since methylation distribution varies among tissues, we trained the model to incorporate similar patterns, significantly improves accuracy with the median R-squared value increasing to 0.728. Additonally, functional analysis reveals Exp2RM's ability to capture coefficient genes in relevant biological processes. This study emphasizes the importance of tissue-specific methylation distribution in enhancing prediction accuracy and provides insights into the functional implications of methylation sites.

Characterization of a Novel N4-Methylcytosine Restriction-Modification System in Deinococcus radiodurans

Deinococcus radiodurans is an extremophilic microorganism that possesses a unique DNA damage repair system, conferring a strong resistance to radiation, desiccation, oxidative stress, and chemical damage. Recently, we discovered that D. radiodurans possesses an N4-methylation (m4C) methyltransferase called M.DraR1, which recognizes the 5'-CCGCGG-3' sequence and methylates the second cytosine. Here, we revealed its cognate restriction endonuclease R.DraR1 and recognized that it is the only endonuclease specially for non-4C-methylated 5'-CCGCGG-3' sequence so far. We designated the particular m4C R.DraR1-M.DraR1 as the DraI R-M system. Bioinformatics searches displayed the rarity of the DraI R-M homologous system. Meanwhile, recombination and transformation efficiency experiments demonstrated the important role of the DraI R-M system in response to oxidative stress. In addition, in vitro activity experiments showed that R.DraR1 could exceptionally cleave DNA substrates with a m5C-methlated 5'-CCGCGG-3' sequence instead of its routine activity, suggesting that this particular R-M component possesses a broader substrate choice. Furthermore, an imbalance of the DraI R-M system led to cell death through regulating genes involved in the maintenance of cell survival such as genome stability, transporter, and energy production. Thus, our research revealed a novel m4C R-M system that plays key roles in maintaining cell viability and defending foreign DNA in D. radiodurans.

Adapting Nanopore Sequencing Basecalling Models for Modification Detection via Incremental Learning and Anomaly Detection

We leverage machine learning approaches to adapt nanopore sequencing basecallers for nucleotide modification detection. We first apply the incremental learning technique to improve the basecalling of modification-rich sequences, which are usually of high biological interests. With sequence backbones resolved, we further run anomaly detection on individual nucleotides to determine their modification status. By this means, our pipeline promises the single-molecule, single-nucleotide and sequence context-free detection of modifications. We benchmark the pipeline using control oligos, further apply it in the basecalling of densely-modified yeast tRNAs and E.coli genomic DNAs, the cross-species detection of N6-methyladenosine (m6A) in mammalian mRNAs, and the simultaneous detection of N1-methyladenosine (m1A) and m6A in human mRNAs. Our IL-AD workflow is available at: https://github.com/wangziyuan66/IL-AD.

Lokatt: a hybrid DNA nanopore basecaller with an explicit duration hidden Markov model and a residual LSTM network

BACKGROUND

Basecalling long DNA sequences is a crucial step in nanopore-based DNA sequencing protocols. In recent years, the CTC-RNN model has become the leading basecalling model, supplanting preceding hidden Markov models (HMMs) that relied on pre-segmenting ion current measurements. However, the CTC-RNN model operates independently of prior biological and physical insights.

RESULTS

We present a novel basecaller named Lokatt: explicit duration Markov model and residual-LSTM network. It leverages an explicit duration HMM (EDHMM) designed to model the nanopore sequencing processes. Trained on a newly generated library with methylation-free Ecoli samples and MinION R9.4.1 chemistry, the Lokatt basecaller achieves basecalling performances with a median single read identity score of 0.930, a genome coverage ratio of 99.750%, on par with existing state-of-the-art structure when trained on the same datasets.

CONCLUSION

Our research underlines the potential of incorporating prior knowledge into the basecalling processes, particularly through integrating HMMs and recurrent neural networks. The Lokatt basecaller showcases the efficacy of a hybrid approach, emphasizing its capacity to achieve high-quality basecalling performance while accommodating the nuances of nanopore sequencing. These outcomes pave the way for advanced basecalling methodologies, with potential implications for enhancing the accuracy and efficiency of nanopore-based DNA sequencing protocols.

Comparison of methylation estimates obtained via MinION nanopore sequencing and sanger bisulfite sequencing in the TRPA1 promoter region

BACKGROUND

Bisulfite sequencing has long been considered the gold standard for measuring DNA methylation at single CpG resolution. However, in recent years several new approaches like nanopore sequencing have been developed due to hints for a partial error-proneness of bisulfite sequencing. Since these errors were shown to be sequence-specific, we aimed to verify the methylation data of a particular region of the TRPA1 promoter from our previous studies obtained by bisulfite sequencing.

METHODS

We compared methylation rates determined by direct bisulfite sequencing and nanopore sequencing following Cas9-mediated PCR-free enrichment.

RESULTS

We could show that CpG methylation levels above 20% corroborate well with our previous data. Within the range between 0 and 20% methylation, however, Sanger sequencing data have to be interpreted cautiously, at least in the investigated region of interest (TRPA1 promotor region).

CONCLUSION

Based on the investigation of the TRPA1- region as an example, the present work can help in choosing the right method out of the two current main approaches for methylation analysis for different individual settings regarding many factors like cohort size, costs and prerequisites that should be fulfilled for each method. All in all, both methods have their raison d'être. Furthermore, the present paper contains and illustrates some important basic information and explanation of how guide RNAs should be located for an optimal outcome in Cas9 mediated PCR free target enrichment.

Getting the chronological age out of DNA: using insights of age-dependent DNA methylation for forensic DNA applications

BACKGROUND

DNA analysis for forensic investigations has a long tradition with important developments and optimizations since its first application. Traditionally, short tandem repeats analysis has been the most powerful method for the identification of individuals. However, in addition, epigenetic changes, i.e., DNA methylation, came into focus of forensic DNA research. Chronological age prediction is one promising application to allow for narrowing the pool of possible individuals who caused a trace, as well as to support the identification of unknown bodies and for age verification of living individuals.

OBJECTIVE

This review aims to provide an overview of the current knowledge, possibilities, and (current) limitations about DNA methylation-based chronological age prediction with emphasis on forensic application.

METHODS

The development, implementation and application of age prediction tools requires a deep understanding about the biological background, the analysis methods, the age-dependent DNA methylation markers, as well as the mathematical models for age prediction and their evaluation. Furthermore, additional influences can have an impact. Therefore, the literature was evaluated in respect to these diverse topics.

CONCLUSION

The numerous research efforts in recent years have led to a rapid change in our understanding of the application of DNA methylation for chronological age prediction, which is now on the way to implementation and validation. Knowledge of the various aspects leads to a better understanding and allows a more informed interpretation of DNAm quantification results, as well as the obtained results by the age prediction tools.

Genome-wide DNA methylation landscape of four Chinese populations and epigenetic variation linked to Tibetan high-altitude adaptation

DNA methylation (DNAm) is one of the major epigenetic mechanisms in humans and is important in diverse cellular processes. The variation of DNAm in the human population is related to both genetic and environmental factors. However, the DNAm profiles have not been investigated in the Chinese population of diverse ethnicities. Here, we performed double-strand bisulfite sequencing (DSBS) for 32 Chinese individuals representing four major ethnic groups including Han Chinese, Tibetan, Zhuang, and Mongolian. We identified a total of 604,649 SNPs and quantified DNAm at more than 14 million CpGs in the population. We found global DNAm-based epigenetic structure is different from the genetic structure of the population, and ethnic difference only partially explains the variation of DNAm. Surprisingly, non-ethnic-specific DNAm variations showed stronger correlation with the global genetic divergence than these ethnic-specific DNAm. Differentially methylated regions (DMRs) among these ethnic groups were found around genes in diverse biological processes. Especially, these DMR-genes between Tibetan and non-Tibetans were enriched around high-altitude genes including EPAS1 and EGLN1, suggesting DNAm alteration plays an important role in high-altitude adaptation. Our results provide the first batch of epigenetic maps for Chinese populations and the first evidence of the association of epigenetic changes with Tibetans' high-altitude adaptation.

A Survey of Archaeal Restriction-Modification Systems

When compared with bacteria, relatively little is known about the restriction-modification (RM) systems of archaea, particularly those in taxa outside of the haloarchaea. To improve our understanding of archaeal RM systems, we surveyed REBASE, the restriction enzyme database, to catalog what is known about the genes and activities present in the 519 completely sequenced archaeal genomes currently deposited there. For 49 (9.4%) of these genomes, we also have methylome data from Single-Molecule Real-Time (SMRT) sequencing that reveal the target recognition sites of the active m6A and m4C DNA methyltransferases (MTases). The gene-finding pipeline employed by REBASE is trained primarily on bacterial examples and so will look for similar genes in archaea. Nonetheless, the organizational structure and protein sequence of RM systems from archaea are highly similar to those of bacteria, with both groups acquiring systems from a shared genetic pool through horizontal gene transfer. As in bacteria, we observe numerous examples of "persistent" DNA MTases conserved within archaeal taxa at different levels. We experimentally validated two homologous members of one of the largest "persistent" MTase groups, revealing that methylation of C(m5C)WGG sites may play a key epigenetic role in Crenarchaea. Throughout the archaea, genes encoding m6A, m4C, and m5C DNA MTases, respectively, occur in approximately the ratio 4:2:1.

DARESOME enables concurrent profiling of multiple DNA modifications with restriction enzymes in single cells and cell-free DNA

5-Methylcytosine (5mC) and 5-hydroxymethylcytosine (5hmC) are the most abundant DNA modifications that have important roles in gene regulation. Detailed studies of these different epigenetic marks aimed at understanding their combined effects and dynamic interconversion are, however, hampered by the inability of current methods to simultaneously measure both modifications, particularly in samples with limited quantities. We present DNA analysis by restriction enzyme for simultaneous detection of multiple epigenomic states (DARESOME), an assay based on modification-sensitive restriction digest and sequential tag ligation that can concurrently perform quantitative profiling of unmodified cytosine, 5mC, and 5hmC in CCGG sites genome-wide. DARESOME reveals the opposing roles of 5mC and 5hmC in gene expression regulation as well as their interconversion during aging in mouse brain. Implementation of DARESOME in single cells demonstrates pronounced 5hmC strand bias that reflects the semiconservative replication of DNA. Last, we showed that DARESOME enables integrative genomic, 5mC, and 5hmC profiling of cell-free DNA that uncovered multiomics cancer signatures in liquid biopsy.

Regulation of plant epigenetic memory in response to cold and heat stress: towards climate resilient agriculture

Plants have evolved to adapt and grow in hot and cold climatic conditions. Some also adapt to daily and seasonal temperature changes. Epigenetic modifications play an important role in regulating plant tolerance under such conditions. DNA methylation and post-translational modifications of histone proteins influence gene expression during plant developmental stages and under stress conditions, including cold and heat stress. While short-term modifications are common, some modifications may persist and result in stress memory that can be inherited by subsequent generations. Understanding the mechanisms of epigenomes responding to stress and the factors that trigger stress memory is crucial for developing climate-resilient agriculture, but such an integrated view is currently limited. This review focuses on the plant epigenetic stress memory during cold and heat stress. It also discusses the potential of machine learning to modify stress memory through epigenetics to develop climate-resilient crops.

Inhibition of DNA and RNA methylation disturbs root development of moso bamboo

DNA methylation (5mC) and N6-methyladenosine (m6A) are two important epigenetics regulators, which have a profound impact on plant growth development. Phyllostachys edulis (P. edulis) is one of the fastest spreading plants due to its well-developed root system. However, the association between 5mC and m6A has seldom been reported in P. edulis. In particular, the connection between m6A and several post-transcriptional regulators remains uncharacterized in P. edulis. Here, our morphological and electron microscope observations showed the phenotype of increased lateral root under RNA methylation inhibitor (DZnepA) and DNA methylation inhibitor (5-azaC) treatment. RNA epitranscriptome based on Nanopore direct RNA sequencing revealed that DZnepA treatment exhibits significantly decreased m6A level in the 3'-untranslated region (3'-UTR), which was accompanied by increased gene expression, full-length ratio, higher proximal poly(A) site usage and shorter poly(A) tail length. DNA methylation levels of CG and CHG were reduced in both coding sequencing and transposable element upon 5-azaC treatment. Cell wall synthesis was impaired under methylation inhibition. In particular, differentially expressed genes showed a high percentage of overlap between DZnepA and 5-azaC treatment, which suggested a potential correlation between two methylations. This study provides preliminary information for a better understanding of the link between m6A and 5mC in root development of moso bamboo.

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.

Long-Read Sequencing Reveals Extensive DNA Methylations in Human Gut Phagenome Contributed by Prevalently Phage-Encoded Methyltransferases

DNA methylation plays a crucial role in the survival of bacteriophages (phages), yet the understanding of their genome methylation remains limited. In this study, DNA methylation patterns are analyzed in 8848 metagenome-assembled high-quality phages from 104 fecal samples using single-molecule real-time sequencing. The results demonstrate that 97.60% of gut phages exhibit methylation, with certain factors correlating with methylation densities. Phages with higher methylation densities appear to have potential viability advantages. Strikingly, more than one-third of the phages possess their own DNA methyltransferases (MTases). Increased MTase copies are associated with higher genome methylation densities, specific methylation motifs, and elevated prevalence of certain phage groups. Notably, the majority of these MTases share close homology with those encoded by gut bacteria, suggesting their exchange during phage-bacterium interactions. Furthermore, these MTases can be employed to accurately predict phage-host relationships. Overall, the findings indicate the widespread utilization of DNA methylation by gut DNA phages as an evasion mechanism against host defense systems, with a substantial contribution from phage-encoded MTases.

Mapping epigenetic modifications by sequencing technologies

The "epigenetics" concept was first described in 1942. Thus far, chemical modifications on histones, DNA, and RNA have emerged as three important building blocks of epigenetic modifications. Many epigenetic modifications have been intensively studied and found to be involved in most essential biological processes as well as human diseases, including cancer. Precisely and quantitatively mapping over 100 [1], 17 [2], and 160 [3] different known types of epigenetic modifications in histone, DNA, and RNA is the key to understanding the role of epigenetic modifications in gene regulation in diverse biological processes. With the rapid development of sequencing technologies, scientists are able to detect specific epigenetic modifications with various quantitative, high-resolution, whole-genome/transcriptome approaches. Here, we summarize recent advances in epigenetic modification sequencing technologies, focusing on major histone, DNA, and RNA modifications in mammalian cells.

Long-Read DNA Sequencing: Recent Advances and Remaining Challenges

DNA sequencing has revolutionized medicine over recent decades. However, analysis of large structural variation and repetitive DNA, a hallmark of human genomes, has been limited by short-read technology, with read lengths of 100-300 bp. Long-read sequencing (LRS) permits routine sequencing of human DNA fragments tens to hundreds of kilobase pairs in size, using both real-time sequencing by synthesis and nanopore-based direct electronic sequencing. LRS permits analysis of large structural variation and haplotypic phasing in human genomes and has enabled the discovery and characterization of rare pathogenic structural variants and repeat expansions. It has also recently enabled the assembly of a complete, gapless human genome that includes previously intractable regions, such as highly repetitive centromeres and homologous acrocentric short arms. With the addition of protocols for targeted enrichment, direct epigenetic DNA modification detection, and long-range chromatin profiling, LRS promises to launch a new era of understanding of genetic diversity and pathogenic mutations in human populations.

Shaping eukaryotic epigenetic systems by horizontal gene transfer

DNA methylation constitutes one of the pillars of epigenetics, relying on covalent bonds for addition and/or removal of chemically distinct marks within the major groove of the double helix. DNA methyltransferases, enzymes which introduce methyl marks, initially evolved in prokaryotes as components of restriction-modification systems protecting host genomes from bacteriophages and other invading foreign DNA. In early eukaryotic evolution, DNA methyltransferases were horizontally transferred from bacteria into eukaryotes several times and independently co-opted into epigenetic regulatory systems, primarily via establishing connections with the chromatin environment. While C5-methylcytosine is the cornerstone of plant and animal epigenetics and has been investigated in much detail, the epigenetic role of other methylated bases is less clear. The recent addition of N4-methylcytosine of bacterial origin as a metazoan DNA modification highlights the prerequisites for foreign gene co-option into the host regulatory networks, and challenges the existing paradigms concerning the origin and evolution of eukaryotic regulatory systems.

DNA hypomethylation silences anti-tumor immune genes in early prostate cancer and CTCs

Cancer is characterized by hypomethylation-associated silencing of large chromatin domains, whose contribution to tumorigenesis is uncertain. Through high-resolution genome-wide single-cell DNA methylation sequencing, we identify 40 core domains that are uniformly hypomethylated from the earliest detectable stages of prostate malignancy through metastatic circulating tumor cells (CTCs). Nested among these repressive domains are smaller loci with preserved methylation that escape silencing and are enriched for cell proliferation genes. Transcriptionally silenced genes within the core hypomethylated domains are enriched for immune-related genes; prominent among these is a single gene cluster harboring all five CD1 genes that present lipid antigens to NKT cells and four IFI16-related interferon-inducible genes implicated in innate immunity. The re-expression of CD1 or IFI16 murine orthologs in immuno-competent mice abrogates tumorigenesis, accompanied by the activation of anti-tumor immunity. Thus, early epigenetic changes may shape tumorigenesis, targeting co-located genes within defined chromosomal loci. Hypomethylation domains are detectable in blood specimens enriched for CTCs.

Y chromosome sequence and epigenomic reconstruction across human populations

Recent advances in long-read sequencing technologies have allowed the generation and curation of more complete genome assemblies, enabling the analysis of traditionally neglected chromosomes, such as the human Y chromosome (chrY). Native DNA was sequenced on a MinION Oxford Nanopore Technologies sequencing device to generate genome assemblies for seven major chrY human haplogroups. We analyzed and compared the chrY enrichment of sequencing data obtained using two different selective sequencing approaches: adaptive sampling and flow cytometry chromosome sorting. We show that adaptive sampling can produce data to create assemblies comparable to chromosome sorting while being a less expensive and time-consuming technique. We also assessed haplogroup-specific structural variants, which would be otherwise difficult to study using short-read sequencing data only. Finally, we took advantage of this technology to detect and profile epigenetic modifications among the considered haplogroups. Altogether, we provide a framework to study complex genomic regions with a simple, fast, and affordable methodology that could be applied to larger population genomics datasets.

Navigating the pitfalls of mapping DNA and RNA modifications

Chemical modifications to nucleic acids occur across the kingdoms of life and carry important regulatory information. Reliable high-resolution mapping of these modifications is the foundation of functional and mechanistic studies, and recent methodological advances based on next-generation sequencing and long-read sequencing platforms are critical to achieving this aim. However, mapping technologies may have limitations that sometimes lead to inconsistent results. Some of these limitations are technical in nature and specific to certain types of technology. Here, however, we focus on common (yet not always widely recognized) pitfalls that are shared among frequently used mapping technologies and discuss strategies to help technology developers and users mitigate their effects. Although the emphasis is primarily on DNA modifications, RNA modifications are also discussed.

Targeted adaptive long-read sequencing for discovery of complex phased variants in inherited retinal disease patients

Inherited retinal degenerations (IRDs) are a heterogeneous group of predominantly monogenic disorders with over 300 causative genes identified. Short-read exome sequencing is commonly used to genotypically diagnose patients with clinical features of IRDs, however, in up to 30% of patients with autosomal recessive IRDs, one or no disease-causing variants are identified. Furthermore, chromosomal maps cannot be reconstructed for allelic variant discovery with short-reads. Long-read genome sequencing can provide complete coverage of disease loci and a targeted approach can focus sequencing bandwidth to a genomic region of interest to provide increased depth and haplotype reconstruction to uncover cases of missing heritability. We demonstrate that targeted adaptive long-read sequencing on the Oxford Nanopore Technologies (ONT) platform of the USH2A gene from three probands in a family with the most common cause of the syndromic IRD, Usher Syndrome, resulted in greater than 12-fold target gene sequencing enrichment on average. This focused depth of sequencing allowed for haplotype reconstruction and phased variant identification. We further show that variants obtained from the haplotype-aware genotyping pipeline can be heuristically ranked to focus on potential pathogenic candidates without a priori knowledge of the disease-causing variants. Moreover, consideration of the variants unique to targeted long-read sequencing that are not covered by short-read technology demonstrated higher precision and F1 scores for variant discovery by long-read sequencing. This work establishes that targeted adaptive long-read sequencing can generate targeted, chromosome-phased data sets for identification of coding and non-coding disease-causing alleles in IRDs and can be applicable to other Mendelian diseases.

6mA DNA Methylation on Genes in Plants Is Associated with Gene Complexity, Expression and Duplication

N6-methyladenine (6mA) DNA methylation has emerged as an important epigenetic modification in eukaryotes. Nevertheless, the evolution of the 6mA methylation of homologous genes after species and after gene duplications remains unclear in plants. To understand the evolution of 6mA methylation, we detected the genome-wide 6mA methylation patterns of four lotus plants (Nelumbo nucifera) from different geographic origins by nanopore sequencing and compared them to patterns in Arabidopsis and rice. Within lotus, the genomic distributions of 6mA sites are different from the widely studied 5mC methylation sites. Consistently, in lotus, Arabidopsis and rice, 6mA sites are enriched around transcriptional start sites, positively correlated with gene expression levels, and preferentially retained in highly and broadly expressed orthologs with longer gene lengths and more exons. Among different duplicate genes, 6mA methylation is significantly more enriched and conserved in whole-genome duplicates than in local duplicates. Overall, our study reveals the convergent patterns of 6mA methylation evolution based on both lineage and duplicate gene divergence, which underpin their potential role in gene regulatory evolution in plants.

Adaptation of Commensal Escherichia coli in Tomato Fruits: Motility, Stress, Virulence

Food contamination can be a serious concern for public health because it can be related to the severe spreading of pathogens. This is a main issue, especially in the case of fresh fruits and vegetables; indeed, they have often been associated with gastrointestinal outbreak events, due to contamination with pathogenic bacteria. However, little is known about the physiological adaptation and bacterial response to stresses encountered in the host plant. Thus, this work aimed to investigate the adaptation of a commensal E. coli strain while growing in tomato pericarp. Pre-adapted and non-adapted cells were compared and used to contaminate tomatoes, demonstrating that pre-adaptation boosted cell proliferation. DNA extracted from pre-adapted and non-adapted cells was sequenced, and their methylation profiles were compared. Hence, genes involved in cell adhesion and resistance against toxic compounds were identified as genes involved in adaptation, and their expression was compared in these two experimental conditions. Finally, pre-adapted and non-adapted E. coli were tested for their ability to resist the presence of toxic compounds, demonstrating that adaptation exerted a protective effect. In conclusion, this work provides new information about the physiological adaptation of bacteria colonizing the tomato fruit pericarp.

Assessing Readability of an 8-Letter Expanded Deoxyribonucleic Acid Alphabet with Nanopores

Chemists have now synthesized new kinds of DNA that add nucleotides to the four standard nucleotides (guanine, adenine, cytosine, and thymine) found in standard Terran DNA. Such "artificially expanded genetic information systems" are today used in molecular diagnostics; to support directed evolution to create medically useful receptors, ligands, and catalysts; and to explore issues related to the early evolution of life. Further applications are limited by the inability to directly sequence DNA containing nonstandard nucleotides. Nanopore sequencing is well-suited for this purpose, as it does not require enzymatic synthesis, amplification, or nucleotide modification. Here, we take the first steps to realize nanopore sequencing of an 8-letter "hachimoji" expanded DNA alphabet by assessing its nanopore signal range using the MspA (Mycobacterium smegmatis porin A) nanopore. We find that hachimoji DNA exhibits a broader signal range in nanopore sequencing than standard DNA alone and that hachimoji single-base substitutions are distinguishable with high confidence. Because nanopore sequencing relies on a molecular motor to control the motion of DNA, we then assessed the compatibility of the Hel308 motor enzyme with nonstandard nucleotides by tracking the translocation of single Hel308 molecules along hachimoji DNA, monitoring the enzyme kinetics and premature enzyme dissociation from the DNA. We find that Hel308 is compatible with hachimoji DNA but dissociates more frequently when walking over C-glycoside nucleosides, compared to N-glycosides. C-glycocide nucleosides passing a particular site within Hel308 induce a higher likelihood of dissociation. This highlights the need to optimize nanopore sequencing motors to handle different glycosidic bonds. It may also inform designs of future alternative DNA systems that can be sequenced with existing motors and pores.

BIND&MODIFY: a long-range method for single-molecule mapping of chromatin modifications in eukaryotes

Epigenetic modifications of histones are associated with development and pathogenesis of disease. Existing approaches cannot provide insights into long-range interactions and represent the average chromatin state. Here we describe BIND&MODIFY, a method using long-read sequencing for profiling histone modifications and transcription factors on individual DNA fibers. We use recombinant fused protein A-M.EcoGII to tether methyltransferase M.EcoGII to protein binding sites to label neighboring regions by methylation. Aggregated BIND&MODIFY signal matches bulk ChIP-seq and CUT&TAG. BIND&MODIFY can simultaneously measure histone modification status, transcription factor binding, and CpG 5mC methylation at single-molecule resolution and also quantifies correlation between local and distal elements.

Hybridization-based CpG methylation level detection using methyl-CpG-binding domain-fused luciferase

Hypermethylation of tumor-suppressor genes and global hypomethylation, which is related to methylation level at the retroelement, have been recognized as features of the cancer genome. In this study, we developed a hybridization-based CpG methylation level detection method using methyl-CpG-binding domain-fused firefly luciferase (MBD-Fluc). In this method, methylated probe oligonucleotides were used to capture target oligonucleotides. Fully methylated and hemimethylated double-stranded DNA (dsDNA) was formed by hybridization of the methylated captured oligonucleotides with methylated or unmethylated target oligonucleotides, respectively. MBD-Fluc specifically binds to fully methylated dsDNA but not to hemimethylated dsDNA; therefore, methylated target oligonucleotides can be detected by measuring the luciferase activity of the bound MBD-Fluc. Using the corresponding methylated probe oligonucleotides, the CpG methylation levels of SEPT9, BRCA1, and long interspersed nuclear element-1 (LINE-1) oligonucleotides were quantified. Moreover, we demonstrated that the emission detection signal was not affected by the methylation state of the overhang region of the target oligonucleotide, which was not hybridized to the probe oligonucleotide, indicating that methylated CpG of the target region could be accurately detected. Unmethylated-CpG-binding domain-fused luciferases and 5-hydroxymethyl-CpG-binding domain-fused luciferases have been constructed, suggesting that other modified bases can be detected by the same platform.

Single-cell and long-read sequencing to enhance modelling of splicing and cell-fate determination

Single-cell sequencing technologies have revolutionised the life sciences and biomedical research. Single-cell sequencing provides high-resolution data on cell heterogeneity, allowing high-fidelity cell type identification, and lineage tracking. Computational algorithms and mathematical models have been developed to make sense of the data, compensate for errors and simulate the biological processes, which has led to breakthroughs in our understanding of cell differentiation, cell-fate determination and tissue cell composition. The development of long-read (a.k.a. third-generation) sequencing technologies has produced powerful tools for investigating alternative splicing, isoform expression (at the RNA level), genome assembly and the detection of complex structural variants (at the DNA level). In this review, we provide an overview of the recent advancements in single-cell and long-read sequencing technologies, with a particular focus on the computational algorithms that help in correcting, analysing, and interpreting the resulting data. Additionally, we review some mathematical models that use single-cell and long-read sequencing data to study cell-fate determination and alternative splicing, respectively. Moreover, we highlight the emerging opportunities in modelling cell-fate determination that result from the combination of single-cell and long-read sequencing technologies.

Detection of nucleotide modifications in bacteria and bacteriophages: Strengths and limitations of current technologies and software

RNA and DNA modifications occur in eukaryotes and prokaryotes, as well as in their viruses, and serve a wide range of functions, from gene regulation to nucleic acid protection. Although the first nucleotide modification was discovered almost 100 years ago, new and unusual modifications are still being described. Nucleotide modifications have also received more attention lately because of their increased significance, but also because new sequencing approaches have eased their detection. Chiefly, third generation sequencing platforms PacBio and Nanopore offer direct detection of modified bases by measuring deviations of the signals. These unusual modifications are especially prevalent in bacteriophage genomes, the viruses of bacteria, where they mostly appear to protect DNA against degradation from host nucleases. In this Opinion article, we highlight and discuss current approaches to detect nucleotide modifications, including hardwares and softwares, and look onward to future applications, especially for studying unusual, rare, or complex genome modifications in bacteriophages. The ability to distinguish between several types of nucleotide modifications may even shed new light on metagenomic studies.

An Introduction to Nanopore Sequencing: Past, Present, and Future Considerations

There has been significant progress made in the field of nanopore biosensor development and sequencing applications, which address previous limitations that restricted widespread nanopore use. These innovations, paired with the large-scale commercialization of biological nanopore sequencing by Oxford Nanopore Technologies, are making the platforms a mainstay in contemporary research laboratories. Equipped with the ability to provide long- and short read sequencing information, with quick turn-around times and simple sample preparation, nanopore sequencers are rapidly improving our understanding of unsolved genetic, transcriptomic, and epigenetic problems. However, there remain some key obstacles that have yet to be improved. In this review, we provide a general introduction to nanopore sequencing principles, discussing biological and solid-state nanopore developments, obstacles to single-base detection, and library preparation considerations. We present examples of important clinical applications to give perspective on the potential future of nanopore sequencing in the field of molecular diagnostics.