Development of a Secure Web-Based Medical Imaging Analysis Platform: The AWESOMME Project.

Precision medicine research benefits from machine learning in the creation of robust models adapted to the processing of patient data. This applies both to pathology identification in images, i.e., annotation or segmentation, and to computer-aided diagnostic for classification or prediction. It comes with the strong need to exploit and visualize large volumes of images and associated medical data. The work carried out in this paper follows on from a main case study piloted in a cancer center. It proposes an analysis pipeline for patients with osteosarcoma through segmentation, feature extraction and application of a deep learning model to predict response to treatment. The main aim of the AWESOMME project is to leverage this work and implement the pipeline on an easy-to-access, secure web platform. The proposed WEB application is based on a three-component architecture: a data server, a heavy computation and authentication server and a medical imaging web-framework with a user interface. These existing components have been enhanced to meet the needs of security and traceability for the continuous production of expert data. It innovates by covering all steps of medical imaging processing (visualization and segmentation, feature extraction and aided diagnostic) and enables the test and use of machine learning models. The infrastructure is operational, deployed in internal production and is currently being installed in the hospital environment. The extension of the case study and user feedback enabled us to fine-tune functionalities and proved that AWESOMME is a modular solution capable to analyze medical data and share research algorithms with in-house clinicians.

Diffusion-Weighted MRI of the Liver in Patients With Chronic Liver Disease: A Comparative Study Between Different Fitting Approaches and Diffusion Models.

BACKGROUND: Diffusion-weighted imaging (DWI) has been considered for chronic liver disease (CLD) characterization. Grading of liver fibrosis is important for disease management. PURPOSE: To investigate the relationship between DWI's parameters and CLD-related features (particularly regarding fibrosis assessment). STUDY TYPE: Retrospective. SUBJECTS: Eighty-five patients with CLD (age: 47.9 ± 15.5, 42.4% females). FIELD STRENGTH/SEQUENCE: 3-T, spin echo-echo planar imaging (SE-EPI) with 12 b-values (0-800 s/mm2 ). ASSESSMENT: Several models statistical models, stretched exponential model, and intravoxel incoherent motion were simulated. The corresponding parameters (Ds , σ, DDC, α, f, D, D*) were estimated on simulation and in vivo data using the nonlinear least squares (NLS), segmented NLS, and Bayesian methods. The fitting accuracy was analyzed on simulated Rician noised DWI. In vivo, the parameters were averaged from five central slices entire liver to compare correlations with histological features (inflammation, fibrosis, and steatosis). Then, the differences between mild (F0-F2) or severe (F3-F6) groups were compared respecting to statistics and classification. A total of 75.3% of patients used to build various classifiers (stratified split strategy and 10-folders cross-validation) and the remaining for testing. STATISTICAL TESTS: Mean squared error, mean average percentage error, spearman correlation, Mann-Whitney U-test, receiver operating characteristic (ROC) curve, area under ROC curve (AUC), sensitivity, specificity, accuracy, precision. A P-value <0.05 was considered statistically significant. RESULTS: In simulation, the Bayesian method provided the most accurate parameters. In vivo, the highest negative significant correlation (Ds , steatosis: r = -0.46, D*, fibrosis: r = -0.24) and significant differences (Ds , σ, D*, f) were observed for Bayesian fitted parameters. Fibrosis classification was performed with an AUC of 0.92 (0.91 sensitivity and 0.70 specificity) with the aforementioned diffusion parameters based on the decision tree method. DATA CONCLUSION: These results indicate that Bayesian fitted parameters may provide a noninvasive evaluation of fibrosis with decision tree. EVIDENCE LEVEL: 1 TECHNICAL EFFICACY: Stage 1.

Optimal control in a magnetization-prepared rapid acquisition gradient-echo sequence.

