Neuroimaging meta regression for coordinate based meta analysis data with a spatial model.

Coordinate-based meta-analysis combines evidence from a collection of neuroimaging studies to estimate brain activation. In such analyses, a key practical challenge is to find a computationally efficient approach with good statistical interpretability to model the locations of activation foci. In this article, we propose a generative coordinate-based meta-regression (CBMR) framework to approximate a smooth activation intensity function and investigate the effect of study-level covariates (e.g. year of publication, sample size). We employ a spline parameterization to model the spatial structure of brain activation and consider four stochastic models for modeling the random variation in foci. To examine the validity of CBMR, we estimate brain activation on 20 meta-analytic datasets, conduct spatial homogeneity tests at the voxel level, and compare the results to those generated by existing kernel-based and model-based approaches. Keywords: generalized linear models; meta-analysis; spatial statistics; statistical modeling.

Partitioning variance in cortical morphometry into genetic, environmental, and subject-specific components.

The relative contributions of genetic variation and experience in shaping the morphology of the adolescent brain are not fully understood. Using longitudinal data from 11,665 subjects in the ABCD Study, we fit vertex-wise variance components including family effects, genetic effects, and subject-level effects using a computationally efficient framework. Variance in cortical thickness and surface area is largely attributable to genetic influence, whereas sulcal depth is primarily explained by subject-level effects. Our results identify areas with heterogeneous distributions of heritability estimates that have not been seen in previous work using data from cortical regions. We discuss the biological importance of subject-specific variance and its implications for environmental influences on cortical development and maturation.

Leveraging the adolescent brain cognitive development study to improve behavioral prediction from neuroimaging in smaller replication samples.

Neuroimaging is a popular method to map brain structural and functional patterns to complex human traits. Recently published observations cast doubt upon these prospects, particularly for prediction of cognitive traits from structural and resting state functional magnetic resonance imaging (MRI). We leverage baseline data from thousands of children in the Adolescent Brain Cognitive DevelopmentSM Study to inform the replication sample size required with univariate and multivariate methods across different imaging modalities to detect reproducible brain-behavior associations. We demonstrate that by applying multivariate methods to high-dimensional brain imaging data, we can capture lower dimensional patterns of structural and functional brain architecture that correlate robustly with cognitive phenotypes and are reproducible with only 41 individuals in the replication sample for working memory-related functional MRI, and ~ 100 subjects for structural and resting state MRI. Even with 100 random re-samplings of 100 subjects in discovery, prediction can be adequately powered with 66 subjects in replication for multivariate prediction of cognition with working memory task functional MRI. These results point to an important role for neuroimaging in translational neurodevelopmental research and showcase how findings in large samples can inform reproducible brain-behavior associations in small sample sizes that are at the heart of many research programs and grants.

FEMA: Fast and efficient mixed-effects algorithm for large sample whole-brain imaging data.

The linear mixed-effects model (LME) is a versatile approach to account for dependence among observations. Many large-scale neuroimaging datasets with complex designs have increased the need for LME; however LME has seldom been used in whole-brain imaging analyses due to its heavy computational requirements. In this paper, we introduce a fast and efficient mixed-effects algorithm (FEMA) that makes whole-brain vertex-wise, voxel-wise, and connectome-wide LME analyses in large samples possible. We validate FEMA with extensive simulations, showing that the estimates of the fixed effects are equivalent to standard maximum likelihood estimates but obtained with orders of magnitude improvement in computational speed. We demonstrate the applicability of FEMA by studying the cross-sectional and longitudinal effects of age on region-of-interest level and vertex-wise cortical thickness, as well as connectome-wide functional connectivity values derived from resting state functional MRI, using longitudinal imaging data from the Adolescent Brain Cognitive DevelopmentSM Study release 4.0. Our analyses reveal distinct spatial patterns for the annualized changes in vertex-wise cortical thickness and connectome-wide connectivity values in early adolescence, highlighting a critical time of brain maturation. The simulations and application to real data show that FEMA enables advanced investigation of the relationships between large numbers of neuroimaging metrics and variables of interest while considering complex study designs, including repeated measures and family structures, in a fast and efficient manner. The source code for FEMA is available via: https://github.com/cmig-research-group/cmig_tools/.

NHS Health Check attendance is associated with reduced multiorgan disease risk: a matched cohort study in the UK Biobank.

