KAS-seq: genome-wide sequencing of single-stranded DNA by N3-kethoxal-assisted labeling

Genome-wide characterization of single-stranded DNA in rice

Single-stranded DNA (ssDNA) is essential for various DNA-templated processes in both eukaryotes and prokaryotes. However, comprehensive characterizations of ssDNA still lag in plants compared to nonplant systems. Here, we conducted in situ S1-sequencing, with starting gDNA ranging from 5 µg to 250 ng, followed by comprehensive characterizations of ssDNA in rice (Oryza sativa L.). We found that ssDNA loci were substantially associated with a subset of non-B DNA structures and functional genomic loci. Subtypes of ssDNA loci had distinct epigenetic features. Importantly, ssDNA may act alone or partly coordinate with non-B DNA structures, functional genomic loci, or epigenetic marks to actively or repressively modulate gene transcription, which is genomic region dependent and associated with the distinct accumulation of RNA Pol II. Moreover, distinct types of ssDNA had differential impacts on the activities and evolution of transposable elements (TEs) (especially common or conserved TEs) in the rice genome. Our study showcases an antibody-independent technique for characterizing non-B DNA structures or functional genomic loci in plants. It lays the groundwork and fills a crucial gap for further exploration of ssDNA, non-B DNA structures, or functional genomic loci, thereby advancing our understanding of their biology in plants.

PTPN2 copper-sensing relays copper level fluctuations into EGFR/CREB activation and associated CTR1 transcriptional repression

Fluxes in human copper levels recently garnered attention for roles in cellular signaling, including affecting levels of the signaling molecule cyclic adenosine monophosphate. We herein apply an unbiased temporal evaluation of the signaling and whole genome transcriptional activities modulated by copper level fluctuations to identify potential copper sensor proteins responsible for driving these activities. We find that fluctuations in physiologically relevant copper levels modulate EGFR signal transduction and activation of the transcription factor CREB. Both intracellular and extracellular assays support Cu1+ inhibition of the EGFR phosphatase PTPN2 (and potentially PTPN1)-via ligation to the PTPN2 active site cysteine side chain-as the underlying mechanism. We additionally show i) copper supplementation drives weak transcriptional repression of the copper importer CTR1 and ii) CREB activity is inversely correlated with CTR1 expression. In summary, our study reveals PTPN2 as a physiological copper sensor and defines a regulatory mechanism linking feedback control of copper stimulated EGFR/CREB signaling and CTR1 expression.

Quantitative analysis of cis-regulatory elements in transcription with KAS-ATAC-seq

Cis-regulatory elements (CREs) are pivotal in orchestrating gene expression throughout diverse biological systems. Accurate identification and in-depth characterization of functional CREs are crucial for decoding gene regulation networks during cellular processes. In this study, we develop Kethoxal-Assisted Single-stranded DNA Assay for Transposase-Accessible Chromatin with Sequencing (KAS-ATAC-seq) to quantitatively analyze the transcriptional activity of CREs. A main advantage of KAS-ATAC-seq lies in its precise measurement of ssDNA levels within both proximal and distal ATAC-seq peaks, enabling the identification of transcriptional regulatory sequences. This feature is particularly adept at defining Single-Stranded Transcribing Enhancers (SSTEs). SSTEs are highly enriched with nascent RNAs and specific transcription factors (TFs) binding sites that define cellular identity. Moreover, KAS-ATAC-seq provides a detailed characterization and functional implications of various SSTE subtypes. Our analysis of CREs during mouse neural differentiation demonstrates that KAS-ATAC-seq can effectively identify immediate-early activated CREs in response to retinoic acid (RA) treatment. Our findings indicate that KAS-ATAC-seq provides more precise annotation of functional CREs in transcription. Future applications of KAS-ATAC-seq would help elucidate the intricate dynamics of gene regulation in diverse biological processes.

Cytoskeletal gene alterations linked to sorafenib resistance in hepatocellular carcinoma

BACKGROUND

Although sorafenib has been consistently used as a first-line treatment for advanced hepatocellular carcinoma (HCC), most patients will develop resistance, and the mechanism of resistance to sorafenib needs further study.

METHODS

Using KAS-seq technology, we obtained the ssDNA profiles within the whole genome range of SMMC-7721 cells treated with sorafenib for differential analysis. We then intersected the differential genes obtained from the analysis of hepatocellular carcinoma patients in GSE109211 who were ineffective and effective with sorafenib treatment, constructed a PPI network, and obtained hub genes. We then analyzed the relationship between the expression of these genes and the prognosis of hepatocellular carcinoma patients.

RESULTS

In this study, we identified 7 hub ERGs (ACTB, CFL1, ACTG1, ACTN1, WDR1, TAGLN2, HSPA8) related to drug resistance, and these genes are associated with the cytoskeleton.