This article proposes a numerical framework to determine the optimal magnetization preparation in a three-dimensional magnetization-prepared rapid gradient-echo (MP-RAGE) sequence to obtain the best achievable contrast between target tissues based on differences in their relaxation times. The benefit lies in the adaptation of the algorithm of optimal control, GRAdient Ascent Pulse Engineering (GRAPE), to the optimization of magnetization preparation in a cyclic sequence without full recovery between each cycle. This numerical approach optimizes magnetization preparation of an arbitrary number of radio frequency pulses to enhance contrast, taking into account the establishment of a steady state in the longitudinal component of the magnetization. The optimal control preparation offers an optimized mixed T 1 / T 2 contrast in this traditional T 1 -weighted sequence. To show the versatility of the proposed method, numerical and in vitro results are described. Examples of contrasts acquired on brain regions of a healthy volunteer are presented for potential applications at 3 T.

Alteration of skeletal muscle energy metabolism assessed by 31 P MRS in clinical routine: Part 2. Clinical application.

In this second part of a two-part paper, we intend to demonstrate the impact of the previously proposed advanced quality control pipeline. To understand its benefit and challenge the proposed methodology in a real scenario, we chose to compare the outcome when applying it to the analysis of two patient populations with significant but highly different types of fatigue: COVID-19 and multiple sclerosis (MS). 31 P-MRS was performed on a 3 T clinical MRI, in 19 COVID-19 patients, 38 MS patients, and 40 matched healthy controls. Dynamic acquisitions using an MR-compatible ergometer ran over a rest (40 s), exercise (2 min), and a recovery phase (6 min). Long and short TR acquisitions were also made at rest for T1 correction. The advanced data quality control pipeline presented in Part 1 is applied to the selected patient cohorts to investigate its impact on clinical outcomes. We first used power and sample size analysis to estimate objectively the impact of adding the quality control score (QCS). Then, comparisons between patients and healthy control groups using the validated QCS were performed using unpaired t tests or Mann-Whitney tests (p < 0.05). The application of the QCS resulted in increased statistical power, changed the values of several outcome measures, and reduced variability (standard deviation). A significant difference was found between the T1PCr and T1Pi values of MS patients and healthy controls. Furthermore, the use of a fixed correction factor led to systematically higher estimated concentrations of PCr and Pi than when using individually corrected factors. We observed significant differences between the two patient populations and healthy controls for resting [PCr]-MS only, [Pi ], [ADP], [H2 PO4 - ], and pH-COVID-19 only, and post-exercise [PCr], [Pi ], and [H2 PO4 - ]-MS only. The dynamic indicators τPCr , τPi , ViPCr , and Vmax were reduced for COVID-19 and MS patients compared with controls. Our results show that QCS in dynamic 31 P-MRS studies results in smaller data variability and therefore impacts study sample size and power. Although QCS resulted in discarded data and therefore reduced the acceptable data and subject numbers, this rigorous and unbiased approach allowed for proper assessment of muscle metabolites and metabolism in patient populations. The outcomes include an increased metabolite T1 , which directly affects the T1 correction factor applied to the amplitudes of the metabolite, and a prolonged τPCr , indicating reduced muscle oxidative capacity for patients with MS and COVID-19.

Alteration of skeletal muscle energy metabolism assessed by phosphorus-31 magnetic resonance spectroscopy in clinical routine, part 1: Advanced quality control pipeline.

