Basic concepts and applications of functional magnetic resonance imaging for radiotherapy of prostate cancer.

Recently, the interest to integrate magnetic resonance imaging (MRI) in radiotherapy for prostate cancer has increased considerably. MRI can contribute in all steps of the radiotherapy workflow from diagnosis, staging, and target definition to treatment follow-up. Of particular interest is the ability of MRI to provide a wide range of functional measures. The complexity of MRI as an imaging modality combined with the growing interest of the application to prostate cancer radiotherapy, emphasize the need for dedicated education within the radiation oncology community. In this context, an overview of the most common as well as a few upcoming functional MR imaging techniques is presented: the basic methodology and measurement is described, the link between the functional measures and the underlying biology is established, and finally relevant applications of functional MRI useful for prostate cancer radiotherapy are given.

Introduction

Due to its superior ability to define soft tissue structures, magnetic resonance imaging (MRI) is today the preferred imaging modality for several anatomical locations such as the brain, the vertebral column, the abdomen and the pelvis. This inherent ability makes MRI of interest for direct implementation in radiotherapy. Detailed anatomical description using a uniform standard nomenclature and content such as provided by The Prostate Imaging – Reporting and Data System Version 2 (PI-RADS v.2) [1] in prostate cancer improves consistency in pretreatment staging of patients and thereby patient selection [2]. MRI for target definition and organ at risk identification also has clear advantages over computed tomography (CT) in many cases, especially for the prostate [3], [4], [5]. These advantages has led to the development of MRI-only treatment planning, where the CT-simulation is excluded and the attenuation data are extracted from MR images [6]. Additionally, MRI offers many possibilities for treatment follow-up including adaptive adjustments during the course of the treatment [7].

Beside the tissue contrast and image resolution, of special interest is the ability of MRI to enable more or less quantitative functional information, such as diffusion, perfusion, and MR spectroscopy acting as imaging biomarkers for relevant tissue changes related to cancer. In this review, “functional MRI” is used for all functional MR imaging methods. Historically, the concept of functional MRI was introduced as a description for a special method sensitive to brain activity during stimulation with different paradigms and abbreviated fMRI. In this report functional MRI is used in a broader sense meaning all functional MR methods, abbreviated FMRI (with an uppercase F). The use of FMRI offers possible advantages in the clinical setting. FMRI may be used for identification of intraprostatic lesions for delineation of boost volumes which have been tested in a randomized clinical trial [8]. The advantages of FMRI may also be of use for identification of malignant lymph nodes and possibly for identification of local relapses after prostatectomy.

FMRI techniques address several important tumor specific characteristics, such as cell density, microvessel structure, tissue perfusion, and oxygenation. This information can be utilized in the diagnostic process, for treatment selection, treatment planning and follow-up. Thereby, FMRI may serve as a tool towards individualized radiotherapy [9]. In line with this increased interest and the growing need for FMRI studies, we present a tutorial review for applications of FMRI to radiotherapy of prostate cancer. It includes the major established FMRI techniques, but also some upcoming methods, which may gain clinical interest in the near future.

Method

Each functional method is characterized by a description of the basic measurement methodology and the link between the functional parameter and underlying biology is established. At the end of each section relevant applications of the functional method to prostate radiotherapy are given. The expectations are that the present review will serve as a base of knowledge for the discussions between the radiation oncologist, the radiologist, and the radiotherapy- and MR-physicists.

A search in the database PubMed (Medline) was performed to capture the most relevant functional MRI techniques for radiotherapy. The following search criteria was used:

