Detection and Characterization of Methylated Circulating Tumor DNAs in Gastric Cancer

Gastric cancer is the 5th most common disease in the world and the 4th most common cause of death. It is diagnosed through esophagogastroduodenoscopy with biopsy, however there are limitations in finding lesions in the early stages. Recently, research has been actively conducted to use liquid biopsy to diagnose various cancers, including gastric cancer. Various substances derived from cancer are reflected in the blood. By analyzing these substances, it was expected that not only the presence or absence of cancer but also the type of cancer can be diagnosed. However, the amount of these substances is extremely small, and even these have various variables depending on the characteristics of the individual or the characteristics of the cancer. To overcome these, we collected methylated DNA fragments using MeDIP and compared them with normal plasma to characterize gastric cancer tissue or patient’s plasma. And we attempted to diagnose gastric cancer using the characteristics of cancer reflected in the blood through the cancer tissue and patient’s plasma. As a result, we confirmed that the consistency of common methylated fragments between tissue and plasma was approximately 41.2%, and found the possibility of diagnosing and characterizing cancer using the characteristics of the fragments through SFR and 5'end-motif analysis.

Keywords:

Subject: Biology and Life Sciences - Biochemistry and Molecular Biology

1. Introduction

1.1. Gastric Cancer and Liquid Biopsy

Despite advances in diagnosis and treatment, cancer is the second leading cause of death worldwide [1]. In particular, gastric cancer is the fifth most common and fourth most death, causing 1.1 million new cases and 7.6 hundred thousand deaths each year in worldwide [2,3]. Recently, the incidence of gastric cancer has been increasing in Korea due to changes in eating habits and various environments [4,5]. The gold standard of diagnosis for gastric cancer is esophagogastroduodenoscopy with biopsy. However, the diagnosis can be often delayed in early stage cancer. In particular, liquid biopsy research using blood is being conducted recently for early diagnosis of various diseases, including cancer [6,7]. In various diseases, including cancer, various characteristics derived from the lesion are reflected in the blood. Direct indicators include exosomes, metabolic materials, cfDNAs, and CTCs. Although it exists in extremely small quantities, analyzing its characteristics has made it possible to diagnose the disease with minimally invasive methods [8,9,10]. Therefore, for a more accurate diagnosis, the unique characteristics of each cancer type must be identified. Representative successful cases include EGFR, KRAS, and BRAF mutations in lung cancer [11,12,13,14]. Recently, omics research on various characteristics has reached a level where it is possible to diagnose not only cancer but also types of cancer [15,16,17]. However, the uniqueness of many types of cancer is still unclear.

1.2. Methylation

Changes in gene sequence cause functional changes, and in fact, many diseases, including cancer, are caused by small changes in genes [18]. It is known that it is sometimes caused by abnormal changes in expression levels rather than functional changes in genes. These epigenetic studies are widely applied in the diagnosis and prediction of cancer. Human cancer is closely associated with changes in methylation [19]. In particular, tissues sometimes have specific methylation profiles [20,21]. Based on these characteristics, it is possible to diagnose cancer and its type using the methylation pattern of cancer-derived substances present in the blood.

1.3. Fragmentomics (Size, 5’end-Motif)

cfDNA in the blood comes from various tissues and cells [22,23]. Recently, it has been revealed that fragments derived from various cancers are shorter than those derived from normal cells due to different biological mechanisms, and the possibility of diagnosis using this is being discussed [24,25,26,27].

All living organisms contain nucleases that can interact with nucleic acids and hydrolyze phosphodiester bonds. Nucleases may also exhibit a cleavage preference for single- or double-stranded nucleic acids. Some nucleases are structure or sequence-specific. In particular, there are nucleases that are predominant and specific to the cancer type or tissue [28,29,30,31]. In gastric cancer, FEN1, APE1, XPF/XPG, MRN complex, and DNase1 are typically known [30].

Many liquid biopsy studies are focused on diagnosing cancer by analyzing the characteristics of cancer-derived cell free nucleic acids in the blood. However, it is still true that in the early stages of cancer, the amount present is very small and a lot of effort and time are required for accurate analysis. In this study, in order to clarify the characteristics of cancer-derived fragments, hypermethylated fragments were collected and their length and 5'end-motif sequences were analyzed. We also confirmed whether this approach was consistent with results from tissue.

