Resampling as a cluster validation technique in fMRI

Evidence of transgenerational effects on autism spectrum disorder using multigenerational space-time cluster detection

Background

Transgenerational epigenetic risks associated with complex health outcomes, such as autism spectrum disorder (ASD), have attracted increasing attention. Transgenerational environmental risk exposures with potential for epigenetic effects can be effectively identified using space-time clustering. Specifically applied to ancestors of individuals with disease outcomes, space-time clustering characterized for vulnerable developmental stages of growth can provide a measure of relative risk for disease outcomes in descendants.

Objectives

(1) Identify space-time clusters of ancestors with a descendent with a clinical ASD diagnosis and matched controls. (2) Identify developmental windows of ancestors with the highest relative risk for ASD in descendants. (3) Identify how the relative risk may vary through the maternal or paternal line.

Methods

Family pedigrees linked to residential locations of ASD cases in Utah have been used to identify space-time clusters of ancestors. Control family pedigrees of none-cases based on age and sex have been matched to cases 2:1. The data have been categorized by maternal or paternal lineage at birth, childhood, and adolescence. A total of 3957 children, both parents, and maternal and paternal grandparents were identified. Bernoulli space-time binomial relative risk (RR) scan statistic was used to identify clusters. Monte Carlo simulation was used for statistical significance testing.

Results

Twenty statistically significant clusters were identified. Thirteen increased RR (> 1.0) space-time clusters were identified from the maternal and paternal lines at a p-value < 0.05. The paternal grandparents carry the greatest RR (2.86–2.96) during birth and childhood in the 1950’s–1960, which represent the smallest size clusters, and occur in urban areas. Additionally, seven statistically significant clusters with RR < 1 were relatively large in area, covering more rural areas of the state.

Conclusion

This study has identified statistically significant space-time clusters during critical developmental windows that are associated with ASD risk in descendants. The geographic space and time clusters family pedigrees with over 3 + generations, which we refer to as a person’s geographic legacy, is a powerful tool for studying transgenerational effects that may be epigenetic in nature. Our novel use of space-time clustering can be applied to any disease where family pedigree data is available.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12942-022-00313-4.

Distinct neurometabolic profiles are evident in the anterior cingulate of young people with major psychiatric disorders

Currently, there are no validated neurobiological methods for distinguishing different pathophysiological pathways in young patients presenting in the early phases of major psychiatric disorders. Hence, treatments are delivered simply on the basis of their possible effects on nonspecific symptom constructs such as depression, cognitive change or psychotic symptoms. In this study, the ratios (relative to creatine) of key metabolites (N-acetyl aspartate, myoinositol, glutamate and glutathione) were measured with proton magnetic resonance spectroscopy ((1)H-MRS) within the anterior cingulate cortex of 88 young persons presenting with major mood or psychotic symptoms. We derived empirically (using a cluster analytical technique) three subgroups of subjects on the basis of their patterns of in vivo brain biochemistry. The three subgroups were distinguished (from each other) by all the four metabolites, in particular, glutathione and glutamate. By contrast, the groups could not be distinguished by differences in terms of other demographic, functional or clinical measures. We propose that this (1)H-MRS-based subclassification system could be used as the basis for much more specific tests of novel intervention strategies (notably, antioxidant and glutamatergic therapies) early in the course of major psychiatric disorders.

Improved temporal clustering analysis method applied to whole-brain data in acupuncture fMRI study

Temporal clustering analysis (TCA) has been proposed as a method for detecting the brain responses of a functional magnetic resonance imaging (fMRI) time series when the time and location of activation are completely unknown. But TCA is not suitable for treating the time series of the whole brain due to the existence of many inactive pixels. In theory, active pixels are located only in gray matter (GM). In this study, SPM2 was used to segment functional images into GM, white matter and cerebrospinal fluid, and only the pixels in GM were considered. Thus, most of inactive pixels are deleted, so that the sensitivity of TCA is greatly improved in the analysis of the whole brain. The same set of acupuncture fMRI data was treated using both conventional TCA and modified TCA (MTCA) for comparing their analytical ability. The results clearly show a significant improvement in the sensitivity achieved by MTCA.

Analyzing consistency of independent components: an fMRI illustration

Independent component analysis (ICA) is a powerful data-driven signal processing technique. It has proved to be helpful in, e.g., biomedicine, telecommunication, finance and machine vision. Yet, some problems persist in its wider use. One concern is the reliability of solutions found with ICA algorithms, resulting from the stochastic changes each time the analysis is performed. The consistency of the solutions can be analyzed by clustering solutions from multiple runs of bootstrapped ICA. Related methods have been recently published either for analyzing algorithmic stability or reducing the variability. The presented approach targets the extraction of additional information related to the independent components, by focusing on the nature of the variability. Practical implications are illustrated through a functional magnetic resonance imaging (fMRI) experiment.

