Search

Abstract Patients with congenital heart disease often have cardiac anatomy that deviates significantly from normal, frequently requiring multiple heart surgeries. Image segmentation from a preoperative cardiovascular magnetic resonance (CMR) scan would enable creation of patient-specific 3D surface models of the heart, which have potential to improve surgical planning, enable surgical simulation, and allow automatic computation of quantitative metrics of heart function. However, there is no publicly available CMR dataset for whole-heart segmentation in patients with congenital heart disease. Here, we release the HVSMR-2.0 dataset, comprising 60 CMR scans alongside manual segmentation masks of the 4 cardiac chambers and 4 great vessels. The images showcase a wide range of heart defects and prior surgical interventions. The dataset also includes masks of required and optional extents of the great vessels, enabling fairer comparisons across algorithms. Detailed diagnoses for each subject are also provided. By releasing HVSMR-2.0, we aim to encourage development of robust segmentation algorithms and clinically relevant tools for congenital heart disease.

IntroductionAs the life expectancy of children with congenital heart disease (CHD) is rapidly increasing and the adult population with CHD is growing, there is an unmet need to improve clinical workflow and efficiency of analysis. Cardiovascular magnetic resonance (CMR) is a noninvasive imaging modality for monitoring patients with CHD. CMR exam is based on multiple breath-hold 2-dimensional (2D) cine acquisitions that should be precisely prescribed and is expert and institution dependent. Moreover, 2D cine images have relatively thick slices, which does not allow for isotropic delineation of ventricular structures. Thus, development of an isotropic 3D cine acquisition and automatic segmentation method is worthwhile to make CMR workflow straightforward and efficient, as the present work aims to establish.MethodsNinety-nine patients with many types of CHD were imaged using a non-angulated 3D cine CMR sequence covering the whole-heart and great vessels. Automatic supervised and semi-supervised deep-learning-based methods were developed for whole-heart segmentation of 3D cine images to separately delineate the cardiac structures, including both atria, both ventricles, aorta, pulmonary arteries, and superior and inferior vena cavae. The segmentation results derived from the two methods were compared with the manual segmentation in terms of Dice score, a degree of overlap agreement, and atrial and ventricular volume measurements.ResultsThe semi-supervised method resulted in a better overlap agreement with the manual segmentation than the supervised method for all 8 structures (Dice score 83.23 ± 16.76% vs. 77.98 ± 19.64%; P-value ≤0.001). The mean difference error in atrial and ventricular volumetric measurements between manual segmentation and semi-supervised method was lower (bias ≤ 5.2 ml) than the supervised method (bias ≤ 10.1 ml).DiscussionThe proposed semi-supervised method is capable of cardiac segmentation and chamber volume quantification in a CHD population with wide anatomical variability. It accurately delineates the heart chambers and great vessels and can be used to accurately calculate ventricular and atrial volumes throughout the cardiac cycle. Such a segmentation method can reduce inter- and intra- observer variability and make CMR exams more standardized and efficient.

