ORIGINAL RESEARCH article

Front. Oncol., 10 February 2023

Sec. Hematologic Malignancies

Volume 13 - 2023 | https://doi.org/10.3389/fonc.2023.1109715

In-depth circulating tumor DNA sequencing for prognostication and monitoring in natural killer/T-cell lymphomas

  • 1. Department of Laboratory Medicine, Yonsei University College of Medicine, Seoul, Republic of Korea

  • 2. Department of Laboratory Medicine and Genetics, Samsung Medical Center, Sungkyunkwan University School of Medicine, Seoul, Republic of Korea

  • 3. Dxome Co. Ltd, 8, Seongnam-si, Gyeonggi-do, Republic of Korea

  • 4. Division of Haematology and Oncology, Department of Medicine, Samsung Medical Center, Sungkyunkwan University School of Medicine, Seoul, Republic of Korea

  • 5. Department of Health Sciences and Technology, Samsung Advanced Institute for Health Sciences and Technology, Sungkyunkwan University, Seoul, Republic of Korea

  • 6. Department of Laboratory Medicine, Korea Worker’s​ Compensation & Welfare Service, Taebaek Hospital, Taebaek-si, Gangwon-do, Republic of Korea

Abstract

Background:

Epstein–Barr virus (EBV) quantitation and current imaging modalities are used for diagnosis and disease monitoring in Extranodal NK/T cell lymphoma (ENKTL) but have limitations. Thus, we explored the utility of circulating tumor DNA (ctDNA) as a diagnostic biomarker.

Methods:

Through in-depth sequencing of 118 blood samples collected longitudinally at different time points from 45 patients, we examined the mutational profile of each sample, estimated its impact on the clinical outcome, and assessed its role as a biomarker in comparison with EBV DNA quantitation.

Results:

The ctDNA concentration was correlated with treatment response, stage, and EBV DNA quantitation. The detection rate of ctDNA mutation was 54.5%, with BCOR (21%) being the most commonly mutated gene in newly diagnosed patients; TP53 mutation (33%) was the most prevalent in patients that experienced a relapse. Additionally, patients in complete remission exhibited a rapid clearance of ENKTL-related somatic mutations, while relapsed patients frequently presented with persisting or emerging mutations. We detected ctDNA mutations in EBV-negative patients (50%) and mutation clearance in EBV-positive patients in remission, suggesting ctDNA genotyping as an efficient complementary monitoring method for ENKTL. Additionally, mutated DDX3X (PFS HR, 8.26) in initial samples predicted poor outcome.

Conclusion:

Our results suggest that ctDNA analysis can be used to genotype at diagnosis and estimate the tumor burden in patients with ENKTL. Furthermore, ctDNA dynamics indicate the potential use of testing it to monitor therapeutic responses and develop new biomarkers for precision ENKTL therapy.

Background

Extranodal natural killer/T-cell lymphoma (ENKTL) is an aggressive subtype of non-Hodgkin lymphoma and the Epstein–Barr virus (EBV), which is believed to play an important role in lymphomagenesis (1), is usually detected in malignant cells. Despite an improvement in treatment strategies including non-anthracycline-based chemotherapy, the disease recurs or progresses in more than 50% of patients with ENKTL (2, 3).

Considering the high probability of relapse even after complete resolution of tumors, the detection of minimal residual disease using imaging studies is fundamental for predicting relapse and developing risk-adapted treatment strategies. However, current imaging diagnostic modalities have low sensitivity and specificity, making it difficult to detect early recurrence or progression and identify dynamic tumor responses and time-wise evolution of clones (4). Blood EBV-DNA copy numbers are widely used as an indirect biomarker in ENKTL because they can be correlated with tumor burden. However, not all patients are EBV-positive, and EBV DNA can be released from necrotic tumor cells rather than living cells (5, 6). Therefore, a new diagnostic method is needed to overcome the limitations of current monitoring modalities.

Several studies have reported various somatic gene mutations associated with ENKTL (7–9); however, a molecular classification based on genetic alterations has not yet been established for these tumor cells. A recent multi-omics study suggested new molecular ENKTL subtypes according to mutation profiles including TP53 and MGA mutations and showed differing outcomes and survival rates among these subtypes (10). Thus, genetic evaluation of tumor biopsy samples prior to treatment may help predict the prognosis of ENKTL. However, a biopsy is not always possible, and the predictive value of archived tissue samples may be influenced by intra-patient heterogeneity, mainly according to the site and timing for tumor tissue sample collection.

Liquid biopsy using circulating tumor DNA (ctDNA) in plasma may be a proxy for the biological nature of tumor tissues. Recent advances in NGS technology can be applied to analyze ctDNA and facilitate cancer diagnosis, surveillance, and prognosis (11). Thus, ctDNA has been suggested as a biomarker, with several studies showing significant associations between ctDNA and relapse in diffuse large B-cell lymphoma (DLBCL), mantle cell lymphoma, and Hodgkin’s lymphoma (12–14).

In this study, we aimed to explore the clinical relevance of ctDNA at the time of diagnosis, relapse, and ENKTL progression using archived plasma samples from a prospective cohort study. The concordance of plasma ctDNA mutations with those of tumor biopsies was analyzed using 171 genes selected according to therapeutic, prognostic, and diagnostic properties (7, 15–18) (Additional file: Table S1). In addition, the effectiveness of ctDNA as a biomarker was compared with that of EBV DNA in patients with ENKTL using longitudinal blood samples serially collected during follow-up.

