Publications

Patients with fibrous dysplasia (FD) often present with craniofacial lesions that affect the trigeminal nerve system. Debilitating pain, headache, and migraine are frequently experienced by FD patients with poor prognosis, while some individuals with similar bone lesions are asymptomatic. The clinical and biological factors that contribute to the etiopathogenesis of pain in craniofacial FD are largely unknown. We present two adult females with comparable craniofacial FD lesion size and location, as measured by 18F-sodium fluoride positron emission tomography/computed tomography (PET/CT), yet their respective pain phenotypes differed significantly. Over 4 weeks, the average pain reported by Patient A was 0.4/0-10 scale. Patient B reported average pain of 7.8/0-10 scale distributed across the entire skull and left facial region. Patient B did not experience pain relief from analgesics or more aggressive treatments (denosumab). In both patients, evaluation of trigeminal nerve divisions (V1, V2, and V3) with CT and magnetic resonance imaging (MRI) revealed nerve compression and displacement with more involvement of the left trigeminal branches relative to the right. First-time employment of diffusion MRI and tractography suggested reduced apparent fiber density within the cisternal segment of the trigeminal nerve, particularly for Patient B and in the left hemisphere. These cases highlight heterogeneous clinical presentation and neurobiological properties in craniofacial FD and also, the disconnect between peripheral pathology and pain severity. We hypothesize that a detailed phenotypic characterization of patients that incorporates an advanced imaging approach probing the trigeminal system may provide enhanced insights into the variable experiences with pain in craniofacial FD.

Almost 25% of all older adults experience difficulty walking. Mobility difficulties for older adults are more pronounced when they perform a simultaneous cognitive task while walking (i.e., dual task walking). Although it is known that aging results in widespread brain atrophy, few studies have integrated across more than one neuroimaging modality to comprehensively examine the structural neural correlates that may underlie dual task walking in older age. We collected spatiotemporal gait data during single and dual task walking for 37 young (18-34 years) and 23 older adults (66-86 years). We also collected T 1-weighted and diffusion-weighted MRI scans to determine how brain structure differs in older age and relates to dual task walking. We addressed two aims: (1) to characterize age differences in brain structure across a range of metrics including volumetric, surface, and white matter microstructure; and (2) to test for age group differences in the relationship between brain structure and the dual task cost (DTcost) of gait speed and variability. Key findings included widespread brain atrophy for the older adults, with the most pronounced age differences in brain regions related to sensorimotor processing. We also found multiple associations between regional brain atrophy and greater DTcost of gait speed and variability for the older adults. The older adults showed a relationship of both thinner temporal cortex and shallower sulcal depth in the frontal, sensorimotor, and parietal cortices with greater DTcost of gait. Additionally, the older adults showed a relationship of ventricular volume and superior longitudinal fasciculus free-water corrected axial and radial diffusivity with greater DTcost of gait. These relationships were not present for the young adults. Stepwise multiple regression found sulcal depth in the left precentral gyrus, axial diffusivity in the superior longitudinal fasciculus, and sex to best predict DTcost of gait speed, and cortical thickness in the superior temporal gyrus to best predict DTcost of gait variability for older adults. These results contribute to scientific understanding of how individual variations in brain structure are associated with mobility function in aging. This has implications for uncovering mechanisms of brain aging and for identifying target regions for mobility interventions for aging populations.

We present a novel weakly-supervised framework for classifying whole slide images (WSIs). WSIs, due to their gigapixel resolution, are commonly processed by patch-wise classification with patch-level labels. However, patch-level labels require precise annotations, which is expensive and usually unavailable on clinical data. With image-level labels only, patch-wise classification would be sub-optimal due to inconsistency between the patch appearance and image-level label. To address this issue, we posit that WSI analysis can be effectively conducted by integrating information at both high magnification (local) and low magnification (regional) levels. We auto-encode the visual signals in each patch into a latent embedding vector representing local information, and down-sample the raw WSI to hardware-acceptable thumbnails representing regional information. The WSI label is then predicted with a Dual-Stream Network (DSNet), which takes the transformed local patch embeddings and multi-scale thumbnail images as inputs and can be trained by the image-level label only. Experiments conducted on three large-scale public datasets demonstrate that our method outperforms all recent state-of-the-art weakly-supervised WSI classification methods.

