Single-base Resolution Research Articles

Single-Base Resolution Detection of N6-Methyladenosine in RNA by Adenosine Deamination Sequencing.

N6-Methyladenosine (m6A) is one of the most abundant and prevalent natural modifications occurring in diverse RNA species. m6A plays a wide range of roles in physiological and pathological processes. Revealing the functions of m6A relies on the faithful detection of individual m6A sites in RNA. However, developing a simple method for the single-base resolution detection of m6A is still a challenging task. Herein, we report an adenosine deamination sequencing (AD-seq) technique for the facile detection of m6A in RNA at single-base resolution. The AD-seq approach capitalizes on the selective deamination of adenosine, but not m6A, by the evolved tRNA adenosine deaminase (TadA) variant of TadA8e or the dimer protein of TadA-TadA8e. In AD-seq, adenosine is deaminated by TadA8e or TadA-TadA8e to form inosine, which pairs with cytidine and is read as guanosine in sequencing. m6A resists deamination due to the interference of the methyl group at the N6 position of adenosine. Thus, the m6A base pairs with thymine and is still read as adenosine in sequencing. The differential readouts from A and m6A in sequencing can achieve the single-base resolution detection of m6A in RNA. Application of the proposed AD-seq successfully identified individual m6A sites in Escherichia coli 23S rRNA. Taken together, the proposed AD-seq allows simple and cost-effective detection of m6A at single-base resolution in RNA, which provides a valuable tool to decipher the functions of m6A in RNA.

Precise characterization of somatic complex structural variations from tumor/control paired long-read sequencing data with nanomonsv.

We present our novel software, nanomonsv, for detecting somatic structural variations (SVs) using tumor and matched control long-read sequencing data with a single-base resolution. The current version of nanomonsv includes two detection modules, Canonical SV module, and Single breakend SV module. Using tumor/control paired long-read sequencing data from three cancer and their matched lymphoblastoid lines, we demonstrate that Canonical SV module can identify somatic SVs that can be captured by short-read technologies with higher precision and recall than existing methods. In addition, we have developed a workflow to classify mobile element insertions while elucidating their in-depth properties, such as 5' truncations, internal inversions, as well as source sites for 3' transductions. Furthermore, Single breakend SV module enables the detection of complex SVs that can only be identified by long-reads, such as SVs involving highly-repetitive centromeric sequences, and LINE1- and virus-mediated rearrangements. In summary, our approaches applied to cancer long-read sequencing data can reveal various features of somatic SVs and will lead to a better understanding of mutational processes and functional consequences of somatic SVs.

Genome-wide DNA methylation sequencing identifies epigenetic perturbations in the upper airways under long-term exposure to moderate levels of ambient air pollution

While the link between exposure to high levels of ambient particulate matter (PM) and increased incidences of respiratory and cardiovascular diseases is widely recognized, recent epidemiological studies have shown that low PM concentrations are equally associated with adverse health effects. As DNA methylation is one of the main mechanisms by which cells regulate and stabilize gene expression, changes in the methylome could constitute early indicators of dysregulated signaling pathways. So far, little is known about PM-associated DNA methylation changes in the upper airways, the first point of contact between airborne pollutants and the human body. Here, we focused on cells of the upper respiratory tract and assessed their genome-wide DNA methylation pattern to explore exposure-associated early regulatory changes. Using a mobile epidemiological laboratory, nasal lavage samples were collected from a cohort of 60 adults that lived in districts with records of low (Simmerath) or moderate (Stuttgart) PM10 levels in Germany. PM10 concentrations were verified by particle measurements on the days of the sample collection and genome-wide DNA methylation was determined by enzymatic methyl sequencing at single-base resolution. We identified 231 differentially methylated regions (DMRs) between moderately and lowly PM10 exposed individuals. A high proportion of DMRs overlapped with regulatory elements, and DMR target genes were involved in pathways regulating cellular redox homeostasis and immune response. In addition, we found distinct changes in DNA methylation of the HOXA gene cluster whose methylation levels have previously been linked to air pollution exposure but also to carcinogenesis in several instances. The findings of this study suggest that regulatory changes in upper airway cells occur at PM10 levels below current European thresholds, some of which may be involved in the development of air pollution-related diseases.

