Detecting DNA cytosine methylation using nanopore sequencing

Stress-induced epigenetic effects driven by maternal lactation in dairy cattle: a comethylation network approach

Epigenetic marks do not follow the Mendelian laws of inheritance. The environment can alter the epigenotype of an individual when exposed to different external stressors. In lactating cows, the first stages of gestation overlap with the lactation peak, creating a negative energy balance that is difficult to overcome with diet. This negative energy balance could affect early embryo development that must compete with the mammary tissue for nutrients. We hypothesize that the methylation profiles of calves born to nonlactating heifers are different from those of calves born to lactating cows. We found 50,277 differentially methylated cytosines and 2,281 differentially methylated regions between these two groups of animals. A comethylation network was constructed to study the correlation between the phenotypes of the mothers and the epigenome of the calves, revealing 265 regions associated with the phenotypes. Our study revealed the presence of DMCs and DMRs in calves gestated by heifers and lactating cows, which were linked to the dam's lactation and the calves' ICAP and milk EBV. Gene-specific analysis highlighted associations with vasculature and organ morphogenesis and cell communication and signalling. These finding support the hypothesis that calves gestated by nonlactating mothers have a different methylation profile than those gestated by lactating cows.

Nanopore signal deviations from pseudouridine modifications in RNA are sequence-specific: quantification requires dedicated synthetic controls

Chemical modifications to mRNA respond dynamically to environmental cues and are important modulators of gene expression. Nanopore direct RNA sequencing has been applied for assessing the presence of pseudouridine (ψ) modifications through basecalling errors and signal analysis. These approaches strongly depend on the sequence context around the modification, and the occupancies derived from these measurements are not quantitative. In this work, we combine direct RNA sequencing of synthetic RNAs bearing site-specific modifications and supervised machine learning models (ModQuant) to achieve near-analytical, site-specific ψ quantification. Our models demonstrate that the ionic current signal features important for accurate ψ classification are sequence dependent and encompass information extending beyond n + 2 and n - 2 nucleotides from the ψ site. This is contradictory to current models, which assume that accurate ψ classification can be achieved with signal information confined to the 5-nucleotide k-mer window (n + 2 and n - 2 nucleotides from the ψ site). We applied our models to quantitatively profile ψ occupancy in five mRNA sites in datasets from seven human cell lines, demonstrating conserved and variable sites. Our study motivates a wider pipeline that uses ground-truth RNA control sets with site-specific modifications for quantitative profiling of RNA modifications. The ModQuant pipeline and guide are freely available at https://github.com/wanunulab/ModQuant .

A phenome-wide association study of methylated GC-rich repeats identifies a GCC repeat expansion in AFF3 associated with intellectual disability

GC-rich tandem repeat expansions (TREs) are often associated with DNA methylation, gene silencing and folate-sensitive fragile sites, and underlie several congenital and late-onset disorders. Through a combination of DNA-methylation profiling and tandem repeat genotyping, we identified 24 methylated TREs and investigated their effects on human traits using phenome-wide association studies in 168,641 individuals from the UK Biobank, identifying 156 significant TRE-trait associations involving 17 different TREs. Of these, a GCC expansion in the promoter of AFF3 was associated with a 2.4-fold reduced probability of completing secondary education, an effect size comparable to several recurrent pathogenic microdeletions. In a cohort of 6,371 probands with neurodevelopmental problems of suspected genetic etiology, we observed a significant enrichment of AFF3 expansions compared with controls. With a population prevalence that is at least fivefold higher than the TRE that causes fragile X syndrome, AFF3 expansions represent a major cause of neurodevelopmental delay.

Visualizing epigenetic modifications and their spatial proximities in single cells using three DNA-encoded amplifying FISH imaging strategies: BEA-FISH, PPDA-FISH and Cell-TALKING

Epigenetic modifications and spatial proximities of nucleic acids and proteins play important roles in regulating physiological processes and disease progression. Currently available cell imaging methods, such as fluorescence in situ hybridization (FISH) and immunofluorescence, struggle to detect low-abundance modifications and their spatial proximities. Here we describe a step-by-step protocol for three DNA-encoded amplifying FISH-based imaging strategies to overcome these challenges for varying applications: base-encoded amplifying FISH (BEA-FISH), pairwise proximity-differentiated amplifying FISH (PPDA-FISH) and cellular macromolecules-tethered DNA walking indexing (Cell-TALKING). They all use the similar core principle of DNA-encoded amplification, which transforms different nonsequence molecular features into unique DNA barcodes for in situ rolling circle amplification and FISH analysis. This involves three key reactions in fixed cell samples: target labeling, DNA encoding and rolling circle amplification imaging. Using this protocol, these three imaging strategies achieve in situ counting of low-abundance modifications alone, the pairwise proximity-differentiated visualization of two modifications and the exploration of multiple modifications around one protein (one-to-many proximity), respectively. Low-abundance modifications, including 5-hydroxymethylcytosine, 5-formylcytosine, 5-hydroxymethyluracil and 5-formyluracil, are clearly visualized in single cells. Various combinatorial patterns of nucleic acid modifications and/or histone modifications are found. The whole protocol takes ~2-4 d to complete, depending on different imaging applications.

Genome assembly in the telomere-to-telomere era

Genome sequences largely determine the biology and encode the history of an organism, and de novo assembly - the process of reconstructing the genome sequence of an organism from sequencing reads - has been a central problem in bioinformatics for four decades. Until recently, genomes were typically assembled into fragments of a few megabases at best, but now technological advances in long-read sequencing enable the near-complete assembly of each chromosome - also known as telomere-to-telomere assembly - for many organisms. Here, we review recent progress on assembly algorithms and protocols, with a focus on how to derive near-telomere-to-telomere assemblies. We also discuss the additional developments that will be required to resolve remaining assembly gaps and to assemble non-diploid genomes.

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

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

Targeting transposable elements in cancer: developments and opportunities

