Independent component analysis (ICA) applied to dynamic oxygen-enhanced MRI (OE-MRI) for robust functional lung imaging at 3 T.

PURPOSE: Dynamic lung oxygen-enhanced MRI (OE-MRI) is challenging due to the presence of confounding signals and poor signal-to-noise ratio, particularly at 3 T. We have created a robust pipeline utilizing independent component analysis (ICA) to automatically extract the oxygen-induced signal change from confounding factors to improve the accuracy and sensitivity of lung OE-MRI. METHODS: Dynamic OE-MRI was performed on healthy participants using a dual-echo multi-slice spoiled gradient echo sequence at 3 T and cyclical gas delivery. ICA was applied to each echo within a thoracic mask. The ICA component relating to the oxygen-enhancement signal was automatically identified using correlation analysis. The oxygen-enhancement component was reconstructed, and the percentage signal enhancement (PSE) was calculated. The lung PSE of current smokers was compared with nonsmokers; scan-rescan repeatability, ICA pipeline repeatability, and reproducibility between two vendors were assessed. RESULTS: ICA successfully extracted a consistent oxygen-enhancement component for all participants. Lung tissue and oxygenated blood displayed the opposite oxygen-induced signal enhancements. A significant difference in PSE was observed between the lungs of current smokers and nonsmokers. The scan-rescan repeatability and the ICA pipeline repeatability were good. CONCLUSION: The developed pipeline demonstrated sensitivity to the signal enhancements of the lung tissue and oxygenated blood at 3 T. The difference in lung PSE between current smokers and nonsmokers indicates a likely sensitivity to lung function alterations that may be seen in mild pathology, supporting future use of our methods in patient studies.

Feasibility of dynamic T2 *-based oxygen-enhanced lung MRI at 3T.

PURPOSE: To demonstrate proof-of-concept of a T2 *-sensitized oxygen-enhanced MRI (OE-MRI) method at 3T by assessing signal characteristics, repeatability, and reproducibility of dynamic lung OE-MRI metrics in healthy volunteers. METHODS: We performed sequence-specific simulations for protocol optimisation and acquired free-breathing OE-MRI data from 16 healthy subjects using a dual-echo RF-spoiled gradient echo approach at 3T across two institutions. Non-linear registration and tissue density correction were applied. Derived metrics included percent signal enhancement (PSE), ∆R2 * and wash-in time normalized for breathing rate (τ-nBR). Inter-scanner reproducibility and intra-scanner repeatability were evaluated using intra-class correlation coefficient (ICC), repeatability coefficient, reproducibility coefficient, and Bland-Altman analysis. RESULTS: Simulations and experimental data show negative contrast upon oxygen inhalation, due to substantial dominance of ∆R2 * at TE > 0.2 ms. Density correction improved signal fluctuations. Density-corrected mean PSE values, aligned with simulations, display TE-dependence, and an anterior-to-posterior PSE reduction trend at TE1 . ∆R2 * maps exhibit spatial heterogeneity in oxygen delivery, featuring anterior-to-posterior R2 * increase. Mean T2 * values across 32 scans were 0.68 and 0.62 ms for pre- and post-O2 inhalation, respectively. Excellent or good agreement emerged from all intra-, inter-scanner and inter-rater variability tests for PSE and ∆R2 *. However, ICC values for τ-nBR demonstrated limited agreement between repeated measures. CONCLUSION: Our results demonstrate the feasibility of a T2 *-weighted method utilizing a dual-echo RF-spoiled gradient echo approach, simultaneously capturing PSE, ∆R2 * changes, and oxygen wash-in during free-breathing. The excellent or good repeatability and reproducibility on intra- and inter-scanner PSE and ∆R2 * suggest potential utility in multi-center clinical applications.

The ISMRM Open Science Initiative for Perfusion Imaging (OSIPI): Results from the OSIPI-Dynamic Contrast-Enhanced challenge.

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.

Ultrafast Charge Transfer Dynamics in a Slip-Stacked Donor-Acceptor-Acceptor System.

