Viscosities of implantable biomaterials in vocal fold augmentation surgery

Vocal fold vibration depends critically on the viscoelasticity of vocal fold tissues. For instance, phonation threshold pressure, a measure of the "ease" of phonation, has been shown to be directly related to the viscosity of the vibrating mucosa. Various implantable biomaterials have been used in vocal fold augmentation surgery, with implantation sites sometimes close to or inside the mucosa. Yet their viscosities or other mechanical properties are seldom known. This study attempts to provide data on viscosities of commonly used phonosurgical biomaterials. Using a parallel-plate rotational rheometer, oscillatory shear experiments were performed on implantable polytetrafluoroethylene (Teflon or Polytef; Mentor Inc., Hingham, MA), collagen (Zyderm; Collagen Corp., Palo Alto, CA), glutaraldehyde crosslinked (GAX) collagen (Phonagel or Zyplast; Collagen Corp.), absorbable gelatin (Gelfoam; Upjohn Co., Kalamazoo, MI), and human abdominal subcutaneous fat. Samples of human vocal fold mucosal tissues were also tested. Under sinusoidal oscillatory shear at 10 Hz and at 37 degrees C, the dynamic viscosity was 116 Pascal-seconds (Pa-s) for polytetrafluoroethylene, 21 Pa-s for gelatin, 8-13 Pa-s for the two types of collagen, 3 Pa-s for fat, and 1 to 3 Pa-s for vocal fold mucosa. Results extrapolated to 100 Hz also show similar differences among the biomaterials, but all values are an order of magnitude lower because of the typical inverse frequency relation (shear thinning effect) for polymeric and biologic materials. The data suggest that the use of fat for vocal fold augmentation may be more conducive to the "ease" of phonation because of its relatively low viscosity, which is closest to physiologic levels. This implication is probably the most relevant in predicting initial outcome of the postoperative voice before there is any significant assimilation (e.g., resorption and fibrosis) of the implanted biomaterial.

The importance of hyaluronic acid in vocal fold biomechanics

OBJECTIVE

This study examined the influence of hyaluronic acid (HA) on the biomechanical properties of the human vocal fold cover (the superficial layer of the lamina propria).

STUDY DESIGN

Vocal fold tissues were freshly excised from 5 adult male cadavers and were treated with bovine testicular hyaluronidase to selectively remove HA from the lamina propria extracellular matrix (ECM). Linear viscoelastic shear properties (elastic shear modulus and dynamic viscosity) of the tissue samples before and after enzymatic treatment were quantified as a function of frequency (0.01 to 15 Hz) by a parallel-plate rotational rheometer at 37 degrees C.

RESULTS

On removing HA from the vocal fold ECM, the elastic shear modulus (G' ) or stiffness of the vocal fold cover decreased by an average of around 35%, while the dynamic viscosity (eta') increased by 70% at higher frequencies (>1 Hz).

CONCLUSION

The results suggested that HA plays an important role in determining the biomechanical properties of the vocal fold cover. As a highly hydrated glycosaminoglycan in the vocal fold ECM, it likely contributes to the maintenance of an optimal tissue viscosity that may facilitate phonation, and an optimal tissue stiffness that may be important for vocal fundamental frequency control.

SIGNIFICANCE

HA has been proposed as a potential bioimplant for the surgical repair of vocal fold ECM defects (eg, vocal fold scarring and sulcus vocalis). Our results suggested that such clinical use may be potentially optimal for voice production from a biomechanical perspective.

Viscoelastic shear properties of human vocal fold mucosa: measurement methodology and empirical results

