Skip to main content

ORIGINAL RESEARCH article

Front. Oncol., 14 November 2019
Sec. Radiation Oncology
This article is part of the Research Topic Magnetic Resonance Imaging for Radiation Therapy View all 17 articles

Predictive Values of MRI and PET Derived Quantitative Parameters for Patterns of Failure in Both p16+ and p16– High Risk Head and Neck Cancer

\nYue Cao,,
Yue Cao1,2,3*Madhava AryalMadhava Aryal1Pin LiPin Li4Choonik LeeChoonik Lee1Matthew Schipper,Matthew Schipper1,4Peter G. HawkinsPeter G. Hawkins1Christina Chapman,Christina Chapman1,5Dawn OwenDawn Owen1Aleksandar F. DragovicAleksandar F. Dragovic1Paul SwiecickiPaul Swiecicki6Keith CasperKeith Casper7Francis WordenFrancis Worden6Theodore S. LawrenceTheodore S. Lawrence1Avraham EisbruchAvraham Eisbruch1Michelle MierzwaMichelle Mierzwa1
  • 1Department of Radiation Oncology, University of Michigan, Ann Arbor, MI, United States
  • 2Department of Radiology, University of Michigan, Ann Arbor, MI, United States
  • 3Department of Biomedical Engineering, University of Michigan, Ann Arbor, MI, United States
  • 4Department of Biostatistics, University of Michigan, Ann Arbor, MI, United States
  • 5Department of Radiation Oncology, VA Ann Arbor Healthcare System, Ann Arbor, MI, United States
  • 6Department of Internal Medicine, University of Michigan, Ann Arbor, MI, United States
  • 7Department of Otolaryngology, University of Michigan, Ann Arbor, MI, United States

Purpose: FDG-PET adds to clinical factors, such tumor stage and p16 status, in predicting local (LF), regional (RF), and distant failure (DF) in poor prognosis locally advanced head and neck cancer (HNC) treated with chemoradiation. We hypothesized that MRI-based quantitative imaging (QI) metrics could add to clinical predictors of treatment failure more significantly than FDG-PET metrics.

Materials and methods: Fifty four patients with poor prognosis HNCs who were enrolled in an IRB approved prospective adaptive chemoradiotherapy trial were analyzed. MRI-derived gross tumor volume (GTV), blood volume (BV), and apparent diffusion coefficient (ADC) pre-treatment and mid-treatment (fraction 10), as well as pre-treatment FDG PET metrics, were analyzed in primary and individual nodal tumors. Cox proportional hazards models for prediction of LRF and DF free survival were used to test the additional value of QI metrics over dominant clinical predictors.

Results: The mean ADC pre-RT and its change rate mid-treatment were significantly higher and lower in p16– than p16+ primary tumors, respectively. A Cox model identified that high mean ADC pre-RT had a high hazard for LF and RF in p16– but not p16+ tumors (p = 0.015). Most interesting, persisting subvolumes of low BV (TVbv) in primary and nodal tumors mid-treatment had high-risk for DF (p < 0.05). Also, total nodal GTV mid-treatment, mean/max SUV of FDG in all nodal tumors, and total nodal TLG were predictive for DF (p < 0.05). When including clinical stage (T4/N3) and total nodal GTV in the model, all nodal PET parameters had a p-value of >0.3, and only TVbv of primary tumors had a p-value of 0.06.

Conclusion: MRI-defined biomarkers, especially persisting subvolumes of low BV, add predictive value to clinical variables and compare favorably with FDG-PET imaging markers. MRI could be well-integrated into the radiation therapy workflow for treatment planning, response assessment, and adaptive therapy.

Introduction

Locoregional failure (LRF) remains a clinical challenge for poor prognosis locally advanced squamous cell carcinoma of the head and neck (HNSCC) treated with definitive chemoradiation therapy (CRT) (1). It is important to identify imaging markers of LRF that identify patients and tumor subvolumes that may benefit from intensified locoregional therapy in the form of radiation boost, targeted systemic therapy, or surgical intervention.

We and others have been developing prognostic and predictive imaging markers of PET and MRI for LRF, distant metastases, progression free survival (PFS), and overall survival (OS) (2-20). Retrospective studies of pre-treatment FDG-PET that quantify cellular glucose metabolism have identified metabolic tumor volume (MTV), total lesion glycolysis (TLG), and mean/max standard uptake value (SUV) in MTV as prognostic for LRF, PFS, and OS in HNSCC (26). Furthermore, FDG-PET has been incorporated into standard of care work-up and follow-up for HNSCC (7, 8). Functional MRI incorporating diffusion and perfusion parameters is an emerging advanced imaging modality in HNSCC. In particular, apparent diffusion coefficient (ADC) correlates with locoregional and distant progression (911). Poorly perfused and low oxygenation tumors have been shown to be associated with LRF and worse survival outcomes (1218).

Despite this progress, it has been difficult to determine which imaging biomarkers should be used to individualize treatment for the patients with locally advanced HNSCC. Most head and neck cancer imaging studies to date include heterogeneous populations of various disease sites, stages, and prognosis. Few imaging studies investigate how p16 status affects imaging parameters pre- and mid-treatment. With respect to ADC in particular, no study to date has evaluated ADC changes during RT for p16+ vs. p16– tumors. A single study investigated ADC differences between HPV+ and HPV– HNSCC, including only 6 HPV+ patients (8%), and found that pre-treatment ADC in HPV+ HNSCC patients was significantly lower than in HPV– patients (19). Furthermore, at the tumor and subtumor level, there is no report on imaging biomarker differences between tumors with local, regional, or distant failure as site of first failure compared to disease free patients. This is an important issue, as it would help stratify the patients for local or systemic intensified or de-intensified therapy. Finally, poorly perfused tumor subvolumes are largely spatially distinct from areas of high FDG uptake and high restricted water diffusion in the same patients, and the spatial correlation between high glucose metabolism and high restricted water diffusion varies greatly from patient to patient (20, 21). These studies question whether both FDG PET and MRI biomarkers are necessary to guide adaptive RT in HNSCC.

This study aimed to (1) investigate p16+ effects on imaging parameters and their early response rates; (2) assess differences between imaging biomarkers of tumors with local, regional or distant progression and those with no evidence disease (NED), and (3) compare the predictive values of MRI and PET biomarkers. We hypothesized that p16+ status could affect imaging biomarkers and their early response rates, and MRI-based QI metrics could add to clinical predictors of treatment failure more significantly than FDG-PET metrics for local, regional and distant failure.