Photoexcitation of molecular electron donor and/or acceptor chromophore aggregates can greatly affect their charge-transfer dynamics. Excitonic coupling not only alters the energy landscape in the excited state but may also open new photophysical pathways, such as symmetry-breaking charge separation (SB-CS). Here, we investigate the impact of excitonic coupling on a covalent donor-acceptor-acceptor system comprising a perylene donor (Per) and two perylenediimide (PDI) acceptor chromophores in which the three components are π-stacked in a geometry that is slipped along their long axes (Per-PDI2). Following selective photoexcitation of PDI, femtosecond transient absorption data for Per-PDI2 is compared to that for the single-donor, single-acceptor Per-PDI system, and the PDI2 dimer, which both have the same interchromophore geometry as Per-PDI2. The data show that electron transfer from Per to the lower exciton state of the PDI dimer is slower than that of the single PDI acceptor system. This is due to the lower free energy of the reaction for charge separation because of the electronic stabilization afforded by the excitonic coupling between the PDIs. While PDI2 was shown previously to undergo ultrafast SB-CS, the strong π-π electronic interaction of Per with the adjacent PDI in Per-PDI2 breaks the electronic symmetry of the PDI dimer, resulting in the oxidation of Per rather than SB-CS. These results show that the electronic coupling between molecules designed to accept charges produced by SB-CS in molecular dimers and the chromophores comprising the dimer must be balanced to favor SB-CS.

Local density dependent potentials for an underlying van der Waals equation of state: A simulation and density functional theory analysis.

There is an ever increasing use of local density dependent potentials in the mesoscale modeling of complex fluids. Questions remain, though, about the dependence of the thermodynamic and structural properties of such systems on the cutoff distance used to calculate these local densities. These questions are particularly acute when it comes to the stability and structure of the vapor/liquid interface. In this article, we consider local density dependent potentials derived from an underlying van der Waals equation of state. We use simulation and density functional theory to examine how the bulk thermodynamic and interfacial properties vary with the cutoff distance, rc, used to calculate the local densities. We show quantitatively how the simulation results for bulk thermodynamic properties and vapor-liquid equilibrium approach the van der Waals limit as rc increases and demonstrate a scaling law for the radial distribution function in the large rc limit. We show that the vapor-liquid interface is stable with a well-defined surface tension and that the interfacial density profile is oscillatory, except for temperatures close to critical. Finally, we show that in the large rc limit, the interfacial tension is proportional to rc and, therefore, unlike the bulk thermodynamic properties, does not approach a constant value as rc increases. We believe that these results give new insights into the properties of local density dependent potentials, in particular their unusual interfacial behavior, which is relevant for modeling complex fluids in soft matter.

Steps on the Path to Clinical Translation: A workshop by the British and Irish Chapter of the ISMRM.

The British and Irish Chapter of the International Society for Magnetic Resonance in Medicine (BIC-ISMRM) held a workshop entitled "Steps on the path to clinical translation" in Cardiff, UK, on 7th September 2022. The aim of the workshop was to promote discussion within the MR community about the problems and potential solutions for translating quantitative MR (qMR) imaging and spectroscopic biomarkers into clinical application and drug studies. Invited speakers presented the perspectives of radiologists, radiographers, clinical physicists, vendors, imaging Contract/Clinical Research Organizations (CROs), open science networks, metrologists, imaging networks, and those developing consensus methods. A round-table discussion was held in which workshop participants discussed a range of questions pertinent to clinical translation of qMR imaging and spectroscopic biomarkers. Each group summarized their findings via three main conclusions and three further questions. These questions were used as the basis of an online survey of the broader UK MR community.

Coherent Vibronic Wavepackets Show Structure-Directed Charge Flow in Host-Guest Donor-Acceptor Complexes.