A standard method for the empirical rheological characterization of viscoelastic materials was adopted to measure the viscoelastic shear properties of human vocal-fold mucosal tissues (the superficial layer of lamina propria). A parallel-plate rotational rheometer was employed to measure shear deformation of viscoelastic tissue samples, which were deformed between two rigid circular plates rotating in small-amplitude sinusoidal oscillations. Elastic and viscous shear moduli of the samples were then quantified as a function of oscillation frequency (0.01 to 15 Hz) based on shear stresses and strains recorded by the rheometer. Data were obtained from 15 excised human larynges (10 male and 5 female). Results showed that the elastic shear modulus mu and the damping ratio zeta of human vocal-fold mucosa were relatively constant across the range of frequencies observed, while the dynamic viscosity eta decreased monotonically with frequency (i.e., shear thinning). Intersubject differences in mu and eta as large as an order of magnitude were observed, part of which may reflect age-related and gender-related differences. Some molecular interpretations of the findings are discussed.

Microstructural characterization of normal and malignant human prostate tissue with vascular, extracellular, and restricted diffusion for cytometry in tumours magnetic resonance imaging

OBJECTIVE

The aim of this study was to demonstrate the feasibility of the recently introduced Vascular, Extracellular, and Restricted Diffusion for Cytometry in Tumours (VERDICT) framework for imaging prostate cancer with diffusion-weighted magnetic resonance imaging (DW-MRI) within a clinical setting.

MATERIALS AND METHODS

The VERDICT framework is a noninvasive microstructure imaging technique that combines an in-depth diffusion MRI acquisition with a mathematical model to estimate and map microstructural tissue parameters such as cell size and density and vascular perfusion. In total, 8 patients underwent 3-T MRI using 9 different b values (100-3000 s/mm). All patients were imaged before undergoing biopsy. Experiments with VERDICT analyzed DW-MRI data from patients with histologically confirmed prostate cancer in areas of cancerous and benign peripheral zone tissue. For comparison, we also fitted commonly used diffusion models such as the apparent diffusion coefficient (ADC), the intravoxel incoherent motion (IVIM), and the kurtosis model. We also investigated correlations of ADC and kurtosis with VERDICT parameters to gain some biophysical insight into the various parameter values.

RESULTS

Eight patients had prostate cancer in the peripheral zone, with Gleason score 3 + 3 (n = 1), 3 + 4 (n = 6), and 4 + 3 (n = 1). The VERDICT model identified a significant increase in the intracellular and vascular volume fraction estimates in cancerous compared with benign peripheral zone, as well as a significant decrease in the volume of the extracellular-extravascular space (EES) (P = 0.05). This is in agreement with manual segmentation of the biopsies for prostate tissue component analysis, which found proliferation of epithelium, loss of surrounding stroma, and an increase in vasculature. The standard ADC and kurtosis parameters were also significantly different (P = 0.05) between tissue types. There was no significant difference in any of the IVIM parameters (P = 0.11 to 0.29). The VERDICT parametric maps from voxel-by-voxel fitting clearly differentiated cancer from benign regions. Kurtosis and ADC parameters correlated most strongly with VERDICT's intracellular volume fraction but also moderately with the EES and vascular fractions.

CONCLUSIONS

The VERDICT model distinguished tumor from benign areas, while revealing differences in microstructure descriptors such as cellular, vascular, and EES fractions. The parameters of ADC and kurtosis models also discriminated between cancer and benign regions. However, VERDICT provides more specific information that disentangles the various microstructural features underlying the changes in ADC and kurtosis. These results highlight the clinical potential of the VERDICT framework and motivate the construction of a shorter, clinically viable imaging protocol to enable larger trials leading to widespread translation of the method.

Hyaluronic acid (with fibronectin) as a bioimplant for the vocal fold mucosa

OBJECTIVES

To measure the viscoelastic shear properties of hyaluronic acid, with and without fibronectin, and to compare them with those of the human vocal fold mucosa and other phonosurgical biomaterials.

METHODS

Viscoelastic shear properties of various implantable biomaterials (Teflon, gelatin, collagen, fat, hyaluronic acid, and hyaluronic acid with fibronectin) were measured with a parallel-plate rotational rheometer. Elastic and viscous shear properties were quantified as a function of oscillation frequency (0.01-15 Hz) at 37 degrees C.

RESULTS

