Detecting DNA cytosine methylation using nanopore sequencing

Constricted Apertures for Dynamic Trapping and Micro-/Nanoscale Discrimination Based on Recapture Kinetics

Sensing via analyte passage through a constricted aperture is a powerful and robust technology which is being utilized broadly, from DNA sequencing to single virus and cell characterization. Micro- and nanoscale structures typically translocate a constricted aperture, or pore, using electrophoretic force. In the present work, we explore the advances in metrology which can be achieved through rapid directional switching of hydrodynamic forces. Interestingly, multipass measurements of microscale and nanoscale structures achieve cell discrimination. We explore this cell-discrimination phenomenon as well as other features of hydrodynamic focusing such as dynamic trapping and discrete interval sensing.

Transcriptional and epi-transcriptional dynamics of SARS-CoV-2 during cellular infection

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) uses subgenomic RNA (sgRNA) to produce viral proteins for replication and immune evasion. We apply long-read RNA and cDNA sequencing to in vitro human and primate infection models to study transcriptional dynamics. Transcription-regulating sequence (TRS)-dependent sgRNA upregulates earlier in infection than TRS-independent sgRNA. An abundant class of TRS-independent sgRNA consisting of a portion of open reading frame 1ab (ORF1ab) containing nsp1 joins to ORF10, and the 3' untranslated region (UTR) upregulates at 48 h post-infection in human cell lines. We identify double-junction sgRNA containing both TRS-dependent and -independent junctions. We find multiple sites at which the SARS-CoV-2 genome is consistently more modified than sgRNA and that sgRNA modifications are stable across transcript clusters, host cells, and time since infection. Our work highlights the dynamic nature of the SARS-CoV-2 transcriptome during its replication cycle.

Transcriptional Silencing of ALDH2 Confers a Dependency on Fanconi Anemia Proteins in Acute Myeloid Leukemia

Hundreds of genes become aberrantly silenced in acute myeloid leukemia (AML), with most of these epigenetic changes being of unknown functional consequence. Here, we demonstrate how gene silencing can lead to an acquired dependency on the DNA repair machinery in AML. We make this observation by profiling the essentiality of the ubiquitination machinery in cancer cell lines using domain-focused CRISPR screening, which revealed Fanconi anemia (FA) proteins UBE2T and FANCL as unique dependencies in AML. We demonstrate that these dependencies are due to a synthetic lethal interaction between FA proteins and aldehyde dehydrogenase 2 (ALDH2), which function in parallel pathways to counteract the genotoxicity of endogenous aldehydes. We show DNA hypermethylation and silencing of ALDH2 occur in a recurrent manner in human AML, which is sufficient to confer FA pathway dependency. Our study suggests that targeting of the ubiquitination reaction catalyzed by FA proteins can eliminate ALDH2-deficient AML. SIGNIFICANCE: Aberrant gene silencing is an epigenetic hallmark of human cancer, but the functional consequences of this process are largely unknown. In this study, we show how an epigenetic alteration leads to an actionable dependency on a DNA repair pathway through the disabling of genetic redundancy.This article is highlighted in the In This Issue feature, p. 2113.

Is Oxford Nanopore sequencing ready for analyzing complex microbiomes?