Materials and methods

Patients and study design

This study included patients with ENKTL who were prospectively enrolled at Samsung Medical Center (Seoul, South Korea) between March 2017 and December 2020 (prospective cohort study for lymphoma: NCT03117036). The inclusion criteria were radiological and pathological diagnoses of ENKTL according to the WHO classification of tumors with hematopoietic and lymphoid tissue criteria (19): 1) aged ≥18 years at the time of diagnosis; relapsed or newly diagnosed with primary tumor invasion and/or regional lymph node involvement, 2) willing to offer peripheral blood samples prior to treatment, during treatment, and during surveillance, and 3) agreement to provide written informed consent to participate in this study. A total of 73 patients fulfilled the inclusion criteria and were enrolled in the study. First, we evaluated the ctDNA analysis platform with 34 cross-sectional archived samples obtained from 28 patients at various disease states including remission (n = 22), relapse or progression (n = 7), partial response to therapy (n = 4), and at diagnosis (n = 1). Second, the utility of ctDNA testing was evaluated in the study cohort, consisting of 34 newly diagnosed and 11 relapsed patients. Serial blood samples collected from baseline (pre-treatment or patients in relapse prior to the state of second-line chemotherapy) to the time of completion of planned treatment, disease relapse, or progression were used to compare the efficacy of ctDNA testing to that of EBV DNA for prediction of treatment outcome (Figure 1A).

Figure 1

Clinicopathological features

Clinicopathological data of the patients were retrieved from electronic medical records and comprised age, sex, date of diagnosis, date of follow-up, prognostic index for natural killer cell lymphoma-EBV (PINK-E), International Prognostic Index score, Eastern Cooperative Oncology Group performance score, pathological information on ENKTL tissue, available NGS tissue biopsy results, and relevant laboratory findings including whole blood EBV quantitation and lactate dehydrogenase level. The Ann Arbor staging system was used for stage assessment. Imaging studies including CT, PET/CT, and magnetic resonance imaging were used to assess stage, extranodal involvement, and recurrence. The survival status of all enrolled patients was updated regularly, and the last update regarding survival and disease status was completed in January 2022.

Sequencing and data analysis

DNA was extracted from plasma and matched FFPE tissue samples. DNA from peripheral blood mononuclear cells (PBMC) was extracted using the QIAamp DNA Mini Kit (QIAGEN, Hilden, Germany), quantified using the TapeStation 4150 (Agilent Technologies, Santa Clara, CA, USA), and fragmented to a mean size of 250 base pairs using a sonicator (Diagnode, Lausanne, Switzerland). DNA in plasma was extracted using Magnetic Circulating DNA Maxi Reagent (Dxome, Sungnam, South Korea).

The size of ctDNA was measured using the TapeStation 4150. The concentration of isolated DNA was measured using a Qubit 4.0 Fluorometer (Thermo Fisher Scientific, Waltham, MA, USA). A total of 30 ng of plasma DNA was used for library preparation; in samples with low DNA quantities, the entire DNA extract was used. Library preparation was performed using the DxSeq Library prep reagent (Dxome, Sungnam, South Korea) according to the manufacturer’s instructions. Targeting an average sequencing depth of 35,000× per sample, sequencing was done on the NovaSeq 6000 platform (Illumina, San Diego, CA, USA). For each sample, PBMCs were sequenced as germline-matched controls using identical panel and library kits, targeting an average depth >2,500×.

Genotyping of ctDNA interpretation

We used the PiSeq and DxSeq softwares (Dxome) for variant identification and annotation, respectively, and additionally searched mutation databases and literature. Variants detected in both ctDNA and PBMCs were considered germline variants or variants of clonal hematopoiesis of indeterminate potential (CHIP). All somatic variants were categorized into Tier I/II, Tier III, or Tier IV according to the recommendation from AMP/ASCO/CAP guideline (20). All variants were visually inspected using the Integrative Genome Viewer (Broad Institute of MIT and Harvard, Cambridge, MA, USA) to exclude false positives. Owing to the lower sequencing coverage of PBMC (2,500×) compared to that of ctDNA (35,000×), potential CHIP variants may be missed in PBMC samples. Therefore, ctDNA-only variants identified in DNMT3A, TET2, and ASXL1 with allele frequencies below 1% were also considered to be CHIP. The variant allele frequency (VAF; %) was calculated as the percentage of sequencing reads of a specific DNA variant divided by the overall coverage at that locus.

Statistical analysis

Statistical analyses were performed using MedCalc version 18.2.1 (MedCalc Software; Mariakerke, Belgium). For continuous data, the Shapiro–Wilk test was used to detect departures from normality; variables were compared using the Mann–Whitney U test. The correlation between ctDNA concentration and EBV quantitation was determined using Spearman’s correlation analysis. The chi-square test or Fisher’s exact test was used to assess differences in the detection of ctDNA between groups. Progression-free survival (PFS) was measured from baseline (before the start of chemotherapy) to the date of progression, death, or the last follow-up. Overall survival (OS) was measured from the date of the initial sample draw prior to the start of chemotherapy to the date of death or the last follow-up. Survival analysis was performed using the Kaplan–Meier method. Statistical significance was defined as P < 0.05.

