Non-Contrast-Enhanced Functional Lung MRI to Evaluate Treatment Response of Allergic Bronchopulmonary Aspergillosis in Patients With Cystic Fibrosis: A Pilot Study

BACKGROUND

Allergic bronchopulmonary aspergillosis (ABPA) in cystic fibrosis (CF) patients is associated with severe lung damage and requires specific therapeutic management. Repeated imaging is recommended to both diagnose and follow-up response to treatment of ABPA in CF. However, high risk of cumulative radiation exposure requires evaluation of free-radiation techniques in the follow-up of CF patients with ABPA.

PURPOSE

To evaluate whether Fourier decomposition (FD) functional lung MRI can detect response to treatment of ABPA in CF patients.

STUDY TYPE

Retrospective longitudinal.

POPULATION

Twelve patients (7M, median-age:14 years) with CF and ABPA with pre- and post-treatment MRI.

FIELD STRENGTH/SEQUENCE

2D-balanced-steady-state free-precession (bSSFP) sequence with FD at 1.5T.

ASSESSMENT

Ventilation-weighted (V) and perfusion-weighted (Q) maps were obtained after FD processing of 2D-coronal bSSFP time-resolved images acquired before and 3-9 months after treatment. Defects extent was assessed on the functional maps using a qualitative semi-quantitative score (0 = absence/negligible, 1 = <50%, 2 = >50%). Mean and coefficient of variation (CV) of the ventilation signal-intensity (VSI) and the perfusion signal-intensity (QSI) were calculated. Measurements were performed independently by three readers and averaged. Inter-reader reproducibility of the measurements was assessed. Pulmonary function tests (PFTs) were performed within 1 week of both MRI studies as markers of the airflow-limitation severity.

STATISTICAL TESTS

Comparisons of medians were performed using the paired Wilcoxon-test. Reproducibility was assessed using intraclass correlation coefficient (ICC). Correlations between MRI and PFT parameters were assessed using the Spearman-test (rho correlation-coefficient). A P-value <0.05 was considered as significant.

RESULTS

Defects extent on both V and Q maps showed a significant reduction after ABPA treatment (4.25 vs. 1.92 for V-defect-score and 5 vs. 2.75 for Q-defect-score). VSI_mean was significantly increased after treatment (280 vs. 167). Qualitative analyses reproducibility showed an ICC > 0.90, while the ICCs of the quantitative measurements was almost perfect (>0.99). Changes in VSI_cv and QSI_cv before and after treatment correlated inversely with changes of FEV1%p (rho = -0.68 for both).

DATA CONCLUSION

Non-contrast-enhanced FD lung MRI has potential to reproducibly assess response to treatment of ABPA in CF patients and correlates with PFT obstructive parameters.

EVIDENCE LEVEL

4 TECHNICAL EFFICACY: Stage 3.

Sea-ice melt determines seasonal phytoplankton dynamics and delimits the habitat of temperate Atlantic taxa as the Arctic Ocean atlantifies

The Arctic Ocean is one of the regions where anthropogenic environmental change is progressing most rapidly and drastically. The impact of rising temperatures and decreasing sea ice on Arctic marine microbial communities is yet not well understood. Microbes form the basis of food webs in the Arctic Ocean, providing energy for larger organisms. Previous studies have shown that Atlantic taxa associated with low light are robust to more polar conditions. We compared to which extent sea ice melt influences light-associated phytoplankton dynamics and biodiversity over two years at two mooring locations in the Fram Strait. One mooring is deployed in pure Atlantic water, and the second in the intermittently ice-covered Marginal Ice Zone. Time-series analysis of amplicon sequence variants abundance over a 2-year period, allowed us to identify communities of co-occurring taxa that exhibit similar patterns throughout the annual cycle. We then examined how alterations in environmental conditions affect the prevalence of species. During high abundance periods of diatoms, polar phytoplankton populations dominated, while temperate taxa were weakly represented. Furthermore, we found that polar pelagic and ice-associated taxa, such as Fragilariopsis cylindrus and Melosira arctica, were more common in Atlantic conditions, while temperate taxa, such as Odontella aurita and Proboscia alata, were less abundant under polar conditions. This suggests that sea ice melt may act as a barrier to the northward expansion of temperate phytoplankton, preventing their dominance in regions still strongly influenced by polar conditions. Our findings highlight the complex interactions between sea ice melt, phytoplankton dynamics, and biodiversity in the Arctic.

Chest magnetic resonance imaging in cystic fibrosis: technique and clinical benefits

Cystic fibrosis (CF) is one of the most common inherited and life-shortening pulmonary diseases in the Caucasian population. With the widespread introduction of newborn screening and the development of modulator therapy, tremendous advances have been made in recent years both in diagnosis and therapy. Since paediatric CF patients tend to be younger and have lower morbidity, the type of imaging modality that should be used to monitor the disease is often debated. Computed tomography (CT) is sensitive to many pulmonary pathologies, but radiation exposure limits its use, especially in children and adolescents. Conventional pulmonary magnetic resonance imaging (MRI) is a valid alternative to CT and, in most cases, provides sufficient information to guide treatment. Given the expected widespread availability of sequences with ultra-short echo times, there will be even fewer reasons to perform CT for follow-up of patients with CF. This review aims to provide an overview of the process and results of monitoring CF with MRI, particularly for centres not specialising in the disease.