Methods

Patients

Imaging analysis was performed on 54 patients [median age of 61 years; 7 females; 31 p16+ (57%)] with advanced HNSCC who were enrolled in a randomized phase II clinical trial between March 2014 and January 2018 (Table 1). The trial was approved by the Institutional Review Board of the University of Michigan, including a parallel imaging study to investigate the predictive values of QI metrics for tumor progression. Written consent was obtained from all enrolled patients. Eligibility included patients with p16+ T4/N3 squamous cell carcinoma of oropharynx or locally advanced p16– HNSCC if planned to undergo definitive CRT. All patients were evaluated for p16 status by immunohistochemistry. After completion of CRT, patients were followed up every 2–3 months per standard care for oncologic outcomes as well as toxicity. Tumor recurrences were scored as LF, RF, or DF, or a combination thereof.

TABLE 1
www.frontiersin.org

Table 1. Patient characteristics.

MRI and PET Acquisition

Patients underwent FDG-PET/CT scans pre-RT within 4 weeks of RT as a part of standard care. Clinical FDG-PET/CT scans were performed on various PET scanners by following the standard clinical protocol (22).

MRI scans were acquired pre-RT (within 2 weeks) and at fraction 10 (20 Gy) as a part of the protocol. All MRI scans were acquired on a 3T scanner (Skyra, Siemens Healthineers), including anatomic, diffusion weighted (DW), and DCE T1-weighted imaging series. All patients were scanned in the treatment position using an individual-patient immobilization 5-point mask and bite block or aquaplast mold as required for treatment. DW images were acquired with spatial resolution of ~1.2 × 1.2 × 4.8 mm and b-values of 50 and 800 s/mm2 by either a 2D spin-echo single shot echo-planar pulse sequence or a readout segmentation of long variable echo-trains (RESOLVE) pulse sequence that reduced geometric distortion (23). Sixty T1-weighted DCE image volumes were acquired using a 3D gradient echo pulse sequence in a sagittal orientation with voxel size ~1.5 × 1.5 × 2.5 mm during an injection of one standard dose of Gd-DTPA. Post-Gd T1-weighted images were acquired in the axial plane with spatial resolution of 0.875 × 0.875 × 3.3 mm by a 2D fast spin echo sequence with fat saturation.

Image Analysis and Registration

Blood volume (BV) maps were quantified from DCE-MRI using the modified Tofts model implemented in an in-house imFIAT Analysis Tool, which was validated using a digital reference object (24). ADC maps were calculated from DW images with b-values of 50 and 800 to mitigate the perfusion effect by using in-house software that was technically validated in a QIN collaborative project (25). Since using the individual-patient immobilization devices reduced gross movement of head and neck during scanning dramatically, BV and ADC maps were reformatted to match voxel-by-voxel of post-Gd T1-weighted images acquired in the same session using coordinates in DICOM headers. SUV of FDG-PET was calculated. Pre-RT FDG-PET/CT and mid-treatment MR images were co-registered to pre-RT post-Gd T1-weighted images using rigid body transformation and mutual information. Target displacement errors, including image mis-registration and geometric distortion in ADC maps, between image series were assessed and reported previously (20). Reproducibility of BV maps was 16%, which was reported previously (26).

Tumor Volumes and Subvolumes

Gross tumor volume (GTV) of primary and nodal disease was contoured individually on post-Gd T1-weighted images by treating attending head and neck radiation oncologists and reviewed by the trial PI (MM). For this cohort of patients with locally advanced HNSCC, gross cystic or necrotic regions and tumor invasion into blood vessels occurred in many tumors, and therefore were excluded from the GTVs for following analyses of quantitative image (QI) metrics by applying simple thresholds. For the ADC analysis, a threshold of >2.7 × 10−3 mm2/s (10% below free water diffusion) was used to exclude gross necrosis and blood vessels, and a threshold of <0.0001 × 10−3 mm2/s was used to exclude air. Then, a low BV subvolume of the GTV (TVBV) was created using a threshold of BV <7.64 ml/100 g reported previously based upon a histogram analysis (16). The low ADC subvolume of the GTV (TVADC) was defined as ADC <1.2 × 10−3 mm2/s based on an ADC-histogram analysis (20), which is also consistent with the mean ADC reported by others (21). A MTV was defined as FDG SUV >50% of a value averaged over 4 voxels with maximum SUVs (MTV50).

Quantitative Imaging Metrics

QI metrics in tumor volumes and their mid-treatment changes were analyzed for prediction of LF, RF, and DF. Tumor volume metrics included GTV, TVBV, TVADC, MTV50. Mean values of ADC and BV in GTV excluding blood vessels and necrosis, mean and max SUVs in MTV50, and TLG of MTV50 were calculated for each primary or nodal tumor as well as for all tumors in each patient.

Treatment

The patients were randomized to a standard arm of RT (70 Gy in 35 fractions) or an experimental arm. In the experimental arm, a union of the persisting TVBV pre-RT to 2 weeks and persisting TVADC pre-RT to 2 weeks received 2.5Gy per fraction for the last 15 of 35 fractions. If the union of persisting subvolumes pre-RT to 2 weeks was <1 cc, the patient was entered into an observation arm and treated by the standard RT (70 Gy in 35 fractions). Patients were planned to receive weekly cisplatin 40 mg/m2, and patients considered to be cisplatin ineligible were treated with weekly carboplatin AUC2.

Statistical Analysis