Results

ctDNA NGS panel evaluation for gene alterations detection

The concordance between tissue-based and ctDNA-based genotyping was analysed in 11 cases. Seven Tier I/II mutations were observed in both plasma and paired tissue, while another seven mutations were found only in tissue and six were found only in plasma (Figure 1B and Additional file: Table S2).

A total of 34 blood samples (average sequencing depth: 21,937.7×, Additional file: Table S3) from 28 patients at different time points were also examined. Most samples (33/34, 97.1%) were obtained during treatment or relapse. Oncogenic or druggable Tier I/II mutations were detected in 10 of 12 samples (83.3%) in patients with persistent disease or disease progression, while such mutations were not found in 22 samples from patients who achieved complete response. TP53 was the most frequently mutated gene (5/12, 42.0%) (Figure 1C). Altogether, these data indicate that our ctDNA analysis platform may be a useful tool for ENKTL diagnosis and monitoring.

Clinical characteristics of the study cohort

In the study cohort, 45 patients were enrolled. The median age at the time of enrolment was 55 years (range, 22–80 years); 64.4% (29/45) of patients were male (Table 1). The majority of patients was in stage I/II (71.1%), and 31 (68.9%) among them had a detectable level of EBV DNA in blood (Table 1). According to the PINK-E risk stratification, 21 (46.7%) patients belonged to intermediate- or high-risk groups (Table 1). Newly diagnosed patients with ENKTL at stage I/II involving the nasal area received concurrent chemo radiotherapy followed by two cycles of VIDL (etoposide, ifosfamide, dexamethasone, and L-asparaginase) chemotherapy. Advanced-stage patients received SMILE (steroid, methotrexate, ifosfamide, L-asparaginase, and etoposide) chemotherapy, whereas gemcitabine-containing chemotherapy or immune checkpoint inhibitors were used for relapsed or refractory disease.

Table 1

CharacteristicPatients, n (%)
SexMale29 (64.4)
Female16 (35.6)
Age (year, median age: 55, range: 22 - 80)≤ 6031 (68.9)
> 6014 (31.1)
Disease statusRelapsed11 (24.4)
Newly diagnosed34 (75.6)
ECOG PS0-143(95.6)
22 (4.4)
Ann Arbor StageI-II32 (71.1)
III-IV13 (28.9)
Primary tumor localizationNasal30 (66.7)
Non-nasal15 (33.3)
Bone marrow involvementPresence6 (15.6)
Absence39 (84.4)
B symptomsPresence7 (15.6)
Absence38 (84.4)
Serum LDHElevated26 (57.8)
Not elevated19 (42.2)
Blood EBV DNADetected31 (68.9)
Not detected14 (31.1)
Serum β2 microglobulin< 2·8 mg/L25 (55.6)
≥ 2·8 mg/L20 (44.4)
IPI riskLow/Low-Intermediate25 (55.6)/10 (22.2)
High-Intermediate/High7 (15.6)/3 (6.7)
PINK-E riskLow24 (53.3)
Intermediate/High21 (46.7)

Characteristics of 45 ENKTL patients in the study.

EBV, Epstein-Barr virus; ECOG PS, Eastern Cooperative Oncology Group performance status score; IPI, International Prognostic Index; LDH, lactate dehydrogenase; PINK-E, prognostic index for natural killer cell lymphoma–Epstein-Barr virus.

ctDNA concentrations as a tumor burden indicators in ENKTL

A total of 118 plasma samples including 45 at baseline, 56 during treatment and surveillance, and 17 at the time of progression were collected (average sequencing depth: 25,076.9×, Additional file: Table S3). The concentrations of total DNA isolated from plasma were highly variable (range, 0.39–511.30 ng/mL), and the median concentration of plasma DNA in pre-treatment or relapsed/refractory samples was 10.66 ng/mL, which was significantly higher than that in complete remission samples (4.71 ng/mL; Mann–Whitney test, P = 0.001). Compared to patients in early-stage disease (stage I/II), those in advanced-stage disease (stage III/IV) had higher plasma DNA concentration (26.55 vs. 7.97 ng/mL; P = 0.003) and blood EBV DNA titer (115,224 vs. 442 copies/mL; P = 0.002) at initial diagnosis prior to treatment (Figures 1D, E). Accordingly, the EBV DNA titer in blood was significantly associated with ctDNA concentration (R2 =0.558; P <0.0001).

Genotyping of ctDNA at diagnosis

Baseline ctDNA testing of initial samples from 34 newly diagnosed showed a detection rate of 54.5%. ctDNA mutations were more frequently observed in patients with advanced-stage disease than in those with early-stage disease (Figure 1F); the presence of a mutation did not correlate with detection of EBV DNA in blood (P=0.827). BCOR (21%) was the most commonly mutated gene, followed by TP53 (18%), DDX3X (15%), and KMT2D (12%) (Figure 1G). The median number of mutations per sample was 1.0 (range; 0–6). The median VAF of ctDNA mutations was 2.7% (interquartile range [IQR]; 0.3%–21.4%), and the majority (79.2%) had allele frequencies exceeding 1.0%.