Transposable elements (TEs), comprising nearly 50% of the human genome, have transitioned from being perceived as "genomic junk" to key players in cancer progression. Contemporary research links TE regulatory disruptions with cancer development, underscoring their therapeutic potential. Advances in long-read sequencing, computational analytics, single-cell sequencing, proteomics, and CRISPR-Cas9 technologies have enriched our understanding of TEs' clinical implications, notably their impact on genome architecture, gene regulation, and evolutionary processes. In cancer, TEs, including long interspersed element-1 (LINE-1), Alus, and long terminal repeat (LTR) elements, demonstrate altered patterns, influencing both tumorigenic and tumor-suppressive mechanisms. TE-derived nucleic acids and tumor antigens play critical roles in tumor immunity, bridging innate and adaptive responses. Given their central role in oncology, TE-targeted therapies, particularly through reverse transcriptase inhibitors and epigenetic modulators, represent a novel avenue in cancer treatment. Combining these TE-focused strategies with existing chemotherapy or immunotherapy regimens could enhance efficacy and offer a new dimension in cancer treatment. This review delves into recent TE detection advancements, explores their multifaceted roles in tumorigenesis and immune regulation, discusses emerging diagnostic and therapeutic approaches centered on TEs, and anticipates future directions in cancer research.

DNA methylation governs the sensitivity of repeats to restriction by the HUSH-MORC2 corepressor

The human silencing hub (HUSH) complex binds to transcripts of LINE-1 retrotransposons (L1s) and other genomic repeats, recruiting MORC2 and other effectors to remodel chromatin. How HUSH and MORC2 operate alongside DNA methylation, a central epigenetic regulator of repeat transcription, remains largely unknown. Here we interrogate this relationship in human neural progenitor cells (hNPCs), a somatic model of brain development that tolerates removal of DNA methyltransferase DNMT1. Upon loss of MORC2 or HUSH subunit TASOR in hNPCs, L1s remain silenced by robust promoter methylation. However, genome demethylation and activation of evolutionarily-young L1s attracts MORC2 binding, and simultaneous depletion of DNMT1 and MORC2 causes massive accumulation of L1 transcripts. We identify the same mechanistic hierarchy at pericentromeric α-satellites and clustered protocadherin genes, repetitive elements important for chromosome structure and neurodevelopment respectively. Our data delineate the epigenetic control of repeats in somatic cells, with implications for understanding the vital functions of HUSH-MORC2 in hypomethylated contexts throughout human development.

An intronic copy number variation in Syntaxin 17 determines speed of greying and melanoma incidence in Grey horses

The Greying with age phenotype in horses involves loss of hair pigmentation whereas skin pigmentation is not reduced, and a predisposition to melanoma. The causal mutation was initially reported as a duplication of a 4.6 kb intronic sequence in Syntaxin 17. The speed of greying varies considerably among Grey horses. Here we demonstrate the presence of two different Grey alleles, G2 carrying two tandem copies of the duplicated sequence and G3 carrying three. The latter is by far the most common allele, probably due to strong selection for the striking white phenotype. Our results reveal a remarkable dosage effect where the G3 allele is associated with fast greying and high incidence of melanoma whereas G2 is associated with slow greying and low incidence of melanoma. The copy number expansion transforms a weak enhancer to a strong melanocyte-specific enhancer that underlies hair greying (G2 and G3) and a drastically elevated risk of melanoma (G3 only). Our direct pedigree-based observation of the origin of a G2 allele from a G3 allele by copy number contraction demonstrates the dynamic evolution of this locus and provides the ultimate evidence for causality of the copy number variation of the 4.6 kb intronic sequence.

Deciphering DNA Methylation in Gestational Diabetes Mellitus: Epigenetic Regulation and Potential Clinical Applications

Gestational diabetes mellitus (GDM) represents a prevalent complication during pregnancy, exerting both short-term and long-term impacts on maternal and offspring health. This review offers a comprehensive outline of DNA methylation modifications observed in various maternal and offspring tissues affected by GDM, emphasizing the intricate interplay between DNA methylation dynamics, gene expression, and the pathogenesis of GDM. Furthermore, it explores the influence of environmental pollutants, maternal nutritional supplementation, and prenatal gut microbiota on GDM development through alterations in DNA methylation profiles. Additionally, this review summarizes recent advancements in DNA methylation-based diagnostics and predictive models in early GDM detection and risk assessment for subsequent type 2 diabetes. These insights contribute significantly to our understanding of the epigenetic mechanisms underlying GDM development, thereby enhancing maternal and fetal health outcomes and advocating further efforts in this field.

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

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

GCphase: an SNP phasing method using a graph partition and error correction algorithm

BACKGROUND

The utilization of long reads for single nucleotide polymorphism (SNP) phasing has become popular, providing substantial support for research on human diseases and genetic studies in animals and plants. However, due to the complexity of the linkage relationships between SNP loci and sequencing errors in the reads, the recent methods still cannot yield satisfactory results.

RESULTS

In this study, we present a graph-based algorithm, GCphase, which utilizes the minimum cut algorithm to perform phasing. First, based on alignment between long reads and the reference genome, GCphase filters out ambiguous SNP sites and useless read information. Second, GCphase constructs a graph in which a vertex represents alleles of an SNP locus and each edge represents the presence of read support; moreover, GCphase adopts a graph minimum-cut algorithm to phase the SNPs. Next, GCpahse uses two error correction steps to refine the phasing results obtained from the previous step, effectively reducing the error rate. Finally, GCphase obtains the phase block. GCphase was compared to three other methods, WhatsHap, HapCUT2, and LongPhase, on the Nanopore and PacBio long-read datasets. The code is available from https://github.com/baimawjy/GCphase .

CONCLUSIONS

Experimental results show that GCphase under different sequencing depths of different data has the least number of switch errors and the highest accuracy compared with other methods.

The Ruminant Telomere-to-Telomere (RT2T) Consortium

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

CGG/CCG Repeat Expansions in LOC642361/NUTM2B-AS1 in Thai Patients With Oculopharyngodistal Myopathy

Objectives

This study characterizes oculopharyngodistal myopathy in 4 Thai patients from 3 families with CGG/CCG repeat expansion in LOC642361/NUTM2B-AS1.

Methods

