Dual Contrast - Magnetic Resonance Fingerprinting (DC-MRF): A Platform for Simultaneous Quantification of Multiple MRI Contrast Agents

Injectable Magnetic Resonance Imaging (MRI) contrast agents have been widely used to provide critical assessments of disease for both clinical and basic science imaging research studies. The scope of available MRI contrast agents has expanded over the years with the emergence of molecular imaging contrast agents specifically targeted to biological markers. Unfortunately, synergistic application of more than a single molecular contrast agent has been limited by MRI’s ability to only dynamically measure a single agent at a time. In this study, a new Dual Contrast - Magnetic Resonance Fingerprinting (DC – MRF) methodology is described that can detect and independently quantify the local concentration of multiple MRI contrast agents following simultaneous administration. This “multi-color” MRI methodology provides the opportunity to monitor multiple molecular species simultaneously and provides a practical, quantitative imaging framework for the eventual clinical translation of molecular imaging contrast agents.

measured dynamically before and after contrast agent administration allowing tracking of the agent's distribution. These studies often incorporate non-targeted contrast agents (e.g., scrambled peptides) as controls to verify the in vivo molecular specificity of the targeting moiety. Due to the difficulty of previously-developed contrast enhanced MRI strategies to uniquely identify two MRI contrast agents administered simultaneously, the targeted and untargeted contrast agents must be studied in separate imaging sessions and likely in separate animal cohorts [15][16][17] . This significant limitation can result in experimental bias due to phenotypic variation. As such, the development of a "multi-color" MRI methodology to independently monitor simultaneously-administered targeted and control MRI contrast agents, and potentially multiple targeted MRI contrast agents, would significantly improve preclinical molecular MRI studies. This "multi-color" MRI capability would also provide a robust pathway for clinical translation of molecular MRI contrast agents by providing the capability to simultaneously compare the biodistribution of a molecular imaging agent with a conventional clinical MRI agent.
The primary limiting factor in the simultaneous assessment of multiple paramagnetic MRI contrast agents is that the unique identification of each agent is challenging. MRI contrast agents directly impact both the T1 and T2 relaxation times according to well-established concentration-dependent linear relationships to their magnetic relaxivities (r 1 and r 2 ) shown in equations (1a) and (1b) below 18 : is the concentration of imaging agent A; T1 0 and T2 0 are the pre-contrast T1 and T2 relaxation times of the tissue; T1 and T2 are the post-contrast T1 and T2 relaxation times; and r 1A and r 2A are the magnetic relaxivities of contrast agent A. Therefore, while an individual MRI contrast agent is typically more sensitive to a particular relaxation parameter (i.e., Gd-chelates for enhancement in T1-weighted imaging acquisitions), each paramagnetic MRI contrast agent still impacts both the T1 and T2 relaxation times. This important factor limits the capability of MRI to independently assess simultaneously-administered contrast agents (e.g., a Gd-based T1 agent and an iron-based T2 agent). The Magnetic Resonance Fingerprinting (MRF) methodology has recently been developed to simultaneously generate inherently co-registered T1 and T2 relaxation time maps in both patients [19][20][21] and animal models 22,23 . MRF uses a unique acquisition and quantification strategy that combines a priori acquisition parameter variation with a dictionary-based pattern matching algorithm to obtain quantitative assessments of multiple imaging parameters simultaneously. Importantly, MRF has been shown to provide quantitative T1 and T2 maps in 10-50 seconds per imaging slice providing the opportunity to dynamically generate quantitative maps of these two important MRI parameters simultaneously [24][25][26][27] . In this study, we demonstrate that the framework for simultaneous T1 and T2 assessments provided by MRF can be used to analytically quantify the local concentration of two different MRI contrast agents present at the same time. Herein, we describe a straightforward multiple contrast agent relaxation model (equations (2a) and (2b) in METHODS) that can be used in combination with the rapid, multi-parametric MRF strategy to independently calculate inherently co-registered concentration maps for two MRI contrast agents. These initial in vitro results represent a proof-of-concept study to: (1) validate the multiple contrast agent relaxation model using a 60 MHz magnetic relaxometer; and (2) demonstrate the application of the rapid MRF method to enable simultaneous calculation of concentration maps for two different paramagnetic MRI contrast agents on a clinical 3 T MRI scanner. Overall, these in vitro studies suggest an imaging framework for future in vivo MRI studies to simultaneously quantify multiple contrast agents. Figure 1 shows the individual magnetic relaxivity plots (R1 and R2 vs concentration) for the gadolinium and manganese contrast agent phantoms obtained from the 60 MHz relaxometer using inversion recovery spin-echo acquisitions for T1 relaxation time measurements followed sequentially by a Carr-Purcell-Meiboom-Gill (CPMG) acquisition for T2 relaxation time measurements. A list of in vitro phantoms and their respective concentrations of contrast agents is shown in Supplementary Table S1. The T1 measurements were repeated twice and the T2 assessments were repeated three times to ensure consistency. The relaxivity data resulted in significant linear correlations as expected from equations (1a) and (1b) (R 2 > 0.993, p < 0.0001). The 60 MHz magnetic relaxivities at 37 °C for the gadolinium agent (slope of the linear correlation lines) were 0.0040 mM −1 ms −1 (r 1 ) and 0.0048 mM −1 ms −1 (r 2 ). The corresponding magnetic relaxivities for the manganese agent were 0.0054 mM −1 ms −1 (r 1 ) and 0.0652 mM −1 ms −1 (r 2 ). The magnetic relaxivity values for the 60 MHz relaxometer, as well as the relaxivities measured at 3 T, are shown in Table 1. The mean T1 and T2 relaxation times for the deionized water phantom with no contrast agent (deionized water alone) were 4250 ms (T1 0 ) and 2760 ms (T2 0 ), respectively. NMR relaxivity measurements are in reasonable agreement with previously reported results for both the Gd 28 and Mn 29,30 contrast agents used here.

