Search PubMed⌕ Search

SEARCH · Search PubMed

Results for “nanopore sequencing”

Search indexed PubMed citations on genomics, clinical trials, systematic reviews and public health. Explore titles, authors and supplied subject terms, then open the PubMed record.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 73 records · Page 4Linked to original sources

Genomic characterization of KPC-2 and NDM coproducing carbapenem-resistant Klebsiella pneumoniae in a hospital: discovery of ST1869 clone and a novel hybrid plasmid.

UNLABELLED: To characterize the plasmid architecture and molecular background of KPC-NDM coproducing carbapenem-resistant Klebsiella pneumoniae (KN-CRKP) in a South China hospital. Five KN-CRKP isolates were collected, including three from one patient. All underwent Illumina sequencing; two (ST11 and ST1869) additionally had Nanopore sequencing. Antimicrobial susceptibility testing strain sequence types, conjugation assays, resistance gene profiling, plasmid typing, genetic structure comparison, core-genome single nucleotide polymorphisms (SNPs) analysis, and plasmid clustering were performed. All isolates exhibited an imipenem minimum inhibitory concentration (MIC) of ≥128 µg/mL and harbored multiple resistance genes. One isolate (1/5) belonged to ST1869 and co-harbored blaKPC-2 and blaNDM-5. The blaNDM-5-carrying plasmid was a novel IncI1/X3 fusion plasmid that also carried blaCMY-42. Unlike several IncX3 plasmids carrying blaNDM in publicly available KN-CRKP genomes from South China, this IncI1/X3 hybrid lacked a complete conjugative transfer system. ST11 was the predominant clone (4/5), co-harboring blaKPC-2 and blaNDM-1. A rare genetic structure, ΔISKpn6-blaKPC-2-ISKpn28, was identified on IncFII plasmids carrying blaKPC-2. Plasmid clustering analysis of 126 comparative KN-CRKP genomes showed diverse sequence types and plasmid backgrounds associated with the KPC/NDM co-production pattern. The observed plasmid diversity and structural variation in KN-CRKP support continued genomic surveillance, with particular attention to the ST1869 clone, the novel IncI1/X3 hybrid plasmid harboring blaNDM-5 and blaCMY-42, and the rare "ΔISKpn6-blaKPC-2-ISKpn28" genetic structure. Expanded genomic data on KN-CRKP are needed to further elucidate its resistance mechanisms and plasmid evolutionary trajectories. IMPORTANCE: The co-production of KPC and NDM carbapenemases in Klebsiella pneumoniae poses a formidable threat to clinical antimicrobial therapy, as these enzymes confer resistance to virtually all β-lactam agents, including carbapenems. Here, we report novel genomic features of KN-CRKP in South China, including the emergence of the ST1869 clone, a unique IncI1/X3 hybrid plasmid harboring blaNDM-5 and blaCMY-42, and the rare ΔISKpn6-blaKPC-2-ISKpn28 genetic structure. These findings substantially expand current understanding of plasmid evolution and resistance gene dissemination in this region. The identification of diverse resistance mechanisms and clonal backgrounds supports enhanced genomic surveillance and infection-control awareness for pan-resistant Enterobacterales.

Plasmids↗

Portable metagenomics for preventive surveillance and outbreak control in livestock and poultry: Pathogen detection, resistome profiling, and antimicrobial stewardship.

Conventional diagnostics for livestock and poultry outbreaks commonly rely on culture or targeted PCR panels, which may be too slow or too narrow to guide early control decisions. Portable metagenomics, particularly real-time nanopore sequencing, offers a route to broad pathogen detection, antimicrobial-resistance gene profiling, and outbreak investigation within an integrated workflow. This implementation-focused review evaluates how near-point-of-care metagenomics may support preventive veterinary medicine through earlier detection, surveillance, cohorting, biosecurity decisions, and antimicrobial stewardship. We synthesize sample-to-answer workflows for enteric and respiratory disease in food-producing animals, including sampling, nucleic-acid extraction, host depletion or target enrichment, library preparation, sequencing, bioinformatics, quality control, and interpretation. Applications in calf diarrhea, bovine respiratory disease, poultry outbreaks, mastitis, and resistome monitoring are considered alongside the central limitation that detection alone does not establish causation. Pathogen and resistance-gene signals must therefore be interpreted with clinical signs, lesions, epidemiology, controls, and confirmatory testing. We also propose a minimum reporting checklist, intended as a practical framework rather than a validated consensus standard. Portable metagenomics is not a replacement for conventional diagnostics, but appropriately validated workflows can reduce uncertainty during time-sensitive outbreaks and support more judicious antimicrobial use.

