New Tools For Intermodal Analysis And Association Testing In Neuroimaging

In the field of neuroimage analysis two key goals are to understand the association of a high- dimensional imaging variable with a phenotype, and to understand relationships between several high-dimensional imaging variables. Several recent studies have shown that the standard “mass- univariate” methods to test an association of an image with a phenotype have inflated type 1 error rates due to invalid assumptions. Here, we propose two new methods to perform association testing in neuroimaging and illustrate the method in two stages of the lifespan. The first is a para- metric bootstrap testing procedure that estimates the joint distribution of test statistical parametric map in order to control the voxel-wise family-wise error rate (FWER). We illustrate the method by identifying sex differences in nonlinear developmental trajectories of cerebral blood flow through adolescence using the Philadelphia Neurodevelopmental Cohort. The second testing procedure is a generalization of Rao’s score test based on projecting the score statistic onto a linear sub- space of a high-dimensional parameter space. The approach provides a way to localize signal in the high-dimensional space by projecting the scores to the subspace where the score test was performed. This allows for inference in the high-dimensional space to be performed on the same degrees of freedom as the score test, effectively reducing the number of comparisons. We illus- trate the method by analyzing a subset of the Alzheimer’s Disease Neuroimaging Initiative dataset. Finally, we propose a new tool to study relationships between neuroimaging modalities that we to describe how the spatial association between cortical thickness and sulcal depth changes in adolescent development

Analysing spatial point patterns in digital pathology: immune cells in high-grade serous ovarian carcinomas

Multiplex immunofluorescence (mIF) imaging technology facilitates the study of the tumour microenvironment in cancer patients. Due to the capabilities of this emerging bioimaging technique, it is possible to statistically analyse, for example, the co-varying location and functions of multiple different types of immune cells. Complex spatial relationships between different immune cells have been shown to correlate with patient outcomes and may reveal new pathways for targeted immunotherapy treatments. This tutorial reviews methods and procedures relating to spatial point patterns for complex data analysis. We consider tissue cells as a realisation of a spatial point process for each patient. We focus on proper functional descriptors for each observation and techniques that allow us to obtain information about inter-patient variation. Ovarian cancer is the deadliest gynaecological malignancy and can resist chemotherapy treatment effective in cancers. We use a dataset of high-grade serous ovarian cancer samples from 51 patients. We examine the immune cell composition (T cells, B cells, macrophages) within tumours and additional information such as cell classification (tumour or stroma) and other patient clinical characteristics. Our analyses, supported by reproducible software, apply to other digital pathology datasets

Faster Family-wise Error Control for Neuroimaging with a Parametric Bootstrap

In neuroimaging, hundreds to hundreds of thousands of tests are performed across a set of brain regions or all locations in an image. Recent studies have shown that the most common family-wise error (FWE) controlling procedures in imaging, which rely on classical mathematical inequalities or Gaussian random field theory, yield FWE rates that are far from the nominal level. Depending on the approach used, the FWER can be exceedingly small or grossly inflated. Given the widespread use of neuroimaging as a tool for understanding neurological and psychiatric disorders, it is imperative that reliable multiple testing procedures are available. To our knowledge, only permutation joint testing procedures have been shown to reliably control the FWER at the nominal level. However, these procedures are computationally intensive due to the increasingly available large sample sizes and dimensionality of the images, and analyses can take days to complete. Here, we develop a parametric bootstrap joint testing procedure. The parametric bootstrap procedure works directly with the test statistics, which leads to much faster estimation of adjusted \emph{p}-values than resampling-based procedures while reliably controlling the FWER in sample sizes available in many neuroimaging studies. We demonstrate that the procedure controls the FWER in finite samples using simulations, and present region- and voxel-wise analyses to test for sex differences in developmental trajectories of cerebral blood flow

Corticosteroids and regional variations in thickness of the human cerebral cortex across the lifespan

International audienceExposures to life stressors accumulate across the lifespan, with possible impact on brain health. Little is known, however, about the mechanisms mediating age-related changes in brain structure. We use a lifespan sample of participants (n = 21 251; 4–97 years) to investigate the relationship between the thickness of cerebral cortex and the expression of the glucocorticoid- and the mineralocorticoid-receptor genes (NR3C1 and NR3C2, respectively), obtained from the Allen Human Brain Atlas. In all participants, cortical thickness correlated negatively with the expression of both NR3C1 and NR3C2 across 34 cortical regions. The magnitude of this correlation varied across the lifespan. From childhood through early adulthood, the profile similarity (between NR3C1/NR3C2 expression and thickness) increased with age. Conversely, both profile similarities decreased with age in late life. These variations do not reflect age-related changes in NR3C1 and NR3C2 expression, as observed in 5 databases of gene expression in the human cerebral cortex (502 donors). Based on the co-expression of NR3C1 (and NR3C2) with genes specific to neural cell types, we determine the potential involvement of microglia, astrocytes, and CA1 pyramidal cells in mediating the relationship between corticosteroid exposure and cortical thickness. Therefore, corticosteroids may influence brain structure to a variable degree throughout life

Studying neuroanatomy using MRI

The study of neuroanatomy using imaging enables key insights into how our brains function, are shaped by genes and environment, and change with development, aging, and disease. Developments in MRI acquisition, image processing, and data modelling have been key to these advances. However, MRI provides an indirect measurement of the biological signals we aim to investigate. Thus, artifacts and key questions of correct interpretation can confound the readouts provided by anatomical MRI. In this review we provide an overview of the methods for measuring macro- and mesoscopic structure and inferring microstructural properties; we also describe key artefacts and confounds that can lead to incorrect conclusions. Ultimately, we believe that, though methods need to improve and caution is required in its interpretation, structural MRI continues to have great promise in furthering our understanding of how the brain works

Interpretable High-Dimensional Inference Via Score Projection With an Application in Neuroimaging

<p>In the fields of neuroimaging and genetics, a key goal is testing the association of a single outcome with a very high-dimensional imaging or genetic variable. Often, summary measures of the high-dimensional variable are created to sequentially test and localize the association with the outcome. In some cases, the associations between the outcome and summary measures are significant, but subsequent tests used to localize differences are underpowered and do not identify regions associated with the outcome. Here, we propose a generalization of Rao’s score test based on projecting the score statistic onto a linear subspace of a high-dimensional parameter space. The approach provides a way to localize signal in the high-dimensional space by projecting the scores to the subspace where the score test was performed. This allows for inference in the high-dimensional space to be performed on the same degrees of freedom as the score test, effectively reducing the number of comparisons. Simulation results demonstrate the test has competitive power relative to others commonly used. We illustrate the method by analyzing a subset of the Alzheimer’s Disease Neuroimaging Initiative dataset. Results suggest cortical thinning of the frontal and temporal lobes may be a useful biological marker of Alzheimer’s disease risk. Supplementary materials for this article are available online.</p