The shear properties of hyaluronic acid were relatively close to those of human vocal fold mucosal tissues reported previously. Hyaluronic acid at specific concentrations (0.5%-1%), with or without fibronectin, was found to exhibit viscous shear properties (viscous shear modulus and dynamic viscosity) similar to those of the average male and female vocal fold mucosa.

CONCLUSIONS

According to a theory that establishes the effects of tissue shear properties on vocal fold oscillation, phonation threshold pressure (a measure of the ease of phonation) is directly related to the viscous shear modulus of the vibrating vocal fold mucosa. Therefore, our findings suggest that hyaluronic acid, either by itself or mixed with fibronectin, may be a potentially optimal bioimplant for the surgical management of vocal fold mucosal defects and lamina propria deficiencies (e.g., scarring) from a biomechanical standpoint.

Does the voice improve with vocal hygiene education? A study of some instrumental voice measures in a group of kindergarten teachers

This study investigated the efficacy of a program of vocal hygiene education designed for an occupational group at risk for dysphonia, namely, kindergarten teachers. Twelve female teachers explored concepts and knowledge of vocal abuse and vocal hygiene in a workshop session and attempted to practice vocal hygiene for 2 months. They showed significant voice improvement as assessed instrumentally in terms of three acoustic and electroglottographic parameters. There was no significant change of voice in a control group of 13 teachers. The results suggest that teachers were able to improve their voices if they significantly reduced vocal abuses in daily life and practiced specific strategies to maintain classroom order and reduce the use of voice in teaching.

Dual-phase cardiac diffusion tensor imaging with strain correction

Purpose

In this work we present a dual-phase diffusion tensor imaging (DTI) technique that incorporates a correction scheme for the cardiac material strain, based on 3D myocardial tagging.

Methods

In vivo dual-phase cardiac DTI with a stimulated echo approach and 3D tagging was performed in 10 healthy volunteers. The time course of material strain was estimated from the tagging data and used to correct for strain effects in the diffusion weighted acquisition. Mean diffusivity, fractional anisotropy, helix, transverse and sheet angles were calculated and compared between systole and diastole, with and without strain correction. Data acquired at the systolic sweet spot, where the effects of strain are eliminated, served as a reference.

Results

The impact of strain correction on helix angle was small. However, large differences were observed in the transverse and sheet angle values, with and without strain correction. The standard deviation of systolic transverse angles was significantly reduced from 35.9±3.9° to 27.8°±3.5° (p<0.001) upon strain-correction indicating more coherent fiber tracks after correction. Myocyte aggregate structure was aligned more longitudinally in systole compared to diastole as reflected by an increased transmural range of helix angles (71.8°±3.9° systole vs. 55.6°±5.6°, p<0.001 diastole). While diastolic sheet angle histograms had dominant counts at high sheet angle values, systolic histograms showed lower sheet angle values indicating a reorientation of myocyte sheets during contraction.

Conclusion

An approach for dual-phase cardiac DTI with correction for material strain has been successfully implemented. This technique allows assessing dynamic changes in myofiber architecture between systole and diastole, and emphasizes the need for strain correction when sheet architecture in the heart is imaged with a stimulated echo approach.

Further studies of phonation threshold pressure in a physical model of the vocal fold mucosa

This paper reports results of further experimentation on a previously developed physical model of the vocal-fold mucosa [Titze et al., J. Acoust. Soc. Am. 97, 3080-3084 (1995)]. The effects of vocal-fold thickness, epithelial membrane thickness, and prephonatory glottal geometry on phonation threshold pressure were studied. Phonation threshold pressures in the range of 0.13 to 0.34 kPa were observed for an 11-mm-thick vocal fold with a 70-micron-thick "epithelial" membrane for different "mucosal" fluid viscosities. Higher threshold pressure was always obtained for thinner vocal folds and thicker membranes. In another set of experiments, lowest offset threshold pressure was obtained for a rectangular or a near-rectangular prephonatory glottis (with a glottal convergence angle within about +/- 3 degrees). It ranged from 0.07 to 0.23 kPa for different glottal half-widths between 2.0 and 6.0 mm. The threshold for more convergent or divergent glottal geometries was consistently higher. This finding only partially agrees with previous analytical work which predicts a lowest threshold for a divergent glottis. The discrepancy between theory and data is likely to be associated with flow separation from a divergent glottis.

