Multimodal Machine Learning in Image-Based and Clinical Biomedicine: Survey and Prospects.

Machine learning (ML) applications in medical artificial intelligence (AI) systems have shifted from traditional and statistical methods to increasing application of deep learning models. This survey navigates the current landscape of multimodal ML, focusing on its profound impact on medical image analysis and clinical decision support systems. Emphasizing challenges and innovations in addressing multimodal representation, fusion, translation, alignment, and co-learning, the paper explores the transformative potential of multimodal models for clinical predictions. It also highlights the need for principled assessments and practical implementation of such models, bringing attention to the dynamics between decision support systems and healthcare providers and personnel. Despite advancements, challenges such as data biases and the scarcity of "big data" in many biomedical domains persist. We conclude with a discussion on principled innovation and collaborative efforts to further the mission of seamless integration of multimodal ML models into biomedical practice.

A snapshot of medical physics practice patterns.

A large number of surveys have been sent to the medical physics community addressing many clinical topics for which the medical physicist is, or may be, responsible. Each survey provides an insight into clinical practice relevant to the medical physics community. The goal of this study was to create a summary of these surveys giving a snapshot of clinical practice patterns. Surveys used in this study were created using SurveyMonkey and distributed between February 6, 2013 and January 2, 2018 via the MEDPHYS and MEDDOS listserv groups. The format of the surveys included questions that were multiple choice and free response. Surveys were included in this analysis if they met the following criteria: more than 20 responses, relevant to radiation therapy physics practice, not single-vendor specific, and formatted as multiple-choice questions (i.e., not exclusively free-text responses). Although the results of free response questions were not explicitly reported, they were carefully reviewed, and the responses were considered in the discussion of each topic. Two-hundred and fifty-two surveys were available, of which 139 passed the inclusion criteria. The mean number of questions per survey was 4. The mean number of respondents per survey was 63. Summaries were made for the following topics: simulation, treatment planning, electron treatments, linac commissioning and quality assurance, setup and treatment verification, IMRT and VMAT treatments, SRS/SBRT, breast treatments, prostate treatments, brachytherapy, TBI, facial lesion treatments, clinical workflow, and after-hours/emergent treatments. We have provided a coherent overview of medical physics practice according to surveys conducted over the last 5 yr, which will be instructive for medical physicists.

Cost-effective immobilization for whole brain radiation therapy.

To investigate the inter- and intra-fraction motion associated with the use of a low-cost tape immobilization technique as an alternative to thermoplastic immobilization masks for whole-brain treatments. The results of this study may be of interest to clinical staff with severely limited resources (e.g., in low-income countries) and also when treating patients who cannot tolerate standard immobilization masks. Setup reproducibility of eight healthy volunteers was assessed for two different immobilization techniques. (a) One strip of tape was placed across the volunteer's forehead and attached to the sides of the treatment table. (b) A second strip was added to the first, under the chin, and secured to the table above the volunteer's head. After initial positioning, anterior and lateral photographs were acquired. Volunteers were positioned five times with each technique to allow calculation of inter-fraction reproducibility measurements. To estimate intra-fraction reproducibility, 5-minute anterior and lateral videos were taken for each technique per volunteer. An in-house software was used to analyze the photos and videos to assess setup reproducibility. The maximum intra-fraction displacement for all volunteers was 2.8 mm. Intra-fraction motion increased with time on table. The maximum inter-fraction range of positions for all volunteers was 5.4 mm. The magnitude of inter-fraction and intra-fraction motion found using the "1-strip" and "2-strip" tape immobilization techniques was comparable to motion restrictions provided by a thermoplastic mask for whole-brain radiotherapy. The results suggest that tape-based immobilization techniques represent an economical and useful alternative to the thermoplastic mask.

Association of graph-based spatial features with overall survival status of glioblastoma patients.