An improved temporal clustering analysis method applied to whole-brain data in fMRI study

Temporal clustering analysis (TCA) has been proposed as a method to detect the brain responses of an fMRI time series when the time and location of the activation are completely unknown. But TCA is still incompetent in dealing with the time series of the whole brain due to the existence of many inactive pixels. If only active pixels are considered, the sensitivity of TCA will be improved greatly and it could be applied to the whole brain. In this study, some modifications were made to TCA to remove inactive pixels, and the applicability of the modified TCA to the whole brain was validated with a set of visual fMRI data. Based on the time series of the modified TCA, activations of the whole brain corresponding to the visual stimulation were detected. Compared with the previous TCA, the modified TCA method shows a significant improvement in the sensitivity to detect activation peaks of the whole brain.

Improved temporal clustering analysis method for detecting multiple response peaks in fMRI

PURPOSE

To develop an improved temporal clustering analysis (TCA) method for detecting multiple active peaks by running the method once.

MATERIALS AND METHODS

Two cases of simulation data and a set of actual fMRI data from nine subjects were used to compare the traditional TCA method with the new method, termed extremum TCA (ETCA). The first case of simulation data simulated event-related activation and block activation in one cerebral area, and the second case simulated event-related activation and block activation in two cerebral areas. An in vivo visual stimulating experiment was performed on a 1.5T MR scanner. All imaging data were processed using both traditional TCA and the new method.

RESULTS

The results of both the simulated and actual fMRI data show that the new method is more sensitive and exact than traditional TCA in detecting multiple response peaks.

CONCLUSION

The new method is effective in detecting multiple activations even when the timing and location of the brain activation are completely unknown.

Graphical-Model-based Morphometric Analysis

We propose a novel method for voxel-based morphometry (VBM), which we call Graphical-Model-based Morphometric Analysis (GAMMA), to identify morphological abnormalities automatically, and to find complex probabilistic associations among voxels in magnetic-resonance images and clinical variables. GAMMA is a fully automatic, nonparametric morphometric-analysis algorithm, with high sensitivity and specificity. It uses a Bayesian network to represent the associations among voxels and the function variable, and uses a contextual-clustering method based on a Markov random field to find clusters in which all voxels have similar associations with the function variable. We use loopy belief propagation to infer the unobserved label field and belief map. As opposed to voxel-based morphometric methods based on general linear models, GAMMA is capable of identifying nonlinear associations among the function variable and voxels. Compared with our previous approach, a Bayesian morphometry algorithm, GAMMA has greater sensitivity, specificity, and computational efficiency.

Controlling the false positive rate in fuzzy clustering using randomization: application to fMRI activation detection

Despite its potential advantages for fMRI analysis, fuzzy C-means (FCM) clustering suffers from limitations such as the need for a priori knowledge of the number of clusters, and unknown statistical significance and instability of the results. We propose a randomization-based method to control the false-positive rate and estimate statistical significance of the FCM results. Using this novel approach, we develop an fMRI activation detection method. The ability of the method in controlling the false-positive rate is shown by analysis of false positives in activation maps of resting-state fMRI data. Controlling the false-positive rate in FCM allows comparison of different fuzzy clustering methods, using different feature spaces, to other fMRI detection methods. In this article, using simulation and real fMRI data, we compare a novel feature space that takes the variability of the hemodynamic response function into account (HRF-based feature space) to the conventional cross-correlation analysis and FCM using the cross-correlation feature space. In both cases, the HRF-based feature space provides a greater sensitivity compared to the cross-correlation feature space and conventional cross-correlation analysis. Application of the proposed method to finger-tapping fMRI data, using HRF-based feature space, detected activation in sub-cortical regions, whereas both of the FCM with cross-correlation feature space and the conventional cross-correlation method failed to detect them.

Fuzzy cluster analysis of high-field functional MRI data