First, we assessed the p16 effect on imaging parameters and parameter change rates at 2 weeks compared to pre-RT using the Mann-Whitney U-test. Secondly, we assessed whether MRI and PET biomarkers had similar predictive values for LRF and DF free survival. For the analysis of LRF, most previous analyses considered either LF, RF, or LRF as an event, of which the model was useful for stratification of the patients but not for stratification of the tumors for intensified adaptive RT. Tumor progression could occur in one or a few treated tumors (primary or nodal tumor) or in none. Therefore, we applied Cox proportional hazards models to individual (primary or nodal) tumors for prediction of failure. The individual tumor failure free rate (ITFFR) was defined from the start of RT to the date of progression of the tested (primary or nodal) tumor. ITFFR times were censored for all tumors from a patient at the earlier of DF, death or last follow-up. Whether primary and nodal tumors can be analyzed together was tested for each imaging parameter. To compare the predictive values of MRI and FDG PET biomarkers, imaging metrics were assessed one at a time in models also including p16 as a co-variable, which is the most important clinical variable for LRF (2729). Distant failure free survival (DFRS) was defined as the time interval from the start of RT to the date of DF. The Cox models were fitted including a single QI metric and clinical stage T4/N3 vs. other (non-T4/N3) as the sole clinical variable (3032), and entering one imaging parameter at a time. Each QI metric was summed up or averaged over all nodal tumors for volume-related or intensity-related metrics, respectively. In the DFFS model, patients were censored at the first occurrence of any local or regional failure, death or last follow-up. If there were any significant differences of imaging parameters between p16– and p16+ tumors, we considered an interaction term in the Cox model or an analysis in different Cox models as appropriate. Since multiple comparisons were made, p-values were corrected using false discovery rate (FDR) control, and corrected p < 0.10 were considered significant. Finally, we assessed if there were any significant differences of imaging biomarkers between the tumors that never progressed, those that demonstrated local or regional progression, and those that were locoregionally controlled but metastasized distantly. This landmark analysis used outcomes at 18 months as a cutoff. The tumors were excluded from the analysis if the tumor had local or regional progression after 18 months or the tumor had no progression but the follow-up was shorter than 18 months. As the data were not Gaussian distributed, non-parametric tests were used: Kruskal-Wallis test for the three-group comparison and Wilcoxon rank test for the comparison between local or regional failure and NED. The p-values were corrected with FDR control, and <0.1 were considered as significant. Since 37% of the patients received higher doses, we tested the dose effect before performing the proposed analyses.

Results

Treatment Failure

This cohort of 54 patients with locally advanced HNSCC had large primary GTVs with a median value of 60.5 cc (range: 10.2–595.2 cc; SD: 86.8 cc; Table 1), which was several times greater than most reported studies (26, 911, 33). Eleven patients (20%) (3 p16+) have had local recurrence. Nine patients (17%) (2 p16+) have had regional recurrence, including one patient (p16–) who failed regionally at two separate treated lymph node locations, and 2 (1 p16– and 1 p16+) who had RF at the locations of non-enlarged/non-FDG avid nodes before RT. Fourteen patients (7 p16+) had distant failure with or without local and regional failure. All cases with LF or RF alone were confirmed pathologically, and distant metastases were diagnosed pathologically or by overt radiographic presentation. Twelve patients have died of HNC (3 p16+), and one patient died cancer-free of other causes. For the patients who did not have progression at the time of analysis, median follow-up was 24 months (range: 10–58 months).

Effects of p16 on Imaging Parameters and Change Rates

We found that both baseline ADC and ADC change after radiation were significantly different between p16+ and p16– primary tumors. The p16– primary tumors had significantly greater mean ADCs pre-RT [1.48 ± 0.05(SEM) μm2/ms], and significantly smaller rates of increase after 10 fractions of RT (10.0% ± 1.2%) than p16+ primary tumors (1.34 ± 0.04 μm2/ms and 21.2% ± 3.1%, p = 0.04, and p = 0.009, respectively). However, there was no significant difference in mean ADC between p16– and p16+ nodal tumors pre-RT or at 2 weeks as well as ADC increased rates (p > 0.7), see Figure 1. Pre-RT GTVs of p16– primary tumors (75 ± 12.1 cc) as well as change rates at 2 weeks (−16.2% ± 3.9%) were similar to p16+ ones (79.2 ± 18.9 cc, and −16.7% ± 3.3%, respectively). Mean GTVs as well as change rates at 2 weeks for p16– and p16+ nodal tumors were not significantly different (p > 0.5), 24.1 ± 8.1 cc and 21.1 ± 5.9 cc of GTVs and −22.4% + 6.7% and −16.5% + 4.5% of change rates for respective p16– and p16+ nodal tumors Also, there was no significant difference in other imaging parameters between p16+ and p16– primary or nodal tumors (p > 0.1). Examples of images are shown in Figure 2.

FIGURE 1
www.frontiersin.org

Figure 1. (A) Mean ADC in p16– and p16+ primary tumors pre-RT and after 10 fractions (2 weeks) of radiation therapy. (B) Mean ADC in p16– and p16+ nodal tumors pre-RT and after 10 fractions (2 weeks) of radiation therapy. (C) Mean ADC change rates in p16– and p16+ primary and nodal tumors after 10 fractions of radiation therapy compared to pre-treatment. *p < 0.05, **p < 0.01.

FIGURE 2
www.frontiersin.org

Figure 2. Post Gd T1 weighted images (left), ADC (second left), blood volume (second right), and SUV (right) FDG PET pre-RT (top) and after 10 fractions of radiation therapy (bottom). GTV: magenta; low BV subvolume: yellow; low ADC subvolume: cyan; MTV: red. Note persistent low BV and low ADC subvolumes after 10 fractions of radiation therapy (yellow arrows).

Predictive Values of MRI and PET Imaging Parameters for Local and Regional Progression

First, we did not detect significant difference in local and regional control rates between two-dose arms yet so that the patients who received different doses were analyzed together. For prediction of local progression, mean ADC pre-RT of primary tumors was the only parameter found significant in a univariate Cox model. Since there was no significant difference in mean ADC between primary and nodal tumors, we combined primary and nodal tumors in a single model (53 primary tumors and 82 nodal tumors). For prediction of ITFFR, considering the p16 effect on mean ADC of primary tumors, the Cox model included p16 status, pre-RT mean ADC, and the interaction of pre-RT mean ADC and p16 status. We found that p16 had a significant effect on tumor control (HR p16+ vs. p16– of 0.21, p = 0.005), pre-RT mean ADC had a significant effect in p16– tumors (HR per 1 SD increase in ADC = 1.9, p = 0.015) but no effect in p16+ tumors (HR = 1.0, p = 1.0). The interaction between p16 status and ADC was not statistically significant (p = 0.24, Table 2).

TABLE 2
www.frontiersin.org

Table 2. Cox model of mean ADC effects.

