Gaussian mixture model-based unsupervised nucleotide modification number detection using nanopore-sequencing readouts

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.

Effects of DNA Methylation on Gene Expression and Phenotypic Traits in Cattle: A Review

Gene expression in cells is determined by the epigenetic state of chromatin. Therefore, the study of epigenetic changes is very important to understand the regulatory mechanism of genes at the molecular, cellular, tissue and organ levels. DNA methylation is one of the most studied epigenetic modifications, which plays an important role in maintaining genome stability and ensuring normal growth and development. Studies have shown that methylation levels in bovine primordial germ cells, the rearrangement of methylation during embryonic development and abnormal methylation during placental development are all closely related to their reproductive processes. In addition, the application of bovine male sterility and assisted reproductive technology is also related to DNA methylation. This review introduces the principle, development of detection methods and application conditions of DNA methylation, with emphasis on the relationship between DNA methylation dynamics and bovine spermatogenesis, embryonic development, disease resistance and muscle and fat development, in order to provide theoretical basis for the application of DNA methylation in cattle breeding in the future.

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.

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.

Direct identification of A-to-I editing sites with nanopore native RNA sequencing

Inosine is a prevalent RNA modification in animals and is formed when an adenosine is deaminated by the ADAR family of enzymes. Traditionally, inosines are identified indirectly as variants from Illumina RNA-sequencing data because they are interpreted as guanosines by cellular machineries. However, this indirect method performs poorly in protein-coding regions where exons are typically short, in non-model organisms with sparsely annotated single-nucleotide polymorphisms, or in disease contexts where unknown DNA mutations are pervasive. Here, we show that Oxford Nanopore direct RNA sequencing can be used to identify inosine-containing sites in native transcriptomes with high accuracy. We trained convolutional neural network models to distinguish inosine from adenosine and guanosine, and to estimate the modification rate at each editing site. Furthermore, we demonstrated their utility on the transcriptomes of human, mouse and Xenopus. Our approach expands the toolkit for studying adenosine-to-inosine editing and can be further extended to investigate other RNA modifications.

Concerted modification of nucleotides at functional centers of the ribosome revealed by single-molecule RNA modification profiling

Nucleotides in RNA and DNA are chemically modified by numerous enzymes that alter their function. Eukaryotic ribosomal RNA (rRNA) is modified at more than 100 locations, particularly at highly conserved and functionally important nucleotides. During ribosome biogenesis, modifications are added at various stages of assembly. The existence of differently modified classes of ribosomes in normal cells is unknown because no method exists to simultaneously evaluate the modification status at all sites within a single rRNA molecule. Using a combination of yeast genetics and nanopore direct RNA sequencing, we developed a reliable method to track the modification status of single rRNA molecules at 37 sites in 18 S rRNA and 73 sites in 25 S rRNA. We use our method to characterize patterns of modification heterogeneity and identify concerted modification of nucleotides found near functional centers of the ribosome. Distinct, undermodified subpopulations of rRNAs accumulate upon loss of Dbp3 or Prp43 RNA helicases, suggesting overlapping roles in ribosome biogenesis. Modification profiles are surprisingly resistant to change in response to many genetic and acute environmental conditions that affect translation, ribosome biogenesis, and pre-mRNA splicing. The ability to capture single-molecule RNA modification profiles provides new insights into the roles of nucleotide modifications in RNA function.

Detecting RNA modification using direct RNA sequencing: A systematic review

Post-transcriptional RNA modifications are involved in a range of important cellular processes, including the regulation of gene expression and fine-tuning of the functions of RNA molecules. To decipher the context-specific functions of these post-transcriptional modifications, it is crucial to accurately determine their transcriptomic locations and modification levels under a given cellular condition. With the newly emerged sequencing technology, especially nanopore direct RNA sequencing, different RNA modifications can be detected simultaneously with a single molecular level resolution. Here we provide a systematic review of 15 published RNA modification prediction tools based on direct RNA sequencing data, including their computational models, input–output formats, supported modification types, and reported performances. Finally, we also discussed the potential challenges and future improvements of nanopore sequencing-based methods for RNA modification detection.

Applications and potentials of nanopore sequencing in the (epi)genome and (epi)transcriptome era

