Methplotlib: analysis of modified nucleotides from nanopore sequencing

Figeno: multi-region genomic figures with long-read support

SUMMARY

The vast amount of publicly available genomic data requires analysis and visualization tools. Here, we present figeno, an application for generating publication-quality FIgures for GENOmics. Figeno particularly focuses on multi-region views across genomic breakpoints and on long reads with base modifications. In addition, we support epigenomic data including ATAC-seq, ChIP-seq or HiC, as well as whole genome sequencing data with copy numbers and structural variants.

AVAILABILITY AND IMPLEMENTATION

Figeno is available as a python package with both a command line and graphical user interface. It can be installed via PyPI and the source code is available at https://github.com/CompEpigen/figeno.

Moving toward the Inclusion of Epigenomics in Bacterial Genome Evolution: Perspectives and Challenges

The universality of DNA methylation as an epigenetic regulatory mechanism belongs to all biological kingdoms. However, while eukaryotic systems have been the primary focus of DNA methylation studies, the molecular mechanisms in prokaryotes are less known. Nevertheless, DNA methylation in prokaryotes plays a pivotal role in many cellular processes such as defense systems against exogenous DNA, cell cycle dynamics, and gene expression, including virulence. Thanks to single-molecule DNA sequencing technologies, genome-wide identification of methylated DNA is becoming feasible on a large scale, providing the possibility to investigate more deeply the presence, variability, and roles of DNA methylation. Here, we present an overview of the multifaceted roles of DNA methylation in prokaryotes and suggest research directions and tools which can enable us to better understand the contribution of DNA methylation to prokaryotic genome evolution and adaptation. In particular, we emphasize the need to understand the presence and role of transgenerational inheritance, as well as the impact of epigenomic signatures on adaptation and genome evolution. Research directions and the importance of novel computational tools are underlined.

Near telomere-to-telomere genome assemblies of two Chlorella species unveil the composition and evolution of centromeres in green algae

BACKGROUND

Centromeres play a crucial and conserved role in cell division, although their composition and evolutionary history in green algae, the evolutionary ancestors of land plants, remains largely unknown.

RESULTS

We constructed near telomere-to-telomere (T2T) assemblies for two Trebouxiophyceae species, Chlorella sorokiniana NS4-2 and Chlorella pyrenoidosa DBH, with chromosome numbers of 12 and 13, and genome sizes of 58.11 Mb and 53.41 Mb, respectively. We identified and validated their centromere sequences using CENH3 ChIP-seq and found that, similar to humans and higher plants, the centromeric CENH3 signals of green algae display a pattern of hypomethylation. Interestingly, the centromeres of both species largely comprised transposable elements, although they differed significantly in their composition. Species within the Chlorella genus display a more diverse centromere composition, with major constituents including members of the LTR/Copia, LINE/L1, and LINE/RTEX families. This is in contrast to green algae including Chlamydomonas reinhardtii, Coccomyxa subellipsoidea, and Chromochloris zofingiensis, in which centromere composition instead has a pronounced single-element composition. Moreover, we observed significant differences in the composition and structure of centromeres among chromosomes with strong collinearity within the Chlorella genus, suggesting that centromeric sequence evolves more rapidly than sequence in non-centromeric regions.

CONCLUSIONS

This study not only provides high-quality genome data for comparative genomics of green algae but gives insight into the composition and evolutionary history of centromeres in early plants, laying an important foundation for further research on their evolution.

CGG repeat expansion in LRP12 in amyotrophic lateral sclerosis

Amyotrophic lateral sclerosis (ALS) is a neurodegenerative disorder characterized by the degeneration of motor neurons. Although repeat expansion in C9orf72 is its most common cause, the pathogenesis of ALS isn't fully clear. In this study, we show that repeat expansion in LRP12, a causative variant of oculopharyngodistal myopathy type 1 (OPDM1), is a cause of ALS. We identify CGG repeat expansion in LRP12 in five families and two simplex individuals. These ALS individuals (LRP12-ALS) have 61-100 repeats, which contrasts with most OPDM individuals with repeat expansion in LRP12 (LRP12-OPDM), who have 100-200 repeats. Phosphorylated TDP-43 is present in the cytoplasm of iPS cell-derived motor neurons (iPSMNs) in LRP12-ALS, a finding that reproduces the pathological hallmark of ALS. RNA foci are more prominent in muscle and iPSMNs in LRP12-ALS than in LRP12-OPDM. Muscleblind-like 1 aggregates are observed only in OPDM muscle. In conclusion, CGG repeat expansions in LRP12 cause ALS and OPDM, depending on the length of the repeat. Our findings provide insight into the repeat length-dependent switching of phenotypes.

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.

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.

Methylartist: Tools for Visualising Modified Bases from Nanopore Sequence Data