BACKGROUND: The NHS Health Check is a preventive programme in the UK designed to screen for cardiovascular risk and to aid in primary disease prevention. Despite its widespread implementation, the effectiveness of the NHS Health Check for longer-term disease prevention is unclear. In this study, we measured the rate of new diagnoses in UK Biobank participants who underwent the NHS Health Check compared with those who did not. METHODS: Within the UK Biobank prospective study, 48,602 NHS Health Check recipients were identified from linked primary care records. These participants were then covariate-matched on an extensive range of socio-demographic, lifestyle, and medical factors with 48,602 participants without record of the check. Follow-up diagnoses were ascertained from health records over an average of 9 years (SD 2 years) including hypertension, diabetes, hypercholesterolaemia, stroke, dementia, myocardial infarction, atrial fibrillation, heart failure, fatty liver disease, alcoholic liver disease, liver cirrhosis, liver failure, acute kidney injury, chronic kidney disease (stage 3 +), cardiovascular mortality, and all-cause mortality. Time-varying survival modelling was used to compare adjusted outcome rates between the groups. RESULTS: In the immediate 2 years after the NHS Health Check, higher diagnosis rates were observed for hypertension, high cholesterol, and chronic kidney disease among health check recipients compared to their matched counterparts. However, in the longer term, NHS Health Check recipients had significantly lower risk across all multiorgan disease outcomes and reduced rates of cardiovascular and all-cause mortality. CONCLUSIONS: The NHS Health Check is linked to reduced incidence of disease across multiple organ systems, which may be attributed to risk modification through earlier detection and treatment of key risk factors such as hypertension and high cholesterol. This work adds important evidence to the growing body of research supporting the effectiveness of preventative interventions in reducing longer-term multimorbidity.

Spatial confidence regions for combinations of excursion sets in image analysis.

The analysis of excursion sets in imaging data is essential to a wide range of scientific disciplines such as neuroimaging, climatology, and cosmology. Despite growing literature, there is little published concerning the comparison of processes that have been sampled across the same spatial region but which reflect different study conditions. Given a set of asymptotically Gaussian random fields, each corresponding to a sample acquired for a different study condition, this work aims to provide confidence statements about the intersection, or union, of the excursion sets across all fields. Such spatial regions are of natural interest as they directly correspond to the questions 'Where do all random fields exceed a predetermined threshold?', or 'Where does at least one random field exceed a predetermined threshold?'. To assess the degree of spatial variability present, our method provides, with a desired confidence, subsets and supersets of spatial regions defined by logical conjunctions (i.e. set intersections) or disjunctions (i.e. set unions), without any assumption on the dependence between the different fields. The method is verified by extensive simulations and demonstrated using task-fMRI data to identify brain regions with activation common to four variants of a working memory task.

A framework for longitudinal latent factor modelling of treatment response in clinical trials with applications to Psoriatic Arthritis and Rheumatoid Arthritis.

OBJECTIVE: Clinical trials involve the collection of a wealth of data, comprising multiple diverse measurements performed at baseline and follow-up visits over the course of a trial. The most common primary analysis is restricted to a single, potentially composite endpoint at one time point. While such an analytical focus promotes simple and replicable conclusions, it does not necessarily fully capture the multi-faceted effects of a drug in a complex disease setting. Therefore, to complement existing approaches, we set out here to design a longitudinal multivariate analytical framework that accepts as input an entire clinical trial database, comprising all measurements, patients, and time points across multiple trials. METHODS: Our framework composes probabilistic principal component analysis with a longitudinal linear mixed effects model, thereby enabling clinical interpretation of multivariate results, while handling data missing at random, and incorporating covariates and covariance structure in a computationally efficient and principled way. RESULTS: We illustrate our approach by applying it to four phase III clinical trials of secukinumab in Psoriatic Arthritis (PsA) and Rheumatoid Arthritis (RA). We identify three clinically plausible latent factors that collectively explain 74.5% of empirical variation in the longitudinal patient database. We estimate longitudinal trajectories of these factors, thereby enabling joint characterisation of disease progression and drug effect. We perform benchmarking experiments demonstrating our method's competitive performance at estimating average treatment effects compared to existing statistical and machine learning methods, and showing that our modular approach leads to relatively computationally efficient model fitting. CONCLUSION: Our multivariate longitudinal framework has the potential to illuminate the properties of existing composite endpoint methods, and to enable the development of novel clinical endpoints that provide enhanced and complementary perspectives on treatment response.

Quality over quantity: powering neuroimaging samples in psychiatry.

Neuroimaging has been widely adopted in psychiatric research, with hopes that these non-invasive methods will provide important clues to the underpinnings and prediction of various mental health symptoms and outcomes. However, the translational impact of neuroimaging has not yet reached its promise, despite the plethora of computational methods, tools, and datasets at our disposal. Some have lamented that too many psychiatric neuroimaging studies have been underpowered with respect to sample size. In this review, we encourage this discourse to shift from a focus on sheer increases in sample size to more thoughtful choices surrounding experimental study designs. We propose considerations at multiple decision points throughout the study design, data modeling and analysis process that may help researchers working in psychiatric neuroimaging boost power for their research questions of interest without necessarily increasing sample size. We also provide suggestions for leveraging multiple datasets to inform each other and strengthen our confidence in the generalization of findings to both population-level and clinical samples. Through a greater emphasis on improving the quality of brain-based and clinical measures rather than merely quantity, meaningful and potentially translational clinical associations with neuroimaging measures can be achieved with more modest sample sizes in psychiatry.

