Highly accelerated non-contrast-enhanced time-resolved 4D MRA using stack-of-stars golden-angle radial acquisition with a self-calibrated low-rank subspace reconstruction.

PURPOSE: To develop a highly accelerated non-contrast-enhanced 4D-MRA technique by combining stack-of-stars golden-angle radial acquisition with a modified self-calibrated low-rank subspace reconstruction. METHODS: A low-rank subspace reconstruction framework was introduced in radial 4D MRA (SUPER 4D MRA) by combining stack-of-stars golden-angle radial acquisition with control-label k-space subtraction-based low-rank subspace modeling. Radial 4D MRA data were acquired and reconstructed using the proposed technique on 12 healthy volunteers and 1 patient with steno-occlusive disease. The performance of SUPER 4D MRA was compared with two temporally constrained reconstruction methods (golden-angle radial sparse parallel [GRASP] and GRASP-Pro) at different acceleration rates in terms of image quality and delineation of blood dynamics. RESULTS: SUPER 4D MRA outperformed the other two reconstruction methods, offering superior image quality with a clear background and detailed delineation of cerebrovascular structures as well as great temporal fidelity in blood flow dynamics. SUPER 4D MRA maintained excellent performance even at higher acceleration rates. CONCLUSIONS: SUPER 4D MRA is a promising technique for highly accelerating 4D MRA acquisition without comprising both temporal fidelity and image quality.

Self-calibrated subspace reconstruction for multidimensional MR fingerprinting for simultaneous relaxation and diffusion quantification.

PURPOSE: To propose a new reconstruction method for multidimensional MR fingerprinting (mdMRF) to address shading artifacts caused by physiological motion-induced measurement errors without navigating or gating. METHODS: The proposed method comprises two procedures: self-calibration and subspace reconstruction. The first procedure (self-calibration) applies temporally local matrix completion to reconstruct low-resolution images from a subset of under-sampled data extracted from the k-space center. The second procedure (subspace reconstruction) utilizes temporally global subspace reconstruction with pre-estimated temporal subspace from low-resolution images to reconstruct aliasing-free, high-resolution, and time-resolved images. After reconstruction, a customized outlier detection algorithm was employed to automatically detect and remove images corrupted by measurement errors. Feasibility, robustness, and scan efficiency were evaluated through in vivo human brain imaging experiments. RESULTS: The proposed method successfully reconstructed aliasing-free, high-resolution, and time-resolved images, where the measurement errors were accurately represented. The corrupted images were automatically and robustly detected and removed. Artifact-free T1, T2, and ADC maps were generated simultaneously. The proposed reconstruction method demonstrated robustness across different scanners, parameter settings, and subjects. A high scan efficiency of less than 20 s per slice has been achieved. CONCLUSION: The proposed reconstruction method can effectively alleviate shading artifacts caused by physiological motion-induced measurement errors. It enables simultaneous and artifact-free quantification of T1, T2, and ADC using mdMRF scans without prospective gating, with robustness and high scan efficiency.

Accelerated 3D multi-echo spin-echo sequence with a subspace constrained reconstruction for whole mouse brain T 2 mapping.

PURPOSE: To accelerate whole-brain quantitative T 2 $$ {\mathrm{T}}_2 $$ mapping in preclinical imaging setting. METHODS: A three-dimensional (3D) multi-echo spin echo sequence was highly undersampled with a variable density Poisson distribution to reduce the acquisition time. Advanced iterative reconstruction based on linear subspace constraints was employed to recover high-quality raw images. Different subspaces, generated using exponential or extended-phase graph (EPG) simulations or from low-resolution calibration images, were compared. The subspace dimension was investigated in terms of T 2 $$ {\mathrm{T}}_2 $$ precision. The method was validated on a phantom containing a wide range of T 2 $$ {\mathrm{T}}_2 $$ and was then applied to monitor metastasis growth in the mouse brain at 4.7T. Image quality and T 2 $$ {\mathrm{T}}_2 $$ estimation were assessed for 3 acceleration factors (6/8/10). RESULTS: The EPG-based dictionary gave robust estimations of a large range of T 2 $$ {\mathrm{T}}_2 $$ . A subspace dimension of 6 was the best compromise between T 2 $$ {\mathrm{T}}_2 $$ precision and image quality. Combining the subspace constrained reconstruction with a highly undersampled dataset enabled the acquisition of whole-brain T 2 $$ {\mathrm{T}}_2 $$ maps, the detection and the monitoring of metastasis growth of less than 500 µ m 3 $$ \mu {\mathrm{m}}^3 $$ . CONCLUSION: Subspace-based reconstruction is suitable for 3D T 2 $$ {\mathrm{T}}_2 $$ mapping. This method can be used to reach an acceleration factor up to 8, corresponding to an acquisition time of 25 min for an isotropic 3D acquisition of 156 µ $$ \mu $$ m on the mouse brain, used here for monitoring metastases growth.