CONCLUSIONS

The cytoskeleton is associated with sorafenib resistance in hepatocellular carcinoma. Using KAS-seq to analyze the early changes in tumor cells treated with drugs is feasible for studying the drug resistance of tumors, which provides reference significance for future research.

Exploring transient global transcriptional changes induced by ascorbic acid revealed via atKAS-seq profiling

Transcription initiates the formation of single-stranded DNA (ssDNA) regions within the genome, delineating transcription bubbles, a highly dynamic genomic process. Kethoxal-assisted single-stranded DNA sequencing (KAS-seq) utilizing N3-kethoxal has emerged as a potent tool for mapping specific guanine positions in ssDNA on a genome-wide scale. However, the original KAS-seq method required the costly Accel-NGS Methyl-seq DNA library kit. This study introduces an optimized iteration of the KAS-seq technique, referred to as adapter-tagged KAS-seq (atKAS-seq), incorporating an adapter tagging strategy. This modification involves integrating sequencing adapters via complementary strand synthesis using random N9 tagging. Additionally, by harnessing the potential of ascorbic acid (ASC), recognized for inducing global epigenetic changes, we employed the atKAS-seq methodology to elucidate critical pathways influenced by short-term, high-dose ASC treatment. Our findings underscore that atKAS-seq enables rapid and precise analyses of transcription dynamics and enhancer activities concurrently. This method offers a streamlined, cost-efficient, and low-input approach, affirming its utility in probing intricate genomic regulatory mechanisms.

G-quadruplex DNA structure is a positive regulator of MYC transcription

DNA structure can regulate genome function. Four-stranded DNA G-quadruplex (G4) structures have been implicated in transcriptional regulation; however, previous studies have not directly addressed the role of an individual G4 within its endogenous cellular context. Using CRISPR to genetically abrogate endogenous G4 structure folding, we directly interrogate the G4 found within the upstream regulatory region of the critical human MYC oncogene. G4 loss leads to suppression of MYC transcription from the P1 promoter that is mediated by the deposition of a de novo nucleosome alongside alterations in RNA polymerase recruitment. We also show that replacement of the endogenous MYC G4 with a different G4 structure from the KRAS oncogene restores G4 folding and MYC transcription. Moreover, we demonstrate that the MYC G4 structure itself, rather than its sequence, recruits transcription factors and histone modifiers. Overall, our work establishes that G4 structures are important features of transcriptional regulation that coordinate recruitment of key chromatin proteins and the transcriptional machinery through interactions with DNA secondary structure, rather than primary sequence.

Trans-vaccenic acid reprograms CD8+ T cells and anti-tumour immunity

Diet-derived nutrients are inextricably linked to human physiology by providing energy and biosynthetic building blocks and by functioning as regulatory molecules. However, the mechanisms by which circulating nutrients in the human body influence specific physiological processes remain largely unknown. Here we use a blood nutrient compound library-based screening approach to demonstrate that dietary trans-vaccenic acid (TVA) directly promotes effector CD8+ T cell function and anti-tumour immunity in vivo. TVA is the predominant form of trans-fatty acids enriched in human milk, but the human body cannot produce TVA endogenously1. Circulating TVA in humans is mainly from ruminant-derived foods including beef, lamb and dairy products such as milk and butter2,3, but only around 19% or 12% of dietary TVA is converted to rumenic acid by humans or mice, respectively4,5. Mechanistically, TVA inactivates the cell-surface receptor GPR43, an immunomodulatory G protein-coupled receptor activated by its short-chain fatty acid ligands6-8. TVA thus antagonizes the short-chain fatty acid agonists of GPR43, leading to activation of the cAMP-PKA-CREB axis for enhanced CD8+ T cell function. These findings reveal that diet-derived TVA represents a mechanism for host-extrinsic reprogramming of CD8+ T cells as opposed to the intrahost gut microbiota-derived short-chain fatty acids. TVA thus has translational potential for the treatment of tumours.

Detection of alternative DNA structures and its implications for human disease

Around 3% of the genome consists of simple DNA repeats that are prone to forming alternative (non-B) DNA structures, such as hairpins, cruciforms, triplexes (H-DNA), four-stranded guanine quadruplexes (G4-DNA), and others, as well as composite RNA:DNA structures (e.g., R-loops, G-loops, and H-loops). These DNA structures are dynamic and favored by the unwinding of duplex DNA. For many years, the association of alternative DNA structures with genome function was limited by the lack of methods to detect them in vivo. Here, we review the recent advancements in the field and present state-of-the-art technologies and methods to study alternative DNA structures. We discuss the limitations of these methods as well as how they are beginning to provide insights into causal relationships between alternative DNA structures, genome function and stability, and human disease.

