Metabolic profiling of murine radiation-induced lung injury with Raman spectroscopy and comparative machine learning.

Radiation-induced lung injury (RILI) is a dose-limiting toxicity for cancer patients receiving thoracic radiotherapy. As such, it is important to characterize metabolic associations with the early and late stages of RILI, namely pneumonitis and pulmonary fibrosis. Recently, Raman spectroscopy has shown utility for the differentiation of pneumonitic and fibrotic tissue states in a mouse model; however, the specific metabolite-disease associations remain relatively unexplored from a Raman perspective. This work harnesses Raman spectroscopy and supervised machine learning to investigate metabolic associations with radiation pneumonitis and pulmonary fibrosis in a mouse model. To this end, Raman spectra were collected from lung tissues of irradiated/non-irradiated C3H/HeJ and C57BL/6J mice and labelled as normal, pneumonitis, or fibrosis, based on histological assessment. Spectra were decomposed into metabolic scores via group and basis restricted non-negative matrix factorization, classified with random forest (GBR-NMF-RF), and metabolites predictive of RILI were identified. To provide comparative context, spectra were decomposed and classified via principal component analysis with random forest (PCA-RF), and full spectra were classified with a convolutional neural network (CNN), as well as logistic regression (LR). Through leave-one-mouse-out cross-validation, we observed that GBR-NMF-RF was comparable to other methods by measure of accuracy and log-loss (p > 0.10 by Mann-Whitney U test), and no methodology was dominant across all classification tasks by measure of area under the receiver operating characteristic curve. Moreover, GBR-NMF-RF results were directly interpretable and identified collagen and specific collagen precursors as top fibrosis predictors, while metabolites with immune and inflammatory functions, such as serine and histidine, were top pneumonitis predictors. Further support for GBR-NMF-RF and the identified metabolite associations with RILI was found as CNN interpretation heatmaps revealed spectral regions consistent with these metabolites.

Stratification of tumour cell radiation response and metabolic signatures visualization with Raman spectroscopy and explainable convolutional neural network.

Reprogramming of cellular metabolism is a driving factor of tumour progression and radiation therapy resistance. Identifying biochemical signatures associated with tumour radioresistance may assist with the development of targeted treatment strategies to improve clinical outcomes. Raman spectroscopy (RS) can monitor post-irradiation biomolecular changes and signatures of radiation response in tumour cells in a label-free manner. Convolutional Neural Networks (CNN) perform feature extraction directly from data in an end-to-end learning manner, with high classification performance. Furthermore, recently developed CNN explainability techniques help visualize the critical discriminative features captured by the model. In this work, a CNN is developed to characterize tumour response to radiotherapy based on its degree of radioresistance. The model was trained to classify Raman spectra of three human tumour cell lines as radiosensitive (LNCaP) or radioresistant (MCF7, H460) over a range of treatment doses and data collection time points. Additionally, a method based on Gradient-Weighted Class Activation Mapping (Grad-CAM) was used to determine response-specific salient Raman peaks influencing the CNN predictions. The CNN effectively classified the cell spectra, with accuracy, sensitivity, specificity, and F1 score exceeding 99.8%. Grad-CAM heatmaps of H460 and MCF7 cell spectra (radioresistant) exhibited high contributions from Raman bands tentatively assigned to glycogen, amino acids, and nucleic acids. Conversely, heatmaps of LNCaP cells (radiosensitive) revealed activations at lipid and phospholipid bands. Finally, Grad-CAM variable importance scores were derived for glycogen, asparagine, and phosphatidylcholine, and we show that their trends over cell line, dose, and acquisition time agreed with previously established models. Thus, the CNN can accurately detect biomolecular differences in the Raman spectra of tumour cells of varying radiosensitivity without requiring manual feature extraction. Finally, Grad-CAM may help identify metabolic signatures associated with the observed categories, offering the potential for automated clinical tumour radiation response characterization.

A postimplant dosimetry simulation framework for robustness evaluation in permanent breast seed implant brachytherapy.

BACKGROUND: Permanent breast seed implant (PBSI) brachytherapy is a promising treatment that has the potential to be widely utilized with increased standardization, optimization, and robustness. Excellent early efficacy and very high patient acceptance were reported, however, to further evaluate and improve planning strategies, a framework to quantify plan robustness to implant uncertainties is necessary. PURPOSE: In this study, we aim to quantify clinical seed displacement using an automated algorithm and develop and validate a PBSI post-implant dosimetry simulation framework to evaluate PBSI plan robustness to implant uncertainties. METHODS AND MATERIALS: Clinical PBSI seed displacements were quantified for 63 consecutive patients. A PBSI simulator was developed in Matlab (2020) by resampling clinical seed displacements and computing a range of possible post-implant dosimetry outcomes under various seed displacement scenarios. Simulations were performed retrospectively on 63 previous clinical plans to evaluate plan robustness to seed displacement. RESULTS: Mean seed displacement for the whole cohort was 10 ± 6 mm. A clinical seed displacement database was established and a user interface was developed for the simulation framework. For all clinical plans, the median (range) value of simulated median ETV V90 in various seed displacement scenarios was 97.8% (87.5-100%). CONCLUSIONS: A PBSI postimplant dosimetry simulation framework was developed and validated. Simulation results showed that the current PTV planning margin is sufficient to provide adequate postimplant dose coverage of ETV. This simulator can be used to evaluate plan robustness to seed displacement and will facilitate future research in improving PBSI planning methods.