From a matched analysis using PBMCs as germline-matched controls, 16 among the 34 newly diagnosed patients harbored 21 variants suspicious for CHIP (Additional file: Table S4). DNMT3A (61.9%, 13/21) was the most frequently mutated gene, followed by KMT2C (14.3%, 3/21). The median VAF of CHIP variants was 0.9%, ranging from 0.1% to 7.2%. Nine (56.3%) among the 16 patients with CHIP variants, harbored Tier I/II mutations other than CHIP. The remaining seven patients with only CHIP variants showed no disease progression.

Comparison of ctDNA mutations between newly diagnosed and relapsed patients

TP53, DDX3X, and BCOR presented the highest mutation frequencies in the newly diagnosed patient group. Of note, patients with advanced-stage disease had more DDX3X and BCOR mutations (P = 0.01 and < 0.001, respectively) than early-stage and relapsed patient groups; thus, at least one genetic variation in the TP53, DDX3X, or BCOR was identified in all patients with advanced-stage disease (Figure 2A).

Figure 2

In 15 patients with relapse, eight (53.4%) harbored Tier I/II mutations (Figure 2B). The detection rate was lower than that for newly diagnosed patients, which might be partly related to the poor quality of samples or low DNA amounts obtained from relapsed patients who were heavily pre-treated prior to sampling. Nevertheless, TP53 (33.0%) was the most commonly mutated gene, followed by POT1 (13.0%).

Feasibility of ctDNA monitoring in combination with EBV quantification

The 45 patients were classified as having advanced- or early-stage disease; the detection rates of baseline ctDNA mutations and EBV quantitation results in each patient were compared (Figure 3A). As expected, 84.6% (11/13) of patients in advanced stage were detected by both methods while EBV- and ctDNA mutation-positive rates in the early stage were low (19/32 and 12/32, respectively) with a low concordance. Whole blood EBV DNA was not detected in the baseline samples from 14 of 45 patients (31.3%), although, seven of 14 patients (one in stage IV and six in stage I/II) harbored ctDNA mutations that reflected their clinical course. This suggests that ctDNA genotyping can be a complementary method to EBV DNA quantification for disease monitoring, especially in EBV-negative patients. Among the 14 patients with EBV-positive and ctDNA mutation -negative results, 13 patients were at early stage, and the median EBV quantitation was 1021.5 copies/mL, which was approximately 1/10 of that for ctDNA-positive cases (median EBV quantitation, 14201 copies/mL).

Figure 3

In seven patients who had ctDNA mutation -positive results at baseline and achieved remission, none had mutations in plasma within 2–4 cycles of chemotherapy whereas EBV DNA was detected from three during monitoring (Figure 3B). This may suggest that EBV quantitation does not fully reflect patient’s disease status and ctDNA genotyping can be more specific.

Longitudinal plasma samples (n = 112) were serially collected from 39 patients in the cohort; Tier I/II mutations were identified in pre-treatment or relapsed/refractory samples (Additional file: Table S5) among 21 patients. A total of 14 patients experienced disease progression or relapse, and the majority (10/14, 71.4%) showed persistently detectable or emerging ctDNA mutations. An exemplary case with disease involving the liver and spleen showed an increase in VAFs of mutations as well as an increase of whole blood EBV DNA titer as the disease progressed (Additional file: Figure S1A). On the other hand, a case with localized disease involving the nasal cavity showed an absence of mutations and a low level of detectable EBV DNA; the VAFs of newly emerging mutations and EBV DNA titer were elevated as new lesions in the spleen appeared (Additional file: Figure S1B). Another exemplary case involving the nasopharynx showed a low VAF at diagnosis and the resolution of mutations after treatment. However, the patient had new intra-abdominal lesions without any lesions in the nasopharynx, with an increase of EBV DNA titer, re-appearance of initial mutations and emergence of new mutations (Figure 4A). Many cases showed frequent changes in the composition of mutations with each genetic mutation indicating the occurrence of clonal evolution (Figures 4B, C).

Figure 4

Prognostic impact of ctDNA mutations

Of the 34 newly diagnosed patients, one was excluded due to follow-up loss immediately after diagnosis, and the remaining 33 patients had a median follow-up period of 21.4 months (IQR; 4.7–25.1 months). Survival analysis showed that DDX3X and BCOR mutations tended to be associated with inferior OS, although this was not statistically significant (Additional file: Figure S2). By contrast, the median PFS of patients with DDX3X mutations (2.2 months) was shorter than that for patients with other gene mutations (25.6 months) or no mutation (25.4 months) (log-rank test; P < 0.0001, Figure 5).

Figure 5

Multivariable Cox regression analysis demonstrated an independent prognostic value of DDX3X mutation for PFS (hazard ratio [HR], 8.256; 95% confidence interval [95%CI], 1.371- 49.722; P = 0.021, Table 2).

Table 2

PFSOS
UnivariateMultivariate*UnivariateMultivariate*
POR95%CIPPOR95%CIP
Age (years)0.618–––0.272–––
BCOR mutation0.477–––0.144––0.288
DDX3X mutation<0.0018.2561.371 to 49.7220.0210.228–––
TP53 mutation0.742–––0.746–––
EBV DNA quantitation (copies/mL)0.201–––0.042––0.162
PINKE: 1 vs ≥20.00611.791.309 to 106.1470.0280.04314.8371.861 to 118.280.011