Repeat-primed PCR analyzed CGG/CCG repeat size in LOC642361/NUTM2B-AS1 in 4 Thai patients suspected of oculopharyngodistal myopathy (OPDM). Clinical records were reviewed for clinicopathologic features.

Results

All patients exhibited strong somatic instabilities of the expanded CGG/CCG repeats, primarily manifesting as oculopharyngeal weakness. Patient 1 had mild finger extensor and intrinsic hand muscle weakness, and although patient 2 lacked limb weakness, both siblings showed electrophysiologic evidence of distal myopathy, indicative of OPDM. Patient 3, the daughter of a sibling with OPDM reported in 2004, lacked limb weakness or leukoencephalopathy on brain MRI. Patient 4, initially misdiagnosed with refractory myasthenia gravis, had generalized muscle weakness.

Discussion

While initially characterized as oculopharyngeal myopathy with leukoencephalopathy (OPML) in a Japanese family, our study suggests a stronger association between CGG/CCG expansion in LOC642361/NUTM2B-AS1 and oculopharyngodistal myopathy (OPDM) rather than OPML. The variable presence or absence of leukoencephalopathy further supports OPDM as the predominant clinical manifestation linked to CGG/CCG expansion in LOC642361/NUTM2B-AS1.

The correlation between CpG methylation and gene expression is driven by sequence variants

Gene promoter and enhancer sequences are bound by transcription factors and are depleted of methylated CpG sites (cytosines preceding guanines in DNA). The absence of methylated CpGs in these sequences typically correlates with increased gene expression, indicating a regulatory role for methylation. We used nanopore sequencing to determine haplotype-specific methylation rates of 15.3 million CpG units in 7,179 whole-blood genomes. We identified 189,178 methylation depleted sequences where three or more proximal CpGs were unmethylated on at least one haplotype. A total of 77,789 methylation depleted sequences (~41%) associated with 80,503 cis-acting sequence variants, which we termed allele-specific methylation quantitative trait loci (ASM-QTLs). RNA sequencing of 896 samples from the same blood draws used to perform nanopore sequencing showed that the ASM-QTL, that is, DNA sequence variability, drives most of the correlation found between gene expression and CpG methylation. ASM-QTLs were enriched 40.2-fold (95% confidence interval 32.2, 49.9) among sequence variants associating with hematological traits, demonstrating that ASM-QTLs are important functional units in the noncoding genome.

Rapid DNA methylation-based classification of pediatric brain tumors from ultrasonic aspirate specimens

BACKGROUND

Although cavitating ultrasonic aspirators are commonly used in neurosurgical procedures, the suitability of ultrasonic aspirator-derived tumor material for diagnostic procedures is still controversial. Here, we explore the feasibility of using ultrasonic aspirator-resected tumor tissue to classify otherwise discarded sample material by fast DNA methylation-based analysis using low pass nanopore whole genome sequencing.

METHODS

Ultrasonic aspirator-derived specimens from pediatric patients undergoing brain tumor resection were subjected to low-pass nanopore whole genome sequencing. DNA methylation-based classification using a neural network classifier and copy number variation analysis were performed. Tumor purity was estimated from copy number profiles. Results were compared to microarray (EPIC)-based routine neuropathological histomorphological and molecular evaluation.

RESULTS

19 samples with confirmed neuropathological diagnosis were evaluated. All samples were successfully sequenced and passed quality control for further analysis. DNA and sequencing characteristics from ultrasonic aspirator-derived specimens were comparable to routinely processed tumor tissue. Classification of both methods was concordant regarding methylation class in 17/19 (89%) cases. Application of a platform-specific threshold for nanopore-based classification ensured a specificity of 100%, whereas sensitivity was 79%. Copy number variation profiles were generated for all cases and matched EPIC results in 18/19 (95%) samples, even allowing the identification of diagnostically or therapeutically relevant genomic alterations.

CONCLUSION

Methylation-based classification of pediatric CNS tumors based on ultrasonic aspirator-reduced and otherwise discarded tissue is feasible using time- and cost-efficient nanopore sequencing.

Ecological genomics in the Northern krill uncovers loci for local adaptation across ocean basins

Krill are vital as food for many marine animals but also impacted by global warming. To learn how they and other zooplankton may adapt to a warmer world we studied local adaptation in the widespread Northern krill (Meganyctiphanes norvegica). We assemble and characterize its large genome and compare genome-scale variation among 74 specimens from the colder Atlantic Ocean and warmer Mediterranean Sea. The 19 Gb genome likely evolved through proliferation of retrotransposons, now targeted for inactivation by extensive DNA methylation, and contains many duplicated genes associated with molting and vision. Analysis of 760 million SNPs indicates extensive homogenizing gene-flow among populations. Nevertheless, we detect signatures of adaptive divergence across hundreds of genes, implicated in photoreception, circadian regulation, reproduction and thermal tolerance, indicating polygenic adaptation to light and temperature. The top gene candidate for ecological adaptation was nrf-6, a lipid transporter with a Mediterranean variant that may contribute to early spring reproduction. Such variation could become increasingly important for fitness in Atlantic stocks. Our study underscores the widespread but uneven distribution of adaptive variation, necessitating characterization of genetic variation among natural zooplankton populations to understand their adaptive potential, predict risks and support ocean conservation in the face of climate change.

DeepBAM: a high-accuracy single-molecule CpG methylation detection tool for Oxford nanopore sequencing

Recent nanopore sequencing system (R10.4) has enhanced base calling accuracy and is being increasingly utilized for detecting CpG methylation state. However, the robustness and universality of the methylation calling model in officially supplied Dorado remains poorly tested. In this study, we obtained heterogeneous datasets from human and plant sources to carry out comprehensive evaluations, which showed that Dorado performed significantly different across datasets. We therefore developed deep neural networks and implemented several optimizations in training a new model called DeepBAM. DeepBAM achieved superior and more stable performances compared with Dorado, including higher area under the ROC curves (98.47% on average and up to 7.36% improvement) and F1 scores (94.97% on average and up to 16.24% improvement) across the datasets. DeepBAM-based whole genome methylation frequencies have achieved >0.95 correlations with BS-seq on four of five datasets, outperforming Dorado in all instances. It enables unraveling allele-specific methylation patterns, including regions of transposable elements. The enhanced performance of DeepBAM paves the way for broader applications of nanopore sequencing in CpG methylation studies.

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

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

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

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