Reconstruction of Raman Spectra of Biochemical Mixtures Using Group and Basis Restricted Non-Negative Matrix Factorization.

Raman spectroscopy is a useful tool for obtaining biochemical information from biological samples. However, interpretation of Raman spectroscopy data in order to draw meaningful conclusions related to the biochemical make up of cells and tissues is often difficult and could be misleading if care is not taken in the deconstruction of the spectral data. Our group has previously demonstrated the implementation of a group- and basis-restricted non-negative matrix factorization (GBR-NMF) framework as an alternative to more widely used dimensionality reduction techniques such as principal component analysis (PCA) for the deconstruction of Raman spectroscopy data as related to radiation response monitoring in both cellular and tissue data. While this method provides better biological interpretability of the Raman spectroscopy data, there are some important factors which must be considered in order to provide the most robust GBR-NMF model. We here evaluate and compare the accuracy of a GBR-NMF model in the reconstruction of three mixture solutions of known concentrations. The factors assessed include the effect of solid versus solutions bases spectra, the number of unconstrained components used in the model, the tolerance of different signal to noise thresholds, and how different groups of biochemicals compare to each other. The robustness of the model was assessed by how well the relative concentration of each individual biochemical in the solution mixture is reflected in the GBR-NMF scores obtained. We also evaluated how well the model can reconstruct original data, both with and without the inclusion of an unconstrained component. Overall, we found that solid bases spectra were generally comparable to solution bases spectra in the GBR-NMF model for all groups of biochemicals. The model was found to be relatively tolerant of high levels of noise in the mixture solutions using solid bases spectra. Additionally, the inclusion of an unconstrained component did not have a significant effect on the deconstruction, on the condition that all biochemicals in the mixture were included as bases chemicals in the model. We also report that some groups of biochemicals achieve a more accurate deconstruction using GBR-NMF than others, likely due to similarity in the individual bases spectra.

Raman spectroscopy and convolutional neural networks for monitoring biochemical radiation response in breast tumour xenografts.

Tumour cells exhibit altered metabolic pathways that lead to radiation resistance and disease progression. Raman spectroscopy (RS) is a label-free optical modality that can monitor post-irradiation biomolecular signatures in tumour cells and tissues. Convolutional Neural Networks (CNN) perform automated feature extraction directly from data, with classification accuracy exceeding that of traditional machine learning, in cases where data is abundant and feature extraction is challenging. We are interested in developing a CNN-based predictive model to characterize clinical tumour response to radiation therapy based on their degree of radiosensitivity or radioresistance. In this work, a CNN architecture is built for identifying post-irradiation spectral changes in Raman spectra of tumour tissue. The model was trained to classify irradiated versus non-irradiated tissue using Raman spectra of breast tumour xenografts. The CNN effectively classified the tissue spectra, with accuracies exceeding 92.1% for data collected 3 days post-irradiation, and 85.0% at day 1 post-irradiation. Furthermore, the CNN was evaluated using a leave-one-out- (mouse, section or Raman map) validation approach to investigate its generalization to new test subjects. The CNN retained good predictive accuracy (average accuracies 83.7%, 91.4%, and 92.7%, respectively) when little to no information for a specific subject was given during training. Finally, the classification performance of the CNN was compared to that of a previously developed model based on group and basis restricted non-negative matrix factorization and random forest (GBR-NMF-RF) classification. We found that CNN yielded higher classification accuracy, sensitivity, and specificity in mice assessed 3 days post-irradiation, as compared with the GBR-NMF-RF approach. Overall, the CNN can detect biochemical spectral changes in tumour tissue at an early time point following irradiation, without the need for previous manual feature extraction. This study lays the foundation for developing a predictive framework for patient radiation response monitoring.

Raman microspectroscopy and machine learning for use in identifying radiation-induced lung toxicity.