Animals↗

A de novo algorithm for allele reconstruction from Oxford nanopore amplicon reads, with application to CYP2D6.

MOTIVATION: The Oxford Nanopore Technologies' sequencing platform offers a path towards bedside genomics, producing long reads that can completely cover a gene of interest, and detect any known or novel variant the gene contains. However, the analysis of these long reads to identify actionable genotypes remains challenging and typically requires customization depending on the target gene. RESULTS: Here, we describe a generic algorithm to accurately reconstruct allele sequences derived from long-reads of amplicon-based data. Rather than calling variants directly from these long-reads, our method takes a "sequence-first" approach, performing an unbiased reconstruction of the underlying amplicon sequences to generate high-confidence reconstructed allele sequences. This is done without user input of the target gene, allowing for any source amplicon to be reconstructed. These high-confidence reconstructed allele sequences are then compared to the genomic reference sequence of the gene to infer the specific diplotype present in the sample. This approach is agnostic towards the number of genes and alleles present and readily detects novel variants. We demonstrate our approach using three independent data sets for CYP2D6, a diverse and complex gene with over 175 known alleles of clinical significance. We show how our approach can accurately recover validated CYP2D6 diplotypes from 20 Coriell samples covering 14 distinct alleles, using different amplicons, flow cell versions, and depths. This includes inferring occurrences of allele duplication events from relative abundances of each allele, a critical factor for ascribing functional effects to a diplotype. Further, we demonstrate our approach's utility for other genomic regions, including HLA. AVAILABILITY: Custom code is available at the following GitHub repository, along with instructions for use and test data: https://github.com/scottdbrown/allele-reconstruction-long-read-amplicon-data. A snapshot of the code at the time of publication is available on Zenodo.org; doi 10.5281/zenodo.19716004. Raw .fastq sequence data for our three sequencing runs is available at the SRA under Bioproject PRJNA1357883 (https://www.ncbi.nlm.nih.gov/bioproject/1357883).

Alleles↗

ModiCal: A Targeted Calibration Workflow for Site-Specific m5C Validation by Nanopore Direct RNA Sequencing.

Accurate identification of RNA 5-methylcytidine (m5C) at the single-nucleotide resolution remains a central challenge in nanopore direct RNA sequencing (DRS). Current global scanning and modification-aware basecalling methods enable transcriptome-wide profiling but often yield high false-positive rates and lack site-specific accuracy. To address this, we repurposed ModiDeC, originally a de novo multimodification classifier, into a targeted, high-precision validation tool for RNA modification sites with prior biochemical knowledge. This was implemented through a three-step calibration workflow that alternates between biochemical and computational modules using the well-characterized m5C2278 site in 25S rRNA as a starting point. Baseline training uses short synthetic RNAs carrying either a methylated or unmodified C2278 as ground truth, followed by IVT-derived calibration and validation in methyltransferase knockout yeast. The baseline model accurately detected the bona fide m5C2278 site but initially produced off-target predictions. Iterative retraining with unmodified IVT signals progressively reduced and ultimately eliminated false positives while maintaining a strong signal at the bona fide site. The final model retained enzyme-dependent detection in wild-type versus knockout yeast and, when explicitly targeted, was also able to detect the second rRNA site, C2870, which remained invisible in the initial analysis. Application to native human prerRNA processing intermediates further resolved two distinct m5C deposition regimes on 28S rRNA, while generalization to dengue virus genomic RNA confirmed that the same calibration logic transfers across diverse RNA contexts. Together, this study establishes a reproducible and transferable framework that integrates biochemical validation with iterative neural network refinement, providing a route toward reliable site-specific m5C confirmation by nanopore direct RNA sequencing.