Direct enzymatic sequencing of 5-methylcytosine at single-base resolution.

5-methylcytosine (5mC) is the most important DNA modification in mammalian genomes. The ideal method for 5mC localization would be both nondestructive of DNA and direct, without requiring inference based on detection of unmodified cytosines. Here we present direct methylation sequencing (DM-Seq), a bisulfite-free method for profiling 5mC at single-base resolution using nanogram quantities of DNA. DM-Seq employs two key DNA-modifying enzymes: a neomorphic DNA methyltransferase and a DNA deaminase capable of precise discrimination between cytosine modification states. Coupling these activities with deaminase-resistant adapters enables accurate detection of only 5mC via a C-to-T transition in sequencing. By comparison, we uncover a PCR-related underdetection bias with the hybrid enzymatic-chemical TET-assisted pyridine borane sequencing approach. Importantly, we show that DM-Seq, unlike bisulfite sequencing, unmasks prognostically important CpGs in a clinical tumor sample by not confounding 5mC with 5-hydroxymethylcytosine. DM-Seq thus offers an all-enzymatic, nondestructive, faithful and direct method for the reading of 5mC alone.

Characteristics of Salvia miltiorrhiza methylome and the regulatory mechanism of DNA methylation in tanshinone biosynthesis.

Salvia miltiorrhiza is a model medicinal plant with significant economic and medicinal value. Its roots produce a group of diterpenoid lipophilic bioactive components, termed tanshinones. Biosynthesis and regulation of tanshinones has attracted widespread interest. However, the methylome of S. miltiorrhiza has not been analysed and the regulatory mechanism of DNA methylation in tanshinone production is largely unknown. Here we report single-base resolution DNA methylomes from roots and leaves. Comparative analysis revealed differential methylation patterns for CG, CHG, and CHH contexts and the association between DNA methylation and the expression of genes and small RNAs. Lowly methylated genes always had higher expression levels and 24-nucleotide sRNAs could be key players in the RdDM pathway in S. miltiorrhiza. DNA methylation variation analysis showed that CHH methylation contributed mostly to the difference. Go enrichment analysis showed that diterpenoid biosynthetic process was significantly enriched for genes with downstream overlapping with hypoCHHDMR in July_root when comparing with those in March_root. Tanshinone biosynthesis-related enzyme genes, such as DXS2, CMK, IDI1, HMGR2, DXR, MDS, CYP76AH1, 2OGD25, and CYP71D373, were less CHH methylated in gene promoters or downstream regions in roots collected in July than those collected in March. Consistently, gene expression was up-regulated in S. miltiorrhiza roots collected in July compared with March and the treatment of DNA methylation inhibitor 5-azacytidine significantly promoted tanshinone production. It suggests that DNA methylation plays a significant regulatory role in tanshinone biosynthesis in S. miltiorrhiza through changing the levels of CHH methylation in promoters or downstreams of key enzyme genes.

The performance of whole genome bisulfite sequencing on DNBSEQ-Tx platform examined by different library preparation strategies