A quantitative output-cost ratio in voice production

A quantitative output-cost ratio (OCR) is proposed for objective use in voice production and is defined as the ratio of the acoustic output intensity to the collision intensity of the vocal folds. Measurement of the OCR is demonstrated in a laboratory experiment using 5 excised larynges and a transducer designed for use on human subjects. Data were gathered at constant fundamental frequency (150 Hz). Subglottal pressure was varied from 1.0 to 1.6 kPa, and glottal width at the vocal processes was varied from a pressed condition to a 2-mm gap. The OCR was plotted as a function of glottal width. With no vocal tract, the excised larynx experiments yielded a broad maxima in the OCR curves, across all subglottal pressure conditions, at about 0.6 mm. Computer simulations indicate that sharper maxima may occur when the influence of the vocal tract is taken into account. The potential clinical utility of the OCR is discussed for treatment of a wide range of voice disorders, including those involving both hyper- and hypoadduction.

Viscoelastic shear properties of human vocal fold mucosa: theoretical characterization based on constitutive modeling

The viscoelastic shear properties of human vocal fold mucosa (cover) were previously measured as a function of frequency [Chan and Titze, J. Acoust. Soc. Am. 106, 2008-2021 (1999)], but data were obtained only in a frequency range of 0.01-15 Hz, an order of magnitude below typical frequencies of vocal fold oscillation (on the order of 100 Hz). This study represents an attempt to extrapolate the data to higher frequencies based on two viscoelastic theories, (1) a quasilinear viscoelastic theory widely used for the constitutive modeling of the viscoelastic properties of biological tissues [Fung, Biomechanics (Springer-Verlag, New York, 1993), pp. 277-292], and (2) a molecular (statistical network) theory commonly used for the rheological modeling of polymeric materials [Zhu et al., J. Biomech. 24, 1007-1018 (1991)]. Analytical expressions of elastic and viscous shear moduli, dynamic viscosity, and damping ratio based on the two theories with specific model parameters were applied to curve-fit the empirical data. Results showed that the theoretical predictions matched the empirical data reasonably well, allowing for parametric descriptions of the data and their extrapolations to frequencies of phonation.

Geometric characterization of the laryngeal cartilage framework for the purpose of biomechanical modeling

Some new anatomic data on the laryngeal cartilage framework have been obtained for the biomechanical modeling of the larynx. This study attempted to define and measure some biomechanically important morphometric features of the laryngeal framework, including both the human and the canine laryngeal frameworks, because the canine larynx has been frequently used as an animal model in gross morphology and in physiological experiments. The larynges of 9 men, 7 women, and 9 dogs were harvested and dissected after death. Linear and angular geometric measurements on the thyroid cartilage, the cricoid cartilage, and the arytenoid cartilage were made with a digital caliper and a protractor, respectively. The results are useful for constructing quantitative biomechanical models of vocal fold vibration and posturing (abduction and adduction), eg, continuum mechanical models and finite-element models of the vocal folds.

Quantitative CEST and MT at 1.5T for monitoring treatment response in glioblastoma: early and late tumor progression during chemoradiation

PURPOSE

Quantitative MRI (qMRI) was performed using a 1.5T protocol that includes a novel chemical exchange saturation transfer/magnetization transfer (CEST/MT) approach. The purpose of this prospective study was to determine if qMRI metrics at baseline, at the 10th and 20th fraction during a 30 fraction/6 week standard chemoradiation (CRT) schedule, and at 1 month following treatment could be an early indicator of response for glioblastoma (GBM).

METHODS