RNA Methylation↗

Complete genome sequence of the Anaplasma phagocytophilum clinical isolate NCH-1.

Anaplasma phagocytophilum is an obligate intracellular gram-negative bacterium and etiologic agent of human granulocytic anaplasmosis. A. phagocytophilum genomic sequencing has historically been performed via short-read platforms. Our optimized bacterial isolation protocol combined with Nanopore sequencing produced a single, closed 1,481,805 bp circular A. phagocytophilum strain NCH-1 chromosome.

Anaplasma phagocytophilum↗

Genome sequence of Staphylococcus epidermidis H1G7 from hilsa (Tenualosa ilisha) gut.

We report the genome of Staphylococcus epidermidis strain H1G7 from hilsa (Tenualosa Ilisha) fish gut. Nanopore sequencing produced a 2.5-Mb assembly (32% GC, two contigs). The genome encodes metal-chelate transport, short-chain fatty acid production, and biosynthetic clusters and lacks virulent factors, revealing genomic features potentially associated with beneficial host interactions.

cyclodipeptide synthase↗

Pan-genome based on chromosome sequences of wild and cultivated Agaricus bisporus.

Agaricus bisporus, one of the most widely cultivated mushrooms around the world, plays an important role in economy and agriculture. In this study, by employing long-reads generated by PacBio and Nanopore sequencing, we assembled six novel high-quality genomes (of which three are telomere-to-telomere assemblies) with sizes 29.6 ~ 30.8 Mb and N50 lengths of 2.5 ~ 2.6 Mb. Combined with public genome data of nine strains, we successfully established a pan-genome of A. bisporus, comprising a total of 14,626 clusters of protein coding genes, of which 50.70%, 7.45%, 24.74% and 17.01% are defined as core, soft core, dispensable, and private clusters, respectively. A total of 5,646 non- redundant structural variants (SVs) were identified among wild and cultivated strains and the genes associated with SV were mapped. This work provides valuable whole-genome sequences and genomic resources across wild and cultivated strains of the most widely cultivated mushroom species for functional analyses of genomes.

Agaricus↗

Long-read sequencing to interrogate strain-level variation among adherent-invasive Escherichia coli isolated from human intestinal tissue.

Adherent-invasive Escherichia coli (AIEC) is a pathovar linked to inflammatory bowel diseases (IBD), especially Crohn's disease, and colorectal cancer. AIEC are genetically diverse, and in the absence of a universal molecular signature, are defined by in vitro functional attributes. The relative ability of difference AIEC strains to colonize, persist, and induce inflammation in an IBD-susceptible host is unresolved. To evaluate strain-level variation among tissue-associated E. coli in the intestines, we develop a long-read sequencing approach to identify AIEC by strain that excludes host DNA. We use this approach to distinguish genetically similar strains and assess their fitness in colonizing the intestine. Here we have assembled complete genomes using long-read nanopore sequencing for a model AIEC strain, NC101, and seven strains isolated from the intestinal mucosa of Crohn's disease and non-Crohn's tissues. We show these strains can colonize the intestine of IBD susceptible mice and induce inflammatory cytokines from cultured macrophages. We demonstrate that these strains can be quantified and distinguished in the presence of 99.5% mammalian DNA and from within a fecal population. Analysis of global genomic structure and specific sequence variation within the ribosomal RNA operon provides a framework for efficiently tracking strain-level variation of closely-related E. coli and likely other commensal/pathogenic bacteria impacting intestinal inflammation in experimental settings and IBD patients.

Animals↗

Assembly and characterization of the first complete mitochondrial genome of Epimedium sagittatum (Sieb. et Zucc.) Maxim (Berberidaceae):an invaluable traditional Chinese medicine.