Exploring the genetic and epigenetic underpinnings of early-onset cancers: Variant prioritization for long read whole genome sequencing from family cancer pedigrees

Despite significant advances in our understanding of genetic cancer susceptibility, known inherited cancer predisposition syndromes explain at most 20% of early-onset cancers. As early-onset cancer prevalence continues to increase, the need to assess previously inaccessible areas of the human genome, harnessing a trio or quad family-based architecture for variant filtration, may reveal further insights into cancer susceptibility. To assess a broader spectrum of variation than can be ascertained by multi-gene panel sequencing, or even whole genome sequencing with short reads, we employed long read whole genome sequencing using an Oxford Nanopore Technology (ONT) PromethION of 3 families containing an early-onset cancer proband using a trio or quad family architecture. Analysis included 2 early-onset colorectal cancer family trios and one quad consisting of two siblings with testicular cancer, all with unaffected parents. Structural variants (SVs), epigenetic profiles and single nucleotide variants (SNVs) were determined for each individual, and a filtering strategy was employed to refine and prioritize candidate variants based on the family architecture. The family architecture enabled us to focus on inapposite variants while filtering variants shared with the unaffected parents, significantly decreasing background variation that can hamper identification of potentially disease causing differences. Candidate d e novo and compound heterozygous variants were identified in this way. Gene expression, in matched neoplastic and pre-neoplastic lesions, was assessed for one trio. Our study demonstrates the feasibility of a streamlined analysis of genomic variants from long read ONT whole genome sequencing and a way to prioritize key variants for further evaluation of pathogenicity, while revealing what may be missing from panel based analyses.

LINE-1 retrotransposons contribute to mouse PV interneuron development

Retrotransposons are mobile DNA sequences duplicated via transcription and reverse transcription of an RNA intermediate. Cis-regulatory elements encoded by retrotransposons can also promote the transcription of adjacent genes. Somatic LINE-1 (L1) retrotransposon insertions have been detected in mammalian neurons. It is, however, unclear whether L1 sequences are mobile in only some neuronal lineages or therein promote neurodevelopmental gene expression. Here we report programmed L1 activation by SOX6, a transcription factor critical for parvalbumin (PV) interneuron development. Mouse PV interneurons permit L1 mobilization in vitro and in vivo, harbor unmethylated L1 promoters and express full-length L1 mRNAs and proteins. Using nanopore long-read sequencing, we identify unmethylated L1s proximal to PV interneuron genes, including a novel L1 promoter-driven Caps2 transcript isoform that enhances neuron morphological complexity in vitro. These data highlight the contribution made by L1 cis-regulatory elements to PV interneuron development and transcriptome diversity, uncovered due to L1 mobility in this milieu.

Viral genome sequencing methods: benefits and pitfalls of current approaches

Whole genome sequencing of viruses provides high-resolution molecular insights, enhancing our understanding of viral genome function and phylogeny. Beyond fundamental research, viral sequencing is increasingly vital for pathogen surveillance, epidemiology, and clinical applications. As sequencing methods rapidly evolve, the diversity of viral genomics applications and catalogued genomes continues to expand. Advances in long-read, single molecule, real-time sequencing methodologies present opportunities to sequence contiguous, haplotype resolved viral genomes in a range of research and applied settings. Here we present an overview of nucleic acid sequencing methods and their applications in studying viral genomes. We emphasise the advantages of different viral sequencing approaches, with a particular focus on the benefits of third-generation sequencing technologies in elucidating viral evolution, transmission networks, and pathogenesis.

Simulation of nanopore sequencing signal data with tunable parameters

In silico simulation of high-throughput sequencing data is a technique used widely in the genomics field. However, there is currently a lack of effective tools for creating simulated data from nanopore sequencing devices, which measure DNA or RNA molecules in the form of time-series current signal data. Here, we introduce Squigulator, a fast and simple tool for simulation of realistic nanopore signal data. Squigulator takes a reference genome, a transcriptome, or read sequences, and generates corresponding raw nanopore signal data. This is compatible with basecalling software from Oxford Nanopore Technologies (ONT) and other third-party tools, thereby providing a useful substrate for development, testing, debugging, validation, and optimization at every stage of a nanopore analysis workflow. The user may generate data with preset parameters emulating specific ONT protocols or noise-free "ideal" data, or they may deterministically modify a range of experimental variables and/or noise parameters to shape the data to their needs. We present a brief example of Squigulator's use, creating simulated data to model the degree to which different parameters impact the accuracy of ONT basecalling and downstream variant detection. This analysis reveals new insights into the nature of ONT data and basecalling algorithms. We provide Squigulator as an open-source tool for the nanopore community.

MethPhaser: methylation-based long-read haplotype phasing of human genomes

The assignment of variants across haplotypes, phasing, is crucial for predicting the consequences, interaction, and inheritance of mutations and is a key step in improving our understanding of phenotype and disease. However, phasing is limited by read length and stretches of homozygosity along the genome. To overcome this limitation, we designed MethPhaser, a method that utilizes methylation signals from Oxford Nanopore Technologies to extend Single Nucleotide Variation (SNV)-based phasing. We demonstrate that haplotype-specific methylations extensively exist in Human genomes and the advent of long-read technologies enabled direct report of methylation signals. For ONT R9 and R10 cell line data, we increase the phase length N50 by 78%-151% at a phasing accuracy of 83.4-98.7% To assess the impact of tissue purity and random methylation signals due to inactivation, we also applied MethPhaser on blood samples from 4 patients, still showing improvements over SNV-only phasing. MethPhaser further improves phasing across HLA and multiple other medically relevant genes, improving our understanding of how mutations interact across multiple phenotypes. The concept of MethPhaser can also be extended to non-human diploid genomes. MethPhaser is available at https://github.com/treangenlab/methphaser .

Twin-based Mendelian Randomization Analyses Highlight Smoking's Effects on Blood DNA Methylation, with Putative Reverse Causation

