Detecting DNA cytosine methylation using nanopore sequencing

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.

Nanopore Single-molecule Analysis of Biomarkers: Providing Possible Clues to Disease Diagnosis

Biomarker detection has attracted increasing interest in recent years due to the minimally or non-invasive sampling process. Single entity analysis of biomarkers is expected to provide real-time and accurate biological information for early disease diagnosis and prognosis, which is critical to the effective disease treatment and is also important in personalized medicine. As an innovative single entity analysis method, nanopore sensing is a pioneering single-molecule detection technique that is widely used in analytical bioanalytical fields. In this review, we overview the recent progress of nanopore biomarker detection as new approaches to disease diagnosis. In highlighted studies, nanopore was focusing on detecting biomarkers of different categories of communicable and noncommunicable diseases, such as pandemic Covid-19, AIDS, cancers, neurologic diseases, etc. Various sensitive and selective nanopore detecting strategies for different types of biomarkers are summarized. In addition, the challenges, opportunities, and direction for future development of nanopore-based biomarker sensors are also discussed.

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.

pycoMeth: a toolbox for differential methylation testing from Nanopore methylation calls

We present pycoMeth, a toolbox to store, manage and analyze DNA methylation calls from long-read sequencing data obtained using the Oxford Nanopore Technologies sequencing platform. Building on a novel, rapid-access, read-level and reference-anchored methylation storage format MetH5, we propose efficient algorithms for haplotype aware, multi-sample consensus segmentation and differential methylation testing. We show that MetH5 is more efficient than existing solutions for storing Oxford Nanopore Technologies methylation calls, and carry out benchmarking for pycoMeth segmentation and differential methylation testing, demonstrating increased performance and sensitivity compared to existing solutions designed for short-read methylation data.

Long-read sequencing of diagnosis and post-therapy medulloblastoma reveals complex rearrangement patterns and epigenetic signatures

Cancer genomes harbor a broad spectrum of structural variants (SVs) driving tumorigenesis, a relevant subset of which escape discovery using short-read sequencing. We employed Oxford Nanopore Technologies (ONT) long-read sequencing in a paired diagnostic and post-therapy medulloblastoma to unravel the haplotype-resolved somatic genetic and epigenetic landscape. We assembled complex rearrangements, including a 1.55-Mbp chromothripsis event, and we uncover a complex SV pattern termed templated insertion (TI) thread, characterized by short (mostly <1 kb) insertions showing prevalent self-concatenation into highly amplified structures of up to 50 kbp in size. TI threads occur in 3% of cancers, with a prevalence up to 74% in liposarcoma, and frequent colocalization with chromothripsis. We also perform long-read-based methylome profiling and discover allele-specific methylation (ASM) effects, complex rearrangements exhibiting differential methylation, and differential promoter methylation in cancer-driver genes. Our study shows the advantage of long-read sequencing in the discovery and characterization of complex somatic rearrangements.

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.

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

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

Flexible and efficient handling of nanopore sequencing signal data with slow5tools

Nanopore sequencing is being rapidly adopted in genomics. We recently developed SLOW5, a new file format with advantages for storage and analysis of raw signal data from nanopore experiments. Here we introduce slow5tools, an intuitive toolkit for handling nanopore data in SLOW5 format. Slow5tools enables lossless data conversion and a range of tools for interacting with SLOW5 files. Slow5tools uses multi-threading, multi-processing, and other engineering strategies to achieve fast data conversion and manipulation, including live FAST5-to-SLOW5 conversion during sequencing. We provide examples and benchmarking experiments to illustrate slow5tools usage, and describe the engineering principles underpinning its performance.

Systematic comparison of tools used for m6A mapping from nanopore direct RNA sequencing

N6-methyladenosine (m6A) has been increasingly recognized as a new and important regulator of gene expression. To date, transcriptome-wide m6A detection primarily relies on well-established methods using next-generation sequencing (NGS) platform. However, direct RNA sequencing (DRS) using the Oxford Nanopore Technologies (ONT) platform has recently emerged as a promising alternative method to study m6A. While multiple computational tools are being developed to facilitate the direct detection of nucleotide modifications, little is known about the capabilities and limitations of these tools. Here, we systematically compare ten tools used for mapping m6A from ONT DRS data. We find that most tools present a trade-off between precision and recall, and integrating results from multiple tools greatly improve performance. Using a negative control could improve precision by subtracting certain intrinsic bias. We also observed variation in detection capabilities and quantitative information among motifs, and identified sequencing depth and m6A stoichiometry as potential factors affecting performance. Our study provides insight into the computational tools currently used for mapping m6A based on ONT DRS data and highlights the potential for further improving these tools, which may serve as the basis for future research.

Profiling the epigenetic landscape of the antigen receptor repertoire: the missing epi-immunogenomics data

High-resolution sequencing methods that capture the epigenetic landscape within the T cell receptor (TCR) gene loci are pivotal for a fundamental understanding of the epigenetic regulatory mechanisms of the TCR repertoire. In our opinion, filling the gaps in our understanding of the epigenetic mechanisms regulating the TCR repertoire will benefit the development of strategies that can modulate the TCR repertoire composition by leveraging the dynamic nature of epigenetic modifications.

Dynamic antagonism between key repressive pathways maintains the placental epigenome