Glioblastoma is the most common malignant brain tumor with less than 15 months median survival. To aid prognosis, there is a need for decision tools that leverage diagnostic modalities such as MRI to inform survival. In this study, we examine higher-order spatial proximity characteristics from habitats and propose two graph-based methods (minimum spanning tree and graph run-length matrix) to characterize spatial heterogeneity over tumor MRI-derived intensity habitats and assess their relationships with overall survival as well as the immune signature status of patients with glioblastoma. A data set of 74 patients was studied based on the availability of post-contrast T1-weighted and T2-weighted fluid attenuated inversion recovery (FLAIR) image data in The Cancer Image Archive (TCIA). We assessed the predictive value of MST- and GRLM-derived features from 2D images for prediction of 12-month survival status and immune signature status of patients with glioblastoma via a receiver operating characteristic curve analysis. For 12-month survival prediction using MST-based method, sensitivity and specificity were 0.82 and 0.79 respectively. For GRLM-based method, sensitivity and specificity were 0.73 and 0.77 respectively. For immune status, sensitivity and specificity were 0.91 and 0.69, respectively, for the GRLM-based method with an immune effector. Our results show that the proposed MST- and GRLM-derived features are predictive of 12-month survival status as well as the immune signature status of patients with glioblastoma. To our knowledge, this is the first application of MST- and GRLM-based proximity analyses for the study of radiologically-defined tumor habitats in glioblastoma.

Low-parameter supervised learning models can discriminate pseudoprogression and true progression in non-perfusion-based MRI.

Discrimination of pseudoprogression and true progression is one challenge to the treatment of malignant gliomas. Although some techniques such as circulating tumor DNA (ctDNA) and perfusion-weighted imaging (PWI) demonstrate promise in distinguishing PsP from TP, we investigate robust and replicable alternatives to distinguish the two entities based on more widely-available media. In this study, we use low-parametric supervised learning techniques based on geographically-weighted regression (GWR) to investigate the utility of both conventional MRI sequences as well as a diffusion-weighted sequence (apparent diffusion coefficient or ADC) in the discrimination of PsP v TP. GWR applied to MRI modality pairs is a unique approach for small sample sizes and is a novel approach in this arena. From our analysis, all modality pairs involving ADC maps, and those involving post-contrast T1 regressed onto T2 showed potential promise. This work on ADC data adds to a growing body of research suggesting the predictive benefits of ADC, and suggests further research on the relationships between post-contrast T1 and T2.Clinical relevance- Few studies have investigated predictive potential of conventional MRI and ADC to detect PsP. Our study adds to the growing research on the topic and presents a new perspective to research by exploiting the utility of ADC in PsP v TP distinction. In addition, our GWR methodology for low-parametric supervised computer vision models demonstrates a unique approach for image processing of small sample sizes.

Clustering-based spatial analysis (CluSA) framework through graph neural network for chronic kidney disease prediction using histopathology images.

Machine learning applied to digital pathology has been increasingly used to assess kidney function and diagnose the underlying cause of chronic kidney disease (CKD). We developed a novel computational framework, clustering-based spatial analysis (CluSA), that leverages unsupervised learning to learn spatial relationships between local visual patterns in kidney tissue. This framework minimizes the need for time-consuming and impractical expert annotations. 107,471 histopathology images obtained from 172 biopsy cores were used in the clustering and in the deep learning model. To incorporate spatial information over the clustered image patterns on the biopsy sample, we spatially encoded clustered patterns with colors and performed spatial analysis through graph neural network. A random forest classifier with various groups of features were used to predict CKD. For predicting eGFR at the biopsy, we achieved a sensitivity of 0.97, specificity of 0.90, and accuracy of 0.95. AUC was 0.96. For predicting eGFR changes in one-year, we achieved a sensitivity of 0.83, specificity of 0.85, and accuracy of 0.84. AUC was 0.85. This study presents the first spatial analysis based on unsupervised machine learning algorithms. Without expert annotation, CluSA framework can not only accurately classify and predict the degree of kidney function at the biopsy and in one year, but also identify novel predictors of kidney function and renal prognosis.