Dynamic mode decomposition of dynamic MRI for assessment of pulmonary ventilation and perfusion

PURPOSE

To introduce dynamic mode decomposition (DMD) as a robust alternative for the assessment of pulmonary functional information from dynamic non-contrast-enhanced acquisitions.

METHODS

Pulmonary fractional ventilation and normalized perfusion maps were obtained using DMD from simulated phantoms as well as in vivo dynamic acquisitions of healthy volunteers at 1.5T. The performance of DMD was compared with conventional Fourier decomposition (FD) and matrix pencil (MP) methods in estimating functional map values. The proposed method was evaluated based on estimated signal amplitude in functional maps across varying number of measurements.

RESULTS

Quantitative assessments performed on phantoms and in vivo measurements indicate that DMD is capable of successfully obtaining pulmonary functional maps. Specifically, compared to FD and MP methods, DMD is able to reduce variations in estimated amplitudes across different number of measurements. This improvement is evident in the fractional ventilation and normalized perfusion maps obtain from phantom simulations with frequency variations and noise, as well as in the maps obtained from in vivo measurements.

CONCLUSIONS

A robust method for accurately estimating pulmonary ventilation and perfusion related signal changes in dynamic acquisitions is presented. The proposed method uses DMD to obtain functional maps reliably, while reducing amplitude variations caused by differences in number of measurements.

Phase-cycled balanced SSFP imaging for non-contrast-enhanced functional lung imaging

PURPOSE

To introduce phase-cycled balanced SSFP (bSSFP) acquisition as an alternative in Fourier decomposition MRI for improved robustness against field inhomogeneities.

METHODS

Series 2D dynamic lung images were acquired in 5 healthy volunteers at 1.5 T and 3 T using bSSFP sequence with multiple RF phase increments and compared with conventional single RF phase increment acquisitions. The approach was evaluated based on functional map homogeneity analysis, while ensuring image and functional map quality by means of SNR and contrast-to-noise ratio analyses.

RESULTS

At both field strengths, functional maps obtained with phase-cycled acquisitions displayed improved robustness against local signal losses compared with single-phase acquisitions. The coefficient of variation (mean ± SD, across volunteers) measured in the ventilation maps resulted in 29.7 ± 2.6 at 1.5 T and 37.5 ± 3.1 at 3 T for phase-cycled acquisitions, compared with 39.9 ± 5.2 at 1.5 T and 49.5 ± 3.7 at 3 T for single-phase acquisitions, indicating a significant improvement ( p < 0.05 $$ p<0.05 $$ ) in ventilation map homogeneity.

CONCLUSIONS

Phase-cycled bSSFP acquisitions improve robustness against field inhomogeneity artifacts and significantly improve ventilation map homogeneity at both field strengths. As such, phase-cycled bSSFP may serve as a robust alternative in lung function assessments.

A dual center and dual vendor comparison study of automated perfusion-weighted phase-resolved functional lung magnetic resonance imaging with dynamic contrast-enhanced magnetic resonance imaging in patients with cystic fibrosis

For sensitive diagnosis and monitoring of pulmonary disease, ionizing radiation-free imaging methods are of great importance. A noncontrast and free-breathing proton magnetic resonance imaging (MRI) technique for assessment of pulmonary perfusion is phase-resolved functional lung (PREFUL) MRI. Since there is no validation of PREFUL MRI across different centers and scanners, the purpose of this study was to compare perfusion-weighted PREFUL MRI with the well-established dynamic contrast-enhanced (DCE) MRI across two centers on scanners from two different vendors. Sixteen patients with cystic fibrosis (CF) (Center 1: 10 patients; Center 2: 6 patients) underwent PREFUL and DCE MRI at 1.5T in the same imaging session. Normalized perfusion-weighted values and perfusion defect percentage (QDP) values were calculated for the whole lung and three central slices (dorsal, central, ventral of the carina). Obtained parameters were compared using Pearson correlation, Spearman correlation, Bland-Altman analysis, Wilcoxon signed-rank test, and Wilcoxon rank-sum test. Moderate-to-strong correlations between normalized perfusion-weighted PREFUL and DCE values were found (posterior slice: r = 0.69, p < 0.01). Spatial overlap of PREFUL and DCE QDP maps showed an agreement of 79.4% for the whole lung. Further, spatial overlap values of Center 1 were not significantly different to those of Center 2 for the three central slices (p > 0.07). The feasibility of PREFUL MRI across two different centers and two different vendors was shown in patients with CF and obtained results were in agreement with DCE MRI.

Validation of Phase-Resolved Functional Lung (PREFUL) Magnetic Resonance Imaging Pulse Wave Transit Time Compared to Echocardiography in Chronic Obstructive Pulmonary Disease

BACKGROUND

Phase-resolved functional lung (PREFUL) magnetic resonance imaging (MRI) pulmonary pulse wave transit time (pPTT) is a contrast agent free, vascular imaging biomarker, but has not been validated in chronic obstructive pulmonary disease (COPD).

PURPOSE

To validate PREFUL with echocardiographic pPTT as a reference standard and to compare arterial/venous pPTT mapping with spirometry and clinical parameters.