Since QI metrics other than mean ADC were significantly different between primary and nodal tumors (p < 0.05), the QI metrics of nodal tumors were tested separately for prediction of regional failure free rates. In Cox models of 82 nodal tumors with p16 status as a co-variate, GTV pre-RT and at 2 weeks, TVBV at 2 weeks, mean and max SUV in MTV50 pre-RT, MTV50 pre-RT, TLG pre-RT, and change in GTV at 2 weeks vs. pre-RT were significant with p < 0.07 with FDR control, see Table 3. It is interesting to note that GTV pre-RT and at 2 weeks as well as mean SUV and TLG pre-RT have the highest c-index (> 0.9). However, MTV50 and TLG as well as TVBV were strongly correlated with GTV pre-RT (range of r between 0.88 and 0.90), suggesting that these metrics are not independent of GTV. The mean and max SUV in MTV50 were strongly correlated each other (r = 0.98) but modestly correlated with GTV pre-RT (range of r between 0.65 and 0.67).

TABLE 3
www.frontiersin.org

Table 3. Cox models for RFFS.

Predictive Values of Imaging Biomarkers for Distant Progression

For prediction of distant progression, Cox models identified that TVBV of primary tumors at 2 weeks, total TVBV of all nodal tumors pre-RT and at 2 weeks, total GTV of all nodal tumors at 2 weeks, mean and max SUV of all nodal MTV50 pre-RT, and TLG pre-RT of all nodal tumors had a nominal p < 0.05 without FDR. With FDR control, total GTV of all nodal tumors at 2 weeks, mean and max SUV of all nodal MTV50 pre-RT, and TLG pre-RT had a p < 0.1, see Table 4. We tested whether the significant predictors could provide any complimentary information to clinical stage of T4/N3 and the sum of all nodal GTVs at 2 weeks for prediction of DF, and found that neither total TVBV, nor mean and max SUV, nor total TLG of all nodal tumors had a p < 0.3, and only TVBV of primary tumors at 2 weeks showed marginally significant (p = 0.06).

TABLE 4
www.frontiersin.org

Table 4. Cox models for DFFS.

Imaging Biomarkers for Differentiation of Tumors With LF (or RF), DF, and NED

For primary tumors, the subvolumes of low BV pre-RT showed a descending trend from LF, DF, and NED with a marginally significant p-value of < 0.06 without FDR control, see Table 5. Figure 3 shows the subvolumes of low BV of primary tumors with LF, DF, and NED pre-RT and at 2 weeks as well as its change rates after 10 factions of RT. Post ad hoc analysis showed that the change rates of low BV subvolume were significant smaller in primary tumors with DF (−0.05% ± 0.16%) than tumors with LF (−0.49 ± 0.08%) and tumors with NED (−0.45 ± 0.09%) with p values of < 0.03 and < 0.015, respectively.

TABLE 5
www.frontiersin.org

Table 5. Differences of imaging biomarkers among tumors with LF (or RF), DF, and NED.

FIGURE 3
www.frontiersin.org

Figure 3. Box and Whisker plots of the subvolumes of low BV in primary tumors with local failure (blue), distant failure (orange), and no evidence disease (gray) pre-RT (A) and after 10 fractions of radiation therapy (B). (C) Plots the change rates of the low BV subvolumes after 10 fractions of radiation compressed to pre-RT.

For nodal tumors, GTV pre-RT and at 2 weeks, the subvolume of low BV pre-RT, mean ADC at 2 weeks, mean BV at 2 weeks, and mean/max SUV of MTV50 pre-RT were different among DF, RF, and NED groups with p < 0.05 without FDR control and p ≤ 0.1 with FDR control, see Table 5. Again, GTV of nodal tumors was a strongest parameter to differentiate the three groups with different outcomes. Regarding the difference between DF and RF groups, only mean BV values at 2 weeks had a p < 0.05 without FDR control but p > 0.1 with FDR control. Figure 4 shows GTVs, the subvolumes of low BV, mean ADC, and mean BV of nodal tumors with RF, DF, and NED pre-RT and at 2 weeks. Figure 5 shows mean SUV, max SUV, and TLG of nodal tumors with RF, DF, and NED pre-RT.

FIGURE 4
www.frontiersin.org

Figure 4. Box and Whisker plots of the GTVs (A,B), the subvolumes of low BV (C,D), mean ADC (E,F), and mean BV (G,H) of nodal tumors with RF (blue), DF (orange), and NED (gray) pre-RT and at 2 weeks.

FIGURE 5
www.frontiersin.org

Figure 5. Box and Whisker plots of the mean SUV (A), max SUV (B), and TLG (C) of nodal tumors with RF (blue), DF (orange), and NED (gray) pre-RT.

Discussion

In this study, we investigated p16 effects on MRI and PET QI metrics, imaging biomarker differences as a function of tumor control (local, regional, or distant), and the predictive values between MRI and PET biomarkers for tumor progression in locally advanced poor prognosis HN cancers. Our cohort of patients had large tumor volumes compared to previously reported literature (26, 911, 33). We found the p16– primary tumors had elevated ADC values pre-RT and low early response rates compared to p16+ tumors; the latter of which has not been previously reported. Also, high mean ADC value pre-RT is a hazard for local and regional failure of p16– tumors. Multiple MRI and PET imaging parameters (including GTV, ADC, BV, SUV, and TLG) predicted RF and DF, but the nodal GTV defined on anatomic MRI was the strongest predictor. Most interesting, we report for the first time that the persistent low BV in primary and nodal tumors during the early course of CRT is associated with high-risk for distant failure. In order to identify patients who may benefit from intensified local therapy in the form of a radiation boost or surgical intervention, or from intensified systemic therapy (30, 34), we analyzed the significant imaging predictors found in Cox modeling for differentiation of the tumors that were controlled compared to those with LF, RF, or DF. The performance of MRI related parameters is stronger than PET parameters. Although PET is a part of standard care, MRI could play an important role from treatment planning, to early response assessment, and boost target definition.

We found a p16 effect on ADC and ADC change rates during the early course of RT. The p16– primary tumors had significantly greater mean ADC values pre-RT and smaller increases in ADC after 2 weeks of CRT than p16+ primary tumors. Furthermore, the p16– tumors from patients with local or regional failure had significantly greater mean ADC values pre-RT and mid-treatment than those from disease free patients. These results are consistent with previous reports that the pre-RT high ADC is negatively prognostic for HN cancers (911). A recent study shows that ADC is significantly and inversely correlated with cell density but also significantly and positively correlated with the percentage area of stroma in laryngeal and hypopharyngeal carcinoma (35). The former finding has been reported previously in animal studies, prostate cancer and lymphomas (3639), and is related to restricted water diffusion due to high cellularity. The latter finding suggests that a large percentage area of stroma in HN cancers is associated with a high ADC. Stroma has been shown to be negatively prognostic in several cancers, to promote tumor growth and invasion, and to potentially protect tumors from delivery of chemotherapy (4045). ADC behaviors in the p16– tumors could be explained by their increased stroma. HPV-related oropharynx cancers are histologically basaloid in histology with significant tumor lymphocytic infiltration, which is associated with improved prognosis (46, 47) and decreased ADC. ADC, although a promising QI metric for differentiation of local and regional failure, and even distant failure, is affected by multiple biologic and physiologic factors, including cell density and stroma as well as cyst and necrosis (In this study, we excluded grossly cystic and necrotic regions for QI metric analysis).