Optimized multi-axis spiral projection MR fingerprinting with subspace reconstruction for rapid whole-brain high-isotropic-resolution quantitative imaging.

PURPOSE: To improve image quality and accelerate the acquisition of 3D MR fingerprinting (MRF). METHODS: Building on the multi-axis spiral-projection MRF technique, a subspace reconstruction with locally low-rank constraint and a modified spiral-projection spatiotemporal encoding scheme called tiny golden-angle shuffling were implemented for rapid whole-brain high-resolution quantitative mapping. Reconstruction parameters such as the locally low-rank regularization parameter and the subspace rank were tuned using retrospective in vivo data and simulated examinations. B0 inhomogeneity correction using multifrequency interpolation was incorporated into the subspace reconstruction to further improve the image quality by mitigating blurring caused by off-resonance effect. RESULTS: The proposed MRF acquisition and reconstruction framework yields high-quality 1-mm isotropic whole-brain quantitative maps in 2 min at better quality compared with 6-min acquisitions of prior approaches. The proposed method was validated to not induce bias in T1 and T2 mapping. High-quality whole-brain MRF data were also obtained at 0.66-mm isotropic resolution in 4 min using the proposed technique, where the increased resolution was shown to improve visualization of subtle brain structures. CONCLUSIONS: The proposed tiny golden-angle shuffling, MRF with optimized spiral-projection trajectory and subspace reconstruction enables high-resolution quantitative mapping in ultrafast acquisition time.

Frequency-modulated SSFP with radial sampling and subspace reconstruction: A time-efficient alternative to phase-cycled bSSFP.

PURPOSE: A novel subspace-based reconstruction method for frequency-modulated balanced steady-state free precession (fmSSFP) MRI is presented. In this work, suitable data acquisition schemes, subspace sizes, and efficiencies for banding removal are investigated. THEORY AND METHODS: By combining a fmSSFP MRI sequence with a 3D stack-of-stars trajectory, scan efficiency is maximized as spectral information is obtained without intermediate preparation phases. A memory-efficient reconstruction routine is implemented by introducing the low-frequency Fourier transform as a subspace which allows for the formulation of a convex reconstruction problem. The removal of banding artifacts is investigated by comparing the proposed acquisition and reconstruction technique to phase-cycled bSSFP MRI. Aliasing properties of different undersampling schemes are analyzed and water/fat separation is demonstrated by reweighting the reconstructed subspace coefficients to generate virtual spectral responses in a post-processing step. RESULTS: A simple root-of-sum-of-squares combination of the reconstructed subspace coefficients yields high-SNR images with the characteristic bSSFP contrast but without banding artifacts. Compared to Golden-Angle trajectories, turn-based sampling schemes were superior in minimizing aliasing across reconstructed subspace coefficients. Water/fat separated images of the human knee were obtained by reweighting subspace coefficients. CONCLUSIONS: The novel subspace-based fmSSFP MRI technique emerges as a time-efficient alternative to phase-cycled bSFFP. The method does not need intermediate preparation phases, offers high SNR and avoids banding artifacts. Reweighting of the reconstructed subspace coefficients allows for generating virtual spectral responses with applications to water/fat separation.

Monitoring Pilot's Mental Workload Using ERPs and Spectral Power with a Six-Dry-Electrode EEG System in Real Flight Conditions.

Recent technological progress has allowed the development of low-cost and highly portable brain sensors such as pre-amplified dry-electrodes to measure cognitive activity out of the laboratory. This technology opens promising perspectives to monitor the "brain at work" in complex real-life situations such as while operating aircraft. However, there is a need to benchmark these sensors in real operational conditions. We therefore designed a scenario in which twenty-two pilots equipped with a six-dry-electrode EEG system had to perform one low load and one high load traffic pattern along with a passive auditory oddball. In the low load condition, the participants were monitoring the flight handled by a flight instructor, whereas they were flying the aircraft in the high load condition. At the group level, statistical analyses disclosed higher P300 amplitude for the auditory target (Pz, P4 and Oz electrodes) along with higher alpha band power (Pz electrode), and higher theta band power (Oz electrode) in the low load condition as compared to the high load one. Single trial classification accuracy using both event-related potentials and event-related frequency features at the same time did not exceed chance level to discriminate the two load conditions. However, when considering only the frequency features computed over the continuous signal, classification accuracy reached around 70% on average. This study demonstrates the potential of dry-EEG to monitor cognition in a highly ecological and noisy environment, but also reveals that hardware improvement is still needed before it can be used for everyday flight operations.

