Recurrence Dynamics of Pathological N2 Non–Small Cell Lung Cancer Based on IASLC Residual Tumor Descriptor
Article information
Abstract
Purpose
This study investigated the recurrence patterns and timing in patients with pathologic N2 (pN2) non-small cell lung cancer (NSCLC) according to the residual tumor (R) descriptor proposed by the International Association for the Study of Lung Cancer (IASLC).
Materials and Methods
From 2004 to 2021, patients with pN2 NSCLC who underwent anatomical resection were analyzed according to the IASLC R criteria using medical records from a single center. Survival analysis was performed using Cox proportional hazards models. Recurrence patterns between complete (R0) and uncertain resections (R[un]) were compared.
Results
In total, 1,373 patients were enrolled in this study: 576 (42.0%) in R0, 286 (20.8%) in R(un), and 511 (37.2%) in R1/R2 according to the IASLC R criteria. The most common reason for R(un) classification was positivity for the highest lymph node (88.8%). In multivariable analysis, the hazard ratios for recurrence in R(un) and R1/R2 compared to R0 were 1.18 (95% confidence interval [CI], 0.96–1.46) and 1.58 (1.31–1.90), respectively. The hazard rate curves displayed similar patterns among groups, peaking at approximately 12 months after surgery. There was a significant difference in distant recurrence patterns between R0 and R(un). Further analysis after stratification with the IASLC N2 descriptor showed significant differences in distant recurrence patterns between R0 and R(un) in patients pN2a1 and pN2a2 disease, but not in those with pN2b disease.
Conclusion
The IASLC R criteria has prognostic relevance in patients with pN2 NSCLC. R(un) is a highly heterogeneous group, and the involvement of the highest mediastinal lymph node can affect distant recurrence patterns.
Introduction
Complete resection (R0) of tumor is one of the most important prognostic factors in the surgical treatment of non–small cell lung cancer (NSCLC). The Union for International Cancer Control (UICC) classification for residual tumors (R) provides an intuitive representation of the completeness of surgical resection and can influence subsequent treatment decisions [1]. However, since the UICC R descriptor is solely determined based on the presence of tumor at the surgical margin, other crucial factors for complete resection, such as the degree of lymph node (LN) evaluation and extranodal extension, are omitted [2-4].
In this situation, the Staging and Prognostic Factors Committee of the International Association for the Study of Lung Cancer (IASLC) recently proposed a new category of R classification based on the expert opinion in 2005 [2], named “uncertain resection [R(un)]” [5,6]. This aimed to provide a more refined classification of resection status, bridging the gap between complete and incomplete resections. To date, the IASLC R classification for survival has been validated through multi-center– or single-center–based studies, demonstrating its prognostic relevance [5-9]. However, these studies have limitations as they lack detailed assessment of recurrence dynamics such as recurrence patterns, timing, and hazard rates according to the R descriptor. The R classification proposed by the IASLC was initially introduced because clinicians continued to encounter ambiguous recurrences despite UICC R0 resection. However, detailed information on cancer recurrence remains unreported. Recently, it has been reported that the recurrence dynamics in NSCLC patients may differ depending on various factors, such as histology and sex [10]. Therefore, to evaluate the prognostic relevance of the R descriptor, it is necessary to understand the recurrence dynamics associated with it.
In this study, we aimed to analyze recurrence dynamics in patients with pathologic N2 (pN2) NSCLC according to the IASLC R classification.
Materials and Methods
1. Patients
We included patients with pathologically confirmed N2 NSCLC who underwent curative intent surgery at the Asan Medical Center, Seoul, South Korea, between January 2004 and December 2021. Clinical information of included patients was extracted from our prospectively collected database. The data of all patients who underwent surgery at our institution were stored in electronic medical records, which were regularly updated by experienced researchers. The exclusion criteria of this study were (1) patients who are either elderly (> 80 years old) or very young (< 18 years old), (2) those with concurrent malignancies, (3) received neoadjuvant treatment before surgery, (4) had wedge resection, (5) presence of distant metastasis, and (6) mortality within 30 days (Fig. 1).
2. Staging workup and treatment
Diagnosis, clinical staging, and surgical strategy were described in a previous study [8]. In summary, the workup for diagnosis and clinical staging was based on imaging modalities such as chest computed tomography (CT), positron emission tomography computed tomography (PET-CT), and brain magnetic resonance imaging. If clinical N2 NSCLC was suspected, mediastinal LN biopsy was selectively performed via mediastinoscopy, endobronchial ultrasound (EBUS), or endoscopic ultrasound. If biopsy-proven N2 metastasis was confirmed, a multidisciplinary team comprising oncologists, pathologists, radiologists, radiation oncologists, and surgeons decided on the treatment approach. If a single N2 metastasis is suspected, upfront surgery is preferred when complete surgical resection can be achieved without the need for pneumonectomy and in the absence of clinical extracapsular extension (ECE). The staging of the patients included in this study was retrospectively determined based on the eighth edition of the American Joint Committee on Cancer Staging Manual [11].
Surgical treatment for NSCLC included the anatomical resection of pulmonary parenchyma and systematic mediastinal lymph node dissection (LND), which refers to the en bloc removal of the ipsilateral mediastinal and hilar LNs, including the surrounding adipose tissue. Specifically, LND includes the 2R, 4R, 7, 8, 9, 10R, 11R, and 12 nodes for right-sided NSCLC, and the 5, 6, 7, 8, 9, 10L, 11L, 12 and selectively 4L (if metastasis was suspected) nodes for left-sided NSCLC. All LNs were designated by the operator during surgery. LN stations were classified according to their anatomical location according to the Mountain-Dresler correction of the American Thoracic Society (until 2008) [12] and the IASLC LN map (since 2008) [13].
Adjuvant chemotherapy was recommended for all patients with pathologically confirmed N2 NSCLC. However, patients with poor performance status or those who expressed no intention of receiving treatment were excluded after consultation with the multidisciplinary team. Systematic adjuvant chemotherapy consisted of a platinum-based regimen administered for four cycles starting from 4 to 6 weeks after surgery. Targeted therapy was widely used as an additional treatment option in cases of recurrence or tumor progression after adjuvant treatment. For adjuvant radiation therapy, a dose of 1.8 Gy per day was administered. Patients who underwent complete resection received a total dose of 50.4 Gy, and patients with positive margins received 55-60 Gy. Patients were considered receiving adjuvant therapy only when the full course of treatment was completed.
3. Follow-up and recurrence
Outpatient-based follow-up was performed at intervals of 3 months for the first 2 years after surgery, every 6 months for the following 3 years, and annually thereafter with chest CT. When recurrence was suspected during surveillance, a systemic imaging workup, including PET-CT, was performed. The date of recurrence was defined as the date of initial detection on the imaging study, and the site of initial recurrence was recorded. After a pathological or clinical diagnosis of recurrence, a treatment plan was established by a multidisciplinary team.
4. Definition of variables
The UICC R classification, which consists of R0 resection with no residual tumor, R1 with microscopic residual tumor, and R2 with macroscopic residual tumor, has demonstrated prognostic value [1]. Subsequently, the IASLC R classification introduced an R(un) between R0 and R1/2, which was validated by several studies. According to the IASLC R classification, for R0 resection, there must be no residual tumor at the resection margin and all of the following: (1) sufficient LN evaluation must be performed (≥ 3 N1 and N2 stations, ≥ 3 N1 and N2 LNs, and LN 7 station evaluated), (2) the highest mediastinal LN must be negative for malignancy, (3) no carcinoma in situ on the bronchial stump, and (4) the pleural lavage test must be negative for malignancy. If any of these conditions is not satisfied, it is classified as R(un) [2]. R1 resection is defined as the presence of microscopic invasion at the resection margin, ECE of LNs at the resection margin of the lung specimen, or those removed separately, as well as malignant pleural or pericardial effusions. R2 resection is defined as the presence of macroscopic remaining tumor after resection, or lack of resection of involved LNs.
The highest mediastinal LN station was defined as the lowest numerically numbered mediastinal LN station evaluated. While systematic LN dissection was routinely performed in our institution, routine dissection of the 4L LN in the case of left-sided NSCLC was not performed. Therefore, when evaluating the 4L LN, 4L was defined as the highest mediastinal LN station, and when not evaluating the 4L LN, 5 was designated. There was ongoing debate regarding the hierarchical order of LN stations due to the variability in LN drainage patterns [14]. However, in most of studies, the lowest numerically numbered LN station was defined as the highest station [5,15].
According to the IASLC recommendation, N2 disease is subclassified into N2a1, N2a2, and N2b. N2a1 is defined as a single N2 metastasis without N1 metastasis (skip metastasis), N2a2 is a single N2 metastasis with N1 metastases, and N2b indicates multiple N2 metastases [16].
The histologic grade was classified according to tumor differentiation and consisted of grade I (well differentiated), grade II (moderately differentiated), and grade III (poorly differentiated), and cases where the grade could not be assessed.
Although the concept of recurrence may be somewhat inconsistent in patients who underwent incomplete resection, in this study, disease progression in patients who underwent incomplete resection was included in the recurrence category.
Locoregional recurrence was defined as recurrence at the surgical site, bronchial stump, ipsilateral lung parenchyma, or mediastinal and hilar LNs. Distant metastasis was referred to as recurrence that extends beyond the locoregional site including metastasis to the pleura, pericardium, other LNs, or distant organs. If locoregional and distant recurrence were detected simultaneously at the time of diagnosis, it was defined as mixed recurrence. Second primary lung cancer was defined according to the criteria of Martini and Melamed [17].
Disease-free survival (DFS) was defined as the time duration from the date of surgery to the date of recurrence, disease progression, death of any cause, or last follow-up. Overall survival (OS) was defined as the time duration from the date of surgery to the date of death of any cause or last follow-up.
5. Statistical analysis
Continuous variables are shown as means and standard deviations, and categorical variables are shown as frequencies and percentages. To compare between IASLC R groups (R0, R[un], and R1/2), analysis of variance for continuous variables and the chi-square tests for categorical variables were used. For comparisons between two groups, Student’s t test or the Mann-Whitney U test for continuous variables and the chi-square tests for categorical variables were applied. DFS and OS were calculated using the Kaplan–Meier method to describe the time to recurrence and death, respectively; these were compared using the log-rank test. Bonferroni correction was used to calculate the p-values in the log-rank test for multiple comparisons of the survival curves. The hazard ratio (HR) of recurrence was evaluated using the life table method, which was the conditional probability of manifesting recurrence within a certain time interval. To display the underlying patterns easier, some instability owing to random variation in the hazard rate estimates was dealt with using the kernel smoothing method. Univariable and multivariable Cox proportional hazards models were applied to identify the prognostic factors for recurrence. After excluding the correlated variables, those with p < 0.05 in the univariable analyses were included in the initial multivariable Cox model. The final multivariable model was selected using the forward stepwise selection method (p ≤ 0.10 for entering the model and p ≤ 0.05 for remaining in the model).
All statistical analyses were performed using R ver. 4.0.2 (The R Foundation for Statistical Computing, Vienna, Austria). p < 0.05 was considered statistically significant.
Results
1. Patient characteristics
A total of 1,373 patients were included in this study. Their mean age was 61.7±10.0 years, and 862 were male (62.8%). According to the UICC R classification, 1,250 (91.0%), 112 (8.2%), and 11 (0.8%) patients were classified as R0, R1, and R2, respectively. According to the IASLC R classification, 576 (42.0%), 286 (20.8%), and 511 (37.2%) were reclassified as R0, R(un), and R1/2, respectively. Among the UICC R0 patients, 576 (46.1%) patients were still IASLC R0 status, while 279 (22.3%) and 395 (31.6%) patients were reclassified as R(un) and R1/2, respectively (S1 Table). Of the 286 patients classified as R(un), 279 were reclassified from UICC R0 and seven were reclassified from UICC R1. Among the patients classified as IASLC R(un) from UICC R0, 25 (9.0%) patients had insufficient LN evaluation, including 19 patients (6.8%) with fewer than three N1 or N2 LN stations and six patients (2.2%) without evaluation of the LN 7 station. The remaining 254 patients (91.0%) were positive for malignancy in the highest mediastinal LN station (S2 Table). Among the patients classified as UICC R1, seven (6.2%) were reassigned to R(un) because of carcinoma in situ on the bronchial resection margin. Furthermore, 395 patients from UICC R0 were reclassified as IASLC R1/2 status because of ECE (S2 Table).
The baseline characteristics of the included patients are shown in Table 1. There were no differences in age (p=0.398), sex (p=0.623), smoking history (p=0.451), comorbidities (p=0.594), forced expiratory volume in 1 second (p=0.650), and diffusing capacity to carbon monoxide (p=0.322) between groups. In addition, there were no differences in the frequency of performing EBUS (p=0.424) and mediastinoscopy (p=0.762), histologic subtype (p=0.659), tumor size (p=0.073), and adjuvant therapy (p=0.095). However, in the R(un), left-sided NSCLC was more common (p=0.001), and a significantly lower number of harvested LNs was observed (p=0.001).
2. Progression outcomes
The mean follow-up duration was 52.88±44.13 months, during which there were 710 cases of recurrence (51.7%) and 667 deaths (48.6%) (S3 Table). According to the IASLC R classification, the 5-year DFS rates were 43.2% for R0, 34.2% for R(un), and 21.1% for R1/2 (p < 0.001). There were significant differences between R0 and R(un) (p < 0.001), as well as between R(un) and R1/2 (p=0.004) (Fig. 2A). The 5-year OS rates were 62.4% for R0, 54.6% for R(un), and 43.4% for R1/2 (p < 0.001). Significant differences were observed between R0 and R(un) (p=0.040), as well as between R(un) and R1/2 (p=0.004) (Fig. 2B).
Univariable and multivariable Cox proportional hazards analyses were performed to identify independent prognostic factors for recurrence in patients with pN2 NSCLC. In the univariable analysis, the IASLC R descriptor was a significant prognostic factor, and the HRs of R(un) and R1/2 compared to R0 were 1.42 (95% confidence interval [CI], 1.16 to 1.73; p=0.001) and 1.88 (95% CI, 1.59 to 2.22; p < 0.001), respectively. In the final multivariable Cox model, the IASLC R descriptor was still a significant prognostic factor, but the prognostic difference between R0 and R(un) became insignificant after adjustment for histologic type, subdivided N2 descriptor, subdivided pathological stage, and the performance of adjuvant therapy (HR, 1.18; 95% CI, 0.96 to 1.46; p=0.123) (Table 2).
3. Recurrence dynamics
Fig. 2C shows the hazard rate for recurrence based on the IASLC R descriptor. Patients with R(un) had a 1.4-fold higher hazard rate of recurrence than those with R0. However, the hazard rate curves displayed similar patterns between the three groups, with a peak at approximately 12 months after surgery.
Subgroup analyses for recurrence sites were performed in R0 and R(un) patients, excluding R1/2 patients who underwent incomplete resection. Between R0 and R(un), there were significant differences in overall recurrence (43.1% vs. 54.5%, p=0.002) and distant recurrence (23.8% vs. 39.2%, p < 0.001), while locoregional (8.0% vs. 6.3%, p=0.451) and mixed (11.3% vs. 9.1%, p=0.385) recurrence were similar. Further analysis was performed after stratification with the IASLC N2 descriptor (Table 3). There was no significant difference in recurrence rate between R0 and R(un) in locoregional and mixed recurrence for patients with N2a1, N2a2, and N2b (Table 3). However, a significant difference was shown in distant recurrence for patients with N2a1 (18.0% vs. 32.7%, p=0.037) and in N2a2 (24.7% vs. 37.9%, p=0.044) (Table 3). On the other hand, a significant difference was not shown in distant recurrence for patients with N2b (31.9% vs. 41.7%, p=0.154) (Table 3).
Multivariable Cox analyses for recurrence were performed according to the recurrence site in R0 and R(un) patients. Several factors were adjusted to identify the prognostic significance of positivity for the highest mediastinal LN node station: age, sex, histology, pathologic stage, the IASLC N2 subclassification, and adjuvant treatment. In the overall recurrence, only the N2 subclassification was identified as an independent prognostic factor. Compared to N2a1, the HR of N2a2 was 1.48 (95% CI, 1.14 to 1.92; p=0.004), and that of N2b was 1.79 (95% CI, 1.34 to 2.40; p < 0.001). In the distant recurrence, both positivity for the highest mediastinal LN station and N2 subclassification were confirmed as independent prognostic factors; the HR of positivity for the highest mediastinal LN station was 1.44 (95% CI, 1.08 to 1.92; p=0.013) (Table 4). However, positivity for the highest mediastinal LN station was not a prognostic factor for locoregional and mixed recurrence.
Discussion
In this study, we evaluated the prognostic relevance and recurrence dynamics according to the R descriptor proposed by the IASLC in patients with pN2 NSCLC. Previous studies have validated the prognostic relevance of the R descriptor for OS [5-9], and this study showed similar results in terms of DFS. Notably, despite having no differences in baseline characteristics among the groups, survival outcomes according to the IASLC R descriptor showed significant differences. However, the difference in DFS between R0 and R(un) was ambiguous after adjustment for N2 subclassification.
Since systematic LND was routinely performed in our institution, only 2.4% (275 of 9,245) of patients who underwent anatomical resection underwent an incomplete LND during the study period. Among the included R(un) patients, 88.8% were assigned to R(un) owing to metastasis in the highest mediastinal LN. Therefore, in this study, the clinical implication of positivity for the highest mediastinal LN could be analyzed by comparing the recurrence dynamics of R0 and R(un).
With regards to the recurrence dynamics between R0 and R(un), there were significant differences in overall recurrence (p=0.002) and distant recurrence (p < 0.001), but not in locoregional (p=0.451) and mixed recurrence (p=0.385) (Table 3). In addition, after adjustment of age, sex, histology, pathologic stage, N2 subclassification, and adjuvant treatment, multivariable analysis according to the recurrence site revealed that positivity for the highest mediastinal LN was an independent prognostic factor in distant recurrence, along with N2 subclassification (HR, 1.44; 95% CI, 1.08 to 1.92; p=0.013). However, it was not a prognostic factor in overall, locoregional, and mixed recurrence (Table 4). Therefore, compared to patients classified as R0, a higher recurrence rate in R(un), determined primarily due to positive highest mediastinal LN, was significantly associated with distant metastasis rather than locoregional metastasis.
After stratification with the N2 subclassification, there were still significant differences in distant recurrence between R0 and R(un) for patients with N2a1 and N2a2 (p=0.037 and p=0.044, respectively), which was not for patients with N2b (p=0.154) (Table 3). Considering that metastasis in the highest mediastinal LN implies a potential risk of metastasis to distant sites along the lymphatics, we think this phenomenon might be related to down-staging from N2b to N2a1 or N2a2 in left-sided patients with positive highest mediastinal LN or from N3 to N2a1 or N2a2 in patients with positive highest mediastinal LN. Even in our institution, where LND is routinely performed, the dissection of 2L and 4L for left-sided tumor is not routinely performed. Wang et al. [18] reported that 4L LN metastasis was observed in 20.9% (29/139) of patients, and those in the 4L dissection group had better OS and DFS than the group that did not. Albeit statistically not significant, difference in distant recurrence between R0 and R(un) was also observed in the N2b group (31.9% vs. 41.7%) (Table 3). It means there could be a potential risk of metastasis to the N3 LN (1R or 1L) or distant sites for R(un) patients.
The prognostic difference between the R0 and R(un) groups was significant in univariable analysis (HR, 1.42; p=0.001) but not in multivariable analysis after adjusting for cofactors such as N2 subclassification (HR, 1.18; p=0.123) (Table 2). Several studies have reported the heterogeneous prognosis of N2 disease [16,19], and result of multivariable analysis in this study also showed significant prognostic differences according to the N2 subclassification (Table 2). Therefore, the prognostic difference between the R0 and R(un) could be due to the proportion of N2 subclassification in the R0 and R(un) included in this study (Table 1). However, results of our study showed that positivity for the highest mediastinal LN was significantly associated with distant recurrence even after various cofactors (HR, 1.44; p=0.013). In contrast, there was no significant association with locoregional recurrence (HR, 0.61; p=0.142) or mixed recurrence (HR, 0.95; p=0.842) (Table 4). The difference in DFS between the R0 and R(un) should be considered in the context of encompassing locoregional, distant, and mixed recurrences. Therefore, despite significant differences in distant recurrence between the two groups, this may not result in a significant difference in DFS. In addition, recently proposed 9th edition N classification of the TNM staging system has shown a clear prognostic discrimination ability by subdividing N2 into N2a and N2b. In contrast, N classification proposed by Huang et al. [19] pointed out limitations of overlaps in survival curve between some subcategories. Therefore, in the future study, validation of the prognostic implications of the positivity for the highest mediastinal LN with larger sample size and the newly proposed 9th edition N classification is needed.
Mountain and Dresler [12] and the Spanish Society of Pneumology and Thoracic Surgery (SEPAR) proposed that the absence of metastasis in the highest mediastinal LN station should be fulfilled to achieve complete resection [20]. Given that achieving complete resection is directly related to the improvement of locoregional recurrence, the fact that positivity for the highest mediastinal LN is highly associated with distant metastasis might be somewhat puzzling for clinicians. For example, if metastasis in the 2R or 4L LN is confirmed during surgery, it is unclear in terms of oncological benefit whether to aggressively perform further dissection of the 1R or 2L LNs. Moreover, even if patients have a positive highest mediastinal LN station, there are currently no additional therapeutic options available for these patients. Nevertheless, these findings should not be interpreted as if the IASLC R(un) status is less relevant in patients with positive highest mediastinal LN. The R descriptors proposed by the IASLC aim to encourage the performance of high-quality surgical staging that is critical for designing and analyzing clinical trials of adjuvant therapies. Consequently, R descriptors should be utilized to help clinicians decide the appropriate adjuvant therapy and perform well-designed clinical trials. Furthermore, based on the findings of this study, patients with positivity for the highest mediastinal LN may have a higher incidence of distant recurrence compared to IASLC R0 patients, suggesting the need for more thorough surveillance for these patients.
There has been a considerable discrepancy regarding the prognostic impact of positivity of the highest mediastinal LN in several studies [5,6,21]. Osarogiagbon et al. [6] reported that positivity for the highest mediastinal LN was not an independent prognostic factor for survival. However, it was from a population-based multi–institutional study conducted in the United States, and the detailed information about cancer recurrence in that study might be insufficient compared to our study. In addition, the majority of their patients (2,004/3,203, 62.6%) underwent incomplete LND, whereas our study group comprises a notably lower proportion of patients who underwent incomplete LND (25/286, 8.7%). Consequently, the predominant reason for the composition of R(un) is due to insufficient LN evaluation, while positivity for the highest LN station accounts for merely 5.8% of cases (119 out of 2,044). Insufficient LN assessment can lead to inaccurate N staging, making it difficult to determine whether the highest mediastinal LN has metastasis. Thus, it might act as a bias for identifying the clinical impact of the highest mediastinal LN on prognosis. Park et al. [21] also reported that positivity for the highest mediastinal LN station was not a prognostic factor in patients with pathologic N2. However, unlike this study, they did not perform the analysis after stratification according to the N2 subclassification, and the results might not be consistent with our results depending on the proportion of patients with N2b. On the other hand, Edwards et al. [5] reported that positivity for the highest mediastinal LN was an independent factor for survival, and that it was more prominent in the N-positive patient group.
There are several limitations in our study. First, because this is a retrospective study performed in a single center, our study has inevitable selection bias. Second, because the LN evaluation strategy in our institution was systematic LND, only 2.4% (275/9,245) of patients underwent insufficient LND during the study period. Therefore, we could not analyze patients with pN0 and pN1 due to the small number of patients with R(un) status. Third, there is still no information about pleural lavage and only seven cases of carcinoma in situ on the bronchus, so we could not perform an analysis of these factors. Fourth, staging-related issues may arise if all patients are restaged using the 8th edition of the UICC TNM staging system. Because this study included data spanning wide durations and encompassing UICC 6th, 7th, and 8th TNM staging systems, adjuvant treatment or clinical decisions might be different according to the staging system working at that time. This might have influenced the recurrence patterns. Finally, we indirectly evaluated the prognostic impact of the highest mediastinal LN positivity through recurrence dynamics. However, whether there is remnant tumor burden in the distal lymphatics or beyond the dissection field in patients positive for the highest mediastinal LN station cannot be confirmed.
In conclusion, the IASLC R criteria has prognostic relevance in patients with pathologic N2 NSCLC. The R(un) is a highly heterogeneous group, and it is possible that the involvement of the highest mediastinal LN involvement affects distant recurrence patterns.
Electronic Supplementary Material
Supplementary materials are available at Cancer Research and Treatment website (https://www.e-crt.org).
Notes
Ethical Statement
This study was approved by the Institutional Review Board of Asan Medical Center, and this center waived the need for informed consent (approval no. 2023-1514, approval date: 28 November 2023).
Author Contributions
Conceived and designed the analysis: Kim IH, Yun JK.
Collected the data: Kim IH, Lee GD.
Contributed data or analysis tools: Kim IH, Lee GD, Choi S, Yun JK.
Performed the analysis: Kim IH, Kim YH, Yun JK.
Wrote the paper: Kim IH, Kim HR, Yun JK.
Writing-review & editing: Kim IH, Lee GD, Choi S, Kim HR, Kim YH, Kim DK, Park SI, Yun JK.
Conflict of Interest
Conflict of interest relevant to this article was not reported.
Acknowledgements
We would like to thank Editage (www.editage.co.kr) for English language editing.