We present a generative probabilistic approach to discovery of disease subtypes determined by the genetic variants. In many diseases, multiple types of pathology may present simultaneously in a patient, making quantification of the disease challenging. Our method seeks common co-occurring image and genetic patterns in a population as a way to model these two different data types jointly. We assume that each patient is a mixture of multiple disease subtypes and use the joint generative model of image and genetic markers to identify disease subtypes guided by known genetic influences. Our model is based on a variant of the so-called topic models that uncover the latent structure in a collection of data. We derive an efficient variational inference algorithm to extract patterns of co-occurrence and to quantify the presence of heterogeneous disease processes in each patient. We evaluate the method on simulated data and illustrate its use in the context of Chronic Obstructive Pulmonary Disease (COPD) to characterize the relationship between image and genetic signatures of COPD subtypes in a large patient cohort.
BACKGROUND: Computed tomography (CT) helps physicians locate and diagnose pathological conditions. In some conditions, having an airway segmentation method which facilitates reconstruction of the airway from chest CT images can help hugely in the assessment of lung diseases. Many efforts have been made to develop airway segmentation algorithms, but methods are usually not optimized to be reliable across different CT scan parameters. METHODS: In this paper, we present a simple and reliable semi-automatic algorithm which can segment tracheal and bronchial anatomy using the open-source 3D Slicer platform. The method is based on a region growing approach where trachea, right and left bronchi are cropped and segmented independently using three different thresholds. The algorithm and its parameters have been optimized to be efficient across different CT scan acquisition parameters. The performance of the proposed method has been evaluated on EXACT’09 cases and local clinical cases as well as on a breathing pig lung phantom using multiple scans and changing parameters. In particular, to investigate multiple scan parameters reconstruction kernel, radiation dose and slice thickness have been considered. Volume, branch count, branch length and leakage presence have been evaluated. A new method for leakage evaluation has been developed and correlation between segmentation metrics and CT acquisition parameters has been considered. RESULTS: All the considered cases have been segmented successfully with good results in terms of leakage presence. Results on clinical data are comparable to other teams’ methods, as obtained by evaluation against the EXACT09 challenge, whereas results obtained from the phantom prove the reliability of the method across multiple CT platforms and acquisition parameters. As expected, slice thickness is the parameter affecting the results the most, whereas reconstruction kernel and radiation dose seem not to particularly affect airway segmentation. CONCLUSION: The system represents the first open-source airway segmentation platform. The quantitative evaluation approach presented represents the first repeatable system evaluation tool for like-for-like comparison between different airway segmentation platforms. Results suggest that the algorithm can be considered stable across multiple CT platforms and acquisition parameters and can be considered as a starting point for the development of a complete airway segmentation algorithm.
There is a clear need to develop non-invasive markers of substantia nigra progression in Parkinson’s disease. We previously found elevated free-water levels in the substantia nigra for patients with Parkinson’s disease compared with controls in single-site and multi-site cohorts. Here, we test the hypotheses that free-water levels in the substantia nigra of Parkinson’s disease increase following 1 year of progression, and that baseline free-water levels in the substantia nigra predict the change in bradykinesia following 1 year. We conducted a longitudinal study in controls (n = 19) and patients with Parkinson’s disease (n = 25). Diffusion imaging and clinical data were collected at baseline and after 1 year. Free-water analyses were performed on diffusion imaging data using blinded, hand-drawn regions of interest in the posterior substantia nigra. A group effect indicated free-water values were increased in the posterior substantia nigra of patients with Parkinson’s disease compared with controls (P = 0.003) and we observed a significant group × time interaction (P 0.05). Free-water values increased for the Parkinson’s disease group after 1 year (P = 0.006), whereas control free-water values did not change. Baseline free-water values predicted the 1 year change in bradykinesia scores (r = 0.74, P 0.001) and 1 year change in Montreal Cognitive Assessment scores (r = -0.44, P = 0.03). Free-water in the posterior substantia nigra is elevated in Parkinson’s disease, increases with progression of Parkinson’s disease, and predicts subsequent changes in bradykinesia and cognitive status over 1 year. These findings demonstrate that free-water provides a potential non-invasive progression marker of the substantia nigra.
RATIONALE: Chronic obstructive pulmonary disease is associated with a worse overall survival in non-small cell lung cancer. Lung emphysema is one component of chronic obstructive pulmonary disease. We hypothesized that emphysema of the tumor region may result in larger tumors and a poorer overall survival. METHODS: We evaluated 304 cases of non-small cell lung cancer from a prospectively enrolled cohort. The lung was divided into equal volumetric thirds (upper, middle, or lower region). Emphysema was defined as percentage of low-attenuation areas less than -950 Hounsfield units (%LAA-950) and measured for each region. Whole-lung %LAA-950 was defined as the emphysema score of the entire lung parenchyma, whereas regional %LAA-950 was the score within that particular region (upper, middle, or lower). The emphysema score of the region in which the tumor occurred was defined as the tumor %LAA-950. Tumor diameter was measured while blinded to characteristics of the lung parenchyma. A proportional hazards model was used to control for multiple factors associated with survival.
Youth football players may incur hundreds of repetitive head impacts (RHI) in one season. Our recent research suggests that exposure to RHI during a critical neurodevelopmental period prior to age 12 may lead to greater later-life mood, behavioral, and cognitive impairments. Here, we examine the relationship between age of first exposure (AFE) to RHI through tackle football and later-life corpus callosum (CC) microstructure using magnetic resonance diffusion tensor imaging (DTI). Forty retired National Football League (NFL) players, ages 40-65, were matched by age and divided into two groups based on their AFE to tackle football: before age 12 or at age 12 or older. Participants underwent DTI on a 3 Tesla Siemens (TIM-Verio) magnet. The whole CC and five subregions were defined and seeded using deterministic tractography. Dependent measures were fractional anisotropy (FA), trace, axial diffusivity, and radial diffusivity. Results showed that former NFL players in the AFE
RATIONALE: Chronic obstructive pulmonary disease (COPD) is defined by the presence of airflow limitation on spirometry, yet subjects with COPD can have marked differences in computed tomography imaging. These differences may be driven by genetic factors. We hypothesized that a genome-wide association study (GWAS) of quantitative imaging would identify loci not previously identified in analyses of COPD or spirometry. In addition, we sought to determine whether previously described genome-wide significant COPD and spirometric loci were associated with emphysema or airway phenotypes. OBJECTIVES: To identify genetic determinants of quantitative imaging phenotypes. METHODS: We performed a GWAS on two quantitative emphysema and two quantitative airway imaging phenotypes in the COPDGene (non-Hispanic white and African American), ECLIPSE (Evaluation of COPD Longitudinally to Identify Predictive Surrogate Endpoints), NETT (National Emphysema Treatment Trial), and GenKOLS (Genetics of COPD, Norway) studies and on percentage gas trapping in COPDGene. We also examined specific loci reported as genome-wide significant for spirometric phenotypes related to airflow limitation or COPD. MEASUREMENTS AND MAIN RESULTS: The total sample size across all cohorts was 12,031, of whom 9,338 were from COPDGene. We identified five loci associated with emphysema-related phenotypes, one with airway-related phenotypes, and two with gas trapping. These loci included previously reported associations, including the HHIP, 15q25, and AGER loci, as well as novel associations near SERPINA10 and DLC1. All previously reported COPD and a significant number of spirometric GWAS loci were at least nominally (P
RATIONALE: Former smoking history and chronic obstructive pulmonary disease (COPD) are potential risk factors for osteoporosis and fractures. Under existing guidelines for osteoporosis screening, women are included but men are not, and only current smoking is considered. OBJECTIVES: To demonstrate the impact of COPD and smoking history on the risk of osteoporosis and vertebral fracture in men and women. METHODS: Characteristics of participants with low volumetric bone mineral density (vBMD) were identified and related to COPD and other risk factors. We tested associations of sex and COPD with both vBMD and fractures adjusting for age, race, body mass index (BMI), smoking, and glucocorticoid use. MEASUREMENTS AND MAIN RESULTS: vBMD by calibrated quantitative computed tomography (QCT), visually scored vertebral fractures, and severity of lung disease were determined from chest CT scans of 3,321 current and ex-smokers in the COPDGene study. Low vBMD as a surrogate for osteoporosis was calculated from young adult normal values. Male smokers had a small but significantly greater risk of low vBMD (2.5 SD below young adult mean by calibrated QCT) and more fractures than female smokers. Low vBMD was present in 58% of all subjects, was more frequent in those with worse COPD, and rose to 84% among subjects with very severe COPD. Vertebral fractures were present in 37% of all subjects and were associated with lower vBMD at each Global Initiative for Chronic Obstructive Lung Disease stage of severity. Vertebral fractures were most common in the midthoracic region. COPD and especially emphysema were associated with both low vBMD and vertebral fractures after adjustment for steroid use, age, pack-years of smoking, current smoking, and exacerbations. Airway disease was associated with higher bone density after adjustment for other variables. Calibrated QCT identified more subjects with abnormal values than the standard dual-energy X-ray absorptiometry in a subset of subjects and correlated well with prevalent fractures. CONCLUSIONS: Male smokers, with or without COPD, have a significant risk of low vBMD and vertebral fractures. COPD was associated with low vBMD after adjusting for race, sex, BMI, smoking, steroid use, exacerbations, and age. Screening for low vBMD by using QCT in men and women who are smokers will increase opportunities to identify and treat osteoporosis in this at-risk population.
BACKGROUND AND PURPOSE: Right Ventricular to Left Ventricular (RV/LV) diameter ratio has been shown to be a prognostic biomarker for patients suffering from acute Pulmonary Embolism (PE). While Computed Tomography Pulmonary Angiography (CTPA) images used to confirm a clinical suspicion of PE do include information of the heart, a numerical RV/LV diameter ratio is not universally reported, likely because of lack in training, inter-reader variability in the measurements, and additional effort by the radiologist. This study designs and validates a completely automated Computer Aided Detection (CAD) system to compute the axial RV/LV diameter ratio from CTPA images so that the RV/LV diameter ratio can be a more objective metric that is consistently reported in patients for whom CTPA diagnoses PE. MATERIALS AND METHODS: The CAD system was designed specifically for RV/LV measurements. The system was tested in 198 consecutive CTPA patients with acute PE. Its accuracy was evaluated using reference standard RV/LV radiologist measurements and its prognostic value was established for 30-day PE-specific mortality and a composite outcome of 30-day PE-specific mortality or the need for intensive therapies. The study was Institutional Review Board (IRB) approved and HIPAA compliant.
Molecular dynamics (MD) simulations give atomically detailed information on structure and dynamics in amphiphilic bilayer systems on timescales up to about 1 μs. The reorientational dynamics of the C-H bonds is conventionally verified by measurements of (13)C or (2)H nuclear magnetic resonance (NMR) longitudinal relaxation rates R1, which are more sensitive to motional processes with correlation times close to the inverse Larmor frequency, typically around 1-10 ns on standard NMR instrumentation, and are thus less sensitive to the 10-1000 ns timescale motion that can be observed in the MD simulations. We propose an experimental procedure for atomically resolved model-free estimation of the C-H bond effective reorientational correlation time τe, which includes contributions from the entire range of all-atom MD timescales and that can be calculated directly from the MD trajectories. The approach is based on measurements of (13)C R1 and R1ρ relaxation rates, as well as (1)H-(13)C dipolar couplings, and is applicable to anisotropic liquid crystalline lipid or surfactant systems using a conventional solid-state NMR spectrometer and samples with natural isotopic composition. The procedure is demonstrated on a fully hydrated lamellar phase of 1-palmitoyl-2-oleoyl-phosphatidylcholine, yielding values of τe from 0.1 ns for the methyl groups in the choline moiety and at the end of the acyl chains to 3 ns for the g1 methylene group of the glycerol backbone. MD simulations performed with a widely used united-atom force-field reproduce the τe-profile of the major part of the acyl chains but underestimate the dynamics of the glycerol backbone and adjacent molecular segments. The measurement of experimental τe-profiles can be used to study subtle effects on C-H bond reorientational motions in anisotropic liquid crystals, as well as to validate the C-H bond reorientation dynamics predicted in MD simulations of amphiphilic bilayers such as lipid membranes.
Whole-body computed tomography (CT) image registration is important for cancer diagnosis, therapy planning and treatment. Such registration requires accounting for large differences between source and target images caused by deformations of soft organs/tissues and articulated motion of skeletal structures. The registration algorithms relying solely on image processing methods exhibit deficiencies in accounting for such deformations and motion. We propose to predict the deformations and movements of body organs/tissues and skeletal structures for whole-body CT image registration using patient-specific non-linear biomechanical modelling. Unlike the conventional biomechanical modelling, our approach for building the biomechanical models does not require time-consuming segmentation of CT scans to divide the whole body into non-overlapping constituents with different material properties. Instead, a Fuzzy C-Means (FCM) algorithm is used for tissue classification to assign the constitutive properties automatically at integration points of the computation grid. We use only very simple segmentation of the spine when determining vertebrae displacements to define loading for biomechanical models. We demonstrate the feasibility and accuracy of our approach on CT images of seven patients suffering from cancer and aortic disease. The results confirm that accurate whole-body CT image registration can be achieved using a patient-specific non-linear biomechanical model constructed without time-consuming segmentation of the whole-body images.