This minireview will discuss the improvements in Oxford Nanopore (Oxford; sequencing technology that make the MinION a viable platform for microbial ecology studies. Specific issues being addressed are the increase in sequence accuracy from 65 to 96.5% during the last 5 years, the ability to obtain a quantifiable/predictive signal from the MinION with respect to target molecule abundance, simple-to-use GUI-based pathways for data analysis and the modest additional equipment needs for sequencing in the field. Coupling these recent improvements with the low capital costs for equipment and the reasonable per sample cost makes MinION sequencing an attractive option for virtually any laboratory.

Multi-omics data integration considerations and study design for biological systems and disease

With the advancement of next-generation sequencing and mass spectrometry, there is a growing need for the ability to merge biological features in order to study a system as a whole. Features such as the transcriptome, methylome, proteome, histone post-translational modifications and the microbiome all influence the host response to various diseases and cancers. Each of these platforms have technological limitations due to sample preparation steps, amount of material needed for sequencing, and sequencing depth requirements. These features provide a snapshot of one level of regulation in a system. The obvious next step is to integrate this information and learn how genes, proteins, and/or epigenetic factors influence the phenotype of a disease in context of the system. In recent years, there has been a push for the development of data integration methods. Each method specifically integrates a subset of omics data using approaches such as conceptual integration, statistical integration, model-based integration, networks, and pathway data integration. In this review, we discuss considerations of the study design for each data feature, the limitations in gene and protein abundance and their rate of expression, the current data integration methods, and microbiome influences on gene and protein expression. The considerations discussed in this review should be regarded when developing new algorithms for integrating multi-omics data.

Low biological fluctuation of mitochondrial CpG and non-CpG methylation at the single-molecule level

Mammalian cytosine DNA methylation (5mC) is associated with the integrity of the genome and the transcriptional status of nuclear DNA. Due to technical limitations, it has been less clear if mitochondrial DNA (mtDNA) is methylated and whether 5mC has a regulatory role in this context. Here, we used bisulfite-independent single-molecule sequencing of native human and mouse DNA to study mitochondrial 5mC across different biological conditions. We first validated the ability of long-read nanopore sequencing to detect 5mC in CpG (5mCpG) and non-CpG (5mCpH) context in nuclear DNA at expected genomic locations (i.e. promoters, gene bodies, enhancers, and cell type-specific transcription factor binding sites). Next, using high coverage nanopore sequencing we found low levels of mtDNA CpG and CpH methylation (with several exceptions) and little variation across biological processes: differentiation, oxidative stress, and cancer. 5mCpG and 5mCpH were overall higher in tissues compared to cell lines, with small additional variation between cell lines of different origin. Despite general low levels, global and single-base differences were found in cancer tissues compared to their adjacent counterparts, in particular for 5mCpG. In conclusion, nanopore sequencing is a useful tool for the detection of modified DNA bases on mitochondria that avoid the biases introduced by bisulfite and PCR amplification. Enhanced nanopore basecalling models will provide further resolution on the small size effects detected here, as well as rule out the presence of other DNA modifications such as oxidized forms of 5mC.

Genome-Wide Methylation Mapping Using Nanopore Sequencing Technology Identifies Novel Tumor Suppressor Genes in Hepatocellular Carcinoma

Downregulation of multiple tumor suppressor genes (TSGs) plays an important role in cancer formation. Recent evidence has accumulated that cancer progression involves genome-wide alteration of epigenetic modifications, which may cause downregulation of the tumor suppressor gene. Using hepatocellular carcinoma (HCC) as a system, we mapped 5-methylcytosine signal at a genome-wide scale using nanopore sequencing technology to identify novel TSGs. Integration of methylation data with gene transcription profile of regenerated liver and primary HCCs allowed us to identify 10 potential tumor suppressor gene candidates. Subsequent validation led us to focus on functionally characterizing one candidate—glucokinase (GCK). We show here that overexpression of GCK inhibits the proliferation of HCC cells via induction of intracellular lactate accumulation and subsequently causes energy crisis due to NAD+ depletion. This suggests GCK functions as a tumor suppressor gene and may be involved in HCC development. In conclusion, these data provide valuable clues for further investigations of the process of tumorigenesis in human cancer.

Discovering multiple types of DNA methylation from bacteria and microbiome using nanopore sequencing

Bacterial DNA methylation occurs at diverse sequence contexts and plays important functional roles in cellular defense and gene regulation. Existing methods for detecting DNA modification from nanopore sequencing data do not effectively support de novo study of unknown bacterial methylomes. In this work, we observed that a nanopore sequencing signal displays complex heterogeneity across methylation events of the same type. To enable nanopore sequencing for broadly applicable methylation discovery, we generated a training dataset from an assortment of bacterial species and developed a method, named nanodisco ( https://github.com/fanglab/nanodisco ), that couples the identification and fine mapping of the three forms of methylation into a multi-label classification framework. We applied it to individual bacteria and the mouse gut microbiome for reliable methylation discovery. In addition, we demonstrated the use of DNA methylation for binning metagenomic contigs, associating mobile genetic elements with their host genomes and identifying misassembled metagenomic contigs.

Whole-genome characterization of chronological age-associated changes in methylome and circular RNAs in moso bamboo (Phyllostachys edulis) from vegetative to floral growth

In mammals, DNA methylation is associated with aging. However, age-related DNA methylation changes during phase transitions largely remain unstudied in plants. Moso bamboo (Phyllostachys edulis) requires a very long time to transition from the vegetative to the floral phase. To comprehensively investigate the association of DNA methylation with aging, we present here single-base-resolution DNA methylation profiles using both high-throughput bisulfite sequencing and single-molecule nanopore-based DNA sequencing, covering the long period of vegetative growth and transition to flowering in moso bamboo. We discovered that CHH methylation gradually accumulates from vegetative to reproductive growth in a time-dependent fashion. Differentially methylated regions, correlating with chronological aging, occurred preferentially at both transcription start sites and transcription termination sites. Genes with CG methylation changes showed an enrichment of Gene Ontology (GO) categories in 'vegetative to reproductive phase transition of meristem'. Combining methylation data with mRNA sequencing revealed that DNA methylation in promoters, introns and exons may have different roles in regulating gene expression. Finally, circular RNA (circRNA) sequencing revealed that the flanking introns of circRNAs are hypermethylated and enriched in long terminal repeat (LTR) retrotransposons. Together, the observations in this study provide insights into the dynamic DNA methylation and circRNA landscapes, correlating with chronological age, which paves the way to study further the impact of epigenetic factors on flowering in moso bamboo.

Targeted nanopore sequencing by real-time mapping of raw electrical signal with UNCALLED

Conventional targeted sequencing methods eliminate many of the benefits of nanopore sequencing, such as the ability to accurately detect structural variants or epigenetic modifications. The ReadUntil method allows nanopore devices to selectively eject reads from pores in real time, which could enable purely computational targeted sequencing. However, this requires rapid identification of on-target reads while most mapping methods require computationally intensive basecalling. We present UNCALLED ( https://github.com/skovaka/UNCALLED ), an open source mapper that rapidly matches streaming of nanopore current signals to a reference sequence. UNCALLED probabilistically considers k-mers that could be represented by the signal and then prunes the candidates based on the reference encoded within a Ferragina-Manzini index. We used UNCALLED to deplete sequencing of known bacterial genomes within a metagenomics community, enriching the remaining species 4.46-fold. UNCALLED also enriched 148 human genes associated with hereditary cancers to 29.6× coverage using one MinION flowcell, enabling accurate detection of single-nucleotide polymorphisms, insertions and deletions, structural variants and methylation in these genes.

Applications of Oxford Nanopore Sequencing in Schizosaccharomyces pombe

Recent years have seen great progresses in third-generation sequencing. New commercial platforms from Oxford Nanopore Technologies (ONT) can generate ultra-long reads from single-molecule nucleic acid fragments of kilobases up to megabases, exceeding the limitation of short reads and dependency on template amplification suffered by the previous generation of sequencing technologies. Moreover, it can detect epigenetic modifications directly, as well as providing all-around field usage, being pocket-sized and low cost. It has already been applied to yeast research in many aspects, such as complete de novo genome assemblies, the phylogeny of large-brewing yeasts, gene isoform identification, and base modification detection. These applications have delivered novel insights into yeast genomic and transcriptomic analysis.

Hybrid genome de novo assembly with methylome analysis of the anaerobic thermophilic subsurface bacterium Thermanaerosceptrum fracticalcis strain DRI-13T

BACKGROUND

There is a dearth of sequenced and closed microbial genomes from environments that exceed > 500 m below level terrestrial surface. Coupled with even fewer cultured isolates, study and understanding of how life endures in the extreme oligotrophic subsurface environments is greatly hindered. Using a de novo hybrid assembly of Illumina and Oxford Nanopore sequences we produced a circular genome with corresponding methylome profile of the recently characterized thermophilic, anaerobic, and fumarate-respiring subsurface bacterium, Thermanaerosceptrum fracticalcis, strain DRI-13T to understand how this microorganism survives the deep subsurface.

RESULTS

The hybrid assembly produced a single circular genome of 3.8 Mb in length with an overall GC content of 45%. Out of the total 4022 annotated genes, 3884 are protein coding, 87 are RNA encoding genes, and the remaining 51 genes were associated with regulatory features of the genome including riboswitches and T-box leader sequences. Approximately 24% of the protein coding genes were hypothetical. Analysis of strain DRI-13T genome revealed: 1) energy conservation by bifurcation hydrogenase when growing on fumarate, 2) four novel bacterial prophages, 3) methylation profile including 76.4% N6-methyladenine and 3.81% 5-methylcytosine corresponding to novel DNA methyltransferase motifs. As well a cluster of 45 genes of unknown protein families that have enriched DNA mCpG proximal to the transcription start sites, and 4) discovery of a putative core of bacteriophage exclusion (BREX) genes surrounded by hypothetical proteins, with predicted functions as helicases, nucleases, and exonucleases.

CONCLUSIONS

The de novo hybrid assembly of strain DRI-13T genome has provided a more contiguous and accurate view of the subsurface bacterium T. fracticalcis, strain DRI-13T. This genome analysis reveals a physiological focus supporting syntrophy, non-homologous double stranded DNA repair, mobility/adherence/chemotaxis, unique methylome profile/recognized motifs, and a BREX defense system. The key to microbial subsurface survival may not rest on genetic diversity, but rather through specific syntrophy niches and novel methylation strategies.

Epigenetically regulated digital signaling defines epithelial innate immunity at the tissue level

