High-resolution 3D whole-heart coronary MRA: a study on the combination of data acquisition in multiple breath-holds and 1D residual respiratory motion compensation.

OBJECT: To study a scan protocol for coronary magnetic resonance angiography based on multiple breath-holds featuring 1D motion compensation and to compare the resulting image quality to a navigator-gated free-breathing acquisition. Image reconstruction was performed using L1 regularized iterative SENSE. MATERIALS AND METHODS: The effects of respiratory motion on the Cartesian sampling scheme were minimized by performing data acquisition in multiple breath-holds. During the scan, repetitive readouts through a k-space center were used to detect and correct the respiratory displacement of the heart by exploiting the self-navigation principle in image reconstruction. In vivo experiments were performed in nine healthy volunteers and the resulting image quality was compared to a navigator-gated reference in terms of vessel length and sharpness. RESULTS: Acquisition in breath-hold is an effective method to reduce the scan time by more than 30 % compared to the navigator-gated reference. Although an equivalent mean image quality with respect to the reference was achieved with the proposed method, the 1D motion compensation did not work equally well in all cases. CONCLUSION: In general, the image quality scaled with the robustness of the motion compensation. Nevertheless, the featured setup provides a positive basis for future extension with more advanced motion compensation methods.

Compressed Sensing Single-Breath-Hold CMR for Fast Quantification of LV Function, Volumes, and Mass.

OBJECTIVES: The purpose of this study was to compare a novel compressed sensing (CS)-based single-breath-hold multislice magnetic resonance cine technique with the standard multi-breath-hold technique for the assessment of left ventricular (LV) volumes and function. BACKGROUND: Cardiac magnetic resonance is generally accepted as the gold standard for LV volume and function assessment. LV function is 1 of the most important cardiac parameters for diagnosis and the monitoring of treatment effects. Recently, CS techniques have emerged as a means to accelerate data acquisition. METHODS: The prototype CS cine sequence acquires 3 long-axis and 4 short-axis cine loops in 1 single breath-hold (temporal/spatial resolution: 30 ms/1.5 × 1.5 mm(2); acceleration factor 11.0) to measure left ventricular ejection fraction (LVEFCS) as well as LV volumes and LV mass using LV model-based 4D software. For comparison, a conventional stack of multi-breath-hold cine images was acquired (temporal/spatial resolution 40 ms/1.2 × 1.6 mm(2)). As a reference for the left ventricular stroke volume (LVSV), aortic flow was measured by phase-contrast acquisition. RESULTS: In 94% of the 33 participants (12 volunteers: mean age 33 ± 7 years; 21 patients: mean age 63 ± 13 years with different LV pathologies), the image quality of the CS acquisitions was excellent. LVEFCS and LVEFstandard were similar (48.5 ± 15.9% vs. 49.8 ± 15.8%; p = 0.11; r = 0.96; slope 0.97; p < 0.00001). Agreement of LVSVCS with aortic flow was superior to that of LVSVstandard (overestimation vs. aortic flow: 5.6 ± 6.5 ml vs. 16.2 ± 11.7 ml, respectively; p = 0.012) with less variability (r = 0.91; p < 0.00001 for the CS technique vs. r = 0.71; p < 0.01 for the standard technique). The intraobserver and interobserver agreement for all CS parameters was good (slopes 0.93 to 1.06; r = 0.90 to 0.99). CONCLUSIONS: The results demonstrated the feasibility of applying the CS strategy to evaluate LV function and volumes with high accuracy in patients. The single-breath-hold CS strategy has the potential to replace the multi-breath-hold standard cardiac magnetic resonance technique.

A comparison of methods for differential expression analysis of RNA-seq data.

BACKGROUND: Finding genes that are differentially expressed between conditions is an integral part of understanding the molecular basis of phenotypic variation. In the past decades, DNA microarrays have been used extensively to quantify the abundance of mRNA corresponding to different genes, and more recently high-throughput sequencing of cDNA (RNA-seq) has emerged as a powerful competitor. As the cost of sequencing decreases, it is conceivable that the use of RNA-seq for differential expression analysis will increase rapidly. To exploit the possibilities and address the challenges posed by this relatively new type of data, a number of software packages have been developed especially for differential expression analysis of RNA-seq data. RESULTS: We conducted an extensive comparison of eleven methods for differential expression analysis of RNA-seq data. All methods are freely available within the R framework and take as input a matrix of counts, i.e. the number of reads mapping to each genomic feature of interest in each of a number of samples. We evaluate the methods based on both simulated data and real RNA-seq data. CONCLUSIONS: Very small sample sizes, which are still common in RNA-seq experiments, impose problems for all evaluated methods and any results obtained under such conditions should be interpreted with caution. For larger sample sizes, the methods combining a variance-stabilizing transformation with the 'limma' method for differential expression analysis perform well under many different conditions, as does the nonparametric SAMseq method.