Exploring Imaging Genetic Markers of Alzheimer's Disease Based on a Novel Nonlinear Correlation Analysis Algorithm.

Alzheimer's disease (AD) is an irreversible neurological disorder characterized by insidious onset. Identifying potential markers in its emergence and progression is crucial for early diagnosis and treatment. Imaging genetics typically merges genetic variables with multiple imaging parameters, employing various association analysis algorithms to investigate the links between pathological phenotypes and genetic variations, and to unearth molecular-level insights from brain images. However, most existing imaging genetics algorithms based on sparse learning assume a linear relationship between genetic factors and brain functions, limiting their ability to discern complex nonlinear correlation patterns and resulting in reduced accuracy. To address these issues, we propose a novel nonlinear imaging genetic association analysis method, Deep Self-Reconstruction-based Adaptive Sparse Multi-view Deep Generalized Canonical Correlation Analysis (DSR-AdaSMDGCCA). This approach facilitates joint learning of the nonlinear relationships between pathological phenotypes and genetic variations by integrating three different types of data: structural magnetic resonance imaging (sMRI), single-nucleotide polymorphism (SNP), and gene expression data. By incorporating nonlinear transformations in DGCCA, our model effectively uncovers nonlinear associations across multiple data types. Additionally, the DSR algorithm clusters samples with identical labels, incorporating label information into the nonlinear feature extraction process and thus enhancing the performance of association analysis. The application of the DSR-AdaSMDGCCA algorithm on real data sets identified several AD risk regions (such as the hippocampus, parahippocampus, and fusiform gyrus) and risk genes (including VSIG4, NEDD4L, and PINK1), achieving maximum classification accuracy with the fewest selected features compared to baseline algorithms. Molecular biology enrichment analysis revealed that the pathways enriched by these top genes are intimately linked to AD progression, affirming that our algorithm not only improves correlation analysis performance but also identifies biologically significant markers.

Integrating multi-omics data of childhood asthma using a deep association model.

Childhood asthma is one of the most common respiratory diseases with rising mortality and morbidity. The multi-omics data is providing a new chance to explore collaborative biomarkers and corresponding diagnostic models of childhood asthma. To capture the nonlinear association of multi-omics data and improve interpretability of diagnostic model, we proposed a novel deep association model (DAM) and corresponding efficient analysis framework. First, the Deep Subspace Reconstruction was used to fuse the omics data and diagnostic information, thereby correcting the distribution of the original omics data and reducing the influence of unnecessary data noises. Second, the Joint Deep Semi-Negative Matrix Factorization was applied to identify different latent sample patterns and extract biomarkers from different omics data levels. Third, our newly proposed Deep Orthogonal Canonical Correlation Analysis can rank features in the collaborative module, which are able to construct the diagnostic model considering nonlinear correlation between different omics data levels. Using DAM, we deeply analyzed the transcriptome and methylation data of childhood asthma. The effectiveness of DAM is verified from the perspectives of algorithm performance and biological significance on the independent test dataset, by ablation experiment and comparison with many baseline methods from clinical and biological studies. The DAM-induced diagnostic model can achieve a prediction AUC of 0.912, which is higher than that of many other alternative methods. Meanwhile, relevant pathways and biomarkers of childhood asthma are also recognized to be collectively altered on the gene expression and methylation levels. As an interpretable machine learning approach, DAM simultaneously considers the non-linear associations among samples and those among biological features, which should help explore interpretative biomarker candidates and efficient diagnostic models from multi-omics data analysis for human complex diseases.

Deep Learning Initialized Compressed Sensing (Deli-CS) in Volumetric Spatio-Temporal Subspace Reconstruction.