To prevent damage to the host or its commensal microbiota, epithelial tissues must match the intensity of the immune response to the severity of a biological threat. Toll-like receptors allow epithelial cells to identify microbe associated molecular patterns. However, the mechanisms that mitigate biological noise in single cells to ensure quantitatively appropriate responses remain unclear. Here we address this question using single cell and single molecule approaches in mammary epithelial cells and primary organoids. We find that epithelial tissues respond to bacterial microbe associated molecular patterns by activating a subset of cells in an all-or-nothing (i.e. digital) manner. The maximum fraction of responsive cells is regulated by a bimodal epigenetic switch that licenses the TLR2 promoter for transcription across multiple generations. This mechanism confers a flexible memory of inflammatory events as well as unique spatio-temporal control of epithelial tissue-level immune responses. We propose that epigenetic licensing in individual cells allows for long-term, quantitative fine-tuning of population-level responses.

Direct Nanopore Sequencing of mRNA Reveals Landscape of Transcript Isoforms in Apicomplexan Parasites

Alternative splicing is a widespread phenomenon in metazoans by which single genes are able to produce multiple isoforms of the gene product. However, this has been poorly characterized in apicomplexans, a major phylum of some of the most important global parasites. Efforts have been hampered by atypical transcriptomic features, such as the high AU content of Plasmodium RNA, but also the limitations of short-read sequencing in deciphering complex splicing events. In this study, we utilized the long read direct RNA sequencing platform developed by Oxford Nanopore Technologies to survey the alternative splicing landscape of Toxoplasma gondii and Plasmodium falciparum. We find that while native RNA sequencing has a reduced throughput, it allows us to obtain full-length or nearly full-length transcripts with comparable quantification to Illumina sequencing. By comparing these data with available gene models, we find widespread alternative splicing, particularly intron retention, in these parasites. Most of these transcripts contain premature stop codons, suggesting that in these parasites, alternative splicing represents a pathway to transcriptomic diversity, rather than expanding proteomic diversity. Moreover, alternative splicing rates are comparable between parasites, suggesting a shared splicing machinery, despite notable transcriptomic differences between the parasites. This study highlights a strategy in using long-read sequencing to understand splicing events at the whole-transcript level and has implications in the future interpretation of transcriptome sequencing studies.

IMPORTANCE We have used a novel nanopore sequencing technology to directly analyze parasite transcriptomes. The very long reads of this technology reveal the full-length genes of the parasites that cause malaria and toxoplasmosis. Gene transcripts must be processed in a process called splicing before they can be translated to protein. Our analysis reveals that these parasites very frequently only partially process their gene products, in a manner that departs dramatically from their human hosts.

The GGC repeat expansion in NOTCH2NLC is associated with oculopharyngodistal myopathy type 3

Oculopharyngodistal myopathy (OPDM) is an adult-onset neuromuscular disease characterized by progressive ocular, facial, pharyngeal and distal limb muscle involvement. Trinucleotide repeat expansions in LRP12 or GIPC1 were recently reported to be associated with OPDM. However, a significant portion of OPDM patients have unknown genetic causes. In this study, long-read whole-genome sequencing and repeat-primed PCR were performed and we identified GGC repeat expansions in the NOTCH2NLC gene in 16.7% (4/24) of a cohort of Chinese OPDM patients, designated as OPDM type 3 (OPDM3). Methylation analysis indicated that methylation levels of the NOTCH2NLC gene were unaltered in OPDM3 patients, but increased significantly in asymptomatic carriers. Quantitative real-time PCR analysis indicated that NOTCH2NLC mRNA levels were increased in muscle but not in blood of OPDM3 patients. Immunofluorescence on OPDM muscle samples and expressing mutant NOTCH2NLC with (GGC)₆₉ repeat expansions in HEK293 cells indicated that mutant NOTCH2NLC-polyglycine protein might be a major component of intranuclear inclusions, and contribute to toxicity in cultured cells. In addition, two RNA-binding proteins, hnRNP A/B and MBNL1, were both co-localized with p62 in intranuclear inclusions in OPDM muscle samples. These results indicated that a toxic protein gain-of-function mechanism and RNA gain-of-function mechanism may both play a vital role in the pathogenic processes of OPDM3. This study extended the spectrum of NOTCH2NLC repeat expansion-related diseases to a predominant myopathy phenotype presenting as OPDM, and provided evidence for possible pathogenesis of these diseases.

Nanopore sequencing and its application to the study of microbial communities

Since its introduction, nanopore sequencing has enhanced our ability to study complex microbial samples through the possibility to sequence long reads in real time using inexpensive and portable technologies. The use of long reads has allowed to address several previously unsolved issues in the field, such as the resolution of complex genomic structures, and facilitated the access to metagenome assembled genomes (MAGs). Furthermore, the low cost and portability of platforms together with the development of rapid protocols and analysis pipelines have featured nanopore technology as an attractive and ever-growing tool for real-time in-field sequencing for environmental microbial analysis. This review provides an up-to-date summary of the experimental protocols and bioinformatic tools for the study of microbial communities using nanopore sequencing, highlighting the most important and recent research in the field with a major focus on infectious diseases. An overview of the main approaches including targeted and shotgun approaches, metatranscriptomics, epigenomics, and epitranscriptomics is provided, together with an outlook to the major challenges and perspectives over the use of this technology for microbial studies.

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

MOTIVATION

Nucleotide modification status can be decoded from the Oxford Nanopore Technologies nanopore-sequencing ionic current signals. Although various algorithms have been developed for nanopore-sequencing-based modification analysis, more detailed characterizations, such as modification numbers, corresponding signal levels and proportions are still lacking.

RESULTS

We present a framework for the unsupervised determination of the number of nucleotide modifications from nanopore-sequencing readouts. We demonstrate the approach can effectively recapitulate the number of modifications, the corresponding ionic current signal levels, as well as mixing proportions under both DNA and RNA contexts. We further show, by integrating information from multiple detected modification regions, that the modification status of DNA and RNA molecules can be inferred. This method forms a key step of de novo characterization of nucleotide modifications, shedding light on the interpretation of various biological questions.

AVAILABILITY AND IMPLEMENTATION

Modified nanopolish: https://github.com/adbailey4/nanopolish/tree/cigar_output. All other codes used to reproduce the results: https://github.com/hd2326/ModificationNumber.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Cellular hallmarks of aging emerge in the ovary prior to primordial follicle depletion

Decline in ovarian reserve with advancing age is associated with reduced fertility and the emergence of metabolic disturbances, osteoporosis, and neurodegeneration. Recent studies have provided insight into connections between ovarian insufficiency and systemic aging, although the basic mechanisms that promote ovarian reserve depletion remain unknown. Here, we sought to determine if chronological age is linked to changes in ovarian cellular senescence, transcriptomic, and epigenetic mechanisms in a mouse model. Histological assessments and transcriptional analyses revealed the accumulation of lipofuscin aggresomes and senescence-related transcripts (Cdkn1a, Cdkn2a, Pai-1 and Hmgb1) significantly increased with advancing age. Transcriptomic profiling and pathway analyses following RNA sequencing, revealed an upregulation of genes related to pro-inflammatory stress and cell-cycle inhibition, whereas genes involved in cell-cycle progression were downregulated; which could be indicative of senescent cell accumulation. The emergence of these senescence-related markers preceded the dramatic decline in primordial follicle reserve observed. Whole Genome Oxidative Bisulfite Sequencing (WGoxBS) found no genome-wide or genomic context-specific DNA methylation and hydroxymethylation changes with advancing age. These findings suggest that cellular senescence may contribute to ovarian aging, and thus, declines in ovarian follicular reserve. Cell-type-specific analyses across the reproductive lifespan are needed to fully elucidate the mechanisms that promote ovarian insufficiency.