BACKGROUND: Chronic obstructive pulmonary disease (COPD) and idiopathic pulmonary fibrosis (IPF) are characterized by shared exposures and clinical features, but distinct genetic and pathologic features exist. These features have not been well-studied using large-scale gene expression datasets. We hypothesized that there are divergent gene, pathway, and cellular signatures between COPD and IPF. METHODS: We performed RNA-sequencing on lung tissues from individuals with IPF (n = 231) and COPD (n = 377) compared to control (n = 267), defined as individuals with normal spirometry. We grouped the overlapping differential expression gene sets based on direction of expression and examined the resultant sets for genes of interest, pathway enrichment, and cell composition. Using gene set variation analysis, we validated the overlap group gene sets in independent COPD and IPF data sets. RESULTS: We found 5010 genes differentially expressed between COPD and control, and 11,454 genes differentially expressed between IPF and control (1% false discovery rate). 3846 genes overlapped between IPF and COPD. Several pathways were enriched for genes upregulated in COPD and downregulated in IPF; however, no pathways were enriched for genes downregulated in COPD and upregulated in IPF. There were many myeloid cell genes with increased expression in COPD but decreased in IPF. We found that the genes upregulated in COPD but downregulated in IPF were associated with lower lung function in the independent validation cohorts. CONCLUSIONS: We identified a divergent gene expression signature between COPD and IPF, with increased expression in COPD and decreased in IPF. This signature is associated with worse lung function in both COPD and IPF.

BACKGROUND: Despite improvement in treatment strategies of atrial fibrillation (AF), a considerable number of patients still experience recurrence of atrial tachyarrhythmia (ATA) following catheter ablation (CA). This study aimed to investigate the prognostic value of left atrial (LA) deformation analysis in a large group of patients undergoing CA for AF. METHODS: This study included 678 patients with AF. Echocardiography including two-dimensional speckle tracking echocardiography (2DSTE) was performed in all patients prior to CA. Logistic regression analysis was used to assess the association between ATA recurrence and LA strain during reservoir phase (LASr), LA strain during contraction phase (LASct), and LA strain during conduit phase (LAScd). RESULTS: During one-year follow-up, 274 (40%) experienced ATA recurrence. Median age of the included study population was 63.2 years (IQR: 55.5, 69.5) and 485 (72%) were male. Patients with recurrence had lower LASr (22.6% vs. 25.1%, p = 0.001) and LASct (10.7% vs. 12.4%, p < 0.001). No difference in LAScd was observed. After adjusting for potential clinical and echocardiographic confounders LASr (OR = 1.04, CI95% [1.01; 1.07], p = 0.015, per 1% decrease) and LASct (OR = 1.06, CI95% [1.02; 1.11], p = 0.007, per 1% decrease) remained independent predictors of recurrence. However, in patients with a normal-sized LA (LA volume index<34 mL/m2), only LASct remained an independent predictor of recurrence (OR = 1.07, CI95% [1.01; 1.12], p = 0.012, per 1% decrease). CONCLUSION: In patients undergoing CA for AF, LA deformation analysis by 2DSTE could be of use in risk stratification in clinical practice regarding ATA recurrence, even in patients with a normal-sized LA.