MHz Relaxometry: Multiple Contrast Agent Relaxation Model Validation.
Equations (3a) and (3b) (in METHODS) were then used to calculate the gadolinium and manganese concentrations for all of the phantoms (n = 17) using the measured T1 and T2 values for each in vitro sample, the calculated magnetic relaxivities for each contrast agent, and the "non-contrast" T1 0 and T2 0 values for the deionized water samples. Plots of the relaxometry-based concentration estimates against the known concentration in each phantom for both gadolinium and manganese chloride are shown in Fig. 2. Concentration estimates obtained from the 60 MHz relaxometry data shown in Fig. 2 were calculated for both the six samples containing both gadolinium and manganese contrast agents (Supplementary Table S1 Phantoms 11-16) as well as the ten samples containing either gadolinium or manganese used in the relaxometry analysis (Fig. 1, Supplementary Table S1 Phantoms 1-10). One additional sample containing only deionized water was also scanned (Supplementary Table S1 Phantom 17). Pure samples with a single contrast agent (either Gd or Mn) were analyzed to increase the number of concentrations calculated and to verify that the method returned a value of 0 if the agent was not present in solution. This comparison resulted in a significant linear correlation (Pearson Correlation: R 2 > 0.998, two-tailed probability p < 0.0001) for both the gadolinium and manganese contrast agents. Importantly, the slopes of these correlations are both near unity (1.003 (Gd) and 0.975 (Mn), respectively) suggesting that the Relaxivity assessments for (a) gadolinium (Gd) and (b) manganese (Mn) contrast agents from a 60 MHz relaxometer using phantoms containing varying concentrations of a single contrast agent. Slopes of the fitted lines of R1, R2 vs. agent concentration (n = 6 for each agent) were used to determine the relaxivities (r 1 and r 2 ) of the two agents. Pearson correlations resulted in significant correlations of concentration vs. R1 and R2 for both contrast agents (R 2 > 0.993, two-tailed probability p < 0.0001).   Table S1 Phantoms 11-16) are shown as a subset of the data separately ( Supplementary Fig. S1). These initial relaxometric results demonstrate that the multiple contrast agent relaxation model shown in equations (2a) and (2b) is capable of providing accurate estimates for the concentration of gadolinium and manganese-based contrast agents. Importantly, these results appear to be consistent whether there is 0, 1, or 2 contrast agents in the phantom.