(prostate AND cancer AND mri AND method) AND ((radiotherapy) OR (radiation therapy)

method was in each search replaced with the functional MRI method of interest.

For each FMRI-method relevant description of the biological interpretation and the applications for radiotherapy have been searched for in the literature. In this case, the search has not been limited to PubMed, but also EMBASE and “gray search” with GOOGLE has been used. From the identified papers or textbooks, the authors have selected the most appropriate sources for the present purpose, i.e. tutorial review.

Results

There was three FMRI methods that dominated the papers published (n = 159), i.e. diffusion weighted imaging, dynamic contrast enhanced MRI and magnetic resonance spectroscopy. The remaining methods had only 25 hits together. There are several limitations of this search strategy, but it serves to point out the three most relevant FMRI methods for radiotherapy applications to prostate cancer. These three methods will therefore have their own chapters, while the less common methods will be handled in one chapter. The major findings are summarized in Table 1 and Fig. 1.

Table 1

The link between the functional measures and the underlying biology.

DiffusionMRI	Tissue oxygenation MRI
cellularity	hypoxia
extracellular volume
membrane permeability	Magnetization transfer MRI
tortuosity	cellularity
microcirculation in capillary network (IVIM)	fibrosis
hypoxia (IVIM)	HP 13C MRI
microstructural organisation of tissue (DTI)	pyruvate glycolysis metabolites
DCE-MRI	lactate
tumor vasculature	alanine
blood vessel permeability	bicarbonate
extravascular extracellular volume
plasma volume
1H MRS
metabolites
citrate
choline
creatine
polyamines

Fig. 1

A graphical illustration on how the three major functional MRI methods can be used in radiotherapy. A plus sign (+) indicate that the method has an advantage for the particular application.

Diffusion MRI

The measured MR signal is sensitive to spatial motion, including both flow and diffusion of water molecules. Diffusion is random free motion due to thermal energy, and the diffusion rate can be described by the diffusion coefficient. In tissue the diffusion of the water molecules are affected by macromolecules, cell membranes and other tissue microstructures that hinder diffusion. Therefore, since the phenomenon is not described by free diffusion molecules, it is denoted apparent diffusion. Several relevant properties of tumor tissue are known to affect the diffusion, e.g. cellularity, extracellular volume fraction, membrane permeability and tortuosity. Compared to normal tissues, tumors have lower diffusion rate due to the increased cellularity and decreased extracellular volume fraction. Diffusion weighted MRI (DWI) is a technique in which the signal intensity in the image is dependent on the degree of diffusion.

The sensitivity of the pulse sequence to diffusion can be regulated by using linear gradients, which encode the spins spatially to diffusion. The spins will precess at different rates depending on location, resulting in dispersion of the phase and signal loss. The same gradient with opposite direction will refocus the spins assuming they have not moved during the time interval between the pulses. Spins that have moved due to diffusion will not rephase completely, since they have experienced a slightly different gradient strength during the two gradient pulses. As a result, regions with low diffusion show high signal in the diffusion weighted MR-image.

The diffusion can also be quantitatively measured. If at least two diffusion weighted images are obtained with different diffusion weighting, often referred to as b-values, the apparent diffusion coefficient (ADC) can be calculated from a monoexponential fit and presented as an ADC-map. It is important to understand that the choice of b-values will result in ADCs, which reflect different aspects of the diffusion (see below) [10].

DWI has evolved to the most important functional method to be added to T2-weighted (T2w) imaging in an examination of the prostate. The sensitivity and specificity for tumor detection increases significantly when T2w imaging is used in combination with DWI [11]. In addition, ADC can also be used for characterization, tumor grading and local tumor staging [12], [13]. In the peripheral zone of the prostate, DWI is the single most important pulse sequence for detection and staging of prostate cancer, and to a large extent determine the PI-RADS v.2 grade together with T2w imaging [1]. Once diagnosed, DWI can be used to improve the definition of the prostate cancer in supplement to T2w images [14]. Since there is some evidence that the degree of restriction of diffusion of prostate cancer is related to Gleason grade, information from ADC could potentially be used in more detail for treatment planning [15]. For treatment techniques, aiming for higher doses (boost) within the prostate, DWI or ADC in a multiparametric setting, can be of guidance to identify the gross target volume for treatment planning [16]. Potentially, the ADC could steer the dose level within the target. Since detailed delineation would be performed using the data on a voxel level, the quality of the ADC measurement needs to be very high, the images geometrically accurate and image artifacts should be insignificant.

MRI diffusion has been used to study treatment response during radiotherapy of prostate cancer [17], [18], [19]. The ADC within tumor tissue has been shown to increase already during the initial phase of the treatment and continues to rise until the completion of therapy months later. These ADC changes seem to occur before changes in levels of PSA [19]. An early increase in ADC reflects increased water mobility through the loss of membrane integrity or an increase in the proportion of total extracellular fluid due to a decrease in cell size or number. The significant difference in ADC between tumor and benign tissue before radiotherapy disappear after completion of the therapy [17]. Diffusion MRI has great potential as biomarker for monitoring treatment response.

The ADC calculated from a monoexponential fit as described above assumes a normal distributed free diffusion. This is a simplification. The diffusion signal deviates from the monoexponential behavior and it becomes apparent for high b-values (>1000 mm2/s). The diffusion kurtosis model comprises the signal behavior by introducing a dedicated kurtosis parameter. It has been hypothesized that the kurtosis parameter represent the interaction of water molecules with cell membranes and intracellular compounds. More irregular and heterogeneous environments would increase the kurtosis [20]. Diffusion kurtosis imaging (DKI) provides an image representing the kurtosis parameter. DKI has been found to differentiate between benign and malignant prostate tissue [21] and also between low grade and high grade tumors [22]. However, a recent systematic review concluded that DKI did not add any value compared to DWI for clinical assessment of prostate cancer [23].The special contrast mechanism may be of interest for monitor treatment response [20]. This was recently indicated for rectal cancer but is yet to be shown for prostate cancer [24].

For low b-values (<200 s/mm2) there is substantial part of the diffusion signal which is affected the by microcirculation of blood in the capillary network. The phenomenon is referred to as intravoxel incoherent motion (IVIM) [10], [25]. Normally, this component is minimized by selecting b-values > 200 s/mm2, or in fact neglected when ADC is measured. However, it is possible to assess both IVIM and ADC, if measurements with a number of lower b-values are performed. There are now indications that the perfusion fraction calculated from the IVIM measurement corresponds to the parameter Ktrans measured by Dynamic Contrast Enhanced MRI (DCE-MRI), specifically when detecting prostate cancer in the transitional zone and locally recurrent disease [26]. In the future IVIM based analyses may provide further insight in diseases and treatment response. In addition, it has recently been found that IVIM in combination with ADC may have a role for hypoxia assessment (see below) [27].

Another special feature of MRI diffusion is diffusion tensor imaging (DTI). The effect from the diffusion on the signal is specific for the direction of the applied diffusion gradient. The ADC map may look very different dependent upon the orientation of the gradient, a result of anisotropic diffusion, i.e. that the diffusion is higher in certain directions. For the prostate, DTI corresponds to the microstructural organization of the tissue and the cancer may change this structure [13]. DTI may also be used to visualize nerve bundles, which may be a tool although not yet clinically used to plan nerve sparing prostatic treatment and monitor post treatment nerve injury and reinnervation [28]. In the future, DTI may contribute to tumor detection or treatment follow-up.

Dynamic contrast enhanced MRI

There are three different techniques to assess perfusion related parameters with MRI. One method, arterial spin labelling (ASL), is based on endogenous contrast from the blood in combination with a special scheme of radiofrequency pulses. The other two methods relies on an i.v. injection of an exogenous contrast agent (CA). Dynamic susceptibility contrast (DSC) MRI is technique that measures the paramagnetic effects in the vessels of the CA during the first pass, while the CA is still enclosed in the intravascular space. ASL and DSC are almost exclusively used for studies of perfusion in the brain, and are not generally used for characterization of tumors in the prostate. The third method, dynamic contrast enhanced MRI (DCE-MRI), is the main method to characterize the vascular properties of tumors in the body including the prostate, and will be described below.

Assessing certain vascular physiological properties is of special interest in cancer since tumor growth creates an environment, which often differs from normal soft tissue with respect to the vascular structure. During tumor growth new vessels are needed, i.e. angiogenesis, but the new vessels will be of different, often poor structure, having a tortuous topology and express a certain leakiness. In short, DCE-MRI can assess tumor vasculature, blood vessel permeability, and extravascular/extracellular volume fraction. Obviously, all these parameters will be of high interest to characterize tumor tissue, as for example with respect to staging and grading [29].

As indicated above, DCE-MRI relies on an i.v. injection of a contrast agent. The standard techniques use clinically available paramagnetic gadolinium (Gd) chelates. The contrast agent affects the relaxation times T1, T2 and T2*. DCE-MRI is entirely built on the T1 shortening effect of the CA. In order to follow the distribution of the CA, i.e. to characterize its pharmacokinetics, a dynamic T1-weighted acquisition is performed.

The most common model for quantitative data analyzing of the time intensity curve is the Tofts model [30]. The analysis results in the parameter Ktrans, which represents the transfer constant from plasma space to tissue space, often referred to as permeability. The analysis also results in the parameters extravascular-extracellular space fraction volume (ve), plasma volume (vp) and the rate constant kep between the ve and vp. The parameter Ktrans is often the parameter of interest for tumor characterization. The rationale behind using this parameter in oncological examinations is that the endothelium in neovascularized cancerous tissue is more prone to leak fluid into the extravascular space than normal vascular structures. The assumption that neovascularization is a key element in the formation and growth of many cancers including prostate cancer has been the basis for DCE-MRI. Although DCE-MRI was introduced more than two decades ago, there is today much more knowledge about the limitations in differentiating prostatitis from cancer in the peripheral zone and benign prostatic hyperplasia from cancer in the transitional zone [31].

The Tofts and similar quantitative models require information such as T1-map of the tissue of interest and an arterial input function (AIF). The T1-map is needed to convert changes of signal to changes of CA concentration and AIF describes the initial concentration of the CA transferred to the tissue of interest and the model. The need for T1-map and AIF make the acquisition protocol complicated. Additionally, the analysis of the data requires substantial computing efforts. Therefore, several semi-quantitative methods have emerged. The most applied method is the initial area under the CA uptake vs time curve (AUC), which basically reflects the initial dynamics of the contrast agent. The AUC often correlates to Ktrans, but the physiologic interpretation is unclear. Still, AUC and similar semi-quantitative techniques are used since the acquisition protocol and the analysis are simplified.

T2w images have similar signal intensities, i.e. low signal, for the malignant tumors as for benign prostatic hyperplasia, prostatitis and post-biopsy hemorrhage. DCE-MRI can therefore contribute to the diagnosis of prostate cancer by increasing the specificity for tumor detection especially for cancer in the peripheral zone of the prostate. T2w imaging with DCE-MRI has shown to be useful for detection recurrent prostate cancers after biochemical failure [32].

DCE-MRI has also a potential role as a method to improve target definition [33]. One benefit of DCE-MRI is that the measurement can often be performed in 3D with relatively high spatial resolution, which enables detailed anatomical information from the prostate. In multiparametric MR imaging (see below) of the prostate, the application of DCE-MRI add information regarding intra-prostatic extension of prostate cancer [34]. The functional information may hypothetically be used to define sub-volumes with different biological properties within the gross tumor volume. The sub-volume may have inferior permeability, with a potential need for different dose strategy based on a threshold, i.e dose sculpting by contours [35]. Biological information on a voxel level may also be used to steer the dose within the target volume, i.e. dose sculpting by numbers or dose painting [36]. Trials with dose painting based on MR imaging is presently scarce, but one example of this concept is represented by the Dutch FLAME study where the intraprostatic boost is delivered based on multiparametric MRI information [37].

DCE-MRI results in new information useful in prognosis as well as assessment of tumor and normal tissue responses to radiation. Use of DCE-MRI in multiparametric MRI of the prostate may play a role in treatment modality selection, target definition, and therapy individualization, although further validation studies are needed [38].

Magnetic resonance spectroscopy

Magnetic resonance spectroscopy (MRS) takes advantage of that spins may resonate at slightly different frequencies, if their molecular environment – chemical bonding – differ. The term chemical shift comes from the difference in resonance frequency relative to a reference compound. A well-known example is the 3.5-ppm chemical shift of fat relative to water in the proton spectrum. The different compounds give rise to different peaks that together result in the MR spectrum. The area under the peak is proportional to the number of nuclei, i.e. the metabolite concentration.

The spectral information does not include spatial localization per se. The spatial localization can be obtained either by using surface coils or by the use of the gradients in the MR-scanner. The latter technique may range from single voxel spectroscopy to chemical shift imaging (CSI). CSI results in 2D/3D arrays of spectra, from which maps of the individual metabolites can be constructed. CSI has been considered a rather time-consuming technique, but new methods are evolving. There is now evidence that acquisition times even with sub ml volume resolution can be in the order of a conventional anatomical imaging sequence [39].

In 1H (proton) MRS, signals from citrate (Cit), choline-containing compounds (Cho), creatine (Cr) and polyamines (PA) are in general detectable in vivo at 1.5 T and higher field strengths as well as mobile lipids. The dominating signals come from water and lipid triglycerides, which therefore need to be suppressed in order to make the metabolites distinguishable. The citrate level are reduced in prostate cancer and can be used to differentiate between tumor and normal prostate tissue [40]. Similarly, the total choline increases in prostate cancer tissue. Generally, there is no difference in total creatine between tumor or normal tissue [40]. The ratio between the metabolites are often used as a readout, since absolute measurements of the metabolites is very difficult. For prostate cancer localization and characterization, the metabolite ratio (Cho + PA + Cr)/Cit has been used as biomarker for malignancies [40].

When 1H MRS is added to the MR imaging protocol the false negatives can be reduced and the specificity can be improved [41]. CSI is both more technically challenging and demanding than DWI and DCE-MRI and the spatial resolution is inherently lower with this technique. Furthermore, most of the reported knowledge around CSI dates several years back during which time other functional imaging techniques, in particular DWI has significantly improved and further developed and being implemented in clinical routine [31], [42]. Due to its complex nature, CSI is not listed as a mandatory functional technique according PI-RADS v.2 guidelines for multiparametric MRI of prostate cancer [1].

The information from CSI can be used to improve the definition of the target [43]. Additionally, spectroscopic data can be used to define regions within the prostate for boost dosing and targets for image guided radiotherapy [44], [45].

MRS has also a role in recurrent cancers. In patients, which have received radiotherapy the radiation induces changes in the tissue properties. It can therefore be difficult to use MRI to distinguish healthy from malignant tissue. The elevated ratio choline to creatine from MRS can still be used to identify residual or recurrent prostate cancer [46]. However, the most important methods for detecting recurrence in irradiated prostate are DWI and DCE-MRI [47].

New and upcoming MRI technologies

Tissue oxygenation MRI

In many tumors, there is an imbalance between oxygen delivery and consumption, leading to hypoxia. The hypoxic environment is known to promote angiogenesis, malignancy, metastases and genetic instability, and to reduce effectiveness of radiation and chemotherapy, and hypoxia is associated with poor prognosis of several cancers [48]. The oxygenation or degree of hypoxia in the tumor is therefore an important factor, both for designing the therapy, as well as an indicator of the progress. However, non-invasive MRI methods to assess the oxygenation level in vivo are scarce.

One way to assess the oxygenation is to use the endogenous contrast mechanism of oxygen on blood. The hemoglobin and deoxyhemoglobin have different magnetic properties and deoxyhemoglobin is strongly paramagnetic. This difference in magnetic susceptibility can be assessed by the T2* relaxation time, a parameter, which reflects susceptibility differences. This phenomenon is known as blood-oxygen-level-dependent (BOLD) contrast. The difference in T2* due to oxygen exposure reflects the deoxyhemoglobin concentration of blood, i.e. the oxygenation. However, the change in T2* will be affected by the blood oxygen saturation, blood volume and hematocrit since these parameters influence the concentration of deoxyhemoglobin of the blood [48]. Therefore, the results may be hard to interpret, unless T2* and T1 is measured simultaneously [49]. There has been indications that the T2* measured in prostate tumors, which were identified on T2w images, correlated to needle oxygen measurements and could be used to identify hypoxic regions [50]. In addition, exposure to 100% oxygen via inhalation during the measurement may facility hypoxia assessment. The change of T1 or T2* due to the oxygen exposure will be related to the initial oxygenation [51]. In prostate tumors, areas with large oxygen enhancement correlate with high perfusion, while no or little enhancement can be attributed to hypoxia, which thereby was possible to identify [52].

A different method to assess hypoxia is based on diffusion. DWI is used to generate separate images of ADC and IVIM, which correlate to the cell density and fractional blood volume (fBV), respectively [27]. The method has recently been validated in prostate tumors, for which areas with low ADC and low fBV represented areas of hypoxia, while oxygenated areas had high ADC and high fBV [27].

Magnetization transfer

In tissues, it is only the protons in the free water pool, which contribute to the MR signal. Protons in macromolecules (proteins, collagen etc) are tightly bound and have very short T2 and are therefore invisible. However, via magnetization transfer, the bound protons can affect the signal of the free protons. Dedicated RF-pulses are used to saturate the bound proton pool, which will exchange magnetization with the free water. This causes a reduced signal from the free water in tissues in which the magnetization transfer (MT) mechanism is prevalent. Since the extent of signal-decay depends on the exchange rate between free and hydration water, MT can be used to provide an alternative tissue contrast in addition to T1, T2, and proton density. The effect, i.e. the measure of the probed macromolecules, is often referred to as the magnetization transfer ratio (MTR), which is the signal without and with the applied RF-pulses.

MT of prostate cancer has not been the subject for many studies. In one study, suspected tumors in the prostate were identified by T2w imaging and the MTR of these regions were measured [53]. Of the patients in which the malignancy of the tumors verified by histopathology, there was a significant higher MTR compared to the patients with negative biopsies. The MTR from negative biopsies also matched values from healthy controls. The increased MTR in tumor is an indication of increased cellularity and an increased amount of high molecular weight material, which is known to correlate with malignity. Recently, it was also found that MTR was significant higher in tumors compared to normal prostate tissue [54]. However, there was no incremental diagnostic value of MTR over conventional DWI.

Hyperpolarized MRI

Normally the signal from tissues is created by the polarization of the nuclei, caused by the high magnetic field in the MR-scanner. However, polarization can also be created by other physical and chemical processes and to levels which are more than 10 000 times higher than what can be achieved by the magnetic field of clinical MR-scanners [55]. This phenomenon, called hyperpolarization (HP) has been used to enhance the signal from gases such as helium-3 or from low concentrations of carbon-13 (13C). The HP is created outside the MR-scanner and outside the body in a dedicated polarizer. Once the hyperpolarized substance leaves the polarizer and is injected in the body, the signal rapidly decays (∼1 min). In many molecules of interest, a carbon atom can be replaced by a 13C isotope, which subsequently is hyperpolarized. After injection of the hyperpolarized substance, it can be followed by MR-imaging or MR-spectroscopy including the metabolites of the molecule. It is worth to note a significant difference between functional studies performed with positron emission tomography (PET) and HP 13C MRI. In PET only the uptake of the injected tracer can be tracked. Using HP 13C MRI, the injected molecule as well as its metabolites can be tracked separately.

It is difficult to find a hyperpolarized molecule for which the decay rate in-vivo is slow enough to allow metabolic studies. However, one such molecule is 13C-pyruvate, a key component in glycolysis. Thereby it enables HP 13C MRI for interesting application in oncology. From injected HP 13C-puruvate, dynamic MR acquisition of the metabolites lactate, alanine and bicarbonate can be performed. Compared to normal cells, tumors have been proven to have a substantially higher uptake as well as higher metabolic activity of pyruvate [55]. The first clinical study of prostate cancer in man using HP 13C-pyruvate was completed recently [56]. The tumor metabolism in 31 patients with prostate cancer was assessed. The level of 13C-lactate/13C-pyruvate was elevated in regions of biopsy-proven cancer. The technique was also able to detect cancer in regions of the prostate that were previously considered tumor–free after examinations with other imaging modalities. However, these type of findings by HP 13C MRI have not yet been proven by histology of resected prostate glands. HP 13C MRI is a new technique, which offers many possibilities. At present there are twelve on-going clinical trials on HP 13C MRI for prostate cancer (December 2018, www.clinicaltrials.gov).

MR-elastography

MR-elastography (MRE) is a method for measurement of the mechanical properties of tissue. This technique requires dedicated hardware, which is commercially available. The mechanical properties can change dramatically due to pathologically processes, such as cancer, fibrosis or inflammation. MR-elastography can assess the tissue properties in large volumes and the measurement results in quantitative information on stiffness.

Mechanical waves, referred to as shear waves, are generated by an external acoustic driver, actuator, for low frequency vibrations. The shear wave displacements in the tissue can be monitored by motion encoded phase contrast MRI. The temporal and spatial characteristics of the wave-field form the basis for an algorithm that transforms the obtained data to a map of the tissue mechanical properties.

The prostate is challenging for MRE. The organ is centrally located in the body and waves from an external actuator are attenuated in the surrounding tissue. An alternative approach is an intra-cavitary endorectal actuator [57]. MRE of prostate has been validated against histopathology with moderate success [58]. One observation was that cancerous tissue was not always stiffer than normal tissue. However, in stiff regions or nodules the localization matched findings by the reference MR images [59]. There are also indications from MRE measurements that radiotherapy makes the gland stiffer [60]. MRE has only been tried out on a limited number of patients with prostate cancer. At present, it seems that MRE of the prostate is inferior to other MR imaging methods, such as T2w, ADC or DCE-MRI.

Multiparametric MRI

Multiparametric MRI (mpMRI) is defined as the combination of an anatomical T1w and T2w imaging and functional imaging techniques. The functional imaging is usually DWI and DCE-MRI, but MRS has also been included in the concept mpMRI. However, in presence of DWI and DCE-MRI, MRS is considered complex and of less added value and therefore often ignored [61]. In addition, DWI is often preferred to DCE-MRI for characterization of lesions, since enhancement in malign lesions can be hard to distinguish from benign conditions, such as prostatitis and prostatic hyperplasia [62].

A lesion is typically depicted as hypointense region on a T2w image, high intensity on DWI (or low ADC), and enhancement on DCE-MRI. Detailed instructions to interpret the images are available in the PI-RADS v.2 document [1], which unfortunately does not provide any recommendations for definition of radiotherapy targets. The mpMRI can be used for tumor detection, staging, assessment of aggressiveness and treatment monitoring [61], [63], [64]. For radiotherapy, mpMRI can have a significant impact on the choice of the treatment. In two studies based on mpMRI findings, 30% of the patients experienced a change in their initial risk group compared to the clinical tumor categorization [65], [66].

After radical prostatectomy or after radiotherapy mpMRI shows a potential to visualize post treatment changes and to detect localized recurrence [47]. For this case, DCE-MRI is of great value and sensitive for detection of small lesions.

Technical aspects

There are some technical aspects on functional MRI imaging of the prostate when the examination is performed in a radiotherapy setting, especially in treatment position. For technical aspects in a standard diagnostic radiology setting, see recent review by Purysko et al. [62]. A flat tabletop in the MR-scanner may create an extra distance between the receiver coil and the body. If the body is kept in a dedicated mould, i.e. treatment position, it can be difficult to use the standard and optimized receiver coils. Both these adjustments may reduce the signal to noise ratio (SNR). At present, there is an urge to create radiotherapy planning workflow based entirely on MRI and exclude the CT. With no density maps from CT, the corresponding information needs to be provided from MR images. There are several methods presented to generate density maps from MR-images [67], and a few are also available commercially [68], [69], and have recently been validated for prostate cancer radiotherapy planning [4], [70]. If the intended use of the examination is to replace the CT and generate density maps from MR-images, it is very important to keep the outer surface of the body intact. One way to accomplish this is to place a dedicated stiff bridge holder between the body and the coil [4], which again may reduce SNR. Altogether, the adjustments to a radiotherapy setting can result in a lower signal to noise ratio (SNR) compared to a standard setting in a radiology. This will affect the quality of functional measurement negatively. For example, the outcome of an ADC measurement is critically dependent on the SNR of the acquired images. In addition, the geometric distortions needs to be monitored and may be of special interest if the functional images are used in for target volume decisions. Typically, this can be an issue for DWI, a method in which MRI sequences sensitive to distortions are often used [71].

For prostate radiotherapy, fiducial markers are often used for positioning the anatomy of the patient accurately for every fraction during the treatment time course [72]. The markers are made of metal in order to create high contrast on X-ray images. Thereby, the markers also create susceptibility artefacts in the MR-images. For normal T1w or T2w imaging this is a minor issue. Certain MRI sequences are very sensitive to susceptibility effects and metallic fiducial markers may prevent reliable measurements of for example diffusion [73]. Other materials such as radiopaque hydrogel are under development to address this issue [74].

In order to obtain quantitative functional data used for example repeated follow-up or adaptive treatment, it can be of extra value to perform the bowel preparation carefully. DWI is very sensitive to gas and to changes of gas filling during the examination. Similarly, even a small improvement on the image quality from using antiperistaltic agents may be considered, when the data should be used for treatment decisions rather than diagnostic procedures [75].

Discussion

Magnetic resonance imaging has played an important role in the diagnosis and management of prostate cancer since it was clinically introduced. Recently, the interest to apply MRI at all stages in the radiotherapy process, i.e. after the diagnostic procedure has been completed, has increased. This was emphasized in an editorial, in which the undeserved needs also were reviewed [76]. One of these needs, i.e. education of the radiation oncology community on MRI and in this case more specifically FMRI, has been addressed in the present review. The recently introduced MR-linac will even further emphasize the need for efficient MR methods in all stages of the radiotherapy workflow [77].

In order to use FMRI methods appropriately and efficiently, it is important to understand both the underlying biological process in the tumor and risk organs during therapy, and how these processes relate to FMRI. In addition, understanding the limitations of what the FMRI techniques are actually measuring, are important for the interpretation of the results. For example, ADCs, initially increase during radiotherapy, as a sign of response, rise to a level above normal tissue due to necrosis, and thereafter return to an intermediate value. Obviously, this indicates that the time point for the measurement, as well as the region of interest for the analysis, is crucial to the interpretation of the results. Additionally, it is reasonable to assume that even in a single tumor there will be an inhomogeneous structure, with a multitude of physiological changes. This increases the complexity and may require more than one functional MRI readout to characterize the tumor. Therefore, multiparametric imaging with combined FMRI techniques will be of especially high value as well as complementary measurements from other imaging modalities. Although, not yet applied to radiotherapy explicitly, the inhomogeneous structure of the tumors has also been successfully addressed with texture analysis of the image data. Several texture features in ADC and mpMRI were able to differentiate between non-cancerous and cancerous tissue including the transition zone and correlate with Gleason score [78], [79].

A well-recognized problem in manual delineation of target structures is the interobserver variability [80], [81]. New methods are being developed for increasing the constancy of tumor delineation [82]. Of special interest is the latest development on artificial intelligence (AI) or machine learning, which may supersede texture analysis in finding other and valuable features for tumor characterization [83]. This could be of particular interest for applications to the MR-linac, when the information, including the functional measures need to be rapidly extracted from the images, while the patient is waiting for the treatment in the machine.

As indicated in this overview several functional MR-methods are available. Although many of the methods have been available in diagnostic radiology for decades, not many methods have been established as a standard tool in clinical routine or even in clinical research trials. There are several explanations. When the functional methods need to be quantitative, they often become complicated and technically challenging both for manufacturers of equipment and for the operators. The methods may require dedicated analyzing software not supplied or available from the vendors. There are limited established standards for data acquisition and data analyzing. However, it is worth to underline that one of the strengths of the functional methods is that they do result in quantitative data. If these functional measures can be brought into robust and validated biomarkers for the imaging community, their working potential is hard to overestimate. More data under different conditions and during different kinds of treatments are needed to reach this level of confidence.

There are several on-going initiatives to standardize FMRI. In Europe there is the European Imaging Biomarkers Alliance [84] within the European Society of Radiology and in the US Radiological Society of North America started the Quantitative Imaging Biomarkers Alliance, QIBA [85]. Within the latter framework, task groups are working on standardization of for example MR-perfusion, diffusion and flow. There are other initiatives on this subject from the American Association of Physicists in Medicine and the National Cancer Institute, although reports are still pending [86].

The different FMRI techniques are of different maturity. The literature on DWI, DCE-MRI and MRS is vast within oncology with applications in most relevant organs. The applications to radiotherapy are less frequent, but steadily increasing. For these FMRI techniques, most vendors provide dedicated sequences as well as tools for analysis. However, there are few standards and the techniques may not be validated.

The remaining techniques for measurement of hypoxia, MTR, hyperpolarization and elastography are more at an experimental level. The measurement techniques varies and the links between read-outs and the biology are less established. An exception may be the method for hypoxia based on DWI [27]. DWI is a method already implemented in most MRI protocols for prostate cancer and therefore hypoxia measurement with DWI may be easier to provide clinically. In addition, new FMRI methods keep emerging, such as the luminal water imaging. A method with a potential to both detect and grade prostate cancer [87].

In parallel to the rapid development of FMRI for radiotherapy, other imaging modalities are gaining interest. Prostate specific PET methods using tracers such 68Ga-PSMA are promising [88]. Future hybrid imaging protocols combining FMRI and prostate specific PET tracers may further increase the possibility to delineate intraprostatic lesions as well as malignant lymph nodes and local relapses. However, the discussion on modalities besides MRI is out of the scope of the present paper.

In conclusion, the role of MR-methods for functional imaging of the prostate as a tool for characterization of tissue and monitoring of response is promising but not established for radiotherapy applications. In this rapidly expanding area of research, as has been outlined in this review, there are many indications that FMRI could play an important role in many of the steps in the radiotherapy process.

Conflict of interest

None of the authors has any conflicts of interest to declare.