2. Results

2.1. Differentially Methylated Region (DMR)

First, changes in cfDNA methylation in the plasma of cancer patients were compared with the healthy human plasma. Paired Blood and tissue from 22 stomach cancer patients and 40 plasmas from healthy people were used (Table 1). The MEDIPS library was used to determine differentially methylated regions (DMR) [32]. The BH (Benjamini-Hochberg) method was set to be used for p-value adjustment, and areas with a value of 0.01 or less based on the adjusted p-value were selected as DMR. The methylation marker information for each cancer type was used as published in Vrba's paper [33], and the methylation marker information for each tissue was used as published in Moss' paper [34]. Total of 1,541,364 DMRs were obtained by comparing tissue and plasma of cancer patients with healthy human plasma, and the total length was approximately 813 Mb (Figure 1A and Supplementary Table S1). About 146 Mb DMRs were identified between plasma from cancer patients and healthy human plasma, with a ratio of 124 Mb gain (hypermethylation) to 22 Mb loss (hypomethylation) of approximately 6:1. Additionally, DMR coverage of 667Mb was confirmed between cancer patient tissue and healthy human plasma, and hypermethylation and hypomethylation of 528Mb and 138Mb were found, respectively (Figure 1). Between patient and healthy plasma, approximately 3% (9,711/316,507) of DMRs were found in the promoter region. Between cancer tissue and healthy plasma, 2.6% (31,922/1,224,857) of DMRs are located in promoter regions (Supplementary Table S1 and Supplementary Figure S2).

2.2. Commonly Methylated Region (CMR)

Next, we analyzed how much methylated DNA there was in common between tissue and plasma from cancer patients. The MACS-2.2.6 program was used to discover methylation regions [35]. The discovered methylation region was then used to determine the common methylation region (CMR). Methylated regions were processed in 300 bp sections. If the methylated region discovered by the MACS program occupied more than 50% of the 300 bp region, the region was considered as methylated. The region is determined by connecting adjacent windows (Supplementary Figure S3A). Overall, the average of peaks was 439,109 and the average number of CMRs was 213,548 and almost of CMRs were not located in DMR (Supplementary Figure S3B,C). When comparing the tissue and plasma of patients with that of healthy human plasma, the average size of the CMR between tissue and cfDNA was approximately 73 Mb (Figure 2). Interestingly, the average similarity was approximately 41.2% between tissue and cfDNA (Supplementary Figure S3D). This shows that cancer tissue-based methylation changes are also reflected at high levels in ctDNA from the blood.

2.3. Short Fragment Ratio (SFR)

To detect cancer, we used cfDNA obtained from MeDIP to compare the length of cfDNA in the plasma of cancer patients and the plasma of healthy donors. Many previous studies have already shown that shorter length cfDNA exists in patients with various diseases, including cancer. Although ctDNA derived from cancer is only a small portion of the total, we distinguished the difference more clearly by using the proportion of short cfDNA. The insert size estimated from mapping information was used as the length of cfMeDNA (circulating free methylated DNA). The short fragment rate (SFR) was calculated as the proportion of short fragments among all fragments. cfMeDNAs with a length of 100 to 240 bp were set as full fragments, and those with a length of 100 to 160 bp were set as short fragments. To calculate the skewness of the cfMeDNA length distribution, 100 to 240 bp fragments were used. Simply looking at the trend, the proportion of short cfDNA in the patient's plasma was high (Figure 3A,B). Next, for detailed quantification, we calculated the SFR by dividing healthies and patients and compared the differences. In conclusion, a high level of significant SFR difference was observed in cancer patients (P value=1.63E-8) (Figure 3C). Surprisingly, focusing on the CpG island resulted in a clearer distinction (P value=1.97E-11) (Figure 3D). To verify the ability to distinguish between normal and cancer using SFR, the ROC curve was calculated. At the whole genome level, the AUC was 0.909. Interestingly, applying this to the CpG island region raised the AUC to 0.950 (Figure 3E,F). We identified how SFRs differed in the promoter regions of 51 well-known cancer-related genes [36]. In particular, the SFR of the BAP1 gene promoter was statistically significantly higher in cancer patients (Figure 4).