Functional magnetic resonance imaging (fMRI) based on blood-oxygen level dependent (BOLD) contrast today is an established brain research method and quickly gains acceptance for complementary clinical diagnosis. However, neither the basic mechanisms like coupling between neuronal activation and haemodynamic response are known exactly, nor can the various artifacts be predicted or controlled. Thus, modeling functional signal changes is non-trivial and exploratory data analysis (EDA) may be rather useful. In particular, identification and separation of artifacts as well as quantification of expected, i.e. stimulus correlated, and novel information on brain activity is important for both, new insights in neuroscience and future developments in functional MRI of the human brain. After an introduction on fuzzy clustering and very high-field fMRI we present several examples where fuzzy cluster analysis (FCA) of fMRI time series helps to identify and locally separate various artifacts. We also present and discuss applications and limitations of fuzzy cluster analysis in very high-field functional MRI: differentiate temporal patterns in MRI using (a) a test object with static and dynamic parts, (b) artifacts due to gross head motion artifacts. Using a synthetic fMRI data set we quantitatively examine the influences of relevant FCA parameters on clustering results in terms of receiver-operator characteristics (ROC) and compare them with a commonly used model-based correlation analysis (CA) approach. The application of FCA in analyzing in vivo fMRI data is shown for (a) a motor paradigm, (b) data from multi-echo imaging, and (c) a fMRI study using mental rotation of three-dimensional cubes. We found that differentiation of true "neural" from false "vascular" activation is possible based on echo time dependence and specific activation levels, as well as based on their signal time-course. Exploratory data analysis methods in general and fuzzy cluster analysis in particular may help to identify artifacts and add novel and unexpected information valuable for interpretation, classification and characterization of functional MRI data which can be used to design new data acquisition schemes, stimulus presentations, neuro(physio)logical paradigms, as well as to improve quantitative biophysical models.

Improved detection of time windows of brain responses in fMRI using modified temporal clustering analysis

Temporal clustering analysis (TCA) has been proposed recently as a method to detect time windows of brain responses in functional MRI (fMRI) studies when the timing and location of the activation are completely unknown. Modifications to the TCA technique are introduced in this report to further improve the sensitivity in detecting brain activation. The modified TCA is based on the integrated signal intensity of a temporal cluster at each time point, while the original TCA is based only on the size of a temporal cluster at each time point. A temporal cluster at each time point is defined, in both TCA methods, as a group of pixels reaching their maximum (or minimum) values at the same time. Both computer simulation and in vivo fMRI experiments have been performed. Compared with the original TCA, the modified TCA shows a significant improvement in the sensitivity to detect activation peaks for determining time windows of brain responses.

Characterization of BOLD activation in multi-echo fMRI data using fuzzy cluster analysis and a comparison with quantitative modeling

A combination of multiple gradient-echo imaging and exploratory data analysis (EDA), i.e. fuzzy cluster analysis (FCA), is proposed for separation and characterization of BOLD activation in single-shot spiral functional magnetic resonance imaging (fMRI) experiments at 3 T. Differentiation of functional activation using FCA is performed by clustering pixel signal changes (DeltaS) as a function of echo time (TE). Further vascular classification is supported by the localization of activation and the comparison with a single-exponential decay model. In some subjects, an additional indication for large vessels within a voxel was found as oscillation of the fMRI signal difference vs echo time (TE). Such large vessels may be separated from small vessel activation and, therefore, our proposed procedure might prove useful if a more specific functional localization is desired in fMRI. In addition to the signal change DeltaS, DeltaT(2)*/T(2)* is significantly different between activated regions. Averaged over all eight subjects DeltaT(2)* is 1.7 +/- 0.2 ms in ROIs with the highest signal change characterized as containing large vessels, whereas in ROIs corresponding to microvascular environment average DeltaT(2)* values are 0.8 +/- 0.1 ms.

Data mining in brain imaging

Data mining in brain imaging is proving to be an effective methodology for disease prognosis and prevention. This, together with the rapid accumulation of massive heterogeneous data sets, motivates the need for efficient methods that filter, clarify, assess, correlate and cluster brain-related information. Here, we present data mining methods that have been or could be employed in the analysis of brain images. These methods address two types of brain imaging data: structural and functional. We introduce statistical methods that aid the discovery of interesting associations and patterns between brain images and other clinical data. We consider several applications of these methods, such as the analysis of task-activation, lesion-deficit, and structure morphological variability; the development of probabilistic atlases; and tumour analysis. We include examples of applications to real brain data. Several data mining issues, such as that of method validation or verification, are also discussed.

Current awareness

In order to keep subscribers up-to-date with the latest developments in their field, John Wiley & Sons are providing a current awareness service in each issue of the journal. The bibliography contains newly published material in the field of NMR in biomedicine. Each bibliography is divided into 9 sections: 1 Books, Reviews ' Symposia; 2 General; 3 Technology; 4 Brain and Nerves; 5 Neuropathology; 6 Cancer; 7 Cardiac, Vascular and Respiratory Systems; 8 Liver, Kidney and Other Organs; 9 Muscle and Orthopaedic. Within each section, articles are listed in alphabetical order with respect to author. If, in the preceding period, no publications are located relevant to any one of these headings, that section will be omitted.