Implementing a standardized phosphorus-31 magnetic resonance spectroscopy (31 P-MRS) dynamic acquisition protocol to evaluate skeletal muscle energy metabolism and monitor muscle fatigability, while being compatible with various longitudinal clinical studies on diversified patient cohorts, requires a high level of technicality and expertise. Furthermore, processing data to obtain reliable results also demands a great degree of expertise from the operator. In this two-part article, we present an advanced quality control approach for data acquired using a dynamic 31 P-MRS protocol. The aim is to provide decision support to the operator to assist in data processing and obtain reliable results based on objective criteria. We present here, in part 1, an advanced data quality control (QC) approach of a dynamic 31 P-MRS protocol. Part 2 is an impact study that will demonstrate the added value of the QC approach to explore data derived from two clinical populations that experience significant fatigue, patients with coronavirus disease 2019 and multiple sclerosis. In part 1, 31 P-MRS was performed using 3-T clinical MRI in 175 subjects from clinical and healthy control populations conducted in a University Hospital. An advanced data QC score (QCS) was developed using multiple objective criteria. The criteria were based on current recommendations from the literature enriched by new proposals based on clinical experience. The QCS was designed to indicate valid and corrupt data and guide necessary objective data editing to extract as much valid physiological data as possible. Dynamic acquisitions using an MR-compatible ergometer ran over a rest (40 s), exercise (2 min), and a recovery phase (6 min). Using QCS enabled rapid identification of subjects with data anomalies, allowing the user to correct the data series or reject them partially or entirely, as well as identify fully valid datasets. Overall, the use of the QCS resulted in the automatic classification of 45% of the subjects, including 58 participants who had data with no criterion violation and 21 participants with violations that resulted in the rejection of all dynamic data. The remaining datasets were inspected manually with guidance, allowing acceptance of full datasets from an additional 80 participants and recovery phase data from an additional 16 subjects. Overall, more anomalies occurred with patient data (35% of datasets) compared with healthy controls (15% of datasets). In conclusion, the QCS ensures a standardized data rejection procedure and rigorous objective analysis of dynamic 31 P-MRS data obtained from patients. This methodology contributes to efforts made to standardize 31 P-MRS practices that have been underway for a decade, with the goal of making it an empowered tool for clinical research.

Evaluation of deep learning models for quality control of MR spectra.

Purpose: While 3D MR spectroscopic imaging (MRSI) provides valuable spatial metabolic information, one of the hurdles for clinical translation is its interpretation, with voxel-wise quality control (QC) as an essential and the most time-consuming step. This work evaluates the accuracy of machine learning (ML) models for automated QC filtering of individual spectra from 3D healthy control and patient datasets. Methods: A total of 53 3D MRSI datasets from prior studies (30 neurological diseases, 13 brain tumors, and 10 healthy controls) were included in the study. Three ML models were evaluated: a random forest classifier (RF), a convolutional neural network (CNN), and an inception CNN (ICNN) along with two hybrid models: CNN + RF, ICNN + RF. QC labels used for training were determined manually through consensus of two MRSI experts. Normalized and cropped real-valued spectra was used as input. A cross-validation approach was used to separate datasets into training/validation/testing sets of aggregated voxels. Results: All models achieved a minimum AUC of 0.964 and accuracy of 0.910. In datasets from neurological disease and controls, the CNN model produced the highest AUC (0.982), while the RF model achieved the highest AUC in patients with brain tumors (0.976). Within tumor lesions, which typically exhibit abnormal metabolism, the CNN AUC was 0.973 while that of the RF was 0.969. Data quality inference times were on the order of seconds for an entire 3D dataset, offering drastic time reduction compared to manual labeling. Conclusion: ML methods accurately and rapidly performed automated QC. Results in tumors highlights the applicability to a variety of metabolic conditions.

Polyphenol Supplementation Did Not Affect Insulin Sensitivity and Fat Deposition During One-Month Overfeeding in Randomized Placebo-Controlled Trials in Men and in Women.

Two randomized placebo-controlled double-blind paralleled trials (42 men in Lyon, 19 women in Lausanne) were designed to test 2 g/day of a grape polyphenol extract during 31 days of high calorie-high fructose overfeeding. Hyperinsulinemic-euglycemic clamps and test meals with [1,1,1-13C3]-triolein were performed before and at the end of the intervention. Changes in body composition were assessed by dual-energy X-ray absorptiometry (DEXA). Fat volumes of the abdominal region and liver fat content were determined in men only, using 3D-magnetic resonance imaging (MRI) and magnetic resonance spectroscopy (MRS) at 3T. Adipocyte's size was measured in subcutaneous fat biopsies. Bodyweight and fat mass increased during overfeeding, in men and in women. While whole body insulin sensitivity did not change, homeostasis model assessment of insulin resistance (HOMA-IR) and the hepatic insulin resistance index (HIR) increased during overfeeding. Liver fat increased in men. However, grape polyphenol supplementation did not modify the metabolic and anthropometric parameters or counteract the changes during overfeeding, neither in men nor in women. Polyphenol intake was associated with a reduction in adipocyte size in women femoral fat. Grape polyphenol supplementation did not counteract the moderated metabolic alterations induced by one month of high calorie-high fructose overfeeding in men and women. The clinical trials are registered under the numbers NCT02145780 and NCT02225457 at ClinicalTrials.gov and available at https://clinicaltrials.gov/ct2/show/NCT02145780 and https://clinicaltrials.gov/ct2/show/NCT02225457.