T DC-MRF: Simultaneous Assessment of Two Paramagnetic MRI Contrast Agents. Similar
to the relaxometry results above, T1 and T2 relaxation time assessments for gadolinium and manganese containing phantoms (n = 17) were obtained using the MRF method on a clinical 3 T MRI scanner at room temperature. Phantoms scanned contained the same contrast agent concentrations as for the NMR relaxometer experiments described above. This data was acquired using a FISP-MRF acquisition repeated 12 times following repositioning to measure average T1 and T2 values for each phantom. In contrast to the relaxometer studies above, the MRF-based T1 and T2 measurements were obtained simultaneously for all phantoms (n = 17). Figure 3 shows the individual magnetic relaxivity plots (R1 and R2 vs concentration) for the gadolinium and manganese contrast agents obtained from the MRF data averaged over the 12 repeats. The relaxivity data resulted in a significant linear correlation for both the gadolinium agent (Pearson Correlation: R 2 ≥ 0.997, two-tailed probability p < 0.0001; r 1 = 0.0056 mM −1 ms −1 ; r 2 = 0.0076 mM −1 ms −1 ) and the manganese agent (R 2 ≥ 0.999, two-tailed probability p < 0.0001; r 1 = 0.0067 mM −1 ms −1 ; r 2 = 0.1144 mM −1 ms −1 ). These MRF-based relaxivity values compared favorably to relaxivity values obtained from conventional inversion recovery and single-echo spin echo MRI experiments for both gadolinium (r 1 = 0.0051 mM −1 ms −1 ; r 2 = 0.0060 mM −1 ms −1 ) and manganese (r 1 = 0.0068 mM −1 ms −1 ; r 2 = 0.1079 mM −1 ms −1 ). Relaxivity results at 60 MHz and 3 T are summarized in Table 1. These relaxivity results also were in reasonable agreement with literature values 28,29,31 . DC-MRF was also used to obtain T1 and T2 relaxation times for the deionized water sample with no contrast agent (T1 0 = 2897 ms; T2 0 = 946 ms).
Representative DC-MRF maps of estimated gadolinium and manganese concentration calculated on a pixel-by-pixel basis from MRF-based T1 and T2 maps using equations (3a) and (3b) are shown in Fig. 4a for both the phantoms containing only gadolinium contrast agent (Supplementary Table S1 Phantoms 1-5; n = 5), only manganese contrast agent (Supplementary Table S1 Phantoms 6-10; n = 5), both gadolinium and manganese agents (Supplementary Table S1 Phantoms 11-16; n = 6), or solvent alone (Supplementary Table S1 Phantom 17; deionized water, n = 1). Theoretical maps of the known phantom concentrations are shown for comparison in Fig. 4b. The estimated gadolinium and manganese concentration maps are visually consistent with the known concentrations. A quantitative comparison of the mean DC-MRF concentration estimates from the region-of-interest analysis with the known concentrations are shown in Fig. 5 and resulted in significant correlations (Pearson Correlations: Gd: R 2 = 0.9987, two-tailed probability p < 0.0001; Mn: R 2 = 0.998, p < 0.0001). The Slopes of the fitted lines of R1, R2 vs. agent concentration (n = 6 for each agent) were used to determine the relaxivities (r 1 and r 2 ) of the two agents. Pearson correlations resulted in significant correlations of concentration vs. R1 and R2 for the two contrast agents (R 2 ≥ 0.997, two-tailed probability p < 0.0001).
mean DC-MRF concentrations were also in reasonable agreement with actual values as evidenced by the slopes of the correlation lines equal to 0.988 (Gd) and 0.980 (Mn), respectively. Results from the phantoms containing mixtures of both contrast agents are presented separately (Supplementary Fig. S2).
Mean and standard deviations of the gadolinium and manganese concentrations were calculated from the MRF-based concentration maps for the 6 phantoms containing both agents using an ROI analysis. The DC-MRF method resulted in significant differences among all samples for both the gadolinium concentration estimates (two-tailed unpaired Student's t-test, p < 0.01) and the manganese concentration estimates (two-tailed unpaired Student's t-test, p < 0.001). Overall, these results suggest that the DC-MRF methodology provides accurate and precise assessments of two paramagnetic MRI contrast agents at the same time.