BACKGROUND: Epimedium sagittatum (Sieb. et Zucc.) Maxim is an invaluable traditional Chinese medicine plant known for its properties of tonifying kidney yang, strengthening bones and muscles, and dispelling rheumatism. The chloroplast (cp) genome of E. sagittatum have been sequenced, offering critical insights for breeding and phylogenetic research. However, the mitochondrial (mt) genome of E. sagittatum remains uncharacterized, limiting comprehensive insights into its genomic evolution. RESULTS: In this study, we assembled the first complete mt genome of E. sagittatum employing Illumina and Nanopore sequencing technology and subsequently investigated comparative analysis with its closely related species. The mt genome of E. sagittatum was assembled as a multi-branched structure with a length of 339,191 bp, within a GC content of 46.91%. Our annotation results have shown 39 protein-coding genes (PCGs), 22 tRNA genes, three rRNA genes and four pseudogenes in the E. sagittatum mt genome. The analysis of sequence repeats has detected 79 simple sequence repeats (SSRs), 10 tandem repeats and 255 dispersed repeats in the E. sagittatum mt genome. A total of 720 C to U RNA editing sites of the 34 PCGs was predicted in E. sagittatum. The codons exhibited a strong preference for A or U bases in the E. sagittatum mt genome. The analysis of nucleotide diversity (Pi) highlighted differences in genetic variability across the tested genes, with atp9 gene exhibiting the highest genetic variation. Selection pressure analysis showed that most genes were affected by negative selection during evolution, whereas ccmB, rps10, and rps12 underwent positive selection in different plants. Additionally, a Bayesian phylogenetic tree showed that E. sagittatum was closely related to E. wushanense and E. pubescens. In total of 14 homologous fragments totaling 8,954 bp were identified between the cp and mt genomes of E. sagittatum. CONCLUSIONS: This study presents the first assembled and annotated mt genome of E. sagittatum, which provides a valuable genetic resource for the Epimedium genus and lays the foundation for investigating the phylogenetic relationship and genetic variation of this invaluable medicinal plant.

Epimedium↗

Genomic analysis of community-associated multidrug-resistant Klebsiella quasipneumoniae subsp. similipneumoniae and the identification of the ST2059-KL1 clone in the U.S.

UNLABELLED: Klebsiella quasipneumoniae subsp. similipneumoniae is an important member of the K. pneumoniae species complex (KpSC) and is increasingly reported as multidrug-resistant (MDR) in healthcare- and community-associated infections. Since clinical laboratories do not routinely distinguish K. quasipneumoniae subsp. similipneumoniae from K. pneumoniae, national prevalence estimates, particularly for MDR, are lacking. In this study, a total of 2,006 community-associated MDR KpSC isolates were collected from 42 U.S. states, with 30 K. quasipneumoniae subsp. similipneumoniae isolates originating from 12 states identified using whole genome sequencing. All isolates were resistant to ceftriaxone and exhibited high rates of resistance to other antimicrobial agents, including ampicillin-sulbactam (56.7%, 17/30), levofloxacin (75.9%, 22/29), and trimethoprim-sulfamethoxazole (53.3%, 16/30). Notably, five isolates were also carbapenem-resistant. Genomic analysis resolved 10 sequence types (STs), with ST2059 (n = 13) and ST414 (n = 9) predominating. Ceftriaxone resistance in most isolates (90%, 27/30) was conferred by an extended-spectrum β-lactamase gene, predominantly blaCTX-M-15 (73.3%, 22/30); the remaining isolates carried either a carbapenemase (blaKPC-3) or an AmpC β-lactamase (blaCMY-2). Nanopore sequencing identified blaCTX-M-15 harbored on two types of IncFIB(Kpn3) antimicrobial resistance (AMR) plasmids, either with or without the conjugative tra gene cluster. Interestingly, the KL1 locus, associated with canonical hypervirulent K. pneumoniae strains, was detected in all ST2059 isolates. Further analysis of public genomic data showed that the KL1 locus is widely distributed across KpSC. KL1 phylogenetic analyses indicated frequent intrasubspecies recombination but limited intersubspecies exchange of KL1. The identification of the dominant MDR K. quasipneumoniae subsp. similipneumoniae KL1-ST2059 clone in the U.S. underscores the importance of ongoing genomic surveillance. IMPORTANCE: Klebsiella quasipneumoniae subsp. similipneumoniae is an underrecognized member of the Klebsiella pneumoniae species complex that is frequently misidentified in clinical laboratories, leading to an incomplete understanding of its role in antimicrobial resistance. In this study, we used large-scale genomic surveillance of community-associated multidrug-resistant isolates across the U.S. to identify this subspecies as a reservoir of clinically relevant resistance plasmids. Notably, we detected a widely distributed ST2059 lineage carrying the K1 capsular locus, a feature traditionally associated with hypervirulent K. pneumoniae. These findings highlight the convergence of resistance and virulence-associated traits in an overlooked species and underscore the need for genomic surveillance to monitor emerging high-risk lineages in community settings.