The Human Genome Project opened an era of (epi)genomic research, and also provided a platform for the development of new sequencing technologies. During and after the project, several sequencing technologies continue to dominate nucleic acid sequencing markets. Currently, Illumina (short-read), PacBio (long-read), and Oxford Nanopore (long-read) are the most popular sequencing technologies. Unlike PacBio or the popular short-read sequencers before it, which, as examples of the second or so-called Next-Generation Sequencing platforms, need to synthesize when sequencing, nanopore technology directly sequences native DNA and RNA molecules. Nanopore sequencing, therefore, avoids converting mRNA into cDNA molecules, which not only allows for the sequencing of extremely long native DNA and full-length RNA molecules but also document modifications that have been made to those native DNA or RNA bases. In this review on direct DNA sequencing and direct RNA sequencing using Oxford Nanopore technology, we focus on their development and application achievements, discussing their challenges and future perspective. We also address the problems researchers may encounter applying these approaches in their research topics, and how to resolve them.

•

Nanopore-seq can dissect native DNA/RNA molecules from any organisms at unlimited length

•

A wide variety of algorithms greatly increase the accuracy of signal decoding in Nanopore-Seq

•

Nanopore-Seq significantly facilitates genome assembly and structural variant calling, and can simultaneously detect base modifications

•

These advantages ensure its great potentials in future medical and agricultural practices

Towards inferring nanopore sequencing ionic currents from nucleotide chemical structures

The characteristic ionic currents of nucleotide kmers are commonly used in analyzing nanopore sequencing readouts. We present a graph convolutional network-based deep learning framework for predicting kmer characteristic ionic currents from corresponding chemical structures. We show such a framework can generalize the chemical information of the 5-methyl group from thymine to cytosine by correctly predicting 5-methylcytosine-containing DNA 6mers, thus shedding light on the de novo detection of nucleotide modifications.

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

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

Identification of differential RNA modifications from nanopore direct RNA sequencing with xPore

RNA modifications, such as N⁶-methyladenosine (m⁶A), modulate functions of cellular RNA species. However, quantifying differences in RNA modifications has been challenging. Here we develop a computational method, xPore, to identify differential RNA modifications from nanopore direct RNA sequencing (RNA-seq) data. We evaluate our method on transcriptome-wide m⁶A profiling data, demonstrating that xPore identifies positions of m⁶A sites at single-base resolution, estimates the fraction of modified RNA species in the cell and quantifies the differential modification rate across conditions. We apply xPore to direct RNA-seq data from six cell lines and multiple myeloma patient samples without a matched control sample and find that many m⁶A sites are preserved across cell types, whereas a subset exhibit significant differences in their modification rates. Our results show that RNA modifications can be identified from direct RNA-seq data with high accuracy, enabling analysis of differential modifications and expression from a single high-throughput experiment.

DNAscent v2: detecting replication forks in nanopore sequencing data with deep learning

Background

Measuring DNA replication dynamics with high throughput and single-molecule resolution is critical for understanding both the basic biology behind how cells replicate their DNA and how DNA replication can be used as a therapeutic target for diseases like cancer. In recent years, the detection of base analogues in Oxford Nanopore Technologies (ONT) sequencing reads has become a promising new method to supersede existing single-molecule methods such as DNA fibre analysis: ONT sequencing yields long reads with high throughput, and sequenced molecules can be mapped to the genome using standard sequence alignment software.

Results

This paper introduces DNAscent v2, software that uses a residual neural network to achieve fast, accurate detection of the thymidine analogue BrdU with single-nucleotide resolution. DNAscent v2 also comes equipped with an autoencoder that interprets the pattern of BrdU incorporation on each ONT-sequenced molecule into replication fork direction to call the location of replication origins termination sites. DNAscent v2 surpasses previous versions of DNAscent in BrdU calling accuracy, origin calling accuracy, speed, and versatility across different experimental protocols. Unlike NanoMod, DNAscent v2 positively identifies BrdU without the need for sequencing unmodified DNA. Unlike RepNano, DNAscent v2 calls BrdU with single-nucleotide resolution and detects more origins than RepNano from the same sequencing data. DNAscent v2 is open-source and available at https://github.com/MBoemo/DNAscent.

Conclusions

This paper shows that DNAscent v2 is the new state-of-the-art in the high-throughput, single-molecule detection of replication fork dynamics. These improvements in DNAscent v2 mark an important step towards measuring DNA replication dynamics in large genomes with single-molecule resolution. Looking forward, the increase in accuracy in single-nucleotide resolution BrdU calls will also allow DNAscent v2 to branch out into other areas of genome stability research, particularly the detection of DNA repair.

Supplementary Information

The online version contains supplementary material available at (10.1186/s12864-021-07736-6).

DNA Methylation in Solid Tumors: Functions and Methods of Detection