DNA and Histone 3 Lysine 27 methylation typically function as repressive modifications and operate within distinct genomic compartments. In mammals, the majority of the genome is kept in a DNA methylated state, whereas the Polycomb repressive complexes regulate the unmethylated CpG-rich promoters of developmental genes. In contrast to this general framework, the extra-embryonic lineages display non-canonical, globally intermediate DNA methylation levels, including disruption of local Polycomb domains. Here, to better understand this unusual landscape's molecular properties, we genetically and chemically perturbed major epigenetic pathways in mouse trophoblast stem cells. We find that the extra-embryonic epigenome reflects ongoing and dynamic de novo methyltransferase recruitment, which is continuously antagonized by Polycomb to maintain intermediate, locally disordered methylation. Despite its disorganized molecular appearance, our data point to a highly controlled equilibrium between counteracting repressors within extra-embryonic cells, one that can seemingly persist indefinitely without bistable features typically seen for embryonic forms of epigenetic regulation.

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.

A gap-free genome assembly of Chlamydomonas reinhardtii and detection of translocations induced by CRISPR-mediated mutagenesis

Genomic assemblies of the unicellular green alga Chlamydomonas reinhardtii have provided important resources for researchers. However, assembly errors, large gaps, and unplaced scaffolds as well as strain-specific variants currently impede many types of analysis. By combining PacBio HiFi and Oxford Nanopore long-read technologies, we generated a de novo genome assembly for strain CC-5816, derived from crosses of strains CC-125 and CC-124. Multiple methods of evaluating genome completeness and base-pair error rate suggest that the final telomere-to-telomere assembly is highly accurate. The CC-5816 assembly enabled previously difficult analyses that include characterization of the 17 centromeres, rDNA arrays on three chromosomes, and 56 insertions of organellar DNA into the nuclear genome. Using Nanopore sequencing, we identified sites of cytosine (CpG) methylation, which are enriched at centromeres. We analyzed CRISPR-Cas9 insertional mutants in the PF23 gene. Two of the three alleles produced progeny that displayed patterns of meiotic inviability that suggested the presence of a chromosomal aberration. Mapping Nanopore reads from pf23-2 and pf23-3 onto the CC-5816 genome showed that these two strains each carry a translocation that was initiated at the PF23 gene locus on chromosome 11 and joined with chromosomes 5 or 3, respectively. The translocations were verified by demonstrating linkage between loci on the two translocated chromosomes in meiotic progeny. The three pf23 alleles display the expected short-cilia phenotype, and immunoblotting showed that pf23-2 lacks the PF23 protein. Our CC-5816 genome assembly will undoubtedly provide an important tool for the Chlamydomonas research community.

Application of third-generation sequencing to herbal genomics

There is a long history of traditional medicine use. However, little genetic information is available for the plants used in traditional medicine, which limits the exploitation of these natural resources. Third-generation sequencing (TGS) techniques have made it possible to gather invaluable genetic information and develop herbal genomics. In this review, we introduce two main TGS techniques, PacBio SMRT technology and Oxford Nanopore technology, and compare the two techniques against Illumina, the predominant next-generation sequencing technique. In addition, we summarize the nuclear and organelle genome assemblies of commonly used medicinal plants, choose several examples from genomics, transcriptomics, and molecular identification studies to dissect the specific processes and summarize the advantages and disadvantages of the two TGS techniques when applied to medicinal organisms. Finally, we describe how we expect that TGS techniques will be widely utilized to assemble telomere-to-telomere (T2T) genomes and in epigenomics research involving medicinal plants.

Benchmarking of Nanopore R10.4 and R9.4.1 flow cells in single-cell whole-genome amplification and whole-genome shotgun sequencing

Third-generation sequencing can be used in human cancer genomics and epigenomic research. Oxford Nanopore Technologies (ONT) recently released R10.4 flow cell, which claimed an improved read accuracy compared to R9.4.1 flow cell. To evaluate the benefits and defects of R10.4 flow cell for cancer cell profiling on MinION devices, we used the human non-small-cell lung-carcinoma cell line HCC78 to construct libraries for both single-cell whole-genome amplification (scWGA) and whole-genome shotgun sequencing. The R10.4 and R9.4.1 reads were benchmarked in terms of read accuracy, variant detection, modification calling, genome recovery rate and compared with the next generation sequencing (NGS) reads. The results highlighted that the R10.4 outperforms R9.4.1 reads, achieving a higher modal read accuracy of over 99.1%, superior variation detection, lower false-discovery rate (FDR) in methylation calling, and comparable genome recovery rate. To achieve high yields scWGA sequencing in the ONT platform as NGS, we recommended multiple displacement amplification with a modified T7 endonuclease Ⅰ cutting procedure as a promising method. In addition, we provided a possible solution to filter the likely false positive sites among the whole genome region with R10.4 by using scWGA sequencing result as a negative control. Our study is the first benchmark of whole genome single-cell sequencing using ONT R10.4 and R9.4.1 MinION flow cells by clarifying the capacity of genomic and epigenomic profiling within a single flow cell. A promising method for scWGA sequencing together with the methylation calling results can benefit researchers who work on cancer cell genomic and epigenomic profiling using third-generation sequencing.

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.

Detection of DNA Modification Using Nanopore Sequencers

DNA modification is a crucial factor of epigenetic modification and has vital functions for gene regulation and phenotype control. A profound understanding of DNA modification requires precise mapping of the modified bases on genomic DNA. In addition to methods such as bisulfite sequencing and single-molecule real-time (SMRT) sequencing of PacBio sequencers, nanopore sequencers can be also utilized for the detection of DNA modification. Here, I will briefly review the three methods for the detection of DNA modification with nanopore sequencers and introduce a protocol using MinION and Megalodon.