Advancing Biophysics Using DNA Origami

DNA origami enables the bottom-up construction of chemically addressable, nanoscale objects with user-defined shapes and tailored functionalities. As such, not only can DNA origami objects be used to improve existing experimental methods in biophysics, but they also open up completely new avenues of exploration. In this review, we discuss basic biophysical concepts that are relevant for prospective DNA origami users. We summarize biochemical strategies for interfacing DNA origami with biomolecules of interest. We describe various applications of DNA origami, emphasizing the added value or new biophysical insights that can be generated: rulers and positioning devices, force measurement and force application devices, alignment supports for structural analysis for biomolecules in cryogenic electron microscopy and nuclear magnetic resonance, probes for manipulating and interacting with lipid membranes, and programmable nanopores. We conclude with some thoughts on so-far little explored opportunities for using DNA origami in more complex environments such as the cell or even organisms.

Genome-wide integration site detection using Cas9 enriched amplification-free long-range sequencing

The gene and cell therapy fields are advancing rapidly, with a potential to treat and cure a wide range of diseases, and lentivirus-based gene transfer agents are the vector of choice for many investigators. Early cases of insertional mutagenesis caused by gammaretroviral vectors highlighted that integration site (IS) analysis was a major safety and quality control checkpoint for lentiviral applications. The methods established to detect lentiviral integrations using next-generation sequencing (NGS) are limited by short read length, inadvertent PCR bias, low yield, or lengthy protocols. Here, we describe a new method to sequence IS using Amplification-free Integration Site sequencing (AFIS-Seq). AFIS-Seq is based on amplification-free, Cas9-mediated enrichment of high-molecular-weight chromosomal DNA suitable for long-range Nanopore MinION sequencing. This accessible and low-cost approach generates long reads enabling IS mapping with high certainty within a single day. We demonstrate proof-of-concept by mapping IS of lentiviral vectors in a variety of cell models and report up to 1600-fold enrichment of the signal. This method can be further extended to sequencing of Cas9-mediated integration of genes and to in vivo analysis of IS. AFIS-Seq uses long-read sequencing to facilitate safety evaluation of preclinical lentiviral vector gene therapies by providing IS analysis with improved confidence.

Megabase-scale methylation phasing using nanopore long reads and NanoMethPhase

The ability of nanopore sequencing to simultaneously detect modified nucleotides while producing long reads makes it ideal for detecting and phasing allele-specific methylation. However, there is currently no complete software for detecting SNPs, phasing haplotypes, and mapping methylation to these from nanopore sequence data. Here, we present NanoMethPhase, a software tool to phase 5-methylcytosine from nanopore sequencing. We also present SNVoter, which can post-process nanopore SNV calls to improve accuracy in low coverage regions. Together, these tools can accurately detect allele-specific methylation genome-wide using nanopore sequence data with low coverage of about ten-fold redundancy.

A meta-analysis uncovers the first sequence variant conferring risk of Bell's palsy

Bell's palsy is the most common cause of unilateral facial paralysis and is defined as an idiopathic and acute inability to control movements of the facial muscles on the affected side. While the pathogenesis remains unknown, previous studies have implicated post-viral inflammation and resulting compression of the facial nerve. Reported heritability estimates of 4-14% suggest a genetic component in the etiology and an autosomal dominant inheritance has been proposed. Here, we report findings from a meta-analysis of genome-wide association studies uncovering the first unequivocal association with Bell's palsy (rs9357446-A; P = 6.79 × 10^-23, OR = 1.23; N_cases = 4714, N_controls = 1,011,520). The variant also confers risk of intervertebral disc disorders (P = 2.99 × 10^-11, OR = 1.04) suggesting a common pathogenesis in part or a true pleiotropy.

Decoding the epitranscriptional landscape from native RNA sequences

Traditional epitranscriptomics relies on capturing a single RNA modification by antibody or chemical treatment, combined with short-read sequencing to identify its transcriptomic location. This approach is labor-intensive and may introduce experimental artifacts. Direct sequencing of native RNA using Oxford Nanopore Technologies (ONT) can allow for directly detecting the RNA base modifications, although these modifications might appear as sequencing errors. The percent Error of Specific Bases (%ESB) was higher for native RNA than unmodified RNA, which enabled the detection of ribonucleotide modification sites. Based on the %ESB differences, we developed a bioinformatic tool, epitranscriptional landscape inferring from glitches of ONT signals (ELIGOS), that is based on various types of synthetic modified RNA and applied to rRNA and mRNA. ELIGOS is able to accurately predict known classes of RNA methylation sites (AUC > 0.93) in rRNAs from Escherichiacoli, yeast, and human cells, using either unmodified in vitro transcription RNA or a background error model, which mimics the systematic error of direct RNA sequencing as the reference. The well-known DRACH/RRACH motif was localized and identified, consistent with previous studies, using differential analysis of ELIGOS to study the impact of RNA m6A methyltransferase by comparing wild type and knockouts in yeast and mouse cells. Lastly, the DRACH motif could also be identified in the mRNA of three human cell lines. The mRNA modification identified by ELIGOS is at the level of individual base resolution. In summary, we have developed a bioinformatic software package to uncover native RNA modifications.

Potential of whole-genome sequencing-based pharmacogenetic profiling

Pharmacogenetics represents a major driver of precision medicine, promising individualized drug selection and dosing. Traditionally, pharmacogenetic profiling has been performed using targeted genotyping that focuses on common/known variants. Recently, whole-genome sequencing (WGS) is emerging as a more comprehensive short-read next-generation sequencing approach, enabling both gene diagnostics and pharmacogenetic profiling, including rare/novel variants, in a single assay. Using the example of the pharmacogene CYP2D6, we demonstrate the potential of WGS-based pharmacogenetic profiling as well as emphasize the limitations of short-read next-generation sequencing. In the near future, we envision a shift toward long-read sequencing as the predominant method for gene diagnostics and pharmacogenetic profiling, providing unprecedented data quality and improving patient care.

Nondestructive enzymatic deamination enables single-molecule long-read amplicon sequencing for the determination of 5-methylcytosine and 5-hydroxymethylcytosine at single-base resolution