BackgroundWhole-genome bisulfite sequencing (WGBS) technology can provide comprehensive DNA methylation at a single-base resolution on a genome-wide scale, and is considered to be the gold standard for the detection of 5-methylcytosine (5 mC). However, the International Human Epigenome Consortium propose a full DNA methylome should have at least 30 fold redundant coverage of the reference genome from a single biological replicate. Therefore, it remains cost prohibitive for large-scale studies. To find a solution, the DNBSEQ-Tx sequencing was developed that can generate up to 6 Tb data in a single run for projects involving large-scale sequencing. ResultsIn this study, we provided two WGBS library construction methods DNB_PREBSseq and DNB_SPLATseq optimized for the DNBSEQ-Tx sequencer, and demonstrated the performance of these two methods on the DNBSEQ-Tx platform, using the DNA extracted from four different cell lines. We also compared the sequencing data from these two WGBS library construction methods with HeLa cell line data from ENCODE sequenced on Illumina HiSeq X Ten and WGBS data of two other cell lines sequenced on HiSeq2500. Various quality control (QC) analyses such as the base quality scores, methylation-bias (m-bias), and conversion efficiency indicated that the data sequenced on the DNBSEQ-Tx platform met the WGBS-required quality controls. Meanwhile, our data closely resembled the coverage shown by the data generated by the Illumina platform. ConclusionsOur study showed that with our optimized methods, DNBSEQ-Tx could generate high-quality WGBS data with relatively good stability for large-scale WGBS sequencing applications. Thus, we conclude that DNBSEQ-Tx can be used for a wide range of WGBS research.

Understanding and modeling human traits and diseases: Insights from the comparative genomics resources of Zoonomia

Understanding and modeling human traits and diseases: Insights from the comparative genomics resources of Zoonomia

A Photoredox Reaction for the Selective Modification of 5-Carboxycytosine in DNA.

Covalent epigenetic modifications contribute to the regulation of important cellular processes during development and differentiation, and changes in their genomic distribution and frequency are linked to the emergence of genetic disease states. Chemical and enzymatic methods that selectively target the orthogonal chemical functionality of epigenetic markers are central to the study of their distribution and function, and considerable research effort has been focused on the development of nondestructive sequencing approaches which preserve valuable DNA samples. Photoredox catalysis enables transformations with tunable chemoselectivity under mild, biocompatible reaction conditions. We report the reductive decarboxylation of 5-carboxycytosine via a novel iridium-based treatment, which represents the first application of visible-light photochemistry to epigenetic sequencing via direct base conversion. We propose that the reaction involves an oxidative quenching cycle beginning with single-electron reduction of the nucleobase by the photocatalyst, followed by hydrogen atom transfer from a thiol. The saturation of the C5-C6 backbone permits decarboxylation of the nonaromatic intermediate, and hydrolysis of the N4-amine constitutes a conversion from a cytosine derivative to a T-like base. This conversion demonstrates selectivity for 5-carboxycytosine over other canonical or modified nucleoside monomers, and is thereby applied to the sequencing of 5-carboxycytosine within modified oligonucleotides. The photochemistry explored in this study can also be used in conjunction with enzymatic oxidation by TET to profile 5-methylcytosine at single-base resolution. Compared to other base-conversion treatments, the rapid photochemical reaction takes place within minutes, which could provide advantages for high-throughput detection and diagnostic applications.

Partial erosion on under-methylated regions and chromatin reprogramming contribute to oncogene activation in IDH mutant gliomas

BackgroundIDH1/2 hotspot mutations are well known to drive oncogenic mutations in gliomas and are well-defined in the WHO 2021 classification of central nervous system tumors. Specifically, IDH mutations lead to aberrant hypermethylation of under-methylated regions (UMRs) in normal tissues through the disruption of TET enzymes. However, the chromatin reprogramming and transcriptional changes induced by IDH-related hypermethylation in gliomas remain unclear.ResultsHere, we have developed a precise computational framework based on Hidden Markov Model to identify altered methylation states of UMRs at single-base resolution. By applying this framework to whole-genome bisulfite sequencing data from 75 normal brain tissues and 15 IDH mutant glioma tissues, we identified two distinct types of hypermethylated UMRs in IDH mutant gliomas. We named them partially hypermethylated UMRs (phUMRs) and fully hypermethylated UMRs (fhUMRs), respectively. We found that the phUMRs and fhUMRs exhibit distinct genomic features and chromatin states. Genes related to fhUMRs were more likely to be repressed in IDH mutant gliomas. In contrast, genes related to phUMRs were prone to be up-regulated in IDH mutant gliomas. Such activation of phUMR genes is associated with the accumulation of active H3K4me3 and the loss of H3K27me3, as well as H3K36me3 accumulation in gene bodies to maintain gene expression stability. In summary, partial erosion on UMRs was accompanied by locus-specific changes in key chromatin marks, which may contribute to oncogene activation.ConclusionsOur study provides a computational strategy for precise decoding of methylation encroachment patterns in IDH mutant gliomas, revealing potential mechanistic insights into chromatin reprogramming that contribute to oncogenesis.