STUDY TYPE

Prospective.

POPULATION

Twenty-one patients (62% female) with COPD and 44 healthy participants (50% female).

FIELD STRENGTH/SEQUENCE

1.5 T; 2D-spoiled gradient-echo sequence.

ASSESSMENT

Three coronal PREFUL MRI slices, echocardiography, and spirometry including forced expiratory volume in 1 second (FEV1, liter) and predicted defined as FEV1 in% divided by the population average FEV1%, were performed. Pulmonary pulse transit time from the main artery to the microvasculature (PREFUL pPTT), to the right upper lobe vein (PREFUL pPTT_av , echo pPTT_av ), from microvasculature to right upper lobe vein (PREFUL_vein ) and the ratio of PREFUL pPTT to PREFUL pPTT_vein were calculated. Body mass index (BMI), Global Initiative for COPD (GOLD) stage 1-4, disease duration, and cigarette packs smoked per day multiplied by the smoked years (pack years) were computed.

STATISTICAL TESTS

Shapiro-Wilk-test, paired-two-sided-t-tests, Bland-Altman-analysis, coefficient of variation, Pearson ρ were applied, pPTT data were compared between 21 subjects from the 44 healthy subjects who were age- and sex-matched to the COPD cohort, P < 0.05 was considered statistically significant.

RESULTS

PREFUL pPTT_av significantly correlated with echo pPTT_av (ρ = 0.95) with 1.85 msec bias, 95% limits of agreement: 55.94 msec, -52.23 msec in all participants (P = 0.59). In the healthy participants, PREFUL and echo pPTT_av significantly correlated with age (ρ = 0.81, ρ = 0.78), FEV1 (ρ = -0.47, ρ = -0.34) and BMI (ρ = 0.56, ρ = 0.51). In COPD patients, PREFUL pPTT significantly correlated with FEV1 predicted (ρ = -0.59), GOLD (ρ = 0.53), disease duration (ρ = 0.54), and pack years (ρ = 0.49).

DATA CONCLUSION

Arteriovenous PTT measured by PREFUL MRI corresponds precisely to echocardiography and appears to be feasible even in severe COPD.

EVIDENCE LEVEL

1 TECHNICAL EFFICACY: Stage 2.

Perfusion quantification using voxel-wise proton density and median signal decay in PREFUL MRI

PURPOSE

Contrast-free lung MRI based on Fourier decomposition is an attractive method to monitor various lung diseases. However, the accuracy of the current perfusion quantification is limited. In this study, a new approach for perfusion quantification based on voxel-wise proton density and median signal decay toward the steady state for Fourier decomposition-based techniques is proposed called Q_Quantified (Q_Quant ).

METHODS

Twenty patients with chronic obstructive pulmonary disease and 18 patients with chronic thromboembolic pulmonary hypertension received phase-resolved functional lung-MRI (PREFUL) and dynamic contrast-enhanced (DCE)-MRI. Nine healthy participants received phase-resolved functional lung-MRI only. Median values of Q_Quant were compared to a Fourier decomposition perfusion quantification presented by Kjørstad et al (Q_Kjørstad ) and validated toward pulmonary blood flow derived by DCE-MRI (PBF_DCE ). Blood fraction maps determined by the new approach were calculated. Regional and global correlation coefficients were calculated, and Bland-Altman plots were created. Histogram analyses of all cohorts were created.

RESULTS

The introduced parameter Q_Quant showed only 2 mL/min/100 mL mean deviation to PBF_DCE in the patient cohort and showed less bias than Q_Kjørstad . Significant increases of regional correlation with PBF_DCE were achieved (r = 0.3 vs. r = 0.2, P < .01*). The trend of global correlation toward PBF_DCE is not uniform, showing higher values for Q_Kjørstad in the chronic obstructive pulmonary disease cohort than for Q_Quant and vice versa in the chronic thromboembolic pulmonary hypertension cohort. In contrast to Q_Kjørstad , Q_Quant perfusion maps indicate a physiologic dorsoventral gradient in supine position similar to PBF_DCE with similar value distribution in the histograms.

CONCLUSION

We proposed a new approach for perfusion quantification of phase-resolved functional lung measurements. The developed parameter Q_Quant reveals a higher accuracy compared to Q_Kjørstad .

Non-contrast pulmonary perfusion MRI in patients with cystic fibrosis

PURPOSE

This study aimed to assess the feasibility of Self-gated Non-Contrast-Enhanced Functional Lung (SENCEFUL) MRI for detection of pulmonary perfusion deficits in patients with cystic fibrosis.

METHODS

Twenty patients with cystic fibrosis and 20 matched healthy controls underwent SENCEFUL-MRI at 1.5 T with reconstruction of perfusion and perfusion phase maps (i.e. comparable to pulse wave delays). Four blinded readers rated both types of maps separately followed by simultaneous assessment thereof. Perfusion phase data was plotted in histograms and a Peak-to-Offset ratio was calculated for comparison to subjective scoring and correlation (Spearman) to lung function parameters. Sensitivity, specificity and positive and negative predictive values were calculated for subjective scoring and Peak-to-Offset ratios. Intraclass correlation (ICC) was used to assess the interrater agreement.

RESULTS