The predominant methodology for DNA methylation analysis relies on the chemical deamination by sodium bisulfite of unmodified cytosine to uracil to permit the differential readout of methylated cytosines. Bisulfite treatment damages the DNA, leading to fragmentation and loss of long-range methylation information. To overcome this limitation of bisulfite-treated DNA, we applied a new enzymatic deamination approach, termed enzymatic methyl-seq (EM-seq), to long-range sequencing technologies. Our methodology, named long-read enzymatic modification sequencing (LR-EM-seq), preserves the integrity of DNA, allowing long-range methylation profiling of 5-methylcytosine (5mC) and 5-hydroxymethylcytosine (5hmC) over multikilobase length of genomic DNA. When applied to known differentially methylated regions (DMRs), LR-EM-seq achieves phasing of >5 kb, resulting in broader and better defined DMRs compared with that previously reported. This result showed the importance of phasing methylation for biologically relevant questions and the applicability of LR-EM-seq for long-range epigenetic analysis at single-molecule and single-nucleotide resolution.

Decoding the epitranscriptional landscape from native RNA sequences

Traditional epitranscriptomics relies on capturing a single RNA modification by antibody or chemical treatment, combined with short-read sequencing to identify its transcriptomic location. This approach is labor-intensive and may introduce experimental artifacts. Direct sequencing of native RNA using Oxford Nanopore Technologies (ONT) can allow for directly detecting the RNA base modifications, although these modifications might appear as sequencing errors. The percent Error of Specific Bases (%ESB) was higher for native RNA than unmodified RNA, which enabled the detection of ribonucleotide modification sites. Based on the %ESB differences, we developed a bioinformatic tool, epitranscriptional landscape inferring from glitches of ONT signals (ELIGOS), that is based on various types of synthetic modified RNA and applied to rRNA and mRNA. ELIGOS is able to accurately predict known classes of RNA methylation sites (AUC > 0.93) in rRNAs from Escherichiacoli, yeast, and human cells, using either unmodified in vitro transcription RNA or a background error model, which mimics the systematic error of direct RNA sequencing as the reference. The well-known DRACH/RRACH motif was localized and identified, consistent with previous studies, using differential analysis of ELIGOS to study the impact of RNA m6A methyltransferase by comparing wild type and knockouts in yeast and mouse cells. Lastly, the DRACH motif could also be identified in the mRNA of three human cell lines. The mRNA modification identified by ELIGOS is at the level of individual base resolution. In summary, we have developed a bioinformatic software package to uncover native RNA modifications.

MinION Nanopore-based detection of Clavibacter nebraskensis, the corn Goss's wilt pathogen, and bacteriomic profiling of necrotic lesions of naturally-infected leaf samples

The Goss’s bacterial wilt pathogen, Clavibacter nebraskensis, of corn is a candidate A1 quarantine organism; and its recent re-emergence and spread in the USA and Canada is a potential biothreat to the crop. We developed and tested an amplicon-based Nanopore detection system for C. nebraskensis (Cn), targeting a purine permease gene. The sensitivity (1 pg) of this system in mock bacterial communities (MBCs) spiked with serially diluted DNA of C. nebraskensis NCPPB 2581^T is comparable to that of real-time PCR. Average Nanopore reads increased exponentially from 125 (1pg) to about 6000 reads (1000 pg) after a 3-hr run-time, with 99.0% of the reads accurately assigned to C. nebraskensis. Three run-times were used to process control MBCs, Cn-spiked MBCs, diseased and healthy leaf samples. The mean Nanopore reads doubled as the run-time is increased from 3 to 6 hrs while from 6 to 12 hrs, a 20% increment was recorded in all treatments. Cn-spiked MBCs and diseased corn leaf samples averaged read counts of 5,100, 11,000 and 14,000 for the respective run-times, with 99.8% of the reads taxonomically identified as C. nebraskensis. The control MBCs and healthy leaf samples had 47 and 14 Nanopore reads, respectively. 16S rRNA bacteriomic profiles showed that Sphingomonas (22.7%) and Clavibacter (21.2%) were dominant in diseased samples while Pseudomonas had only 3.5% relative abundance. In non-symptomatic leaf samples, however, Pseudomonas (20.0%) was dominant with Clavibacter at 0.08% relative abundance. This discrepancy in Pseudomonas abundance in the samples was corroborated by qPCR using EvaGreen chemistry. Our work outlines a new useful tool for diagnosis of the Goss’s bacterial wilt disease; and provides the first insight on Pseudomonas community dynamics in necrotic leaf lesions.

Adaptation of Human Ribosomal RNA for Nanopore Sequencing of Canonical and Modified Nucleotides

Historically, RNA has been sequenced as cDNA copies derived from reverse transcription of cellular RNA followed by PCR amplification. Recently, RNA sequencing using nanopores has emerged as an alternative. Using this technology, individual cellular RNA strands are read directly as they are driven through nanoscale pores by an applied voltage. The speed of translocation is regulated by a helicase that is loaded onto each RNA strand by an adapter that also facilitates capture by the nanopore electric field. Here we describe a technique for adapting human ribosomal RNA subunits for nanopore sequencing. Using this strategy, a single Oxford Nanopore MinION run delivered 470,907 sequence reads of which 396,048 aligned to ribosomal RNA, with 28S, 18S, 5.8S, and 5S coverage of 6053, 369,472, 16,058, and 4465 reads, respectively. Example alignments that reveal putative nucleotide modifications are provided.

Liquid Chromatography-Mass Spectrometry Analysis of Cytosine Modifications

Liquid chromatography-tandem mass spectrometry (LC-MS/MS) is a widely used technique in the global analysis of epigenetic DNA modifications. The high-resolution chromatographic separation along with sensitive MS detection permits the identification and quantification of deoxyribonucleosides with precision and reliability. Although there have been tremendous advances in LC and MS instrumentation in recent years, sample preparation has not experienced a similar rate of development and is often a bottleneck to chemical analysis. Here we present a protocol for identification and quantification of cytosine modifications that combines a robust and efficient method to generate single nucleosides from genomic DNA samples followed by direct LC-MS/MS analysis.

Mapping DNA Methylation in Mammals: The State of the Art

A complete understanding of the dynamics and function of cytosine modifications in mammalian biology is lacking. Central to achieving this understanding is the availability of techniques that permit sensitive and specific genome-wide mapping of DNA modifications in mammalian DNA. The last decade has seen the development of a vast arsenal of novel profiling approaches enabling epigeneticists to tackle research questions that were previously out of reach. Here, we review the techniques currently available for profiling DNA modifications in mammals, discuss their strengths and weaknesses, and speculate on the future direction of DNA modification profiling technologies.

Epigenetics of Hepatic Insulin Resistance

Insulin resistance (IR) is largely recognized as a unifying feature that underlies metabolic dysfunction. Both lifestyle and genetic factors contribute to IR. Work from recent years has demonstrated that the epigenome may constitute an interface where different signals may converge to promote IR gene expression programs. Here, we review the current knowledge of the role of epigenetics in hepatic IR, focusing on the roles of DNA methylation and histone post-translational modifications. We discuss the broad epigenetic changes observed in the insulin resistant liver and its associated pathophysiological states and leverage on the wealth of 'omics' studies performed to discuss efforts in pinpointing specific loci that are disrupted by these changes. We envision that future studies, with increased genomic resolution and larger cohorts, will further the identification of biomarkers of early onset hepatic IR and assist the development of targeted interventions. Furthermore, there is growing evidence to suggest that persistent epigenetic marks may be acquired over prolonged exposure to disease or deleterious exposures, highlighting the need for preventative medicine and long-term lifestyle adjustments to avoid irreversible or long-term alterations in gene expression.