Advances in the oral microbiota and rapid detection of oral infectious diseases

Several studies have shown that the dysregulation of the oral microbiota plays a crucial role in human health conditions, such as dental caries, periodontal disease, oral cancer, other oral infectious diseases, cardiovascular diseases, diabetes, bacteremia, and low birth weight. The use of traditional detection methods in conjunction with rapidly advancing molecular techniques in the diagnosis of harmful oral microorganisms has expanded our understanding of the diversity, location, and function of the microbiota associated with health and disease. This review aimed to highlight the latest knowledge in this field, including microbial colonization; the most modern detection methods; and interactions in disease progression. The next decade may achieve the rapid diagnosis and precise treatment of harmful oral microorganisms.

Factors and Methods for the Detection of Gene Expression Regulation

Gene-expression regulation involves multiple processes and a range of regulatory factors. In this review, we describe the key factors that regulate gene expression, including transcription factors (TFs), chromatin accessibility, histone modifications, DNA methylation, and RNA modifications. In addition, we also describe methods that can be used to detect these regulatory factors.

Studying Whole-Genome Duplication Using Experimental Evolution of Chlamydomonas

In this chapter, we present the use of Chlamydomonas reinhardtii in experiments designed to study the evolutionary impacts of whole genome duplication. We shortly introduce the algal species and depict why it is an excellent model for experimental evolution. Subsequently, we discuss the most relevant steps and methods in the design of a ploidy-related Chlamydomonas experiment. These steps include strain selection, ploidy determination, different methods of making diplo- and polyploid Chlamydomonas cells, replication, culturing conditions, preservation, and the ways to quantify phenotypic and genotypic change.

Portable nanopore-sequencing technology: Trends in development and applications

Sequencing technology is the most commonly used technology in molecular biology research and an essential pillar for the development and applications of molecular biology. Since 1977, when the first generation of sequencing technology opened the door to interpreting the genetic code, sequencing technology has been developing for three generations. It has applications in all aspects of life and scientific research, such as disease diagnosis, drug target discovery, pathological research, species protection, and SARS-CoV-2 detection. However, the first- and second-generation sequencing technology relied on fluorescence detection systems and DNA polymerization enzyme systems, which increased the cost of sequencing technology and limited its scope of applications. The third-generation sequencing technology performs PCR-free and single-molecule sequencing, but it still depends on the fluorescence detection device. To break through these limitations, researchers have made arduous efforts to develop a new advanced portable sequencing technology represented by nanopore sequencing. Nanopore technology has the advantages of small size and convenient portability, independent of biochemical reagents, and direct reading using physical methods. This paper reviews the research and development process of nanopore sequencing technology (NST) from the laboratory to commercially viable tools; discusses the main types of nanopore sequencing technologies and their various applications in solving a wide range of real-world problems. In addition, the paper collates the analysis tools necessary for performing different processing tasks in nanopore sequencing. Finally, we highlight the challenges of NST and its future research and application directions.

Robust methylation-based classification of brain tumours using nanopore sequencing

BACKGROUND

DNA methylation-based classification of cancer provides a comprehensive molecular approach to diagnose tumours. In fact, DNA methylation profiling of human brain tumours already profoundly impacts clinical neuro-oncology. However, current implementation using hybridisation microarrays is time consuming and costly. We recently reported on shallow nanopore whole-genome sequencing for rapid and cost-effective generation of genome-wide 5-methylcytosine profiles as input to supervised classification. Here, we demonstrate that this approach allows us to discriminate a wide spectrum of primary brain tumours.

RESULTS

Using public reference data of 82 distinct tumour entities, we performed nanopore genome sequencing on 382 tissue samples covering 46 brain tumour (sub)types. Using bootstrap sampling in a cohort of 55 cases, we found that a minimum set of 1000 random CpG features is sufficient for high-confidence classification by ad hoc random forests. We implemented score recalibration as a confidence measure for interpretation in a clinical context and empirically determined a platform-specific threshold in a randomly sampled discovery cohort (N = 185). Applying this cut-off to an independent validation series (n = 184) yielded 148 classifiable cases (sensitivity 80.4%) and demonstrated 100% specificity. Cross-lab validation demonstrated robustness with concordant results across four laboratories in 10/11 (90.9%) cases. In a prospective benchmarking (N = 15), the median time to results was 21.1 h.

CONCLUSIONS

In conclusion, nanopore sequencing allows robust and rapid methylation-based classification across the full spectrum of brain tumours. Platform-specific confidence scores facilitate clinical implementation for which prospective evaluation is warranted and ongoing.

Using Nanocompore to Identify RNA Modifications from Direct RNA Nanopore Sequencing Data