Head and Neck Radiation Therapy Patterns of Practice Variability Identified as a Challenge to Real-World Big Data: Results From the Learning from Analysis of Multicentre Big Data Aggregation (LAMBDA) Consortium.

Purpose: Outside of randomized clinical trials, it is difficult to develop clinically relevant evidence-based recommendations for radiation therapy (RT) practice guidelines owing to lack of comprehensive real-world data. To address this knowledge gap, we formed the Learning from Analysis of Multicenter Big Data Aggregation consortium to cooperatively implement RT data standardization, develop software solutions for data analysis, and recommend clinical practice change based on real-world data analyzed. The first phase of this "Big Data" study aimed at characterizing variability in clinical practice patterns of dosimetric data for organs at risk (OARs) that would undermine subsequent use of large-scale, electronically aggregated data to characterize associations with outcomes. Evidence from this study was used as the basis for practical recommendations to improve data quality. Methods and Materials: Dosimetric details of patients with head and neck cancer treated with radiation therapy between 2014 and 2019 were analyzed. Institutional patterns of practice were characterized, including structure nomenclature, volumes, and frequency of contouring. Dose volume histogram (DVH) distributions were characterized and compared with institutional constraints and literature values. Results: Plans for 4664 patients treated to a mean plan dose of 64.4 ± 13.2 Gy in 32 ± 4 fractions were aggregated. Before implementation of TG-263 guidelines in each institution, there was variability in OAR nomenclature across institutions and structures. With evidence from this study, we identified a targeted and practical set of recommendations aimed at improving the quality of real-world data. Conclusions: Quantifying similarities and differences among institutions for OAR structures and DVH metrics is the launching point for next steps to investigate potential relationships between DVH parameters and patient outcomes.

Differentiating intrinsic SERS spectra from a mixture by sampling induced composition gradient and independent component analysis.

By generating a composition gradient on a highly uniform SERS substrate and applying independent component analysis, we demonstrate that one can extract the intrinsic SERS spectrum of individual components from SERS spectra obtained from a two-component mixture.

Unsupervised machine learning for identifying important visual features through bag-of-words using histopathology data from chronic kidney disease.

Pathologists use visual classification to assess patient kidney biopsy samples when diagnosing the underlying cause of kidney disease. However, the assessment is qualitative, or semi-quantitative at best, and reproducibility is challenging. To discover previously unknown features which predict patient outcomes and overcome substantial interobserver variability, we developed an unsupervised bag-of-words model. Our study applied to the C-PROBE cohort of patients with chronic kidney disease (CKD). 107,471 histopathology images were obtained from 161 biopsy cores and identified important morphological features in biopsy tissue that are highly predictive of the presence of CKD both at the time of biopsy and in one year. To evaluate the performance of our model, we estimated the AUC and its 95% confidence interval. We show that this method is reliable and reproducible and can achieve 0.93 AUC at predicting glomerular filtration rate at the time of biopsy as well as predicting a loss of function at one year. Additionally, with this method, we ranked the identified morphological features according to their importance as diagnostic markers for chronic kidney disease. In this study, we have demonstrated the feasibility of using an unsupervised machine learning method without human input in order to predict the level of kidney function in CKD. The results from our study indicate that the visual dictionary, or visual image pattern, obtained from unsupervised machine learning can predict outcomes using machine-derived values that correspond to both known and unknown clinically relevant features.

Stress Testing Pathology Models with Generated Artifacts.

BACKGROUND: Machine learning models provide significant opportunities for improvement in health care, but their "black-box" nature poses many risks. METHODS: We built a custom Python module as part of a framework for generating artifacts that are meant to be tunable and describable to allow for future testing needs. We conducted an analysis of a previously published digital pathology classification model and an internally developed kidney tissue segmentation model, utilizing a variety of generated artifacts including testing their effects. The artifacts simulated were bubbles, tissue folds, uneven illumination, marker lines, uneven sectioning, altered staining, and tissue tears. RESULTS: We found that there is some performance degradation on the tiles with artifacts, particularly with altered stains but also with marker lines, tissue folds, and uneven sectioning. We also found that the response of deep learning models to artifacts could be nonlinear. CONCLUSIONS: Generated artifacts can provide a useful tool for testing and building trust in machine learning models by understanding where these models might fail.