Designing and controlling charge transfer (CT) pathways in organic semiconductors are important for solar energy applications. To be useful, a photogenerated, Coulombically bound CT exciton must further separate into free charge carriers; direct observations of the detailed CT relaxation pathways, however, are lacking. Here, photoinduced CT and relaxation dynamics in three host-guest complexes, where a perylene (Per) electron donor guest is incorporated into two symmetric and one asymmetric extended viologen cyclophane acceptor hosts, are presented. The central ring in the extended viologen is either p-phenylene (ExV2+) or electron-rich 2,5-dimethoxy-p-phenylene (ExMeOV2+), resulting in two symmetric cyclophanes with unsubstituted or methoxy-substituted central rings, ExBox4+ and ExMeOBox4+, respectively, and an asymmetric cyclophane with one of the central viologen rings being methoxylated ExMeOVBox4+. Upon photoexcitation, the asymmetric host-guest ExMeOVBox4+ ⊃ Per complex exhibits directional CT toward the energetically unfavorable methoxylated side due to structural restrictions that facilitate strong interactions between the Per donor and the ExMeOV2+ side. The CT state relaxation pathways are probed using ultrafast optical spectroscopy by focusing on coherent vibronic wavepackets, which are used to identify CT relaxations along charge localization and vibronic decoherence coordinates. Specific low- and high-frequency nuclear motions are direct indicators of a delocalized CT state and the degree of CT character. Our results show that the CT pathway can be controlled by subtle chemical modifications of the acceptor host in addition to illustrating how coherent vibronic wavepackets can be used to probe the nature and time evolution of the CT states.

First-in-human technique translation of oxygen-enhanced MRI to an MR Linac system in patients with head and neck cancer.

BACKGROUND AND PURPOSE: Tumour hypoxia is prognostic in head and neck cancer (HNC), associated with poor loco-regional control, poor survival and treatment resistance. The advent of hybrid MRI - radiotherapy linear accelerator or 'MR Linac' systems - could permit imaging for treatment adaptation based on hypoxic status. We sought to develop oxygen-enhanced MRI (OE-MRI) in HNC and translate the technique onto an MR Linac system. MATERIALS AND METHODS: MRI sequences were developed in phantoms and 15 healthy participants. Next, 14 HNC patients (with 21 primary or local nodal tumours) were evaluated. Baseline tissue longitudinal relaxation time (T1) was measured alongside the change in 1/T1 (termed ΔR1) between air and oxygen gas breathing phases. We compared results from 1.5 T diagnostic MR and MR Linac systems. RESULTS: Baseline T1 had excellent repeatability in phantoms, healthy participants and patients on both systems. Cohort nasal concha oxygen-induced ΔR1 significantly increased (p < 0.0001) in healthy participants demonstrating OE-MRI feasibility. ΔR1 repeatability coefficients (RC) were 0.023-0.040 s-1 across both MR systems. The tumour ΔR1 RC was 0.013 s-1 and the within-subject coefficient of variation (wCV) was 25% on the diagnostic MR. Tumour ΔR1 RC was 0.020 s-1 and wCV was 33% on the MR Linac. ΔR1 magnitude and time-course trends were similar on both systems. CONCLUSION: We demonstrate first-in-human translation of volumetric, dynamic OE-MRI onto an MR Linac system, yielding repeatable hypoxia biomarkers. Data were equivalent on the diagnostic MR and MR Linac systems. OE-MRI has potential to guide future clinical trials of biology guided adaptive radiotherapy.

Two-Dimensional Electronic Spectroscopy Reveals Vibrational Modes Coupled to Charge Transfer in a Julolidine-BODIPY Dyad.

Understanding charge transfer (CT) dynamics in molecular donor-acceptor (D-A) dyads can provide insight into developing efficient D-A molecules for capturing solar energy. Here, we characterize the excited-state evolution of a julolidine-BODIPY (Jul-BD) D-A system with an emissive CT state using time-resolved fluorescence, femtosecond transient absorption, and two-dimensional electronic spectroscopies. Comparison of these results with those from phenyl-BODIPY (Ph-BD) allows us to identify the dynamics at play during CT state formation and its subsequent conversion to either a fully charge-separated or triplet state. Photoexcitation of Jul-BD in tetrahydrofuran results in the formation of an initial emissive CT state that relaxes before fully charge-separating. In contrast, Jul-BD in toluene exhibits similar CT state dynamics, albeit at slower timescales, before decaying to a terminal triplet species. Quantum beat analysis at early times in both solvents shows several vibronic modes, which are corroborated using density functional theory (DFT) calculations. For Ph-BD, a single 220 cm-1 compression mode about the single bond linking the phenyl to BODIPY modulates their orbital overlap. Three active vibronic modes, 147, 174, and 214 cm-1, are found in Jul-BD, regardless of the dielectric constant of the medium. These motions correspond to compression and torsional motions along the single bond joining Jul to BD and are responsible for the evolution of the spontaneous and stimulated emission features in the time-resolved spectroscopic data, which is further supported by time-dependent DFT calculations of the steady-state absorption spectrum of the Jul-BD as a function of increasing D-A dihedral core angle. These findings show how torsional and compression motions can play a pivotal role in intramolecular CT between a D and an A linked by a single bond.