Readers attributed pathological ratings 2.2-3.5 times more frequently to the CF-group. The sensitivity with regard to a correct assignment to CF was similar between ratings (perfusion only vs. perfusions phase only vs. simultaneous assessment: 0.54-0.56), while specificity increased from 0.75 to 0.85 for simultaneous assessment. ICC was 0.77-0.84 for subjective scoring. ROC-analysis of Peak-to-Offset ratios on a mean per-subject basis revealed a sensitivity of 0.75 and specificity of 0.85 (PPV 0.83, NPV 0.77). Functional pulmonary parameters indicative of bronchial obstruction and Peak-to-Offset ratios showed positive correlation (FEV1: 0.77; FEF75: 0.76).

CONCLUSIONS

SENCEFUL-MRI bears the potential for monitoring CF including disease-associated patterns of altered pulmonary perfusion. The proposed Peak-to-Offset ratio derived from pulmonary perfusion phase measurements could represent an objective future marker for perfusion impairment.

Flow Volume Loop and Regional Ventilation Assessment Using Phase-Resolved Functional Lung (PREFUL) MRI: Comparison With 129 Xenon Ventilation MRI and Lung Function Testing

BACKGROUND

Regional flow volume loop ventilation-weighted noncontrast-enhanced proton lung MRI in free breathing has emerged as a novel technique for assessment of regional lung ventilation, but has yet not been validated with ¹²⁹ Xenon MRI (¹²⁹ Xe-MRI), a direct visualization of ventilation in healthy volunteers, cystic fibrosis (CF), and chronic obstructive pulmonary disease (COPD) patients.

PURPOSE

To compare regional ventilation and regional flow volume loops measured by noncontrast-enhanced ventilation-weighted phase-resolved functional lung MRI (PREFUL-MRI) with ¹²⁹ Xe-MRI ventilation imaging and with lung function test parameters.

STUDY TYPE

Retrospective study.

POPULATION

Twenty patients with COPD, eight patients with CF, and six healthy volunteers.

FIELD STRENGTH/SEQUENCE

PREFUL and ¹²⁹ Xe-MRI gradient echo sequences were acquired at 1.5T.

ASSESSMENT

Coronal slices of PREFUL-MRI (free breathing) and ¹²⁹ Xe-MRI (single breath-hold) were acquired on the same day, matched by their ventrodorsal position and coregistered for evaluation. Ventilation defect percentage (VDP) was calculated based on regional ventilation (RV), regional flow volume loops (RFVL), or ¹²⁹ Xe-MRI with two different threshold methods. A combined VDP was calculated for RV and RFVL. Additionally, lung function testing was performed (such as the forced expiratory volume in 1 second [FEV₁ ]) was used.

STATISTICAL TESTS

The obtained parameters were compared using Wilcoxon tests, correlated using Spearman's correlation coefficient (r), and agreement between PREFUL and ¹²⁹ Xe-MRI parameters was assessed using Bland-Altman analysis and Dice coefficients.

RESULTS

VDP measured by PREFUL and ¹²⁹ Xe were significantly correlated with both thresholding techniques (r = 0.62-0.69, P < 0.05 for all) and with lung function test parameters. Combined RV and RFVL PREFUL defect maps correlated with lung function testing (eg, with FEV₁ r = -0.87 P < 0.05), and showed better regional agreement to ¹²⁹ Xe-MRI ventilation defects (Dice coefficient defect 0.413) with significantly higher VDP values (10.2 ± 27.3, P = 0.04) than either PREFUL defect map alone.

DATA CONCLUSION

Combined RV and RFVL PREFUL defect maps likely increase sensitivity to mild airway obstruction with increased VDP values compared to ¹²⁹ Xe-MRI, and correlate strongly with lung function test parameters.

LEVEL OF EVIDENCE

3 TECHNICAL EFFICACY STAGE: 2.

Meteorological rhythms of respiratory and circulatory diseases revealed by Harmonic Analysis

The intricately fluctuating onset of respiratory and circulatory diseases displays rhythms of multi-scaled meteorological conditions due to their sensitivity to weather changes. The intrinsic meteorological rhythms of these diseases are revealed in this bio-meteorological study via Fourier decomposition and harmonic analysis. Daily emergency room (ER) visit data for respiratory and circulatory diseases from three comprehensive hospitals in Haidian district of Beijing, China were used in the analysis. Meteorological data included three temperature metrics, relative humidity, sunshine duration, daily mean air pressure, and wind speed. The Fourier decomposition and harmonic analysis on ER visits and meteorological variables involve frequency, period, and power of all harmonics. The results indicated that: i) for respiratory morbidity, a strong climatic annual rhythm responding to annual temperature change was firstly revealed; its ratio of spectral density was 16-33%. Moreover, significant correlations existed between the high-frequency fluctuations (<30 d) of morbidity and short-term harmonics of humidity and solar duration. High-frequency harmonics of temperature and pressure showed no statistically significant effect. ii) With regard to all types of circulatory morbidity, their annual periodicity was weaker than that of respiratory diseases, whose harmonic energy took a ratio less than 8%. Besides, the power of all high-frequency harmonics of circulatory morbidity accounted for up to 70-90% in the original sequences, and their relationship to many short-term meteorological factors were significant, including the mean and maximum temperatures, wind speed, and solar duration. iii) The weekly rhythm appeared in respiratory ER visits with 15% of harmonic variance but not prominent in circulatory morbidity. In summary, by decomposing the sequence of respiratory and circulatory diseases as well as recognizing their meteorological rhythms, different responses to meteorological conditions on various time scales were identified.