Deep brain stimulation (DBS) is a well-established neurosurgical procedure for movement disorders that is also being explored for treatment-resistant psychiatric conditions. This review highlights important consideration for DBS simulation and data analysis. The literature on DBS has expanded considerably in recent years, and this article aims to identify important trends in the field. During DBS planning, surgery, and follow up sessions, several large data sets are created for each patient, and it becomes clear that any group analysis of such data is a big data analysis problem and has to be handled with care. The aim of this review is to provide an update and overview from a neuroengineering perspective of the current DBS techniques, technical aids, and emerging tools with the focus on patient-specific electric field (EF) simulations, group analysis, and visualization in the DBS domain. Examples are given from the state-of-the-art literature including our own research. This work reviews different analysis methods for EF simulations, tractography, deep brain anatomical templates, and group analysis. Our analysis highlights that group analysis in DBS is a complex multi-level problem and selected parameters will highly influence the result. DBS analysis can only provide clinically relevant information if the EF simulations, tractography results, and derived brain atlases are based on as much patient-specific data as possible. A trend in DBS research is creation of more advanced and intuitive visualization of the complex analysis results suitable for the clinical environment.

In cerebral small vessel disease (CSVD), both white matter hyperintensities (WMH) of presumed vascular origin and the normal-appearing white matter (NAWM) contain microstructural brain alterations on diffusion-weighted MRI (DWI). Contamination of DWI-derived metrics by extracellular free-water can be corrected with free-water (FW) imaging. We investigated the alterations in FW and FW-corrected fractional anisotropy (FA-t) in WMH and surrounding tissue and their association with cerebrovascular risk factors. We analysed 1,000 MRI datasets from the Hamburg City Health Study. DWI was used to generate FW and FA-t maps. WMH masks were segmented on FLAIR and T1-weighted MRI and dilated repeatedly to create 8 NAWM masks representing increasing distance from WMH. Linear models were applied to compare FW and FA-t across WMH and NAWM masks and in association with cerebrovascular risk. Median age was 64 ± 14 years. FW and FA-t were altered 8 mm and 12 mm beyond WMH, respectively. Smoking was significantly associated with FW in NAWM (p = 0.008) and FA-t in WMH (p = 0.008) and in NAWM (p = 0.003) while diabetes and hypertension were not. Further research is necessary to examine whether FW and FA-t alterations in NAWM are predictors for developing WMH.

There are many ways to acquire and process diffusion MRI (dMRI) data for group studies, but it is unknown which maximizes the sensitivity to white matter (WM) pathology. Inspired by this question, we analyzed data acquired for diffusion tensor imaging (DTI) and diffusion kurtosis imaging (DKI) at 3T (3T-DTI and 3T-DKI) and DTI at 7T in patients with systemic lupus erythematosus (SLE) and healthy controls (HC). Parameter estimates in 72 WM tracts were obtained using TractSeg. The impact on the sensitivity to WM pathology was evaluated for the diffusion protocol, the magnetic field strength, and the processing pipeline. Sensitivity was quantified in terms of Cohen’s d for group comparison. Results showed that the choice of diffusion protocol had the largest impact on the effect size. The effect size in fractional anisotropy (FA) across all WM tracts was 0.26 higher when derived by DTI than by DKI and 0.20 higher in 3T compared with 7T. The difference due to the diffusion protocol was larger than the difference due to magnetic field strength for the majority of diffusion parameters. In contrast, the difference between including or excluding different processing steps was near negligible, except for the correction of distortions from eddy currents and motion which had a clearly positive impact. For example, effect sizes increased on average by 0.07 by including motion and eddy correction for FA derived from 3T-DTI. Effect sizes were slightly reduced by the incorporation of denoising and Gibbs-ringing removal (on average by 0.011 and 0.005, respectively). Smoothing prior to diffusion model fitting generally reduced effect sizes. In summary, 3T-DTI in combination with eddy current and motion correction yielded the highest sensitivity to WM pathology in patients with SLE. However, our results also indicated that the 3T-DKI and 7T-DTI protocols used here may be adjusted to increase effect sizes.