The study included 51 newly diagnosed GBM patients. Four regions-of-interest (ROI) were analyzed: (i) the radiation defined clinical target volume (CTV), (ii) radiation defined gross tumor volume (GTV), (iii) enhancing-tumor regions, and (iv) FLAIR-hyperintense regions. Quantitative CEST, MT, T₁ and T₂ parameters were compared between those patients progressing within 6.9 months (early), and those progressing after CRT (late), using mixed modelling. Exploratory predictive modelling was performed to identify significant predictors of early progression using a multivariable LASSO model.

RESULTS

Results were dependent on the specific tumor ROI analyzed and the imaging time point. The baseline CEST asymmetry within the CTV was significantly higher in the early progression cohort. Other significant predictors included the T₂ of the MT pools (for semi-solid at fraction 20 and water at 1 month after CRT), the exchange rate (at fraction 20) and the MGMT methylation status.

CONCLUSIONS

We observe the potential for multiparametric qMRI, including a novel pulsed CEST/MT approach, to show potential in distinguishing early from late progression GBM cohorts. Ultimately, the goal is to personalize therapeutic decisions and treatment adaptation based on non-invasive imaging-based biomarkers.

Characterization and correction of eddy-current artifacts in unipolar and bipolar diffusion sequences using magnetic field monitoring

Diffusion tensor imaging (DTI) of moving organs is gaining increasing attention but robust performance requires sequence modifications and dedicated correction methods to account for system imperfections. In this study, eddy currents in the "unipolar" Stejskal-Tanner and the velocity-compensated "bipolar" spin-echo diffusion sequences were investigated and corrected for using a magnetic field monitoring approach in combination with higher-order image reconstruction. From the field-camera measurements, increased levels of second-order eddy currents were quantified in the unipolar sequence relative to the bipolar diffusion sequence while zeroth and linear orders were found to be similar between both sequences. Second-order image reconstruction based on field-monitoring data resulted in reduced spatial misalignment artifacts and residual displacements of less than 0.43 mm and 0.29 mm (in the unipolar and bipolar sequences, respectively) after second-order eddy-current correction. Results demonstrate the need for second-order correction in unipolar encoding schemes but also show that bipolar sequences benefit from second-order reconstruction to correct for incomplete intrinsic cancellation of eddy-currents.

Spatially varying properties of the vocal ligament contribute to its eigenfrequency response

The vocal ligament is known to have nonlinear variation in geometry, yet this is rarely considered in empirical or computational studies. This paper investigates the effects of a nonlinear variation of the anterior-to-posterior geometry and the corresponding spatial variation in elastic modulus on the fundamental frequency of vibration for the vocal ligament. Uniaxial tensile tests were performed on a vocal ligament specimen dissected from an excised 60-year-old male larynx. Digital image correlation (DIC) was used to obtain the spatial deformation field for the entire ligament specimen. DIC results revealed that the tensile deformation was very heterogeneous, with the least amount of deformation occurring in the region of smallest cross-sectional area. The elastic modulus was calculated locally and was found to be approximately 10 times higher at the midpoint of the vocal ligament than in the anterior and posterior macula flavae regions. Based on the spatially varying material properties obtained, finite element models (isotropic and transversely isotropic) were created to investigate how the effects of varying cross-section, heterogeneous stiffness, and anisotropy could affect the fundamental frequency of vibration. It was found that the spatial cross-section variation and the spatially varying anisotropy (i.e. modulus ratio) are significant to predictions of the vibration characteristics. Fundamental frequencies predicted with a finite element model are discussed in view of rotatory inertia and contribution of transverse shear deformation.

Estimation of viscoelastic shear properties of vocal-fold tissues based on time-temperature superposition