Analysis of prognostic factors for ENKTL patients.

*The variables that entered the multivariate Cox regression model included variables identified by univariate analysis with P <0.15.

PFS, progression-free survival; HR, hazard ratio; CI, confidence interval; PINK, prognostic index of natural killer lymphoma; PINK-E, prognostic index for natural killer cell lymphoma–Epstein-Barr virus; OS, overall survival.

Discussion

The aggressiveness of ENKTL could differ according to the involved sites as well as the EBV DNA titer in blood; patients with tumors in non-nasal areas and high level of EBV DNA at diagnosis frequently show an aggressive clinical course, resulting in poor outcomes even in cases of localized disease. Therefore, an assessment of pre-treatment risk may be fundamental and influence the outcome of primary treatments, especially when most patients show unfavorable outcomes if the first-line treatment fails. Fluorodeoxyglucose (FDG)-PET/CT has been evaluated as a tool for risk prediction in lymphoma; however, its value is limited because ENKTL tumor cells are generally mixed with inflammatory cells, and combined inflammation may increase FDG uptake (21, 22). Furthermore, coagulative necrosis frequently found in ENKTL may also decrease FDG uptake. By contrast, EBV-DNA level in blood reflects the tumor burden and can be used for disease monitoring (6, 23, 24). However, EBV infection is highly prevalent in healthy individuals (25), and EBV DNA is frequently detected in elderly people (26). Therefore, EBV DNA level in blood may not be accurate for estimation of tumor burden in ENKTL patients. Although a previous study suggested that the Deauville score, reflecting FDG uptake and EBV DNA in blood, at the end of treatment could predict risk of relapse in ENKTL, its value is not sufficient to estimate treatment strategies for ENKTL (27).

In this study, we evaluated the feasibility and performance of in-depth cell-free ctDNA sequencing as a tool for prognostication and monitoring of ENKTL. We used a sensitive customized ctDNA NGS panel covering 171 genes, allowing detection of mutations with very low VAF (0.25%). The number of genes in the panel was larger and the detection sensitivity was higher than those in previous studies performed in NK or T-cell lymphomas (28, 29). The high-sensitivity of our platform allowed the detection of mutations repeatedly reported as main genetic alterations in ENKTL such as those in TP53, KRAS, EP300, and MGA, even those existing in very low VAFs.

Previous studies monitoring ENKTL patients have reported that clonal changes were as frequent as in other types of lymphoma (29, 30), suggesting that genomic profiling using proper multi-gene panels is needed both for initial diagnosis as well as monitoring. In this study, the three most frequently mutated genes (BCOR, TP53, and DDX3X) were found in all patients with advanced disease, in agreement with previous studies that have reported that mutations in DDX3X and TP53 were adverse prognostic markers in ENKTL (10). DDX3X, a RNA helicase family gene, and TP53, a tumor suppressor gene, are assumed to have a close relationship with the biological process of ENKTL (8). Loss-of-function mutations of DDX3X are also frequent in B cell lymphomas including Burkitt lymphoma and DLBCL (31). DDX3X mutants exhibited loss of suppressive effects on cell-cycle progression in NK cells and transcriptional activation of pro-proliferative pathways in ENKTL (8, 32). Mutations of DDX3X is an adverse prognostic marker in CHOP-treated ENKTL patients, but this is not the case in ENKTL patients treated with asparaginase (ASP)-based regimens (10). In our study, all patients were treated with ASP-based regimens, yet the DDX3X mutation showed adverse effect on PFS in this study.

The TP53 mutation was more frequent in relapsed patients (33.0%) than in newly diagnosed patients (18.0%). Considering that clonal evolution is frequent in lymphoma, most TP53 mutations detected in relapsed/refractory patients may be from emerging clones during treatment that contribute to tumor relapse, as reported by previous studies in other lymphomas (33, 34). However, further studies including a larger number of serial samples from patients with ENKTL are needed to verify the clonal evolution of TP53.

Herein, we demonstrated that ctDNA genotyping is a useful tool to diagnose ENKTL, estimate tumor burden, monitor treatment response, and detect emerging mutations of residual disease or relapse, and therefore, the assay may be suggested as an option for patient care in combination with conventional tools (Figure 6). Current diagnostic tools for ENKTL include imaging studies and EBV DNA level in blood. However, even if the treatment was very effective in imaging studies, uptake alterations sometimes remained and EBV DNA level could increase without clinical evidence of residuals or progression of ENKTL (26, 35, 36). Therefore, when the clinical features and results of the conventional monitoring tools conflict, a ctDNA analysis may help to estimate the patient’s condition and make therapeutic decision. Specifically, we found that the ctDNA genotyping could be used for tumor burden monitoring in EBV-negative patients. After chemotherapy or concurrent chemoradiotherapy (CCRT) treatments, ctDNA analysis may be performed along with imaging and EBV DNA quantitation. In patients showing discrepancies between imaging and EBV results, ctDNA genotyping may provide additional information for medical decisions. In such cases, if mutation is not detected in plasma sample (negative ctDNA), surveillance monitoring without additional chemotherapy and other causes of elevated EBV could be considered. However, if mutation is detected in plasma sample (positive ctDNA), a short-term follow-up or additional secondary therapy could be planned due to the possibility of residual ENKTL. By contrast, in patients with relapse, ctDNA analysis may indicate the necessity for a change in treatment regimen and allow more accurate prognoses based on changes in the mutation spectrum.