Background: Tumor-related hyperintensities in high b-value diffusion-weighted imaging (DWI) are radiologically important in the workup of gliomas. However, the white matter may also appear as hyperintense, which may conflate interpretation. Purpose: To investigate whether DWI with spherical b-tensor encoding (STE) can be used to suppress white matter and enhance the conspicuity of glioma hyperintensities unrelated to white matter. Materials and Methods: Twenty-five patients with a glioma tumor and at least one pathology-related hyperintensity on DWI underwent conventional MRI at 3 T. The DWI was performed both with linear and spherical tensor encoding (LTE-DWI and STE-DWI). The LTE-DWI here refers to the DWI obtained with conventional diffusion encoding and averaged across diffusion-encoding directions. Retrospectively, the differences in contrast between LTE-DWI and STE-DWI, obtained at a b-value of 2,000 s/mm2, were evaluated by comparing hyperintensities and contralateral normal-appearing white matter (NAWM) both visually and quantitatively in terms of the signal intensity ratio (SIR) and contrast-to-noise ratio efficiency (CNReff). Results: The spherical tensor encoding DWI was more effective than LTE-DWI at suppressing signals from white matter and improved conspicuity of pathology-related hyperintensities. The median SIR improved in all cases and on average by 28%. The median (interquartile range) SIR was 1.9 (1.6 - 2.1) for STE and 1.4 (1.3 - 1.7) for LTE, with a significant difference of 0.4 (0.3 -0.5) (p < 10-4, paired U-test). In 40% of the patients, the SIR was above 2 for STE-DWI, but with LTE-DWI, the SIR was below 2 for all patients. The CNReff of STE-DWI was significantly higher than of LTE-DWI: 2.5 (2 - 3.5) vs. 2.3 (1.7 - 3.1), with a significant difference of 0.4 (-0.1 -0.6) (p < 10-3, paired U-test). The STE improved CNReff in 70% of the cases. We illustrate the benefits of STE-DWI in three patients, where STE-DWI may facilitate an improved radiological description of tumor-related hyperintensity, including one case that could have been missed out if only LTE-DWI was inspected. Conclusion: The contrast mechanism of high b-value STE-DWI results in a stronger suppression of white matter than conventional LTE-DWI, and may, therefore, be more sensitive and specific for assessment of glioma tumors and DWI-hyperintensities.

Extracellular free water (FW) increases are suggested to better provide pathophysiological information in brain aging than conventional biomarkers such as fractional anisotropy. The aim of the present study was to determine the relationship between conventional biomarkers, FW in white matter hyperintensities (WMH), FW in normal appearing white matter (NAWM) and in white matter tracts and executive functions (EF) with a speed component in elderly persons. We examined 226 healthy elderly participants (median age 69.83 years, IQR: 56.99-74.42) who underwent brain MRI and neuropsychological examination. FW in WMH and in NAWM as well as FW corrected diffusion metrics and measures derived from conventional MRI (white matter hyperintensities, brain volume, lacunes) were used in partial correlation (adjusted for age) to assess their correlation with EF with a speed component. Random forest analysis was used to assess the relative importance of these variables as determinants. Lastly, linear regression analyses of FW in white matter tracts corrected for risk factors of cognitive and white matter deterioration, were used to examine the role of specific tracts on EF with a speed component, which were then ranked with random forest regression. Partial correlation analyses revealed that almost all imaging metrics showed a significant association with EF with a speed component (r=-0.213 - 0.266). Random forest regression highlighted FW in WMH and in NAWM as most important among all diffusion and structural MRI metrics. The fornix (R2=0.421, p=0.018) and the corpus callosum (genu (R2=0.418, p=0.021), prefrontal (R2=0.416, p=0.026), premotor (R2=0.418, p=0.021)) were associated with EF with a speed component in tract based regression analyses and had highest variables importance. In a normal aging population FW in WMH and NAWM is more closely related to EF with a speed component than standard DTI and brain structural measures. Higher amounts of FW in the fornix and the frontal part of the corpus callosum leads to deteriorating EF with a speed component.