Encapsulating Semiconductor Quantum Dots in Supramolecular Cages Enables Ultrafast Guest-Host Electron and Vibrational Energy Transfer.

In the field of supramolecular chemistry, host-guest systems have been extensively explored to encapsulate a wide range of substrates, owing to emerging functionalities in nanoconfined space that cannot be achieved in dilute solutions. However, host-guest chemistry is still limited to encapsulation of small guests. Herein, we construct a water-soluble metallo-supramolecular hexagonal prism with a large hydrophobic cavity by anchoring multiple polyethylene glycol chains onto the building blocks. Then, assembled prisms are able to encapsulate quantum dots (QDs) with diameters of less than 5.0 nm. Furthermore, we find that the supramolecular cage around each QD strongly modifies the photophysics of the QD by universally increasing the rates of QD relaxation processes via ultrafast electron and vibrational energy transfer. Taken together, these efforts expand the scope of substrates in host-guest systems and provide a new approach to tune the optical properties of QDs.

Multi-Maintenance Olaparib Therapy in Relapsed, Germline BRCA1/2-Mutant High-Grade Serous Ovarian Cancer (MOLTO): A Phase II Trial.

PURPOSE: A single maintenance course of a PARP inhibitor (PARPi) improves progression-free survival (PFS) in germline BRCA1/2-mutant high-grade serous ovarian cancer (gBRCAm-HGSOC). The feasibility of a second maintenance course of PARPi was unknown. PATIENTS AND METHODS: Phase II trial with two entry points (EP1, EP2). Patients were recruited prior to rechallenge platinum. Patients with relapsed, gBRCAm-HGSOC were enrolled at EP1 if they were PARPi-naïve. Patients enrolled at EP2 had received their first course of olaparib prior to trial entry. EP1 patients were retreated with olaparib after RECIST complete/partial response (CR/PR) to platinum. EP2 patients were retreated with olaparib ± cediranib after RECIST CR/PR/stable disease to platinum and according to the platinum-free interval. Co-primary outcomes were the proportion of patients who received a second course of olaparib and the proportion who received olaparib retreatment for ≥6 months. Functional homologous recombination deficiency (HRD), somatic copy-number alteration (SCNA), and BRCAm reversions were investigated in tumor and liquid biopsies. RESULTS: Twenty-seven patients were treated (EP1 = 17, EP2 = 10), and 19 were evaluable. Twelve patients (63%) received a second course of olaparib and 4 received olaparib retreatment for ≥6 months. Common grade ≥2 adverse events during olaparib retreatment were anemia, nausea, and fatigue. No cases of MDS/AML occurred. Mean duration of olaparib treatment and retreatment differed (12.1 months vs. 4.4 months; P < 0.001). Functional HRD and SCNA did not predict PFS. A BRCA2 reversion mutation was detected in a post-olaparib liquid biopsy. CONCLUSIONS: A second course of olaparib can be safely administered to women with gBRCAm-HGSOC but is only modestly efficacious. See related commentary by Gonzalez-Ochoa and Oza, p. 2563.

Criteria for the translation of radiomics into clinically useful tests.

Computer-extracted tumour characteristics have been incorporated into medical imaging computer-aided diagnosis (CAD) algorithms for decades. With the advent of radiomics, an extension of CAD involving high-throughput computer-extracted quantitative characterization of healthy or pathological structures and processes as captured by medical imaging, interest in such computer-extracted measurements has increased substantially. However, despite the thousands of radiomic studies, the number of settings in which radiomics has been successfully translated into a clinically useful tool or has obtained FDA clearance is comparatively small. This relative dearth might be attributable to factors such as the varying imaging and radiomic feature extraction protocols used from study to study, the numerous potential pitfalls in the analysis of radiomic data, and the lack of studies showing that acting upon a radiomic-based tool leads to a favourable benefit-risk balance for the patient. Several guidelines on specific aspects of radiomic data acquisition and analysis are already available, although a similar roadmap for the overall process of translating radiomics into tools that can be used in clinical care is needed. Herein, we provide 16 criteria for the effective execution of this process in the hopes that they will guide the development of more clinically useful radiomic tests in the future.