Effect of intravenously injected gadolinium-based contrast agents on functional lung parameters derived by PREFUL MRI

PURPOSE

To evaluate the influence of intravenously administered gadolinium-based contrast agents on functional ventilation and perfusion parameters derived by phase-resolved functional lung (PREFUL) MRI.

METHODS

Fourteen participants underwent functional MRI at 1.5T using a 2D spoiled gradient echo sequence during free breathing. Three data sets of PREFUL images were obtained-the 1st data set was acquired in mean 33:46 min (SD = 6:20 min) prior, the 2nd and 3rd data sets 43 and 91 s (both SD = 1.9 s), respectively, after i.v. application of gadobutrol. Full respiratory and cardiac cycles were reconstructed and functional parameters of regional ventilation (RV), perfusion (Q), and quantified perfusion (Q_Quant ) together with perfusion-defected percentages (QDP), ventilation-defected percentages (VDP), and ventilation-perfusion match (VQM) were calculated and compared for systematic differences between the acquired data sets.

RESULTS

RV- and Q-values presented no significant alteration after gadobutrol administration. Consequently, QDP, VDP, and VQ maps were not significantly different. Sørensen-Dice coefficients of QDP and VDP maps between the different series varied up to ±9%. Q_Quant was significantly increased after the application of gadobutrol (1st vs. 2nd series, P = 0.0021; 1st vs. 3rd, P = 0.0188), which can be explained by the velocity-dependent signal in the completely blood-filled voxel (ROI of the aorta) after shortening of T₁ relaxation time (1st vs. 2nd series, P = 0.0003; 1st vs. 3rd series, P = 0.0008).

CONCLUSION

Except for quantified perfusion, all evaluated functional parameters including ventilation- and perfusion-weighted maps derived by PREFUL MRI were independent of gadolinium-based contrast agents, which is important for the design of MRI protocols in future studies.

Free-breathing quantification of regional ventilation derived by phase-resolved functional lung (PREFUL) MRI

PURPOSE

To test the feasibility of regional fully quantitative ventilation measurement in free breathing derived by phase-resolved functional lung (PREFUL) MRI in the supine and prone positions. In addition, the influence of T₂ * relaxation time on ventilation quantification is assessed.

METHODS

Twelve healthy volunteers underwent functional MRI at 1.5 T using a 2D triple-echo spoiled gradient echo sequence allowing for quantitative measurement of T₂ * relaxation time. Minute ventilation (ΔV) was quantified by conventional fractional ventilation (FV) and the newly introduced regional ventilation (VR), which corrects volume errors due to image registration. ΔV_FV versus ΔV_VR and ΔV_VR versus ΔV_VR with T₂ * correction were compared using Bland-Altman plots and correlation analysis. The repeatability and physiological plausibility of all measurements were tested in the supine and prone positions.

RESULTS

On global and regional scales a strong correlation was observed between ΔV_FV versus ΔV_VR and ΔV_VR versus ΔV_VRT2* (r > 0.93); however, regional Bland-Altman analysis showed systematic differences (p < 0.0001). Unlike ΔV_VRT2* , ΔV_VR and ΔV_FV showed expected physiologic anterior-posterior gradients, which decreased in the supine but not in the prone position at second measurement during 3 min in the same position. For all quantification methods a moderate repeatability (coefficient of variation <20%) of ventilation was found.

CONCLUSION

A fully quantified regional ventilation measurement using ΔV_VR in free breathing is feasible and shows physiologically plausible results. In contrast to conventional ΔV_FV, volume errors due to image registration are eliminated with the ΔV_VR approach. However, correction for the T₂ * effect remains challenging.

Nonuniform Fourier-decomposition MRI for ventilation- and perfusion-weighted imaging of the lung

PURPOSE

To improve the robustness of pulmonary ventilation- and perfusion-weighted imaging with Fourier decomposition (FD) MRI in the presence of respiratory and cardiac frequency variations by replacing the standard fast Fourier transform with the more general nonuniform Fourier transform.

THEORY AND METHODS

Dynamic coronal single-slice MRI of the thorax was performed in 11 patients and 5 healthy volunteers on a 1.5T whole-body scanner using a 2D ultra-fast balanced steady-state free-precession sequence with temporal resolutions of 4-9 images/s. For the proposed nonuniform Fourier-decomposition (NUFD) approach, the original signal with variable physiological frequencies that was acquired with constant sampling rate was retrospectively transformed into a signal with (ventilation or perfusion) frequency-adapted sampling rate. For that purpose, frequency tracking was performed with the synchro-squeezed wavelet transform. Ventilation- and perfusion-weighted NUFD amplitude and signal delay maps were generated and quantitatively compared with regularly sampled FD maps based on their signal-to-noise ratio (SNR).

RESULTS