Discussion
In this initial report, we have demonstrated the capability of the Dual Contrast (DC)-MRF method to simultaneously measure the concentration of two paramagnetic MRI contrast agents. Herein, we present initial in vitro relaxometry measurements at 60 MHz to evaluate the proposed multiple contrast agent relaxation model for two paramagnetic MRI contrast agents (Equations (2a) and (2b)). We also show initial in vitro DC-MRF results on a clinical 3 T MRI scanner demonstrating the capability of DC-MRF to independently quantify the local concentration of two paramagnetic MRI contrast agents using simultaneously-measured T1 and T2 relaxation times. Overall, DC-MRF provides an imaging platform that can be used to independently monitor multiple paramagnetic MRI contrast agents with numerous clinical and preclinical molecular imaging applications.
We first validated the multiple contrast agent relaxation model described in equations (2a) and (2b). The significant correlations between the known phantom concentrations and the 60 MHz relaxometric estimates (R 2 > 0.998, Fig. 2) as well as the 3 T DC-MRF estimates (R 2 > 0.998, Figs. 4 and 5) suggests that this model is accurate for the concentration ranges of the two specific MRI contrast agents (Multihance ® , gadobenate dimeglumine and MnCl 2 ) used in this study. These results also demonstrate the ability of the model to accurately measure contrast agent concentrations at multiple field strengths. This straightforward linear relaxation model assumes that the two contrast agents impact the overall T1 and T2 relaxation times independently with minimal interactions between the two contrast agents. While this model may be expected to be valid for moderate agent concentrations, high local concentrations of one or both of the contrast agents could result in deviations from the linear model as the agents compete for interactions with the surrounding water molecules. Implicit in this model also is the assumption that the relaxivities (r 1 and r 2 ) are distinctly different for the two MRI contrast agents. Therefore, this relaxation model may have limitations if used to detect two MRI contrast agents with similar r 1 and r 2 values (e.g., gadopentetate dimeglumine and gadolinium-diethylenetriamine pentaacetic acid (Gd-DTPA)). Furthermore, multiple follow-on studies (in vitro and in vivo) will be needed to thoroughly explore the limitations of the proposed relaxation model. Regardless, these initial results demonstrate that the multiple contrast agent relaxation model can provide an analytical basis to determine the concentration of two different paramagnetic MRI contrast agents.
A key advantage of DC-MRF is that it provides the opportunity to simultaneously and dynamically detect multiple MRI contrast agents. Prior studies have attempted to detect multiple MRI contrast agents in a single scanning session using sequential administration of contrast agents 32 , ratiometric methods to detect the presence of activatable MRI agents 14 , and machine learning 33 . Other groups have utilized chemical exchange saturation transfer (CEST) MRI techniques which can have reduced sensitivity on low-field (≤ 3 T) MRI scanners [34][35][36] . It is important to note that these prior studies primarily incorporated sequential MRI assessments to separately assess the different MRI contrast agents. Sequential T1 and T2 measurements (or any other method that results in both T1 and T2 maps) could be used with equations (3a) and (3b) in place of MRF-based T1 and T2 assessments under the assumption that the concentration of the agent is not appreciably changing during the measurement time. Difficulty with this strategy arises because this assumption may or may not be valid based on the disease state and agents used. The ability of the MRF data to be acquired in as little as 10 seconds provides the opportunity to dynamically assess a wide variety of contrast agents regardless of the pharmacokinetics making it a more general solution with fewer required assumptions. Further, prior studies have shown that MRF is more temporally efficient than other rapid MRI techniques 37,38 , and has been implemented on both clinical and preclinical MRI scanners. An additional benefit of the simultaneous measurement of T1 and T2 provided by MRF is these two relaxation time maps always being co-registered regardless of subject motion. This important feature allows the pixelwise concentration maps shown in Fig. 4b to be calculated without any mismatch errors or utilization of additional co-registration methodology. Therefore, DC-MRF may provide an adaptable, quantitative imaging framework to assess two MRI contrast agents simultaneously for a wide variety of imaging applications.
While these initial in vitro results suggest that DC-MRF can provide accurate assessments of two paramagnetic MRI contrast agents, in vivo imaging studies will be required to more fully evaluate the DC-MRF methodology. A schematic for an eventual in vivo DC-MRF implementation is shown in Supplementary Fig. S3. These in vivo DC-MRF experiments could be conducted similarly to conventional dynamic contrast enhanced MRI studies where images are collected before, during, and after administration of a contrast agent [39][40][41][42] . For in vivo DC-MRF, a pre-contrast MRF scan is performed to obtain baseline T1 and T2 maps (T1 0 and T2 0 in equations (2a) and (2b)). Both MRI contrast agents would then be administered simultaneously as a mixture (pseudo-color orange in syringe and tumor in Supplementary Fig. S3) during dynamic acquisition of MRF-based T1 and T2 maps. These inherently co-registered T1 and T2 maps would be then used with equations (3a) and (3b) to calculate in vivo concentration maps for each contrast agent (Agent A, red; Agent B, yellow; Supplementary Fig. S3).
Importantly, there are several challenges that must be overcome to utilize the DC-MRF methodology for in vivo experiments. First, MRF has been previously shown to be sensitive to inhomogeneities in both the B1 and B0 fields. As variation in B1 and B0 are expected to significantly increase for in vivo experiments, accurate B0 shimming and/ or B1 corrections will be needed to avoid significant errors in T1 and T2 measurements 43 as well as the calculated contrast agent concentrations. Additionally, in vivo relaxivities will likely be different from the in vitro relaxivities shown here due to complex molecular interactions experienced by the contrast agents in vivo. This may be particularly complicated as these interactions can vary considerably between normal tissues and pathologies. Addressing this important limitation will require careful in vivo relaxivity assessments and potentially validation using elemental analysis of excised tissues. The pre-contrast T1 and T2 relaxation times of the tissues and pathologies of interest may also pose specific challenges. For example, if the tissue of interest has a low T2 relaxation time before contrast administration (e.g., liver), resolving the concentration of the two agents may be problematic due to excessive T2 decay. In addition, in vivo studies may be impacted by multiple factors including partial volume effects, flow effect, and magnetization transfer and/or chemical exchange that may benefit from multi-exponential relaxation models instead of the simple mono-exponential models used here 44,45 . Despite these challenges, DC-MRF offers a unique opportunity to expand the portfolio of in vivo contrast-enhanced MRI applications.
In addition to the need for follow-on in vivo studies, the DC-MRF results presented in this initial report have multiple areas for future exploration. As described above, the multiple contrast agent relaxation model described herein provides an analytical solution to independently quantify two MRI contrast agents. It is conceivable that DC-MRF could be applied for three or more contrast agents. However, in that case, equations (2a) and (2b) become underdetermined. As such, using the DC-MRF approach to differentiate three or more agents would require additional measurements and/or alternative numerical or machine learning 33 approaches to resolve these agent concentrations accurately. As described above, another limitation of the DC-MRF methodology is the requirement for two agents with different relaxivities. In this initial proof-of-concept study, we used a clinical gadolinium agent and a manganese chloride based agent with substantial differences in r 1 and r 2 (Figs. 1 and 3). The minimum differences in relaxivity needed to reliably differentiate two contrast agents remains undetermined. Reduced relaxivity differences would also likely limit the ability of the DC-MRF methodology to detect small concentration changes. Additionally, these relaxivities are also known to change as a function of magnetic field strength. Therefore, in vivo validation studies would likely be required for each MRI field strength and for a variety of MRI contrast agents to explore the limitations of DC-MRF.
Although DC-MRF may have numerous applications, this new methodology is particularly well suited to support the development and eventual clinical translation of molecular MRI contrast agents. As described above, DC-MRF would provide the opportunity to detect both a molecularly-targeted contrast agent and a control non-targeted contrast agent at the same time in the same subject. The only constraint in achieving these simultaneous assessments using DC-MRF is that the two contrast agents must have different relaxivities. This constraint may require that the targeted and control contrast agents incorporate different lanthanides (e.g., gadolinium and dysprosium) in order to produce differential relaxivities while also retaining similar pharmacokinetic properties. A similar approach could also be used to allow multiple molecular imaging targets to be assessed simultaneously. For example, two targeted MRI contrast agents could be used to simultaneously track drug delivery and therapeutic impact in cancer treatment (e.g., a therapeutic agent 46 labeled with gadolinium, and a second agent targeted to apoptosis 47 labeled with dysprosium). Importantly, the DC-MRF methodology may also directly aid in the clinical translation of molecular imaging agents by allowing for the direct comparison between a conventional, non-targeted clinical MRI contrast agent (e.g., Multihance ® ) with a molecular imaging agent in a single patient. In this way, DC-MRF would provide the opportunity to efficiently establish the molecular specificity of the molecular contrast agent in heterogeneous human diseases.
In conclusion, we describe a new Dual Contrast -Magnetic Resonance Fingerprinting technique that can be used to independently quantify the local concentration of two paramagnetic MRI contrast agents administered simultaneously. These initial in vitro results validate the proposed multiple contrast agent relaxation model and demonstrate a new quantitative imaging methodology that can be used to simultaneously generate concentration maps for two different MRI contrast agents. Overall, these results suggest a new application for the MRF technology to provide quantitative assessments that lays the foundation for numerous clinical and preclinical multi-agent imaging applications.