DNA methylation, i.e., addition of methyl group to 5′-carbon of cytosine residues in CpG dinucleotides, is an important epigenetic modification regulating gene expression, and thus implied in many cellular processes. Deregulation of DNA methylation is strongly associated with onset of various diseases, including cancer. Here, we review how DNA methylation affects carcinogenesis process and give examples of solid tumors where aberrant DNA methylation is often present. We explain principles of methods developed for DNA methylation analysis at both single gene and whole genome level, based on (i) sodium bisulfite conversion, (ii) methylation-sensitive restriction enzymes, and (iii) interactions of 5-methylcytosine (5mC) with methyl-binding proteins or antibodies against 5mC. In addition to standard methods, we describe recent advances in next generation sequencing technologies applied to DNA methylation analysis, as well as in development of biosensors that represent their cheaper and faster alternatives. Most importantly, we highlight not only advantages, but also disadvantages and challenges of each method.

Instrumental analysis of RNA modifications

Organisms from all domains of life invest a substantial amount of energy for the introduction of RNA modifications into nearly all transcripts studied to date. Instrumental analysis of RNA can focus on the modified residues and reveal the function of these epitranscriptomic marks. Here, we will review recent advances and breakthroughs achieved by NMR spectroscopy, sequencing, and mass spectrometry of the epitranscriptome.

Analysis of RNA Modifications by Second- and Third-Generation Deep Sequencing: 2020 Update

The precise mapping and quantification of the numerous RNA modifications that are present in tRNAs, rRNAs, ncRNAs/miRNAs, and mRNAs remain a major challenge and a top priority of the epitranscriptomics field. After the keystone discoveries of massive m⁶A methylation in mRNAs, dozens of deep sequencing-based methods and protocols were proposed for the analysis of various RNA modifications, allowing us to considerably extend the list of detectable modified residues. Many of the currently used methods rely on the particular reverse transcription signatures left by RNA modifications in cDNA; these signatures may be naturally present or induced by an appropriate enzymatic or chemical treatment. The newest approaches also include labeling at RNA abasic sites that result from the selective removal of RNA modification or the enhanced cleavage of the RNA ribose-phosphate chain (perhaps also protection from cleavage), followed by specific adapter ligation. Classical affinity/immunoprecipitation-based protocols use either antibodies against modified RNA bases or proteins/enzymes, recognizing RNA modifications. In this survey, we review the most recent achievements in this highly dynamic field, including promising attempts to map RNA modifications by the direct single-molecule sequencing of RNA by nanopores.

Quantitative profiling of N6-methyladenosine at single-base resolution in stem-differentiating xylem of Populus trichocarpa using Nanopore direct RNA sequencing

There are no comprehensive methods to identify N6-methyladenosine (m6A) at single-base resolution for every single transcript, which is necessary for the estimation of m6A abundance. We develop a new pipeline called Nanom6A for the identification and quantification of m6A modification at single-base resolution using Nanopore direct RNA sequencing based on an XGBoost model. We validate our method using methylated RNA immunoprecipitation sequencing (MeRIP-Seq) and m6A-sensitive RNA-endoribonuclease-facilitated sequencing (m6A-REF-seq), confirming high accuracy. Using this method, we provide a transcriptome-wide quantification of m6A modification in stem-differentiating xylem and reveal that different alternative polyadenylation (APA) usage shows a different ratio of m6A.

Probably Correct: Rescuing Repeats with Short and Long Reads

Ever since the introduction of high-throughput sequencing following the human genome project, assembling short reads into a reference of sufficient quality posed a significant problem as a large portion of the human genome-estimated 50-69%-is repetitive. As a result, a sizable proportion of sequencing reads is multi-mapping, i.e., without a unique placement in the genome. The two key parameters for whether or not a read is multi-mapping are the read length and genome complexity. Long reads are now able to span difficult, heterochromatic regions, including full centromeres, and characterize chromosomes from "telomere to telomere". Moreover, identical reads or repeat arrays can be differentiated based on their epigenetic marks, such as methylation patterns, aiding in the assembly process. This is despite the fact that long reads still contain a modest percentage of sequencing errors, disorienting the aligners and assemblers both in accuracy and speed. Here, I review the proposed and implemented solutions to the repeat resolution and the multi-mapping read problem, as well as the downstream consequences of reference choice, repeat masking, and proper representation of sex chromosomes. I also consider the forthcoming challenges and solutions with regards to long reads, where we expect the shift from the problem of repeat localization within a single individual to the problem of repeat positioning within pangenomes.