OBJECTIVE: In this work, we explore and develop a method that uses Raman spectroscopy to measure and differentiate radiation induced toxicity in murine lungs with the goal of setting the foundation for a predictive disease model. METHODS: Analysis of Raman tissue data is achieved through a combination of techniques. We first distinguish between tissue measurements and air pockets in the lung by using group and basis restricted non-negative matrix factorization. We then analyze the tissue spectra using sparse multinomial logistic regression to discriminate between fibrotic gradings. Model validation is achieved by splitting the data into a training set containing 70% of the data and a test set with the remaining 30%; classification accuracy is used as the performance metric. We also explore several other potential classification tasks wherein the response considered is the grade of pneumonitis and fibrosis sickness. RESULTS: A classification accuracy of 91.6% is achieved on the test set of fibrotic gradings, illustrating the ability of Raman measurements to detect differing levels of fibrotic disease among the murine lungs. It is also shown via further modeling that coarser consideration of fibrotic grading via binning (ie. 'Low', 'Medium', 'High') does not degrade performance. Finally, we consider preliminary models for pneumonitis discrimination using the same methodologies.

Radiation treatment response and hypoxia biomarkers revealed by machine learning assisted Raman spectroscopy in tumour cells and xenograft tissues.

Recent advancements in anatomical imaging of tumours as treatment targets have led to improvements in RT. However, it is unlikely that improved anatomical imaging alone will be the sole driver for new advances in personalised RT. Biochemically based radiobiological information is likely to be required for next-generation improvements in the personalisation of radiotherapy dose prescriptions to individual patients. In this paper, we use Raman spectroscopy (RS), an optical technique, to monitor individual biochemical response to radiation within a tumour microenvironment. We spatially correlate individual biochemical responses to augmentatively derived hypoxic maps within the tumour microenvironment. Furthermore, we pair RS with a data analytical framework combining (i) group and basis restricted non-negative matrix factorization (GBR-NMF), (ii) a random forest (RF) classifier, (iii) and a feature metric importance calculation method, Shapley Additive exPlanations (SHAP), in order to ascertain the relative importance of individual biochemicals in describing the overall biological response as observed with RS. The current study found that the GBR-NMF-RF-SHAP model helped identify a wide range of radiation response biomarkers and hypoxia indicators (e.g., glycogen, lipids, DNA, amino acids) in H460 human lung cancer cells and H460 xenografts. Correlations between the hypoxic regions and Raman chemical biomarkers (e.g., glycogen, alanine, and arginine) were also identified in H460 xenografts. To summarize, GBR-NMF-RF-SHAP combined with RS can be applied to monitor the RT-induced biochemical response within cellular and tissue environments. Individual biochemicals were identified that (i) contributed to overall biological response to radiation, and (ii) spatially correlated with hypoxic regions of the tumour. RS combined with our analytical pipeline shows promise for further understanding of individual biochemical dynamics in radiation response for use in cancer therapy.

Prediction of disease progression indicators in prostate cancer patients receiving HDR-brachytherapy using Raman spectroscopy and semi-supervised learning: a pilot study.