KAS-Analyzer: a novel computational framework for exploring KAS-seq data

Motivation

Kethoxal-assisted ssDNA sequencing (KAS-seq) is rapidly gaining popularity as a robust and effective approach to study the nascent dynamics of transcriptionally engaged RNA polymerases through profiling of genome-wide single-stranded DNA (ssDNA). Its latest variant, spKAS-seq, a strand-specific version of KAS-seq, has been developed to map genome-wide R-loop structures by detecting imbalances of ssDNA on two strands. However, user-friendly, open-source computational tools tailored for KAS-seq data are still lacking.

Results

Here, we introduce KAS-Analyzer, the first comprehensive computational framework aimed at streamlining and enhancing the analysis and interpretation of KAS-seq and spKAS-seq data. In addition to standard analyses, KAS-Analyzer offers many novel tools specifically designed for KAS-seq data, including, but not limited to: calculation of transcription-related metrics, identification of single-stranded transcribing (SST) enhancers, high-resolution mapping of R-loops, and differential RNA polymerase activity analysis. We provided a detailed overview of KAS-seq data and its diverse applications through the implementation of KAS-Analyzer. Using the example time-course KAS-seq datasets, we further showcase the robust capabilities of KAS-Analyzer for investigating dynamic transcriptional regulatory programs in response to UVB radiation.

Availability and implementation

KAS-Analyzer is available at https://github.com/Ruitulyu/KAS-Analyzer.

Adrenergic and mesenchymal signatures are identifiable in cell-free DNA and correlate with metastatic disease burden in children with neuroblastoma

Background

Cell free DNA (cfDNA) profiles of 5-hydroxymethylcytosine (5-hmC), an epigenetic marker of open chromatin and active gene expression, are correlated with metastatic disease burden in patients with neuroblastoma. Neuroblastoma tumors are comprised of adrenergic (ADRN) and mesenchymal (MES) cells, and the relative abundance of each in tumor biopsies has prognostic implications. We hypothesized that ADRN and MES specific signatures could be quantified in cfDNA 5-hmC profiles and would augment the detection of metastatic burden in patients with neuroblastoma.

Methods

We previously performed an integrative analysis to identify ADRN and MES specific genes (n=373 and n=159, respectively). Purified DNA from cell lines was serial diluted with healthy donor cfDNA. Using Gene Set Variation Analysis (GSVA), ADRN and MES signatures were optimized. We then quantified signature scores, and our prior neuroblastoma signature, in cfDNA from 84 samples from 46 high-risk patients including 21 patients with serial samples.

Results

Samples from patients with higher metastatic burden had increased GSVA scores for both ADRN and MES gene signatures (p < 0.001). While ADRN and MES signature scores tracked together in serially collected samples, we identified instances of patients with increases in either MES or ADRN score at relapse.

Conclusions

While it is feasible to identify ADRN and MES signatures using 5-hmC profiles of cfDNA from neuroblastoma patients and correlate these signatures to metastatic burden, additional data are needed to determine the optimal strategies for clinical implementation. Prospective evaluation in larger cohorts is ongoing.

R-loop-dependent promoter-proximal termination ensures genome stability

The proper regulation of transcription is essential for maintaining genome integrity and executing other downstream cellular functions1,2. Here we identify a stable association between the genome-stability regulator sensor of single-stranded DNA (SOSS)3 and the transcription regulator Integrator-PP2A (INTAC)4-6. Through SSB1-mediated recognition of single-stranded DNA, SOSS-INTAC stimulates promoter-proximal termination of transcription and attenuates R-loops associated with paused RNA polymerase II to prevent R-loop-induced genome instability. SOSS-INTAC-dependent attenuation of R-loops is enhanced by the ability of SSB1 to form liquid-like condensates. Deletion of NABP2 (encoding SSB1) or introduction of cancer-associated mutations into its intrinsically disordered region leads to a pervasive accumulation of R-loops, highlighting a genome surveillance function of SOSS-INTAC that enables timely termination of transcription at promoters to constrain R-loop accumulation and ensure genome stability.

PTPN2 copper-sensing rapidly relays copper level fluctuations into EGFR/CREB activation and associated CTR1 transcriptional repression