Short echo time dual-frequency MR Elastography with Optimal Control RF pulses.

Magnetic Resonance Elastography (MRE) quantifies the mechanical properties of tissues, typically applying motion encoding gradients (MEG). Multifrequency results allow better characterizations of tissues using data usually acquired through sequential monofrequency experiments. High frequencies are difficult to reach due to slew rate limitations and low frequencies induce long TEs, yielding magnitude images with low SNR. We propose a novel strategy to perform simultaneous multifrequency MRE in the absence of MEGs: using RF pulses designed via the Optimal Control (OC) theory. Such pulses control the spatial distribution of the MRI magnetization phase so that the resulting transverse magnetization reproduces the phase pattern of an MRE acquisition. The pulse is applied with a constant gradient during the multifrequency mechanical excitation to simultaneously achieve slice selection and motion encoding. The phase offset sampling strategy can be adapted according to the excitation frequencies to reduce the acquisition time. Phantom experiments were run to compare the classical monofrequency MRE to the OC based dual-frequency MRE method and showed excellent agreement between the reconstructed shear storage modulus G'. Our method could be applied to simultaneously acquire low and high frequency components, which are difficult to encode with the classical MEG MRE strategy.

Spurious phase correction in rapid metabolic imaging.

IDEAL-type magnetic resonance spectroscopic imaging (MRSI) sequences require the acquisition of several datasets using optimized sampling in the time domain to reconstruct metabolite maps. Each unitary scan consists of a selective slice (2D) or slab (3D) excitation followed by an evolution time and then the acquisition of the spatially encoded signal. It is critical that the phase variation during the evolution time for each scan is only dependent on chemical shifts. In this paper, we described the apparition of spurious phase due to either the transmit or the receive frequency. The presence of this unwanted phase depends on (i) where the commutation between these two frequencies is performed and (ii) how it is done, as there are two phase commutation modes: continuous and coherent. We present the correction needed in function of the different cases. It appears that some solutions are universal. However, it is critical to know which case is implemented on the MRI scanner, which is not always easy information to have. We illustrated several cases with our preclinical MRI by using the IDEAL spiral method on a 13C phantom.

Direct Comparison of Bayesian and Fermi Deconvolution Approaches for Myocardial Blood Flow Quantification: In silico and Clinical Validations.

Cardiac magnetic resonance myocardial perfusion imaging can detect coronary artery disease and is an alternative to single-photon emission computed tomography or positron emission tomography. However, the complex, non-linear MR signal and the lack of robust quantification of myocardial blood flow have hindered its widespread clinical application thus far. Recently, a new Bayesian approach was developed for brain imaging and evaluation of perfusion indexes (Kudo et al., 2014). In addition to providing accurate perfusion measurements, this probabilistic approach appears more robust than previous approaches, particularly due to its insensitivity to bolus arrival delays. We assessed the performance of this approach against a well-known and commonly deployed model-independent method based on the Fermi function for cardiac magnetic resonance myocardial perfusion imaging. The methods were first evaluated for accuracy and precision using a digital phantom to test them against the ground truth; next, they were applied in a group of coronary artery disease patients. The Bayesian method can be considered an appropriate model-independent method with which to estimate myocardial blood flow and delays. The digital phantom comprised a set of synthetic time-concentration curve combinations generated with a 2-compartment exchange model and a realistic combination of perfusion indexes, arterial input dynamics, noise and delays collected from the clinical dataset. The myocardial blood flow values estimated with the two methods showed an excellent correlation coefficient (r 2 > 0.9) under all noise and delay conditions. The Bayesian approach showed excellent robustness to bolus arrival delays, with a similar performance to Fermi modeling when delays were considered. Delays were better estimated with the Bayesian approach than with Fermi modeling. An in vivo analysis of coronary artery disease patients revealed that the Bayesian approach had an excellent ability to distinguish between abnormal and normal myocardium. The Bayesian approach was able to discriminate not only flows but also delays with increased sensitivity by offering a clearly enlarged range of distribution for the physiologic parameters.