Volunteers and patients showed statistically significant increases of SNR in frequency-adapted NUFD results compared to regularly sampled FD results. For ventilation data, the mean SNR increased by 43.4 % ± 25.3 % and 24.4 % ± 31.9 % in volunteers and patients, respectively; for perfusion data, SNR increased by 93.0 % ± 36.1 % and 75.6 % ± 62.8 % . Two patients showed perfusion signal in pulmonary areas with NUFD that could not be imaged with FD.

CONCLUSION

This study demonstrates that using nonuniform Fourier transform in combination with frequency tracking can significantly increase SNR and reduce frequency overlaps by collecting the signal intensity onto single frequency bins.

A framework for Fourier-decomposition free-breathing pulmonary 1 H MRI ventilation measurements

PURPOSE

To develop a rapid Fourier decomposition (FD) free-breathing pulmonary ¹ H MRI (FDMRI) image processing and biomarker pipeline for research use.

METHODS

We acquired MRI in 20 asthmatic subjects using a balanced steady-state free precession (bSSFP) sequence optimized for ventilation imaging. 2D ¹ H MRI series were segmented by enforcing the spatial similarity between adjacent images and the right-to-left lung volume-ratio. The segmented lung series were co-registered using a coarse-to-fine deformable registration framework that used dual optimization techniques. All pairwise registrations were implemented in parallel and FD was performed to generate 2D ventilation-weighted maps and ventilation-defect-percent (VDP). Lung segmentation and registration accuracy were evaluated by comparing algorithm and manual lung-masks, deformed manual lung-masks, and fiducials in the moving and fixed images using Dice-similarity-coefficient (DSC), mean-absolute-distance (MAD), and target-registration-error (TRE). The relationship of FD-VDP and ³ He-VDP was evaluated using the Pearson-correlation-coefficient (r) and Bland Altman analysis. Algorithm reproducibility was evaluated using the coefficient-of-variation (CoV) and intra-class-correlation-coefficient (ICC) for segmentation, registration, and FD-VDP components.

RESULTS

For lung segmentation, there was a DSC of 95 ± 1.5% and MAD of 2.3 ± 0.5 mm, and for registration there was a DSC of 97 ± 0.8%, MAD of 1.6 ± 0.4 mm and TRE of 3.6 ± 1.2 mm. Reproducibility for segmentation DSC (CoV/ICC = 0.5%/0.92), registration TRE (CoV/ICC = 0.4%/0.98), and FD-VDP (Cov/ICC = 3.9%/0.97) was high. The pipeline required 10 min/subject. FD-VDP was correlated with ³ He-VDP (r = 0.69, P < 0.001) although there was a bias toward lower FD-VDP (bias = -4.9%).

CONCLUSIONS

We developed and evaluated a pipeline that provides a rapid and precise method for FDMRI ventilation maps.

Free-breathing Functional Pulmonary MRI: Response to Bronchodilator and Bronchoprovocation in Severe Asthma

RATIONALE AND OBJECTIVES

Ventilation heterogeneity is a hallmark feature of asthma. Our objective was to evaluate ventilation heterogeneity in patients with severe asthma, both pre- and post-salbutamol, as well as post-methacholine (MCh) challenge using the lung clearance index, free-breathing pulmonary ¹H magnetic resonance imaging (FDMRI), and inhaled-gas MRI ventilation defect percent (VDP).

MATERIALS AND METHODS

Sixteen severe asthmatics (49 ± 10 years) provided written informed consent to an ethics board-approved protocol. Spirometry, plethysmography, and multiple breath nitrogen washout to measure the lung clearance index were performed during a single visit within 15 minutes of MRI. Inhaled-gas MRI and FDMRI were performed pre- and post-bronchodilator to generate VDP. For asthmatics with forced expiratory volume in 1 second (FEV₁) >70%_predicted, MRI was also performed before and after MCh challenge. Wilcoxon signed-rank tests, Spearman correlations, and a repeated-measures analysis of variance were performed.

RESULTS

Hyperpolarized ³He (P = .02) and FDMRI (P = .02) VDP significantly improved post-salbutamol and for four asthmatics who could perform MCh (n = 4). ³He and FDMRI VDP significantly increased at the provocative concentration of MCh, resulting in a 20% decrease in FEV₁ (PC₂₀) and decreased post-bronchodilator (P = .02), with a significant difference between methods (P = .01). FDMRI VDP was moderately correlated with ³He VDP (ρ = .61, P = .01), but underestimated VDP relative to ³He VDP (-6 ± 9%). Whereas ³He MRI VDP was significantly correlated with the lung clearance index, FDMRI was not (ρ = .49, P = .06).

CONCLUSIONS

FDMRI VDP generated in free-breathing asthmatic patients was correlated with static inspiratory breath-hold ³He MRI VDP but underestimated VDP relative to ³He MRI VDP. Although less sensitive to salbutamol and MCh, FDMRI VDP may be considered for asthma patient evaluations at centers without inhaled-gas MRI.

Feasibility of quantitative regional ventilation and perfusion mapping with phase-resolved functional lung (PREFUL) MRI in healthy volunteers and COPD, CTEPH, and CF patients

PURPOSE

In this feasibility study, a phase-resolved functional lung imaging postprocessing method for extraction of dynamic perfusion (Q) and ventilation (V) parameters using a conventional 1H lung MRI Fourier decomposition acquisition is introduced.