2.4. 5’end-Motif

Short fragments occur depending on the type of DNase acting. In particular, it is known to show characteristic differences when DNA is fragmented within cancer cells. We classified the collected short cfDNAs using the sequences of their 5’end-motifs. The length of the motif used in the 5’end-motif frequency analysis is 4 bp. In other previous studies, abnormal activation of DNase I was reported in gastric cancer [30]. As a result of the action of DNase I, thiamine end motif is generated, and in our results, T-end SFRs accounted for a statistically significantly high proportion in top 10 rank (Table 2). The t-test was used in the statistical test for frequency differences between groups, and p-value adjustment was performed by the Benjamini-Hochberg procedure.

3. Discussion

We hypothesized that there would be a clearer difference between cancer and normal if methylation changes were limited to cfDNA. In this study, hypermethylated cfDNA predominates in both cancer and healthy donor plasma due to MeDIP capture. Hypermethylation is more prevalent in cancer genomes. As our results show, MeDIP had the effect of making these characteristics more evident. It is well known that many diseases, including cancer, are related to methylation change. Moreover, methylation change related with age, gender, species, and organ specific factor [37]. However, to be free from these factors, individual cancer-derived substances were compared with a pool of 40 healthy donor plasma results.

Of course, comparing within the same type of specimen can yield more accurate results than comparing plasma and tissue. Nevertheless, there are many studies attempting to determine whether the characteristics of cancer tissue are reflected in blood for liquid biopsy. We obtained common DMR between cancer tissue and plasma of cancer patients compare with plasma of healthy people. Although a more in-depth study of the obtained results was not conducted, we believe it will be good material for future research.

Because methylation differences were well reflected using MeDIP, we expected to discover CMRs within DMRs. Extreme DMR was also set between cancer and normal plasma for dramatic results. Unfortunately, CMR seems to have little to do with DMR. This result is expected to have occurred because a very small amount of cfDNA including ctDNA was used as a control. Nonetheless, surprisingly, there was a match rate of about 40% between cancer tissue and cancer patient plasma. This is the result of the characteristics of cancer tissue being reflected in the blood. Sometimes, many tissues have tissue-specific mutations and methylation patterns depending on the type of cancer [38,39]. It is expected that the CMR may contain clues that can predict not only the presence or absence of cancer but also the type of cancer.

Short fragment DNA provides useful clues for ctDNA classification using liquid biopsy. Especially in the early stages, ctDNA accounts for less than 1% of cfDNA [40,41]. Therefore, there are limitations in clearly distinguishing cancer from normal with only these small changes. We hypothesized that cfDNA derived from cancer could be more clearly identified through methylation capture. Surprisingly, when our method used SFR to distinguish cancer from normal, the AUC of ROC was 0.909. Additionally, when limited to CpG islands, the AUC increased to 0.95. This shows that SFR can more clearly distinguish between cancer and normal conditions under certain conditions such as methylation. We applied methylation capture and CpG island piecemeal, but we expect that better discrimination ability will be achieved in the future.

All living organisms have enzymes that break down nucleic acids. And DNases and RNases act differently depending on the type of tissue or disease. The action of these different DNases results in shorter length of cancer-derived cfDNA. In particular, all nucleic acids exposed due to apoptosis or necrosis are fragmented by various DNases and RNases, and some are exuded into the blood. Verification is necessary to clearly reveal the direct basis for the relationship between DNase and 5’ end motif in the results presented in our study. Nevertheless, our study shows that cancer diagnosis is possible by utilizing the characteristics of ctDNA in cfDNA. In particular, if the DNase action is more clearly identified for each cancer type, it is expected that cancer can be diagnosed using only the short sequence of the fragment (only 4 bp).