Harmonic wideband simultaneous dual-frequency MR Elastography.

Magnetic resonance elastography (MRE) is used to non-invasively quantify viscoelastic properties of tissues based on the measurement of propagation characteristics of shear waves. Because some of these viscoelastic parameters show a frequency dependence, multifrequency analysis allows us to measure the wave propagation dispersion, leading to a better characterization of tissue properties. Conventionally, motion encoding gradients (MEGs) oscillating at the same frequency as the mechanical excitation encode motion. Hence, multifrequency data is usually obtained by sequentially repeating monochromatic wave excitations experiments at different frequencies. The result is that the total acquisition time is multiplied by a factor corresponding to the number of repetitions of monofrequency experiments, which is a major limitation of multifrequency MRE. In order to make it more accessible, a novel single-shot harmonic wideband dual-frequency MRE method is proposed. Two superposed shear waves of different frequencies are simultaneously generated and propagate in a sample. Trapezoidal oscillating MEGs are used to encode mechanical vibrations having frequencies that are an odd multiple of the MEG frequency. The number of phase offsets is optimized to reduce the acquisition time. For this purpose, a sampling method not respecting the Shannon theorem is used to produce a controlled temporal aliasing that allows us to encode both frequencies without any additional examination time. Phantom experiments were run to compare conventional monofrequency MRE with the single-shot dual-frequency MRE method and showed excellent agreement between the reconstructed shear storage moduli G'. In addition, dual-frequency MRE yielded an increased signal-to-noise ratio compared with conventional monofrequency MRE acquisitions when encoding the high frequency component. The novel proposed multifrequency MRE method could be applied to simultaneously acquire more than two frequency components, reducing examination time. Further studies are needed to confirm its applicability in preclinical and clinical models.

Automatic myocardial ischemic lesion detection on magnetic resonance perfusion weighted imaging prior perfusion quantification: A pre-modeling strategy.

Even if cardiovascular magnetic resonance (CMR) perfusion imaging has proven its relevance for visual detection of ischemia, myocardial blood flow (MBF) quantification at the voxel observation scale remains challenging. Integration of an automated segmentation step, prior to perfusion index estimation, might be a significant reconstruction component that could allow sustainable assumptions and constraint enlargement prior to advanced modeling. Current clustering techniques, such as bullseye representation or manual delineation, are not designed to discriminate voxels belonging to the lesion from healthy areas. Hence, the resulting average time-intensity curve, which is assumed to represent the dynamic contrast enhancement inside of a lesion, might be contaminated by voxels with perfectly healthy microcirculation. This study introduces a hierarchical lesion segmentation approach based on time-intensity curve features that considers the spatial particularities of CMR myocardial perfusion. A first k-means clustering approach enables this method to perform coarse clustering, which is refined by a novel spatiotemporal region-growing (STRG) segmentation, thus ensuring spatial and time-intensity curve homogeneity. Over a cohort of 30 patients, myocardial blood flow (MBF) measured in voxels of lesion regions detected with STRG was significantly lower than in regions drawn manually (mean difference = 0.14, 95% CI [0.07, 0.2]) and defined with the bullseye template (mean difference = 0.25, 95% CI [0.17, 0.36]). Over the 90 analyzed slices, the median Dice score calculated against the ground truth ranged between 0.62 and 0.67, the inclusion coefficients ranged between 0.62 and 0.76 and the centroid distances ranged between 0.97 and 3.88 mm. Therefore, though these metrics highlight spatial differences, they could not be used as an index to evaluate the accuracy and performance of the method, which can only be attested by the variability of the MBF clinical index.