RNA modifications can alter the behavior of RNA molecules depending on where they are located on the strands. Traditionally, RNA modifications have been detected and characterized by biophysical assays, mass spectrometry, or specific next-generation sequencing techniques, but are limited to specific modifications or are low throughput. Nanopore is a platform capable of sequencing RNA strands directly, which permits transcriptome-wide detection of RNA modifications. RNA modifications alter the nanopore raw signal relative to the canonical form of the nucleotide, and several software tools detect these signal alterations. One such tool is Nanocompore, which compares the ionic current features between two different experimental conditions (i.e., with and without RNA modifications) to detect RNA modifications. Nanocompore is not limited to a single type of RNA modification, has a high specificity for detecting RNA modifications, and does not require model training. To use Nanocompore, the following steps are needed: (i) the data must be basecalled and aligned to the reference transcriptome, then the raw ionic current signals are aligned to the sequences and transformed into a Nanocompore-compatible format; (ii) finally, the statistical testing is conducted on the transformed data and produces a table of p-value predictions for the positions of the RNA modifications. These steps can be executed with several different methods, and thus we have also included two alternative protocols for running Nanocompore. Once the positions of RNA modifications are determined by Nanocompore, users can investigate their function in various metabolic pathways. © 2023 The Authors. Current Protocols published by Wiley Periodicals LLC. Basic Protocol: RNA modification detection by Nanocompore Alternate Protocol 1: RNA modification detection by Nanocompore with f5c Alternate Protocol 2: RNA modification detection by Nanocompore using Nextflow.

Genome-wide evaluation of the effect of short tandem repeat variation on local DNA methylation

Short tandem repeats (STRs) contribute significantly to genetic diversity in humans, including disease-causing variation. Although the effect of STR variation on gene expression has been extensively assessed, their impact on epigenetics has been poorly studied and limited to specific genomic regions. Here, we investigated the hypothesis that some STRs act as independent regulators of local DNA methylation in the human genome and modify risk of common human traits. To address these questions, we first analyzed two independent data sets comprising PCR-free whole-genome sequencing (WGS) and genome-wide DNA methylation levels derived from whole-blood samples in 245 (discovery cohort) and 484 individuals (replication cohort). Using genotypes for 131,635 polymorphic STRs derived from WGS using HipSTR, we identified 11,870 STRs that associated with DNA methylation levels (mSTRs) of 11,774 CpGs (Bonferroni P < 0.001) in our discovery cohort, with 90% successfully replicating in our second cohort. Subsequently, through fine-mapping using CAVIAR we defined 585 of these mSTRs as the likely causal variants underlying the observed associations (fm-mSTRs) and linked a fraction of these to previously reported genome-wide association study signals, providing insights into the mechanisms underlying complex human traits. Furthermore, by integrating gene expression data, we observed that 12.5% of the tested fm-mSTRs also modulate expression levels of nearby genes, reinforcing their regulatory potential. Overall, our findings expand the catalog of functional sequence variants that affect genome regulation, highlighting the importance of incorporating STRs in future genetic association analysis and epigenetics data for the interpretation of trait-associated variants.

Computational Methods for Single-cell DNA Methylome Analysis

Dissecting intercellular epigenetic differences is key to understanding tissue heterogeneity. Recent advances in single-cell DNA methylome profiling have presented opportunities to resolve this heterogeneity at the maximum resolution. While these advances enable us to explore frontiers of chromatin biology and better understand cell lineage relationships, they pose new challenges in data processing and interpretation. This review surveys the current state of computational tools developed for single-cell DNA methylome data analysis. We discuss critical components of single-cell DNA methylome data analysis, including data preprocessing, quality control, imputation, dimensionality reduction, cell clustering, supervised cell annotation, cell lineage reconstruction, gene activity scoring, and integration with transcriptome data. We also highlight unique aspects of single-cell DNA methylome data analysis and discuss how techniques common to other single-cell omics data analyses can be adapted to analyze DNA methylomes. Finally, we discuss existing challenges and opportunities for future development.

Comparison of Single Molecule, Real-Time Sequencing and Nanopore Sequencing for Analysis of the Size, End-Motif, and Tissue-of-Origin of Long Cell-Free DNA in Plasma

BACKGROUND

Recent studies using single molecule, real-time (SMRT) sequencing revealed a substantial population of analyzable long cell-free DNA (cfDNA) in plasma. Potential clinical utilities of such long cfDNA in pregnancy and cancer have been demonstrated. However, the performance of different long-read sequencing platforms for the analysis of long cfDNA remains unknown.

METHODS

Size biases of SMRT sequencing by Pacific Biosciences (PacBio) and nanopore sequencing by Oxford Nanopore Technologies (ONT) were evaluated using artificial mixtures of sonicated human and mouse DNA of different sizes. cfDNA from plasma samples of pregnant women at different trimesters, hepatitis B carriers, and patients with hepatocellular carcinoma were sequenced with the 2 platforms.

RESULTS

Both platforms showed biases to sequence longer (1500 bp vs 200 bp) DNA fragments, with PacBio showing a stronger bias (5-fold overrepresentation of long fragments vs 2-fold in ONT). Percentages of cfDNA fragments 500 bp were around 6-fold higher in PacBio compared with ONT. End motif profiles of cfDNA from PacBio and ONT were similar, yet exhibited platform-dependent patterns. Tissue-of-origin analysis based on single-molecule methylation patterns showed comparable performance on both platforms.

CONCLUSIONS

SMRT sequencing generated data with higher percentages of long cfDNA compared with nanopore sequencing. Yet, a higher number of long cfDNA fragments eligible for the tissue-of-origin analysis could be obtained from nanopore sequencing due to its much higher throughput. When analyzing the size and end motif of cfDNA, one should be aware of the analytical characteristics and possible biases of the sequencing platforms being used.

Vitamin C boosts DNA demethylation in TET2 germline mutation carriers

BACKGROUND