Although the genes we used in our analysis were limited to 22, these genes are already well-known cancer-related genes. In particular, BAP1 is a well-known TSG, and its role in gastric cancer has also been identified. This study did not conduct experiments to support the function of BAP1 in gastric cancer. However, shorter fragmentation of TSGs and oncogenes may not be unrelated to the development of cancer. Characteristics of the fragments, such as differences in methylation or length of these genes in the blood, may be the result of their origin from actual cancer tissue. In other words, identification of these characteristics can provide a clear basis for cancer diagnosis using liquid biopsy. Unfortunately, gastric was not ranked in classification using tissue-specific hypermethylation markers (data not shown). Perhaps methylation capture may simplify features of tissue, including stomach, that we were unaware of.

In conclusion, we initiated this study out of curiosity to see whether known properties of ctDNA could be further elucidated through methylation. Although there were limitations in obtaining normal controls, we believe that our approach offers new possibilities for future gastric cancer diagnosis and characterization using liquid biopsy. We plan to conduct large-scale clinical research using the methods identified in this study. In addition, by securing a strong normal control group for making clarify characteristics of tissue and plasma, and applying various analysis methods to cfDNA, including ctDNA, it is expected that we will be able to get one step closer to cancer diagnosis using liquid biopsy.

4. Materials and Methods

4.1. 40 Healthy Donor Plasmas, 22 Pairs of GC Tissue and Blood, DNA Extraction

In this study, tissue and blood were collected in pairs from 22 patients (IRB : CUH2022-8-012). And for comparative study and normalization, 40 healthy human plasmas were used. cfDNA was collected using by Apostile cfDNA prep kit (Beckman Coulter). Genomic DNA was extracted from the tissue using a DNeasy kit (Qiagen), and cut into sizes of less than 200bp using Bioruptor Pico (Diagenode). Size and quality of all nucleic acid was measured using TapeStation 4150 (Agilent Technologies) and Qubit 4 Fluorometer (Thermofisher Scientific).

4.2. cfMeDIP Seq

For the MeDIP seq of cfDNA and fragmented gDNA, we followed previous other methods [42]. The sequencing data was generated using MGI G400 sequencer. Each output of sequencing produced secured more than 2 times depth.

4.3. Data Processing and Mapping

The Cutadapt-2.4 program was used to select high-quality sequence data [43]. Adapter sequences and low-quality bases from the back of the original sequence were removed to obtain a high-quality sequence. The high-quality sequences were mapped to the hg38 human genome and the mem command of the BWA-0.7.15 program was used for the mapping task [44]. Duplicate sequence data were removed using the MarkDuplicates command in the Picard-2.7.1 program. Flowcharts and tools used in data processing such as DMR, CMR, and SFR are available in Supplementary Figure S1.

Supplementary Materials

The following supporting information can be downloaded at the website of this paper posted on Preprints.org.

Author Contributions

Seung Young Seo: Conceptualization (Equal); Data curation (Equal); Formal analysis (Equal); Resources (Equal); Supervision (Equal); Writing-review & editing (Equal). Sang Hee Youn: Conceptualization (Equal); Data curation (Equal); Formal analysis (Equal); Resources (Equal); Writing-review & editing (Supporting). Jin-Han Bae: Conceptualization (Equal); Data curation (Lead); Formal analysis (Lead); Methodology (Lead); Software (Lead); Validation (Lead); Visualization (Lead); Writing-original draft (Lead). Sung-Hun Lee: Conceptualization (Equal); Data curation (Equal); Formal analysis (Equal); Supervision (Equal); Writing-review & editing (Equal). Sun Young Lee: Conceptualization (Equal); Funding acquisition (Lead); Investigation (Lead) Project administration (Lead); Resources (Lead); Supervision (Lead); Writing-review & editing (Equal).

Funding

This paper was supported by Fund of Biomedical Research Institute, Jeonbuk National University Hospital (CUH2022-0006).

Institutional Review Board Statement

The study was approved by the Ethics Committee/Institutional Review Board (IRB) of Jeonbuk National University Hospital (IRB : CUH2022-8-012). In addition, prior to the commencement of the study, written informed consent was obtained from all participants.

Data Availability Statement

All data are available upon reasonable request to the corresponding author of the article.

Conflicts of Interest statement

The authors have no conflicts of interest related to this work.

References