Novel molecular requirements for CRISPR RNA-guided transposition.

CRISPR-associated transposases(CASTs) direct DNA integration downstream of target sites using the RNA-guided DNA binding activity of nuclease-deficient CRISPR-Cas systems. Transposition relies on several key protein-protein and protein-DNA interactions, but little is known about the explicit sequence requirements governing efficient transposon DNA integration activity. Here, we exploit pooled library screening and high-throughput sequencing to reveal novel sequence determinants during transposition by the Type I-F Vibrio cholerae CAST system(VchCAST). On the donor DNA, large transposon end libraries revealed binding site nucleotide preferences for the TnsB transposase, as well as an additional conserved region that encoded a consensus binding site for integration host factor (IHF). Remarkably, we found that VchCAST requires IHF for efficient transposition, thus revealing a novel cellular factor involved in CRISPR-associated transpososome assembly. On the target DNA, we uncovered preferred sequence motifs at the integration site that explained previously observed heterogeneity with single-base pair resolution. Finally, we exploited our library data to design modified transposon variants that enable in-frame protein tagging. Collectively, our results provide new clues about the assembly and architecture of the paired-end complex formed between TnsB and the transposon DNA, and inform the design of custom payload sequences for genome engineering applications with CAST systems.

The methylation landscape of giga-genome and the epigenetic timer of age in Chinese pine

Epigenetics has been revealed to play a crucial role in the long-term memory in plants. However, little is known about whether the epigenetic modifications occur with age progressively in conifers. Here, we present the single-base resolution DNA methylation landscapes of the 25-gigabase Chinese pine (Pinus tabuliformis) genome at different ages. The result shows that DNA methylation is closely coupled with the regulation of gene transcription. The age-dependent methylation profile with a linearly increasing trend is the most significant pattern of DMRs between ages. Two segments at the five-prime end of the first ultra-long intron in DAL1, a conservative age biomarker in conifers, shows a gradual decline of CHG methylation as the age increased, which is highly correlated with its expression profile. Similar high correlation is also observed in nine other age marker genes. Our results suggest that DNA methylation serves as an important epigenetic signature of developmental age in conifers.

Abstract 6018: Sequencing of 5-hydroxymethylcytosine to single base resolution