Accurate regulation of DNA methylation is necessary for normal cells to differentiate, develop and function. TET2 catalyzes stepwise DNA demethylation in hematopoietic cells. Mutations in the TET2 gene predispose to hematological malignancies by causing DNA methylation overload and aberrant epigenomic landscape. Studies on mice and cell lines show that the function of TET2 is boosted by vitamin C. Thus, by strengthening the demethylation activity of TET2, vitamin C could play a role in the prevention of hematological malignancies in individuals with TET2 dysfunction. We recently identified a family with lymphoma predisposition where a heterozygous truncating germline mutation in TET2 segregated with nodular lymphocyte-predominant Hodgkin lymphoma. The mutation carriers displayed a hypermethylation pattern that was absent in the family members without the mutation.

METHODS

In a clinical trial of 1 year, we investigated the effects of oral 1 g/day vitamin C supplementation on DNA methylation by analyzing genome-wide DNA methylation and gene expression patterns from the family members.

RESULTS

We show that vitamin C reinforces the DNA demethylation cascade, reduces the proportion of hypermethylated loci and diminishes gene expression differences between TET2 mutation carriers and control individuals.

CONCLUSIONS

These results suggest that vitamin C supplementation increases DNA methylation turnover and provide a basis for further work to examine the potential benefits of vitamin C supplementation in individuals with germline and somatic TET2 mutations.

TRIAL REGISTRATION

This trial was registered at EudraCT with reference number of 2018-000155-41 (01.04.2019).

Parent-of-origin detection and chromosome-scale haplotyping using long-read DNA methylation sequencing and Strand-seq

Hundreds of loci in human genomes have alleles that are methylated differentially according to their parent of origin. These imprinted loci generally show little variation across tissues, individuals, and populations. We show that such loci can be used to distinguish the maternal and paternal homologs for all human autosomes without the need for the parental DNA. We integrate methylation-detecting nanopore sequencing with the long-range phase information in Strand-seq data to determine the parent of origin of chromosome-length haplotypes for both DNA sequence and DNA methylation in five trios with diverse genetic backgrounds. The parent of origin was correctly inferred for all autosomes with an average mismatch error rate of 0.31% for SNVs and 1.89% for insertions or deletions (indels). Because our method can determine whether an inherited disease allele originated from the mother or the father, we predict that it will improve the diagnosis and management of many genetic diseases.

Dynamic, adaptive sampling during nanopore sequencing using Bayesian experimental design

Nanopore sequencers can select which DNA molecules to sequence, rejecting a molecule after analysis of a small initial part. Currently, selection is based on predetermined regions of interest that remain constant throughout an experiment. Sequencing efforts, thus, cannot be re-focused on molecules likely contributing most to experimental success. Here we present BOSS-RUNS, an algorithmic framework and software to generate dynamically updated decision strategies. We quantify uncertainty at each genome position with real-time updates from data already observed. For each DNA fragment, we decide whether the expected decrease in uncertainty that it would provide warrants fully sequencing it, thus optimizing information gain. BOSS-RUNS mitigates coverage bias between and within members of a microbial community, leading to improved variant calling; for example, low-coverage sites of a species at 1% abundance were reduced by 87.5%, with 12.5% more single-nucleotide polymorphisms detected. Such data-driven updates to molecule selection are applicable to many sequencing scenarios, such as enriching for regions with increased divergence or low coverage, reducing time-to-answer.

Quantifying propagation of DNA methylation and hydroxymethylation with iDEMS

DNA methylation is a critical epigenetic mark in mammalian cells. Many aspects of DNA methylation maintenance have been characterized; however, the exact kinetics of post-replicative methylation maintenance remain a subject of debate. Here we develop isolation of DNA by 5-ethynyl-deoxyuridine labelling for mass spectrometry (iDEMS), a highly sensitive, quantitative mass spectrometry-based method for measuring DNA modifications on metabolically labelled DNA. iDEMS reveals an unexpectedly hemi-methylated landscape on nascent DNA. Combining iDEMS with metabolic labelling reveals that methylation maintenance is outpaced by cell division in mouse embryonic stem cells. Our approach shows that hydroxymethylation is perpetually asymmetric between sister strands in favour of the parental, template strand. iDEMS can be coupled with immunoprecipitation of chromatin proteins, revealing features of DNA methylation-histone modification crosstalk and suggesting a model for interplay between methylation and nucleosome assembly. iDEMS therefore elucidates long-standing questions about DNA modification propagation and provides an important orthogonal technology to understanding this process in dynamic cellular contexts.

Comprehensive variant discovery in the era of complete human reference genomes

Advances in long-read sequencing technologies have broadened our understanding of genetic variation in the human population, uncovered new complex structural variants and offered an opportunity to elucidate new variant associations with disease.

Phasing DNA Methylation

Haplotyping enables the study of allele-specific events. Heterozygous variants, primarily single nucleotide variants (SNVs), enable the assignment of the paternal and maternal origin of the chromosomes and are widely employed to phase sequencing reads to their haplotype of origin. Certain long-read technologies enable the detection of both the DNA sequence and DNA modifications. These long reads and their inherent methylation information are suitable for genome-wide haplotyping and allele-specific DNA methylation analysis. Here, we describe the workflow to phase reads and DNA methylation using nanopore sequencing.

Post-bisulfite Adaptor Tagging Based on an ssDNA Ligation Technique (tPBAT)