Nagai, H.; Kim, Y.H. Cancer prevention from the perspective of global cancer burden patterns. Journal of thoracic disease 2017, 9, 448. [Google Scholar] [CrossRef]
Sung, H.; Ferlay, J.; Siegel, R.L.; Laversanne, M.; Soerjomataram, I.; Jemal, A.; Bray, F. Global cancer statistics 2020: GLOBOCAN estimates of incidence and mortality worldwide for 36 cancers in 185 countries. CA: a cancer journal for clinicians 2021, 71, 209–249. [Google Scholar] [CrossRef] [PubMed]
Arnold, M.; Rutherford, M.J.; Bardot, A.; Ferlay, J.; Andersson, T.M.; Myklebust, T.Å.; Tervonen, H.; Thursfield, V.; Ransom, D.; Shack, L. Progress in cancer survival, mortality, and incidence in seven high-income countries 1995–2014 (ICBP SURVMARK-2): a population-based study. The Lancet Oncology 2019, 20, 1493–1505. [Google Scholar] [CrossRef] [PubMed]
Park, S.H.; Kang, M.J.; Yun, E.H.; Jung, K.-W. Epidemiology of gastric cancer in Korea: trends in incidence and survival based on Korea Central Cancer Registry Data (1999–2019). Journal of Gastric Cancer 2022, 22, 160. [Google Scholar] [CrossRef] [PubMed]
Hong, S.; Won, Y.-J.; Lee, J.J.; Jung, K.-W.; Kong, H.-J.; Im, J.-S.; Seo, H.G. Cancer statistics in Korea: incidence, mortality, survival, and prevalence in 2018. Cancer Research and Treatment: Official Journal of Korean Cancer Association 2021, 53, 301. [Google Scholar] [CrossRef] [PubMed]
Ma, S.; Zhou, M.; Xu, Y.; Gu, X.; Zou, M.; Abudushalamu, G.; Yao, Y.; Fan, X.; Wu, G. Clinical application and detection techniques of liquid biopsy in gastric cancer. Molecular Cancer 2023, 22, 7. [Google Scholar] [CrossRef] [PubMed]
Lone, S.N.; Nisar, S.; Masoodi, T.; Singh, M.; Rizwan, A.; Hashem, S.; El-Rifai, W.; Bedognetti, D.; Batra, S.K.; Haris, M. Liquid biopsy: a step closer to transform diagnosis, prognosis and future of cancer treatments. Molecular cancer 2022, 21, 79. [Google Scholar] [CrossRef] [PubMed]
Soda, N.; Clack, K.; Shiddiky, M.J. Recent advances in liquid biopsy technologies for cancer biomarker detection. Sensors & Diagnostics 2022, 1, 343–375. [Google Scholar] [CrossRef]
Alix-Panabières, C.; Pantel, K. Clinical applications of circulating tumor cells and circulating tumor DNA as liquid biopsy. Cancer discovery 2016, 6, 479–491. [Google Scholar] [CrossRef]
Diehl, F.; Schmidt, K.; Choti, M.A.; Romans, K.; Goodman, S.; Li, M.; Thornton, K.; Agrawal, N.; Sokoll, L.; Szabo, S.A. Circulating mutant DNA to assess tumor dynamics. Nature medicine 2008, 14, 985–990. [Google Scholar] [CrossRef]
Fujii, H.; Nagakura, H.; Kobayashi, N.; Kubo, S.; Tanaka, K.; Watanabe, K.; Horita, N.; Hara, Y.; Nishikawa, M.; Miura, K. Liquid biopsy for detecting epidermal growth factor receptor mutation among patients with non-small cell lung cancer treated with afatinib: a multicenter prospective study. BMC cancer 2022, 22, 1035. [Google Scholar] [CrossRef] [PubMed]
Nigro, M.C.; Marchese, P.V.; Deiana, C.; Casadio, C.; Galvani, L.; Di Federico, A.; De Giglio, A. Clinical utility and application of liquid biopsy genotyping in lung cancer: a comprehensive review. Lung Cancer: Targets and Therapy 2023, 11–25. [Google Scholar] [CrossRef] [PubMed]
Nacchio, M.; Sgariglia, R.; Gristina, V.; Pisapia, P.; Pepe, F.; De Luca, C.; Migliatico, I.; Clery, E.; Greco, L.; Vigliar, E. KRAS mutations testing in non-small cell lung cancer: The role of Liquid biopsy in the basal setting. Journal of thoracic disease 2020, 12, 3836. [Google Scholar] [CrossRef] [PubMed]
Fernandes, G.; Rodrigues, A.; Matos, C.; Barata, F.; Cirnes, L.; Ferreira, L.; Lopes, J.A.; Felizardo, M.; Fidalgo, P.; Brito, U. Liquid biopsy in the management of advanced lung cancer: Implementation and practical aspects. Cancer Treatment and Research Communications 2023, 36, 100725. [Google Scholar] [CrossRef] [PubMed]
Di Sario, G.; Rossella, V.; Famulari, E.S.; Maurizio, A.; Lazarevic, D.; Giannese, F.; Felici, C. Enhancing clinical potential of liquid biopsy through a multi-omic approach: A systematic review. Frontiers in Genetics 2023, 14, 1152470. [Google Scholar] [CrossRef] [PubMed]
Amelio, I.; Bertolo, R.; Bove, P.; Buonomo, O.C.; Candi, E.; Chiocchi, M.; Cipriani, C.; Di Daniele, N.; Ganini, C.; Juhl, H. Liquid biopsies and cancer omics. Cell Death Discovery 2020, 6, 131. [Google Scholar] [CrossRef] [PubMed]
Chen, G.; Zhang, J.; Fu, Q.; Taly, V.; Tan, F. Integrative analysis of multi-omics data for liquid biopsy. British journal of cancer 2023, 128, 505–518. [Google Scholar] [CrossRef] [PubMed]
Dulbecco, R. A turning point in cancer research: sequencing the human genome. Science 1986, 231, 1055–1056. [Google Scholar] [CrossRef]
Jones, P.A.; Baylin, S.B. The epigenomics of cancer. Cell 2007, 128, 683–692. [Google Scholar] [CrossRef]
Zhou, J.; Sears, R.L.; Xing, X.; Zhang, B.; Li, D.; Rockweiler, N.B.; Jang, H.S.; Choudhary, M.N.; Lee, H.J.; Lowdon, R.F. Tissue-specific DNA methylation is conserved across human, mouse, and rat, and driven by primary sequence conservation. BMC genomics 2017, 18, 1–17. [Google Scholar] [CrossRef]
Ghosh, S.; Yates, A.J.; Frühwald, M.C.; Miecznikowski, J.C.; Plass, C.; Smiraglia, D. Tissue specific DNA methylation of CpG islands in normal human adult somatic tissues distinguishes neural from non-neural tissues. Epigenetics 2010, 5, 527–538. [Google Scholar] [CrossRef]
Sun, K.; Jiang, P.; Chan, K.A.; Wong, J.; Cheng, Y.K.; Liang, R.H.; Chan, W.-k.; Ma, E.S.; Chan, S.L.; Cheng, S.H. Plasma DNA tissue mapping by genome-wide methylation sequencing for noninvasive prenatal, cancer, and transplantation assessments. Proceedings of the National Academy of Sciences 2015, 112, E5503–E5512. [Google Scholar] [CrossRef] [PubMed]
Li, W.; Li, Q.; Kang, S.; Same, M.; Zhou, Y.; Sun, C.; Liu, C.-C.; Matsuoka, L.; Sher, L.; Wong, W.H. CancerDetector: ultrasensitive and non-invasive cancer detection at the resolution of individual reads using cell-free DNA methylation sequencing data. Nucleic acids research 2018, 46, e89–e89. [Google Scholar] [CrossRef] [PubMed]
Mouliere, F.; Rosenfeld, N. Circulating tumor-derived DNA is shorter than somatic DNA in plasma. Proceedings of the National Academy of Sciences 2015, 112, 3178–3179. [Google Scholar] [CrossRef] [PubMed]
Chen, E.; Cario, C.L.; Leong, L.; Lopez, K.; Márquez, C.P.; Chu, C.; Li, P.S.; Oropeza, E.; Tenggara, I.; Cowan, J. Cell-free DNA concentration and fragment size as a biomarker for prostate cancer. Scientific reports 2021, 11, 5040. [Google Scholar] [CrossRef] [PubMed]
Lo, Y.D.; Han, D.S.; Jiang, P.; Chiu, R.W. Epigenetics, fragmentomics, and topology of cell-free DNA in liquid biopsies. Science 2021, 372, eaaw3616. [Google Scholar] [CrossRef] [PubMed]
Cristiano, S.; Leal, A.; Phallen, J.; Fiksel, J.; Adleff, V.; Bruhm, D.C.; Jensen, S.Ø.; Medina, J.E.; Hruban, C.; White, J.R. Genome-wide cell-free DNA fragmentation in patients with cancer. Nature 2019, 570, 385–389. [Google Scholar] [CrossRef] [PubMed]
Desai, N.A.; Shankar, V. Single-strand-specific nucleases. FEMS microbiology reviews 2003, 26, 457–491. [Google Scholar] [CrossRef] [PubMed]
Balian, A.; Hernandez, F.J. Nucleases as molecular targets for cancer diagnosis. Biomarker Research 2021, 9, 1–16. [Google Scholar] [CrossRef]
Yang, W. Nucleases: diversity of structure, function and mechanism. Quarterly reviews of biophysics 2011, 44, 1–93. [Google Scholar] [CrossRef]
Zhou, Z.; Ma, M.-J.L.; Chan, R.W.; Lam, W.J.; Peng, W.; Gai, W.; Hu, X.; Ding, S.C.; Ji, L.; Zhou, Q. Fragmentation landscape of cell-free DNA revealed by deconvolutional analysis of end motifs. Proceedings of the National Academy of Sciences 2023, 120, e2220982120. [Google Scholar] [CrossRef] [PubMed]
Lienhard, M.; Grimm, C.; Morkel, M.; Herwig, R.; Chavez, L. MEDIPS: genome-wide differential coverage analysis of sequencing data derived from DNA enrichment experiments. Bioinformatics 2014, 30, 284–286. [Google Scholar] [CrossRef]
Vrba, L.; Futscher, B.W. A suite of DNA methylation markers that can detect most common human cancers. Epigenetics 2018, 13, 61–72. [Google Scholar] [CrossRef] [PubMed]
Moss, J.; Magenheim, J.; Neiman, D.; Zemmour, H.; Loyfer, N.; Korach, A.; Samet, Y.; Maoz, M.; Druid, H.; Arner, P. Comprehensive human cell-type methylation atlas reveals origins of circulating cell-free DNA in health and disease. Nature communications 2018, 9, 5068. [Google Scholar] [CrossRef] [PubMed]
Zhang, Y.; Liu, T.; Meyer, C.A.; Eeckhoute, J.; Johnson, D.S.; Bernstein, B.E.; Nusbaum, C.; Myers, R.M.; Brown, M.; Li, W. Model-based analysis of ChIP-Seq (MACS). Genome biology 2008, 9, 1–9. [Google Scholar] [CrossRef] [PubMed]
Lee, S.H.; Park, C.K.; Lee, S.Y.; Choi, C.M. Treatment outcomes and safety of afatinib in advanced squamous cell lung cancer progressed after platinum-based doublet chemotherapy and immunotherapy (SPACE study). Thoracic Cancer 2021, 12, 1264–1268. [Google Scholar] [CrossRef]
Boks, M.P.; Derks, E.M.; Weisenberger, D.J.; Strengman, E.; Janson, E.; Sommer, I.E.; Kahn, R.S.; Ophoff, R.A. The relationship of DNA methylation with age, gender and genotype in twins and healthy controls. PloS one 2009, 4, e6767. [Google Scholar] [CrossRef]
Lakshminarasimhan, R.; Liang, G. The role of DNA methylation in cancer. DNA Methyltransferases-Role and Function 2016, 151–172. [Google Scholar] [CrossRef] [PubMed]
Nishiyama, A.; Nakanishi, M. Navigating the DNA methylation landscape of cancer. Trends in Genetics 2021, 37, 1012–1027. [Google Scholar] [CrossRef]
Stejskal, P.; Goodarzi, H.; Srovnal, J.; Hajdúch, M.; van’t Veer, L.J.; Magbanua, M.J.M. Circulating tumor nucleic acids: biology, release mechanisms, and clinical relevance. Molecular cancer 2023, 22, 15. [Google Scholar] [CrossRef]
Neumann, M.H.; Bender, S.; Krahn, T.; Schlange, T. ctDNA and CTCs in liquid biopsy–current status and where we need to progress. Computational and structural biotechnology journal 2018, 16, 190–195. [Google Scholar] [CrossRef] [PubMed]
Shen, S.Y.; Burgener, J.M.; Bratman, S.V.; De Carvalho, D.D. Preparation of cfMeDIP-seq libraries for methylome profiling of plasma cell-free DNA. Nature protocols 2019, 14, 2749–2780. [Google Scholar] [CrossRef] [PubMed]
Martin, M. Cutadapt removes adapter sequences from high-throughput sequencing reads. EMBnet. journal 2011, 17, 10–12. [Google Scholar] [CrossRef]
Li, H.; Durbin, R. Fast and accurate short read alignment with Burrows–Wheeler transform. bioinformatics 2009, 25, 1754–1760. [Google Scholar] [CrossRef]