The structure, function and evolution of a complete human chromosome 8

The complete assembly of each human chromosome is essential for understanding human biology and evolution^1,2. Here we use complementary long-read sequencing technologies to complete the linear assembly of human chromosome 8. Our assembly resolves the sequence of five previously long-standing gaps, including a 2.08-Mb centromeric α-satellite array, a 644-kb copy number polymorphism in the β-defensin gene cluster that is important for disease risk, and an 863-kb variable number tandem repeat at chromosome 8q21.2 that can function as a neocentromere. We show that the centromeric α-satellite array is generally methylated except for a 73-kb hypomethylated region of diverse higher-order α-satellites enriched with CENP-A nucleosomes, consistent with the location of the kinetochore. In addition, we confirm the overall organization and methylation pattern of the centromere in a diploid human genome. Using a dual long-read sequencing approach, we complete high-quality draft assemblies of the orthologous centromere from chromosome 8 in chimpanzee, orangutan and macaque to reconstruct its evolutionary history. Comparative and phylogenetic analyses show that the higher-order α-satellite structure evolved in the great ape ancestor with a layered symmetry, in which more ancient higher-order repeats locate peripherally to monomeric α-satellites. We estimate that the mutation rate of centromeric satellite DNA is accelerated by more than 2.2-fold compared to the unique portions of the genome, and this acceleration extends into the flanking sequence.

New Twists in Detecting mRNA Modification Dynamics

Modified nucleotides in mRNA are an essential addition to the standard genetic code of four nucleotides in animals, plants, and their viruses. The emerging field of epitranscriptomics examines nucleotide modifications in mRNA and their impact on gene expression. The low abundance of nucleotide modifications and technical limitations, however, have hampered systematic analysis of their occurrence and functions. Selective chemical and immunological identification of modified nucleotides has revealed global candidate topology maps for many modifications in mRNA, but further technical advances to increase confidence will be necessary. Single-molecule sequencing introduced by Oxford Nanopore now promises to overcome such limitations, and we summarize current progress with a particular focus on the bioinformatic challenges of this novel sequencing technology.

Intraoperative DNA methylation classification of brain tumors impacts neurosurgical strategy

BACKGROUND

Brain tumor surgery must balance the benefit of maximal resection against the risk of inflicting severe damage. The impact of increased resection is diagnosis-specific. However, the precise diagnosis is typically uncertain at surgery due to limitations of imaging and intraoperative histomorphological methods. Novel and accurate strategies for brain tumor classification are necessary to support personalized intraoperative neurosurgical treatment decisions. Here, we describe a fast and cost-efficient workflow for intraoperative classification of brain tumors based on DNA methylation profiles generated by low coverage nanopore sequencing and machine learning algorithms.

METHODS

We evaluated 6 independent cohorts containing 105 patients, including 50 pediatric and 55 adult patients. Ultra-low coverage whole-genome sequencing was performed on nanopore flow cells. Data were analyzed using copy number variation and ad hoc random forest classifier for the genome-wide methylation-based classification of the tumor.

RESULTS

Concordant classification was obtained between nanopore DNA methylation analysis and a full neuropathological evaluation in 93 of 105 (89%) cases. The analysis demonstrated correct diagnosis in 6/6 cases where frozen section evaluation was inconclusive. Results could be returned to the operating room at a median of 97 min (range 91-161 min). Precise classification of the tumor entity and subtype would have supported modification of the surgical strategy in 12 out of 20 patients evaluated intraoperatively.

CONCLUSION

Intraoperative nanopore sequencing combined with machine learning diagnostics was robust, sensitive, and rapid. This strategy allowed DNA methylation-based classification of the tumor to be returned to the surgeon within a timeframe that supports intraoperative decision making.

Targeted methods for epigenetic age predictions in mice

Age-associated DNA methylation reflects aspect of biological aging—therefore epigenetic clocks for mice can elucidate how the aging process in this model organism is affected by specific treatments or genetic background. Initially, age-predictors for mice were trained for genome-wide DNA methylation profiles and we have recently described a targeted assay based on pyrosequencing of DNA methylation at only three age-associated genomic regions. Here, we established alternative approaches using droplet digital PCR (ddPCR) and barcoded bisulfite amplicon sequencing (BBA-seq). At individual CG dinucleotides (CpGs) the correlation of DNA methylation with chronological age was slightly higher for pyrosequencing and ddPCR as compared to BBA-seq. On the other hand, BBA-seq revealed that neighboring CpGs tend to be stochastically modified at murine age-associated regions. Furthermore, the binary sequel of methylated and non-methylated CpGs in individual reads can be used for single-read predictions, which may reflect heterogeneity in epigenetic aging. In comparison to C57BL/6 mice the single-read age-predictions using BBA-seq were also accelerated in the shorter-lived DBA/2 mice, and in C57BL/6 mice with a lifespan quantitative trait locus of DBA/2 mice. Taken together, we describe alternative targeted methods for epigenetic age predictions that provide new perspectives for aging-intervention studies in mice.

A comprehensive evaluation of long read error correction methods

BACKGROUND

Third-generation single molecule sequencing technologies can sequence long reads, which is advancing the frontiers of genomics research. However, their high error rates prohibit accurate and efficient downstream analysis. This difficulty has motivated the development of many long read error correction tools, which tackle this problem through sampling redundancy and/or leveraging accurate short reads of the same biological samples. Existing studies to asses these tools use simulated data sets, and are not sufficiently comprehensive in the range of software covered or diversity of evaluation measures used.

RESULTS

In this paper, we present a categorization and review of long read error correction methods, and provide a comprehensive evaluation of the corresponding long read error correction tools. Leveraging recent real sequencing data, we establish benchmark data sets and set up evaluation criteria for a comparative assessment which includes quality of error correction as well as run-time and memory usage. We study how trimming and long read sequencing depth affect error correction in terms of length distribution and genome coverage post-correction, and the impact of error correction performance on an important application of long reads, genome assembly. We provide guidelines for practitioners for choosing among the available error correction tools and identify directions for future research.

CONCLUSIONS

Despite the high error rate of long reads, the state-of-the-art correction tools can achieve high correction quality. When short reads are available, the best hybrid methods outperform non-hybrid methods in terms of correction quality and computing resource usage. When choosing tools for use, practitioners are suggested to be careful with a few correction tools that discard reads, and check the effect of error correction tools on downstream analysis. Our evaluation code is available as open-source at https://github.com/haowenz/LRECE .

Simulation of Nanopore Sequencing Signals Based on BiGRU