Standardised lesion segmentation for imaging biomarker quantitation: a consensus recommendation from ESR and EORTC.

BACKGROUND: Lesion/tissue segmentation on digital medical images enables biomarker extraction, image-guided therapy delivery, treatment response measurement, and training/validation for developing artificial intelligence algorithms and workflows. To ensure data reproducibility, criteria for standardised segmentation are critical but currently unavailable. METHODS: A modified Delphi process initiated by the European Imaging Biomarker Alliance (EIBALL) of the European Society of Radiology (ESR) and the European Organisation for Research and Treatment of Cancer (EORTC) Imaging Group was undertaken. Three multidisciplinary task forces addressed modality and image acquisition, segmentation methodology itself, and standards and logistics. Devised survey questions were fed via a facilitator to expert participants. The 58 respondents to Round 1 were invited to participate in Rounds 2-4. Subsequent rounds were informed by responses of previous rounds. RESULTS/CONCLUSIONS: Items with ≥ 75% consensus are considered a recommendation. These include system performance certification, thresholds for image signal-to-noise, contrast-to-noise and tumour-to-background ratios, spatial resolution, and artefact levels. Direct, iterative, and machine or deep learning reconstruction methods, use of a mixture of CE marked and verified research tools were agreed and use of specified reference standards and validation processes considered essential. Operator training and refreshment were considered mandatory for clinical trials and clinical research. Items with a 60-74% agreement require reporting (site-specific accreditation for clinical research, minimal pixel number within lesion segmented, use of post-reconstruction algorithms, operator training refreshment for clinical practice). Items with ≤ 60% agreement are outside current recommendations for segmentation (frequency of system performance tests, use of only CE-marked tools, board certification of operators, frequency of operator refresher training). Recommendations by anatomical area are also specified.

Bayesian methods provide a practical real-world evidence framework for evaluating the impact of changes in radiotherapy.

PURPOSE: Retrospective studies have identified a link between the average set-up error of lung cancer patients treated with image-guided radiotherapy (IGRT) and survival. The IGRT protocol was subsequently changed to reduce the action threshold. In this study, we use a Bayesian approach to evaluate the clinical impact of this change to practice using routine 'real-world' patient data. METHODS AND MATERIALS: Two cohorts of NSCLC patients treated with IGRT were compared: pre-protocol change (N = 780, 5 mm action threshold) and post-protocol change (N = 411, 2 mm action threshold). Survival models were fitted to each cohort and changes in the hazard ratios (HR) associated with residual set-up errors was assessed. The influence of using an uninformative and a skeptical prior in the model was investigated. RESULTS: Following the reduction of the action threshold, the HR for residual set-up error towards the heart was reduced by up to 10%. Median patient survival increased for patients with set-up errors towards the heart, and remained similar for patients with set-up errors away from the heart. Depending on the prior used, a residual hazard ratio may remain. CONCLUSIONS: Our analysis found a reduced hazard of death and increased survival for patients with residual set-up errors towards versus away from the heart post-protocol change. This study demonstrates the value of a Bayesian approach in the assessment of technical changes in radiotherapy practice and supports the consideration of adopting this approach in further prospective evaluations of changes to clinical practice.

Clinical translation of quantitative magnetic resonance imaging biomarkers - An overview and gap analysis of current practice.