Figure 1. Size and composition of DMR obtained by comparing tissue and plasma of gastric cancer patients with healthy donor plasma. (A) The graph shows the total number and size of DMRs obtained by comparing the gastric cancer patient's plasma or tissue to healthy donor plasma. (B) Red and blue bars represent the proportion of methylation gain and loss among total DMRs, respectively, when compared to heathy donor plasma. (C) and (D) The graph shows the total size for methylation gained and methylation lost.

Figure 2. Size and distribution of CMR obtained by comparing tissue and plasma of gastric cancer patients. Red and black bars represent DMRs found only in the patient's plasma and tissue, respectively. The green bar shows commonly presented methylation region in the patient's plasma and tissue.

Figure 3. Distribution of SFR and ability to distinguish between gastric cancer and normal. (A) The graph shows the length and distribution of cfDNA obtained from the plasma of cancer patients (red) and healthy people (black). (B) SFRs were displayed for each individual and divided into genomic regions (grey) and CpG islands (red and black). (C) and (D) The graph shows the SFR between the two groups of cancer patients and healthy people in genomic and CpG island region, respectively. (E) and (F) The graph shows ROC curve and AUC score of SFR in genomic and CpG island region, respectively. ROC, receiver operating characteristic; AUC, area under the ROC curve.

Figure 4. Differences in SFR within cancer-related genes. The graph shows the difference in SFR in the plasma of 22 patients (red) and 40 healthy people (blue) for each promoter of 51 cancer-related genes.