Figure 6

There were some limitations in our study. First, some of the tissue samples were analyzed by panels that targeted different genes; therefore, only genes targeted by both panels could be compared. Second, EBV DNA quantitation measurement was done in the whole blood samples. Studies reported that the plasma is a preferable sample than whole blood in measuring EBV DNA as a marker of EBV diseases (6, 37, 38). However, comprehensive recommendations and international guidelines are still lacking (39), and the institute where this study was conducted used whole blood samples for EBV quantitation in clinical practice that we used the whole blood sample. Thus, tumor burden analysis of ctDNA mutation and plasma EBV DNA quantitation might be needed in the future study. Third, although this study enrolled a larger number of patients with longitudinal samples than previous studies (16, 29), the sample size may still be too small to establish a firm conclusion. The topic should be further explored through a larger, prospective clinical trial. Nevertheless, ctDNA mutations showed dynamic properties during treatment, suggesting that they could be excellent indicators of disease course.

In conclusion, ctDNA testing at diagnosis may confirm the presence of a mutation that can be subsequently traced in blood samples during monitoring in patients with ENKTL. Testing it could act as an excellent marker providing additional information compared to that by imaging and EBV DNA quantitation and contribute to making effective medical decisions.

Statements

Data availability statement

The original contributions presented in the study are included in the article/Supplementary Material. Further inquiries can be directed to the author ().

Ethics statement

The studies involving human participants were reviewed and approved by Ethics Committee of Samsung Medical Center. The patients/participants provided their written informed consent to participate in this study.

Author contributions

JK and H-YK designed the study, interpreted data, and wrote the manuscript; ZC, SH, and HJ performed the molecular studies; WK, SK, H-YK interpreted data; SY, collected data; SS designed the study and contributed to manuscript revision; SK, JC, and S-TL provided study materials and contributed to manuscript revision.

Funding

This research was supported by grants provided by National Research Foundation of Korea funded by the Korean government (NRF-2021R1A2C1007531 and NRF-2017R1D1A1B04028149).

Acknowledgments

The authors thank Medical Illustration & Design, part of the Medical Research Support Services of Yonsei University College of Medicine, for all artistic support related to this work.

Conflict of interest

ZC, SH, HJ, JC, and S-TL are employees of Dxome co., Ltd.

The remaining authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest

Publisher’s note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Supplementary material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fonc.2023.1109715/full#supplementary-material