The low BV in primary tumors and persisting during the early course of RT have reported previously to be associated with LF (1217). However, there is no report that the low BV and its low response rate in HNSCC during the early course of RT is associated with DF. The subvolumes of low BV in primary tumors show a descending trend from LF, to DF and NED. The response rate of low BV could be used to differentiate the tumor at high-risk for LF or DF from NED, and thereby adapting intensified local or systematic therapy for the patients with different progression risks.

Pretreatment FDG QI metrics, including MTV, TLG and mean/max, have been reported to be correlated with PFS and OS in the patients with HN cancers treated with CRT (36). We found that the high mean/max SUV and large TLG in nodal tumors were risk factors for nodal failure, and that the sum of TLG over all nodal MTVs was a negative prognostic factor for DFFS, which is consistent with several previous reports (24). Although TLG accounts for both the size and SUV of MTV, we found that nodal TLG was strongly correlated with MRI-defined GTV, and nodal GTV was the strongest predictor for RF in our study. For prediction of RF and DF, several other MRI parameters (including GTV, ADC, and BV) perform as well as FDG PET related parameters. When including T4/N3 and total nodal GTV in the Cox model, no other imaging parameters including PET were found to be significant. Finally, there were no FDG PET related parameters that could predict LF.

Radiomics analysis of CT and PET features is another area of imaging analysis that could provide complimentary information to the present study. Radiomics analysis that extracts the large amounts of quantitative textural features from CT, PET, and MRI has been investigated for the prediction of local control, PFS, and OS in head and neck cancers (4852). Through the feature selection and reduction processes, a small number of features have been found to have prognostic or predictive value. These features include general categories of statistical energy, shape compactness, gray level non-homogeneity, and gray level non-uniformity. These features may represent different tumor phenotypes. However, it is hard to link the feature to tumor physiology, pathology and biology. Furthermore, radiomics approaches require a large amount of high quality image data, and high-throughput.

A limitation of the present analysis includes RT boost of tumor subvolumes with persistent low BV and low ADC on our clinical trial. This could affect QI metrics that are identified for prediction of treatment failure. We will perform this analysis on patients who are on the standard treatment arm when the trial is completed and the data have matured. Nevertheless, we found that persistent low BV in primary and nodal tumors carries a high-risk for nodal and distant failure, the low response rate of low BV has a high-risk for distant failure, and the low response rate of ADC is for p16– primary tumors. MRI derived biomarkers perform at least as well as FDG PET defined ones. As MRI based planning is already well-integrated into radiation therapy, our findings suggest that MRI based response assessment will be a valuable guide in adaptive radiation therapy.

Data Availability Statement

The image data that were collected in the patients with at least 10 months follow-up in the study are included in the manuscript/supplementary files.

Ethics Statement

The studies involving human participants were reviewed and approved by the institute review board of University Michigan. The patients/participants provided their written informed consent to participate in this study.

Author Contributions

All authors contributed significantly for study design, patient enrollment, image acquisition and analysis, statistic analysis, data interpretation, or writing.

Funding

This work was supported by NIH/NCI grants U01CA183848 and RO1CA184153.

Conflict of Interest

YC is co-owner of a US patent of No. 61/656,323, which is entitled “The subvolume identification for prediction of treatment outcome”.

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.

References

1. Jemal A, Bray F, Center MM, Ferlay J, Ward E, Forman D. Global cancer statistics. CA Cancer J Clin. (2011) 61:69–90. doi: 10.3322/caac.20107

PubMed Abstract | CrossRef Full Text | Google Scholar

2. Chotchutipan T, Rosen BS, Hawkins PG, Lee JY, Saripalli AL, Thakkar D, et al. Volumetric (18) F-FDG-PET parameters as predictors of locoregional failure in low-risk HPV-related oropharyngeal cancer after definitive chemoradiation therapy. Head Neck. (2019) 41:366–73. doi: 10.1002/hed.25505

PubMed Abstract | CrossRef Full Text | Google Scholar

3. Lin YC, Chen SW, Hsieh TC, Yen KY, Yang SN, Wang YC, et al. Risk stratification of metastatic neck nodes by CT and PET in patients with head and neck cancer receiving definitive radiotherapy. J Nucl Med. (2015) 56:183–9. doi: 10.2967/jnumed.114.148023

PubMed Abstract | CrossRef Full Text | Google Scholar

4. Pollom EL, Song J, Durkee BY, Aggarwal S, Bui T, von Eyben R, et al. Prognostic value of midtreatment FDG-PET in oropharyngeal cancer. Head Neck. (2016) 38:1472–8. doi: 10.1002/hed.24454

PubMed Abstract | CrossRef Full Text | Google Scholar

5. Schwartz DL, Harris J, Yao M, Rosenthal DI, Opanowski A, Levering A, et al. Metabolic tumor volume as a prognostic imaging-based biomarker for head-and-neck cancer: pilot results from Radiation Therapy Oncology Group protocol 0522. Int J Radiat Oncol Biol Phys. (2015) 91:721–9. doi: 10.1016/j.ijrobp.2014.12.023

PubMed Abstract | CrossRef Full Text | Google Scholar

6. Tang C, Murphy JD, Khong B, La TH, Kong C, Fischbein NJ, et al. Validation that metabolic tumor volume predicts outcome in head-and-neck cancer. Int J Radiat Oncol Biol Phys. (2012) 83:1514–20. doi: 10.1016/j.ijrobp.2011.10.023

PubMed Abstract | CrossRef Full Text | Google Scholar

7. Gupta T, Master Z, Kannan S, Agarwal JP, Ghsoh-Laskar S, Rangarajan V, et al. Diagnostic performance of post-treatment FDG PET or FDG PET/CT imaging in head and neck cancer: a systematic review and meta-analysis. Eur J Nucl Med Mol Imaging. (2011) 38:2083–95. doi: 10.1007/s00259-011-1893-y