Drug Resistance, Multiple, Bacterial↗

Genetic heterogeneity in patients with enlarged vestibular aqueduct and Pendred syndrome.

BACKGROUND: Pathogenic variants in the SLC26A4 gene, encoding for Cl-/HCO3- and I- anion transporter pendrin, are associated with non-syndromic hearing loss with enlarged vestibular aqueduct (NSEVA) and Pendred syndrome (PDS). In the Caucasian population, up to 75% of patients fail to identify a genetic cause through biallelic mutations in the SLC26A4 gene. The CEVA haplotype could therefore play an important role in the diagnostics of NSEVA. The aim of the study was to determine the genetic etiology of hearing loss with EVA or with fully developed PDS in 37 probands and the functional characterization of novel variants identified in the SLC26A4 gene. METHODS: To determine the genetic etiology, Sanger sequencing, WES and KASP genotyping assay were used. Functional characterization of SLC26A4 variants c.140G>A (p.R47Q), c.415G>A (p.G139R), c.441G>A (p.M147I), c.481T>A (p.F161I), c.1589A>C (p.Y530S) and c.2260del (p.D754Ifs*5) involved determination of iodide influx, total and plasma membrane pendrin expression level and subcellular localization of pendrin by confocal imaging. The nanopore sequencing of nasopharyngeal swab samples was performed to confirm the pathogenic effect of potential splice site variant c.415G>A. RESULTS: Biallelic variants in the SLC26A4 gene (M2 genotype) were identified in ten probands and a complete CEVA haplotype was confirmed in three probands harbouring SLC26A4 monoallelic variants (M1 genotype). Fifteen variants in the SLC26A4 gene were identified in total, three of which are novel. The functional characterization of the novel variants and variants which were not yet functionally characterized confirmed the pathogenic potential of five out of six tested variants (p.G139R, p.M147I, p.Y530S, p.D754Ifs*5, and p.F161I). Analysis of nasopharyngeal swab samples confirmed exon 4 skipping due to novel variant SLC26A4:c.415G>A. Probands with biallelic SLC26A4 variants had significantly larger thyroid volume per m2 of body surface area than subjects with monoallelic SLC26A4 variants and the CEVA haplotype. CONCLUSIONS: The genetic aetiology was determined in 13 out of 37 probands (35%), seven manifested with PDS and six with NSEVA. The present study highlights the importance of functional testing to confirm the pathogenicity of SLC26A4 variants and the phenotype-genotype correlation in SLC26A4-related disorders.

Humans↗

Coexistence of carbapenemase and hypervirulence-associated genes among Klebsiella pneumoniae high-risk clones in Hungary.