Chemical-Shift-Encoded Magnetic Resonance Imaging and Spectroscopy to Reveal Immediate and Long-Term Multi-Organs Composition Changes of a 14-Days Periodic Fasting Intervention: A Technological and Case Report.

Objectives: The aim of this study was to investigate the feasibility of measuring the effects of a 14-day Periodic Fasting (PF) intervention (<200 cal) on multi-organs of primary interest (liver, visceral/subcutaneous/bone marrow fat, muscle) using non-invasive advanced magnetic resonance spectroscopic (MRS) and imaging (MRI) methods. Methods: One subject participated in a 14-day PF under daily supervision of nurses and specialized physicians, ingesting a highly reduced intake: 200 Kcal/day coupled with active walking and drinking at least 3 L of liquids/day. The fasting was preceded by a 7-day pre-fasting vegetarian period and followed by 14 days of stepwise reintroduction of food. The longitudinal study collected imaging and biological data before the fast, at peak fasting, and 7 days, 1 month, and 4 months after re-feeding. Body fat mass in the trunk, abdomen, and thigh, liver and muscle mass, were respectively computed using advanced MRI and MRS signal modeling. Fat fraction, MRI relativity index T2* and susceptibility (Chi), as well as Fatty acid composition, were calculated at all-time points. Results: A decrease in body weight (BW: -9.5%), quadriceps muscle volume (-3.2%), Subcutaneous and Visceral Adipose Tissue (SAT -34.4%; VAT -20.8%), liver fat fraction (PDFF = 1.4 vs. 2.6 % at baseline) but increase in Spine Bone Marrow adipose tissue (BMAT) associated with a 10% increase in global adiposity fraction (PDFF: 54.4 vs. 50.9%) was observed. Femoral BMAT showed minimal changes compared to spinal level, with a slight decrease (-3.1%). Interestingly, fatty acid (FA) pattern changes differed depending on the AT locations. In muscle, all lipids increased after fasting, with a greater increase of intramyocellular lipid (IMCL: from 2.7 to 6.3 mmol/kg) after fasting compared to extramyocellular lipid (EMCL: from 6.2 to 9.5 mmol/kg) as well as Carnosine (6.9 to 8.1 mmol/kg). Heterogenous and reverse changes were also observed after re-feeding depending on the organ. Conclusion: These results suggest that investigating the effects of a 14-day PF intervention using advanced MRI and MRS is feasible. Quantitative MR indexes are a crucial adjunct to further understanding the effective changes in multiple crucial organs especially liver, spin, and muscle, differences between adipose tissue composition and the interplay that occurs during periodic fasting.

3D Chemical Shift-Encoded MRI for Volume and Composition Quantification of Abdominal Adipose Tissue During an Overfeeding Protocol in Healthy Volunteers.