Data-driven topo-climatic mapping with machine learning methods

Automatic environmental monitoring networks enforced by wireless communication technologies provide large and ever increasing volumes of data nowadays. The use of this information in natural hazard research is an important issue. Particularly useful for risk assessment and decision making are the spatial maps of hazard-related parameters produced from point observations and available auxiliary information. The purpose of this article is to present and explore the appropriate tools to process large amounts of available data and produce predictions at fine spatial scales. These are the algorithms of machine learning, which are aimed at non-parametric robust modelling of non-linear dependencies from empirical data. The computational efficiency of the data-driven methods allows producing the prediction maps in real time which makes them superior to physical models for the operational use in risk assessment and mitigation. Particularly, this situation encounters in spatial prediction of climatic variables (topo-climatic mapping). In complex topographies of the mountainous regions, the meteorological processes are highly influenced by the relief. The article shows how these relations, possibly regionalized and non-linear, can be modelled from data using the information from digital elevation models. The particular illustration of the developed methodology concerns the mapping of temperatures (including the situations of Föhn and temperature inversion) given the measurements taken from the Swiss meteorological monitoring network. The range of the methods used in the study includes data-driven feature selection, support vector algorithms and artificial neural networks.

Handheld-assisted field data collection for occupational risk assessment.

A computerized handheld procedure is presented in this paper. It is intended as a database complementary tool, to enhance prospective risk analysis in the field of occupational health. The Pendragon forms software (version 3.2) has been used to implement acquisition procedures on Personal Digital Assistants (PDAs) and to transfer data to a computer in an MS-Access format. The data acquisition strategy proposed relies on the risk assessment method practiced at the Institute of Occupational Health Sciences (IST). It involves the use of a systematic hazard list and semi-quantitative risk assessment scales. A set of 7 modular forms has been developed to cover the basic need of field audits. Despite the minor drawbacks observed, the results obtained so far show that handhelds are adequate to support field risk assessment and follow-up activities. Further improvements must still be made in order to increase the tool effectiveness and field adequacy.

Improved temporal resolution for functional studies with reduced number of segments with three-dimensional echo planar imaging.