Post-bisulfite adaptor tagging (PBAT) is a concept that enables the preparation of an efficient sequencing library from bisulfite-treated DNA, and it also means the protocol implemented the concept. Although the previous PBAT or rPBAT was sensitive enough for single-cell methylome analysis, the protocol had several drawbacks owing to the repeated random priming reactions. To resolve these problems, we developed a unique single-strand DNA ligation technique, termed TACS ligation, and established a new protocol called tPBAT. With tPBAT, the data quality improved, with a longer insert and higher mapping rate than that obtained with rPBAT. In addition, paired-end sequencing and indexing were supported by the default. In this chapter, the tPBAT protocol is introduced, and a thorough description of its application to small samples is provided.

cvlr: finding heterogeneously methylated genomic regions using ONT reads

Summary

Nanopore reads encode information on the methylation status of cytosines in CpG dinucleotides. The length of the reads makes it comparatively easy to look at patterns consisting of multiple loci; here, we exploit this property to search for regions where one can define subpopulations of molecules based on methylation patterns. As an example, we run our clustering algorithm on known imprinted genes; we also scan chromosome 15 looking for windows corresponding to heterogeneous methylation. Our software can also compute the covariance of methylation across these regions while keeping into account the mixture of different types of reads.

Availability and implementation

https://github.com/EmanueleRaineri/cvlr.

Contact

simon.heath@cnag.crg.eu.

Supplementary information

Supplementary data are available at Bioinformatics Advances online.

Targeted Nanopore Resequencing and Methylation Analysis of LINE-1 Retrotransposons

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

Methylation biomarkers for early cancer detection and diagnosis: Current and future perspectives

The increase in recent scientific studies on cancer biomarkers has brought great new insights into the field. Moreover, novel technological breakthroughs such as long read sequencing and microarrays have enabled high throughput profiling of many biomarkers, while advances in bioinformatic tools have made the possibility of developing highly reliable and accurate biomarkers a reality. These changes triggered renewed interest in biomarker research and provided tremendous opportunities for enhancing cancer management and improving early disease detection. DNA methylation alterations are known to accompany and contribute to carcinogenesis, making them promising biomarkers for cancer, namely due to their stability, frequency and accessibility in bodily fluids. The advent of newer minimally invasive experimental methods such as liquid biopsies provide the perfect setting for methylation-based biomarker development and application. Despite their huge potential, accurate and robust biomarkers for the conclusive diagnosis of most cancer types are still not routinely used, hence a strong need for sustained research in this field is still needed. This review provides a brief exposition of current methylation biomarkers for cancer diagnosis and early detection, including markers already in clinical use as well as various upcoming ones. It also outlines how recent big data and novel technologies will revolutionise the next generation of cancer tests in supplementing or replacing currently existing invasive techniques.

Nanopore Sequencing to Identify Transposable Element Insertions and Their Epigenetic Modifications

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

Nanopore Direct RNA Sequencing Data Processing and Analysis Using MasterOfPores

This chapter describes MasterOfPores v.2 (MoP2), an open-source suite of pipelines for processing and analyzing direct RNA Oxford Nanopore sequencing data. The MoP2 relies on the Nextflow DSL2 framework and Linux containers, thus enabling reproducible data analysis in transcriptomic and epitranscriptomic studies. We introduce the key concepts of MoP2 and provide a step-by-step fully reproducible and complete example of how to use the workflow for the analysis of S. cerevisiae total RNA samples sequenced using MinION flowcells. The workflow starts with the pre-processing of raw FAST5 files, which includes basecalling, read quality control, demultiplexing, filtering, mapping, estimation of per-gene/transcript abundances, and transcriptome assembly, with support of the GPU computing for the basecalling and read demultiplexing steps. The secondary analyses of the workflow focus on the estimation of RNA poly(A) tail lengths and the identification of RNA modifications. The MoP2 code is available at https://github.com/biocorecrg/MOP2 and is distributed under the MIT license.

Efficient real-time selective genome sequencing on resource-constrained devices

BACKGROUND

Third-generation nanopore sequencers offer selective sequencing or "Read Until" that allows genomic reads to be analyzed in real time and abandoned halfway if not belonging to a genomic region of "interest." This selective sequencing opens the door to important applications such as rapid and low-cost genetic tests. The latency in analyzing should be as low as possible for selective sequencing to be effective so that unnecessary reads can be rejected as early as possible. However, existing methods that employ a subsequence dynamic time warping (sDTW) algorithm for this problem are too computationally intensive that a massive workstation with dozens of CPU cores still struggles to keep up with the data rate of a mobile phone-sized MinION sequencer.

RESULTS

In this article, we present Hardware Accelerated Read Until (HARU), a resource-efficient hardware-software codesign-based method that exploits a low-cost and portable heterogeneous multiprocessor system-on-chip platform with on-chip field-programmable gate arrays (FPGA) to accelerate the sDTW-based Read Until algorithm. Experimental results show that HARU on a Xilinx FPGA embedded with a 4-core ARM processor is around 2.5× faster than a highly optimized multithreaded software version (around 85× faster than the existing unoptimized multithreaded software) running on a sophisticated server with a 36-core Intel Xeon processor for a SARS-CoV-2 dataset. The energy consumption of HARU is 2 orders of magnitudes lower than the same application executing on the 36-core server.

CONCLUSIONS

HARU demonstrates that nanopore selective sequencing is possible on resource-constrained devices through rigorous hardware-software optimizations. The source code for the HARU sDTW module is available as open source at https://github.com/beebdev/HARU, and an example application that uses HARU is at https://github.com/beebdev/sigfish-haru.