Bayesian Lesion Estimation with a Structured Spike-and-Slab Prior.

Neural demyelination and brain damage accumulated in white matter appear as hyperintense areas on T2-weighted MRI scans in the form of lesions. Modeling binary images at the population level, where each voxel represents the existence of a lesion, plays an important role in understanding aging and inflammatory diseases. We propose a scalable hierarchical Bayesian spatial model, called BLESS, capable of handling binary responses by placing continuous spike-and-slab mixture priors on spatially-varying parameters and enforcing spatial dependency on the parameter dictating the amount of sparsity within the probability of inclusion. The use of mean-field variational inference with dynamic posterior exploration, which is an annealing-like strategy that improves optimization, allows our method to scale to large sample sizes. Our method also accounts for underestimation of posterior variance due to variational inference by providing an approximate posterior sampling approach based on Bayesian bootstrap ideas and spike-and-slab priors with random shrinkage targets. Besides accurate uncertainty quantification, this approach is capable of producing novel cluster size based imaging statistics, such as credible intervals of cluster size, and measures of reliability of cluster occurrence. Lastly, we validate our results via simulation studies and an application to the UK Biobank, a large-scale lesion mapping study with a sample size of 40,000 subjects.

Decentralized Mixed Effects Modeling in COINSTAC.

Performing group analysis on magnetic resonance imaging (MRI) data with linear mixed-effects (LME) models is challenging due to its large dimensionality and inherent multi-level covariance structure. In addition, as large-scale collaborative projects become commonplace in neuroimaging, data must increasingly be stored and analyzed from different locations. In such settings, substantial overhead can occur in terms of data transfer and coordination between participating research groups. In some cases, data cannot be pooled together due to privacy or regulatory concerns. In this work, we propose a decentralized LME model to perform a large-scale analysis of data from different collaborations without data pooling. This method is efficient as it overcomes the hurdles of data sharing and has lower bandwidth and memory requirements for analysis than the centralized modeling approach. We evaluate our model using features extracted from structural magnetic resonance imaging (sMRI) data. Results highlight gray matter reductions in the temporal lobe/insula and medial frontal regions in schizophrenia, consistent with prior studies. Our analysis also demonstrates that decentralized LME models achieve similar performance compared to the models trained with all the data in one location. We also implement the decentralized LME approach in COINSTAC, an open source, decentralized platform for federating neuroimaging analysis, providing an easy to use tool for dissemination to the neuroimaging community.

Ventricular volume asymmetry as a novel imaging biomarker for disease discrimination and outcome prediction.

Aims: Disruption of the predictable symmetry of the healthy heart may be an indicator of cardiovascular risk. This study defines the population distribution of ventricular asymmetry and its relationships across a range of prevalent and incident cardiorespiratory diseases. Methods and results: The analysis includes 44 796 UK Biobank participants (average age 64.1 ± 7.7 years; 51.9% women). Cardiovascular magnetic resonance (CMR) metrics were derived using previously validated automated pipelines. Ventricular asymmetry was expressed as the ratio of left and right ventricular (LV and RV) end-diastolic volumes. Clinical outcomes were defined through linked health records. Incident events were those occurring for the first time after imaging, longitudinally tracked over an average follow-up time of 4.75 ± 1.52 years. The normal range for ventricular symmetry was defined in a healthy subset. Participants with values outside the 5th-95th percentiles of the healthy distribution were classed as either LV dominant (LV/RV > 112%) or RV dominant (LV/RV < 80%) asymmetry. Associations of LV and RV dominant asymmetry with vascular risk factors, CMR features, and prevalent and incident cardiovascular diseases (CVDs) were examined using regression models, adjusting for vascular risk factors, prevalent diseases, and conventional CMR measures. Left ventricular dominance was linked to an array of pre-existing vascular risk factors and CVDs, and a two-fold increased risk of incident heart failure, non-ischaemic cardiomyopathies, and left-sided valvular disorders. Right ventricular dominance was associated with an elevated risk of all-cause mortality. Conclusion: Ventricular asymmetry has clinical utility for cardiovascular risk assessment, providing information that is incremental to traditional risk factors and conventional CMR metrics.

Two common and distinct forms of variation in human functional brain networks.