PURPOSE: To introduce a new k-space traversal strategy for segmented three-dimensional echo planar imaging (3D EPI) that encodes two partitions per radiofrequency excitation, effectively reducing the number excitations used to acquire a 3D EPI dataset by half. METHODS: The strategy was evaluated in the context of functional MRI applications for: image quality compared with segmented 3D EPI, temporal signal-to-noise ratio (tSNR) (the ability to detect resting state networks compared with multislice two-dimensional (2D) EPI and segmented 3D EPI, and temporal resolution (the ability to separate cardiac- and respiration-related fluctuations from the desired blood oxygen level-dependent signal of interest). RESULTS: Whole brain images with a nominal voxel size of 2 mm isotropic could be acquired with a temporal resolution under half a second using traditional parallel imaging acceleration up to 4× in the partition-encode direction and using novel data acquisition speed-up of 2× with a 32-channel coil. With 8× data acquisition speed-up in the partition-encode direction, 3D reduced excitations (RE)-EPI produced acceptable image quality without introduction of noticeable additional artifacts. Due to increased tSNR and better characterization of physiological fluctuations, the new strategy allowed detection of more resting state networks compared with multislice 2D-EPI and segmented 3D EPI. CONCLUSION: 3D RE-EPI resulted in significant increases in temporal resolution for whole brain acquisitions and in improved physiological noise characterization compared with 2D-EPI and segmented 3D EPI. Magn Reson Med 72:786-792, 2014. © 2013 Wiley Periodicals, Inc.

Human Primary Auditory Cortex Follows the Shape of Heschl's Gyrus.

The primary auditory cortex (PAC) is central to human auditory abilities, yet its location in the brain remains unclear. We measured the two largest tonotopic subfields of PAC (hA1 and hR) using high-resolution functional MRI at 7 T relative to the underlying anatomy of Heschl's gyrus (HG) in 10 individual human subjects. The data reveals a clear anatomical-functional relationship that, for the first time, indicates the location of PAC across the range of common morphological variants of HG (single gyri, partial duplications, and complete duplications). In 20/20 individual hemispheres, two primary mirror-symmetric tonotopic maps were clearly observed with gradients perpendicular to HG. PAC spanned both divisions of HG in cases of partial and complete duplications (11/20 hemispheres), not only the anterior division as commonly assumed. Specifically, the central union of the two primary maps (the hA1-R border) was consistently centered on the full Heschl's structure: on the gyral crown of single HGs and within the sulcal divide of duplicated HGs. The anatomical-functional variants of PAC appear to be part of a continuum, rather than distinct subtypes. These findings significantly revise HG as a marker for human PAC and suggest that tonotopic maps may have shaped HG during human evolution. Tonotopic mappings were based on only 16 min of fMRI data acquisition, so these methods can be used as an initial mapping step in future experiments designed to probe the function of specific auditory fields.

Repositioning precision of coronary arteries measured on X-ray angiography and its implications for coronary MR angiography.

BACKGROUND: To test the hypothesis that intervals with superior beat-to-beat coronary artery repositioning precision exist in the cardiac cycle, to design a coronary MR angiography (MRA) methodology in response, and to ascertain its performance. METHODS: Coronary repositioning precision in consecutive heartbeats was measured on x-ray coronary angiograms of 17 patients and periods with the highest repositioning precision were identified. In response, the temporal order of coronary MRA pulse sequence elements required modification and the T2 -prep now follows (T2 -post) rather than precedes the imaging part of the sequence. The performance of T2 -post was quantitatively compared (signal-to-noise [SNR], contrast-to-noise [CNR], vessel sharpness) to that of T2 -prep in vivo. RESULTS: Coronary repositioning precision is <1 mm at peak systole and in mid diastole. When comparing systolic T2 -post to diastolic T2 -prep, CNR and vessel sharpness remained unchanged (both P = NS) but SNR for muscle and blood increased by 104% and 36% (both P < 0.05), respectively. CONCLUSION: Windows with improved coronary repositioning precision exist in the cardiac cycle: one in peak systole and one in mid diastole. Peak-systolic imaging necessitates a re-design of conventional coronary MRA pulse sequences and leads to image quality very similar to that of conventional mid-diastolic data acquisition but improved SNR. J. Magn. Reson. Imaging 2015;41:1251-1258. © 2014 Wiley Periodicals, Inc.

Simultaneous B(0)- and B(1)+-map acquisition for fast localized shim, frequency, and RF power determination in the heart at 3 T.

High-field (>or=3 T) cardiac MRI is challenged by inhomogeneities of both the static magnetic field (B(0)) and the transmit radiofrequency field (B(1)+). The inhomogeneous B fields not only demand improved shimming methods but also impede the correct determination of the zero-order terms, i.e., the local resonance frequency f(0) and the radiofrequency power to generate the intended local B(1)+ field. In this work, dual echo time B(0)-map and dual flip angle B(1)+-map acquisition methods are combined to acquire multislice B(0)- and B(1)+-maps simultaneously covering the entire heart in a single breath hold of 18 heartbeats. A previously proposed excitation pulse shape dependent slice profile correction is tested and applied to reduce systematic errors of the multislice B(1)+-map. Localized higher-order shim correction values including the zero-order terms for frequency f(0) and radiofrequency power can be determined based on the acquired B(0)- and B(1)+-maps. This method has been tested in 7 healthy adult human subjects at 3 T and improved the B(0) field homogeneity (standard deviation) from 60 Hz to 35 Hz and the average B(1)+ field from 77% to 100% of the desired B(1)+ field when compared to more commonly used preparation methods.

Respiratory motion artifact suppression in diffusion-weighted MR imaging of the spine.

Diffusion-weighted spin-echo imaging of the spine has been successfully implemented for differentiation of benign fracture edema and tumor infiltration of the vertebral body. Nevertheless, this technique still suffers from insufficient image quality in numerous patients due to motion artifacts. The aim of this study was to investigate the impact of variable respiratory motion artifact suppression techniques on image quality in diffusion-weighted spin-echo imaging of the spine. In addition to phase-encoding reordering, a newly implemented right hemi-diaphragmaitc navigator for respiratory gating was used. Subjective and objective image quality parameters were compared. Respiratory motion artifact suppression has a major impact on image quality in diffusion-weighted imaging of the spine. Phase-encoding reordering does not enhance image quality while right hemi-diaphragmatic respiratory navigator gating significantly improves image quality at the cost of data acquisition time. Navigator gating should be used if standard spin-echo diffusion-weighted imaging demonstrates insufficient image quality.

Self-Navigation with Compressed Sensing for 2D Translational Motion Correction in Free-Breathing Coronary MRI: A Feasibility Study.

PURPOSE: Respiratory motion correction remains a challenge in coronary magnetic resonance imaging (MRI) and current techniques, such as navigator gating, suffer from sub-optimal scan efficiency and ease-of-use. To overcome these limitations, an image-based self-navigation technique is proposed that uses "sub-images" and compressed sensing (CS) to obtain translational motion correction in 2D. The method was preliminarily implemented as a 2D technique and tested for feasibility for targeted coronary imaging. METHODS: During a 2D segmented radial k-space data acquisition, heavily undersampled sub-images were reconstructed from the readouts collected during each cardiac cycle. These sub-images may then be used for respiratory self-navigation. Alternatively, a CS reconstruction may be used to create these sub-images, so as to partially compensate for the heavy undersampling. Both approaches were quantitatively assessed using simulations and in vivo studies, and the resulting self-navigation strategies were then compared to conventional navigator gating. RESULTS: Sub-images reconstructed using CS showed a lower artifact level than sub-images reconstructed without CS. As a result, the final image quality was significantly better when using CS-assisted self-navigation as opposed to the non-CS approach. Moreover, while both self-navigation techniques led to a 69% scan time reduction (as compared to navigator gating), there was no significant difference in image quality between the CS-assisted self-navigation technique and conventional navigator gating, despite the significant decrease in scan time. CONCLUSIONS: CS-assisted self-navigation using 2D translational motion correction demonstrated feasibility of producing coronary MRA data with image quality comparable to that obtained with conventional navigator gating, and does so without the use of additional acquisitions or motion modeling, while still allowing for 100% scan efficiency and an improved ease-of-use. In conclusion, compressed sensing may become a critical adjunct for 2D translational motion correction in free-breathing cardiac imaging with high spatial resolution. An expansion to modern 3D approaches is now warranted.

Precursor ion scans for the targeted detection of stable-isotope-labeled peptides.

Stable isotope labels are routinely introduced into proteomes for quantification purposes. Full labeling of cells in varying biological states, followed by sample mixing, fractionation and intensive data acquisition, is used to obtain accurate large-scale quantification of total protein levels. However, biological processes often affect only a small group of proteins for a short time, resulting in changes that are difficult to detect against the total proteome background. An alternative approach could be the targeted analysis of the proteins synthesized in response to a given biological stimulus. Such proteins can be pulse-labeled with a stable isotope by metabolic incorporation of 'heavy' amino acids. In this study we investigated the specific detection and identification of labeled proteins using acquisition methods based on Precursor Ion Scans (PIS) on a triple-quadrupole ion trap mass spectrometer. PIS-based methods were set to detect unique immonium ions originating from labeled peptides. Different labels and methods were tested in standard mixtures to optimize performance. We showed that, in comparison with an untargeted analysis on the same instrument, the approach allowed a several-fold increase in the specificity of detection of labeled proteins over unlabeled ones. The technique was applied to the identification of proteins secreted by human cells into growth media containing bovine serum proteins, allowing the preferential detection of labeled cellular proteins over unlabeled bovine ones. However, compared with untargeted acquisitions on two different instruments, the PIS-based strategy showed some limitations in sensitivity. We discuss possible perspectives of the technique.

Practical signal-to-noise ratio quantification for sensitivity encoding: Application to coronary MR angiography.

Purpose: To develop and evaluate a practical method for the quantification of signal-to-noise ratio (SNR) on coronary MR angiograms (MRA) acquired with parallel imaging.Materials and Methods: To quantify the spatially varying noise due to parallel imaging reconstruction, a new method has been implemented incorporating image data acquisition followed by a fast noise scan during which radio-frequency pulses, cardiac triggering and navigator gating are disabled. The performance of this method was evaluated in a phantom study where SNR measurements were compared with those of a reference standard (multiple repetitions). Subsequently, SNR of myocardium and posterior skeletal muscle was determined on in vivo human coronary MRA.Results: In a phantom, the SNR measured using the proposed method deviated less than 10.1% from the reference method for small geometry factors (<= 2). In vivo, the noise scan for a 10 min coronary MRA acquisition was acquired in 30 s. Higher signal and lower SNR, due to spatially varying noise, were found in myocardium compared with posterior skeletal muscle.Conclusion: SNR quantification based on a fast noise scan is a validated and easy-to-use method when applied to three-dimensional coronary MRA obtained with parallel imaging as long as the geometry factor remains low.

Free-breathing renal magnetic resonance angiography with steady-state free-precession and slab-selective spin inversion combined with radial k-space sampling and water-selective excitation.

The impact of radial k-space sampling and water-selective excitation on a novel navigator-gated cardiac-triggered slab-selective inversion prepared 3D steady-state free-precession (SSFP) renal MR angiography (MRA) sequence was investigated. Renal MRA was performed on a 1.5-T MR system using three inversion prepared SSFP approaches: Cartesian (TR/TE: 5.7/2.8 ms, FA: 85 degrees), radial (TR/TE: 5.5/2.7 ms, FA: 85 degrees) SSFP, and radial SSFP combined with water-selective excitation (TR/TE: 9.9/4.9 ms, FA: 85 degrees). Radial data acquisition lead to significantly reduced motion artifacts (P < 0.05). SNR and CNR were best using Cartesian SSFP (P < 0.05). Vessel sharpness and vessel length were comparable in all sequences. The addition of a water-selective excitation could not improve image quality. In conclusion, radial k-space sampling reduces motion artifacts significantly in slab-selective inversion prepared renal MRA, while SNR and CNR are decreased. The addition of water-selective excitation could not improve the lower CNR in radial scanning.