Oxford Nanopore sequencing is an important sequencing technology, which reads the nucleotide sequence by detecting the electrical current signal changes when DNA molecule is forced to pass through a biological nanopore. The research on signal simulation of nanopore sequencing is highly desirable for method developments of nanopore sequencing applications. To improve the simulation accuracy, we propose a novel signal simulation method based on Bi-directional Gated Recurrent Units (BiGRU). In this method, the signal processing model based on BiGRU is built to replace the traditional low-pass filter to post-process the ground-truth signal calculated by the input nucleotide sequence and nanopore sequencing pore model. Gaussian noise is then added to the filtered signal to generate the final simulated signal. This method can accurately model the relation between ground-truth signal and real-world sequencing signal through experimental sequencing data. The simulation results reveal that the proposed method utilizing the powerful learning ability of the neural network can generate the simulated signal that is closer to the real-world sequencing signal in the time and frequency domains than the existing simulation method.

Impact of lossy compression of nanopore raw signal data on basecalling and consensus accuracy

Motivation

Nanopore sequencing provides a real-time and portable solution to genomic sequencing, enabling better assembly, structural variant discovery and modified base detection than second generation technologies. The sequencing process generates a huge amount of data in the form of raw signal contained in fast5 files, which must be compressed to enable efficient storage and transfer. Since the raw data is inherently noisy, lossy compression has potential to significantly reduce space requirements without adversely impacting performance of downstream applications.

Results

We explore the use of lossy compression for nanopore raw data using two state-of-the-art lossy time-series compressors, and evaluate the tradeoff between compressed size and basecalling/consensus accuracy. We test several basecallers and consensus tools on a variety of datasets at varying depths of coverage, and conclude that lossy compression can provide 35–50% further reduction in compressed size of raw data over the state-of-the-art lossless compressor with negligible impact on basecalling accuracy (⁠≲0.2% reduction) and consensus accuracy (⁠≲0.002% reduction). In addition, we evaluate the impact of lossy compression on methylation calling accuracy and observe that this impact is minimal for similar reductions in compressed size, although further evaluation with improved benchmark datasets is required for reaching a definite conclusion. The results suggest the possibility of using lossy compression, potentially on the nanopore sequencing device itself, to achieve significant reductions in storage and transmission costs while preserving the accuracy of downstream applications.

Availabilityand implementation

The code is available at https://github.com/shubhamchandak94/lossy_compression_evaluation.

Analytical validity of nanopore sequencing for rapid SARS-CoV-2 genome analysis

Viral whole-genome sequencing (WGS) provides critical insight into the transmission and evolution of Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2). Long-read sequencing devices from Oxford Nanopore Technologies (ONT) promise significant improvements in turnaround time, portability and cost, compared to established short-read sequencing platforms for viral WGS (e.g., Illumina). However, adoption of ONT sequencing for SARS-CoV-2 surveillance has been limited due to common concerns around sequencing accuracy. To address this, here we perform viral WGS with ONT and Illumina platforms on 157 matched SARS-CoV-2-positive patient specimens and synthetic RNA controls, enabling rigorous evaluation of analytical performance. We report that, despite the elevated error rates observed in ONT sequencing reads, highly accurate consensus-level sequence determination was achieved, with single nucleotide variants (SNVs) detected at >99% sensitivity and >99% precision above a minimum ~60-fold coverage depth, thereby ensuring suitability for SARS-CoV-2 genome analysis. ONT sequencing also identified a surprising diversity of structural variation within SARS-CoV-2 specimens that were supported by evidence from short-read sequencing on matched samples. However, ONT sequencing failed to accurately detect short indels and variants at low read-count frequencies. This systematic evaluation of analytical performance for SARS-CoV-2 WGS will facilitate widespread adoption of ONT sequencing within local, national and international COVID-19 public health initiatives.

Simultaneous profiling of chromatin accessibility and methylation on human cell lines with nanopore sequencing

Probing epigenetic features on DNA has tremendous potential to advance our understanding of the phased epigenome. In this study, we use nanopore sequencing to evaluate CpG methylation and chromatin accessibility simultaneously on long strands of DNA by applying GpC methyltransferase to exogenously label open chromatin. We performed nanopore sequencing of nucleosome occupancy and methylome (nanoNOMe) on four human cell lines (GM12878, MCF-10A, MCF-7 and MDA-MB-231). The single-molecule resolution allows footprinting of protein and nucleosome binding, and determination of the combinatorial promoter epigenetic signature on individual molecules. Long-read sequencing makes it possible to robustly assign reads to haplotypes, allowing us to generate a fully phased human epigenome, consisting of chromosome-level allele-specific profiles of CpG methylation and chromatin accessibility. We further apply this to a breast cancer model to evaluate differential methylation and accessibility between cancerous and noncancerous cells.

Inducible cell-specific mouse models for paired epigenetic and transcriptomic studies of microglia and astroglia

Epigenetic regulation of gene expression occurs in a cell type-specific manner. Current cell-type specific neuroepigenetic studies rely on cell sorting methods that can alter cell phenotype and introduce potential confounds. Here we demonstrate and validate a Nuclear Tagging and Translating Ribosome Affinity Purification (NuTRAP) approach for temporally controlled labeling and isolation of ribosomes and nuclei, and thus RNA and DNA, from specific central nervous system cell types. Analysis of gene expression and DNA modifications in astrocytes or microglia from the same animal demonstrates differential usage of DNA methylation and hydroxymethylation in CpG and non-CpG contexts that corresponds to cell type-specific gene expression. Application of this approach in LPS treated mice uncovers microglia-specific transcriptome and epigenome changes in inflammatory pathways that cannot be detected with tissue-level analysis. The NuTRAP model and the validation approaches presented can be applied to any brain cell type for which a cell type-specific cre is available.

Ana Chucair-Elliott, Sarah Ocaňas et al. present a NuTRAP approach for simultaneous analysis of transcript expression and DNA modifications in two specific mouse brain cell types, astrocytes and microglia. They further apply this approach to identify molecular changes in microglia following LPS treatment and identify both transcriptomic and epigenomic alterations not apparent in tissue-level analyses.

A hybrid correcting method considering heterozygous variations by a comprehensive probabilistic model

Background

The emergence of the third generation sequencing technology, featuring longer read lengths, has demonstrated great advancement compared to the next generation sequencing technology and greatly promoted the biological research. However, the third generation sequencing data has a high level of the sequencing error rates, which inevitably affects the downstream analysis. Although the issue of sequencing error has been improving these years, large amounts of data were produced at high sequencing errors, and huge waste will be caused if they are discarded. Thus, the error correction for the third generation sequencing data is especially important. The existing error correction methods have poor performances at heterozygous sites, which are ubiquitous in diploid and polyploidy organisms. Therefore, it is a lack of error correction algorithms for the heterozygous loci, especially at low coverages.

Results