The cortex has a characteristic layout with specialized functional areas forming distributed large-scale networks. However, substantial work shows striking variation in this organization across people, which relates to differences in behavior. While most previous work treats individual differences as linked to boundary shifts between the borders of regions, here we show that cortical 'variants' also occur at a distance from their typical position, forming ectopic intrusions. Both 'border' and 'ectopic' variants are common across individuals, but differ in their location, network associations, properties of subgroups of individuals, activations during tasks, and prediction of behavioral phenotypes. Border variants also track significantly more with shared genetics than ectopic variants, suggesting a closer link between ectopic variants and environmental influences. This work argues that these two dissociable forms of variation-border shifts and ectopic intrusions-must be separately accounted for in the analysis of individual differences in cortical systems across people.

Feasibility of multiorgan risk prediction with routinely collected diagnostics: a prospective cohort study in the UK Biobank.

OBJECTIVES: Despite rising rates of multimorbidity, existing risk assessment tools are mostly limited to a single outcome of interest. This study tests the feasibility of producing multiple disease risk estimates with at least 70% discrimination (area under the receiver operating curve, AUROC) within the time and information constraints of the existing primary care health check framework. DESIGN: Observational prospective cohort study SETTING: UK Biobank. PARTICIPANTS: 228 240 adults from the UK population. INTERVENTIONS: None. MAIN OUTCOME MEASURES: Myocardial infarction, atrial fibrillation, heart failure, stroke, all-cause dementia, chronic kidney disease, fatty liver disease, alcoholic liver disease, liver cirrhosis and liver failure. RESULTS: Using a set of predictors easily gathered at the standard primary care health check (such as the National Health Service Health Check), we demonstrate that it is feasible to simultaneously produce risk estimates for multiple disease outcomes with AUROC of 70% or greater. These predictors can be entered once into a single form and produce risk scores for stroke (AUROC 0.727, 95% CI 0.713 to 0.740), all-cause dementia (0.823, 95% CI 0.810 to 0.836), myocardial infarction (0.785, 95% CI 0.775 to 0.795), atrial fibrillation (0.777, 95% CI 0.768 to 0.785), heart failure (0.828, 95% CI 0.818 to 0.838), chronic kidney disease (0.774, 95% CI 0.765 to 0.783), fatty liver disease (0.766, 95% CI 0.753 to 0.779), alcoholic liver disease (0.864, 95% CI 0.835 to 0.894), liver cirrhosis (0.763, 95% CI 0.734 to 0.793) and liver failure (0.746, 95% CI 0.695 to 0.796). CONCLUSIONS: Easily collected diagnostics can be used to assess 10-year risk across multiple disease outcomes, without the need for specialist computing or invasive biomarkers. Such an approach could increase the utility of existing data and place multiorgan risk information at the fingertips of primary care providers, thus creating opportunities for longer-term multimorbidity prevention. Additional work is needed to validate whether these findings would hold in a larger, more representative cohort outside the UK Biobank.

MRI economics: Balancing sample size and scan duration in brain wide association studies.

A pervasive dilemma in neuroimaging is whether to prioritize sample size or scan duration given fixed resources. Here, we systematically investigate this trade-off in the context of brain-wide association studies (BWAS) using resting-state functional magnetic resonance imaging (fMRI). We find that total scan duration (sample size × scan duration per participant) robustly explains individual-level phenotypic prediction accuracy via a logarithmic model, suggesting that sample size and scan duration are broadly interchangeable. The returns of scan duration eventually diminish relative to sample size, which we explain with principled theoretical derivations. When accounting for fixed costs associated with each participant (e.g., recruitment, non-imaging measures), we find that prediction accuracy in small-scale BWAS might benefit from much longer scan durations (>50 min) than typically assumed. Most existing large-scale studies might also have benefited from smaller sample sizes with longer scan durations. Both logarithmic and theoretical models of the relationships among sample size, scan duration and prediction accuracy explain well-predicted phenotypes better than poorly-predicted phenotypes. The logarithmic and theoretical models are also undermined by individual differences in brain states. These results replicate across phenotypic domains (e.g., cognition and mental health) from two large-scale datasets with different algorithms and metrics. Overall, our study emphasizes the importance of scan time, which is ignored in standard power calculations. Standard power calculations inevitably maximize sample size at the expense of scan duration. The resulting prediction accuracies are likely lower than would be produced with alternate designs, thus impeding scientific discovery. Our empirically informed reference is available for future study design: WEB_APPLICATION_LINK.

Erratum: Separating Clinical and Subclinical Depression by Big Data Informed Structural Vulnerability Index and Its impact on Cognition: ENIGMA Dot Product.

In the PSB article published in Biocomputing 2022: Proceedings of the Pacific Symposium, pp. 133-143; doi: 10.1142/9789811250477_0013 (https://www.ncbi.nlm.nih.gov/pmc/articles/PMC8719281/), the following author name is missing: Si Gao MS