PubMed Abstract | CrossRef Full Text | Google Scholar

8. Isles MG, McConkey C, Mehanna HM. A systematic review and meta-analysis of the role of positron emission tomography in the follow up of head and neck squamous cell carcinoma following radiotherapy or chemoradiotherapy. Clin Otolaryngol. (2008) 33:210–22. doi: 10.1111/j.1749-4486.2008.01688.x

CrossRef Full Text | Google Scholar

9. Hatakenaka M, Nakamura K, Yabuuchi H, Shioyama Y, Matsuo Y, Ohnishi K, et al. Pretreatment apparent diffusion coefficient of the primary lesion correlates with local failure in head-and-neck cancer treated with chemoradiotherapy or radiotherapy. Int J Radiat Oncol Biol Phys. (2011) 81:339–45. doi: 10.1016/j.ijrobp.2010.05.051

CrossRef Full Text | Google Scholar

10. Kim S, Loevner L, Quon H, Sherman E, Weinstein G, Kilger A, et al. Diffusion-weighted magnetic resonance imaging for predicting and detecting early response to chemoradiation therapy of squamous cell carcinomas of the head and neck. Clin Cancer Res. (2009) 15:986–94. doi: 10.1158/1078-0432.CCR-08-1287

PubMed Abstract | CrossRef Full Text | Google Scholar

11. Lambrecht M, Van Calster B, Vandecaveye V, De Keyzer F, Roebben I, Hermans R, et al. Integrating pretreatment diffusion weighted MRI into a multivariable prognostic model for head and neck squamous cell carcinoma. Radiother Oncol. (2014) 110:429–34. doi: 10.1016/j.radonc.2014.01.004

PubMed Abstract | CrossRef Full Text | Google Scholar

12. Bernstein JM, Kershaw LE, Withey SB, Lowe NM, Homer JJ, Slevin NJ, et al. Tumor plasma flow determined by dynamic contrast-enhanced MRI predicts response to induction chemotherapy in head and neck cancer. Oral Oncol. (2015) 51:508–13. doi: 10.1016/j.oraloncology.2015.01.013

PubMed Abstract | CrossRef Full Text | Google Scholar

13. Cao Y, Popovtzer A, Li D, Chepeha DB, Moyer JS, Prince ME, et al. Early prediction of outcome in advanced head-and-neck cancer based on tumor blood volume alterations during therapy: a prospective study. Int J Radiat Oncol Biol Phys. (2008) 72:1287–90. doi: 10.1016/j.ijrobp.2008.08.024

PubMed Abstract | CrossRef Full Text | Google Scholar

14. Chawla S, Kim S, Loevner LA, Hwang WT, Weinstein G, Chalian A, et al. Prediction of disease-free survival in patients with squamous cell carcinomas of the head and neck using dynamic contrast-enhanced MR imaging. AJNR Am J Neuroradiol. (2011) 32:778–84. doi: 10.3174/ajnr.A2376

PubMed Abstract | CrossRef Full Text | Google Scholar

15. Hermans R, Meijerink M, Van den Bogaert W, Rijnders A, Weltens C, Lambin P. Tumor perfusion rate determined noninvasively by dynamic computed tomography predicts outcome in head-and-neck cancer after radiotherapy. Int J Radiat Oncol Biol Phys. (2003) 57:1351–6. doi: 10.1016/S0360-3016(03)00764-8

PubMed Abstract | CrossRef Full Text | Google Scholar

16. Wang P, Popovtzer A, Eisbruch A, Cao Y. An approach to identify, from DCE MRI, significant subvolumes of tumors related to outcomes in advanced head-and-neck cancer. Med Phys. (2012) 39:5277–85. doi: 10.1118/1.4737022

PubMed Abstract | CrossRef Full Text | Google Scholar

17. Brizel DM, Dodge RK, Clough RW, Dewhirst MW. Oxygenation of head and neck cancer: changes during radiotherapy and impact on treatment outcome. Radiother Oncol. (1999) 53:113–7. doi: 10.1016/S0167-8140(99)00102-4

PubMed Abstract | CrossRef Full Text | Google Scholar

18. Stadler P, Becker A, Feldmann HJ, Hansgen G, Dunst J, Wurschmidt F, et al. Influence of the hypoxic subvolume on the survival of patients with head and neck cancer. Int J Radiat Oncol Biol Phys. (1999) 44:749–54. doi: 10.1016/S0360-3016(99)00115-7

PubMed Abstract | CrossRef Full Text | Google Scholar

19. Driessen JP, van Bemmel AJ, van Kempen PM, Janssen LM, Terhaard CH, Pameijer FA, et al. Correlation of human papillomavirus status with apparent diffusion coefficient of diffusion-weighted MRI in head and neck squamous cell carcinomas. Head Neck. (2016) 38 (Suppl 1):E613–8. doi: 10.1002/hed.24051

PubMed Abstract | CrossRef Full Text | Google Scholar

20. Teng F, Aryal M, Lee J, Lee C, Shen X, Hawkins PG, et al. Adaptive boost target definition in high-risk head and neck cancer based on multi-imaging risk biomarkers. Int J Radiat Oncol Biol Phys. (2018) 102:969–77. doi: 10.1016/j.ijrobp.2017.12.269

PubMed Abstract | CrossRef Full Text | Google Scholar

21. Houweling AC, Wolf AL, Vogel WV, Hamming-Vrieze O, van Vliet-Vroegindeweij C, van de Kamer JB, et al. FDG-PET and diffusion-weighted MRI in head-and-neck cancer patients: implications for dose painting. Radiother Oncol. (2013) 106:250–4. doi: 10.1016/j.radonc.2013.01.003

PubMed Abstract | CrossRef Full Text | Google Scholar

22. Daisne JF, Sibomana M, Bol A, Doumont T, Lonneux M, Gregoire V. Tri-dimensional automatic segmentation of PET volumes based on measured source-to-background ratios: influence of reconstruction algorithms. Radiother Oncol. (2003) 69:247–50. doi: 10.1016/S0167-8140(03)00270-6

PubMed Abstract | CrossRef Full Text | Google Scholar