Introduction: Spatio-temporal MRI methods enable whole-brain multi-parametric mapping at ultra-fast acquisition times through efficient k-space encoding, but can have very long reconstruction times, which limit their integration into clinical practice. Deep learning (DL) is a promising approach to accelerate reconstruction, but can be computationally intensive to train and deploy due to the large dimensionality of spatio-temporal MRI. DL methods also need large training data sets and can produce results that don't match the acquired data if data consistency is not enforced. The aim of this project is to reduce reconstruction time using DL whilst simultaneously limiting the risk of deep learning induced hallucinations, all with modest hardware requirements. Methods: Deep Learning Initialized Compressed Sensing (Deli-CS) is proposed to reduce the reconstruction time of iterative reconstructions by "kick-starting" the iterative reconstruction with a DL generated starting point. The proposed framework is applied to volumetric multi-axis spiral projection MRF that achieves whole-brain T1 and T2 mapping at 1-mm isotropic resolution for a 2-minute acquisition. First, the traditional reconstruction is optimized from over two hours to less than 40 minutes while using more than 90% less RAM and only 4.7 GB GPU memory, by using a memory-efficient GPU implementation. The Deli-CS framework is then implemented and evaluated against the above reconstruction. Results: Deli-CS achieves comparable reconstruction quality with 50% fewer iterations bringing the full reconstruction time to 20 minutes. Conclusion: Deli-CS reduces the reconstruction time of subspace reconstruction of volumetric spatio-temporal acquisitions by providing a warm start to the iterative reconstruction algorithm.

Quantitative longitudinal mapping of radiation-treated prostate cancer using MR fingerprinting with radial acquisition and subspace reconstruction.

MR fingerprinting (MRF) enables fast multiparametric quantitative imaging with a single acquisition and has been shown to improve diagnosis of prostate cancer. However, most prostate MRF studies were performed with spiral acquisitions that are sensitive to B0 inhomogeneities and consequent blurring. In this work, a radial MRF acquisition with a novel subspace reconstruction technique was developed to enable fast T1/T2 mapping in the prostate in under 4 min. The subspace reconstruction exploits the extensive temporal correlations in the MRF dictionary to pre-compute a low dimensional space for the solution and thus reduce the number of radial spokes to accelerate the acquisition. Iterative reconstruction with the subspace model and additional regularization of the signal representation in the subspace is performed to minimize the number of spokes and maintain matching quality and SNR. Reconstruction accuracy was assessed using the ISMRM NIST phantom. In-vivo validation was performed on two healthy subjects and two prostate cancer patients undergoing radiation therapy. The longitudinal repeatability was quantified using the concordance correlation coefficient (CCC) in one of the healthy subjects by repeated scans over 1 year. One prostate cancer patient was scanned at three time points, before initiating therapy and following brachytherapy and external beam radiation. Changes in the T1/T2 maps obtained with the proposed method were quantified. The prostate, peripheral and transitional zones, and visible dominant lesion were delineated for each study, and the statistics and distribution of the quantitative mapping values were analyzed. Significant image quality improvements compared with standard reconstruction methods were obtained with the proposed subspace reconstruction method. A notable decrease in the spread of the T1/T2 values without biasing the estimated mean values was observed with the subspace reconstruction and agreed with reported literature values. The subspace reconstruction enabled visualization of small differences in T1/T2 values in the tumor region within the peripheral zone. Longitudinal imaging of a volunteer subject yielded CCC of 0.89 for MRF T1, and 0.81 for MRF T2 in the prostate gland. Longitudinal imaging of the prostate patient confirmed the feasibility of capturing radiation treatment related changes. This work is a proof-of-concept for a high resolution and fast quantitative mapping using golden-angle radial MRF combined with a subspace reconstruction technique for longitudinal treatment response assessment in subjects undergoing radiation treatment.

Cramér-Rao Bound Optimized Subspace Reconstruction in Quantitative MRI.

We extend the traditional framework for estimating subspace bases that maximize the preserved signal energy to additionally preserve the Cramér-Rao bound (CRB) of the biophysical parameters and, ultimately, improve accuracy and precision in the quantitative maps. To this end, we introduce an approximate compressed CRB based on orthogonalized versions of the signal's derivatives with respect to the model parameters. This approximation permits singular value decomposition (SVD)-based minimization of both the CRB and signal losses during compression. Compared to the traditional SVD approach, the proposed method better preserves the CRB across all biophysical parameters with negligible cost to the preserved signal energy, leading to reduced bias and variance of the parameter estimates in simulation. In vivo, improved accuracy and precision are observed in two quantitative neuroimaging applications, permitting the use of smaller basis sizes in subspace reconstruction and offering significant computational savings.

A Novel Longitudinal Phenotype-Genotype Association Study Based on Deep Feature Extraction and Hypergraph Models for Alzheimer's Disease.