Empirical data on the viscoelastic shear properties of human vocal-fold mucosa (cover) were recently reported at relatively low frequency (0.01-15 Hz). For the data to become relevant to voice production, attempts have been made to parametrize and extrapolate the data to higher frequencies using constitutive modeling [Chan and Titze, J. Acoust. Soc. Am. 107, 565-580 (2000)]. This study investigated the feasibility of an alternative approach for data extrapolation, namely the principle of time-temperature superposition (TTS). TTS is a hybrid theoretical-empirical approach widely used by rheologists to estimate the viscoelastic properties of polymeric systems at time or frequency scales not readily accessible experimentally. It is based on the observation that for many polymers, the molecular configurational changes that occur in a given time scale at a low temperature correspond to those that occur in a shorter time scale at a higher temperature. Using a rotational rheometer, the elastic shear modulus (G') and viscous shear modulus (G'') of vocal-fold cover (superficial layer of lamina propria) tissue samples were measured at 0.01-15 Hz at relatively low temperatures (5 degrees-37 degrees C). Data were empirically shifted according to TTS, yielding composite "master curves" for predicting the magnitude of the shear moduli at higher frequencies at 37 degrees C. Results showed that TTS may be a feasible approach for estimating the viscoelastic shear properties of vocal-fold tissues at frequencies of phonation (on the order of 100-1000 Hz).

Accuracy and precision of apparent diffusion coefficient measurements on a 1.5 T MR-Linac in central nervous system tumour patients

BACKGROUND AND PURPOSE

MRI linear accelerators (MR-Linacs) may allow treatment adaptation to be guided by quantitative MRI including diffusion-weighted imaging (DWI). The aim of this study was to evaluate the accuracy and precision of apparent diffusion coefficient (ADC) measurements from DWI on a 1.5 T MR-Linac in patients with central nervous system (CNS) tumours through comparison with a diagnostic scanner.

MATERIALS AND METHODS

CNS patients were treated using a 1.5 T Elekta Unity MR-Linac. DWI was acquired during MR-Linac treatment and on a Philips Ingenia 1.5 T. The agreement between the two scanners on median ADC over the gross tumour/clinical target volumes (GTV/CTV) and in brain regions (white/grey matter, cerebrospinal fluid (CSF)) was computed. Repeated scans were used to estimate ADC repeatability. Daily changes in ADC over the GTV of high-grade gliomas were characterized from MR-Linac scans.

RESULTS

DWI from 59 patients was analyzed. MR-Linac ADC measurements showed a small bias relative to Ingenia measurements in white matter, grey matter, GTV, and CTV (bias: -0.05 ± 0.03, -0.08 ± 0.05, -0.1 ± 0.1, -0.08 ± 0.07 μm²/ms). ADC differed substantially in CSF (bias: -0.5 ± 0.3 μm²/ms). The repeatability of MR-Linac ADC over white/grey matter was similar to previous reports (coefficients of variation for median ADC: 1.4%/1.8%). MR-Linac ADC changes in the GTV were detectable.

CONCLUSIONS

It is possible to obtain ADC measurements in the brain on a 1.5 T MR-Linac that are comparable to those of diagnostic-quality scanners. This technical validation study adds to the foundation for future studies that will correlate brain tumour ADC with clinical outcomes.

High-resolution diffusion tensor imaging of the human kidneys using a free-breathing, multi-slice, targeted field of view approach

Fractional anisotropy (FA) obtained by diffusion tensor imaging (DTI) can be used to image the kidneys without any contrast media. FA of the medulla has been shown to correlate with kidney function. It is expected that higher spatial resolution would improve the depiction of small structures within the kidney. However, the achievement of high spatial resolution in renal DTI remains challenging as a result of respiratory motion and susceptibility to diffusion imaging artefacts. In this study, a targeted field of view (TFOV) method was used to obtain high-resolution FA maps and colour-coded diffusion tensor orientations, together with measures of the medullary and cortical FA, in 12 healthy subjects. Subjects were scanned with two implementations (dual and single kidney) of a TFOV DTI method. DTI scans were performed during free breathing with a navigator-triggered sequence. Results showed high consistency in the greyscale FA, colour-coded FA and diffusion tensors across subjects and between dual- and single-kidney scans, which have in-plane voxel sizes of 2 × 2 mm(2) and 1.2 × 1.2 mm(2) , respectively. The ability to acquire multiple contiguous slices allowed the medulla and cortical FA to be quantified over the entire kidney volume. The mean medulla and cortical FA values were 0.38 ± 0.017 and 0.21 ± 0.019, respectively, for the dual-kidney scan, and 0.35 ± 0.032 and 0.20 ± 0.014, respectively, for the single-kidney scan. The mean FA between the medulla and cortex was significantly different (p < 0.001) for both dual- and single-kidney implementations. High-spatial-resolution DTI shows promise for improving the characterization and non-invasive assessment of kidney function.

Osseointegrated implants for single tooth replacement in general practice: a 1-year report from a multicentre prospective study

The single tooth implant is a treatment option for the replacement of missing single teeth and in many cases is the treatment of choice. It is, however, an expensive treatment requiring a co-ordinated approach to the surgical and restorative aspects of treatment. In this study, a group of dentists in general and restorative dental practice and with no previous experience in implant surgery underwent an intensive training course in all aspects of implant treatment for single teeth. Using a system of simplified instrumentation with a strict adherence to protocol, the group installed and restored single tooth implants ad modum Brånemark in a wide range of clinical situations. At the one year follow-up period following crown insertion, the success rate of treatment compared favourably with results reported from centres using the specialist team approach to treatment. The results of this study indicate that further consideration should be given to the training of general dentists so that improved delivery of dental health care can be provided at a more economic level.

High grade glioma radiation therapy on a high field 1.5 Tesla MR-Linac - workflow and initial experience with daily adapt-to-position (ATP) MR guidance: A first report

Purpose

This study reports the workflow and initial clinical experience of high grade glioma (HGG) radiotherapy on the 1.5 T MR-Linac (MRL), with a focus on the temporal variations of the tumor and feasibility of multi-parametric image (mpMRI) acquisition during routine treatment workflow.

Materials and methods

Ten HGG patients treated with radiation within the first year of the MRL's clinical operation, between October 2019 and August 2020, were identified from a prospective database. Workflow timings were recorded and online adaptive plans were generated using the Adapt-To-Position (ATP) workflow. Temporal variation within the FLAIR hyperintense region (FHR) was assessed by the relative FHR volumes (n = 281 contours) and migration distances (maximum linear displacement of the volume). Research mpMRIs were acquired on the MRL during radiation and changes in selected functional parameters were investigated within the FHR.

Results

All patients completed radiotherapy to a median dose of 60 Gy (range, 54-60 Gy) in 30 fractions (range, 30-33), receiving a total of 287 fractions on the MRL. The mean in-room time per fraction with or without post-beam research imaging was 42.9 minutes (range, 25.0-69.0 minutes) and 37.3 minutes (range, 24.0-51.0 minutes), respectively. Three patients (30%) required re-planning between fractions 9 to 12 due to progression of tumor and/or edema identified on daily MRL imaging. At the 10, 20, and 30-day post-first fraction time points 3, 3, and 4 patients, respectively, had a FHR volume that changed by at least 20% relative to the first fraction. Research mpMRIs were successfully acquired on the MRL. The median apparent diffusion coefficient (ADC) within the FHR and the volumes of FLAIR were significantly correlated when data from all patients and time points were pooled (R=0.68, p<.001).

Conclusion

We report the first clinical series of HGG patients treated with radiotherapy on the MRL. The ATP workflow and treatment times were clinically acceptable, and daily online MRL imaging triggered adaptive re-planning for selected patients. Acquisition of mpMRIs was feasible on the MRL during routine treatment workflow. Prospective clinical outcomes data is anticipated from the ongoing UNITED phase 2 trial to further refine the role of MR-guided adaptive radiotherapy.

Chemical exchange saturation transfer MRI in central nervous system tumours on a 1.5 T MR-Linac

PURPOSE

To describe the implementation and initial results of using Chemical Exchange Saturation Transfer (CEST) for monitoring patients with central nervous system (CNS) tumours treated using a 1.5 tesla MR-guided radiotherapy system.

METHODS

CNS patients were treated with up to 30 fractions (total dose up to 60 Gy) using a 1.5 T Elekta Unity MR-Linac. CEST scans were obtained in 54 subjects at one or more time points during treatment. CEST metrics, including the amide magnetization transfer ratio (MTR_Amide), nuclear Overhauser effect (NOE) MTR (MTR_NOE) and asymmetry, were quantified in phantoms and CNS patients. The signal was investigated between tumour and white matter, across time, and across disease categories including high- and low-grade tumours.

RESULTS

The gross tumour volume (GTV) exhibited lower MTR_Amide and MTR_NOE and higher asymmetry compared to contralateral normal appearing white matter. Signal changes in the GTV during fractionated radiotherapy were observed. There were differences between high- and low-grade tumours, with higher CEST asymmetry associated with higher grade disease.

CONCLUSION

CEST MRI using a 1.5 T MR-Linac was demonstrated to be feasible for in vivo imaging of CNS tumours. CEST images showed tumour/white-matter contrast, temporal CEST signal changes, and associations with tumour grade. These results show promise for the eventual goal of using metabolic imaging to inform the design of adaptive radiotherapy protocols.

Using the robust principal component analysis algorithm to remove RF spike artifacts from MR images

Brief bursts of RF noise during MR data acquisition (“k‐space spikes”) cause disruptive image artifacts, manifesting as stripes overlaid on the image. RF noise is often related to hardware problems, including vibrations during gradient‐heavy sequences, such as diffusion‐weighted imaging. In this study, we present an application of the Robust Principal Component Analysis (RPCA) algorithm to remove spike noise from k‐space. Methods: Corrupted k‐space matrices were decomposed into their low‐rank and sparse components using the RPCA algorithm, such that spikes were contained within the sparse component and artifact‐free k‐space data remained in the low‐rank component. Automated center refilling was applied to keep the peaked central cluster of k‐space from misclassification in the sparse component. Results: This algorithm was demonstrated to effectively remove k‐space spikes from four data types under conditions generating spikes: (i) mouse heart T₁ mapping, (ii) mouse heart cine imaging, (iii) human kidney diffusion tensor imaging (DTI) data, and (iv) human brain DTI data. Myocardial T₁ values changed by 86.1 ± 171 ms following despiking, and fractional anisotropy values were recovered following despiking of DTI data. Conclusion: The RPCA despiking algorithm will be a valuable postprocessing method for retrospectively removing stripe artifacts without affecting the underlying signal of interest. Magn Reson Med 75:2517–2525, 2016. © 2015 The Authors. Magnetic Resonance in Medicine published by Wiley Periodicals, Inc. on behalf of International Society for Magnetic Resonance in Medicine. This is an open access article under the terms of the Creative Commons Attribution License, which permits use, distribution and reproduction in any medium, provided the original work is properly cited.

Removing undersampling artifacts in DCE-MRI studies using independent components analysis

In breast MRI mammography both high temporal resolution and high spatial resolution have been shown to be important in improving specificity. Adaptive methods such as projection reconstruction time-resolved imaging of contrast kinetics (PR-TRICKS) allow images to be reconstructed at various temporal and spatial resolutions from the same data set. The main disadvantage is that the undersampling, which is necessary to produce high temporal resolution images, leads to the presence of streak artifacts in the images. We present a novel method of removing these artifacts using independent components analysis (ICA) and demonstrate that this results in a significant improvement in image quality for both simulation studies and for patient dynamic contrast-enhanced (DCE)-MRI images. We also investigate the effect of artifacts on two quantitative measures of contrast enhancement. Using simulation studies we demonstrate that streak artifacts lead to pronounced periodic oscillations in pixel concentration curves which, in turn, lead to increased errors and introduce bias into heuristic measurements. ICA filtering significantly reduces this bias and improves accuracy. Pharmacokinetic modeling was more robust and there was no evidence of bias due to the presence of streak artifacts. ICA filtering did not significantly reduce the errors in the estimated pharmacokinetic parameters; however, the chi-squared error was greatly reduced after ICA filtering.