Theranostics 2020; 10(19):8677-8690. doi:10.7150/thno.47298 This issue
1. Department of Surgery and Cancer, Imperial College London, United Kingdom.
2. Department of Radiology/Nuclear Medicine, Imperial College Healthcare NHS Trust, London, United Kingdom.
3. Department of Surgery and Cancer, Imperial College Healthcare NHS Trust, London, United Kingdom.
4. Department of Medicine and Integrated Care, Imperial College Healthcare NHS Trust, London, United Kingdom.
Purpose: The spatio-molecular distribution of choline and its metabolites in tumors is highly heterogeneous. Due to regulation of choline metabolism by hypoxic transcriptional signaling and other survival factors, we envisage that detection of such heterogeneity in patient tumors could provide the basis for advanced localized therapy. However, non-invasive methods to assess this phenomenon in patients are limited. We investigated such heterogeneity in Non-Small Cell Lung Cancer (NSCLC) with [18F]fluoromethyl-(1,2-2H4) choline ([18F]D4-FCH) and positron emission tomography/computed tomography (PET/CT).
Experimental design: [18F]D4-FCH (300.5±72.9MBq [147.60-363.6MBq]) was administered intravenously to 17 newly diagnosed NSCLC patients. PET/CT scans were acquired concurrently with radioactive blood sampling to permit mathematical modelling of blood-tissue transcellular rate constants. Comparisons were made with biopsy-derived choline kinase-α (CHKα) expression and diagnostic [18F]fluorodeoxyglucose ([18F]FDG) scans.
Results: Oxidation of [18F]D4-FCH to [18F]D4-fluorobetaine was suppressed (48.58±0.31% parent at 60 min) likely due to the deuterium isotope effect embodied within the design of the radiotracer. Early (5 min) and late (60 min) images showed specific uptake of tracer in all 51 lesions (tumors, lymph nodes and metastases) from 17 patients analyzed. [18F]D4-FCH-derived uptake (SUV60max) in index primary lesions (n=17) ranged between 2.87-10.13; lower than that of [18F]FDG PET [6.89-22.64]. Mathematical modelling demonstrated net irreversible uptake of [18F]D4-FCH at steady-state, and parametric mapping of the entire tumor showed large intratumorally heterogeneity in radiotracer retention, which is likely to have influenced correlations with biopsy-derived CHKα expression.
Conclusions: [18F]D4-FCH is detectable in NSCLC with large intratumorally heterogeneity, which could be exploited in the future for targeting localized therapy.
Keywords: [18F] Fluoromethyl-(1, 2-2H4)-choline, Positron Emission Tomography, heterogeneity choline metabolism, choline kinase alpha, lung cancer.
Choline radiotracers, including [11C]choline, [18F]fluoromethylcholine ([18F]FCH) and [18F]fluoroethylcholine ([18F]FEC) are approved for use in detecting recurrent prostate cancer in patients with rising PSA despite having received definitive treatment together with investigational use in cancers of the brain, lung, liver and bladder . Evolving literature also suggests potential utility in response monitoring . Tumor PET images represent transport, phosphorylation and retention of choline metabolites in tumor cells . [11C]choline and [18F]FCH are rapidly oxidized to [11C]betaine/[18F]fluorobetaine by choline oxidase mainly in kidney and liver tissues, with metabolites detectable in plasma soon after injection of the radiotracer [4-7]. To preclude interpretation of complex imaging signals, we developed a more metabolically stable fluorocholine analogue, [18F]fluoromethyl-(1, 2-2H4)choline ([18F]D4-FCH) based on the deuterium isotope effect [6, 8, 9].
Radiotherapy, alone or in combination with systemic therapies are the mainstay of treatment in patients diagnosed with stage I-IV NSCLC, and with advances in radiotherapy planning and delivery, a remarkable improvement in local tumor control has been seen [10, 11]. It is hypothesized that imaging approaches for detecting the multiple factors that influence radioresistance would make further impact. For instance, hypoxia/hypoxia signaling is known to remarkably influence radioresistance, which has led to investigation of pre-therapy [18F]fluoromisonidazole PET ([18F]FMISO PET) to identify aggressive hypoxic components of tumors which can be targeted with a radiation boost with aim of reducing recurrence . In view of logistic issues with the latter e.g. 2-3 h to allow for (a relatively low) imaging contrast, other strategies are warranted. Choline transport and metabolism are upregulated in cancer cells and serve as potential targets in PET imaging using radiolabeled choline tracers. Choline-related metabolites including phosphocholine display a highly heterogeneous spatio-molecular distribution in tumors by high resolution mass spectrometry imaging, and expression of choline kinase-α (CHKα) is regulated transcriptionally by hypoxia . In early-stage NSCLC CHKα mRNA expression is prognostic indicating that aspects of tumor aggressiveness regulate the heterogeneous choline metabolic phenotype . The most widely used radiotracer for tumor imaging is [18F]fluorodeoxyglucose ([18F]FDG) which measures glucose utilization. [18F]FDG does not overlap with hypoxia measured by [18F]FMISO PET at the global or voxel level.
Having developed [18F]D4-FCH for late (60 min) imaging of choline transport and phosphorylation in pre-clinical models and demonstrated safety/appropriate organ dosimetry in healthy volunteers, the objective of this study was to investigate imaging characteristics of the tracer in tumors and metastases, assess tumor heterogeneity and further confirm whether simplified measures of radiotracer uptake were associated with more complex kinetic variables so as to support use of such simplified measures; particular for standard late (60 min) imaging where, for [11C]choline and [18F]FCH, parent radiotracer levels in blood are either very low or undetectable. The setting for this first human study is non-small cell lung cancer (NSCLC). The choice is predicated on previous work suggesting that CHKα expression correlates with survival ; and routine use of [18F]FDG PET in this setting, to enable comparison.
A prospective non-randomized feasibility study was conducted in patients with newly diagnosed localized, locally advanced or metastatic NSCLC who were fit for surgery, radiotherapy or systemic therapy. Patients with histologically confirmed NSCLC were recruited from Respiratory and Oncology clinics at Imperial College NHS Healthcare Trust over a period of 3 years (2014-2017). All patients included were ≥18 yr., good performance status (ECOG ≥2), and had clinically acceptable blood parameters. All had at least one measurable site of disease ≥ 2cm (primary tumor or lymph node). Excluded from study were patients participating in other clinical trials with investigational drugs within 30 days prior to study entry, pregnant or lactating women and patients diagnosed with an invasive malignancy within the last 5 years (excluding basal cell carcinoma). Ethical approval was granted for this study by the West London and GTAC NRES committee, and the Medicine and Healthcare Products Regulatory Agency (MHRA). All patients provided written informed consent in accordance with the Declaration of Helsinki and the administration of radioactivity was approved by the Administration of Radioactive Substances Advisory Committee (ARSAC).
[18F]D4-FCH was synthesized from its precursor as previously described by our group, 2011 [16, 17]. As part of standard investigations, patients underwent routine imaging with contrast enhanced computed tomography (CECT) of chest and upper abdomen and a staging [18F]FDG PET/CT. [18F]D4-FCH PET/CT was performed after routine standard investigations and prior to any therapy. Images were acquired on a Siemens Biograph 6 TruePoint PET/CT scanner (with TrueV; extended field of view) with 21.6-cm axial and 60.5-cm transaxial fields of view. A thoracic CT scan (130 kV; exposure, 30 effective mAs; pitch, 1.5; slice thickness, 5 mm; rotation time, 0.6 s; resulting in an effective dose of 2.5 mSv) was performed for attenuation correction and anatomical localization. This was then followed by tracer injection (a maximum dose of 370MBq of [18F]D4-FCH) and dynamic imaging (single-bed position over thorax and mediastinum) for 66 min followed by a half-body static (mid-thigh to vertex) attenuation CT and PET scan (3 min per bed position). Discrete venous bloods and plasma (at 5, 10, 15, 30, 60 min post injection) were obtained for radioactivity counting and metabolite analysis, as previously described [16, 17] and applied in oncological studies [18, 19]. The time course of metabolite and parent fraction measured directly by HPLC (experimental data) for D4 Choline are summarized in Figure S1.
All dynamic [18F]D4-FCH PET/CT data were reconstructed using the ordered-subsets expectation maximization algorithm (OSEM - 3 iterations and 21 subsets) with corrections for dead time, scatter, attenuation, and radioactive decay. For each patient, volumes of interest (VOIs) for both [18F]D4-FCH and [18F]FDG PET/ CT were generated using Hermes (Hermes Diagnostics, Stockholm, Sweden) by a single investigator (SD) onto fused image datasets (PET/CT). Volumes of interests (VOIs) were drawn manually over the whole tumor, in normal background organs (2 cm fixed spheres in normal lung tissue, vertebral body and muscle), in the aorta (1cm fixed sphere) and in lymph nodes, when present. An additional VOI was drawn around the aorta (1cm fixed sphere) for an image derived calculation of the arterial input function. The parent plasma input function was generated by normalizing to plasma-over-blood ratio and correcting for metabolites. In addition to that, to account for the presence of labelled metabolites, the metabolite time activity curve in plasma was used as input of the double input spectral analysis . PET data was quantified with the semi-quantitative Standardized Uptake Value (SUV) and tracer kinetics was studied with quantitative methods: spectral, graphical and compartmental analyses. The work by Boellaard's lab for [18F]FCH  indicated discordance between SUV and kinetic variables, which warranted investigation of tracer kinetics for [18F]D4-FCH. PET data were quantified with the semi-quantitative Standardized Uptake Value (SUV) and tracer kinetics were studied with quantitative methods: spectral, graphical and compartmental analyses. The analysis started with the implementation of Spectral Analysis (SA) which, without any a priori assumption on the kinetics of the tracer, provides a “kinetic spectrum” representing the functional processes in which the tracer is involved. Hence, from this spectrum it is possible to obtain a complete description of tracer kinetics. It is also possible to identify the number and the type of compartments necessary for the data modelling- necessary information for the following graphical and compartmental analyses.
SUV: PET quantification is routinely performed using the SUV, a semi-quantitative index computed as the raw image counts normalized by the injected dose, and generally, the patient's body weight. From SUV tissue radioactivity curves, we evaluated mean and maximum values of SUV and areas under the SUV curve (SUVAUC; from zero to time t) at 5, 30 and 60 min (SUV5mean, SUV30mean, SUV60mean, SUV5max, SUV30max, SUV60max, SUVAUC5mean, SUVAUC30mean, SUVAUC60mean, SUVAUC5max, SUVAUC30max, SUVAUC60max). We also evaluated mean and maximum SUV for the later static [18F]D4-FCH scan (SUV60mean_s, SUV60max_s).
Spectral analysis (SA): Without any a priori assumption on the kinetics of the tracer, SA provided the K1 (K1_SA), the influx rate constant, and the Ki (Ki_SA), net uptake of tracer in tissue. Spectral analysis provided also the impulse response function (IRF) of the compartmental models used in PET, evaluated as an analytical sum of exponentials . SA was conducted for each patient choosing a grid of 100 components in the range of (0.0001, 1) [23, 24]. A double input SA (DI-SA) was also implemented to account for the contribution of the radiometabolites. This approach is based on the method implemented by Tomasi et al. [20, 25].
Graphical analysis: Graphical methods are the simplest approaches for dynamic PET data quantification as they are based on a linearization of compartmental model's differential equations. Patlak plot was used to estimate Ki (Ki_P), the net irreversible uptake rate constant, which quantifies the rate at which the tracer is irreversibly trapped . The Logan plot is the counterpart of the Patlak plot for reversible radiotracers. This method was used to estimate the single macroparameter Vt (Vt_L), the total volume of distribution of the tracer in the tissue .
Compartmental analysis: To investigate the best kinetic model, methods by Fan and colleagues , were employed with the MICK software (modelling-Input-function-Compartmental-Kinetics) and Matlab (Mathworks, R2015b). Three different kinetic models were tested: the reversible 1-tissue 2k (1TCM2k) model, the 2-tissue 3k (2TCM4k) model, and the irreversible 2-tissue 3k model (for which k implies the rate constant for tracer for different kinetic compartments).
Parametric maps and distribution kernels: A Kernel function was used to model the voxel distribution profiles of each parameter and to obtain skewness, kurtosis and ratio (skewness/kurtosis) values for each parameter.
The model providing the best fits to the lesion time-activity curves was selected on the basis of the Akaike information criterion (AIC) :
where m denotes number of parameters, n equals the sum of degrees of freedom and number of parameters, and wrss denotes weighted residual sum of squares .
The Kolmogorov-Smirnov test was used to assess the normal distribution of each parameter distribution. Differences among tumor and healthy tissues were tested with the Wilcoxon Sum Rank test and corrected for multiple comparisons using Bonferroni approach. Spearman's rank test was used to verify the correlation between parameters. All statistical tests were run in Matlab (Mathworks, R2018b).
CHKα immunohistochemistry was conducted as previously described studies using a primary polyclonal, human anti-CHKα antibody as per manufacturer's instructions (Sigma-Aldrich, HPA024153,1:20 dilution); the intensity of cytoplasmic and nuclear staining were scored (score 1: low intensity, 2: moderate and 3: high) [31, 32]. The relationship between PET uptake parameters (SUV60mean and SUV60max) and CHKα expression, determined by immunohistochemistry on lung tumor tissue and lymph node specimens, was established. In addition, commercially available tissue microarrays (TMAs; US Biomax Inc.) were obtained (HLug-Ade090Lym-01; HLug-Squ090Lym-01; HLug-Ade150Sur-01; HLug-Squ150Sur-01) and staining performed as above.
Seventeen patients with confirmed NSCLC were enrolled to the study. All patients fulfilled the inclusion criteria and provided written informed consent. Patient characteristics and primary treatment modalities are shown (Table 1). The mean age ± sd were 64 ± 11y; age range 39-84 y; weight ± sd, 68.9 ± 15kg; weight range, 49.7-93.3 kg), respectively. Nine patients were diagnosed with squamous cell carcinoma and 8 had adenocarcinoma. Twelve out of 17 patients had nodal metastases diagnosed on mediastinoscopy and confirmed as positive on tumor and nodal resection on surgical histopathology. Five patients eligible for surgery underwent a lobectomy or pneumonectomy, or wedge resection with hilar and mediastinal lymph node sampling, or en bloc resection. Primary systemic treatment consisted of platinum doublet chemotherapy (gemcitabine and carboplatin, pemetrexed and cisplatin or carboplatin). Other treatment options included radical radiotherapy, chemoradiotherapy, and palliative high-dose radiotherapy. Treatment decisions were made by the lung cancer multidisciplinary team. Of the adenocarcinoma patients tested, 6/6 were negative for EGFR (epidermal growth factor receptor) mutation, 2/3 were ALK (Anaplastic lymphoma kinase gene) positive, and 1/5 were KRAS (Kristen rat sarcoma viral oncogene) mutant. Median follow-up for all patients was 15 months (range 3-47 months). Eleven patients out of 17 died of lung cancer at time of last follow-up. Local and distant recurrence were identified in 9 patients.
|Patient ID||Histology||TNM||Primary treatment|
|2||Squamous||T4N2M1b||Chemotherapy, high dose palliative RT|
|7||Squamous||T1bN1M0||Surgery and adjuvant chemotherapy|
Radiochemical purity of [18F]D4-FCH was 100% on completion of synthesis with a mean (± sd) specific activity of 61 ± 76 GBq/μmol (range, 21-353 GBq/μmol) and pH of 5.35 ± 0.45 (range 4.50 -5.81). The mean (± sd) dose injected in all patients was 300.5 ± 72. 9MBq (range 147.60 - 363.6MBq). [18F]D4-FCH was well tolerated by all patients with no immediate or delayed complications observed. As per healthy volunteers , the radiotracer was stable after injection into patients with only 48.6% of the parent converting mainly to a single metabolite ([18F]D4-fluorobetaine) after 60 min. Figure 1A shows the chemical structure of [18F]D4-FCH, the typical radio-high-performance liquid chromatograms (radio-HPLC, 1200 series system; Agilent) for a representative patient (Figure 1B), and the mean percentage metabolite and parent tracer fraction over time (Figure 1C).
Structure and metabolic stability of [18F]D4-FCH. (A) Chemical structure of [18F]fluoromethyl-(1,2-2H4) choline ([18F]D4-FCH). (B) Typical high-performance liquid chromatogram (HPLC) of [18F]D4- FCH in plasma at 2.5, 5, 10 and 30 min time-points in a patient. Red arrows point to parent tracer and blue arrows point to metabolite peak (betaine) of [18F]D4-FCH. Y-axis scale adjusted to highlight metabolite peak. CPM = counts per minute. (C) Mean percentage parent tracer and metabolite fraction over time (min).
Images acquired with [18F]D4-FCH PET, and comparison to [18F]FDG PET are shown (Figures 2A and B). Mean decay-corrected time activity curves (TACs) for tumor, lymph nodes and normal tissues are shown in Figure 2C (SUV at 5, 30 and 60 min were similar, Table 2; Table S1). Physiological uptake of [18F]D4-FCH PET was noted in the salivary glands, liver, kidneys, pancreas and bladder. All primary tumors and involved lymph nodes were visible on [18F]D4-FCH PET/CT. NSCLC lesions were avid for both [18F]D4-FCH PET and [18F]FDG in patients who underwent [18F]FDG PET/CT as part of their staging work-up. [18F]D4-FCH derived uptake (SUVmax) in index primary lesions (n=17) ranged between 2.87 and 10.13; lower than that of [18F]FDG PET/CT (6.89 and 22.64) (Table 2; Table S1). On analysis of [18F]D4-FCH PET/CT (static) alone, the SUVmean and SUVmax values were higher in the lymph nodes and metastases compared with primary tumor in 7 out of 17 patients. Corresponding nodal lesions on [18F]FDG PET/CT (where available) showed in general, higher SUV values but in 1 patient (patient 3), SUV was found to be lower than that of [18F]D4-FCH PET/CT. The noted [18F]D4-FCH PET uptake within the mediastinal lymph nodes was negative for malignancy on correlation with surgical histopathology. It is unclear whether relevant or not, but on further investigation this patient was noted to be a heavy smoker. Unlike prostate cancer where low level symmetrical [18F]FCH uptake is often seen in reactive nodes , this was not seen with [18F]D4-FCH PET. The same number of lesions were noted with [18F]FDG PET/CT apart from 1 patient where [18F]FDG demonstrated uptake of a lower vertebral metastasis (patient 6).
The behavior of the tracer was investigated in a stepwise manner. With spectral analysis the kinetics of [18F]D4-FCH was described by a kinetic spectrum, which allowed the identification of the number and type of compartments necessary for the following graphical and compartmental analyses. On the basis of the AIC, standard and double input spectral analyses  were comparable (AIC of standard SA = 4.98 ± 34.76, AIC of DI-SA = 7.17 ± 34.61) with slightly better fits provided by DI-SA. SA revealed irreversible trapping component of [18F]D4-FCH in the tissue in 8 patients (Figure 3A, Figure S2). This was further confirmed by the fractional Retention of Tracer (FRT) calculated as IRF at 60 minutes relative to 1 min (p > 0.05) (Figure 3B) and by compartmental analysis (fitted to a two-tissue irreversible model - 2TCM3k - p < 0.05). Both Patlak and Logan plots are justified by the result of the spectral analysis. In fact, while the spectrum revealed a low frequency slow (irreversible) component, this did not occur in the entire cohort. SA revealed also a reversible trapping component in 3 patients. For this reason, both Patlak and Logan plot were implemented. The best fit was provided by Patlak plot (AIC = -57.00 ± 25.39), followed by Logan plot (AIC = 101.06 ± 37.22). Mean SA K1 (K1_SA), SA Ki (Ki_SA), SA parent tracer volume of distribution (VTpar), SA metabolite volume of distribution (VTmetab), DI-SA Ki, Patlak Ki (Ki_P), and Logan Vt (Vt_L) in primary tumors and healthy lung are summarized in Table 3. SA revealed irreversible trapping component of [18F]D4-FCH in the tissue in 8 patients (Figure 3A, Figure S2).
Patient characteristics and PET parameters in primary tumor.
|Patient||Sex||TNM||Histology||[18F]D4-FCH Dynamic Scan||[18F]D4-FCH Static Scan||[18F]FDG Static Scan|
|SUV 5 min||SUV 30 min||SUV 60 min||SUV 60 min||SUV 60 min|
M, Male; F, Female
Results of quantitative analyses of dynamic PET data. Mean (± sd) values for Patlak Ki (Ki_P) [mL/cm3/min], SA K1 (K1_SA) [mL/cm3/min], SA Ki (Ki_SA) [mL/cm3/min], Logan Vt (Vt_L) [mL/cm3] and the compartmental model Ki (Ki_CM) [mL/cm3/min] in primary tumor and healthy lung. A double input SA was used to investigate the volumes of distribution of the parent and of the metabolite in the lesions. VT_par and VT_metab are the volumes of distribution for the parent and the metabolite, respectively. Parent and metabolite both contribute to the final PET signal.
|Single Input SA||Double Input SA|
|Primary tumor||0.05 (±0.04)||6.38 (±7.56)||8E-05 (±9E-05)||0.26 (±0.31)||0.99 (±0.70)||0.37 (±0.63)||8.75E-05 (±1E-04)|
|Healthy lung||0.0 2 (±0.22)||2.87 (±2.32)||4E-05 (±5E-05)||0.07 (±0.99)|
We previously demonstrated in mouse tumors that [18F]D4-FCH is phosphorylated to [18F]D4-phosphocholine or oxidized to [18F]D4-fluorobetaine, (with the latter occurring mainly in organs like liver or kidneys)[9, 25]. However, we did not conduct a head-to-head comparison between [18F]FCH and [18F]D4-FCH, and to date, only rodent studies have been conducted in a head-to-head manner. DI-SA was used to investigate the volumes of distribution of the parent and of the metabolite in the tissue. Table 3 shows the two volumes of distribution and the Ki_SA in primary tumors and lymph nodes. In most lesions with measurable trapping component, the main contributor to the distribution volume was the parent fraction; the higher contribution of metabolites to certain tumors with low tracer retention was also evident.
Diagnostic biopsy tissues from primary tumor were available in 10 out of 17 patients, 6 of whom were positive for CHKα protein expression (Figure 4A and B) (lymph node biopsies were not available). No correlation was seen between PET uptake variables and CHKα protein expression, e.g., SUVmax and CHKα expression (p = 0.968). Of interest, two out of the 6 patients, who demonstrated higher [18F]D4-FCH uptake in lymph node metastases than primary tumor (patient 2 and 10), were negative for CHKα expression indicating tissue heterogeneity. To affirm that staining heterogeneity was not just pertinent to our prospective study, we investigated CHKα protein expression in commercially available TMAs (US Biomax Inc.). Sixty cores of squamous cell carcinoma and 60 cores of adenocarcinoma with no survival data, as well as 75 cores of each histopathology sub-type with corresponding survival data were available. Analysis of these tissue demonstrated substantial heterogeneity (~50% staining) of CHKα protein expression in lung lesions of both histologies (Figure 4C) supporting data from our prospective study. 48.95% of cores (35/69 squamous cell and 35/75 adenocarcinoma) showed cytoplasmic staining for CHKα. CHKα expression with staining intensity +2 was found in all T2a/b, N0/N1 and M0 cases. No correlation with survival and stage of disease was seen.
Initially we performed histogram analysis of 3D parametric maps of Ki_P, K1_SA and Vt_L generated from Graphical and Spectral methods for tumor, normal lung (as reference tissue), lymph node metastases and vertebra (Figure 5A). Statistical differences between tumor and healthy tissues (entire volume-of-interest) were tested with the Wilcoxon rank sum test and were found to be significant (p < 0.05) for Ki_P, K1_SA, Ki_SA and Vt_L. To interrogate spatial relationships, we constructed histograms of the data. Figure 5B shows the histogram analysis K1_SA and Ki_SA for tumor, normal lung, lymph node metastases and vertebra in a representative patient. In the majority of analyzed cases, the voxel distribution profiles were non-Gaussian and modelling was performed with a Kernel function . From the results of distribution analysis, the flattened tumor distributions profiles (Figure 5B) testify the wider range of values covered by parametric maps (theoretically from higher skewness, lower kurtosis and/or higher skewness/kurtosis ratio relative to normal lung). Notably, a statistically significant lower kurtosis was seen for the K1_SA kernel for CHKα negative tumors compared to CHKα positive tumors (Figure 5C and D; Table S2), indicating perhaps that one reason for CHKα negativity may be due to heterogenous sampling.
The imaging protocol could be simplified using semi-quantitative methods instead of a full kinetic analysis. Results from SUV analyses and the quantification of dynamic PET data were compared. In particular, we compared Ki_P, K1_SA, Ki_SA, Ki_CM and Vt_L to the mean and maximum SUV at 5 and 60 minutes ( SUV5mean, SUV5max, SUV60mean, SUV60max), and to the SUVAUC mean and maximum evaluated at 5 and 60 minutes (SUVAUC_5mean, SUVAUC_5max, SUVAUC_60mean, SUCAUC_60max). Results of the linear regression analysis are summarized in Table S3. Commonly used simplified methods showed moderate correlation to Ki_P and K1_SA (Ki_P vs. SUV60mean, R2 = 0.5308; Ki_P vs. SUVAUC_60mean, R2 = 0.5246; K1_SA vs. SUV60mean, R2 = 0.4617; K1_SA vs. SUVAUC_60mean, R2 = 0.4502). The lowest correlations resulted from the comparison between SUVs and Ki_SA, Ki_CM and Vt_L.
[18F]D4-FCH and [18F]FDG Tissue uptake. (A) Axial views on CT, fused PET/CT and PET for diagnostic [18F]FDG and [18F]D4-FCH in a patient with squamous cell carcinoma of the right lung (patient 3 - see Table 2); (B) similar images in a patient with an adenocarcinoma of the right lung (patient 13 - see Table 2), SUV intensity threshold 0-10; and; (C) Mean decay-corrected time-activity curves (TACs) for normal tissue (normal lung, bone, including bone marrow, and muscle) and tumor and lymph nodes in dynamic PET imaging.
Spectral analysis of [18F]D4-FCH in lung tumor. (A) The unit impulse response function (IRF) of tumor compared with reference tissue (normal lung), lymph node metastases and vertebra (values of the fractional retention of the tracer (FRT) are shown) and spectrum of the kinetic components obtained using spectral analysis. (B) Box-and-whisker plots of FRT showing the highest retention component in tumor and lymph node metastases. Less retention was seen in the reference tissue (normal lung) and vertebra. Node refers to involved lymph node metastases.
Because most patients had diagnostic [18F]FDG scans, we also assessed the correlation between [18F]D4-FCH and [18F]FDG SUVs (in the case of [18F]D4-FCH, both early - ~5 min and late ~60 min images from the dynamic scan). Correlations were weak when the two static scans (R2 = 0.1928 and R2 = 0.2478 for SUVmean and SUVmax, respectively) were compared, and only moderate when [18F]FDG SUVs were compared to [18F]D4-FCH SUVs evaluated from the dynamic scan (R2 = 0.4063 and R2 = 0.5080 for SUVmean and SUVmax, respectively).
The deuterium isotope effect involving reduced tunneling of choline substrates to the catabolic (choline oxidase) enzyme site has been proved in this study to stabilize [18F]D4-FCH in systemic circulation, thus permitting substrate transport and retention to be determined at early and late timepoints. By using [18F]D4-FCH PET/CT, we report that substantial heterogeneity for choline substrate transport and retention exists in lung tumors and might confound biopsy measurements of choline pathway protein expression.
This is the first study of [18F]D4-FCH in patients, having previously asserted safety and established dosimetry for human use . We confirmed the slow catabolism of [18F]D4-FCH to [18F]D4-fluorobetaine in patients using radio-HPLC. This is in contrast to [11C]choline and [18F]fluorocholine, which are rapidly catabolized [17, 21] leading to lower proportion of available substrate for tumor uptake for non-deuteriated analogues[35, 36]. Previous studies with choline tracers have attempted to overcome the rapid catabolism by scanning at 5-15 min. We demonstrated an ability to conduct scanning consistently at early (5 min) or late (60 min) timepoints; it is appreciated that differences in perfusion could affect interpretation of early scan protocols. Consequently, in all primary lung lesions uptake with [18F]D4-FCH PET/CT was demonstrated in this study. Using mathematical modelling we were able to interpret the trapping of [18F]D4-FCH as contributed by substrate transport and conversion to phosphocholine.
The availability of routine [18F]FDG PET/CT scans in 12 of the 17 patients allowed us to establish to what extent the two tracers were used by lung tumors. While lesion detection was comparable between the two tracers, uptake intensity (SUV) was overall higher with [18F]FDG (2.87-10.13 vs 6.89-22.64). Studies in lung cancer using [11C]choline PET, have shown variable results. When comparing [11C]choline, typically evaluated at 5-15 min, and [18F]FDG, typically evaluated at 60 min for detection of mediastinal lymph node metastases, Hara et al reported excellent performance (100% sensitivity for both tracers in 29 patients) prior to thoracotomy and ipsilateral node dissection . Other studies have reported a lower diagnostic sensitivity for [11C]choline compared to [18F]FDG PET in detecting nodal disease, but superior sensitivity for detecting of metastatic disease[38, 39]. However, it is important to note that the above studies were PET-only (rather than PET/CT), and in the study by Hara and colleagues, sensitivity may have been overestimated by the low SUV criteria used for lesion detection .
Choline kinase alpha (CHKα) expression in NSCLC. (A) Immunohistochemistry panel showing CHKα expression in lung tissue of patients with NSCLC. Examples are shown of CHKα negative, CHKα positive +1, +2 and +3 intensities on staining. Images are x200 magnification. (B) CHKα expression in lung tissue was found to be negative in four patients (Neg, patients 2, 10,12 and 14), +1 in four patients (patients 4, 5, 6 and 7), +2 in patient 3 and +2/+3 in patient 9. (C) Commercial TMAs (US Biomax Inc.): CHKα was expressed in both primary lesions and lymph nodes in 50% of squamous cell tissue and 63.3% of adenocarcinoma cores. 26.7% and 23.3% of cores expressed CHKα only in the primary squamous cell and adenocarcinoma tissue, respectively. 3.3% and 6.7% of cores expressed CHKα only in the squamous cell and adenocarcinoma lymph nodes, respectively. In both squamous and adenocarcinoma primary lesion and lymph nodes, no expression was found in 20% and 6.7% of cores, respectively.
Parametric analysis of [18F]D4-FCH. (A) Spectral analysis at voxel level resulted with the K1 parametric map. Patlak and Logan methods resulted with 3D parametric maps of Ki and Vt, respectively. (B) Kernel Distribution analysis of K1_SA (on the left) and Ki_SA (on the right) of a representative patient. Tumor, reference tissue (normal lung) and vertebra are shown in this patient as an example. No lymph node metastases were present. On the y axis, the total number of data points per bin. The kernel distributions illustrate the heterogeneity of tumor tissue compared with normal tissue. Normal lung, shown in green, has a more homogenous peaked distribution, as expected. (C) Global kernel distributions of parametric maps (Spectral Analysis K1 on the left and SA_Ki on the right) evaluated in healthy lung (in blue, n = 10), lesions which showed a negative CHKα expression (in red, n = 3) and lesions which showed a positive CHKα expression (in green, n = 7). On the y axis, N is the total number of data points per bin. (D) Kurtosis and Skewness values of Spectral Analysis K1 map in patients with negative (n = 3) and positive (n = 7) CHKα expression. Statistically significant difference was found between the K1_SA kurtosis of CHKα- and CHKα+ lesions (p < 0.05).
The high spatio-molecular distribution of choline metabolites in mouse tumors , and aforementioned resistance phenotype regulating/associated with choline metabolism supports exploitation of this radiotracer for localized therapy planning. Choline tracers have indeed been used to assess target volume, support focal dose escalation decisions and for monitoring response to external beam radiotherapy , particularly in the setting of radical prostate radiotherapy with concurrent androgen deprivation . Thus, we investigated whether similar spatio-molecular distribution of choline metabolites exists and found that radiotracer transfer from blood to tissue, (described by K1) and radiotracer retention (described by Ki; mainly due to labelled phosphocholine formation) were highly heterogeneous in distribution. This provides the opportunity to test in a prospective manner, the role of [18F]D4-FCH PET/CT for planning focal dose escalation to “highest-risk” tumor regions.
Increased CHKα expression has been described in several tumor types [42-44]. The role of CHKα mRNA in malignancy and clinical relevance as a prognostic factor in early-stage lung cancer has been shown in a previous study by Ramirez-Molina et al , where majority of patients were stage 1b and tumors were all surgically excised. They demonstrated that patients with higher CHKα expression (1.91 times higher than that of healthy tissues) had a worse clinical outcome than those with low CHKα expression . We did not find any significant statistical correlation between [18F]D4-FCH and CHKα expression (p = 0.968). While this could be interpreted as a lack of influence of CHKα on substrate retention, our data disputes this hypothesis. a) Mathematical modelling supported a role for substrate trapping due to formation of phosphocholine and not transport, b) Voxel-by-voxel parametric maps indicated existence of high spatio-molecular heterogeneity particularly in 'CHKα negative' suggesting potential biopsy sampling challenges in these tumors might have contributed to this lack of correlation. A larger number of patients would have provided more confidence for this latter hypothesis. Nonetheless external validation in a larger lung cancer tissue bank demonstrated that this phenotype of variable CHKα protein expression in lung cancer was not only pertinent to our prospective study, with only ~50% of cases staining positive for the protein.
Following the study of Boellaard's lab, which employed pharmacokinetic modeling to validate the use of simplified methods for quantification of [18F]FCH , the correlation between SUV and kinetic parameters was also tested in order to establish whether simplified measures could be used instead of a full kinetic modelling. In agreement with Verwer et al. , we obtained poor correlations between both SUV and SUVAUC and Ki_SA, Ki_CM and Vt_L. We obtained higher (but moderate) correlations between mean SUV and SUVAUC evaluated at 60 min and Ki_P and K1_SA. Our results showed that a simpler imaging protocol could not directly replace a full kinetic quantification as SUV does not take into account possible influences by metabolite formation or blood volume. However, the results showed that longer-scan SUV measures give higher correlations (60 min compared to 5 min) and supports SUV60mean compared to SUV60max .
This is an exploratory/ proof of concept study. Based on previous exploratory studies in our group, we estimated that 25 patients would be sufficient to provide initial estimates of the primary and secondary outcome measures. Over the 3-year period in running of this study, the limitations have been multifactorial. Both logistics and difficulties in enrollment have had an impact, with delay in completion of study and reaching target accrual. We recruited 17 patients, and all were included in the final analyses. Subsequently, data analysis was hampered with limitations in diagnostic tissue availability and the lack of [18F]FDG PET/CT in a few patients during staging, whereby no direct comparison with [18F]D4-FCH could be made.
Despite low [18F]D4-FCH uptake in NSCLC patients, there remains the question as to whether PET imaging using choline may still have utility in clinical practice. [18F]FDG PET/CT showed higher SUV than [18F]D4-FCH PET/CT , thus it is unlikely that [18F]D4-FCH will replace [18F]FDG for staging in lung cancer. However, FDG uptake is more homogeneous as it is utilized by all viable cells. We show here that [18F]D4-FCH uptake is heterogeneous, demonstrating potential use in treatment planning (Figure S3). The quantitative measurements provided by [18F]D4-FCH PET/CT uptake in this study may be complementary alongside [18F]FDG PET/CT in primary staging and management. Previous studies have explored the use of metabolic imaging to directly compare dose-escalation plans with standard uniform dose-prescriptions [46-48]. Delivery of an integrated simultaneous boost guided by metabolic imaging may provide the spatial mapping sought in radiotherapy planning.
In conclusion, we have shown late (60 min p.i) imaging with the new choline tracer, [18F]D4-FCH PET/CT designed according to the deuterium isotope effect principles, is feasible and that its intratumorally distribution is highly heterogeneous warranting future studies using this radiotracer for potential radiotherapy dose-delivery and treatment response monitoring.
[18F]D4-FCH: [18F]fluoromethyl-(1,2-2H4) choline; [18F]FCH: [18F]fluoromethylcholine; [18F]FDG: [18F]fluorodeoxyglucose; [18F]FEC: [18F]fluoroethylcholine; 1TCM2k: 1-tissue 2k compartmental model; 2TCM4k: 2-tissue 3k reversible compartmental model; AIC: akaike information criterion; ARSAC: administration of radioactive substances advisory committee; CECT: contrast enhanced computed tomography; CHKα: choline kinase-α; CT: computed tomography; DI-SA: double input spectral analysis; EGFR: epidermal growth factor receptor; FMISO: [18F]fluoromisonidazole; HPLC: high-performance liquid chromatography; IRF: impulse response function; KRAS: kristen rat sarcoma viral oncogene; MHRA: medicine and healthcare products regulatory agency; NSCLC: non-small cell lung cancer; OSEM: ordered-subsets expectation maximization algorithm; PET: positron emission tomography; SA: spectral analysis; SUV: standardized uptake value; TMA: tissue microarrays; VOI: volumes of interest.
Supplementary figures and tables.
The authors would like to thank all patients who selflessly contributed to this research. We would also like to thank Invicro Ltd., Histopathology department at Imperial College Healthcare NHS trust, Tissue bank, Imperial College, London and the West London Cancer Research Network. In addition, the authors would like to thank Dr Louis Allott for providing the image of the chemical structure of [18F]D4-FCH and advice on the manuscript.
Study funded by the UK Medical Research Council award MR/N020782/1. The authors would like to acknowledge support from Imperial Experimental Cancer Medicines Centre and Imperial NIHR Biomedical Research Centre.
Imperial College London is a member of the Cancer Research UK funded National Cancer Imaging Translational Accelerator award (NCITA; C2536/A28680).
All procedures performed in studies involving human participants were in accordance with the ethical standards of the institutional research committee and with the 1964 Helsinki declaration and its later amendments or comparable ethical standards.
The authors have declared that no competing interest exists.
1. Talbot J-N, Fartoux L, Balogova S, Nataf V, Kerrou K, Gutman F. et al. Detection of hepatocellular carcinoma with PET/CT: a prospective comparison of 18F-fluorocholine and 18F-FDG in patients with cirrhosis or chronic liver disease. J Nucl Med. 2010;51:1699-706
2. Challapalli A, Aboagye EO. Positron Emission Tomography Imaging of Tumor Cell Metabolism and Application to Therapy Response Monitoring. Front Oncol. 2016;6:44
3. DeGrado TR, Coleman RE, Wang S, Baldwin SW, Orr MD, Robertson CN. et al. Synthesis and evaluation of 18F-labeled choline as an oncologic tracer for positron emission tomography: initial findings in prostate cancer. Cancer Res. 2001;61:110-7
4. Bansal A, Shuyan W, Hara T, Harris RA, Degrado TR. Biodisposition and metabolism of [(18)F]fluorocholine in 9L glioma cells and 9L glioma-bearing fisher rats. Eur J Nucl Med Mol Imaging. 2008;35:1192-203
5. Roivainen A, Forsback S, Grönroos T, Lehikoinen P, Kähkönen M, Sutinen E. et al. Blood metabolism of [methyl-11 C] choline; implications for in vivo imaging with positron emission tomography. Eur J Nucl Med. 2000;27:25-32
6. Smith G, Zhao Y, Leyton J, Shan B, Perumal M, Turton D. et al. Radiosynthesis and pre-clinical evaluation of [18F] fluoro-[1, 2-2H4] choline. Nucl Med Biol. 2011;38:39-51
7. Verwer EE, Oprea-Lager DE, van den Eertwegh AJ, van Moorselaar RJ, Windhorst AD, Schwarte LA. et al. Quantification of 18F-fluorocholine kinetics in patients with prostate cancer. J Nucl Med. 2015;56:365-71
8. Gadda G. pH and deuterium kinetic isotope effects studies on the oxidation of choline to betaine-aldehyde catalyzed by choline oxidase. Biochim Biophys Acta. 2003;1650:4-9
9. Leyton J, Smith G, Zhao Y, Perumal M, Nguyen Q-D, Robins E. et al. [18F] fluoromethyl-[1, 2-2H4]-choline: a novel radiotracer for imaging choline metabolism in tumors by positron emission tomography. Cancer Res. 2009;69:7721-8
10. Cheng M, Jolly S, Quarshie WO, Kapadia N, Vigneau FD, Kong FS. Modern Radiation Further Improves Survival in Non-Small Cell Lung Cancer: An Analysis of 288,670 Patients. J Cancer. 2019;10:168-77
11. Vera P, Mihailescu S-D, Lequesne J, Modzelewski R, Bohn P, Hapdey S. et al. Radiotherapy boost in patients with hypoxic lesions identified by 18 F-FMISO PET/CT in non-small-cell lung carcinoma: can we expect a better survival outcome without toxicity?[RTEP5 long-term follow-up]. Eur J Nucl Med Mol Imaging. 2019;46:1448-56
12. Watanabe S, Inoue T, Okamoto S, Magota K, Takayanagi A, Sakakibara-Konishi J. et al. Combination of FDG-PET and FMISO-PET as a treatment strategy for patients undergoing early-stage NSCLC stereotactic radiotherapy. EJNMMI Res. 2019;9:104
13. Glunde K, Shah T, Winnard PT, Raman V, Takagi T, Vesuna F. et al. Hypoxia regulates choline kinase expression through hypoxia-inducible factor-1α signaling in a human prostate cancer model. Cancer Res. 2008;68:172-80
14. de Molina AR, Sarmentero-Estrada J, Belda-Iniesta C, Tarón M, de Molina VR, Cejas P. et al. Expression of choline kinase alpha to predict outcome in patients with early-stage non-small-cell lung cancer: a retrospective study. Lancet Oncol. 2007;8:889-97
15. Kroenke M, Hirata K, Gafita A, Watanabe S, Okamoto S, Magota K. et al. Voxel based comparison and texture analysis of 18F-FDG and 18F-FMISO PET of patients with head-and-neck cancer. PLoS One. 2019 14
16. Smith G, Zhao Y, Leyton J, Shan B, Nguyen QD, Perumal M. et al. Radiosynthesis and pre-clinical evaluation of [(18)F]fluoro-[1,2-(2)H(4)]choline. Nucl Med Biol. 2011;38:39-51
17. Challapalli A, Sharma R, Hallett WA, Kozlowski K, Carroll L, Brickute D. et al. Biodistribution and Radiation Dosimetry of Deuterium-Substituted 18F-Fluoromethyl-[1, 2-2H4]Choline in Healthy Volunteers. J Nucl Med. 2014;55:256-63
18. Sharma R, Valls PO, Inglese M, Dubash S, Chen M, Gabra H. et al. [(18)F]Fluciclatide PET as a biomarker of response to combination therapy of pazopanib and paclitaxel in platinum-resistant/refractory ovarian cancer. Eur J Nucl Med Mol Imaging. 2020;47:1239-51
19. Sharma R, Inglese M, Dubash S, Lu H, Pinato DJ, Sanghera C. et al. Monitoring response to transarterial chemoembolization in hepatocellular carcinoma using (18)F-Fluorothymidine Positron Emission Tomography. J Nucl Med. 2020
20. Tomasi G, Kimberley S, Rosso L, Aboagye E, Turkheimer F. Double-input compartmental modeling and spectral analysis for the quantification of positron emission tomography data in oncology. Phys Med Biol. 2012;57:1889
21. Verwer EE, Oprea-Lager DE, van den Eertwegh A, van Moorselaar R, Windhorst AD, Schwarte LA. et al. Quantification of 18F-fluorocholine kinetics in patients with prostate cancer. J Nucl Med. 2015;56:365-71
22. Veronese M, Rizzo G, Turkheimer FE, Bertoldo A. SAKE: a new quantification tool for positron emission tomography studies. Comput Methods Programs Biomed. 2013;111:199-213
23. Bertoldo A, Vicini P, Sambuceti G, Lammertsma AA, Parodi O, Cobelli C. Evaluation of compartmental and spectral analysis models of [18F]FDG kinetics for heart and brain studies with PET. IEEE Trans Biomed Eng. 1998;45:1429-48
24. Turkheimer F, Moresco RM, Lucignani G, Sokoloff L, Fazio F, Schmidt K. The use of spectral analysis to determine regional cerebral glucose utilization with positron emission tomography and [18F]fluorodeoxyglucose: theory, implementation, and optimization procedures. J Cereb Blood Flow Metab. 1994;14:406-22
25. Witney TH, Alam IS, Turton DR, Smith G, Carroll L, Brickute D. et al. Evaluation of deuterated 18F-and 11C-labeled choline analogs for cancer detection by positron emission tomography. Clin Cancer Res. 2012;18:1063-72
26. Patlak CS, Blasberg RG, Fenstermacher JD. Graphical evaluation of blood-to-brain transfer constants from multiple-time uptake data. J Cereb Blood Flow Metab. 1983;3:1-7
27. Logan J, Fowler JS, Volkow ND, Wolf AP, Dewey SL, Schlyer DJ. et al. Graphical analysis of reversible radioligand binding from time-activity measurements applied to [N-11C-methyl]-(-)-cocaine PET studies in human subjects. J Cereb Blood Flow Metab. 1990;10:740-7
28. Fan Z, Calsolaro V, Atkinson RA, Femminella GD, Waldman A, Buckley C. et al. Flutriciclamide (18F-GE180) PET: First-in-Human PET Study of Novel Third-Generation In vivo Marker of Human Translocator Protein. J Nucl Med. 2016;57:1753-9
29. Glatting G, Kletting P, Reske SN, Hohl K, Ring C. Choosing the optimal fit function: comparison of the Akaike information criterion and the F-test. Med Phys. 2007;34:4285-92
30. Turkheimer FE, Hinz R, Cunningham VJ. On the undecidability among kinetic models: from model selection to model averaging. J Cereb Blood Flow Metab. 2003;23:490-8
31. Contractor K, Challapalli A, Barwick T, Winkler M, Hellawell G, Hazell S. et al. Use of [11C]choline PET-CT as a noninvasive method for detecting pelvic lymph node status from prostate cancer and relationship with choline kinase expression. Clin Cancer Res. 2011;17:7673-83
32. Challapalli A, Trousil S, Hazell S, Kozlowski K, Gudi M, Aboagye EO. et al. Exploiting altered patterns of choline kinase-alpha expression on human prostate tissue to prognosticate prostate cancer. J Clin Pathol. 2015;68:703-9
33. Rietbergen DD, van der Hiel B, Vogel W, Stokkel MP. Mediastinal lymph node uptake in patients with prostate carcinoma on F18-choline PET/CT. Nucl Med Commun. 2011;32:1143-7
34. Veronese M, Bertoldo A, Bishu S, Unterman A, Tomasi G, Smith CB. et al. A spectral analysis approach for determination of regional rates of cerebral protein synthesis with the L-[1-11C] leucine PET method. J Cereb Blood Flow Metab. 2010;30:1460-76
35. Hara T, Kosaka N, Shinoura N, Kondo T. PET imaging of brain tumor with [methyl-11C]choline. J Nucl Med. 1997;38:842-7
36. Roivainen A, Forsback S, Gronroos T, Lehikoinen P, Kahkonen M, Sutinen E. et al. Blood metabolism of [methyl-11C]choline; implications for in vivo imaging with positron emission tomography. Eur J Nucl Med. 2000;27:25-32
37. Hara T, Inagaki K, Kosaka N, Morita T. Sensitive detection of mediastinal lymph node metastasis of lung cancer with 11C-choline PET. J Nucl Med. 2000;41:1507-13
38. Khan N, Oriuchi N, Zhang H, Higuchi T, Tian M, Inoue T. et al. A comparative study of 11C-choline PET and [18F]fluorodeoxyglucose PET in the evaluation of lung cancer. Nucl Med Commun. 2003;24:359-66
39. Pieterman RM, Que TH, Elsinga PH, Pruim J, van Putten JW, Willemsen AT. et al. Comparison of (11)C-choline and (18)F-FDG PET in primary diagnosis and staging of patients with thoracic cancer. J Nucl Med. 2002;43:167-72
40. Verma V, Choi JI, Sawant A, Gullapalli RP, Chen W, Alavi A. et al. Use of PET and other functional imaging to guide target delineation in radiation oncology. Seminars in radiation oncology: Elsevier. 2018 p. 171-7
41. Challapalli A, Barwick T, Tomasi G, O'Doherty M, Contractor K, Stewart S. et al. Exploring the potential of [11C] choline-PET/CT as a novel imaging biomarker for predicting early treatment response in prostate cancer. Nucl Med Commun. 2014;35:20-9
42. Ramirez de Molina A, Gutierrez R, Ramos MA, Silva JM, Silva J, Bonilla F. et al. Increased choline kinase activity in human breast carcinomas: clinical evidence for a potential novel antitumor strategy. Oncogene. 2002;21:4317-22
43. Ramirez de Molina A, Rodriguez-Gonzalez A, Gutierrez R, Martinez-Pineiro L, Sanchez J, Bonilla F. et al. Overexpression of choline kinase is a frequent feature in human tumor-derived cell lines and in lung, prostate, and colorectal human cancers. Biochem Biophys Res Commun. 2002;296:580-3
44. Ramirez de Molina A, Sarmentero-Estrada J, Belda-Iniesta C, Taron M, Ramirez de Molina V, Cejas P. et al. Expression of choline kinase alpha to predict outcome in patients with early-stage non-small-cell lung cancer: a retrospective study. Lancet Oncol. 2007;8:889-97
45. Adams MC, Turkington TG, Wilson JM, Wong TZ. A systematic review of the factors affecting accuracy of SUV measurements. American Journal of Roentgenology. 2010;195:310-20
46. Wang X, Cui H, Gong G, Fu Z, Zhou J, Gu J. et al. Computational delineation and quantitative heterogeneity analysis of lung tumor on 18F-FDG PET for radiation dose-escalation. Sci Rep. 2018;8:10649
47. Thureau S, Dubray B, Modzelewski R, Bohn P, Hapdey S, Vincent S. et al. FDG and FMISO PET-guided dose escalation with intensity-modulated radiotherapy in lung cancer. Radiat Oncol. 2018;13:208
48. Yap ML, Sun A, Higgins J, Clarke K, Marshall A, Becker N. et al. Adaptive Dose Escalation using Serial Four-dimensional Positron Emission Tomography/Computed Tomography Scans during Radiotherapy for Locally Advanced Non-small Cell Lung Cancer. Clin Oncol (R Coll Radiol). 2016;28(e):199-205
Corresponding author: Professor Eric O Aboagye, Department of Surgery and Cancer, Imperial College London, Hammersmith Hospital, Du Cane Rd., London W120NN, U.K. E-mail: eric.aboagyeac.uk (020 3313 3759).