PURPOSE: This overview of the current landscape of quantitative magnetic resonance imaging biomarkers (qMR IBs) aims to support the standardisation of academic IBs to assist their translation to clinical practice. METHODS: We used three complementary approaches to investigate qMR IB use and quality management practices within the UK: 1) a literature search of qMR and quality management terms during 2011-2015 and 2016-2020; 2) a database search for clinical research studies using qMR IBs during 2016-2020; and 3) a survey to ascertain the current availability and quality management practices for clinical MRI scanners and associated equipment at research institutions across the UK. RESULTS: The analysis showed increased use of all qMR methods between the periods 2011-2015 and 2016-2020 and diffusion-tensor MRI and volumetry to be popular methods. However, the "translation ratio" of journal articles to clinical research studies was higher for qMR methods that have evidence of clinical translation via a commercial route, such as fat fraction and T2 mapping. The number of journal articles citing quality management terms doubled between the periods 2011-2015 and 2016-2020; although, its proportion relative to all journal articles only increased by 3.0%. The survey suggested that quality assurance (QA) and quality control (QC) of data acquisition procedures are under-reported in the literature and that QA/QC of acquired data/data analysis are under-developed and lack consistency between institutions. CONCLUSIONS: We summarise current attempts to standardise and translate qMR IBs, and conclude by outlining the ideal quality management practices and providing a gap analysis between current practice and a metrological standard.

Radiomic Features From Diffusion-Weighted MRI of Retroperitoneal Soft-Tissue Sarcomas Are Repeatable and Exhibit Change After Radiotherapy.

Background: Size-based assessments are inaccurate indicators of tumor response in soft-tissue sarcoma (STS), motivating the requirement for new response imaging biomarkers for this rare and heterogeneous disease. In this study, we assess the test-retest repeatability of radiomic features from MR diffusion-weighted imaging (DWI) and derived maps of apparent diffusion coefficient (ADC) in retroperitoneal STS and compare baseline repeatability with changes in radiomic features following radiotherapy (RT). Materials and Methods: Thirty patients with retroperitoneal STS received an MR examination prior to treatment, of whom 23/30 were investigated in our repeatability analysis having received repeat baseline examinations and 14/30 patients were investigated in our post-treatment analysis having received an MR examination after completing pre-operative RT. One hundred and seven radiomic features were extracted from the full manually delineated tumor region using PyRadiomics. Test-retest repeatability was assessed using an intraclass correlation coefficient (baseline ICC), and post-radiotherapy variance analysis (post-RT-IMS) was used to compare the change in radiomic feature value to baseline repeatability. Results: For the ADC maps and DWI images, 101 and 102 features demonstrated good baseline repeatability (baseline ICC > 0.85), respectively. Forty-three and 2 features demonstrated both good baseline repeatability and a high post-RT-IMS (>0.85), respectively. Pearson correlation between the baseline ICC and post-RT-IMS was weak (0.432 and 0.133, respectively). Conclusions: The ADC-based radiomic analysis shows better test-retest repeatability compared with features derived from DWI images in STS, and some of these features are sensitive to post-treatment change. However, good repeatability at baseline does not imply sensitivity to post-treatment change.

Habitat Imaging of Tumors Enables High Confidence Sub-Regional Assessment of Response to Therapy.

Imaging biomarkers are used in therapy development to identify and quantify therapeutic response. In oncology, use of MRI, PET and other imaging methods can be complicated by spatially complex and heterogeneous tumor micro-environments, non-Gaussian data and small sample sizes. Linear Poisson Modelling (LPM) enables analysis of complex data that is quantitative and can operate in small data domains. We performed experiments in 5 mouse models to evaluate the ability of LPM to identify responding tumor habitats across a range of radiation and targeted drug therapies. We tested if LPM could identify differential biological response rates. We calculated the theoretical sample size constraints for applying LPM to new data. We then performed a co-clinical trial using small data to test if LPM could detect multiple therapeutics with both improved power and reduced animal numbers compared to conventional t-test approaches. Our data showed that LPM greatly increased the amount of information extracted from diffusion-weighted imaging, compared to cohort t-tests. LPM distinguished biological response rates between Calu6 tumors treated with 3 different therapies and between Calu6 tumors and 4 other xenograft models treated with radiotherapy. A simulated co-clinical trial using real data detected high precision per-tumor treatment effects in as few as 3 mice per cohort, with p-values as low as 1 in 10,000. These findings provide a route to simultaneously improve the information derived from preclinical imaging while reducing and refining the use of animals in cancer research.