References

  • 1

    Montes-MojarroIAFendFQuintanilla-MartinezL. EBV and the pathogenesis of NK/T cell lymphoma. Cancers (Basel) (2021) 13(6):1414-32. doi: 10.3390/cancers13061414

  • 2

    WangLLiLRZhangLWangJW. The landscape of new drugs in extranodal NK/T-cell lymphoma. Cancer Treat Rev (2020) 89:102065. doi: 10.1016/j.ctrv.2020.102065

  • 3

    FoxCPCivalleroMKoYHManniMSkrypetsTPileriSet al. Survival outcomes of patients with extranodal natural-killer T-cell lymphoma: a prospective cohort study from the international T-cell project. Lancet Haematol (2020) 7(4):e284–e94. doi: 10.1016/S2352-3026(19)30283-2

  • 4

    TrotmanJLuminariSBoussettaSVersariADupuisJTychyjCet al. Prognostic value of PET-CT after first-line therapy in patients with follicular lymphoma: a pooled analysis of central scan review in three multicentre studies. Lancet Haematol (2014) 1(1):e17–27. doi: 10.1016/S2352-3026(14)70008-0

  • 5

    ItoYKimuraHMaedaYHashimotoCIshidaFIzutsuKet al. Pretreatment EBV-DNA copy number is predictive of response and toxicities to SMILE chemotherapy for extranodal NK/T-cell lymphoma, nasal type. Clin Cancer Res (2012) 18(15):4183–90. doi: 10.1158/1078-0432.CCR-12-1064

  • 6

    SuzukiRYamaguchiMIzutsuKYamamotoGTakadaKHarabuchiYet al. Prospective measurement of Epstein-Barr virus-DNA in plasma and peripheral blood mononuclear cells of extranodal NK/T-cell lymphoma, nasal type. Blood (2011) 118(23):6018–22. doi: 10.1182/blood-2011-05-354142

  • 7

    ChoiSGoJHKimEKLeeHLeeWMChoCSet al. Mutational analysis of extranodal NK/T-cell lymphoma using targeted sequencing with a comprehensive cancer panel. Genomics Inform (2016) 14(3):78–84. doi: 10.5808/GI.2016.14.3.78

  • 8

    JiangLGuZHYanZXZhaoXXieYYZhangZGet al. Exome sequencing identifies somatic mutations of DDX3X in natural killer/T-cell lymphoma. Nat Genet (2015) 47(9):1061–6. doi: 10.1038/ng.3358

  • 9

    LeeSParkHYKangSYKimSJHwangJLeeSet al. Genetic alterations of JAK/STAT cascade and histone modification in extranodal NK/T-cell lymphoma nasal type. Oncotarget (2015) 6(19):17764. doi: 10.18632/oncotarget.3776

  • 10

    XiongJCuiBWWangNDaiYTZhangHWangCFet al. Genomic and transcriptomic characterization of natural killer T cell lymphoma. Cancer Cell (2020) 37(3):403–19.e6. doi: 10.1016/j.ccell.2020.02.005

  • 11

    DiazLAJr.BardelliA. Liquid biopsies: genotyping circulating tumor DNA. J Clin Oncol (2014) 32(6):579–86. doi: 10.1200/JCO.2012.45.2011

  • 12

    CamusVViennotMLequesneJViaillyPJBohersEBessiLet al. Targeted genotyping of circulating tumor DNA for classical Hodgkin lymphoma monitoring: a prospective study. Haematologica (2021) 106(1):154–62. doi: 10.3324/haematol.2019.237719

  • 13

    LakhotiaRMelaniCDunleavyKPittalugaSSabaNLindenbergLet al. Circulating tumor DNA predicts therapeutic outcome in mantle cell lymphoma. Blood Adv (2022) 6(8):2667–80. doi: 10.1182/bloodadvances.2021006397

  • 14

    HerreraAFTracySCroftBOpatSRayJLovejoyAFet al. Risk profiling of patients with relapsed/refractory diffuse large b-cell lymphoma by measuring circulating tumor DNA. Blood Adv (2022) 6(6):1651–60. doi: 10.1182/bloodadvances.2021006415

  • 15

    Arzuaga-MendezJPrieto-FernandezELopez-LopezEMartin-GuerreroIGarcia-RuizJCGarcia-OradA. Cell-free DNA as a biomarker in diffuse large b-cell lymphoma: A systematic review. Crit Rev Oncol Hematol (2019) 139:7–15. doi: 10.1016/j.critrevonc.2019.04.013

  • 16

    LiQZhangWLiJXiongJLiuJChenTet al. Plasma circulating tumor DNA assessment reveals KMT2D as a potential poor prognostic factor in extranodal NK/T-cell lymphoma. biomark Res (2020) 8:27. doi: 10.1186/s40364-020-00205-4

  • 17

    de MelSSoonGSMokYChungTHJeyasekharanADChngWJet al. The genomics and molecular biology of natural Killer/T-cell lymphoma: Opportunities for translation. Int J Mol Sci (2018) 19(7):1931. doi: 10.3390/ijms19071931

  • 18

    WatataniYSatoYMiyoshiHSakamotoKNishidaKGionYet al. Molecular heterogeneity in peripheral T-cell lymphoma, not otherwise specified revealed by comprehensive genetic profiling. Leukemia (2019) 33(12):2867–83. doi: 10.1038/s41375-019-0473-1

  • 19

    SwerdlowSHCampoEHarrisNLJaffeESPileriSASteinHet al. WHO classification of tumours of haematopoietic and lymphoid tissues. 4th edition. Lyon, France:International Agency of Research on Cancer (2017).

  • 20

    LiMMDattoMDuncavageEJKulkarniSLindemanNIRoySet al. Standards and guidelines for the interpretation and reporting of sequence variants in cancer: A joint consensus recommendation of the association for molecular pathology, American society of clinical oncology, and college of American pathologists. J Mol Diagn (2017) 19(1):4–23. doi: 10.1016/j.jmoldx.2016.10.002

  • 21

    KhandaniAHDunphyCHMeteesatienPDufaultDLIvanovicMSheaTC. Glut1 and Glut3 expression in lymphoma and their association with tumor intensity on 18F-fluorodeoxyglucose positron emission tomography. Nucl Med Commun (2009) 30(8):594–601. doi: 10.1097/MNM.0b013e32832cc295

  • 22

    GuoRXuPChengSLinMZhongHLiWet al. Comparison of nasopharyngeal MR, 18 f-FDG PET/CT, and 18 f-FDG PET/MR for local detection of natural Killer/T-cell lymphoma, nasal type. Front Oncol (2020) 10:576409. doi: 10.3389/fonc.2020.576409

  • 23

    HaverkosBMPanZGruAAFreudAGRabinovitchRXu-WelliverMet al. Extranodal NK/T cell lymphoma, nasal type (ENKTL-NT): An update on epidemiology, clinical presentation, and natural history in north American and European cases. Curr Hematol Malig Rep (2016) 11(6):514–27. doi: 10.1007/s11899-016-0355-9

  • 24

    LeiKIChanLYChanW-YJohnsonPJLoYD. Diagnostic and prognostic implications of circulating cell-free Epstein-Barr virus DNA in natural killer/T-cell lymphoma. Clin Cancer Res (2002) 8(1):29–34.

  • 25

    DowdJBPalermoTBriteJMcDadeTWAielloA. Seroprevalence of Epstein-Barr virus infection in US children ages 6-19, 2003-2010. PloS One (2013) 8(5):e64921. doi: 10.1371/journal.pone.0064921

  • 26

    TrottierHButeauCRobitailleNDuvalMTucciMLacroixJet al. Transfusion-related Epstein-Barr virus infection among stem cell transplant recipients: a retrospective cohort study in children. Transfusion (2012) 52(12):2653–63. doi: 10.1111/j.1537-2995.2012.03611.x

  • 27

    KimSJChoiJYHyunSHKiCSOhDAhnYCet al. Risk stratification on the basis of deauville score on PET-CT and the presence of Epstein-Barr virus DNA after completion of primary treatment for extranodal natural killer/T-cell lymphoma, nasal type: a multicentre, retrospective analysis. Lancet Haematol (2015) 2(2):e66–74. doi: 10.1016/S2352-3026(15)00002-2

  • 28

    KimSJKimYJYoonSERyuKJParkBParkDet al. Circulating tumor DNA-based genotyping and monitoring for predicting disease relapses of patients with peripheral T-cell lymphomas. Cancer Res Treat (2023) 55(1):291–303. doi: 10.4143/crt.2022.017

  • 29

    QiFCaoZChenBChaiYLinJYeJet al. Liquid biopsy in extranodal NK/T-cell lymphoma: a prospective analysis of cell-free DNA genotyping and monitoring. Blood Adv (2021) 5(11):2505–14. doi: 10.1182/bloodadvances.2020001637

  • 30

    Chang-TsuYTai-ChungHChi-ChenCYi-KuangCJia-HueiTJau-YuLet al. Clonal sequence tracking reveals TET2-mutated extranodal NK/T-cell lymphoma disseminated independent of Epstein Barr virus. Haematologica (2019) 104(10):e489–e92. doi: 10.3324/haematol.2018.215327

  • 31

    GongCKrupkaJAGaoJGrigoropoulosNFGiotopoulosGAsbyRet al. Sequential inverse dysregulation of the RNA helicases DDX3X and DDX3Y facilitates MYC-driven lymphomagenesis. Mol Cell (2021) 81(19):4059–75 e11. doi: 10.1016/j.molcel.2021.07.041

  • 32

    de MelSHueSSJeyasekharanADChngWJNgSB. Molecular pathogenic pathways in extranodal NK/T cell lymphoma. J Hematol Oncol (2019) 12(1):33. doi: 10.1186/s13045-019-0716-7

  • 33

    EliasCFlorenceCPaoloGUlrichJSarkaPRichardRet al. TP53 aberrations in chronic lymphocytic leukemia: an overview of the clinical implications of improved diagnostics. Haematologica (2018) 103(12):1956–68. doi: 10.3324/haematol.2018.187583

  • 34

    BerendsenMRStevensWBCvan den BrandMvan KriekenJHScheijenB. Molecular genetics of relapsed diffuse Large b-cell lymphoma: Insight into mechanisms of therapy resistance. Cancers (Basel) (2020) 12(12):3553-78. doi: 10.3390/cancers12123553

  • 35

    XueW-QWangT-MHuangJ-WZhangJ-BHeY-QWuZ-Yet al. A comprehensive analysis of genetic diversity of EBV reveals potential high-risk subtypes associated with nasopharyngeal carcinoma in China. Virus Evol (2021) 7(1):veab010. doi: 10.1093/ve/veab010

  • 36

    JahrSHentzeHEnglischSHardtDFackelmayerFOHeschR-Det al. DNA Fragments in the blood plasma of cancer patients: quantitations and evidence for their origin from apoptotic and necrotic cells. Cancer Res (2001) 61(4):1659–65.

  • 37

    HaJYChoHSungHJungARLeeYSLeeSWet al. Superiority of Epstein-Barr virus DNA in the plasma over whole blood for prognostication of extranodal NK/T cell lymphoma. Front Oncol (2020) 10:594692. doi: 10.3389/fonc.2020.594692

  • 38

    KanakryJAHegdeAMDurandCMMassieABGreerAEAmbinderRFet al. The clinical significance of EBV DNA in the plasma and peripheral blood mononuclear cells of patients with or without EBV diseases. Blood (2016) 127(16):2007–17. doi: 10.1182/blood-2015-09-672030

  • 39

    LeeAWMLeeVHFNgWTStrojanPSabaNFRinaldoAet al. A systematic review and recommendations on the use of plasma EBV DNA for nasopharyngeal carcinoma. Eur J Cancer (2021) 153:109–22. doi: 10.1016/j.ejca.2021.05.022