Simultaneous profiling of histone modifications and DNA methylation via nanopore sequencing

The interplay between histone modifications and DNA methylation drives the establishment and maintenance of the cellular epigenomic landscape, but it remains challenging to investigate the complex relationship between these epigenetic marks across the genome. Here we describe a nanopore-sequencing-based-method, nanoHiMe-seq, for interrogating the genome-wide localization of histone modifications and DNA methylation from single DNA molecules. nanoHiMe-seq leverages a nonspecific methyltransferase to exogenously label adenine bases proximal to antibody-targeted modified nucleosomes in situ. The labelled adenines and the endogenous methylated CpG sites are simultaneously detected on individual nanopore reads using a hidden Markov model, which is implemented in the nanoHiMe software package. We demonstrate the utility, robustness and sensitivity of nanoHiMe-seq by jointly profiling DNA methylation and histone modifications at low coverage depths, concurrently determining phased patterns of DNA methylation and histone modifications, and probing the intrinsic connectivity between these epigenetic marks across the genome.

Epitranscriptomics in parasitic protists: Role of RNA chemical modifications in posttranscriptional gene regulation

"Epitranscriptomics" is the new RNA code that represents an ensemble of posttranscriptional RNA chemical modifications, which can precisely coordinate gene expression and biological processes. There are several RNA base modifications, such as N6-methyladenosine (m6A), 5-methylcytosine (m5C), and pseudouridine (Ψ), etc. that play pivotal roles in fine-tuning gene expression in almost all eukaryotes and emerging evidences suggest that parasitic protists are no exception. In this review, we primarily focus on m6A, which is the most abundant epitranscriptomic mark and regulates numerous cellular processes, ranging from nuclear export, mRNA splicing, polyadenylation, stability, and translation. We highlight the universal features of spatiotemporal m6A RNA modifications in eukaryotic phylogeny, their homologs, and unique processes in 3 unicellular parasites-Plasmodium sp., Toxoplasma sp., and Trypanosoma sp. and some technological advances in this rapidly developing research area that can significantly improve our understandings of gene expression regulation in parasites.

Crossing Bacterial Genomic Features and Methylation Patterns with MeStudio: An Epigenomic Analysis Tool

DNA methylation is one of the most observed epigenetic modifications. It is present in eukaryotes and prokaryotes and is related to several biological phenomena, including gene flow and adaptation to environmental conditions. The widespread use of third-generation sequencing technologies allows direct and easy detection of genome-wide methylation profiles, offering increasing opportunities to understand and exploit the epigenomic landscape of individuals and populations. Here, we present a pipeline named MeStudio, with the aim of analyzing and combining genome-wide methylation profiles with genomic features. Outputs report the presence of DNA methylation in coding sequences (CDSs) and noncoding sequences, including both intergenic sequences and sequences upstream of the CDS. We apply this novel tool, showing the usage and performance of MeStudio, on a set of single-molecule real-time sequencing outputs from strains of the bacterial species Sinorhizobium meliloti.

Identification and Mechanism of Action of the Global Secondary Metabolism Regulator SaraC in Stereum hirsutum

DNA methylation is an important factor in the regulation of gene expression. In analyzing genomic data of Stereum hirsutum FP-91666, we found a hypothetical bifunctional transcription regulator/O6Meguanine-DNA methyltransferase (named SaraC), which is widely present in both bacteria and fungi, and confirmed that its function in bacteria is mainly for DNA reparation. In this paper, we confirmed that SaraC has the function of DNA binding and demethylation through surface plasma resonance and reaction experiments in vitro. Then, we achieved the overexpression of SaraC (OES) in S. hirsutum, sequenced the methylation and transcription levels of the whole-genome, and further conducted untargeted metabolomics analyses of the OES transformants and the wild type (WT). The results confirmed that the overall-methylation levels of the transformants were significantly downregulated, and various genes related to secondary metabolism were upregulated. Through comparative untargeted metabolomic analyses, it showed that OES SA6 transformant produced a greater number of hybrid polyketides, and we identified 2 novel hybrid polyketides from the fermentation products of SA6. Our results show that overexpression SaraC can effectively stimulate the expression of secondary-metabolism-related genes, which could be a broad-spectrum tool for discovery of metabolites due to its cross-species conservation. IMPORTANCE Fungi are one of the important sources of active compounds. However, in fungi, most of the secondary metabolic biosynthetic gene clusters are weakly expressed or silenced under conventional culture conditions. How to efficiently excavate potential new compounds contained in fungi is becoming a research hot spot in the world. In this study, we found a DNA demethylation protein (SaraC) and confirmed that it is a global secondary metabolism regulator in Stereum hirsutum FP-91666. In the past, SaraC-like proteins were mainly regarded as DNA repair proteins, but our findings proved that it will be a powerful tool for mining secondary metabolites for overexpression of SaraC, which can effectively stimulate the expression of genes related to secondary metabolism.

Integrated chromatin accessibility and DNA methylation analysis to reveal the critical epigenetic modification and regulatory mechanism in gonadal differentiation of the sequentially hermaphroditic fish, Monopterus albus

BACKGROUND

Monopterus albus is a hermaphroditic and economically farmed fish that undergoes sex reversal from ovary to testis via ovotestis during gonadal development. The epigenetic changes that are associated with gonadal development in this species remain unclear.

METHODS