This work combines Raman spectroscopy (RS) with supervised learning methods-group and basis restricted non-negative matrix factorisation (GBR-NMF) and linear discriminant analysis (LDA)-to aid in the prediction of clinical indicators of disease progression in a cohort of 9 patients receiving high dose rate brachytherapy (HDR-BT) as the primary treatment for intermediate risk (D'Amico) prostate adenocarcinoma. The combination of Raman spectroscopy and GBR-NMF-sparseLDA modelling allowed for the prediction of the following clinical information; Gleason score, cancer of the prostate risk assessment (CAPRA) score of pre-treatment biopsies and a Ki67 score of < 3.5% or > 3.5% in post treatment biopsies. The three clinical indicators of disease progression investigated in this study were predicted using a single set of Raman spectral data acquired from each individual biopsy, obtained pre HDR-BT treatment. This work highlights the potential of RS, combined with supervised learning, as a tool for the prediction of multiple types of clinically relevant information to be acquired simultaneously using pre-treatment biopsies, therefore opening up the potential for avoiding the need for multiple immunohistochemistry (IHC) staining procedures (H&E, Ki67) and blood sample analysis (PSA) to aid in CAPRA scoring.

Raman spectroscopy and supervised learning as a potential tool to identify high-dose-rate-brachytherapy induced biochemical profiles of prostate cancer.

High-dose-rate-brachytherapy (HDR-BT) is an increasingly attractive alternative to external beam radiation-therapy for patients with intermediate risk prostate cancer. Despite this, no bio-marker based method currently exists to monitor treatment response, and the changes which take place at the biochemical level in hypo-fractionated HDR-BT remain poorly understood. The aim of this pilot study is to assess the capability of Raman spectroscopy (RS) combined with principal component analysis (PCA) and random-forest classification (RF) to identify radiation response profiles after a single dose of 13.5 Gy in a cohort of nine patients. We here demonstrate, as a proof-of-concept, how RS-PCA-RF could be utilised as an effective tool in radiation response monitoring, specifically assessing the importance of low variance PCs in complex sample sets. As RS provides information on the biochemical composition of tissue samples, this technique could provide insight into the changes which take place on the biochemical level, as result of HDR-BT treatment.

Group and Basis Restricted Non-Negative Matrix Factorization and Random Forest for Molecular Histotype Classification and Raman Biomarker Monitoring in Breast Cancer.

Raman spectroscopy is a non-invasive optical technique that can be used to investigate biochemical information embedded in cells and tissues exposed to ionizing radiation used in cancer therapy. Raman spectroscopy could potentially be incorporated in personalized radiation treatment design as a tool to monitor radiation response in at the metabolic level. However, tracking biochemical dynamics remains challenging for Raman spectroscopy. Here we developed a novel analytical framework by combining group and basis restricted non-negative matrix factorization and random forest (GBR-NMF-RF). This framework can monitor radiation response profiles in different molecular histotypes and biochemical dynamics in irradiated breast cancer cells. Five subtypes of; human breast cancer (MCF-7, BT-474, MDA-MB-230, and SK-BR-3) and normal cells derived from human breast tissue (MCF10A) which had been exposed to ionizing radiation were tested in this framework. Reference Raman spectra of 20 biochemicals were collected and used as the constrained Raman biomarkers in the GBR-NMF-RF framework. We obtained scores for individual biochemicals corresponding to the contribution of each Raman reference spectrum to each spectrum obtained from the five cell types. A random forest classifier was then fitted to the chemical scores for performing molecular histotype classifications (HER2, PR, ER, Ki67, and cancer versus non-cancer) and assessing the importance of the Raman biochemical basis spectra for each classification test. Overall, the GBR-NMF-RF framework yields classification results with high accuracy (>97%), high sensitivity (>97%), and high specificity (>97%). Variable importance calculated in the random forest model indicated high contributions from glycogen and lipids (cholesterol, phosphatidylserine, and stearic acid) in molecular histotype classifications.

Raman spectroscopy and group and basis-restricted non negative matrix factorisation identifies radiation induced metabolic changes in human cancer cells.

This work combines single cell Raman spectroscopy (RS) with group and basis restricted non-negative matrix factorisation (GBR-NMF) to identify individual biochemical changes associated with radiation exposure in three human cancer cell lines. The cell lines analysed were derived from lung (H460), breast (MCF7) and prostate (LNCaP) tissue and are known to display varying degrees of radio sensitivity due to the inherent properties of each cell type. The GBR-NMF approach involves the deconstruction of Raman spectra into component biochemical bases using a library of Raman spectra of known biochemicals present in the cells. Subsequently, scores are obtained on each of these bases which can be directly correlated with the contribution of each chemical to the overall Raman spectrum. We validated GBR-NMF through the correlation of GBR-NMF-derived glycogen scores with scores that were previously observed using principal component analysis (PCA). Phosphatidylcholine, glucose, arginine and asparagine showed a distinct differential score pattern between radio-resistant and radio-sensitive cell types. In summary, the GBR-NMF approach allows for the monitoring of individual biochemical radiation-response dynamics previously unattainable with more traditional PCA-based approaches.

Monitor Ionizing Radiation-Induced Cellular Responses with Raman Spectroscopy, Non-Negative Matrix Factorization, and Non-Negative Least Squares.

Radiation therapy (RT) is one of the most commonly prescribed cancer treatments. New tools that can accurately monitor and evaluate individual patient responses would be a major advantage and lend to the implementation of personalized treatment plans. In this study, Raman spectroscopy (RS) was applied to examine radiation-induced cellular responses in H460, MCF7, and LNCaP cancer cell lines across different dose levels and times post-irradiation. Previous Raman data analysis was conducted using principal component analysis (PCA), which showed the ability to extract biological information of glycogen. In the current studies, the use of non-negative matrix factorization (NMF) allowed for the discovery of multiplexed biological information, specifically uncovering glycogen-like and lipid-like component bases. The corresponding scores of glycogen and previously unidentified lipids revealed the content variations of these two chemicals in the cellular data. The NMF decomposed glycogen and lipid-like bases were able to separate the cancer cell lines into radiosensitive and radioresistant groups. A further lipid phenotype investigation was also attempted by applying non-negative least squares (NNLS) to the lipid-like bases decomposed individually from three cell lines. Qualitative differences found in lipid weights for each lipid-like basis suggest the lipid phenotype differences in the three tested cancer cell lines. Collectively, this study demonstrates that the application of NMF and NNLS on RS data analysis to monitor ionizing radiation-induced cellular responses can yield multiplexed biological information on bio-response to RT not revealed by conventional chemometric approaches.