BACKGROUND: Overweight and obesity are major worldwide health concerns characterized by an abnormal accumulation of fat in adipose tissue (AT) and liver. PURPOSE: To evaluate the volume and the fatty acid (FA) composition of the subcutaneous adipose tissue (SAT) and the visceral adipose tissue (VAT) and the fat content in the liver from 3D chemical-shift-encoded (CSE)-MRI acquisition, before and after a 31-day overfeeding protocol. STUDY TYPE: Prospective and longitudinal study. SUBJECTS: Twenty-one nonobese healthy male volunteers. FIELD STRENGTH/SEQUENCE: A 3D spoiled-gradient multiple echo sequence and STEAM sequence were performed at 3T. ASSESSMENT: AT volume was automatically segmented on CSE-MRI between L2 to L4 lumbar vertebrae and compared to the dual-energy X-ray absorptiometry (DEXA) measurement. CSE-MRI and MR spectroscopy (MRS) data were analyzed to assess the proton density fat fraction (PDFF) in the liver and the FA composition in SAT and VAT. Gas chromatography-mass spectrometry (GC-MS) analyses were performed on 13 SAT samples as a FA composition countermeasure. STATISTICAL TESTS: Paired t-test, Pearson's correlation coefficient, and Bland-Altman plots were used to compare measurements. RESULTS: SAT and VAT volumes significantly increased (P < 0.001). CSE-MRI and DEXA measurements were strongly correlated (r = 0.98, P < 0.001). PDFF significantly increased in the liver (+1.35, P = 0.002 for CSE-MRI, + 1.74, P = 0.002 for MRS). FA composition of SAT and VAT appeared to be consistent between localized-MRS and CSE-MRI (on whole segmented volume) measurements. A significant difference between SAT and VAT FA composition was found (P < 0.001 for CSE-MRI, P = 0.001 for MRS). MRS and CSE-MRI measurements of the FA composition were correlated with the GC-MS results (for ndb: rMRS/GC-MS = 0.83 P < 0.001, rCSE-MRI/GC-MS = 0.84, P = 0.001; for nmidb: rMRS/GC-MS = 0.74, P = 0.006, rCSE-MRI/GC-MS = 0.66, P = 0.020) DATA CONCLUSION: The follow-up of liver PDFF, volume, and FA composition of AT during an overfeeding diet was demonstrated through different methods. The CSE-MRI sequence associated with a dedicated postprocessing was found reliable for such quantification. LEVEL OF EVIDENCE: 1 Technical Efficacy: Stage 2 J. Magn. Reson. Imaging 2019;49:1587-1599.

A simplified framework to optimize MRI contrast preparation.

PURPOSE: This article proposes a rigorous optimal control framework for the design of preparation schemes that optimize MRI contrast based on relaxation time differences. METHODS: Compared to previous optimal contrast preparation schemes, a drastic reduction of the optimization parameter number is performed. The preparation scheme is defined as a combination of several block pulses whose flip angles, phase terms and inter-pulse delays are optimized to control the magnetization evolution. RESULTS: The proposed approach reduces the computation time of B 0 -robust preparation schemes to around a minute (whereas several hours were required with previous schemes), with negligible performance loss. The chosen parameterization allows to formulate the total preparation duration as a constraint, which improves the overall compromise between contrast performance and preparation time. Simulation, in vitro and in vivo results validate this improvement, illustrate the straightforward applicability of the proposed approach, and point out its flexibility in terms of achievable contrasts. Major improvement is especially achieved for short-T2 enhancement, as shown by the acquisition of a non-trivial contrast on a rat brain, where a short-T2 white matter structure (corpus callosum) is enhanced compared to surrounding gray matter tissues (hippocampus and neocortex). CONCLUSIONS: This approach proposes key advances for the design of optimal contrast preparation sequences, that emphasize their ability to generate non-standard contrasts, their potential benefit in a clinical context, and their straightforward applicability on any MR system.

Constant gradient elastography with optimal control RF pulses.

This article presents a new motion encoding strategy to perform magnetic resonance elastography (MRE). Instead of using standard motion encoding gradients, a tailored RF pulse is designed to simultaneously perform selective excitation and motion encoding in presence of a constant gradient. The RF pulse is designed with a numerical optimal control algorithm, in order to obtain a magnetization phase distribution that depends on the displacement characteristics inside each voxel. As a consequence, no post-excitation encoding gradients are required. This offers numerous advantages, such as reducing eddy current artifacts, and relaxing the constraint on the gradients maximum switch rate. It also allows to perform MRE with ultra-short TE acquisition schemes, which limits T2 decay and optimizes signal-to-noise ratio. The pulse design strategy is developed and analytically analyzed to clarify the encoding mechanism. Finally, simulations, phantom and ex vivo experiments show that phase-to-noise ratios are improved when compared to standard MRE encoding strategies.

Comparison of MRI-derived vs. traditional estimations of fatty acid composition from MR spectroscopy signals.