In this article, we propose a error correction method, named QIHC. QIHC is a hybrid correction method, which needs both the next generation and third generation sequencing data. QIHC greatly enhances the sensitivity of identifying the heterozygous sites from sequencing errors, which leads to a high accuracy on error correction. To achieve this, QIHC established a set of probabilistic models based on Bayesian classifier, to estimate the heterozygosity of a site and makes a judgment by calculating the posterior probabilities. The proposed method is consisted of three modules, which respectively generates a pseudo reference sequence, obtains the read alignments, estimates the heterozygosity the sites and corrects the read harboring them. The last module is the core module of QIHC, which is designed to fit for the calculations of multiple cases at a heterozygous site. The other two modules enable the reads mapping to the pseudo reference sequence which somehow overcomes the inefficiency of multiple mappings that adopt by the existing error correction methods.

Conclusions

To verify the performance of our method, we selected Canu and Jabba to compare with QIHC in several aspects. As a hybrid correction method, we first conducted a groups of experiments under different coverages of the next-generation sequencing data. QIHC is far ahead of Jabba on accuracy. Meanwhile, we varied the coverages of the third generation sequencing data and compared performances again among Canu, Jabba and QIHC. QIHC outperforms the other two methods on accuracy of both correcting the sequencing errors and identifying the heterozygous sites, especially at low coverage. We carried out a comparison analysis between Canu and QIHC on the different error rates of the third generation sequencing data. QIHC still performs better. Therefore, QIHC is superior to the existing error correction methods when heterozygous sites exist.

Enhancer hijacking determines extrachromosomal circular MYCN amplicon architecture in neuroblastoma

MYCN amplification drives one in six cases of neuroblastoma. The supernumerary gene copies are commonly found on highly rearranged, extrachromosomal circular DNA (ecDNA). The exact amplicon structure has not been described thus far and the functional relevance of its rearrangements is unknown. Here, we analyze the MYCN amplicon structure using short-read and Nanopore sequencing and its chromatin landscape using ChIP-seq, ATAC-seq and Hi-C. This reveals two distinct classes of amplicons which explain the regulatory requirements for MYCN overexpression. The first class always co-amplifies a proximal enhancer driven by the noradrenergic core regulatory circuit (CRC). The second class of MYCN amplicons is characterized by high structural complexity, lacks key local enhancers, and instead contains distal chromosomal fragments harboring CRC-driven enhancers. Thus, ectopic enhancer hijacking can compensate for the loss of local gene regulatory elements and explains a large component of the structural diversity observed in MYCN amplification.

MYCN amplification is common in neuroblastomas. Here the authors analyse the MYCN amplicon structure and its epigenetic regulation by integrating short- and longread genomic and epigenomic data and find two classes of MYCN amplicons in neuroblastomas, one driven by local enhancers and the other by hijacking of distal regulatory elements.

Nanopore Sequencing Enables Comprehensive Transposable Element Epigenomic Profiling

Transposable elements (TEs) drive genome evolution and are a notable source of pathogenesis, including cancer. While CpG methylation regulates TE activity, the locus-specific methylation landscape of mobile human TEs has to date proven largely inaccessible. Here, we apply new computational tools and long-read nanopore sequencing to directly infer CpG methylation of novel and extant TE insertions in hippocampus, heart, and liver, as well as paired tumor and non-tumor liver. As opposed to an indiscriminate stochastic process, we find pronounced demethylation of young long interspersed element 1 (LINE-1) retrotransposons in cancer, often distinct to the adjacent genome and other TEs. SINE-VNTR-Alu (SVA) retrotransposons, including their internal tandem repeat-associated CpG island, are near-universally methylated. We encounter allele-specific TE methylation and demethylation of aberrantly expressed young LINE-1s in normal tissues. Finally, we recover the complete sequences of tumor-specific LINE-1 insertions and their retrotransposition hallmarks, demonstrating how long-read sequencing can simultaneously survey the epigenome and detect somatic TE mobilization.

High quality genome assemblies of Mycoplasma bovis using a taxon-specific Bonito basecaller for MinION and Flongle long-read nanopore sequencing

BACKGROUND

Implementation of Third-Generation Sequencing approaches for Whole Genome Sequencing (WGS) all-in-one diagnostics in human and veterinary medicine, requires the rapid and accurate generation of consensus genomes. Over the last years, Oxford Nanopore Technologies (ONT) released various new devices (e.g. the Flongle R9.4.1 flow cell) and bioinformatics tools (e.g. the in 2019-released Bonito basecaller), allowing cheap and user-friendly cost-efficient introduction in various NGS workflows. While single read, overall consensus accuracies, and completeness of genome sequences has been improved dramatically, further improvements are required when working with non-frequently sequenced organisms like Mycoplasma bovis. As an important primary respiratory pathogen in cattle, rapid M. bovis diagnostics is crucial to allow timely and targeted disease control and prevention. Current complete diagnostics (including identification, strain typing, and antimicrobial resistance (AMR) detection) require combined culture-based and molecular approaches, of which the first can take 1-2 weeks. At present, cheap and quick long read all-in-one WGS approaches can only be implemented if increased accuracies and genome completeness can be obtained.

RESULTS

Here, a taxon-specific custom-trained Bonito v.0.1.3 basecalling model (custom-pg45) was implemented in various WGS assembly bioinformatics pipelines. Using MinION sequencing data, we showed improved consensus accuracies up to Q45.2 and Q46.7 for reference-based and Canu de novo assembled M. bovis genomes, respectively. Furthermore, the custom-pg45 model resulted in mean consensus accuracies of Q45.0 and genome completeness of 94.6% for nine M. bovis field strains. Improvements were also observed for the single-use Flongle sequencer (mean Q36.0 accuracies and 80.3% genome completeness).

CONCLUSIONS

These results implicate that taxon-specific basecalling of MinION and single-use Flongle Nanopore long reads are of great value to be implemented in rapid all-in-one WGS tools as evidenced for Mycoplasma bovis as an example.

Sensitive alignment using paralogous sequence variants improves long-read mapping and variant calling in segmental duplications

The ability to characterize repetitive regions of the human genome is limited by the read lengths of short-read sequencing technologies. Although long-read sequencing technologies such as Pacific Biosciences (PacBio) and Oxford Nanopore Technologies can potentially overcome this limitation, long segmental duplications with high sequence identity pose challenges for long-read mapping. We describe a probabilistic method, DuploMap, designed to improve the accuracy of long-read mapping in segmental duplications. It analyzes reads mapped to segmental duplications using existing long-read aligners and leverages paralogous sequence variants (PSVs)—sequence differences between paralogous sequences—to distinguish between multiple alignment locations. On simulated datasets, DuploMap increased the percentage of correctly mapped reads with high confidence for multiple long-read aligners including Minimap2 (74.3–90.6%) and BLASR (82.9–90.7%) while maintaining high precision. Across multiple whole-genome long-read datasets, DuploMap aligned an additional 8–21% of the reads in segmental duplications with high confidence relative to Minimap2. Using DuploMap-aligned PacBio circular consensus sequencing reads, an additional 8.9 Mb of DNA sequence was mappable, variant calling achieved a higher F₁ score and 14 713 additional variants supported by linked-read data were identified. Finally, we demonstrate that a significant fraction of PSVs in segmental duplications overlaps with variants and adversely impacts short-read variant calling.