Abstract DNA methylation is an epigenetic regulator of gene expression with important functions in development and diseases such as cancer. The modified cytosines, 5-methylcytosine (5mC) and 5-hydroxymethylcytosine (5hmC) are identified using Illumina sequencing. Libraries are generated using either NEBNext® EM-seq™, an enzyme-based workflow, or by using sodium bisulfite conversion. These methods, however, cannot differentiate between 5mC and 5hmC. Distinguishing these modifications is important as there is increased interest in 5hmCs role in regulating gene expression. Methods currently exist to enable discrimination of 5mC and 5hmC, for example, oxBS-seq and TAB-seq. These methods still rely upon sodium bisulfite conversion and have reduced data quality due to fragmentation and loss of DNA. Here we describe an enzymatic method that enables specific detection of 5hmC. 5hmC is detected using two enzymatic steps. Firstly, 5hmCs are glucosylated, which protects them from subsequent deamination by APOBEC. In contrast, cytosines and 5mCs are deaminated resulting in their conversion to uracil and thymine, respectively. During Illumina sequencing deaminated cytosine and 5mC are read as thymine and the glucosylated 5hmC are read as cytosine. Subtractive analysis of 5hmC data from EM-seq data, which detects both 5mC and 5hmC, permits the identification of individual 5mC and 5hmC sites. 5hmC data were generated for inputs of 0.1 ng to 200 ng DNA isolated from E14 mouse embryonic stem cells and human brain. The 5hmC libraries had similar characteristics to EM-seq libraries, including expected insert sizes due to intact DNA molecules, low duplication rates and minimal GC bias. T4147 phage DNA was used as an internal control, as all cytosines are 5-hydroxymethylated, with 98-99% of cytosines identified as 5hmC. 5mC and 5hmC levels were also profiled during E14 cell differentiation for 10 days. Interestingly, 5hmC levels decreased whereas 5mC levels increased particularly during the first five days of differentiation. LC-MS/MS quantification of this same DNA mirrored the changes observed by sequencing. The ability to discriminate between 5mC and 5hmC will provide key insights into the role of these cytosine modifications in development and disease. Citation Format: Daniel J. Evanich, V. K. Chaithanya Ponnaluri, Vaishnavi Panchapakesa, Ariel Erijman, Matthew A. Campbell, Nan Dai, Bradley W. Langhorst, Romualdas Vaisvila, Louise Williams. Sequencing of 5-hydroxymethylcytosine to single base resolution [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2023; Part 1 (Regular and Invited Abstracts); 2023 Apr 14-19; Orlando, FL. Philadelphia (PA): AACR; Cancer Res 2023;83(7_Suppl):Abstract nr 6018.

Primer exchange reaction-amplified protein-nucleic acid interactions for ultrasensitive and specific microRNA detection

Protein-nucleic acid interactions are not only fundamental to genetic regulation and cellular metabolism, but molecular basis to artificial biosensors. However, such interactions are generally weak and dynamic, making their specific and sensitive quantitative detection challenging. By using primer exchange reaction (PER)-amplified protein-nucleic acid interactions, we here design a universal and ultrasensitive electrochemical sensor to quantify microRNAs (miRNAs) in blood. This PER-miR sensor leverages specific recognition between S9.6 antibodies and miRNA/DNA hybrids to couple with PER-derived multi-enzyme catalysis for ultrasensitive miRNA detection. Surface binding kinetic analysis shows a rational Kd (8.9 nM) between the miRNA/DNA heteroduplex and electrode-attached S9.6 antibody. Based on such a favorable affinity, the programmable PER amplification enables the sensor to detect target miRNAs with sensitivity up to 90.5 aM, three orders of magnitude higher than that without PER in routine design, and with specificity of single-base resolution. Furthermore, the PER-miR sensor allows detecting multiple miRNAs in parallel, measuring target miRNA in lysates across four types of cell lines, and differentiating tumor patients from healthy individuals by directly analyzing the human blood samples (n = 40). These advantages make the sensor a promising tool to enable quantitative sensing of biomolecular interactions and precision diagnostics.

Nucleic Acid Amplification by Template-Dominated Click Chemistry for Ultrasensitive DNA/RNA Detection on an Electrochemical Readout Platform.

As an enzyme-free exponential nucleic acid amplification method, the click chemistry-mediated ligation chain reaction (ccLCR) has shown great prospects in the molecular diagnosis. However, the current optics-based ccLCR is challenged by remarkable nonspecific amplification, severely hindering its future application. This study demonstrated that the severe nonspecific amplification was generated probably due to high random collision in the high DNA probe concentration (μM level). To solve this hurdle, a nucleic acid template-dominated ccLCR was constructed using nM-level DNA probes and read on an electrochemical platform (cc-eLCR). Under the optimal conditions, the proposed cc-eLCR detected a low-level nucleic acid target (1 fM) with a single-base resolution. Furthermore, this assay was applied to detect the target of interest in cell extracts with a satisfactory result. The proposed cc-eLCR offers huge possibility for click chemistry-mediated enzyme-free exponential nucleic acid amplification in the application of medical diagnosis and biomedical research.

DNA methylation variation is crucial to restore adventitious rooting ability during in vitro shoot culture-induced rejuvenation in apple rootstock.

In vitro shoot culture has been widely used for restoring adventitious rooting ability in rooting recalcitrant woody perennial species for the past few decades, but its molecular mechanism is largely uncovered. DNA methylation is an essential epigenetic mark that participates in many biological processes. Recent reports suggested a role of DNA methylation in vitro culture in plants. In this study, we characterized the single-base resolution DNA methylome and transcriptome of adult and in vitro shoot culture-induced rejuvenation cuttings of apple rootstock M9T337. We found a global decrease in DNA methylation during rejuvenation, which may be correlated with increased expression of DNA demethylase genes and decreased expression of DNA methyltransferase genes. We additionally documented DNA hypomethylation in 'T337'_R in gene protomer associated with higher transcript levels of several adventitious rooting-related genes. The application of a DNA methylation inhibitor (5-azacytidine) enhanced the adventitious rooting ability and the expression level of adventitious rooting-related genes, such as, MdANT, MdMPK3, MdABCB21, MdCDC48, MdKIN8B, pri-MdMIR156a5 and pri-MdMIR156a12. Together, the DNA hypomethylation is critical for the rejuvenation-dependent adventitious rooting ability in apple rootstock. In addition, increased DNA methylation was also found in thousands of genes in 'T337'_R. We additionally documented that DNA hypermethylation is required for inhibition of adventitious rooting-repressed genes, such as MdGAD5a, encoding glutamate decarboxylase, which can catalyze glutamate decarboxylated to form γ-aminobutyric acid (GABA). Our results revealed that in vitro shoot culture-dependent DNA methylation variation plays important roles in adventitious rooting in apple rootstock.

Localized assembly for long reads enables genome-wide analysis of repetitive regions at single-base resolution in human genomes

BackgroundLong-read sequencing technologies have the potential to overcome the limitations of short reads and provide a comprehensive picture of the human genome. However, the characterization of repetitive sequences by reconstructing genomic structures at high resolution solely from long reads remains difficult. Here, we developed a localized assembly method (LoMA) that constructs highly accurate consensus sequences (CSs) from long reads.MethodsWe developed LoMA by combining minimap2, MAFFT, and our algorithm, which classifies diploid haplotypes based on structural variants and CSs. Using this tool, we analyzed two human samples (NA18943 and NA19240) sequenced with the Oxford Nanopore sequencer. We defined target regions in each genome based on mapping patterns and then constructed a high-quality catalog of the human insertion solely from the long-read data.ResultsThe assessment of LoMA showed a high accuracy of CSs (error rate < 0.3%) compared with raw data (error rate > 8%) and superiority to a previous study. The genome-wide analysis of NA18943 and NA19240 identified 5516 and 6542 insertions (≥ 100 bp), respectively. Most insertions (~ 80%) were derived from tandem repeats and transposable elements. We also detected processed pseudogenes, insertions in transposable elements, and long insertions (> 10 kbp). Finally, our analysis suggested that short tandem duplications are associated with gene expression and transposons.ConclusionsOur analysis showed that LoMA constructs high-quality sequences from long reads with substantial errors. This study revealed the true structures of the insertions with high accuracy and inferred the mechanisms for the insertions, thus contributing to future human genome studies. LoMA is available at our GitHub page: https://github.com/kolikem/loma.

EBS-seq: enrichment-based method for accurate analysis of 5-hydroxymethylcytosine at single-base resolution

BackgroundA growing body of research has emphasized 5-hydroxymethylcytosine (5hmC) as an important epigenetic mark. High-resolution methods to detect 5hmC require high sequencing depth and are therefore expensive. Many studies have used enrichment-based methods to detect 5hmC; however, conventional enrichment-based methods have limited resolution. To overcome these limitations, we developed EBS-seq, a cost-efficient method for 5hmC detection with single-base resolution that combines the advantages of high-resolution methods and enrichment-based methods.ResultsEBS-seq uses selective labeling of 5hmC, deamination of cytosine and 5-methylcytosine, pull-down of labeled 5hmC, and C-to-T conversion during DNA amplification. Using this method, we profiled 5hmC in HEK293T cells and two colorectal cancer samples. Compared with conventional enrichment-based 5hmC detection, EBS-seq improved 5hmC signals by localizing them at single-base resolution. Furthermore, EBS-seq was able to determine 5hmC levels in CpG-dense regions where distortion of signals can occur, such as CpG islands and CpG shores. Comparing EBS-seq and conventional high-resolution 5hmC detection by ACE-seq, we showed that EBS-seq is more effective at finding 5hmC sites. Using EBS-seq, we found strong associations between gene expression and gene-body 5hmC content in both HEK293T cells and colorectal cancer samples.ConclusionsEBS-seq is a reliable and cost-efficient method for 5hmC detection because it simultaneously enriches 5hmC-containing DNA fragments and localizes 5hmC signals at single-base resolution. This method is a promising choice for 5hmC detection in challenging clinical samples with low 5hmC levels, such as cancer tissues.

Abstract P5-06-05: Whole-genome bisulfite sequencing of single circulating tumor cells identifies cellular methylation heterogeneity in metastatic breast cancer

Abstract Background: Although different patterns and changes in DNA methylation play an important role in cancer progression and tumor subtype differentiation, methylation remains relatively less understood in the context of tumor heterogeneity. Intra-tumor heterogeneity influences chemotherapy response, resistance development, and metastatic progression, and can be identified with liquid biopsy, a non-invasive approach to monitor cancer progression and provide predictive information. In this study, we sequenced circulating tumor cells (CTCs) to interrogate whole-genome methylation at single-cell level and single-base resolution. Methods: We enumerated CTCs in 7.5ml of whole blood samples from over 130 metastatic breast cancer patients using CellSearch, and selected blood samples with more than 20 CTCs to be further processed by DEPArray. After single-cell library construction using the Swift Accel-NGS Adaptase Module, whole-genome bisulfite sequencing (WGBS) was performed. The sequencing data were then aligned to the reference genome, and the methylation information was extracted by Bismark. We compared methylation profiles between CTC and WBC samples from each patient to identify differentially methylated regions (DMRs) using Methylkit. Multiple-comparison was conducted by SMART2 to identify DMRs within CTCs. Heatmap was plotted based on the regional methylation rates of DMRs. CpG and genic annotations were performed by annotatr. The pathway and function enrichment analyses (KEGG and GO) were conducted using clusterProfiler. Genomic region-based enrichment was conducted using LOLA. t-SNE clustering was used to investigate intra-patient heterogeneity. Results: A total of 376 cells from nine patients that passed our selection criteria were isolated for WGBS. The mean sequencing depth of all single cells was 0.8×, and an average mapping rate of 54% was achieved. We first compared CTCs and WBCs, and observed a global hypo-methylation pattern in CTCs (average methylation rate 68.8% in CTCs vs. 76.7% in WBCs). Moreover, approximately 2,000 highly differentially methylated regions were found in each patient, with hundreds of hypo- or hyper-methylated genes related to these DMRs. Hypo-methylated DMRs showed genomic region-based enrichment in breast-, prostate-, and fibroblast-related region sets. Then, we investigated methylation heterogeneity within CTCs, where t-SNE clustering identified four different subsets in one representative patient. About 1,500 hypo-methylated DMRs were found, differentiating those four subgroups. KEGG pathway analysis indicated enrichment in the Rap1 signaling pathway and focal adhesion. GO enrichment analysis highlighted the regulation of GTPase activity and membrane potential. Conclusions: Our study identified differentially methylated regions in CTCs of metastatic breast cancer patients, and demonstrated intra-patient heterogeneity based on cellular methylation information. Future studies are warranted to validate our findings and explore their biological mechanisms and clinical relevance. Citation Format: Rui Luo, Weelic Chong, Zhenchao Zhang, Maysa Abu-Khalaf, Daniel Silver, Frederick Fellin, Rebecca Jaslow, AnaMaria Lopez, Terrence Cescon, Kevan Ip, Ronald Myers, Qiang Wei, Bingshan Li, Chun Wang, Hushan Yang. Whole-genome bisulfite sequencing of single circulating tumor cells identifies cellular methylation heterogeneity in metastatic breast cancer [abstract]. In: Proceedings of the 2022 San Antonio Breast Cancer Symposium; 2022 Dec 6-10; San Antonio, TX. Philadelphia (PA): AACR; Cancer Res 2023;83(5 Suppl):Abstract nr P5-06-05.

The Transcriptome-Wide Mapping of 2-Methylthio-N6-isopentenyladenosine at Single-Base Resolution.

Hundreds of modified bases have been identified and enzymatically modified to transfer RNAs (tRNAs) to regulate RNA function in various organisms. 2-Methylthio-N6-isopentenyladenosine (ms2i6A), a hypermodified base found at tRNA position 37, exists in both prokaryotes and eukaryotes. ms2i6A is traditionally identified by separating and digesting each tRNA from total RNA using RNA mass spectrometry. A transcriptome-wide and single-base resolution method that enables absolute mapping of ms2i6A along with analysis of its distribution in different RNAs is lacking. Here, through chemoselective methylthio group bioconjugation, we introduce a new approach (redox activated chemical tagging sequencing, ReACT-seq) to detect ms2i6A transcriptome-wide at single-base resolution. Using the chemoselectivity between the methylthio group and oxaziridine group, ms2i6A is bio-orthogonally tagged with an azide group without interference of canonical nucleotides, advancing enrichment of methylthio group modified RNAs prior to sequencing. ReACT-seq was demonstrated on nine known tRNAs and proved to be highly accurate, and the reverse transcription stop (RT-stop) character enables ReACT-seq detection at single-base resolution. In addition, ReACT-seq identified that the modification of ms2i6A is conservative and may not exist in other RNAs.

Universal and highly accurate detection of circulating tumor DNA mutation in non-small cell lung cancer based on CRISPR/Cas12a system

Circulating tumor DNA (ctDNA) is a promising biomarker for real-time, minimally invasive diagnostics and monitoring in patients with non-small cell lung cancer (NSCLC), especially when representative tissue biopsies are not available. However, the current methods for ctDNA detection are cumbersome and expensive. While advanced CRISPR/Cas-based assays offer advantages of simplicity, low cost and high sensitivity, their application for ctDNA detection is restricted by the requirement of a protospacer adjacent motif (PAM) near the mutation site and off-target cleavages (i.e., false-positive results) due to the extreme similarities between the mutant and wild-type sequences, especially single nucleotide variants. Herein, we propose a novel strategy comprising recombinase polymerase amplification (RPA) and CRISPR/Cas12a to detect ctDNA with high universality and accuracy. The use of artificially inserted PAMs by modified RPA primers or suboptimal PAMs unlocks the PAM restriction; introducing single- or double-base mismatches in CRISPR RNA effectively reduces the off-target effects and improves the specificity to single-base resolution. Under optimized conditions, this method detected ctDNA mutations with a limit of detection at 100 aM and identified mutations down to 0.02% variant allele frequency in 50 min, requiring only isothermal control. We successfully applied this method to multiple clinical samples of NSCLC and the results were validated using real-time polymerase chain reaction analysis. In summary, we established a rapid, sensitive, universal and highly accurate method for ctDNA detection that has great potential application in the early diagnosis, therapy guidance and prognosis prediction of NSCLC.