Epigenome-wide association studies (EWAS) aim to identify differentially methylated loci associated with complex traits and disorders. EWAS of cigarette smoking shows some of the most widespread DNA methylation (DNAm) associations in blood. However, traditional EWAS cannot differentiate between causation and confounding, leading to ambiguity in etiological interpretations. Here, we apply an integrated approach combining Mendelian Randomization and twin-based Direction-of-Causation analyses (MR-DoC) to examine causality underlying smoking-associated blood DNAm changes in the Netherlands Twin Register (N=2577). Evidence across models suggests that current smoking's causal effects on DNAm likely drive many of the previous EWAS findings, implicating functional pathways relevant to several adverse health outcomes of smoking, including hemopoiesis, cell- and neuro-development, and immune regulation. Additionally, we find evidence of potential reverse causal influences at some DNAm sites, with 17 of these sites enriched for gene regulatory functional elements in the brain. The top three sites with evidence of DNAm's effects on smoking annotate to genes involved in G protein-coupled receptor signaling (GNG7, RGS3) and innate immune response (SLC15A4), elucidating potential biological risk factors for smoking. This study highlights the utility of integrating genotypic and DNAm measures in twin cohorts to clarify the causal relationships between health behaviors and blood DNAm.

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

Bacteria rely on DNA methylation for restriction-modification systems and epigenetic control of gene expression. Here, we use direct detection of methylated bases by nanopore sequencing to monitor global DNA methylation in Alphaproteobacteria, where use of this technique has not yet been reported. One representative of this order, Caulobacter crescentus, relies on DNA methylation to control cell cycle progression, but it is unclear whether other members of this order, such as Brucella abortus, depend on the same systems. We addressed these questions by first measuring CcrM-dependent DNA methylation in Caulobacter and showing excellent correlation between nanopore-based detection and previously published results. We then directly measure the impact of Lon-mediated CcrM degradation on the epigenome, verifying that loss of Lon results in pervasive methylation. We also show that the AlkB demethylase has no global impact on DNA methylation during normal growth. Next, we report on the global DNA methylation in B. abortus for the first time and find that CcrM-dependent methylation is reliant on Lon but impacts the two chromosomes differently. Finally, we explore the impact of the MucR transcription factor, known to compete with CcrM methylation, on the Brucella methylome and share the results with a publicly available visualization package. Our work demonstrates the utility of nanopore-based sequencing for epigenome measurements in Alphaproteobacteria and reveals new features of CcrM-dependent methylation in a zoonotic pathogen.IMPORTANCEDNA methylation plays an important role in bacteria, maintaining genome integrity and regulating gene expression. We used nanopore sequencing to directly measure methylated bases in Caulobacter crescentus and Brucella abortus. In Caulobacter, we showed that stabilization of the CcrM methyltransferase upon loss of the Lon protease results in prolific methylation and discovered that the putative methylase AlkB is unlikely to have a global physiological effect. We measured genome-wide methylation in Brucella for the first time, revealing a similar role for CcrM in cell-cycle methylation but a more complex regulation by the Lon protease than in Caulobacter. Finally, we show how the virulence factor MucR impacts DNA methylation patterns in Brucella.

Long-read sequencing for brain tumors

Brain tumors and genomics have a long-standing history given that glioblastoma was the first cancer studied by the cancer genome atlas. The numerous and continuous advances through the decades in sequencing technologies have aided in the advanced molecular characterization of brain tumors for diagnosis, prognosis, and treatment. Since the implementation of molecular biomarkers by the WHO CNS in 2016, the genomics of brain tumors has been integrated into diagnostic criteria. Long-read sequencing, also known as third generation sequencing, is an emerging technique that allows for the sequencing of longer DNA segments leading to improved detection of structural variants and epigenetics. These capabilities are opening a way for better characterization of brain tumors. Here, we present a comprehensive summary of the state of the art of third-generation sequencing in the application for brain tumor diagnosis, prognosis, and treatment. We discuss the advantages and potential new implementations of long-read sequencing into clinical paradigms for neuro-oncology patients.

Decoding chromosome organization using CheC-PLS: chromosome conformation by proximity labeling and long-read sequencing

Genomic approaches have provided detailed insight into chromosome architecture. However, commonly deployed techniques do not preserve connectivity-based information, leaving large-scale genome organization poorly characterized. Here, we developed CheC-PLS: a proximity-labeling technique that indelibly marks, and then decodes, protein-associated sites. CheC-PLS tethers dam methyltransferase to a protein of interest, followed by Nanopore sequencing to identify methylated bases - indicative of in vivo proximity - along reads >100kb. As proof-of-concept we analyzed, in budding yeast, a cohesin-based meiotic backbone that organizes chromatin into an array of loops. Our data recapitulates previously obtained association patterns, and, importantly, exposes variability between cells. Single read data reveals cohesin translocation on DNA and, by anchoring reads onto unique regions, we define the internal organization of the ribosomal DNA locus. Our versatile technique, which we also deployed on isolated nuclei with nanobodies, promises to illuminate diverse chromosomal processes by describing the in vivo conformations of single chromosomes.

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

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

Direct transposition of native DNA for sensitive multimodal single-molecule sequencing

Concurrent readout of sequence and base modifications from long unamplified DNA templates by Pacific Biosciences of California (PacBio) single-molecule sequencing requires large amounts of input material. Here we adapt Tn5 transposition to introduce hairpin oligonucleotides and fragment (tagment) limiting quantities of DNA for generating PacBio-compatible circular molecules. We developed two methods that implement tagmentation and use 90-99% less input than current protocols: (1) single-molecule real-time sequencing by tagmentation (SMRT-Tag), which allows detection of genetic variation and CpG methylation; and (2) single-molecule adenine-methylated oligonucleosome sequencing assay by tagmentation (SAMOSA-Tag), which uses exogenous adenine methylation to add a third channel for probing chromatin accessibility. SMRT-Tag of 40 ng or more human DNA (approximately 7,000 cell equivalents) yielded data comparable to gold standard whole-genome and bisulfite sequencing. SAMOSA-Tag of 30,000-50,000 nuclei resolved single-fiber chromatin structure, CTCF binding and DNA methylation in patient-derived prostate cancer xenografts and uncovered metastasis-associated global epigenome disorganization. Tagmentation thus promises to enable sensitive, scalable and multimodal single-molecule genomics for diverse basic and clinical applications.