Radiomics feature robustness as measured using an MRI phantom.

Radiomics involves high-throughput extraction of large numbers of quantitative features from medical images and analysis of these features to predict patients' outcome and support clinical decision-making. However, radiomics features are sensitive to several factors, including scanning protocols. The purpose of this study was to investigate the robustness of magnetic resonance imaging (MRI) radiomics features with various MRI scanning protocol parameters and scanners using an MRI radiomics phantom. The variability of the radiomics features with different scanning parameters and repeatability measured using a test-retest scheme were evaluated using the coefficient of variation and intraclass correlation coefficient (ICC) for both T1- and T2-weighted images. For variability measures, the features were categorized into three groups: large, intermediate, and small variation. For repeatability measures, the average T1- and T2-weighted image ICCs for the phantom (0.963 and 0.959, respectively) were higher than those for a healthy volunteer (0.856 and 0.849, respectively). Our results demonstrated that various radiomics features are dependent on different scanning parameters and scanners. The radiomics features with a low coefficient of variation and high ICC for both the phantom and volunteer can be considered good candidates for MRI radiomics studies. The results of this study will assist current and future MRI radiomics studies.

Discriminating pseudoprogression and true progression in diffuse infiltrating glioma using multi-parametric MRI data through deep learning.

Differentiating pseudoprogression from true tumor progression has become a significant challenge in follow-up of diffuse infiltrating gliomas, particularly high grade, which leads to a potential treatment delay for patients with early glioma recurrence. In this study, we proposed to use a multiparametric MRI data as a sequence input for the convolutional neural network with the recurrent neural network based deep learning structure to discriminate between pseudoprogression and true tumor progression. In this study, 43 biopsy-proven patient data identified as diffuse infiltrating glioma patients whose disease progressed/recurred were used. The dataset consists of five original MRI sequences; pre-contrast T1-weighted, post-contrast T1-weighted, T2-weighted, FLAIR, and ADC images as well as two engineered sequences; T1post-T1pre and T2-FLAIR. Next, we used three CNN-LSTM models with a different set of sequences as input sequences to pass through CNN-LSTM layers. We performed threefold cross-validation in the training dataset and generated the boxplot, accuracy, and ROC curve, AUC from each trained model with the test dataset to evaluate models. The mean accuracy for VGG16 models ranged from 0.44 to 0.60 and the mean AUC ranged from 0.47 to 0.59. For CNN-LSTM model, the mean accuracy ranged from 0.62 to 0.75 and the mean AUC ranged from 0.64 to 0.81. The performance of the proposed CNN-LSTM with multiparametric sequence data was found to outperform the popular convolutional CNN with a single MRI sequence. In conclusion, incorporating all available MRI sequences into a sequence input for a CNN-LSTM model improved diagnostic performance for discriminating between pseudoprogression and true tumor progression.

A PET Radiomics Model to Predict Refractory Mediastinal Hodgkin Lymphoma.

First-order radiomic features, such as metabolic tumor volume (MTV) and total lesion glycolysis (TLG), are associated with disease progression in early-stage classical Hodgkin lymphoma (HL). We hypothesized that a model incorporating first- and second-order radiomic features would more accurately predict outcome than MTV or TLG alone. We assessed whether radiomic features extracted from baseline PET scans predicted relapsed or refractory disease status in a cohort of 251 patients with stage I-II HL who were managed at a tertiary cancer center. Models were developed and tested using a machine-learning algorithm. Features extracted from mediastinal sites were highly predictive of primary refractory disease. A model incorporating 5 of the most predictive features had an area under the curve (AUC) of 95.2% and total error rate of 1.8%. By comparison, the AUC was 78% for both MTV and TLG and was 65% for maximum standardize uptake value (SUVmax). Furthermore, among the patients with refractory mediastinal disease, our model distinguished those who were successfully salvaged from those who ultimately died of HL. We conclude that our PET radiomic model may improve upfront stratification of early-stage HL patients with mediastinal disease and thus contribute to risk-adapted, individualized management.