INTRODUCTION: The composition of fatty acids in the body is gaining increasing interest, and can be followed up noninvasively by quantitative magnetic resonance spectroscopy (MRS). However, current MRS quantification methods have been shown to provide different quantitative results in terms of lipid signals, with possible varying outcomes for a given biological examination. Quantitative magnetic resonance imaging using multigradient echo sequence (MGE-MRI) has recently been added to MRS approaches. In contrast, these methods fit the undersampled magnetic resonance temporal signal with a simplified model function (expressing the triglyceride [TG] spectrum with only three TG parameters), specific implementations and prior knowledge. In this study, an adaptation of an MGE-MRI method to MRS lipid quantification is proposed. METHODS: Several versions of the method - with time data fully or undersampled, including or excluding the spectral peak T2 knowledge in the fitting - were compared theoretically and on Monte Carlo studies with a time-domain, peak-fitting approach. Robustness, repeatability and accuracy were also inspected on in vitro oil acquisitions and test-retest in vivo subcutaneous adipose tissue acquisitions, adding results from the reference LCModel method. RESULTS: On simulations, the proposed method provided TG parameter estimates with the smallest variability, but with a possible bias, which was mitigated by fitting on undersampled data and considering peak T2 values. For in vitro measurements, estimates for all approaches were correlated with theoretical values and the best concordance was found for the usual MRS method (LCModel and peak fitting). Limited in vivo test-retest variability was found (4.1% for PUFAindx, 0.6% for MUFAindx and 3.6% for SFAindx), as for LCModel (7.6% for PUFAindx, 7.8% for MUFAindx and 3.0% for SFAindx). CONCLUSION: This study shows that fitting the three TG parameters directly on MRS data is one valuable solution to circumvent the poor conditioning of the MRS quantification problem.

In vivo MRS for the assessment of mouse colon using a dedicated endorectal coil: initial findings.

Inflammatory bowel disease is a common group of inflammation conditions that can affect the colon and the rectum. These pathologies require a careful follow-up of patients to prevent the development of colorectal cancer. Currently, conventional endoscopy is used to depict alterations of the intestinal walls, and biopsies are performed on suspicious lesions for further analysis (histology). MRS enables the in vivo analysis of biochemical content of tissues (i.e. without removing any samples). Combined with dedicated endorectal coils (ERCs), MRS provides new ways of characterizing alterations of tissues. An MRS in vivo protocol was specifically set up on healthy mice and on mice chemically treated to induce colitis. Acquisitions were performed on a 4.7 T system using a linear volume birdcage coil for the transmission of the B1 magnetic field, and a dedicated ERC was used for signal reception. Colon-wall complex, lumen and visceral fat were assessed on healthy and treated mice with voxel sizes ranging from 0.125 µL to 2 µL while keeping acquisition times below 3 min. The acquired spectra show various biochemical contents such as α- and ß-methylene but also glycerol backbone and diacyl. Choline was detected in tumoral regions. Visceral fat regions display a high lipid content with no water, whereas colon-wall complex exhibits both high lipid and high water contents. To the best of our knowledge, this is the first time that in vivo MRS using an ERC has been performed in the assessment of colon walls and surrounding structures. It provides keys for the in vivo characterization of small local suspicious lesions and offers complementary solutions to biopsies.

Optimal control design of preparation pulses for contrast optimization in MRI.

This work investigates the use of MRI radio-frequency (RF) pulses designed within the framework of optimal control theory for image contrast optimization. The magnetization evolution is modeled with Bloch equations, which defines a dynamic system that can be controlled via the application of the Pontryagin Maximum Principle (PMP). This framework allows the computation of optimal RF pulses that bring the magnetization to a given state to obtain the desired contrast after acquisition. Creating contrast through the optimal manipulation of Bloch equations is a new way of handling contrast in MRI, which can explore the theoretical limits of the system. Simulation experiments carried out on-resonance quantify the contrast improvement when compared to standard T1 or T2 weighting strategies. The use of optimal pulses is also validated for the first time in both in vitro and in vivo experiments on a small-animal 4.7T MR system. Results demonstrate their robustness to static field inhomogeneities as well as the fact that they can be embedded in standard imaging sequences without affecting standard parameters such as slice selection or echo type. In vivo results on rat and mouse brains illustrate the ability of optimal contrast pulses to create non-trivial contrasts on well-studied structures (white matter versus gray matter).