METHODS

Time series of coronal gradient-echo MR images with a temporal resolution of 288 to 324 ms of two healthy volunteers, one patient with chronic thromboembolic hypertension, one patient with cystic fibrosis, and one patient with chronic obstructive pulmonary disease were acquired at 1.5 T. Using a sine model to estimate cardiac and respiratory phases of each image, all images were sorted to reconstruct full cardiac and respiratory cycles. Time to peak (TTP), V/Q maps, and fractional ventilation flow-volume loops were calculated.

RESULTS

For the volunteers, homogenous ventilation and perfusion TTP maps (V-TTP, Q-TTP) were obtained. The chronic thromboembolic hypertension patient showed increased perfusion TTP in hypoperfused regions in visual agreement with dynamic contrast-enhanced MRI, which improved postpulmonary endaterectomy surgery. Cystic fibrosis and chronic obstructive pulmonary disease patients showed a pattern of increased V-TTP and Q-TTP in regions of hypoventilation and decreased perfusion. Fractional ventilation flow-volume loops of the chronic obstructive pulmonary disease patient were smaller in comparison with the healthy volunteer, and showed regional differences in visual agreement with functional small airways disease and emphysema on CT.

CONCLUSIONS

This study shows the feasibility of phase-resolved functional lung imaging to gain quantitative information regarding regional lung perfusion and ventilation without the need for ultrafast imaging, which will be advantageous for future clinical translation. Magn Reson Med 79:2306-2314, 2018. © 2017 International Society for Magnetic Resonance in Medicine.

Robust power spectral estimation for EEG data

BACKGROUND

Typical electroencephalogram (EEG) recordings often contain substantial artifact. These artifacts, often large and intermittent, can interfere with quantification of the EEG via its power spectrum. To reduce the impact of artifact, EEG records are typically cleaned by a preprocessing stage that removes individual segments or components of the recording. However, such preprocessing can introduce bias, discard available signal, and be labor-intensive. With this motivation, we present a method that uses robust statistics to reduce dependence on preprocessing by minimizing the effect of large intermittent outliers on the spectral estimates.

NEW METHOD

Using the multitaper method (Thomson, 1982) as a starting point, we replaced the final step of the standard power spectrum calculation with a quantile-based estimator, and the Jackknife approach to confidence intervals with a Bayesian approach. The method is implemented in provided MATLAB modules, which extend the widely used Chronux toolbox.

RESULTS

Using both simulated and human data, we show that in the presence of large intermittent outliers, the robust method produces improved estimates of the power spectrum, and that the Bayesian confidence intervals yield close-to-veridical coverage factors.

COMPARISON TO EXISTING METHOD

The robust method, as compared to the standard method, is less affected by artifact: inclusion of outliers produces fewer changes in the shape of the power spectrum as well as in the coverage factor.

CONCLUSION

In the presence of large intermittent outliers, the robust method can reduce dependence on data preprocessing as compared to standard methods of spectral estimation.

Matrix pencil decomposition of time-resolved proton MRI for robust and improved assessment of pulmonary ventilation and perfusion

PURPOSE

To present an improved and robust method of pulmonary function assessment from time-resolved proton MRI using a matrix pencil (MP) method in combination with a linear least squares analysis.

METHODS

Simulations of the signal time course in lung parenchyma were performed to compare the accuracy of Fourier decomposition (FD) and MP methods for the estimation of respiratory and cardiac amplitudes. Series of two-dimensional time-resolved lung images were acquired in healthy volunteers at 1.5 T using ultra-fast steady-state free precession. Qualitative lung ventilation- and perfusion-weighted images as well as a quantitative map of fractional ventilation, perfusion, and blood arrival time were calculated using the proposed MP method and compared with the contemporary FD technique. A region-of-interest analysis was performed on the quantitative data.

RESULTS

The signal analysis performed using MP decomposition resulted in reduced variability of the estimated respiratory and cardiac amplitudes in comparison with FD for both simulated and in vivo data.

CONCLUSION

MP decomposition provides an automatic, robust, and more accurate estimation of amplitudes of respiratory and cardiac signal modulations in the lung parenchyma than the contemporary FD technique. Magn Reson Med 77:336-342, 2017. © 2016 Wiley Periodicals, Inc.

Reproducibility of fractional ventilation derived by Fourier decomposition after adjusting for tidal volume with and without an MRI compatible spirometer

PURPOSE

To reduce the influence of tidal volume on fractional ventilation (FV) derived by Fourier decomposition (FD).

METHODS

Twelve volunteers were examined on a 1.5 Tesla scanner. Spoiled gradient echo imaging of coronal and sagittal slices of the lung were performed. The tidal volume variations between different acquisitions were studied by reproducibility and repeatability measurements. To adjust the FV derived by FD for tidal volume differences between the measurements, during all acquisitions, the lung volume changes were measured by a spirometer and used to calculate a global FV parameter. As an alternative, using the FD data, the lung area changes were calculated and used for the adjustment.

RESULTS