Table 1. Information of specimens in study.

Group	Ethnicity	Tissue	Sample	Capture	# Samples	# Individual
Gastric Cancer	Korean	Plasma	cfDNA	MeDIP	22	22
Gastric Cancer	Korean	Cancer	gDNA	MeDIP	22	22
Healthy	Hispanic	Plasma	cfDNA	MeDIP	40	40

Table 2. List of top 10 significant 5’ End-Motif.

Motif	Average Frequency (Rate)		P-value	Adjusted P-value
Motif	Cancer Patients	Healthy Individuals	P-value	Adjusted P-value
TCAA	0.002575	0.001606	1.35E-15	1.85E-13
TGCG	0.001420	0.000867	1.44E-15	1.85E-13
TACG	0.000666	0.000449	2.25E-15	1.92E-13
AGCA	0.003536	0.004193	4.81E-15	3.08E-13
TCCT	0.004531	0.003161	2.64E-14	1.35E-12
TCAG	0.002890	0.001785	3.16E-14	1.35E-12
AGCC	0.003252	0.003959	1.42E-13	5.18E-12
TCAC	0.003681	0.002555	1.81E-13	5.79E-12
TCCG	0.000936	0.000623	2.76E-13	7.85E-12
ATCA	0.001637	0.002306	7.14E-13	1.83E-11

Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

© 2024 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (http://creativecommons.org/licenses/by/4.0/).

Copyright: This open access article is published under a Creative Commons CC BY 4.0 license, which permit the free download, distribution, and reuse, provided that the author and preprint are cited in any reuse.

MDPI Initiatives

Important Links

Choose an area of interest and we will send you notifications of new preprints at your preferred frequency.

Disclaimer