23. Zhao M, Liu Z, Sha Y, Wang S, Ye X, Pan Y, et al. Readout-segmented echo-planar imaging in the evaluation of sinonasal lesions: a comprehensive comparison of image quality in single-shot echo-planar imaging. Magn Reson Imaging. (2016) 34:166–72. doi: 10.1016/j.mri.2015.10.010

PubMed Abstract | CrossRef Full Text | Google Scholar

24. Cao Y, Li D, Shen Z, Normolle D. Sensitivity of quantitative metrics derived from DCE MRI and a pharmacokinetic model to image quality and acquisition parameters. Acad Radiol. (2010) 17:468–78. doi: 10.1016/j.acra.2009.10.021

PubMed Abstract | CrossRef Full Text | Google Scholar

25. Newitt DC, Malyarenko D, Chenevert TL, Quarles CC, Bell L, Fedorov A, et al. Multisite concordance of apparent diffusion coefficient measurements across the NCI Quantitative Imaging Network. J Med Imaging. (2018) 5:011003. doi: 10.1117/1.JMI.5.1.011003

PubMed Abstract | CrossRef Full Text | Google Scholar

26. Aryal M, Lee C, Hawkins PG, Chapman C, Eisbruch A, Mierzwa M, et al. Real-time quantitative assessment of accuracy and precision of blood volume derived from DCE MRI in individual patients during a clinical trial. Tomography. (2019) 5:61–7. doi: 10.18383/j.tom.2018.00029

PubMed Abstract | CrossRef Full Text | Google Scholar

27. Chung CH, Zhang Q, Kong CS, Harris J, Fertig EJ, Harari PM, et al. p16 protein expression and human papillomavirus status as prognostic biomarkers of nonoropharyngeal head and neck squamous cell carcinoma. J Clin Oncol. (2014) 32:3930–8. doi: 10.1200/JCO.2013.54.5228

PubMed Abstract | CrossRef Full Text | Google Scholar

28. Lassen P, Lacas B, Pignon JP, Trotti A, Zackrisson B, Zhang Q, et al. Prognostic impact of HPV-associated p16–expression and smoking status on outcomes following radiotherapy for oropharyngeal cancer: the MARCH-HPV project. Radiother Oncol. (2018) 126:107–15. doi: 10.1016/j.radonc.2017.10.018

PubMed Abstract | CrossRef Full Text | Google Scholar

29. Rosenthal DI, Harari PM, Giralt J, Bell D, Raben D, Liu J, et al. Association of human papillomavirus and p16 status with outcomes in the IMCL-9815 phase III registration trial for patients with locoregionally advanced oropharyngeal squamous cell carcinoma of the head and neck treated with radiotherapy with or without cetuximab. J Clin Oncol. (2016) 34:1300–8. doi: 10.1200/JCO.2015.62.5970

CrossRef Full Text | Google Scholar

30. Cohen EE, Karrison TG, Kocherginsky M, Mueller J, Egan R, Huang CH, et al. Phase III randomized trial of induction chemotherapy in patients with N2 or N3 locally advanced head and neck cancer. J Clin Oncol. (2014) 32:2735–43. doi: 10.1200/JCO.2013.54.6309

CrossRef Full Text | Google Scholar

31. Hawkins PG, Mierzwa ML, Bellile E, Jackson WC, Malloy KM, Chinn SB, et al. Impact of American Joint Committee on Cancer Eighth Edition clinical stage and smoking history on oncologic outcomes in human papillomavirus-associated oropharyngeal squamous cell carcinoma. Head Neck. (2019) 41:857–64. doi: 10.1002/hed.25336

PubMed Abstract | CrossRef Full Text | Google Scholar

32. Spector ME, Chinn SB, Bellile E, Gallagher KK, Ibrahim M, Vainshtein J, et al. Matted nodes as a predictor of distant metastasis in advanced-stage III/IV oropharyngeal squamous cell carcinoma. Head Neck. (2016) 38:184–90. doi: 10.1002/hed.23882

PubMed Abstract | CrossRef Full Text | Google Scholar

33. Martens RM, Noij DP, Koopman T, Zwezerijnen B, Heymans M, de Jong MC, et al. Predictive value of quantitative diffusion-weighted imaging and 18-F-FDG-PET in head and neck squamous cell carcinoma treated by (chemo)radiotherapy. Eur J Radiol. (2019) 113:39–50. doi: 10.1016/j.ejrad.2019.01.031

PubMed Abstract | CrossRef Full Text | Google Scholar

34. Lin TA, Garden AS, Elhalawani H, Elgohari B, Jethanandani A, Ng SP, et al. Radiographic retropharyngeal lymph node involvement in human papillomavirus-associated oropharyngeal carcinoma: patterns of involvement and impact on patient outcomes. Cancer. (2019) 125:1536–46. doi: 10.1002/cncr.31944

CrossRef Full Text | Google Scholar

35. Driessen JP, Caldas-Magalhaes J, Janssen LM, Pameijer FA, Kooij N, Terhaard CH, et al. Diffusion-weighted MR imaging in laryngeal and hypopharyngeal carcinoma: association between apparent diffusion coefficient and histologic findings. Radiology. (2014) 272:456–63. doi: 10.1148/radiol.14131173

PubMed Abstract | CrossRef Full Text | Google Scholar

36. Barajas RF Jr, Rubenstein JL, Chang JS, Hwang J, Cha S. Diffusion-weighted MR imaging derived apparent diffusion coefficient is predictive of clinical outcome in primary central nervous system lymphoma. AJNR Am J Neuroradiol. (2010) 31:60–6. doi: 10.3174/ajnr.A1750

PubMed Abstract | CrossRef Full Text | Google Scholar

37. Gibbs P, Liney GP, Pickles MD, Zelhof B, Rodrigues G, Turnbull LW. Correlation of ADC and T2 measurements with cell density in prostate cancer at 3.0 Tesla. Invest Radiol. (2009) 44:572–6. doi: 10.1097/RLI.0b013e3181b4c10e

PubMed Abstract | CrossRef Full Text | Google Scholar

38. Wang XZ, Wang B, Gao ZQ, Liu JG, Liu ZQ, Niu QL, et al. Diffusion-weighted imaging of prostate cancer: correlation between apparent diffusion coefficient values and tumor proliferation. J Magn Reson Imaging. (2009) 29:1360–6. doi: 10.1002/jmri.21797

PubMed Abstract | CrossRef Full Text | Google Scholar

