ABSTRACT
PURPOSE: K trans $$ {K}^{\mathrm{trans}} $$ has often been proposed as a quantitative imaging biomarker for diagnosis, prognosis, and treatment response assessment for various tumors. None of the many software tools for K trans $$ {K}^{\mathrm{trans}} $$ quantification are standardized. The ISMRM Open Science Initiative for Perfusion Imaging-Dynamic Contrast-Enhanced (OSIPI-DCE) challenge was designed to benchmark methods to better help the efforts to standardize K trans $$ {K}^{\mathrm{trans}} $$ measurement. METHODS: A framework was created to evaluate K trans $$ {K}^{\mathrm{trans}} $$ values produced by DCE-MRI analysis pipelines to enable benchmarking. The perfusion MRI community was invited to apply their pipelines for K trans $$ {K}^{\mathrm{trans}} $$ quantification in glioblastoma from clinical and synthetic patients. Submissions were required to include the entrants' K trans $$ {K}^{\mathrm{trans}} $$ values, the applied software, and a standard operating procedure. These were evaluated using the proposed OSIP I gold $$ \mathrm{OSIP}{\mathrm{I}}_{\mathrm{gold}} $$ score defined with accuracy, repeatability, and reproducibility components. RESULTS: Across the 10 received submissions, the OSIP I gold $$ \mathrm{OSIP}{\mathrm{I}}_{\mathrm{gold}} $$ score ranged from 28% to 78% with a 59% median. The accuracy, repeatability, and reproducibility scores ranged from 0.54 to 0.92, 0.64 to 0.86, and 0.65 to 1.00, respectively (0-1 = lowest-highest). Manual arterial input function selection markedly affected the reproducibility and showed greater variability in K trans $$ {K}^{\mathrm{trans}} $$ analysis than automated methods. Furthermore, provision of a detailed standard operating procedure was critical for higher reproducibility. CONCLUSIONS: This study reports results from the OSIPI-DCE challenge and highlights the high inter-software variability within K trans $$ {K}^{\mathrm{trans}} $$ estimation, providing a framework for ongoing benchmarking against the scores presented. Through this challenge, the participating teams were ranked based on the performance of their software tools in the particular setting of this challenge. In a real-world clinical setting, many of these tools may perform differently with different benchmarking methodology.
Subject(s)
Contrast Media , Magnetic Resonance Imaging , Humans , Reproducibility of Results , Magnetic Resonance Imaging/methods , Software , AlgorithmsABSTRACT
BACKGROUND: Diffusion-weighted imaging (DWI) is commonly used to detect prostate cancer, and a major clinical challenge is differentiating aggressive from indolent disease. PURPOSE: To compare 14 site-specific parametric fitting implementations applied to the same dataset of whole-mount pathologically validated DWI to test the hypothesis that cancer differentiation varies with different fitting algorithms. STUDY TYPE: Prospective. POPULATION: Thirty-three patients prospectively imaged prior to prostatectomy. FIELD STRENGTH/SEQUENCE: 3 T, field-of-view optimized and constrained undistorted single-shot DWI sequence. ASSESSMENT: Datasets, including a noise-free digital reference object (DRO), were distributed to the 14 teams, where locally implemented DWI parameter maps were calculated, including mono-exponential apparent diffusion coefficient (MEADC), kurtosis (K), diffusion kurtosis (DK), bi-exponential diffusion (BID), pseudo-diffusion (BID*), and perfusion fraction (F). The resulting parametric maps were centrally analyzed, where differentiation of benign from cancerous tissue was compared between DWI parameters and the fitting algorithms with a receiver operating characteristic area under the curve (ROC AUC). STATISTICAL TEST: Levene's test, P < 0.05 corrected for multiple comparisons was considered statistically significant. RESULTS: The DRO results indicated minimal discordance between sites. Comparison across sites indicated that K, DK, and MEADC had significantly higher prostate cancer detection capability (AUC range = 0.72-0.76, 0.76-0.81, and 0.76-0.80 respectively) as compared to bi-exponential parameters (BID, BID*, F) which had lower AUC and greater between site variation (AUC range = 0.53-0.80, 0.51-0.81, and 0.52-0.80 respectively). Post-processing parameters also affected the resulting AUC, moving from, for example, 0.75 to 0.87 for MEADC varying cluster size. DATA CONCLUSION: We found that conventional diffusion models had consistent performance at differentiating prostate cancer from benign tissue. Our results also indicated that post-processing decisions on DWI data can affect sensitivity and specificity when applied to radiological-pathological studies in prostate cancer. LEVEL OF EVIDENCE: 1 TECHNICAL EFFICACY: Stage 3.
Subject(s)
Diffusion Magnetic Resonance Imaging , Prostatic Neoplasms , Diffusion Magnetic Resonance Imaging/methods , Humans , Male , Prospective Studies , Prostatic Neoplasms/diagnostic imaging , ROC Curve , Retrospective Studies , Sensitivity and SpecificityABSTRACT
PURPOSE: Current state-of-the-art models for estimating the pharmacokinetic parameters do not account for intervoxel movement of the contrast agent (CA). We introduce an optimal mass transport (OMT) formulation that naturally handles intervoxel CA movement and distinguishes between advective and diffusive flows. METHOD: Ten patients with head and neck squamous cell carcinoma (HNSCC) were enrolled in the study between June 2014 and October 2015 and underwent DCE MRI imaging prior to beginning treatment. The CA tissue concentration information was taken as the input in the data-driven OMT model. The OMT approach was tested on HNSCC DCE data that provides quantitative information for forward flux ( ΦF ) and backward flux ( ΦB ). OMT-derived ΦF was compared with the volume transfer constant for CA, Ktrans , derived from the Extended Tofts Model (ETM). RESULTS: The OMT-derived flows showed a consistent jump in the CA diffusive behavior across the images in accordance with the known CA dynamics. The mean forward flux was 0.0082 ± 0.0091 ( min-1 ) whereas the mean advective component was 0.0052 ± 0.0086 ( min-1 ) in the HNSCC patients. The diffusive percentages in forward and backward flux ranged from 8.67% to 18.76% and 12.76% to 30.36%, respectively. The OMT model accounts for intervoxel CA movement and results show that the forward flux ( ΦF ) is comparable with the ETM-derived Ktrans . CONCLUSIONS: This is a novel data-driven study based on optimal mass transport principles applied to patient DCE imaging to analyze CA flow in HNSCC.
Subject(s)
Carcinoma, Squamous Cell/diagnostic imaging , Contrast Media/pharmacokinetics , Diffusion Magnetic Resonance Imaging , Head and Neck Neoplasms/diagnostic imaging , Carcinoma, Squamous Cell/virology , Gadolinium DTPA/pharmacokinetics , Head and Neck Neoplasms/virology , Humans , Kinetics , Models, Theoretical , Papillomavirus Infections/diagnostic imaging , Reproducibility of Results , Retrospective Studies , Treatment OutcomeABSTRACT
Patients with asymptomatic kidney stones have a high rate of progression to becoming symptomatic kidney stones when followed for several years. Small kidney stones are often found incidentally on imaging when evaluating patients for kidney donation, and there is a concern that after nephrectomy, the donor may become symptomatic and incur damage to the remaining kidney. We reviewed kidney donors at our institution with asymptomatic stones and surveyed them several years after donation to see if the stones became clinically active.
Subject(s)
Incidental Findings , Kidney Calculi/diagnosis , Kidney Transplantation/methods , Living Donors , Nephrectomy , Asymptomatic Diseases , Disease Progression , Female , Follow-Up Studies , Humans , Kidney Calculi/complications , Kidney Transplantation/adverse effects , Male , Nephrectomy/adverse effects , Risk Factors , Time Factors , Tomography, X-Ray Computed , Treatment OutcomeABSTRACT
Unraveling the relationship between molecular signatures in the brain and their functional, architectonic, and anatomic correlates is an important neuroscientific goal. It is still not well understood whether the diversity demonstrated by histological studies in the human brain is reflected in the spatial patterning of whole brain transcriptional profiles. Using genome-wide maps of transcriptional distribution of the human brain by the Allen Brain Institute, we test the hypothesis that gene expression profiles are specific to anatomically described brain regions. In this work, we demonstrate that this is indeed the case by showing that gene similarity clusters appear to respect conventional basal-cortical and caudal-rostral gradients. To fully investigate the causes of this observed spatial clustering, we test a connectionist hypothesis that states that the spatial patterning of gene expression in the brain is simply reflective of the fiber tract connectivity between brain regions. We find that although gene expression and structural connectivity are not determined by each other, they do influence each other with a high statistical significance. This implies that spatial diversity of gene expressions is a result of mainly location-specific features but is influenced by neuronal connectivity, such that like cellular species preferentially connects with like cells.
Subject(s)
Brain/anatomy & histology , Brain/physiology , Gene Expression , Atlases as Topic , Brain Mapping , Databases, Genetic , Diffusion Tensor Imaging , Humans , Magnetic Resonance Imaging , Microarray Analysis , Models, Genetic , Models, Neurological , Neural Pathways/anatomy & histology , Neural Pathways/physiologyABSTRACT
There is a need to develop user-friendly imaging tools estimating robust quantitative biomarkers (QIBs) from multiparametric (mp)MRI for clinical applications in oncology. Quantitative metrics derived from (mp)MRI can monitor and predict early responses to treatment, often prior to anatomical changes. We have developed a vendor-agnostic, flexible, and user-friendly MATLAB-based toolkit, MRI-Quantitative Analysis and Multiparametric Evaluation Routines ("MRI-QAMPER", current release v3.0), for the estimation of quantitative metrics from dynamic contrast-enhanced (DCE) and multi-b value diffusion-weighted (DW) MR and MR relaxometry. MRI-QAMPER's functionality includes generating numerical parametric maps from these methods reflecting tumor permeability, cellularity, and tissue morphology. MRI-QAMPER routines were validated using digital reference objects (DROs) for DCE and DW MRI, serving as initial approval stages in the National Cancer Institute Quantitative Imaging Network (NCI/QIN) software benchmark. MRI-QAMPER has participated in DCE and DW MRI Collaborative Challenge Projects (CCPs), which are key technical stages in the NCI/QIN benchmark. In a DCE CCP, QAMPER presented the best repeatability coefficient (RC = 0.56) across test-retest brain metastasis data, out of ten participating DCE software packages. In a DW CCP, QAMPER ranked among the top five (out of fourteen) tools with the highest area under the curve (AUC) for prostate cancer detection. This platform can seamlessly process mpMRI data from brain, head and neck, thyroid, prostate, pancreas, and bladder cancer. MRI-QAMPER prospectively analyzes dose de-escalation trial data for oropharyngeal cancer, which has earned it advanced NCI/QIN approval for expanded usage and applications in wider clinical trials.
Subject(s)
Multiparametric Magnetic Resonance Imaging , Prostatic Neoplasms , Male , Humans , Contrast Media , Prostatic Neoplasms/diagnostic imaging , Prostatic Neoplasms/pathology , Medical Oncology , BiomarkersABSTRACT
We describe a dataset from patients who received ablative radiation therapy for locally advanced pancreatic cancer (LAPC), consisting of computed tomography (CT) and cone-beam CT (CBCT) images with physician-drawn organ-at-risk (OAR) contours. The image datasets (one CT for treatment planning and two CBCT scans at the time of treatment per patient) were collected from 40 patients. All scans were acquired with the patient in the treatment position and in a deep inspiration breath-hold state. Six radiation oncologists delineated the gastrointestinal OARs consisting of small bowel, stomach and duodenum, such that the same physician delineated all image sets belonging to the same patient. Two trained medical physicists further edited the contours to ensure adherence to delineation guidelines. The image and contour files are available in DICOM format and are publicly available from The Cancer Imaging Archive ( https://doi.org/10.7937/TCIA.ESHQ-4D90 , Version 2). The dataset can serve as a criterion standard for evaluating the accuracy and reliability of deformable image registration and auto-segmentation algorithms, as well as a training set for deep-learning-based methods.
Subject(s)
Pancreatic Neoplasms , Radiotherapy Planning, Computer-Assisted , Humans , Cone-Beam Computed Tomography/methods , Image Processing, Computer-Assisted/methods , Pancreatic Neoplasms/diagnostic imaging , Pancreatic Neoplasms/radiotherapy , Radiotherapy Planning, Computer-Assisted/methods , Reproducibility of Results , Tomography, X-Ray ComputedABSTRACT
Objective.Delineating swallowing and chewing structures aids in radiotherapy (RT) treatment planning to limit dysphagia, trismus, and speech dysfunction. We aim to develop an accurate and efficient method to automate this process.Approach.CT scans of 242 head and neck (H&N) cancer patients acquired from 2004 to 2009 at our institution were used to develop auto-segmentation models for the masseters, medial pterygoids, larynx, and pharyngeal constrictor muscle using DeepLabV3+. A cascaded framework was used, wherein models were trained sequentially to spatially constrain each structure group based on prior segmentations. Additionally, an ensemble of models, combining contextual information from axial, coronal, and sagittal views was used to improve segmentation accuracy. Prospective evaluation was conducted by measuring the amount of manual editing required in 91 H&N CT scans acquired February-May 2021.Main results. Medians and inter-quartile ranges of Dice similarity coefficients (DSC) computed on the retrospective testing set (N = 24) were 0.87 (0.85-0.89) for the masseters, 0.80 (0.79-0.81) for the medial pterygoids, 0.81 (0.79-0.84) for the larynx, and 0.69 (0.67-0.71) for the constrictor. Auto-segmentations, when compared to two sets of manual segmentations in 10 randomly selected scans, showed better agreement (DSC) with each observer than inter-observer DSC. Prospective analysis showed most manual modifications needed for clinical use were minor, suggesting auto-contouring could increase clinical efficiency. Trained segmentation models are available for research use upon request viahttps://github.com/cerr/CERR/wiki/Auto-Segmentation-models.Significance.We developed deep learning-based auto-segmentation models for swallowing and chewing structures in CT and demonstrated its potential for use in treatment planning to limit complications post-RT. To the best of our knowledge, this is the only prospectively-validated deep learning-based model for segmenting chewing and swallowing structures in CT. Segmentation models have been made open-source to facilitate reproducibility and multi-institutional research.
Subject(s)
Deep Learning , Head and Neck Neoplasms , Deglutition , Humans , Mastication , Organs at Risk , Radiotherapy Planning, Computer-Assisted/methods , Reproducibility of Results , Retrospective Studies , Tomography, X-Ray Computed/methodsABSTRACT
PURPOSE: Contouring inconsistencies are known but understudied in clinical radiation therapy trials. We applied auto-contouring to the Radiation Therapy Oncology Group (RTOG) 0617 dose escalation trial data. We hypothesized that the trial heart doses were higher than reported due to inconsistent and insufficient heart segmentation. We tested our hypothesis by comparing doses between deep-learning (DL) segmented hearts and trial hearts. METHODS AND MATERIALS: The RTOG 0617 data were downloaded from The Cancer Imaging Archive; the 442 patients with trial hearts and dose distributions were included. All hearts were resegmented using our DL pipeline and quality assured to meet the requirements for clinical implementation. Dose (V5%, V30%, and mean heart dose) was compared between the 2 sets of hearts (Wilcoxon signed-rank test). Each dose metric was associated with overall survival (Cox proportional hazards). Lastly, 18 volume similarity metrics were assessed for the hearts and correlated with |DoseDL - DoseRTOG0617| (linear regression; significance: P ≤ .0028; corrected for 18 tests). RESULTS: Dose metrics were significantly higher for DL hearts compared with trial hearts (eg, mean heart dose: 15 Gy vs 12 Gy; P = 5.8E-16). All 3 DL heart dose metrics were stronger overall survival predictors than those of the trial hearts (median, P = 2.8E-5 vs 2.0E-4). Thirteen similarity metrics explained |DoseDL - DoseRTOG0617|; the axial distance between the 2 centers of mass was the strongest predictor (CENTAxial; median, R2 = 0.47; P = 6.1E-62). CENTAxial agreed with the qualitatively identified inconsistencies in the superior direction. The trial's qualitative heart contouring score was not correlated with |DoseDL - DoseRTOG0617| (median, R2 = 0.01; P = .02) or with any of the similarity metrics (median, Rs = 0.13 [range, -0.22 to 0.31]). CONCLUSIONS: Using a coherent heart definition, as enabled through our open-source DL algorithm, the trial heart doses in RTOG 0617 were found to be significantly higher than previously reported, which may have led to an even more rapid mortality accumulation. Auto-segmentation is likely to reduce contouring and dose inconsistencies and increase the quality of clinical RT trials.
Subject(s)
Clinical Trials, Phase III as Topic , Deep Learning , Heart/diagnostic imaging , Heart/radiation effects , Organs at Risk/diagnostic imaging , Organs at Risk/radiation effects , Adult , Aged , Aged, 80 and over , Algorithms , Carcinoma, Non-Small-Cell Lung/therapy , Confidence Intervals , Female , Humans , Linear Models , Lung Neoplasms/therapy , Male , Middle Aged , Proportional Hazards Models , Radiotherapy Dosage , Statistics, Nonparametric , Tomography, X-Ray ComputedABSTRACT
The present study aims to monitor longitudinal changes in simulated tumor interstitial fluid pressure (IFP) and velocity (IFV) values using dynamic contrast-enhanced (DCE)-MRI-based computational fluid modeling (CFM) in pancreatic ductal adenocarcinoma (PDAC) patients. Nine PDAC patients underwent MRI, including DCE-MRI, on a 3-Tesla MRI scanner at pre-treatment (TX (0)), after the first fraction of stereotactic body radiotherapy (SBRT, (D1-TX)), and six weeks post-TX (D2-TX). The partial differential equation of IFP formulated from the continuity equation, incorporating the Starling Principle of fluid exchange, Darcy velocity, and volume transfer constant (Ktrans), was solved in COMSOL Multiphysics software to generate IFP and IFV maps. Tumor volume (Vt), Ktrans, IFP, and IFV values were compared (Wilcoxon and Spearman) between the time- points. D2-TX Ktrans values were significantly different from pre-TX and D1-TX (p < 0.05). The D1-TX and pre-TX mean IFV values exhibited a borderline significant difference (p = 0.08). The IFP values varying <3.0% between the three time-points were not significantly different (p > 0.05). Vt and IFP values were strongly positively correlated at pre-TX (ρ = 0.90, p = 0.005), while IFV exhibited a strong negative correlation at D1-TX (ρ = -0.74, p = 0.045). Vt, Ktrans, IFP, and IFV hold promise as imaging biomarkers of early response to therapy in PDAC.
ABSTRACT
BACKGROUND AND PURPOSE: Reducing trismus in radiotherapy for head and neck cancer (HNC) is important. Automated deep learning (DL) segmentation and automated planning was used to introduce new and rarely segmented masticatory structures to study if trismus risk could be decreased. MATERIALS AND METHODS: Auto-segmentation was based on purpose-built DL, and automated planning used our in-house system, ECHO. Treatment plans for ten HNC patients, treated with 2 Gy × 35 fractions, were optimized (ECHO0). Six manually segmented OARs were replaced with DL auto-segmentations and the plans re-optimized (ECHO1). In a third set of plans, mean doses for auto-segmented ipsilateral masseter and medial pterygoid (MIMean, MPIMean), derived from a trismus risk model, were implemented as dose-volume objectives (ECHO2). Clinical dose-volume criteria were compared between the two scenarios (ECHO0 vs. ECHO1; ECHO1 vs. ECHO2; Wilcoxon signed-rank test; significance: p < 0.01). RESULTS: Small systematic differences were observed between the doses to the six auto-segmented OARs and their manual counterparts (median: ECHO1 = 6.2 (range: 0.4, 21) Gy vs. ECHO0 = 6.6 (range: 0.3, 22) Gy; p = 0.007), and the ECHO1 plans provided improved normal tissue sparing across a larger dose-volume range. Only in the ECHO2 plans, all patients fulfilled both MIMean and MPIMean criteria. The population median MIMean and MPIMean were considerably lower than those suggested by the trismus model (ECHO0: MIMean = 13 Gy vs. ≤42 Gy; MPIMean = 29 Gy vs. ≤68 Gy). CONCLUSIONS: Automated treatment planning can efficiently incorporate new structures from DL auto-segmentation, which results in trismus risk sparing without deteriorating treatment plan quality. Auto-planning and deep learning auto-segmentation together provide a powerful platform to further improve treatment planning.
ABSTRACT
BACKGROUND AND PURPOSE: To determine the ability of diffusion-weighted imaging (DWI) and dynamic contrast-enhanced magnetic resonance imaging (DCE-MRI) to predict long-term response of brain metastases prior to and within 72 hours of stereotactic radiosurgery (SRS). METHODS: In this prospective pilot study, multiple b-value DWI and T1-weighted DCE-MRI were performed in patients with brain metastases before and within 72 hours following SRS. Diffusion-weighted images were analyzed using the monoexponential and intravoxel incoherent motion (IVIM) models. DCE-MRI data were analyzed using the extended Tofts pharmacokinetic model. The parameters obtained with these methods were correlated with brain metastasis outcomes according to modified Response Assessment in Neuro-Oncology Brain Metastases criteria. RESULTS: We included 25 lesions from 16 patients; 16 patients underwent pre-SRS MRI and 12 of 16 patients underwent both pre- and early (within 72 hours) post-SRS MRI. The perfusion fraction (f) derived from IVIM early post-SRS was higher in lesions demonstrating progressive disease than in lesions demonstrating stable disease, partial response, or complete response (q = .041). Pre-SRS extracellular extravascular volume fraction, ve , and volume transfer coefficient, Ktrans , derived from DCE-MRI were higher in nonresponders versus responders (q = .041). CONCLUSIONS: Quantitative DWI and DCE-MRI are feasible imaging methods in the pre- and early (within 72 hours) post-SRS evaluation of brain metastases. DWI- and DCE-MRI-derived parameters demonstrated physiologic changes (tumor cellularity and vascularity) and offer potentially useful biomarkers that can predict treatment response. This allows for initiation of alternate therapies within an effective time window that may help prevent disease progression.
Subject(s)
Brain Neoplasms/diagnostic imaging , Brain Neoplasms/radiotherapy , Diffusion Magnetic Resonance Imaging , Radiosurgery , Adult , Aged , Contrast Media , Disease Progression , Humans , Male , Middle Aged , Pilot Projects , Prospective StudiesABSTRACT
We aimed to assess longitudinal changes in quantitative imaging metric values obtained from diffusion-weighted (DW-) and dynamic contrast-enhanced magnetic resonance imaging (DCE)-MRI at pre-treatment (TX[0]), immediately after the first fraction of stereotactic body radiotherapy (D1-TX[1]), and 6 weeks post-TX (Post-TX[2]) in patients with pancreatic ductal adenocarcinoma. Ten enrolled patients (n = 10) underwent DW- and DCE-MRI examinations on a 3.0 T scanner. The apparent diffusion coefficient, ADC (mm2/s), was derived from DW imaging data using a monoexponential model. The tissue relaxation rate, R1t, time-course data were fitted with a shutter-speed model, which provides estimates of the volume transfer constant, Ktrans (min-1), extravascular extracellular volume fraction, ve , and mean lifetime of intracellular water protons, τ i (seconds). Wilcoxon rank-sum test compared the mean values, standard deviation, skewness, kurtosis, and relative percentage (r, %) changes (Δ) in ADC, Ktrans, ve , and τ i values between the magnetic resonance examinations. rADCΔ2-0 values were significantly greater than rADCΔ1-0 values (P = .009). rKtransΔ2-0 values were significantly lower than rKtransΔ1-0 values (P = .048). rveΔ2-1 and rveΔ2-0 values were significantly different (P = .016). rτ iΔ2-1 values were significantly lower than rτ iΔ2-0 values (P = .008). For group comparison, the pre-TX mean and kurtosis of ADC (P = .18 and P = .14), skewness and kurtosis of Ktrans values (P = .14 for both) showed a leaning toward significant difference between patients who experienced local control (n = 2) and failed early (n = 4). DW- and DCE-MRI-derived quantitative metrics could be useful biomarkers to evaluate longitudinal changes to stereotactic body radiotherapy in patients with pancreatic ductal adenocarcinoma.
Subject(s)
Magnetic Resonance Imaging , Pancreatic Neoplasms , Radiosurgery , Aged , Antineoplastic Combined Chemotherapy Protocols , Benchmarking , Contrast Media , Female , Humans , Male , Middle Aged , Pancreatic Neoplasms/diagnostic imaging , Pancreatic Neoplasms/radiotherapy , Pancreatic Neoplasms/surgeryABSTRACT
Background: Early imaging-based treatment response assessment of brain metastases following stereotactic radiosurgery (SRS) remains challenging. The aim of this study is to determine whether early (within 12 weeks) intratumoral changes in interstitial fluid pressure (IFP) and velocity (IFV) estimated from computational fluid modeling (CFM) using dynamic contrast-enhanced (DCE) MRI can predict long-term outcomes of lung cancer brain metastases (LCBMs) treated with SRS. Methods: Pre- and post-treatment T1-weighted DCE-MRI data were obtained in 41 patients treated with SRS for intact LCBMs. The imaging response was assessed using RANO-BM criteria. For each lesion, extravasation of contrast agent measured from Extended Tofts pharmacokinetic Model (volume transfer constant, Ktrans) was incorporated into a computational fluid model to estimate tumor IFP and IFV. Estimates of mean IFP and IFV and heterogeneity (skewness and kurtosis) were calculated for each lesion from pre- and post-SRS imaging. The Wilcoxon rank-sum test was utilized to assess for significant differences in IFP, IFV, and IFP/IFV change (Δ) between response groups. Results: Fifty-three lesions from 41 patients were included. Median follow-up time after SRS was 11 months. The objective response (OR) rate (partial or complete response) was 79%, with 21% demonstrating stable disease (SD) or progressive disease (PD). There were significant response group differences for multiple posttreatment and Δ CFM parameters: post-SRS IFP skewness (mean -0.405 vs. -0.691, p = 0.022), IFP kurtosis (mean 2.88 vs. 3.51, p = 0.024), and IFV mean (5.75e-09 vs. 4.19e-09 m/s, p = 0.027); and Δ IFP kurtosis (mean -2.26 vs. -0.0156, p = 0.017) and IFV mean (1.91e-09 vs. 2.38e-10 m/s, p = 0.013). Posttreatment and Δ thresholds predicted non-OR with high sensitivity (sens): post-SRS IFP skewness (-0.432, sens 84%), kurtosis (2.89, sens 84%), and IFV mean (4.93e-09 m/s, sens 79%); and Δ IFP kurtosis (-0.469, sens 74%) and IFV mean (9.90e-10 m/s, sens 74%). Conclusions: Objective response was associated with lower post-treatment tumor heterogeneity, as represented by reductions in IFP skewness and kurtosis. These results suggest that early post-treatment assessment of IFP and IFV can be used to predict long-term response of lung cancer brain metastases to SRS, allowing a timelier treatment modification.
ABSTRACT
We developed and tested the feasibility of computational fluid modeling (CFM) based on dynamic contrast-enhanced magnetic resonance imaging (DCE-MRI) for quantitative estimation of interstitial fluid pressure (IFP) and velocity (IFV) in patients with head and neck (HN) cancer with locoregional lymph node metastases. Twenty-two patients with HN cancer, with 38 lymph nodes, underwent pretreatment standard MRI, including DCE-MRI, on a 3-Tesla scanner. CFM simulation was performed with the finite element method in COMSOL Multiphysics software. The model consisted of a partial differential equation (PDE) module to generate 3D parametric IFP and IFV maps, using the Darcy equation and Ktrans values (min-1, estimated from the extended Tofts model) to reflect fluid influx into tissue from the capillary microvasculature. The Spearman correlation (ρ) was calculated between total tumor volumes and CFM estimates of mean tumor IFP and IFV. CFM-estimated tumor IFP and IFV mean ± standard deviation for the neck nodal metastases were 1.73 ± 0.39 (kPa) and 1.82 ± 0.9 × (10-7 m/s), respectively. High IFP estimates corresponds to very low IFV throughout the tumor core, but IFV rises rapidly near the tumor boundary where the drop in IFP is precipitous. A significant correlation was found between pretreatment total tumor volume and CFM estimates of mean tumor IFP (ρ = 0.50, P = 0.004). Future studies can validate these initial findings in larger patients with HN cancer cohorts using CFM of the tumor in concert with DCE characterization, which holds promise in radiation oncology and drug-therapy clinical trials.
Subject(s)
Extracellular Fluid , Head and Neck Neoplasms , Magnetic Resonance Spectroscopy , Contrast Media , Extracellular Fluid/physiology , Feasibility Studies , Female , Head and Neck Neoplasms/diagnostic imaging , Head and Neck Neoplasms/physiopathology , Humans , Male , PressureABSTRACT
An open-source library of implementations for deep-learning-based image segmentation and outcomes models based on radiotherapy and radiomics is presented. As oncology treatment planning becomes increasingly driven by automation, such a library of model implementations is crucial to (i) validate existing models on datasets collected at different institutions, (ii) automate segmentation, (iii) create ensembles for improving performance and (iv) incorporate validated models in the clinical workflow. Inclusion of deep-learning-based image segmentation and outcomes models in the same library provides a fully automated and reproduceable pipeline to estimate prognosis. The library was developed with the Computational Environment for Radiological Research (CERR) software platform. Centralizing model implementations in CERR builds upon its rich set of radiotherapy and radiomics tools and caters to the world-wide user base. CERR provides well-validated feature extraction pipelines for radiotherapy dosimetry and radiomics with fine control over the calculation settings, allowing users to select appropriate parameters used in model derivation. Models for automatic image segmentation are distributed via containers, allowing them to be deployed with a variety of scientific computing architectures. The library includes implementations of popular DVH-based models outlined in the Quantitative Analysis of Normal Tissue Effects in the Clinic effort and recently published literature. Radiomics models include features from the Image Biomarker Standardization Initiative and application-specific features found to be relevant across multiple sites and image modalities. The library is distributed as a module within CERR at https://www.github.com/cerr/CERR under the GNU-GPL copyleft with additional restrictions on clinical and commercial use and provision to dual license in future.
Subject(s)
Deep Learning , Image Processing, Computer-Assisted/methods , Reproducibility of ResultsABSTRACT
Quantitative kurtosis phantoms are sought by multicenter clinical trials to establish accuracy and precision of quantitative imaging biomarkers on the basis of diffusion kurtosis imaging (DKI) parameters. We designed and evaluated precision, reproducibility, and long-term stability of a novel isotropic (i)DKI phantom fabricated using four families of chemicals based on vesicular and lamellar mesophases of liquid crystal materials. The constructed iDKI phantoms included negative control monoexponential diffusion materials to independently characterize noise and model-induced bias in quantitative kurtosis parameters. Ten test-retest DKI studies were performed on four scanners at three imaging centers over a six-month period. The tested prototype phantoms exhibited physiologically relevant apparent diffusion, Dapp, and kurtosis, Kapp, parameters ranging between 0.4 and 1.1 (×10-3 mm2/s) and 0.8 and 1.7 (unitless), respectively. Measured kurtosis phantom Kapp exceeded maximum fit model bias (0.1) detected for negative control (zero kurtosis) materials. The material-specific parameter precision [95% CI for Dapp: 0.013-0.022(×10-3 mm2/s) and for Kapp: 0.009-0.076] derived from the test-retest analysis was sufficient to characterize thermal and temporal stability of the prototype DKI phantom through correlation analysis of inter-scan variability. The present study confirms a promising chemical design for stable quantitative DKI phantom based on vesicular mesophase of liquid crystal materials. Improvements to phantom preparation and temperature monitoring procedures have potential to enhance precision and reproducibility for future multicenter iDKI phantom studies.
Subject(s)
Diffusion Magnetic Resonance Imaging/standards , Phantoms, Imaging/standards , Diffusion Magnetic Resonance Imaging/methods , Equipment Design , Humans , Image Interpretation, Computer-Assisted/methods , Liquid Crystals , Reproducibility of Results , TemperatureABSTRACT
The aim of this study was to establish the repeatability measures of quantitative Gaussian and non-Gaussian diffusion metrics using diffusion-weighted imaging (DWI) data from phantoms and patients with head-and-neck and papillary thyroid cancers. The Quantitative Imaging Biomarker Alliance (QIBA) DWI phantom and a novel isotropic diffusion kurtosis imaging phantom were scanned at 3 different sites, on 1.5T and 3T magnetic resonance imaging systems, using standardized multiple b-value DWI acquisition protocol. In the clinical component of this study, a total of 60 multiple b-value DWI data sets were analyzed for test-retest, obtained from 14 patients (9 head-and-neck squamous cell carcinoma and 5 papillary thyroid cancers). Repeatability of quantitative DWI measurements was assessed by within-subject coefficient of variation (wCV%) and Bland-Altman analysis. In isotropic diffusion kurtosis imaging phantom vial with 2% ceteryl alcohol and behentrimonium chloride solution, the mean apparent diffusion (Dapp × 10-3 mm2/s) and kurtosis (Kapp, unitless) coefficient values were 1.02 and 1.68 respectively, capturing in vivo tumor cellularity and tissue microstructure. For the same vial, Dapp and Kapp mean wCVs (%) were ≤1.41% and ≤0.43% for 1.5T and 3T across 3 sites. For pretreatment head-and-neck squamous cell carcinoma, apparent diffusion coefficient, D, D*, K, and f mean wCVs (%) were 2.38%, 3.55%, 3.88%, 8.0%, and 9.92%, respectively; wCVs exhibited a higher trend for papillary thyroid cancers. Knowledge of technical precision and bias of quantitative imaging metrics enables investigators to properly design and power clinical trials and better discern between measurement variability versus biological change.
Subject(s)
Diffusion Magnetic Resonance Imaging/standards , Head and Neck Neoplasms/diagnostic imaging , Phantoms, Imaging , Squamous Cell Carcinoma of Head and Neck/diagnostic imaging , Thyroid Neoplasms/diagnostic imaging , Adult , Aged , Diffusion Magnetic Resonance Imaging/methods , Female , Humans , Image Interpretation, Computer-Assisted/methods , Male , Middle Aged , Reproducibility of ResultsABSTRACT
BACKGROUND: Platelet-activating factor acetylhydrolase 1B1 (LIS1), a critical mediator of neuronal migration in developing brain, is expressed throughout life. However, relatively little is known about LIS1 function in the mature brain. We previously demonstrated that LIS1 involvement in the formation and turnover of synaptic protrusions and synapses of young brain after neuronal migration is complete. Here we examine the requirement for LIS1 to maintain hippocampal circuit function in adulthood. METHODS: Effects of conditional Lis1 inactivation in excitatory pyramidal neurons, starting in juvenile mouse brain, were probed using high-resolution approaches combining mouse genetics, designer receptor exclusively activated by designer drug technology to specifically manipulate CA1 pyramidal neuron excitatory activity, electrophysiology, hippocampus-selective behavioral testing, and magnetic resonance imaging tractography to examine the connectivity of LIS1-deficient neurons. RESULTS: We found progressive excitatory and inhibitory postsynaptic dysfunction as soon as 10 days after conditional inactivation of Lis1 targeting CA1 pyramidal neurons. Surprisingly, by postnatal day 60 it also caused CA1 histological disorganization, with a selective decline in parvalbumin-expressing interneurons and further reduction in inhibitory neurotransmission. Accompanying these changes were behavioral and cognitive deficits that could be rescued by either designer receptor exclusively activated by designer drug-directed specific increases in CA1 excitatory transmission or pharmacological enhancement of gamma-aminobutyric acid transmission. Lagging behind electrophysiological changes was a progressive, selective decline in neural connectivity, affecting hippocampal efferent pathways documented by magnetic resonance imaging tractography. CONCLUSIONS: LIS1 supports synaptic function and plasticity of mature CA1 neurons. Postjuvenile loss of LIS1 disrupts the structure and cellular composition of the hippocampus, its connectivity with other brain regions, and cognition dependent on hippocampal circuits.
Subject(s)
1-Alkyl-2-acetylglycerophosphocholine Esterase/metabolism , Cognition/physiology , Hippocampus/cytology , Microtubule-Associated Proteins/metabolism , Neurons/physiology , Synapses/physiology , 1-Alkyl-2-acetylglycerophosphocholine Esterase/genetics , Animals , Animals, Newborn , Calcium-Calmodulin-Dependent Protein Kinase Type 2/genetics , Calcium-Calmodulin-Dependent Protein Kinase Type 2/metabolism , Cell Movement/genetics , Clonazepam/pharmacology , Cognition/drug effects , Conditioning, Psychological/physiology , Excitatory Postsynaptic Potentials/drug effects , Excitatory Postsynaptic Potentials/genetics , Fear/physiology , GABA Modulators/pharmacology , Hippocampus/diagnostic imaging , Locomotion/genetics , Mice , Mice, Inbred C57BL , Mice, Transgenic , Microtubule-Associated Proteins/genetics , Nerve Tissue Proteins/metabolism , Neuronal Plasticity/drug effects , Neuronal Plasticity/genetics , Neurons/drug effects , Receptors, G-Protein-Coupled/genetics , Receptors, G-Protein-Coupled/metabolism , Recognition, Psychology/physiology , Synapses/drug effectsABSTRACT
Spatiotemporal tau pathology progression is regarded as highly stereotyped within each type of degenerative condition. For instance, AD has a progression of tau pathology consistently beginning in the entorhinal cortex, the locus coeruleus, and other nearby noradrenergic brainstem nuclei, before spreading to the rest of the limbic system as well as the cingulate and retrosplenial cortices. Proposed explanations for the consistent spatial patterns of tau pathology progression, as well as for why certain regions are selectively vulnerable to exhibiting pathology over the course of disease generally focus on transsynaptic spread proceeding via the brain's anatomic connectivity network in a cell-independent manner or on cell-intrinsic properties that might render some cell populations or regions uniquely vulnerable. We test connectivity based explanations of spatiotemporal tau pathology progression and regional vulnerability against cell-intrinsic explanation, using regional gene expression profiles as a proxy. We find that across both exogenously seeded and non-seeded tauopathic mouse models, the connectivity network provides a better explanation than regional gene expression profiles, even when such profiles are limited to specific sets of tau risk-related genes only. Our results suggest that, regardless of the location of pathology initiation, tau pathology progression is well characterized by a model positing entirely cell-type and molecular environment independent transsynaptic spread via the mouse brain's connectivity network. These results further suggest that regional vulnerability to tau pathology is mainly governed by connectivity with regions already exhibiting pathology, rather than by cell-intrinsic factors.