Summary

Methylartist is a consolidated suite of tools for processing, visualizing and analysing nanopore-derived modified base calls. All detectable methylation types (e.g. 5mCpG, 5hmC, 6mA) are supported, enabling integrated study of base pairs when modified naturally or as part of an experimental protocol.

Availability and implementation

Methylartist is implemented in Python and is installable via PyPI and bioconda. Source code and test data are available at https://github.com/adamewing/methylartist.

Supplementary information

Supplementary data are available at Bioinformatics online.

ModPhred: an integrative toolkit for the analysis and storage of nanopore sequencing DNA and RNA modification data

MOTIVATION

DNA and RNA modifications can now be identified using nanopore sequencing. However, we currently lack a flexible software to efficiently encode, store, analyze and visualize DNA and RNA modification data.

RESULTS

Here, we present ModPhred, a versatile toolkit that facilitates DNA and RNA modification analysis from nanopore sequencing reads in a user-friendly manner. ModPhred integrates probabilistic DNA and RNA modification information within the FASTQ and BAM file formats, can be used to encode multiple types of modifications simultaneously, and its output can be easily coupled to genomic track viewers, facilitating the visualization and analysis of DNA and RNA modification information in individual reads in a simple and computationally efficient manner.

AVAILABILITY AND IMPLEMENTATION

ModPhred is available at https://github.com/novoalab/modPhred, is implemented in Python3, and is released under an MIT license. Docker images with all dependencies preinstalled are also provided.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Nanopore sequencing technology, bioinformatics and applications

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

NanoMethViz: An R/Bioconductor package for visualizing long-read methylation data

A key benefit of long-read nanopore sequencing technology is the ability to detect modified DNA bases, such as 5-methylcytosine. The lack of R/Bioconductor tools for the effective visualization of nanopore methylation profiles between samples from different experimental groups led us to develop the NanoMethViz R package. Our software can handle methylation output generated from a range of different methylation callers and manages large datasets using a compressed data format. To fully explore the methylation patterns in a dataset, NanoMethViz allows plotting of data at various resolutions. At the sample-level, we use dimensionality reduction to look at the relationships between methylation profiles in an unsupervised way. We visualize methylation profiles of classes of features such as genes or CpG islands by scaling them to relative positions and aggregating their profiles. At the finest resolution, we visualize methylation patterns across individual reads along the genome using the spaghetti plot and heatmaps, allowing users to explore particular genes or genomic regions of interest. In summary, our software makes the handling of methylation signal more convenient, expands upon the visualization options for nanopore data and works seamlessly with existing methylation analysis tools available in the Bioconductor project. Our software is available at https://bioconductor.org/packages/NanoMethViz.

Rapid Characterization of Complex Killer Cell Immunoglobulin-Like Receptor (KIR) Regions Using Cas9 Enrichment and Nanopore Sequencing

Long-read sequencing approaches have considerably improved the quality and contiguity of genome assemblies. Such platforms bear the potential to resolve even extremely complex regions, such as multigenic immune families and repetitive stretches of DNA. Deep sequencing coverage, however, is required to overcome low nucleotide accuracy, especially in regions with high homopolymer density, copy number variation, and sequence similarity, such as the MHC and KIR gene clusters of the immune system. Therefore, we have adapted a targeted enrichment protocol in combination with long-read sequencing to efficiently annotate complex KIR gene regions. Using Cas9 endonuclease activity, segments of the KIR gene cluster were enriched and sequenced on an Oxford Nanopore Technologies platform. This provided sufficient coverage to accurately resolve and phase highly complex KIR haplotypes. Our strategy eliminates PCR-induced amplification errors, facilitates rapid characterization of large and complex multigenic regions, including its epigenetic footprint, and is applicable in multiple species, even in the absence of a reference genome.

High-quality Arabidopsis thaliana Genome Assembly with Nanopore and HiFi Long Reads

Arabidopsis thaliana is an important and long-established model species for plant molecular biology, genetics, epigenetics, and genomics. However, the latest version of reference genome still contains significant number of missing segments. Here, we report a high-quality and almost complete Col-0 genome assembly with two gaps (Col-XJTU) using combination of Oxford Nanopore Technology ultra-long reads, PacBio high-fidelity long reads, and Hi-C data. The total genome assembly size is 133,725,193 bp, introducing 14.6 Mb of novel sequences compared to the TAIR10.1 reference genome. All five chromosomes of Col-XJTU assembly are highly accurate with consensus quality (QV) scores > 60 (ranging from 62 to 68), which are higher than those of TAIR10.1 reference (QV scores ranging from 45 to 52). We have completely resolved chromosome (Chr) 3 and Chr5 in a telomere-to-telomere manner. Chr4 has been completely resolved except the nucleolar organizing regions, which comprise long repetitive DNA fragments. The Chr1 centromere (CEN1), reportedly around 9 Mb in length, is particularly challenging to assemble due to the presence of tens of thousands of CEN180 satellite repeats. Using the cutting-edge sequencing data and novel computational approaches, we assembled about 4 Mb of sequence for CEN1 and a 3.5-Mb-long CEN2. We investigated the structure and epigenetics of centromeres. We detected four clusters of CEN180 monomers, and found that the centromere-specific histone H3-like protein (CENH3) exhibits a strong preference for CEN180 cluster 3. Moreover, we observed hypomethylation patterns in CENH3-enriched regions. We believe that this high-quality genome assembly, Col-XJTU, would serve as a valuable reference to better understand the global pattern of centromeric polymorphisms, as well as genetic and epigenetic features in plants.