39. Zelhof B, Pickles M, Liney G, Gibbs P, Rodrigues G, Kraus S, et al. Correlation of diffusion-weighted magnetic resonance data with cellularity in prostate cancer. BJU Int. (2009) 103:883–8. doi: 10.1111/j.1464-410X.2008.08130.x

PubMed Abstract | CrossRef Full Text | Google Scholar

40. Neesse A, Michl P, Frese KK, Feig C, Cook N, Jacobetz MA, et al. Stromal biology and therapy in pancreatic cancer. Gut. (2011) 60:861–8. doi: 10.1136/gut.2010.226092

PubMed Abstract | CrossRef Full Text | Google Scholar

41. Mueller MM, Fusenig NE. Friends or foes - bipolar effects of the tumour stroma in cancer. Nat Rev Cancer. (2004) 4:839–49. doi: 10.1038/nrc1477

PubMed Abstract | CrossRef Full Text | Google Scholar

42. Wang K, Ma W, Wang J, Yu L, Zhang X, Wang Z, et al. Tumor-stroma ratio is an independent predictor for survival in esophageal squamous cell carcinoma. J Thorac Oncol. (2012) 7:1457–61. doi: 10.1097/JTO.0b013e318260dfe8

PubMed Abstract | CrossRef Full Text | Google Scholar

43. de Kruijf EM, van Nes JG, van de Velde CJ, Putter H, Smit VT, Liefers GJ, et al. Tumor-stroma ratio in the primary tumor is a prognostic factor in early breast cancer patients, especially in triple-negative carcinoma patients. Breast Cancer Res Treat. (2011) 125:687–96. doi: 10.1007/s10549-010-0855-6

PubMed Abstract | CrossRef Full Text | Google Scholar

44. Huijbers A, Tollenaar RA, v Pelt GW, Zeestraten EC, Dutton S, McConkey CC, et al. The proportion of tumor-stroma as a strong prognosticator for stage II and III colon cancer patients: validation in the VICTOR trial. Ann Oncol. (2013) 24:179–85. doi: 10.1093/annonc/mds246

PubMed Abstract | CrossRef Full Text | Google Scholar

45. Wiseman BS, Werb Z. Stromal effects on mammary gland development and breast cancer. Science. (2002) 296:1046–9. doi: 10.1126/science.1067431

PubMed Abstract | CrossRef Full Text | Google Scholar

46. Ward MJ, Thirdborough SM, Mellows T, Riley C, Harris S, Suchak K, et al. Tumour-infiltrating lymphocytes predict for outcome in HPV-positive oropharyngeal cancer. Br J Cancer. (2014) 110:489–500. doi: 10.1038/bjc.2013.639

PubMed Abstract | CrossRef Full Text | Google Scholar

47. Chakravarthy A, Henderson S, Thirdborough SM, Ottensmeier CH, Su X, Lechner M, et al. Human papillomavirus drives tumor development throughout the head and neck: improved prognosis is associated with an immune response largely restricted to the oropharynx. J Clin Oncol. (2016) 34:4132–41. doi: 10.1200/JCO.2016.68.2955

PubMed Abstract | CrossRef Full Text | Google Scholar

48. Bogowicz M, Leijenaar RTH, Tanadini-Lang S, Riesterer O, Pruschy M, Studer G, et al. Post-radiochemotherapy PET radiomics in head and neck cancer - The influence of radiomics implementation on the reproducibility of local control tumor models. Radiother Oncol. (2017) 125:385–91. doi: 10.1016/j.radonc.2017.10.023

PubMed Abstract | CrossRef Full Text | Google Scholar

49. Ou D, Blanchard P, Rosellini S, Levy A, Nguyen F, Leijenaar RTH, et al. Predictive and prognostic value of CT based radiomics signature in locally advanced head and neck cancers patients treated with concurrent chemoradiotherapy or bioradiotherapy and its added value to Human Papillomavirus status. Oral Oncol. (2017) 71:150–5. doi: 10.1016/j.oraloncology.2017.06.015

PubMed Abstract | CrossRef Full Text | Google Scholar

50. Aerts HJ, Velazquez ER, Leijenaar RT, Parmar C, Grossmann P, Carvalho S, et al. Decoding tumour phenotype by noninvasive imaging using a quantitative radiomics approach. Nat Commun. (2014) 5:4006. doi: 10.1038/ncomms5006

PubMed Abstract | CrossRef Full Text | Google Scholar

51. Jansen JF, Lu Y, Gupta G, Lee NY, Stambuk HE, Mazaheri Y, et al. Texture analysis on parametric maps derived from dynamic contrast-enhanced magnetic resonance imaging in head and neck cancer. World J Radiol. (2016) 8:90–7. doi: 10.4329/wjr.v8.i1.90

PubMed Abstract | CrossRef Full Text | Google Scholar

52. Zhang H, Graham CM, Elci O, Griswold ME, Zhang X, Khan MA, et al. Locally advanced squamous cell carcinoma of the head and neck: CT texture and histogram analysis allow independent prediction of overall survival in patients treated with induction chemotherapy. Radiology. (2013) 269:801–9. doi: 10.1148/radiol.13130110

PubMed Abstract | CrossRef Full Text | Google Scholar

Keywords: MRI, head and neck cancer, radiation therapy, imaging biomarker, adaptive therapy

Citation: Cao Y, Aryal M, Li P, Lee C, Schipper M, Hawkins PG, Chapman C, Owen D, Dragovic AF, Swiecicki P, Casper K, Worden F, Lawrence TS, Eisbruch A and Mierzwa M (2019) Predictive Values of MRI and PET Derived Quantitative Parameters for Patterns of Failure in Both p16+ and p16– High Risk Head and Neck Cancer. Front. Oncol. 9:1118. doi: 10.3389/fonc.2019.01118

Received: 01 August 2019; Accepted: 08 October 2019;
Published: 14 November 2019.

Edited by:

John Varlotto, University of Massachusetts Medical School, United States

Reviewed by:

Vivek Verma, Allegheny General Hospital, United States
Bilgin Kadri Aribas, Bülent Ecevit University, Turkey

Copyright © 2019 Cao, Aryal, Li, Lee, Schipper, Hawkins, Chapman, Owen, Dragovic, Swiecicki, Casper, Worden, Lawrence, Eisbruch and Mierzwa. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Yue Cao, eXVlY2FvJiN4MDAwNDA7bWVkLnVtaWNoLmVkdQ==

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.