Methods
This section provides details on the validation of the multiple contrast agent relaxation model and the MRF acquisition used in combination to simultaneously estimate the concentration of two paramagnetic MRI contrast agents.
Multiple Contrast Agent Relaxation Model. As described in the INTRODUCTION section, a single paramagnetic MRI contrast agent exhibits concentration-dependent T1 and T2 relaxation effects as described by equations (1a) and (1b), respectively. Herein, we are proposing a straightforward linear model to incorporate a second MRI contrast agent B as shown in equations (2a) and (2b): To test the multiple contrast agent relaxation model in equations (2a) and (2b), each of the samples were transferred to 5 mm NMR tubes (Norell, 507-HP-7) and scanned on a Bruker Minispec 60 MHz relaxometer (Bruker Biospin, Billerica, MA). The relaxometer was used to obtain T1 and T2 relaxation time assessments for each sample using an inversion recovery spin echo technique (T1: 7 inversion times) and a multi-echo Carr-Purcell-Meiboom-Gill (CPMG) MRI acquisition (1,000-10,0000 echoes) 48 . All relaxometry experiments were conducted at 37 °C. The relaxometric measurements were repeated twice for T1 measurements and three times for T2 measurements to ensure accurate assessments. These repeated measures were averaged to obtain single T1 and T2 values for each in vitro phantom.
The resulting T1 and T2 values for each phantom were then used to calculate the agent concentration in each tube in a multi-step process. First, the T1 and T2 values for the solvent (deionized water) alone were measured and established as "pre-contrast" relaxation times (T1 0 and T2 0 ). Second, the phantoms containing the individual agents were analyzed to calculate the magnetic relaxivities of both agents (r 1G , r 2G , r 1M , r 2M ) through a linear least-squares fit to the plot of R1 (1/T1 in ms −1 ; Equation (1a)) and R2 (1/T2 in ms −1 ; Equation (1b)) as a function of gadolinium or manganese concentration (in mM) using established methods. The slopes of the resulting fits were used as the magnetic relaxivities for each agent. From these relaxivity results, as well as the "pre-contrast" T1 0 and T2 0 values, the T1 and T2 relaxation times for each sample were used to analytically calculate the gadolinium and manganese contrast agent concentrations for each sample using equations (3a) and (3b). These equations were derived directly from the algebraic solution to equations (2a) and (2b) (with Gd = A and Mn = B): where ΔR1 = 1/T1-1/T1 0 , ΔR2 = 1/T2-1/T2 0 . The agent concentrations calculated from equations (3a) and (3b) were then compared with known concentrations in each phantom using Pearson correlations with a probability of p < 0.05 used as a determination for significance.