Advances in single-cell long-read sequencing technologies

With an increase in accuracy and throughput of long-read sequencing technologies, they are rapidly being assimilated into the single-cell sequencing pipelines. For transcriptome sequencing, these techniques provide RNA isoform-level information in addition to the gene expression profiles. Long-read sequencing technologies not only help in uncovering complex patterns of cell-type specific splicing, but also offer unprecedented insights into the origin of cellular complexity and thus potentially new avenues for drug development. Additionally, single-cell long-read DNA sequencing enables high-quality assemblies, structural variant detection, haplotype phasing, resolving high-complexity regions, and characterization of epigenetic modifications. Given that significant progress has primarily occurred in single-cell RNA isoform sequencing (scRiso-seq), this review will delve into these advancements in depth and highlight the practical considerations and operational challenges, particularly pertaining to downstream analysis. We also aim to offer a concise introduction to complementary technologies for single-cell sequencing of the genome, epigenome and epitranscriptome. We conclude by identifying certain key areas of innovation that may drive these technologies further and foster more widespread application in biomedical science.

Comparative genomics of the closely related fungal genera Cryptococcus and Kwoniella reveals karyotype dynamics and suggests evolutionary mechanisms of pathogenesis

In exploring the evolutionary trajectories of both pathogenesis and karyotype dynamics in fungi, we conducted a large-scale comparative genomic analysis spanning the Cryptococcus genus, encompassing both global human fungal pathogens and nonpathogenic species, and related species from the sister genus Kwoniella. Chromosome-level genome assemblies were generated for multiple species, covering virtually all known diversity within these genera. Although Cryptococcus and Kwoniella have comparable genome sizes (about 19.2 and 22.9 Mb) and similar gene content, hinting at preadaptive pathogenic potential, our analysis found evidence of gene gain (via horizontal gene transfer) and gene loss in pathogenic Cryptococcus species, which might represent evolutionary signatures of pathogenic development. Genome analysis also revealed a significant variation in chromosome number and structure between the 2 genera. By combining synteny analysis and experimental centromere validation, we found that most Cryptococcus species have 14 chromosomes, whereas most Kwoniella species have fewer (11, 8, 5, or even as few as 3). Reduced chromosome number in Kwoniella is associated with formation of giant chromosomes (up to 18 Mb) through repeated chromosome fusion events, each marked by a pericentric inversion and centromere loss. While similar chromosome inversion-fusion patterns were observed in all Kwoniella species with fewer than 14 chromosomes, no such pattern was detected in Cryptococcus. Instead, Cryptococcus species with less than 14 chromosomes showed reductions primarily through rearrangements associated with the loss of repeat-rich centromeres. Additionally, Cryptococcus genomes exhibited frequent interchromosomal translocations, including intercentromeric recombination facilitated by transposons shared between centromeres. Overall, our findings advance our understanding of genetic changes possibly associated with pathogenicity in Cryptococcus and provide a foundation to elucidate mechanisms of centromere loss and chromosome fusion driving distinct karyotypes in closely related fungal species, including prominent global human pathogens.

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

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

Unveiling microbial diversity: harnessing long-read sequencing technology

Long-read sequencing has recently transformed metagenomics, enhancing strain-level pathogen characterization, enabling accurate and complete metagenome-assembled genomes, and improving microbiome taxonomic classification and profiling. These advancements are not only due to improvements in sequencing accuracy, but also happening across rapidly changing analysis methods. In this Review, we explore long-read sequencing's profound impact on metagenomics, focusing on computational pipelines for genome assembly, taxonomic characterization and variant detection, to summarize recent advancements in the field and provide an overview of available analytical methods to fully leverage long reads. We provide insights into the advantages and disadvantages of long reads over short reads and their evolution from the early days of long-read sequencing to their recent impact on metagenomics and clinical diagnostics. We further point out remaining challenges for the field such as the integration of methylation signals in sub-strain analysis and the lack of benchmarks.

Nanopore Current Events Magnifier (nanoCEM): a novel tool for visualizing current events at modification sites of nanopore sequencing