INTRODUCTION: Strains of Klebsiella pneumoniae carrying hypervirulence and carbapenemase genes represent a rapidly emerging global public health threat. Our study aimed to comprehensively characterise the genomics of hypervirulence-associated and carbapenemase genes carrying K. pneumoniae (hv(a)CpKp) isolates in Hungary. MATERIALS AND METHODS: Between January 2022 and April 2024, 89 aerobactin (iucA-D/iutA)-positive non-duplicate carbapenemase-producing K. pneumoniae isolates from 15 Hungarian healthcare institutes underwent short-read (Illumina, MiSeq, NextSeq) whole-genome sequencing, followed by detailed plasmid analysis using long-read sequencing (Nanopore, MinION) in a representative subset of 32 strains. RESULTS: Most isolates (79/89) belonged to the high-risk clone ST147. Hypervirulence-associated (hva) genes-including rmpA/rmpA2, peg344, shiF, iucA-D, and iutA-were universally present, and 59 isolates possessed chromosomally integrated yersiniabactin loci. Most isolates (87/89) carried the bla NDM-1 carbapenemase gene. Hypervirulence-associated genes were most frequently (29/32) associated with IncHI1B/IncFIB(Mar) plasmids. Notably, we identified plasmids carrying both hva and carbapenemase genes-designated as hybrid plasmids-in 13 of 32 strains. The bla NDM-1 was linked to the IS26 transposase and was present in conserved, identical cassettes on all bla NDM-1-carrying plasmids. DISCUSSION/CONCLUSION: Our study identified hv(a)CpKp strains, particularly the ST147 clone, circulating in Hungary. Our findings highlight the need for routine virulence gene monitoring and continuous genomic and plasmid-based surveillance to mitigate the clinical and epidemiological impact of emerging hv(a)CpKp lineages.

Klebsiella pneumoniae↗

Molecular Evolution and Zoonotic Potential of Muju Virus (Orthohantavirus puumalaense) in Craseomys regulus, Republic of Korea.

Orthohantavirus puumalaense causes hemorrhagic fever with renal syndrome in Europe, with Puumala virus (PUUV) as its primary representative. Muju virus (MUJV), harbored by Craseomys regulus, an Arvicolinae rodent species endemic to the Republic of Korea (ROK), is also a genotype of O. puumalaense. However, their genomic diversity and zoonotic potential remain largely unknown. To investigate their prevalence, 185 voles were collected from 23 regions of the ROK between 2012 and 2023. Serological assays detected anti-PUUV immunoglobulin G antibodies in five samples (3.1%), whereas reverse-transcription polymerase chain reaction confirmed MUJV RNA in identical specimens (2.7%). Amplicon-based nanopore sequencing facilitates near-complete genome recovery, enabling high-resolution comparative analysis. Phylogenetic analysis revealed distinct genetic lineages in Gangwon and Jeollabuk Provinces. Evolutionary rate estimates indicated greater sequence divergence in the S and L segments than in the M segment. A zoonotic risk assessment revealed that most MUJV variants exhibited moderate-to-high spillover potential. The molecular detection of MUJV in Cheorwon, Gangwon Province, expands its known geographic range and provides the first molecular evidence of MUJV circulation in this region. These findings highlight the need for continued surveillance and seroprevalence studies of MUJV to assess its potential for human exposure and public health relevance in the ROK.

Animals↗

Field-isolate recombinant tick-borne encephalitis viruses define reporter-stability guidelines for antiviral testing in flaviviruses.

As arthropod-borne viruses continue to threaten populations globally, there is a pressing need for experimental systems that enable rapid antiviral discovery. Reverse-genetics platforms producing recombinant reporter orthoflaviviruses have been developed to address this gap. Here, we present two new recombinant tick-borne encephalitis viruses (TBEVrec) generated on a European-subtype Haselmühl Tiho1 isolate backbone. A reporter gene, either eGFP or Nluc, was inserted in the capsid-coding region of the genome downstream of the capsid RNA regulatory signal and separated from the complete viral polyprotein by a 2A self-cleaving peptide. TBEVrec was better rescued using the circular polymerase extension reaction (CPER) than with the infectious subgenomic amplicon (ISA) method. TBEVrec replicated efficiently in relevant human cell lines, with comparable replication to wild-type TBEV in a neuronal cell line and moderately reduced titers and RNA levels in immune-derived cell lines. Using either eGFP or Nluc, we illustrate how TBEVrec enabled high-content RNAi screening, highlighting Nucleolin and PRKD1 as potential TBEV host factors, and drug testing on a benchtop plate reader. Nanopore sequencing of the eGFP insert revealed that the reporter is excised without affecting flanking regions. Comparative analysis of eGFP and Nluc further shows that this instability is time- and cell type-dependent, and that Nluc is comparatively more stable. From these observations, we outline safeguards and design principles that are broadly applicable both to the rescue of existing constructs and to the design of future recombinant reporter virus platforms.