The human thalamus is a highly connected brain structure, which is key for the control of numerous functions and is involved in several neurological disorders. Recently, neuroimaging studies have increasingly focused on the volume and connectivity of the specific nuclei comprising this structure, rather than looking at the thalamus as a whole. However, accurate identification of cytoarchitectonically designed histological nuclei on standard in vivo structural MRI is hampered by the lack of image contrast that can be used to distinguish nuclei from each other and from surrounding white matter tracts. While diffusion MRI may offer such contrast, it has lower resolution and lacks some boundaries visible in structural imaging. In this work, we present a Bayesian segmentation algorithm for the thalamus. This algorithm combines prior information from a probabilistic atlas with likelihood models for both structural and diffusion MRI, allowing segmentation of 25 thalamic labels per hemisphere informed by both modalities. We present an improved probabilistic atlas, incorporating thalamic nuclei identified from histology and 45 white matter tracts surrounding the thalamus identified in ultra-high gradient strength diffusion imaging. We present a family of likelihood models for diffusion tensor imaging, ensuring compatibility with the vast majority of neuroimaging datasets that include diffusion MRI data. The use of these diffusion likelihood models greatly improves identification of nuclear groups versus segmentation based solely on structural MRI. Dice comparison of 5 manually identifiable groups of nuclei to ground truth segmentations show improvements of up to 10 percentage points. Additionally, our chosen model shows a high degree of reliability, with median test-retest Dice scores above 0.85 for four out of five nuclei groups, whilst also offering improved detection of differential thalamic involvement in Alzheimer’s disease (AUROC 81.98%). The probabilistic atlas and segmentation tool will be made publicly available as part of the neuroimaging package FreeSurfer (https://freesurfer.net/fswiki/ThalamicNucleiDTI).

Hypoplastic left heart syndrome (HLHS) is a severe congenital heart defect in which the right ventricle and associated tricuspid valve (TV) alone support the circulation. TV failure is thus associated with heart failure, and the outcome of TV valve repair are currently poor. 3D echocardiography (3DE) can generate high-quality images of the valve, but segmentation is necessary for precise modeling and quantification. There is currently no robust methodology for rapid TV segmentation, limiting the clinical application of these technologies to this challenging population. We utilized a Fully Convolutional Network (FCN) to segment tricuspid valves from transthoracic 3DE. We trained on 133 3DE image-segmentation pairs and validated on 28 images. We then assessed the effect of varying inputs to the FCN using Mean Boundary Distance (MBD) and Dice Similarity Coefficient (DSC). The FCN with the input of an annular curve achieved a median DSC of 0.86 [IQR: 0.81–0.88] and MBD of 0.35 [0.23–0.4] mm for the merged segmentation and an average DSC of 0.77 [0.73–0.81] and MBD of 0.6 [0.44–0.74] mm for individual TV leaflet segmentation. The addition of commissural landmarks improved individual leaflet segmentation accuracy to an MBD of 0.38 [0.3–0.46] mm. FCN-based segmentation of the tricuspid valve from transthoracic 3DE is feasible and accurate. The addition of an annular curve and commissural landmarks improved the quality of the segmentations with MBD and DSC within the range of human inter-user variability. Fast and accurate FCN-based segmentation of the tricuspid valve in HLHS may enable rapid modeling and quantification, which in the future may inform surgical planning. We are now working to deploy this network for public use.

Objective: To personalize the prognostication of post-stroke outcome using MRI-detected cerebrovascular pathology, we sought to investigate the association between the excessive white matter hyperintensity (WMH) burden unaccounted for by the traditional stroke risk profile of individual patients and their long-term functional outcomes after a stroke.Methods: We included 890 patients who survived after an acute ischemic stroke from the MRI-Genetics Interface Exploration (MRI-GENIE) study, for whom data on vascular risk factors (VRFs), including age, sex, atrial fibrillation, diabetes mellitus, hypertension, coronary artery disease, smoking, prior stroke history, as well as acute stroke severity, 3- to−6-month modified Rankin Scale score (mRS), WMH, and brain volumes, were available. We defined the unaccounted WMH (uWMH) burden via modeling of expected WMH burden based on the VRF profile of each individual patient. The association of uWMH and mRS score was analyzed by linear regression analysis. The odds ratios of patients who achieved full functional independence (mRS < 2) in between trichotomized uWMH burden groups were calculated by pair-wise comparisons.Results: The expected WMH volume was estimated with respect to known VRFs. The uWMH burden was associated with a long-term functional outcome (β = 0.104, p < 0.01). Excessive uWMH burden significantly reduced the odds of achieving full functional independence after a stroke compared to the low and average uWMH burden [OR = 0.4, 95% CI: (0.25, 0.63), p < 0.01 and OR = 0.61, 95% CI: (0.42, 0.87), p < 0.01, respectively].Conclusion: The excessive amount of uWMH burden unaccounted for by the traditional VRF profile was associated with worse post-stroke functional outcomes. Further studies are needed to evaluate a lifetime brain injury reflected in WMH unrelated to the VRF profile of a patient as an important factor for stroke recovery and a plausible indicator of brain health.

Most existing algorithms for automatic 3D morphometry of human brain MRI scans are designed for data with near-isotropic voxels at approximately 1 mm resolution, and frequently have contrast constraints as well-typically requiring T1-weighted images (e.g., MP-RAGE scans). This limitation prevents the analysis of millions of MRI scans acquired with large inter-slice spacing in clinical settings every year. In turn, the inability to quantitatively analyze these scans hinders the adoption of quantitative neuro imaging in healthcare, and also precludes research studies that could attain huge sample sizes and hence greatly improve our understanding of the human brain. Recent advances in convolutional neural networks (CNNs) are producing outstanding results in super-resolution and contrast synthesis of MRI. However, these approaches are very sensitive to the specific combination of contrast, resolution and orientation of the input images, and thus do not generalize to diverse clinical acquisition protocols – even within sites. In this article, we present SynthSR, a method to train a CNN that receives one or more scans with spaced slices, acquired with different contrast, resolution and orientation, and produces an isotropic scan of canonical contrast (typically a 1 mm MP-RAGE). The presented method does not require any preprocessing, beyond rigid coregistration of the input scans. Crucially, SynthSR trains on synthetic input images generated from 3D segmentations, and can thus be used to train CNNs for any combination of contrasts, resolutions and orientations without high-resolution real images of the input contrasts. We test the images generated with SynthSR in an array of common downstream analyses, and show that they can be reliably used for subcortical segmentation and volumetry, image registration (e.g., for tensor-based morphometry), and, if some image quality requirements are met, even cortical thickness morphometry. The source code is publicly available at https://github.com/BBillot/SynthSR.

ObjectiveNeuroimaging measurements of brain structural integrity are thought to be surrogates for brain health, but precise assessments require dedicated advanced image acquisitions. By means of quantitatively describing conventional images, radiomic analyses hold potential for evaluating brain health. We sought to: (1) evaluate radiomics to assess brain structural integrity by predicting white matter hyperintensities burdens (WMH) and (2) uncover associations between predictive radiomic features and clinical phenotypes.MethodsWe analyzed a multi-site cohort of 4,163 acute ischemic strokes (AIS) patients with T2-FLAIR MR images with total brain and WMH segmentations. Radiomic features were extracted from normal-appearing brain tissue (brain mask–WMH mask). Radiomics-based prediction of personalized WMH burden was done using ElasticNet linear regression. We built a radiomic signature of WMH with stable selected features predictive of WMH burden and then related this signature to clinical variables using canonical correlation analysis (CCA).ResultsRadiomic features were predictive of WMH burden (R2 = 0.855 ± 0.011). Seven pairs of canonical variates (CV) significantly correlated the radiomics signature of WMH and clinical traits with respective canonical correlations of 0.81, 0.65, 0.42, 0.24, 0.20, 0.15, and 0.15 (FDR-corrected p-valuesCV1–6 < 0.001, p-valueCV7 = 0.012). The clinical CV1 was mainly influenced by age, CV2 by sex, CV3 by history of smoking and diabetes, CV4 by hypertension, CV5 by atrial fibrillation (AF) and diabetes, CV6 by coronary artery disease (CAD), and CV7 by CAD and diabetes.ConclusionRadiomics extracted from T2-FLAIR images of AIS patients capture microstructural damage of the cerebral parenchyma and correlate with clinical phenotypes, suggesting different radiographical textural abnormalities per cardiovascular risk profile. Further research could evaluate radiomics to predict the progression of WMH and for the follow-up of stroke patients’ brain health.

Background: Magnetic resonance imaging (MRI) serves as a cornerstone in defining stroke phenotype and etiological subtype through examination of ischemic stroke lesion appearance and is therefore an essential tool in linking genetic traits and stroke. Building on baseline MRI examinations from the centralized and structured radiological assessments of ischemic stroke patients in the Stroke Genetics Network, the results of the MRI-Genetics Interface Exploration (MRI-GENIE) study are described in this work.Methods: The MRI-GENIE study included patients with symptoms caused by ischemic stroke (N = 3,301) from 12 international centers. We established and used a structured reporting protocol for all assessments. Two neuroradiologists, using a blinded evaluation protocol, independently reviewed the baseline diffusion-weighted images (DWIs) and magnetic resonance angiography images to determine acute lesion and vascular occlusion characteristics.Results: In this systematic multicenter radiological analysis of clinical MRI from 3,301 acute ischemic stroke patients according to a structured prespecified protocol, we identified that anterior circulation infarcts were most prevalent (67.4%), that infarcts in the middle cerebral artery (MCA) territory were the most common, and that the majority of large artery occlusions 0 to 48 h from ictus were in the MCA territory. Multiple acute lesions in one or several vascular territories were common (11%). Of 2,238 patients with unilateral DWI lesions, 52.6% had left-sided infarct lateralization (P = 0.013 for χ2 test).Conclusions: This large-scale analysis of a multicenter MRI-based cohort of AIS patients presents a unique imaging framework facilitating the relationship between imaging and genetics for advancing the knowledge of genetic traits linked to ischemic stroke.

Abstract Fetal health is critically dependent on placental function, especially placental transport of oxygen from mother to fetus. When fetal growth is compromised, placental insufficiency must be distinguished from modest genetic growth potential. If placental insufficiency is present, the physician must trade off the risk of prolonged fetal exposure to placental insufficiency against the risks of preterm delivery. Current ultrasound methods to evaluate the placenta are indirect and insensitive. We propose to use Blood-Oxygenation-Level-Dependent (BOLD) MRI with maternal hyperoxia to quantitatively assess mismatch in placental function in seven monozygotic twin pairs naturally matched for genetic growth potential. In-utero BOLD MRI time series were acquired at 29 to 34 weeks gestational age. Maps of oxygen Time-To-Plateau (TTP) were obtained in the placentas by voxel-wise fitting of the time series. Fetal brain and liver volumes were measured based on structural MR images. After delivery, birth weights were obtained and placental pathological evaluations were performed. Mean placental TTP negatively correlated with fetal liver and brain volumes at the time of MRI as well as with birth weights. Mean placental TTP positively correlated with placental pathology. This study demonstrates the potential of BOLD MRI with maternal hyperoxia to quantify regional placental function in vivo.

White matter hyperintensity (WMH) burden is a critically important cerebrovascular phenotype linked to prediction of diagnosis and prognosis of diseases, such as acute ischemic stroke (AIS). However, current approaches to its quantification on clinical MRI often rely on time intensive manual delineation of the disease on T2 fluid attenuated inverse recovery (FLAIR), which hinders high-throughput analyses such as genetic discovery.In this work, we present a fully automated pipeline for quantification of WMH in clinical large-scale studies of AIS. The pipeline incorporates automated brain extraction, intensity normalization and WMH segmentation using spatial priors. We first propose a brain extraction algorithm based on a fully convolutional deep learning architecture, specifically designed for clinical FLAIR images. We demonstrate that our method for brain extraction outperforms two commonly used and publicly available methods on clinical quality images in a set of 144 subject scans across 12 acquisition centers, based on dice coefficient (median 0.95; inter-quartile range 0.94–0.95; p