Reproducibility analysis of unadjusted coronal FV showed a determination coefficient of R2 = 71% and an intraclass correlation coefficient of ICC = 93%. Differences in the measurements could be ascribed to different tidal volumes. Area adjusted values exhibited an increased R2 of 84% and a higher ICC of 97%. For the coronal middle slice/sagittal slices in free breathing, the inter-volunteer coefficient of variation was reduced from 0.23/0.28 (unadjusted) to 0.16/0.20 (spirometer) or 0.12/0.13 (area).

CONCLUSION

The calculation of lung area changes is sufficient to increase the reproducibility of FV in a volunteer cohort avoiding the need for an MRI compatible spirometer. Magn Reson Med 76:1542-1550, 2016. © 2015 International Society for Magnetic Resonance in Medicine.

Free-breathing pulmonary 1H and Hyperpolarized 3He MRI: comparison in COPD and bronchiectasis

RATIONALE AND OBJECTIVES

In this proof-of-concept demonstration, we aimed to quantitatively and qualitatively compare pulmonary ventilation abnormalities derived from Fourier decomposition of free-breathing (1)H magnetic resonance imaging (FDMRI) to hyperpolarized (3)He MRI in subjects with chronic obstructive pulmonary disease (COPD) and bronchiectasis.

MATERIALS AND METHODS

All subjects provided written informed consent to a protocol approved by a local research ethics board and Health, Canada, and they underwent MRI, computed tomography (CT), spirometry, and plethysmography during a single 2-hour visit. Semiautomated segmentation was used to generate ventilation defect measurements derived from FDMRI and (3)He MRI, and these were compared using analysis of variance and Pearson correlations.

RESULTS

Twenty-six subjects were evaluated including 12 COPD subjects (67 ± 9 years) and 14 bronchiectasis subjects (70 ± 11 years). For COPD subjects, FDMRI and (3)He MRI ventilation defect percent (VDP) was 7 ± 6% and 24 ± 14%, respectively (P < .001; bias = -16 ± 9%). In COPD subjects, FDMRI was significantly correlated with (3)He MRI VDP (r = .88; P = .0001), (3)He MRI apparent diffusion coefficient (r = .71; P < .05), airways resistance (r = .60; P < .05), and RA950 (r = .80; P < .01). In subjects with bronchiectasis, FDMRI VDP (5 ± 3%) and (3)He MRI VDP (18 ± 9%) were significantly different (P < .001) and not correlated (P > .05). The Dice similarity coefficient (DSC) for FDMRI and (3)He MRI ventilation was 86 ± 7% for COPD and 86 ± 4% for bronchiectasis subjects (P > .05); the DSC for FDMRI ventilation defects and CT RA950 was 19 ± 20% in COPD and 2 ± 3% in bronchiectasis subjects (P < .01).

CONCLUSIONS

FDMRI and (3)He MRI VDP were strongly related in COPD but not in bronchiectasis subjects. In COPD only, FDMRI ventilation defects were spatially related with (3)He ventilation defects and emphysema.

Performance of perfusion-weighted Fourier decomposition MRI for detection of chronic pulmonary emboli

PURPOSE

To evaluate the test performance of perfusion-weighted Fourier-decomposition (pw-FD) magnetic resonance imaging (MRI) in comparison to dynamic contrast-enhanced (DCE)-MRI as a reference standard in patients with known or suspected chronic pulmonary embolism (PE).

MATERIALS AND METHODS

In 64 patients, chronic PE was ruled out or confirmed by DCE-MRI using a time-resolved angiography with stochastic trajectories (TWIST) sequence in one breath-hold. Pw-FD-MRI was performed using a 2D fast low-angle shot (FLASH) sequence in free-breathing. After a nonrigid image registration, FD was applied to generate pw-images. Lungs were scored by two radiologists (2 and 12 years of lung MRI experience) visually for each lobe and segment for hypoperfused areas. For intra- and interobserver variability, the MR images were analyzed 2 months after the first analysis, blinded to the results of the first reader.

RESULTS

PE was diagnosed by DCE-MRI in 39 patients. For the pw-FD MRI sensitivity, specificity, accuracy, and positive and negative predictive value for diagnosis of PE were 100%, 95%, 98%, 98%, and 100% on a per-patient basis, 94%, 94%, 94%, 95%, 94% on a per-lobe basis, and 82%, 92%, 88%, 88%, 88% on a segmental basis, respectively. Detection of subsegmental and segmental hypoperfusion using pw-FD MRI showed a moderate agreement with DCE-MRI (kappa of 0.68; 95% confidence interval: 0.64; 0.72).

CONCLUSION

Pw-FD of the lung is a feasible test to diagnose chronic PE on a per-patient level during free-breathing without the use of ionizing radiation or contrast agents.

Exact results for amplitude spectra of fitness landscapes

Starting from fitness correlation functions, we calculate exact expressions for the amplitude spectra of fitness landscapes as defined by Stadler [1996. Landscapes and their correlation functions. J. Math. Chem. 20, 1] for common landscape models, including Kauffman's NK-model, rough Mount Fuji landscapes and general linear superpositions of such landscapes. We further show that correlations decaying exponentially with the Hamming distance yield exponentially decaying spectra similar to those reported recently for a model of molecular signal transduction. Finally, we compare our results for the model systems to the spectra of various experimentally measured fitness landscapes. We claim that our analytical results should be helpful when trying to interpret empirical data and guide the search for improved fitness landscape models.