DEMARCATE: Density-based magnetic resonance image clustering for assessing tumor heterogeneity in cancer.

Tumor heterogeneity is a crucial area of cancer research wherein inter- and intra-tumor differences are investigated to assess and monitor disease development and progression, especially in cancer. The proliferation of imaging and linked genomic data has enabled us to evaluate tumor heterogeneity on multiple levels. In this work, we examine magnetic resonance imaging (MRI) in patients with brain cancer to assess image-based tumor heterogeneity. Standard approaches to this problem use scalar summary measures (e.g., intensity-based histogram statistics) that do not adequately capture the complete and finer scale information in the voxel-level data. In this paper, we introduce a novel technique, DEMARCATE (DEnsity-based MAgnetic Resonance image Clustering for Assessing Tumor hEterogeneity) to explore the entire tumor heterogeneity density profiles (THDPs) obtained from the full tumor voxel space. THDPs are smoothed representations of the probability density function of the tumor images. We develop tools for analyzing such objects under the Fisher-Rao Riemannian framework that allows us to construct metrics for THDP comparisons across patients, which can be used in conjunction with standard clustering approaches. Our analyses of The Cancer Genome Atlas (TCGA) based Glioblastoma dataset reveal two significant clusters of patients with marked differences in tumor morphology, genomic characteristics and prognostic clinical outcomes. In addition, we see enrichment of image-based clusters with known molecular subtypes of glioblastoma multiforme, which further validates our representation of tumor heterogeneity and subsequent clustering techniques.

Associating spatial diversity features of radiologically defined tumor habitats with epidermal growth factor receptor driver status and 12-month survival in glioblastoma: methods and preliminary investigation.

We analyzed the spatial diversity of tumor habitats, regions with distinctly different intensity characteristics of a tumor, using various measurements of habitat diversity within tumor regions. These features were then used for investigating the association with a 12-month survival status in glioblastoma (GBM) patients and for the identification of epidermal growth factor receptor (EGFR)-driven tumors. T1 postcontrast and T2 fluid attenuated inversion recovery images from 65 GBM patients were analyzed in this study. A total of 36 spatial diversity features were obtained based on pixel abundances within regions of interest. Performance in both the classification tasks was assessed using receiver operating characteristic (ROC) analysis. For association with 12-month overall survival, area under the ROC curve was 0.74 with confidence intervals [0.630 to 0.858]. The sensitivity and specificity at the optimal operating point ([Formula: see text]) on the ROC were 0.59 and 0.75, respectively. For the identification of EGFR-driven tumors, the area under the ROC curve (AUC) was 0.85 with confidence intervals [0.750 to 0.945]. The sensitivity and specificity at the optimal operating point ([Formula: see text]) on the ROC were 0.76 and 0.83, respectively. Our findings suggest that these spatial habitat diversity features are associated with these clinical characteristics and could be a useful prognostic tool for magnetic resonance imaging studies of patients with GBM.

Spatial Habitat Features Derived from Multiparametric Magnetic Resonance Imaging Data Are Associated with Molecular Subtype and 12-Month Survival Status in Glioblastoma Multiforme.