Nanopore sequencing technologies have enabled the direct detection of base modifications in DNA or RNA molecules. Despite these advancements, the tools for visualizing electrical current, essential for analyzing base modifications, are often lacking in clarity and compatibility with diverse nanopore pipelines. Here, we present Nanopore Current Events Magnifier (nanoCEM, https://github.com/lrslab/nanoCEM), a Python command-line tool designed to facilitate the identification of DNA/RNA modification sites through enhanced visualization and statistical analysis. Compatible with the four preprocessing methods including 'f5c resquiggle', 'f5c eventalign', 'Tombo' and 'move table', nanoCEM is applicable to RNA and DNA analysis across multiple flow cell types. By utilizing rescaling techniques and calculating various statistical features, nanoCEM provides more accurate and comparable visualization of current events, allowing researchers to effectively observe differences between samples and showcase the modified sites.

An epigenetic toolbox for conservation biologists

Ongoing climatic shifts and increasing anthropogenic pressures demand an efficient delineation of conservation units and accurate predictions of populations' resilience and adaptive potential. Molecular tools involving DNA sequencing are nowadays routinely used for these purposes. Yet, most of the existing tools focusing on sequence-level information have shortcomings in detecting signals of short-term ecological relevance. Epigenetic modifications carry valuable information to better link individuals, populations, and species to their environment. Here, we discuss a series of epigenetic monitoring tools that can be directly applied to various conservation contexts, complementing already existing molecular monitoring frameworks. Focusing on DNA sequence-based methods (e.g. DNA methylation, for which the applications are readily available), we demonstrate how (a) the identification of epi-biomarkers associated with age or infection can facilitate the determination of an individual's health status in wild populations; (b) whole epigenome analyses can identify signatures of selection linked to environmental conditions and facilitate estimating the adaptive potential of populations; and (c) epi-eDNA (epigenetic environmental DNA), an epigenetic-based conservation tool, presents a non-invasive sampling method to monitor biological information beyond the mere presence of individuals. Overall, our framework refines conservation strategies, ensuring a comprehensive understanding of species' adaptive potential and persistence on ecologically relevant timescales.

DNA modifications: Biomarkers for the exposome?

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

Single-base resolution quantitative genome methylation analysis in the model bacterium Helicobacter pylori by enzymatic methyl sequencing (EM-Seq) reveals influence of strain, growth phase, and methyl homeostasis

BACKGROUND

Bacterial epigenetics is a rapidly expanding research field. DNA methylation by diverse bacterial methyltransferases (MTases) contributes to genomic integrity and replication, and many recent studies extended MTase function also to global transcript regulation and phenotypic variation. Helicobacter pylori is currently one of those bacterial species which possess the highest number and the most variably expressed set of DNA MTases. Next-generation sequencing technologies can directly detect DNA base methylation. However, they still have limitations in their quantitative and qualitative performance, in particular for cytosine methylation.

RESULTS

As a complementing approach, we used enzymatic methyl sequencing (EM-Seq), a technology recently established that has not yet been fully evaluated for bacteria. Thereby, we assessed quantitatively, at single-base resolution, whole genome cytosine methylation for all methylated cytosine motifs in two different H. pylori strains and isogenic MTase mutants. EM-Seq reliably detected both m5C and m4C methylation. We demonstrated that three different active cytosine MTases in H. pylori provide considerably different levels of average genome-wide single-base methylation, in contrast to isogenic mutants which completely lost specific motif methylation. We found that strain identity and changed environmental conditions, such as growth phase and interference with methyl donor homeostasis, significantly influenced quantitative global and local genome-wide methylation in H. pylori at specific motifs. We also identified significantly hyper- or hypo-methylated cytosines, partially linked to overlapping MTase target motifs. Notably, we revealed differentially methylated cytosines in genome-wide coding regions under conditions of methionine depletion, which can be linked to transcript regulation.

CONCLUSIONS

This study offers new knowledge on H. pylori global and local genome-wide methylation and establishes EM-Seq for quantitative single-site resolution analyses of bacterial cytosine methylation.

Examining chromatin heterogeneity through PacBio long-read sequencing of M.EcoGII methylated genomes: an m6A detection efficiency and calling bias correcting pipeline

Recent studies have combined DNA methyltransferase footprinting of genomic DNA in nuclei with long-read sequencing, resulting in detailed chromatin maps for multi-kilobase stretches of genomic DNA from one cell. Theoretically, nucleosome footprints and nucleosome-depleted regions can be identified using M.EcoGII, which methylates adenines in any sequence context, providing a high-resolution map of accessible regions in each DNA molecule. Here, we report PacBio long-read sequence data for budding yeast nuclei treated with M.EcoGII and a bioinformatic pipeline which corrects for three key challenges undermining this promising method. First, detection of m6A in individual DNA molecules by the PacBio software is inefficient, resulting in false footprints predicted by random gaps of seemingly unmethylated adenines. Second, there is a strong bias against m6A base calling as AT content increases. Third, occasional methylation occurs within nucleosomes, breaking up their footprints. After correcting for these issues, our pipeline calculates a correlation coefficient-based score indicating the extent of chromatin heterogeneity within the cell population for every gene. Although the population average is consistent with that derived using other techniques, we observe a wide range of heterogeneity in nucleosome positions at the single-molecule level, probably reflecting cellular chromatin dynamics.

The Applications of Nanopore Sequencing Technology in Animal and Human Virus Research

In recent years, an increasing number of viruses have triggered outbreaks that pose a severe threat to both human and animal life, as well as caused substantial economic losses. It is crucial to understand the genomic structure and epidemiology of these viruses to guide effective clinical prevention and treatment strategies. Nanopore sequencing, a third-generation sequencing technology, has been widely used in genomic research since 2014. This technology offers several advantages over traditional methods and next-generation sequencing (NGS), such as the ability to generate ultra-long reads, high efficiency, real-time monitoring and analysis, portability, and the ability to directly sequence RNA or DNA molecules. As a result, it exhibits excellent applicability and flexibility in virus research, including viral detection and surveillance, genome assembly, the discovery of new variants and novel viruses, and the identification of chemical modifications. In this paper, we provide a comprehensive review of the development, principles, advantages, and applications of nanopore sequencing technology in animal and human virus research, aiming to offer fresh perspectives for future studies in this field.

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

BACKGROUND

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

METHODS

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

RESULTS

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

CONCLUSIONS

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

Prediction of m6A and m5C at single-molecule resolution reveals a transcriptome-wide co-occurrence of RNA modifications

The epitranscriptome embodies many new and largely unexplored functions of RNA. A significant roadblock hindering progress in epitranscriptomics is the identification of more than one modification in individual transcript molecules. We address this with CHEUI (CH3 (methylation) Estimation Using Ionic current). CHEUI predicts N6-methyladenosine (m6A) and 5-methylcytosine (m5C) in individual molecules from the same sample, the stoichiometry at transcript reference sites, and differential methylation between any two conditions. CHEUI processes observed and expected nanopore direct RNA sequencing signals to achieve high single-molecule, transcript-site, and stoichiometry accuracies in multiple tests using synthetic RNA standards and cell line data. CHEUI's capability to identify two modification types in the same sample reveals a co-occurrence of m6A and m5C in individual mRNAs in cell line and tissue transcriptomes. CHEUI provides new avenues to discover and study the function of the epitranscriptome.

DNA Methylation Analysis of Growth Differences between Upright and Inverted Cuttings of Populus yunnanensis

DNA methylation is an important mechanism for epigenetic modifications that have been shown to be associated with responses to plant development. Previous studies found that inverted Populus yunnanensis cuttings were still viable and could develop into complete plants. However, the growth status of inverted cuttings was weaker than that of upright cuttings, and the sprouting time of inverted cuttings was later than that of upright cuttings. There is currently no research on DNA methylation patterns in inverted cuttings of Populus yunnanensis. In this study, we detected genome-wide methylation patterns of stem tips of Populus yunnanensis at the early growth stage and the rapid growth stage by Oxford Nanopore Technologies (ONT) methylation sequencing. We found that the methylation levels of CpG, CHG, CHH, and 6mA were 41.34%, 33.79%, 17.27%, and 12.90%, respectively, in the genome of inverted poplar cuttings, while the methylation levels of the four methylation types were higher in the genome of upright poplar cuttings than in inverted cuttings, 41.90%, 34.57%, 18.09%, and 14.11%, suggesting important roles for DNA methylation in poplar cells. In all comparison groups, CpG-type methylation genes in the Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway were annotated to pathways associated with carbon metabolism, ribosome biogenesis in eukaryotes, glycolysis/gluconeogenesis, pyruvate metabolism, and mRNA detection pathways, suggesting that different biological processes are activated in upright and inverted cuttings. The results show that methylation genes are commonly present in the poplar genome, but only a few of them are involved in the regulation of expression in the growth and development of inverted cuttings. From this, we screened the DET2 gene for significant differences in methylation levels in upright or inverted cuttings. The DET2 gene is a key gene in the Brassinolide (BRs) synthesis pathway, and BRs have an important influence on the growth and development process of plants. These results provide important clues for studying DNA methylation patterns in P. yunnanensis.

MASTR-seq: Multiplexed Analysis of Short Tandem Repeats with sequencing

More than 60 human disorders have been linked to unstable expansion of short tandem repeat (STR) tracts. STR length and the extent of DNA methylation is linked to disease pathology and can be mosaic in a cell type-specific manner in several repeat expansion disorders. Mosaic phenomenon have been difficult to study to date due to technical bias intrinsic to repeat sequences and the need for multi-modal measurements at single-allele resolution. Nanopore long-read sequencing accurately measures STR length and DNA methylation in the same single molecule but is cost prohibitive for studies assessing a target locus across multiple experimental conditions or patient samples. Here, we describe MASTR-seq, M ultiplexed A nalysis of S hort T andem R epeats, for cost-effective, high-throughput, accurate, multi-modal measurements of DNA methylation and STR genotype at single-allele resolution. MASTR-seq couples long-read sequencing, Cas9-mediated target enrichment, and PCR-free multiplexed barcoding to achieve a >ten-fold increase in on-target read mapping for 8-12 pooled samples in a single MinION flow cell. We provide a detailed experimental protocol and computational tools and present evidence that MASTR-seq quantifies tract length and DNA methylation status for CGG and CAG STR loci in normal-length and mutation-length human cell lines. The MASTR-seq protocol takes approximately eight days for experiments and one additional day for data processing and analyses.

Key points

We provide a protocol for MASTR-seq: M ultiplexed A nalysis of S hort T andem R epeats using Cas9-mediated target enrichment and PCR-free, multiplexed nanopore sequencing. MASTR-seq achieves a >10-fold increase in on-target read proportion for highly repetitive, technically inaccessible regions of the genome relevant for human health and disease.MASTR-seq allows for high-throughput, efficient, accurate, and cost-effective measurement of STR length and DNA methylation in the same single allele for up to 8-12 samples in parallel in one Nanopore MinION flow cell.

Advancements in long-read genome sequencing technologies and algorithms

The recent advent of long read sequencing technologies, such as Pacific Biosciences (PacBio) and Oxford Nanopore technology (ONT), have led to substantial improvements in accuracy and computational cost in sequencing genomes. However, de novo whole-genome assembly still presents significant challenges related to the quality of the results. Pursuing de novo whole-genome assembly remains a formidable challenge, underscored by intricate considerations surrounding computational demands and result quality. As sequencing accuracy and throughput steadily advance, a continuous stream of innovative assembly tools floods the field. Navigating this dynamic landscape necessitates a reasonable choice of sequencing platform, depth, and assembly tools to orchestrate high-quality genome reconstructions. This comprehensive review delves into the intricate interplay between cutting-edge long read sequencing technologies, assembly methodologies, and the ever-evolving field of genomics. With a focus on addressing the pivotal challenges and harnessing the opportunities presented by these advancements, we provide an in-depth exploration of the crucial factors influencing the selection of optimal strategies for achieving robust and insightful genome assemblies.

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

Introduction

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

Methods

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

Results

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

Discussion

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

Restriction modification systems in archaea: A panoramic outlook

Restriction modification (RM) systems are one of the ubiquitous yet primitive defense responses employed by bacteria and archaea with the primary role of safeguarding themselves against invading bacteriophages. Protection of the host occurs by the cleavage of the invading foreign DNA via restriction endonucleases with concomitant methylation of host DNA with the aid of a methyltransferase counterpart. RM systems have been extensively studied in bacteria, however, in the case of archaea there are limited reports of RM enzymes that are investigated to date owing to their inhospitable growth demands. This review aims to broaden the knowledge about what is known about the diversity of RM systems in archaea and encapsulate the current knowledge on restriction and modification enzymes characterized in archaea so far and the role of RM systems in the milieu of archaeal biology.

Detection of ribonucleotides embedded in DNA by Nanopore sequencing

Ribonucleotides represent the most common non-canonical nucleotides found in eukaryotic genomes. The sources of chromosome-embedded ribonucleotides and the mechanisms by which unrepaired rNMPs trigger genome instability and human pathologies are not fully understood. The available sequencing technologies only allow to indirectly deduce the genomic location of rNMPs. Oxford Nanopore Technologies (ONT) may overcome such limitation, revealing the sites of rNMPs incorporation in genomic DNA directly from raw sequencing signals. We synthesized two types of DNA molecules containing rNMPs at known or random positions and we developed data analysis pipelines for DNA-embedded ribonucleotides detection by ONT. We report that ONT can identify all four ribonucleotides incorporated in DNA by capturing rNMPs-specific alterations in nucleotide alignment features, current intensity, and dwell time. We propose that ONT may be successfully employed to directly map rNMPs in genomic DNA and we suggest a strategy to build an ad hoc basecaller to analyse native genomes.