We produced DNA methylome, transcriptome, and chromatin accessibility maps of the key stages of gonad development: ovary, ovotestis, and testis. The expression of the key candidate genes was detected using qRT-PCR and in situ hybridization and the methylation levels were analysed using bisulphite sequencing PCR. Promoter activity and regulation were assessed using dual-luciferase reporter assays.

RESULTS

Gonadal development exhibits highly dynamic transcriptomic, DNA methylation, and chromatin accessibility changes. We found that DNA methylation status, especially of the transcription start site, was significantly negatively correlated with gene expression while chromatin accessibility exhibited no correlation with gene expression during gonadal development. The epigenetic signatures revealed many novel regulatory elements and genes involved in sex reversal, which were validated. DNA methylation detection and site mutation of plastin-2 promoter, as a candidate gene, revealed that DNA methylation could impact the binding of transcription factor dmrt1 and foxl2 through methylation and demethylation to regulate plastin-2 expression during gonadal development.

CONCLUSIONS

These data provide novel insights into epigenetic modification and help elucidate the potential molecular mechanism by which dynamic modification of DNA methylation plays a crucial role in gonadal development.

Mammalian DNA methylome dynamics: mechanisms, functions and new frontiers

DNA methylation is a highly conserved epigenetic modification that plays essential roles in mammalian gene regulation, genome stability and development. Despite being primarily considered a stable and heritable epigenetic silencing mechanism at heterochromatic and repetitive regions, whole genome methylome analysis reveals that DNA methylation can be highly cell-type specific and dynamic within proximal and distal gene regulatory elements during early embryonic development, stem cell differentiation and reprogramming, and tissue maturation. In this Review, we focus on the mechanisms and functions of regulated DNA methylation and demethylation, highlighting how these dynamics, together with crosstalk between DNA methylation and histone modifications at distinct regulatory regions, contribute to mammalian development and tissue maturation. We also discuss how recent technological advances in single-cell and long-read methylome sequencing, along with targeted epigenome-editing, are enabling unprecedented high-resolution and mechanistic dissection of DNA methylome dynamics.

Single-Cell DNA Methylation Analysis in Cancer

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

Species-specific basecallers improve actual accuracy of nanopore sequencing in plants

BACKGROUND

Long-read sequencing platforms offered by Oxford Nanopore Technologies (ONT) allow native DNA containing epigenetic modifications to be directly sequenced, but can be limited by lower per-base accuracies. A key step post-sequencing is basecalling, the process of converting raw electrical signals produced by the sequencing device into nucleotide sequences. This is challenging as current basecallers are primarily based on mixtures of model species for training. Here we utilise both ONT PromethION and higher accuracy PacBio Sequel II HiFi sequencing on two plants, Phebalium stellatum and Xanthorrhoea johnsonii, to train species-specific basecaller models with the aim of improving per-base accuracy. We investigate sequencing accuracies achieved by ONT basecallers and assess accuracy gains by training single-species and species-specific basecaller models. We also evaluate accuracy gains from ONT's improved flowcells (R10.4, FLO-PRO112) and sequencing kits (SQK-LSK112). For the truth dataset for both model training and accuracy assessment, we developed highly accurate, contiguous diploid reference genomes with PacBio Sequel II HiFi reads.

RESULTS

Basecalling with ONT Guppy 5 and 6 super-accurate gave almost identical results, attaining read accuracies of 91.96% and 94.15%. Guppy's plant-specific model gave highly mixed results, attaining read accuracies of 91.47% and 96.18%. Species-specific basecalling models improved read accuracy, attaining 93.24% and 95.16% read accuracies. R10.4 sequencing kits also improve sequencing accuracy, attaining read accuracies of 95.46% (super-accurate) and 96.87% (species-specific).

CONCLUSIONS

The use of a single mixed-species basecaller model, such as ONT Guppy super-accurate, may be reducing the accuracy of nanopore sequencing, due to conflicting genome biology within the training dataset and study species. Training of single-species and genome-specific basecaller models improves read accuracy. Studies that aim to do large-scale long-read genotyping would primarily benefit from training their own basecalling models. Such studies could use sequencing accuracy gains and improving bioinformatics tools to improve study outcomes.

5-Methyl-cytosine stabilizes DNA but hinders DNA hybridization revealed by magnetic tweezers and simulations

5-Methyl-cytosine (5mC) is one of the most important DNA modifications and plays versatile biological roles. It is well known that 5mC stabilizes DNA duplexes. However, it remains unclear how 5mC affects the kinetics of DNA melting and hybridization. Here, we studied the kinetics of unzipping and rezipping using a 502-bp DNA hairpin by single-molecule magnetic tweezers. Under constant loading rates, 5mC increases the unzipping force but counterintuitively decreases the rezipping force at various salt and temperature conditions. Under constant forces, the non-methylated DNA hops between metastable states during unzipping and rezipping, which implies low energy barriers. Surprisingly, the 5mC DNA can't rezip after fully unzipping unless much lower forces are applied, where it rezips stochastically in a one-step manner, which implies 5mC kinetically hinders DNA hybridization and high energy barriers in DNA hybridization. All-atom molecular dynamics simulations reveal that the 5mC kinetically hinders DNA hybridization due to steric effects rather than electrostatic effects caused by the additional methyl groups of cytosines. Considering the possible high speed of DNA unzipping and zipping during replication and transcription, our findings provide new insights into the biological roles of 5mC.