One of the most common and aggressive malignant brain tumors is Glioblastoma multiforme. Despite the multimodality treatment such as radiation therapy and chemotherapy (temozolomide: TMZ), the median survival rate of glioblastoma patient is less than 15 months. In this study, we investigated the association between measures of spatial diversity derived from spatial point pattern analysis of multiparametric magnetic resonance imaging (MRI) data with molecular status as well as 12-month survival in glioblastoma. We obtained 27 measures of spatial proximity (diversity) via spatial point pattern analysis of multiparametric T1 post-contrast and T2 fluid-attenuated inversion recovery MRI data. These measures were used to predict 12-month survival status (≤12 or >12 months) in 74 glioblastoma patients. Kaplan-Meier with receiver operating characteristic analyses was used to assess the relationship between derived spatial features and 12-month survival status as well as molecular subtype status in patients with glioblastoma. Kaplan-Meier survival analysis revealed that 14 spatial features were capable of stratifying overall survival in a statistically significant manner. For prediction of 12-month survival status based on these diversity indices, sensitivity and specificity were 0.86 and 0.64, respectively. The area under the receiver operating characteristic curve and the accuracy were 0.76 and 0.75, respectively. For prediction of molecular subtype status, proneural subtype shows highest accuracy of 0.93 among all molecular subtypes based on receiver operating characteristic analysis. We find that measures of spatial diversity from point pattern analysis of intensity habitats from T1 post-contrast and T2 fluid-attenuated inversion recovery images are associated with both tumor subtype status and 12-month survival status and may therefore be useful indicators of patient prognosis, in addition to providing potential guidance for molecularly-targeted therapies in Glioblastoma multiforme.

Purification and characterization of a 7-kDa protein from Clonorchis sinensis adult worms.