DNA METHYLATION PROFILING OF A CNIDARIAN-ALGAL SYMBIOSIS USING NANOPORE SEQUENCING

Symbiosis with protists is common among cnidarians such as corals and sea anemones and is associated with homeostatic and phenotypic changes in the host that could have epigenetic underpinnings, such as methylation of CpG dinucleotides. We leveraged the sensitivity to base modifications of nanopore sequencing to probe the effect of symbiosis with the chlorophyte Elliptochloris marina on methylation in the sea anemone Anthopleura elegantissima. We first validated the approach by comparison of nanopore-derived methylation levels with CpG depletion analysis of a published transcriptome, finding that high methylation levels are associated with CpG depletion as expected. Next, using reads generated exclusively from aposymbiotic anemones, a largely complete draft genome comprising 243 Mb was assembled. Reads from aposymbiotic and symbiotic sea anemones were then mapped to this genome and assessed for methylation using the program Nanopolish, which detects signal disruptions from base modifications as they pass through the nanopore. Based on assessment of 452,841 CpGs for which there was adequate read coverage (approximately 8% of the CpGs in the genome), symbiosis with E. marina was, surprisingly, associated with only subtle changes in the host methylome. However, we did identify one extended genomic region with consistently higher methylation among symbiotic individuals. The region was associated with a DNA polymerase zeta that is noted for its role in translesion synthesis, which opens interesting questions about the biology of this symbiosis. Our study highlights the power and relative simplicity of nanopore sequencing for studies of nucleic acid base modifications in non-model species.

Quantitative profiling of pseudouridylation dynamics in native RNAs with nanopore sequencing

Nanopore RNA sequencing shows promise as a method for discriminating and identifying different RNA modifications in native RNA. Expanding on the ability of nanopore sequencing to detect N⁶-methyladenosine, we show that other modifications, in particular pseudouridine (Ψ) and 2'-O-methylation (Nm), also result in characteristic base-calling 'error' signatures in the nanopore data. Focusing on Ψ modification sites, we detected known and uncovered previously unreported Ψ sites in mRNAs, non-coding RNAs and rRNAs, including a Pus4-dependent Ψ modification in yeast mitochondrial rRNA. To explore the dynamics of pseudouridylation, we treated yeast cells with oxidative, cold and heat stresses and detected heat-sensitive Ψ-modified sites in small nuclear RNAs, small nucleolar RNAs and mRNAs. Finally, we developed a software, nanoRMS, that estimates per-site modification stoichiometries by identifying single-molecule reads with altered current intensity and trace profiles. This work demonstrates that Nm and Ψ RNA modifications can be detected in cellular RNAs and that their modification stoichiometry can be quantified by nanopore sequencing of native RNA.

DNAModAnnot: a R toolbox for DNA modification filtering and annotation

MOTIVATION

Long-read sequencing technologies can be employed to detect and map DNA modifications at the nucleotide resolution on a genome-wide scale. However, published software packages neglect the integration of genomic annotation and comprehensive filtering when analyzing patterns of modified bases detected using Pacific Biosciences (PacBio) or Oxford Nanopore Technologies (ONT) data. Here, we present DNAModAnnot, a R package designed for the global analysis of DNA modification patterns using adapted filtering and visualization tools.

RESULTS

We tested our package using PacBio sequencing data to analyze patterns of the 6-methyladenine (6 mA) in the ciliate Paramecium tetraurelia, in which high 6 mA amounts were previously reported. We found Paramecium tetraurelia 6 mA genome-wide distribution to be similar to other ciliates. We also performed 5-methylcytosine (5mC) analysis in human lymphoblastoid cells using ONT data and confirmed previously known patterns of 5mC. DNAModAnnot provides a toolbox for the genome-wide analysis of different DNA modifications using PacBio and ONT long-read sequencing data.

AVAILABILITY

DNAModAnnot is distributed as a R package available via GitHub (https://github.com/AlexisHardy/DNAModAnnot).

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.