Summary

Keywords

liquid biopsy, Extranodal natural killer/T cell lymphoma, biomarker, CtDNA, molecular biomarker

Citation

Kim JJ, Kim H-Y, Choi Z, Hwang S, Jeong H, Choi JR, Yoon SE, Kim WS, Kim S-H, Kim H-J, Shin S-Y, Lee S-T and Kim SJ (2023) In-depth circulating tumor DNA sequencing for prognostication and monitoring in natural killer/T-cell lymphomas. Front. Oncol. 13:1109715. doi: 10.3389/fonc.2023.1109715

Received

28 November 2022

Accepted

24 January 2023

Published

10 February 2023

Volume

13 - 2023

Edited by

Sylvain Choquet, Hôpitaux Universitaires Pitié Salpêtrière, France

Reviewed by

Mir Basharath Alikhan, NorthShore University HealthSystem, United States; Arnaud Jaccard, Centre Hospitalier Universitaire de Limoges, France

Updates

Copyright

*Correspondence: Sang-Yong Shin, ; Seung-Tae Lee, ; Seok Jin Kim,

†These authors have contributed equally to this work

This article was submitted to Hematologic Malignancies, a section of the journal Frontiers in Oncology

Disclaimer

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.

Outline

Figures

Cite article

Copy to clipboard


Export citation file


Share article

Article metrics