Fluxes in human intra- and extracellular copper levels recently garnered attention for roles in cellular signaling, including affecting levels of the signaling molecule cyclic adenosine monophosphate (cAMP). We herein applied an unbiased temporal evaluation of the whole-genome transcriptional activities modulated by fluctuations in copper levels to identify the copper sensor proteins responsible for driving these activities. We found that fluctuations in physiologically-relevant copper levels rapidly modulate EGFR/MAPK/ERK signal transduction and activation of the transcription factor cAMP response element-binding protein (CREB). Both intracellular and extracellular assays support Cu 1+ inhibition of the EGFR-phosphatase PTPN2 (and potentially the homologous PTPN1)-via direct ligation to the PTPN2 active site cysteine side chain-as the underlying mechanism of copper-stimulated EGFR signal transduction activation. Depletion of copper represses this signaling pathway. We additionally show i ) copper supplementation drives transcriptional repression of the copper importer CTR1 and ii ) CREB activity is inversely correlated with CTR1 expression. In summary, our study reveals PTPN2 as a physiological copper sensor and defines a regulatory mechanism linking feedback control of copper-stimulated MAPK/ERK/CREB-signaling and CTR1 expression, thereby uncovering a previously unrecognized link between copper levels and cellular signal transduction.

5-hydroxymethylcytosine profiling of cell-free DNA identifies bivalent genes that are prognostic of survival in high-risk neuroblastoma

Neuroblastoma is the most common extra-cranial solid tumor in childhood and epigenetic dysregulation is a key driver of this embryonal disease. In cell-free DNA from neuroblastoma patients with high-risk disease, we found increased 5-hydroxymethylcytosine (5-hmC) deposition on Polycomb Repressive Complex 2 (PRC2) target genes, a finding previously described in the context of bivalent genes. As bivalent genes, defined as genes bearing both activating (H3K4me3) and repressive (H3K27me3) chromatin modifications, have been shown to play an important role in development and cancer, we investigated the potential role of bivalent genes in maintaining a de-differentiated state in neuroblastoma and their potential use as a biomarker. We identified 313 genes that bore bivalent chromatin marks, were enriched for mediators of neuronal differentiation, and were transcriptionally repressed across a panel of heterogenous neuroblastoma cell lines. Through gene set variance analysis, we developed a clinically implementable bivalent signature. In three distinct clinical cohorts, low bivalent signature was significantly and independently associated with worse clinical outcome in high-risk neuroblastoma patients. Thus, low expression of bivalent genes is a biomarker of ultra-high-risk disease and may represent a therapeutic opportunity in neuroblastoma.

Noncanonical DNA structures are drivers of genome evolution

In addition to the canonical right-handed double helix, other DNA structures, termed 'non-B DNA', can form in the genomes across the tree of life. Non-B DNA regulates multiple cellular processes, including replication and transcription, yet its presence is associated with elevated mutagenicity and genome instability. These discordant cellular roles fuel the enormous potential of non-B DNA to drive genomic and phenotypic evolution. Here we discuss recent studies establishing non-B DNA structures as novel functional elements subject to natural selection, affecting evolution of transposable elements (TEs), and specifying centromeres. By highlighting the contributions of non-B DNA to repeated evolution and adaptation to changing environments, we conclude that evolutionary analyses should include a perspective of not only DNA sequence, but also its structure.

spKAS-seq reveals R-loop dynamics using low-input materials by detecting single-stranded DNA with strand specificity

R-loops affect transcription and genome stability. Dysregulation of R-loops is related to human diseases. Genome-wide R-loop mapping typically uses the S9.6 antibody or inactive ribonuclease H, both requiring a large number of cells with varying results observed depending on the approach applied. Here, we present strand-specific kethoxal-assisted single-stranded DNA (ssDNA) sequencing (spKAS-seq) to map R-loops by taking advantage of the presence of a ssDNA in the triplex structure. We show that spKAS-seq detects R-loops and their dynamics at coding sequences, enhancers, and other intergenic regions with as few as 50,000 cells. A joint analysis of R-loops and chromatin-bound RNA binding proteins (RBPs) suggested that R-loops can be RBP binding hotspots on the chromatin.

Replication dependent and independent mechanisms of GAA repeat instability

Trinucleotide repeat instability is a driver of human disease. Large expansions of (GAA)n repeats in the first intron of the FXN gene are the cause Friedreich's ataxia (FRDA), a progressive degenerative disorder which cannot yet be prevented or treated. (GAA)n repeat instability arises during both replication-dependent processes, such as cell division and intergenerational transmission, as well as in terminally differentiated somatic tissues. Here, we provide a brief historical overview on the discovery of (GAA)n repeat expansions and their association to FRDA, followed by recent advances in the identification of triplex H-DNA formation and replication fork stalling. The main body of this review focuses on the last decade of progress in understanding the mechanism of (GAA)n repeat instability during DNA replication and/or DNA repair. We propose that the discovery of additional mechanisms of (GAA)n repeat instability can be achieved via both comparative approaches to other repeat expansion diseases and genome-wide association studies. Finally, we discuss the advances towards FRDA prevention or amelioration that specifically target (GAA)n repeat expansions.