CPER↗

Protocol to improve isoform-level quantification of low-abundance transcripts via STALARD pre-amplification.

STALARD (selective target amplification for low-abundance RNA detection) enables isoform-level quantification of low-abundance RNAs using conventional laboratory equipment. Here, we describe steps for RNA isolation, primer design, reverse transcription, selective target amplification, and downstream analysis. The protocol couples selective pre-amplification with a quantitative reverse-transcription PCR (RT-qPCR) readout and optional nanopore sequencing. Using 1 μg input RNA and 12 pre-amplification cycles, STALARD reduces Cq values by approximately 10-12 cycles, bringing the target into a reliably quantifiable range. For complete details on the use and execution of this protocol, please refer to Jeong et al.1.

Gene Expression↗

In vitro reconstitution of chromatin replication recapitulates symmetric histone recycling.

Symmetric histone recycling is vital for maintaining epigenetic inheritance upon eukaryotic DNA replication. Recent genome-wide studies have uncovered key determinants of this process, but how these factors collectively support parental histone transfer remains incompletely understood. Here, we successfully reconstitute histone recycling with 24 purified proteins and analyze the products digested by Micrococcal nuclease with Repli-pore-seq, the newly developed pipeline combining nanopore sequencing and deep-learning-based classification. As a result, we identify histones symmetrically recycled as tetrasomes or hexasomes on nucleosome-favorable sequences. We also observe the discordance of the recycled position between lagging and leading strands on the GC-rich DNA sequences. Moreover, removal of Pol δ, Pol32, Dpb3/4, Ctf4, Csm3/Tof1, or Mrc1 disrupts the balance of histone recycling between the two daughter strands, whereas removal of Ctf4, Csm3/Tof1, or Mrc1 additionally alters the positions at which histones were recycled. Furthermore, addition of the lagging-strand maturation factors Fen1 and Cdc9 enhances histone recycling to the lagging strand. These findings provide critical insights into the molecular players and mechanisms underlying symmetric histone recycling.

Histones↗

Fluctuating DNA methylation tracks cancer evolution at clinical scale.

Cancer development and response to treatment are evolutionary processes1,2, but characterizing evolutionary dynamics at a clinically meaningful scale has remained challenging3. Here we develop a new methodology called EVOFLUx, based on natural DNA methylation barcodes fluctuating over time4, that quantitatively infers evolutionary dynamics using only a bulk tumour methylation profile as input. We apply EVOFLUx to 1,976 well-characterized lymphoid cancer samples spanning a broad spectrum of diseases and show that initial tumour growth rate, malignancy age and epimutation rates vary by orders of magnitude across disease types. We measure that subclonal selection occurs only infrequently within bulk samples and detect occasional examples of multiple independent primary tumours. Clinically, we observe faster initial tumour growth in more aggressive disease subtypes, and that evolutionary histories are strong independent prognostic factors in two series of chronic lymphocytic leukaemia. Using EVOFLUx for phylogenetic analyses of aggressive Richter-transformed chronic lymphocytic leukaemia samples detected that the seed of the transformed clone existed decades before presentation. Orthogonal verification of EVOFLUx inferences is provided using additional genetic data, including long-read nanopore sequencing, and clinical variables. Collectively, we show how widely available, low-cost bulk DNA methylation data precisely measure cancer evolutionary dynamics, and provides new insights into cancer biology and clinical behaviour.

Humans↗