Traditional image genetics primarily uses linear models to investigate the relationship between brain image data and genetic data for Alzheimer's disease (AD) and does not take into account the dynamic changes in brain phenotype and connectivity data across time between different brain areas. In this work, we proposed a novel method that combined Deep Subspace reconstruction with Hypergraph-Based Temporally-constrained Group Sparse Canonical Correlation Analysis (DS-HBTGSCCA) to discover the deep association between longitudinal phenotypes and genotypes. The proposed method made full use of dynamic high-order correlation between brain regions. In this method, the deep subspace reconstruction technique was applied to retrieve the nonlinear properties of the original data, and hypergraphs were used to mine the high-order correlation between two types of rebuilt data. The molecular biological analysis of the experimental findings demonstrated that our algorithm was capable of extracting more valuable time series correlation from the real data obtained by the AD neuroimaging program and finding AD biomarkers across multiple time points. Additionally, we used regression analysis to verify the close relationship between the extracted top brain areas and top genes and found the deep subspace reconstruction approach with a multi-layer neural network was helpful in enhancing clustering performance.

NEAR: An artifact removal pipeline for human newborn EEG data.

Electroencephalography (EEG) is arising as a valuable method to investigate neurocognitive functions shortly after birth. However, obtaining high-quality EEG data from human newborn recordings is challenging. Compared to adults and older infants, datasets are typically much shorter due to newborns' limited attentional span and much noisier due to non-stereotyped artifacts mainly caused by uncontrollable movements. We propose Newborn EEG Artifact Removal (NEAR), a pipeline for EEG artifact removal designed explicitly for human newborns. NEAR is based on two key steps: 1) A novel bad channel detection tool based on the Local Outlier Factor (LOF), a robust outlier detection algorithm; 2) A parameter calibration procedure for adapting to newborn EEG data the algorithm Artifacts Subspace Reconstruction (ASR), developed for artifact removal in mobile adult EEG. Tests on simulated data showed that NEAR outperforms existing methods in removing representative newborn non-stereotypical artifacts. NEAR was validated on two developmental populations (newborns and 9-month-old infants) recorded with two different experimental designs (frequency-tagging and ERP). Results show that NEAR artifact removal successfully reproduces established EEG responses from noisy datasets, with a higher statistical significance than the one obtained by existing artifact removal methods. The EEGLAB-based NEAR pipeline is freely available at https://github.com/vpKumaravel/NEAR.

Automated preprocessing and phase-amplitude coupling analysis of scalp EEG discriminates infantile spasms from controls during wakefulness.

OBJECTIVE: Delta-gamma phase-amplitude coupling in EEG is useful for localizing epileptic sources and to evaluate severity in children with infantile spasms. We (1) develop an automated EEG preprocessing pipeline to clean data using artifact subspace reconstruction (ASR) and independent component (IC) analysis (ICA) and (2) evaluate delta-gamma modulation index (MI) as a method to distinguish children with epileptic spasms (cases) from normal controls during sleep and awake. METHODS: Using 400 scalp EEG datasets (200 sleep, 200 awake) from 100 subjects, we calculated MI after applying high-pass and line-noise filters (Clean 0), and after ASR followed by either conservative (Clean 1) or stringent (Clean 2) artifactual IC rejection. Classification of cases and controls using MI was evaluated with Receiver Operating Characteristics (ROC) to obtain area under curve (AUC). RESULTS: The artifact rejection algorithm reduced raw signal variance by 29-45% and 38-60% for Clean 1 and Clean 2, respectively. MI derived from sleep data, with or without preprocessing, robustly classified the groups (all AUC > 0.98). In contrast, group classification using MI derived from awake data was successful only after Clean 2 (AUC = 0.85). CONCLUSIONS: We have developed an automated EEG preprocessing pipeline to perform artifact rejection and quantify delta-gamma modulation index.

Review of EEG-based pattern classification frameworks for dyslexia.

Dyslexia is a disability that causes difficulties in reading and writing despite average intelligence. This hidden disability often goes undetected since dyslexics are normal and healthy in every other way. Electroencephalography (EEG) is one of the upcoming methods being researched for identifying unique brain activation patterns in dyslexics. The aims of this paper are to examine pros and cons of existing EEG-based pattern classification frameworks for dyslexia and recommend optimisations through the findings to assist future research. A critical analysis of the literature is conducted focusing on each framework's (1) data collection, (2) pre-processing, (3) analysis and (4) classification methods. A wide range of inputs as well as classification approaches has been experimented for the improvement in EEG-based pattern classification frameworks. It was uncovered that incorporating reading- and writing-related tasks to experiments used in data collection may help improve these frameworks instead of using only simple tasks, and those unwanted artefacts caused by body movements in the EEG signals during reading and writing activities could be minimised using artefact subspace reconstruction. Further, support vector machine is identified as a promising classifier to be used in EEG-based pattern classification frameworks for dyslexia.