In Vitro DC-MRF Assessments at 3 T.
To determine the capability of the DC-MRF technique to provide quantitative imaging based assessments of contrast agent concentration, we obtained MRF-based T1 and T2 maps of the gadolinium-containing and manganese-containing phantoms evaluated in the relaxometry studies described above. The solutions in these phantoms were taken from the same source solutions as the NMR studies and placed into 15 mL centrifuge tubes (Fisher Scientific, S50712). All MRF studies were conducted on a Siemens Skyra 3 T MRI scanner (Siemens Healthineers, Erlangen, Germany). The MRF acquisition (Siemens Work-In-Progress #881v23) utilized a FISP acquisition kernel designed with a priori variation in both flip angle (FA) and repetition time (TR, baseline TR of 12 ms) to generate MRF signal evolution profiles sensitive to both T1 and T2 relaxation times 25 . The MRF method acquired 3000 images with time-varying contrast generated by the FA and TR variation. This FISP-MRF implementation included a non-selective inversion preparation (inversion time = 21 ms) immediately prior to the FISP-MRF image acquisitions to increase the sensitivity of the MRF signal evolution profiles to T1 relaxation times 22 . The MRF acquisition also incorporated undersampled spiral trajectories as described previously 25 . The acquisition time of one slice was 47 seconds with a FOV of 380 × 380 mm, an image matrix of 352 × 352, and a slice thickness of 5 mm. Measurement of the excitation (B1) field was incorporated to mitigate the effects of inhomogeneous B1 field on the T1 and T2 relaxation time estimates 43 . A fundamental component of the MRF reconstruction is the development of a large dictionary of signal evolution profiles that are subsequently "matched" to the acquired MRF signal evolution profile for each imaging voxel using vector-based inner product comparisons. The MRF dictionary was created as described previously 24 assuming a mono-exponential relaxation model. MRF image reconstruction as well as subsequent dictionary matching were performed on the Siemens Skyra 3 T MRI scanner. Generation of the quantitative MRF-based T1 and T2 relaxation time maps was attained by matching the acquired MRF profiles on a pixel-by-pixel basis to the MRF dictionary of simulated profiles from all logical combinations of T1 (10-100 ms, increment = 10 ms; 100-1000 ms, increment = 20 ms; 1000-2000 ms, increment = 40 ms; 2000-4500 ms, increment = 100 ms) and T2 (2-100 ms, increment = 2 ms; 100-150 ms, increment = 5 ms; 160-300 ms, increment = 10 ms; 300-800 ms, increment = 50 ms; 800-1600 ms, increment = 100 ms; 1600-3000 ms, increment = 200 ms). The T1 and T2 maps were exported for further offline processing in MATLAB (MathWorks, Natick, MA).
Mean T1 and T2 values for each phantom were obtained from the MRF maps using a region of interest (ROI) analysis. Similar to the relaxometric studies above, mean MRF-based T1 and T2 values for the sample with no contrast agent (deionized water only) was used as a measure of T1 0 and T2 0 , respectively. The mean T1 and T2 relaxation times from the phantoms containing only a single agent were used to estimate the magnetic relaxivities for the gadolinium and manganese contrast agents. The T1 and T2 maps from the MRF acquisition were then used to calculate gadolinium and manganese concentration maps using the calculated 3 T relaxivities and equations (3a) and (3b), respectively. An ROI analysis of the maps was then used to calculate a mean gadolinium and manganese concentration value for each phantom. The mean DC-MRF concentration estimates were compared with known values using Pearson correlations. The MRF acquisition was repeated 12 times allowing for sample repositioning as well as scanner adjustments to test the capability of the DC-MRF method to statistically differentiate samples with both contrast agents.
Magnetic relaxivities were also obtained using conventional MRI assessments for comparison with the MRF-based relaxivity assessments. T1 relaxation time assessments were obtained with an inversion recovery spin echo acquisition (8 inversion times), and T2 relaxation times were obtained with a single-echo spin echo acquisition (8 echo times). Fitting was performed with the appropriate mono-exponential model generating T1 and T2 maps. Equations (1a) and (1b) were then used to calculate r 1 and r 2 for the gadolinium and manganese contrast agents, respectively.
SCiENtiFiC REPoRtS | 7: 8431 | DOI:10.1038/s41598-017-08762-9 Statistical Analysis. Data were compared using Pearson correlation and unpaired, two-tailed Student's t-tests when appropriate. In this study, p < 0.05 was used to establish statistical significance and the p-values are reported as less than the largest reasonable, round number. Repeated studies were performed on the same days after repositioning of the specimen and performing scanner adjustments including magnetic field shimming, radiofrequency excitation calibration, and receiver gain adjustments. Data Availability. The data and protocols used in the current study are available from the corresponding author on reasonable request.