A 7-kDa protein was purified from extracts of adult Clonorchis sinensis by a combination of ammonium sulfate precipitation, anion exchange chromatography, cation exchange chromatography, gel-filtration chromatography, and reversed-phase FPLC. The 7-kDa protein exists in the excretory-secretory products of adult C. sinensis, but not in extracts of adult Paragonimus westermani. Also, the 7-kDa protein reacted with the sera of patients with clonorchiasis but not with paragonimiasis or normal human sera. To observe the localization of the 7-kDa protein in the tissue of adult C. sinensis, an immunogold labeling method was followed using anti-7-kDa antibody. The gold particles were observed in the basal layer below the tegumental syncytium, in the interstitial matrix of the parenchyma, and in the content of the uterus. The 7-kDa cDNA was obtained through reverse transcription-polymerase chain reaction using a primer designed from N-terminal sequence analysis. Rapid amplification of cDNA ends (5'-RACE) was used to obtain the complete protein coding sequence. The sequence encodes a 90-amino acid polypeptide. The deduced amino acid sequence of the 7-kDa protein revealed no homology with proteins of different organisms reported so far. These results suggest that the 7-kDa protein is a fluid antigen and may be valuable as a tool for the immunodiagnosis of clonorchiasis.

Expression of cysteine proteinase of Clonorchis sinensis and its use in serodiagnosis of clonorchiasis.

A gene encoding cysteine proteinase from Clonorchis sinensis has been cloned and expressed in Escherichia coli. The cysteine proteinase cDNA fragment was amplified by reverse transcription-polymerase chain reaction using degenerate oligonucleotide primers derived from conserved active site of cysteine proteinases. The 5' and 3' regions of the gene were amplified using rapid amplification of cDNA ends. The cloned gene has an open reading frame of 696 bp and deduced amino acid sequence of 232. Sequence analysis and alignment showed significant homologies with the eukaryotic cysteine proteinases and conservation of the Cys, His, and Asp residues that form the catalytic triad. Analysis of the expressed protein on sodium dodecyl sulfate-polyacrylamide gel electrophoresis showed that the molecular weight of the protein was approximately 28.5 kDa. Proteolytic activity of the expressed protein was inhibited by cysteine proteinase inhibitors such as L-trans-epoxysuccinyl-leucylamide-(4-guanidino)-butane, iodoacetic acid, and leupeptin. The expressed protein showed biochemical properties similar to those of cysteine proteinases of other parasites. The expressed protein strongly reacted with the sera from patients with clonorchiasis but not with the sera from patients with paragonimiasis, fascioliasis, cysticercosis, and sparganosis, or with sera from normal human controls. These results suggest that the expressed protein may be valuable as a specific diagnostic material for the immunodiagnosis of clonorchiasis.

Comparison of analytical and numerical analysis of the reference region model for DCE-MRI.

This study compared three methods for analyzing DCE-MRI data with a reference region (RR) model: a linear least-square fitting with numerical analysis (LLSQ-N), a nonlinear least-square fitting with numerical analysis (NLSQ-N), and an analytical analysis (NLSQ-A). The accuracy and precision of estimating the pharmacokinetic parameter ratios KR and VR, where KR is defined as a ratio between the two volume transfer constants, K(trans,TOI) and K(trans,RR), and VR is the ratio between the two extracellular extravascular volumes, ve,TOI and ve,RR, were assessed using simulations under various signal-to-noise ratios (SNRs) and temporal resolutions (4, 6, 30, and 60s). When no noise was added, the simulations showed that the mean percent error (MPE) for the estimated KR and VR using the LLSQ-N and NLSQ-N methods ranged from 1.2% to 31.6% with various temporal resolutions while the NLSQ-A method maintained a very high accuracy (<1.0×10(-4) %) regardless of the temporal resolution. The simulation also indicated that the LLSQ-N and NLSQ-N methods appear to underestimate the parameter ratios more than the NLSQ-A method. In addition, seven in vivo DCE-MRI datasets from spontaneously occurring canine brain tumors were analyzed with each method. Results for the in vivo study showed that KR (ranging from 0.63 to 3.11) and VR (ranging from 2.82 to 19.16) for the NLSQ-A method were both higher than results for the other two methods (KR ranging from 0.01 to 1.29 and VR ranging from 1.48 to 19.59). A temporal downsampling experiment showed that the averaged percent error for the NLSQ-A method (8.45%) was lower than the other two methods (22.97% for LLSQ-N and 65.02% for NLSQ-N) for KR, and the averaged percent error for the NLSQ-A method (6.33%) was lower than the other two methods (6.57% for LLSQ-N and 13.66% for NLSQ-N) for VR. Using simulations, we showed that the NLSQ-A method can estimate the ratios of pharmacokinetic parameters more accurately and precisely than the NLSQ-N and LLSQ-N methods over various SNRs and temporal resolutions. All simulations were validated with in vivo DCE MRI data.

An analysis of the pharmacokinetic parameter ratios in DCE-MRI using the reference region model.

Dynamic contrast-enhanced magnetic resonance imaging (DCE-MRI) is performed by obtaining sequential MRI images, before, during and after the injection of a contrast agent. It is usually used to observe the exchange of contrast agent between the vascular space and extravascular extracellular space (EES), and provide information about blood volume and microvascular permeability. To estimate the kinetic parameters derived from the pharmacokinetic model, accurate knowledge of the arterial input function (AIF) is very important. However, the AIF is usually unknown, and it remains very difficult to obtain such information noninvasively. In this article, without knowledge of the AIF, we applied a reference region (RR) model to analyze the kinetic parameters. The RR model usually depends on kinetic parameters found in previous studies of a reference region. However, both the assignment of reference region parameters (intersubject variation) and the selection of the reference region itself (intrasubject variation) may confound the results obtained by RR methods. Instead of using literature values for those pharmacokinetic parameters of the reference region, we proposed to use two pharmacokinetic parameter ratios between the tissue of interest (TOI) and the reference region. Specifically, one parameter K(R) is calculated as the ratio between the volume transfer constant K(trans) of the TOI and RR. Similarly, another parameter V(R) is calculated as the ratio between the extravascular extracellular volume fraction v(e) of the TOI and RR. To investigate the consistency of the two ratios, the K(trans) of the RR was varied ranging from 0.1 to 1.0 min(-1), covering the cited literature values. A simulated dataset with different levels of Gaussian noises and an in vivo dataset acquired from five canine brains with spontaneous occurring brain tumors were used to study the proposed ratios. It is shown from both datasets that these ratios are independent of K(trans) of the RR, implying that there is potentially no need to obtain information about literature values from the reference region for future pharmacokinetic modeling and analysis.