A structural covariance network (SCN) has been used successfully in structural magnetic resonance imaging (sMRI) studies. However, most SCNs have been constructed by a unitary marker that is insensitive for discriminating different disease phases. The aim of this study was to devise a novel regional radiomics similarity network (R2SN) that could provide more comprehensive information in morphological network analysis. R2SNs were constructed by computing the Pearson correlations between the radiomics features extracted from any pair of regions for each subject (AAL atlas). We further assessed the small-world property of R2SNs, and we evaluated the reproducibility in different datasets and through test-retest analysis. The relationships between the R2SNs and general intelligence/interregional coexpression of genes were also explored. R2SNs could be replicated in different datasets, regardless of the use of different feature subsets. R2SNs showed high reproducibility in the test-retest analysis (intraclass correlation coefficient > 0.7). In addition, the small-word property (σ > 2) and the high correlation between gene expression (R = 0.29, p < 0.001) and general intelligence were determined for R2SNs. Furthermore, the results have also been repeated in the Brainnetome atlas. R2SNs provide a novel, reliable, and biologically plausible method to understand human morphological covariance based on sMRI.

Gray matter volume and cortical thickness are some of the most popular brain morphological measures of structural magnetic resonance imaging (sMRI). These patterns are important for understanding complex brain cognitive function. However, most of the studies typically analyze single/several anatomical regions independently without considering associations among brain regions. The structural covariance network (SCN) is often used to reconstruct the brain structural network from sMRI and is commonly used to measure the association between regions in the human brain with morphological similarity. However, most of the individual SCNs have been constructed by a unitary marker such as gray volume/cortical thickness with hyposensitivity. We develop a novel, reliable and biologically plausible brain network to understand human morphological covariance based on sMRI.

Structural magnetic resonance imaging (sMRI) plays an important role in neuroscience, including in evaluations of gray matter volume and cortical thickness, which are some of the most popular brain morphological measures. However, most of the studies typically analyze single/several anatomical regions independently without considering associations among brain regions (Alexander-Bloch et al., 2013; Pichet Binette et al., 2020), especially regarding which complex heterogeneous network patterns can be used to characterize the brain by supporting information transformation; these patterns are important for understanding complex brain cognitive function (Bullmore & Sporns, 2012; Alexander-Bloch et al., 2013). Specifically, the structural covariance network (SCN) is often used to reconstruct the brain structural network from sMRI based on the similarity of gray matter morphology (He et al., 2007; Tijms et al., 2012) and is commonly used to measure the association between interregions in the human brain with morphological similarity (He et al., 2007; Montembeault et al., 2012; Tijms et al., 2012; Montembeault et al., 2016; Zhang et al., 2017; Seidlitz et al., 2018; Kong et al., 2019; Spreng et al., 2019).

The SCNs typically consist of nodes and edges, representing the predefined brain regions and the statistical similarity between them based on predefined morphological markers such as volume or thickness (For a review, see Alexander-Bloch et al., 2013). Several methodologies have been introduced for the reconstruction of connectome maps based on sMRI either at the group level [such as independent component analysis (Guo et al., 2015; Pichet Binette et al., 2020) and seed-based technology (Zielinski et al., 2010)] or the individual level. Specifically, SCN based on cortical thickness (He et al., 2007, 2008) or gray matter volume (Yao et al., 2010) was well studied at the group level, whereas the abovementioned biomarkers were also developed to construct brain networks using single or multiple morphological features at the individual level (Tijms et al., 2012; Wee et al., 2013; Kong et al., 2014; Kim et al., 2016). All these methods have been used to investigate network alterations in brain-related diseases (Yao et al., 2010; Zheng et al., 2015; Bethlehem et al., 2017; Yu et al., 2018). Seidlitz and colleagues proposed a morphometric similarity network that captures cortical cytoarchitecture and is linked to individual cognitive performance (Seidlitz et al., 2018), and it has been applied to understand major depressive disorder (Li et al., 2021). Despite the progress in constructing different brain networks, a well-validated and widely accessible model for mapping the brain network architecture of anatomical brain regions in an individual human brain is needed. Radiomics is a powerful, robust method to extract more detailed information (Li et al., 2019), including intensity and texture features from each brain region (Parmar et al., 2015; Gillies et al., 2016; Chaddad et al., 2018). Texture analysis describes a variety of features that quantify the variation in the patterns of intensity, including some imperceptible information to the human visual system (Kassner & Thornhill, 2010). However, there is a literature gap regarding the construction of a radiomics-based similarity network, as well as the associated attributes, which could be a feasible anatomical topological mapping of the individual brain.

In this study, the first aim was to develop a novel regional radiomics similarity network (R2SN) approach. Building on this foundation, the second aim was to explore the reproducibility, small-world property, and biological basis of the R2SNs, including the relationship between R2SN indices and the coexpression of the gene or general fluid intelligence (gF) scores. The results confirmed that R2SNs provide a novel, robust, and biologically plausible model and a new perspective for understanding the human brain; therefore, R2SNs have great promise in further studies.

Subjects

A total of 848 subjects from the Human Connectome Project (HCP, https://www.humanconnectome.org/study/hcp-young-adult/document/) were included in our study; all subjects were cognitively normal controls (NCs; age: 28.82 ± 3.68, sex (M/F): 371/477, fluid intelligence: 16.53 ± 4.86, Flanker inhibitory control and attention (FICA) test score: 112.41 ± 10.06). The HCP was established in 2009 with an overarching objective of studying human brain connectivity and its variability in healthy adults (Van Essen et al., 2012). All HCP subjects were evaluated by 3T MR scanners. In the HCP protocol, fluid intelligence was assessed using a form of Raven’s progressive matrices with 24 items. Detailed subject information can be found at https://www.humanconnectome.org/study/hcp-young-adult/document/1200-subjects-data-release and can also be found in a previous study (Van Essen et al., 2012).

Data Preprocessing and Radiomics Feature Extraction

For each subject, the T1-weighted MRI image was aligned to Montreal Neurological Institute (MNI) space by using a combined linear and nonlinear registration (including N4 bias field correction) and resampled to 1 mm × 1 mm × 1 mm (Xie et al., 2016) (Figure 1A). Then, 47 radiomics features in each brain region were extracted with each region defined in the AAL atlas (a total of 90 regions) (Tzourio-Mazoyer et al., 2002). The radiomics features consisted of 14 intensity features and 33 texture features (Figure 1B). All features were described in the study by Aerts et al. (2014) and implemented as in-house MATLAB scripts (https://github.com/YongLiulab/). The definitions and detailed descriptions of the radiomics features can be found in previous publications (Aerts et al., 2014; Feng et al., 2018; Zhao et al., 2020) and Supporting Information section S01. Redundancy features, defined as features having a high correlation with other features (R > 0.9), were removed before subsequent analysis (Supporting Information section S02). Therefore, a final feature matrix with 25 × 90 for each individual was obtained for further analysis (Figure 1C; Table 1).

Figure 1.

Schematic of the data analysis pipeline. (A) Data preprocessing. (B) The computation of radiomics features in each brain region. (C) Feature matrix of radiomics features. (D) Network construction with Pearson correlation. (E) The reproducibility of the network, (F) the network parameters of the R2SN, (G) the correlation between the R2SN and gene expression, and (H) the correlation between the R2SN and the gF score/FICA test score. gF indicates fluid intelligence; FICA indicates Flanker inhibitory control and attention.

Figure 1.

Schematic of the data analysis pipeline. (A) Data preprocessing. (B) The computation of radiomics features in each brain region. (C) Feature matrix of radiomics features. (D) Network construction with Pearson correlation. (E) The reproducibility of the network, (F) the network parameters of the R2SN, (G) the correlation between the R2SN and gene expression, and (H) the correlation between the R2SN and the gF score/FICA test score. gF indicates fluid intelligence; FICA indicates Flanker inhibitory control and attention.

Close modal
Table 1.

The reserved features after removing superfluous features

Intensity features energy Textural features Autocorrelation IMC1 
kurtosis Cluster prominence Maximum probability 
maximum Cluster Shade Sum entropy 
MAD Cluster Tendency Short run-length emphasis 
minimum Contrast Long run-length emphasis 
skewness Correlation Gray-level nonuniformity 
entropy Energy Low gray-level run-length emphasis 
  Entropy High gray-level run-length emphasis 
  Homogeneity1 High gray-level long run-length emphasis 
Intensity features energy Textural features Autocorrelation IMC1 
kurtosis Cluster prominence Maximum probability 
maximum Cluster Shade Sum entropy 
MAD Cluster Tendency Short run-length emphasis 
minimum Contrast Long run-length emphasis 
skewness Correlation Gray-level nonuniformity 
entropy Energy Low gray-level run-length emphasis 
  Entropy High gray-level run-length emphasis 
  Homogeneity1 High gray-level long run-length emphasis 

R2SN Construction

An individual R2SN was constructed by feature normalization, followed by radiomics similarity matrix establishment. Specifically, feature normalization was implemented by adopting a common min-max feature normalization scheme, whereas the radiomics similarity matrix was established by mapping the individual’s radiomics features into a radiomics similarity matrix of pairwise interregional Pearson’s correlations (Figure 1D). Briefly, for each subject, a 90 × 90 connection matrix was obtained, the node of this network was defined as the region based on the AAL atlas, and the edge was calculated by Pearson’s correlations between interregional radiomics features (Figure 1D). The mean value and the standard deviation (Std) of the R2SNs were computed to estimate the fluctuation of R2SNs in these young NCs (Figure 1E).

The Topological Structure of R2SNs

To explore the topological structure of the brain characterized by R2SNs, a variety of graph-theoretical network parameters were computed, including shortest path length (L), clustering coefficient (C), and small-world property, after binarization of R2SNs by using a threshold ranging from 0.5 to 0.75 (step size = 0.01). Briefly, the small-world index was defined as sigma (σ) = γ/λ, with gamma being (γ) = CR2SN/Crandom, lambda being (λ) = LR2SN/Lrandom, and Crandom and Lrandom being the C and L of the random network (random sampling of edges to yield a matrix with the same number of connections), respectively (Figure 1F). In addition, we computed the neighbor degree, degree, betweenness centrality, and global efficiency of the R2SNs by using a threshold ranging from 0.5 to 0.75 (step size = 0.01). The detailed definitions of these parameters can be found in a previous study (Rubinov & Sporns, 2010) and were computed by using the Brainnetome fMRI Toolkit (BRANT) (Xu et al., 2018) (Figure 1F).

Reproducibility of the Network

We demonstrated three points of view to estimate the reproducibility of R2SNs, including the reproducibility of R2SN among different datasets, the reproducibility of the individual R2SN among different scans, and the robustness of R2SNs to methodological variations. First, the HCP dataset was randomly divided into two subdatasets for 1,000 iterations (424 subjects for each subdataset). The Pearson correlation coefficient of the two subdatasets in each iteration was used to estimate the consistency of the mR2SN (Seidlitz et al., 2018; Jin et al., 2020a; Zhao et al., 2020), which was defined as the mean value of the R2SNs in the specific dataset. Finally, the distribution of 1,000 Pearson correlation coefficients was used to estimate the reproducibility of the R2SNs (Seidlitz et al., 2018; Jin et al., 2020a; Zhao et al., 2020).

Notably, the test-retest analysis was used to estimate the reproducibility in different observers when measuring the same index (Friedman et al., 2008; Leijenaar et al., 2013; Maclaren et al., 2014; Noble et al., 2019). To investigate the reproducibility of the R2SN, a test-retest analysis was conducted using 21 subjects with four images acquired during different visits (https://duke.edu/∼morey005/ScanRescanData/). In this dataset, each subject was scanned on two different days, two scans were conducted 1 hr apart on day 1 (scans 1A and 1B), and two scans were conducted 1 hr apart at a second session 7–9 days later (scans 2A and 2B) (Morey et al., 2010). The intraclass correlation coefficient [ICC; ICC = (BMS − WMS)/BMS] was used to estimate the reproducibility of each edge of the R2SNs, where BMS is the between-subjects mean square and WMS is the within-subject mean square. The ICC has a value between 0 and 1; ICC = 0 indicates no reproducibility, and ICC = 1 indicates absolute reproducibility (Shrout & Fleiss, 1979) (Figure 1E).

We also quantified the robustness of R2SNs to methodological variations, including randomly reducing the number of radiomics features for analysis (i.e., only 20 radiomics features rather than all 25 features were involved as a predefined marker for network construction). The Pearson correlation coefficient was used to estimate the similarity of the mR2SN, which was constructed with 20 radiomics features and 25 radiomics features. The distribution of the Pearson correlation coefficient of 1,000 simulations was performed to assess the robustness of R2SNs to methodological variations (Seidlitz et al., 2018; Jin et al., 2020a; Zhao et al., 2020) (Figure 1E).

In addition, we explored whether a significant correlation could be obtained between the mean connective strength and the size of each node (the size of each node was roughly estimated with the voxel number based on the AAL atlas) (Tzourio-Mazoyer et al., 2002).

The Relationship Between the R2SN and the Gene Similarity Network

To further explore the biological basis of R2SNs, we continued our investigation by computing the relationship between the R2SN and the gene similarity network (GSN). The GSN was constructed with the Allen atlas (https://human.brain-map.org/) (Zeng et al., 2012) and predefined genes (six subjects). The nodes were defined as the Allen atlas and mapped to the AAL brain regions based on MNI coordinates by using “abagen” (https://github.com/rmarkello/abagen) (Arnatkeviciute et al., 2019), and the edges were computed by the Pearson coefficient between the gene expression of any pair of regions. To correct the spatial autocorrelation effect within R2SN and GSN, the Euclidean distance of each pair of brain regions was employed as a concomitant variable. Pearson’s correlation coefficient between the connectivities of the mR2SN and the connectivities of the mGSN (the mean value of the GSN) was calculated to assess the similarity between the R2SN and GSN (Figure 1G). Furthermore, Pearson’s correlation coefficient between the mean connection strengths of nodes of the mR2SN and mGSN was calculated to evaluate the similarity of these two networks.

Association Between the Network Properties of the R2SN and Cognitive Differences

We have investigated the network properties, reproducibility, and biological basis of the R2SN, and we assumed that the individual network properties could represent the individual’s differences in cognitive ability. In the HCP project, fluid intelligence scores were obtained as an index for measuring the subjects’ intrinsic cognitive ability. In addition, the Flanker task measures both a participant’s attention and inhibitory control. Therefore, we focused on the associations of common gF and the FICA test score as the indices to estimate the cognitive difference among young NCs, similar to Seidlitz et al. (2018). On this basis, we performed Pearson’s correlation analysis between gF/FICA scores and connectivity strength and network parameters (the neighbor degree, degree, betweenness centrality, global efficiency, shortest path length, and clustering coefficients) of the individual’s network to explore the relationship between the R2SN and gF/FICA scores in the HCP dataset.

The Replicability of Results in Brainnetome Atlas

To explore the replicability of the results in different atlases, we repeated the above analysis in the Brainnetome atlas (Fan et al., 2016). First, the R2SN was created based on the Brainnetome atlas (N = 246 without cerebellar), and then the reproducibility, small-world properties and the biological basis of R2SN were explored with the same pipeline.

R2SN: Edge Properties

After redundancy removal, 25 radiomics features were reserved for each brain region, and the R2SN was constructed with those predefined features for each subject (a symmetrical matrix with a size of 90 × 90). The connective strength of the mR2SN ranged from −0.56 to 0.99 in the HCP dataset (N = 848), and the Std value in each edge of the R2SN was significantly smaller than the mean value, which ranged from 0.01 to 0.36, but with more than 95% of them ranging from 0.01 to 0.20 (Figure 2A).

Figure 2.

The reproducibility of the R2SN. (A) The distribution of variance in the HCP dataset. (B) The distribution of the correlation coefficients of the R2SN network for each pair of two datasets (the HCP dataset was divided into two parts for 1,000 random repetitions). (C) The distribution of the correlation coefficients of the R2SN network with 20 features (random selection of 20 features from 25 predefined features 1,000 times) and the R2SN created with 25 features. (D) The ICC value of the R2SN with test-retest analysis.

Figure 2.

The reproducibility of the R2SN. (A) The distribution of variance in the HCP dataset. (B) The distribution of the correlation coefficients of the R2SN network for each pair of two datasets (the HCP dataset was divided into two parts for 1,000 random repetitions). (C) The distribution of the correlation coefficients of the R2SN network with 20 features (random selection of 20 features from 25 predefined features 1,000 times) and the R2SN created with 25 features. (D) The ICC value of the R2SN with test-retest analysis.

Close modal

R2SN: Network Properties

The C, L, and sparseness values are shown in Figure 3AC. In addition, the lambda (λ) value showed the ratio of the L of the R2SN and the random network, and the gamma (γ) value showed the ratio of the C of the R2SN and the random network. As a result, the value of γ was significantly larger than 1 (Figure 3D), the value of λ was close to 1 (Figure 3E), and the small-world index sigma (σ) was also significantly larger than 1 by different thresholds of binarization (from 0.5 to 0.75 with a step size = 0.01) (Figure 3F). In addition, the betweenness centrality (Figure 3G), neighbor degree (Figure 3H), degree (Figure 3I) and global efficiency (Figure 3J) of the R2SNs with different thresholds (from 0.5 to 0.75 with a step size = 0.01) are also shown in Figure 3.

Figure 3.

The network parameters of the R2SN with different correlation thresholds (0.5–0.75 and step size = 0.01). (A) The clustering coefficient—the “black” means the R2SN, and the “blue” means the random network; (B) the shortest path length; (C) sparseness, and the small-world parameter, including (D) the gamma value (the ratio of clustering coefficient between R2SN and random network); (E) the lambda value (the ratio of shortest path length between the R2SN and the random network); (F) the sigma value (the ratio of gamma and lambda); (G) the betweenness centrality value; (H) the neighbor degree value; (I) the degree value; and (J) the global efficiency value.

Figure 3.

The network parameters of the R2SN with different correlation thresholds (0.5–0.75 and step size = 0.01). (A) The clustering coefficient—the “black” means the R2SN, and the “blue” means the random network; (B) the shortest path length; (C) sparseness, and the small-world parameter, including (D) the gamma value (the ratio of clustering coefficient between R2SN and random network); (E) the lambda value (the ratio of shortest path length between the R2SN and the random network); (F) the sigma value (the ratio of gamma and lambda); (G) the betweenness centrality value; (H) the neighbor degree value; (I) the degree value; and (J) the global efficiency value.

Close modal

Reproducibility of the R2SN

A series of hypothesis testing was performed to further assess the reproducibility of the R2SN. Most notably, high consistency was found in any two mR2SNs constructed by different datasets (1,000 random simulations), with the Pearson coefficient ranging from 0.9997 to 1 (Figure 2B). In addition, high consistency was also obtained by mR2SNs constructed with a different number of features (20 randomly selected features and all features), and the R value ranged from 0.879 to 0.997 (Figure 2C). More importantly, the R2SN had a high ICC value (ICC > 0.7) within more than a 95% edge by test-retest analysis (Figure 2D). We did not find a significant correlation between the size of the node and the mean connective strength (R = 0.14, p = 0.17).

The Association Between the R2SN and the GSN

For each subject from the Allen dataset, the GSN was constructed based on 15,633 predefined genes (https://human.brain-map.org/). The mean connective strength of each node of the mR2SN and the mGSN are shown in Figure 4A (mR2SN) and Figure 4B (mGSN). We also computed the similarity between the mR2SN and the mGSN (edge-based), and a significant correlation was found between the two networks with R = 0.29 (p < 0.001) (Figure 4C and Figure 4D), meaning that the brain region with high morphometric similarity also tended to have a high transcriptional similarity of the gene. We also computed the average of the off-diagonal elements of a row or column in the radiomics similarity matrix (node-based), and a significant correlation was also found between the mR2SN and the mGSN regarding the mean connective strength of each node (R = 0.32, p = 0.002).

Figure 4.

The correlation between the R2SN and the gene expression network. (A) The mean connective strength of the R2SN, which was mapped to the surface area from the AAL template. (B) The mean connective strength of the GSN, which was mapped to the surface area from the AAL template. The value of the color bar was normalized with the max-min method. (C) The heat map for the R2SN and GSN. Some negative correlations were generated with the GSN, and this phenomenon might have been caused by the deletion of genes in some brain regions. (D) The scatter diagram of the correlation between GSN and R2SN, where the Euclidean distance between each pair of regions of interest was employed as a concomitant variable within the network.

Figure 4.

The correlation between the R2SN and the gene expression network. (A) The mean connective strength of the R2SN, which was mapped to the surface area from the AAL template. (B) The mean connective strength of the GSN, which was mapped to the surface area from the AAL template. The value of the color bar was normalized with the max-min method. (C) The heat map for the R2SN and GSN. Some negative correlations were generated with the GSN, and this phenomenon might have been caused by the deletion of genes in some brain regions. (D) The scatter diagram of the correlation between GSN and R2SN, where the Euclidean distance between each pair of regions of interest was employed as a concomitant variable within the network.

Close modal

The Association Between the R2SN and Cognitive Differences

The fluid intelligence scores ranged from 4 to 24, and the FICA test scores ranged from 89 to 142 for all subjects (https://www.humanconnectome.org/study/hcp-young-adult). The Pearson’s correlation showed that approximately 5% of connections had a significant correlation with the fluid intelligence score (Bonferroni-corrected p < 0.05, with N = 4,005) (Figure 5A). As Figure 5 shows, the neighbor degree of 15 nodes (Figure 5B), the degree of 16 nodes (Figure 5C), the global efficiency of 20 nodes (Figure 5D), and the clustering coefficient of 8 nodes (Figure 5E) were significantly correlated with fluid intelligence (Bonferroni-corrected p < 0.001, with N = 90) (when the threshold = 0.50).

Figure 5.

The correlation between the R2SN and cognitive ability. (A) Connectivity showed a significant correlation with the individuals’ gF score (Bonferroni-corrected p < 0.001); the brain regions in which the (B) neighbor degree value, (C) degree value, (D) global efficiency value, and (E) clustering coefficient showed a significant association with the gF score (Bonferroni-corrected p < 0.001). (F) Connectivity showed a significant correlation with the individuals’ FICA score (Bonferroni-corrected p < 0.001); the brain regions in which the (B) neighbor degree, (C) degree, (D) global efficiency, and (E) clustering coefficient showed a significant association with the FICA score (Bonferroni-corrected p < 0.001). FICA indicates the Flanker inhibitory control and attention test score, and gF indicates the fluid intelligence.

Figure 5.

The correlation between the R2SN and cognitive ability. (A) Connectivity showed a significant correlation with the individuals’ gF score (Bonferroni-corrected p < 0.001); the brain regions in which the (B) neighbor degree value, (C) degree value, (D) global efficiency value, and (E) clustering coefficient showed a significant association with the gF score (Bonferroni-corrected p < 0.001). (F) Connectivity showed a significant correlation with the individuals’ FICA score (Bonferroni-corrected p < 0.001); the brain regions in which the (B) neighbor degree, (C) degree, (D) global efficiency, and (E) clustering coefficient showed a significant association with the FICA score (Bonferroni-corrected p < 0.001). FICA indicates the Flanker inhibitory control and attention test score, and gF indicates the fluid intelligence.

Close modal

Pearson’s correlation showed that approximately 8% of connections had a significant correlation with the FICA test score (Bonferroni-corrected p < 0.05, with N = 4,005) (Figure 5F). In addition, the neighbor degree of 7 nodes (Figure 5G), degree of 7 nodes (Figure 5H), global efficiency of 42 nodes (Figure 5I), and clustering coefficient of 38 nodes (Figure 5J) were significantly correlated with the FICA test score (Bonferroni-corrected p < 0.001, with N = 90) (when the threshold is 0.50). The correlation between gF/FICA and the shortest path length and betweenness centrality is shown in Supporting Information section S03 due to the weaker significance.

The Replicability of the Results in the Brainnetome Atlas

The results can be repeated in the Brainnetome atlas. First, high consistency was found in any two mR2SNs constructed by different datasets or subsets of the feature. More importantly, the R2SN had a high ICC value (ICC > 0.7) within more than 95% of edges by the test-retest analysis (Supporting Information Figure S1). Furthermore, the value of γ was significantly larger than 1, the value of λ was close to 1 and the small-world index sigma (σ) was also significantly larger than 1 by different thresholds of binarization (from 0.5 to 0.75 with a step size = 0.01) (Supporting Information Figure S2). In addition, the relationships between R2SN and gene expression (Supporting Information Figure S3)/fluid intelligence (Supporting Information Figure S4) were also verified in the Brainnetome atlas (Supporting Information S03).

The present study provides a pipeline for transforming radiomics feature maps into a pairwise interregional radiomics similarity matrix at the individual level, which still represents a gap according to the literature. The systematic results confirmed that R2SNs provide a novel, robust, and biologically plausible model for understanding the human brain.

Radiomics is a powerful method to extract more detailed information from brain images, which includes intensity and texture features, and which could improve the ability of the SCN to represent the morphology of the brain as a network with the relativeness of pair of regions (Parmar et al., 2015; Gillies et al., 2016; Chaddad et al., 2018). Texture features can quantify the variations in intensity or patterns, including those features that are imperceptible to the human visual system (Aerts et al., 2014). Numerous studies have emphasized the importance of radiomics (Gillies et al., 2016) and have regarded it as a bridge between imaging and personalized medicine (Lambin et al., 2017). More importantly, radiomics also improve the diagnostic precision, treatment, and prognosis of a tumor (Aerts et al., 2014; Coroller et al., 2015; Huang et al., 2018) as well as the diagnosis of Alzheimer’s disease (Sorensen et al., 2016; Sorensen et al., 2017; Feng et al., 2018; Zhao et al., 2020; Ding et al., 2021). Radiomics is a method that extracts a large number of features of the brain region with high reproducibility (Li et al., 2019) and provides a better measure to evaluate the characteristics of the brain region.

Representing the morphology of the brain as a network has an advantage in that the structure of the brain can be described statistically with tools from graph theory (Tijms et al., 2012). Similar to most of the brain networks (e.g., functional brain networks, diffusion-weighted imaging, and SCNs) (He et al., 2007; Tijms et al., 2012; Heinze et al., 2015; Seidlitz et al., 2018), R2SNs also had complex topologies. In R2SNs, the high-degree nodes were located in the frontal lobe, parietal lobe, and occipital lobe, and low-degree nodes were found in the temporal lobe and subcortical lobe. The nodes with strong connection strengths might work more cooperatively with other brain regions, and the nodes with lower connection strengths might have a more specific function in the brain (Seidlitz et al., 2018). The different connective patterns between a GSN and R2SN might be caused by the difference in microstructure in the brain, such as the histological classification of cortical areas (Solari & Stoner, 2011; Seidlitz et al., 2018).

Notably, R2SN can be replicated in different datasets, demonstrated by a subsampling strategy from the same dataset for 1,000 iterations. R2SNs showed high ICCs (a prominent statistic to measure test-retest reliability) in different visit images, indicating the robustness of the radiomics similarity connectivity under the consensus that reproducibility was the most important property for a novel method in MRI analysis (Stikov et al., 2019). The high consistency of R2SNs, which was constructed with different features, also confirmed the methodological variations (Seidlitz et al., 2018). These properties of R2SNs support the firm foundation for the credibility of the results.

Imaging biomarkers are taken as the cornerstone of the radiology community, and imaging genetics has established heritable phenotypes for quantitative genetics of brain phenotypes. As expected, a high correlation was found between the gene expression network and the R2SN (Zeng et al., 2012), meaning that cortical areas with high morphometric similarity also tended to have high transcriptional similarity (Seidlitz et al., 2018). The structure of the brain region was controlled by gene expression (Warden & Mayfield, 2017), and the variations in intensity or pattern of the brain region can be reflected by radiomics features (Wang et al., 2019). These results further suggested that gene expression can be reflected by R2SNs. In brief, R2SNs have a genetic basis, and the present findings provide the possibility of estimating the risk of a gene by using R2SNs and indicate a certain degree of evidence for genetic disease research, such as that conducted regarding Alzheimer’s disease (Jin et al., 2020b; Zhao et al., 2020).

It is vital to evaluate the association between the individual network architecture and the cognitive ability or psychological functions of the brain (Li et al., 2009; van den Heuvel et al., 2009; Seidlitz et al., 2018). For high-degree hub nodes and the global efficiency of the connectome to be preferentially affected by clinical brain disorders associated with cognitive impairment, the relationship between high-degree nodes and cognitive status should obtain a higher performance in classification, prediction, and so on (Crossley et al., 2014). As we all know that fluid intelligence is a common measure to estimate the cognitive ability among NCs and to identify the cognitive difference among individuals, higher fluid intelligence scores corresponded to more efficient information transfer in the brain (Li et al., 2009). The fluid intelligence score refers to the ability to solve abstract problems that do not depend on acquired knowledge and change with age (Gray et al., 2003; Cole et al., 2012; Kievit et al., 2016). A set of fronto-parietal brain regions has been reported to be associated with fluid intelligence in brain imaging (Duncan et al., 2000; Cole et al., 2012; Woolgar et al., 2018). Comprehensive evidence has indicated that the lingual gyrus, caudate nucleus, rolandic operculum, and frontal lobe play an important role in an individual’s cognitive ability (Preusse et al., 2011; Rhein et al., 2014; Santarnecchi et al., 2015; Zhang et al., 2016). The significant correlation between fluid intelligence and efficiency of the hub regions in R2SNs also indicates that an individual’s cognitive ability is linked to the brain network architecture (Li et al., 2009; van den Heuvel et al., 2009). In addition, the network measures exclude global efficiency, which is also correlated with cognitive ability. Thus, we speculate that the R2SN might be a brain network with a solid biological basis.

This study has some limitations. First, we only studied R2SNs based on the AAL atlas and Brainnetome atlas; therefore, the network properties and their basis need to be further validated by other reputable brain atlases. Second, a unified framework for interpreting these measures and their alterations in different brain diseases is needed. The integration of existing studies shows that different models for brain structural networks do not align uniformly across the brain, and the coupling between structural network models and/or functional connectivity remodeling will help support the underpinnings of functional specialization and cognition. Third, more samples from different independent scanners and more cognition measures may improve the statistical power of the analysis, allowing scientists to explore the neural mechanisms of R2SNs in the future.

R2SN is a network with high reproducibility and a biological basis; thus, an R2SN might serve as an improved, novel method and shed new light on future MRI studies. We assume that an R2SN could provide a powerful technology platform for measuring the anatomical connectome in vivo and be applied for the diagnosis of a variety of diseases in the future.

Data were provided by the Human Connectome Project, WU-Minn Consortium (Principal Investigators: David Van Essen and Kamil Ugurbil; 1U54MH091657) funded by the 16 National Institutes of Health institutes and centers that support the NIH Blueprint for Neuroscience Research and by the McDonnell Center for Systems Neuroscience at Washington University.

Kun Zhao: Conceptualization; Data curation; Formal analysis; Investigation; Methodology; Visualization; Writing – original draft. Qiang Zheng: Funding acquisition; Investigation; Writing – review & editing. Tongtong Che: Writing – review & editing. Martin Dyrba: Writing – review & editing. Qiongling Li: Writing – review & editing. Yanhui Ding: Funding acquisition; Writing – review & editing. Yuanjie Zheng: Funding acquisition; Writing – review & editing. Yong Liu: Conceptualization; Funding acquisition; Writing – original draft; Writing – review & editing. Shuyu Li: Conceptualization; Funding acquisition; Writing – original draft; Writing – review & editing.

Supporting information for this article is available at https://doi.org/10.1162/netn_a_00200. All subjects of this study are downloaded from the Human Connectome Project (HCP, https://www.humanconnectome.org/study/hcp-young-adult/document/), and the scripts of the radiomics features are available at https://github.com/YongLiulab.

Shuyu Li, the National Natural Science Foundation of China, Award ID: 81972160, 81622025. Yong Liu, the National Natural Science Foundation of China, Award ID: 81871438, the Beijing Natural Science Funds for Distinguished Young Scholars, Award ID: JQ20036. Qiang Zheng, the National Natural Science Foundation of China, Award ID: 61802330. Yuanjie Zheng, the National Natural Science Foundation of China, Award ID: 81871508, 61773246, the Taishan Scholar Program of Shandong Province of China, Award ID: TSHW201502038. Yanhui Ding, the Natural Science Foundation of Shandong Province, Award ID: ZR2020MF051.

Structural covariance network (SCN):

A brain network from sMRI based on the similarity of gray matter morphology and is commonly used to measure the association between interregions in the human brain with morphological similarity.

Radiomics:

A powerful, robust method to extract more detailed information from each brain region.

Small world network:

Network with high clustering coefficient and short shortest path length.

Fluid intelligence:

A measure which can reflect the cognitive ability.

Flanker inhibitory control and attention test:

A measure which can reflect both a participant’s attention and inhibitory control.

Shortest path length:

A shortest path between two nodes in a graph is a path with the minimum number of edge.

Sigma (σ):

The small-world index was defined as sigma (σ) = γ/λ, with gamma being (γ) = CR2SN/Crandom, lambda being (λ) = LR2SN/Lrandom, and with Crandom and Lrandom being the C and L of the random network. C: clustering coefficient; L: shortest path length.

Random network:

Network with low clustering coefficient and short shortest path length.

Global efficiency:

The global efficiency is the average inverse shortest path length in the network.

Gene similarity network:

A network based on the similarity of gene expression between each pair of brain regions.

Aerts
,
H. J.
,
Velazquez
,
E. R.
,
Leijenaar
,
R. T.
,
Parmar
,
C.
,
Grossmann
,
P.
,
Carvalho
,
S.
, …
Lambin
,
P.
(
2014
).
Decoding tumour phenotype by noninvasive imaging using a quantitative radiomics approach
.
Nature Communications
,
5
,
4006
. ,
[PubMed]
Alexander-Bloch
,
A.
,
Giedd
,
J. N.
, &
Bullmore
,
E.
(
2013
).
Imaging structural co-variance between human brain regions
.
Nature Reviews Neuroscience
,
14
(
5
),
322
336
. ,
[PubMed]
Arnatkeviciute
,
A.
,
Fulcher
,
B. D.
, &
Fornito
,
A.
(
2019
).
A practical guide to linking brain-wide gene expression and neuroimaging data
.
NeuroImage
,
189
,
353
367
. ,
[PubMed]
Bethlehem
,
R. A. I.
,
Romero-Garcia
,
R.
,
Mak
,
E.
,
Bullmore
,
E. T.
, &
Baron-Cohen
,
S.
(
2017
).
Structural covariance networks in children with autism or ADHD
.
Cerebral Cortex
,
27
(
8
),
4267
4276
. ,
[PubMed]
Bullmore
,
E.
, &
Sporns
,
O.
(
2012
).
The economy of brain network organization
.
Nature Reviews Neuroscience
,
13
(
5
),
336
349
. ,
[PubMed]
Chaddad
,
A.
,
Desrosiers
,
C.
, &
Niazi
,
T.
(
2018
).
Deep radiomic analysis of MRI related to Alzheimer’s disease
.
IEEE Access
,
6
,
58213
58221
.
Cole
,
M. W.
,
Yarkoni
,
T.
,
Repovs
,
G.
,
Anticevic
,
A.
, &
Braver
,
T. S.
(
2012
).
Global connectivity of prefrontal cortex predicts cognitive control and intelligence
.
Journal of Neuroscience
,
32
(
26
),
8988
8999
. ,
[PubMed]
Coroller
,
T. P.
,
Grossmann
,
P.
,
Hou
,
Y.
,
Rios Velazquez
,
E.
,
Leijenaar
,
R. T.
,
Hermann
,
G.
, …
Aerts
,
H. J.
(
2015
).
CT-based radiomic signature predicts distant metastasis in lung adenocarcinoma
.
Radiotherapy and Oncology
,
114
(
3
),
345
350
. ,
[PubMed]
Crossley
,
N. A.
,
Mechelli
,
A.
,
Scott
,
J.
,
Carletti
,
F.
,
Fox
,
P. T.
,
McGuire
,
P.
, &
Bullmore
,
E. T.
(
2014
).
The hubs of the human connectome are generally implicated in the anatomy of brain disorders
.
Brain
,
137
(
Pt 8
),
2382
2395
. ,
[PubMed]
Ding
,
Y.
,
Zhao
,
K.
,
Che
,
T.
,
Du
,
K.
,
Su
,
H.
,
Liu
,
S.
, …
Liu
,
Y.
(
2021
).
Quantitative radiomic features as new biomarkers for Alzheimer’s disease: An amyloid PET study
.
Cerebral Cortex
. ,
[PubMed]
Duncan
,
J.
,
Seitz
,
R. J.
,
Kolodny
,
J.
,
Bor
,
D.
,
Herzog
,
H.
,
Ahmed
,
A.
, …
Emslie
,
H.
(
2000
).
A neural basis for general intelligence
.
Science
,
289
(
5478
),
457
460
. ,
[PubMed]
Fan
,
L.
,
Li
,
H.
,
Zhuo
,
J.
,
Zhang
,
Y.
,
Wang
,
J.
,
Chen
,
L.
, …
Jiang
,
T.
(
2016
).
The Human Brainnetome Atlas: A new brain atlas based on connectional architecture
.
Cerebral Cortex
,
26
(
8
),
3508
3526
. ,
[PubMed]
Feng
,
F.
,
Wang
,
P.
,
Zhao
,
K.
,
Zhou
,
B.
,
Yao
,
H.
,
Meng
,
Q.
, …
Liu
,
Y.
(
2018
).
Radiomic features of hippocampal subregions in Alzheimer’s disease and amnestic mild cognitive impairment
.
Frontiers in Aging Neuroscience
,
10
,
290
. ,
[PubMed]
Friedman
,
L.
,
Stern
,
H.
,
Brown
,
G. G.
,
Mathalon
,
D. H.
,
Turner
,
J.
,
Glover
,
G. H.
, …
Potkin
,
S. G.
(
2008
).
Test-retest and between-site reliability in a multicenter fMRI study
.
Human Brain Mapping
,
29
(
8
),
958
972
. ,
[PubMed]
Gillies
,
R. J.
,
Kinahan
,
P. E.
, &
Hricak
,
H.
(
2016
).
Radiomics: Images are more than pictures, they are data
.
Radiology
,
278
(
2
),
563
577
. ,
[PubMed]
Gray
,
J. R.
,
Chabris
,
C. F.
, &
Braver
,
T. S.
(
2003
).
Neural mechanisms of general fluid intelligence
.
Nature Neuroscience
,
6
(
3
),
316
322
. ,
[PubMed]
Guo
,
X.
,
Wang
,
Y.
,
Guo
,
T.
,
Chen
,
K.
,
Zhang
,
J.
,
Li
,
K.
, …
Yao
,
L.
(
2015
).
Structural covariance networks across healthy young adults and their consistency
.
Journal of Magnetic Resonance Imaging
,
42
(
2
),
261
268
. ,
[PubMed]
He
,
Y.
,
Chen
,
Z.
, &
Evans
,
A.
(
2008
).
Structural insights into aberrant topological patterns of large-scale cortical networks in Alzheimer’s disease
.
Journal of Neuroscience
,
28
(
18
),
4756
4766
. ,
[PubMed]
He
,
Y.
,
Chen
,
Z. J.
, &
Evans
,
A. C.
(
2007
).
Small-world anatomical networks in the human brain revealed by cortical thickness from MRI
.
Cerebral Cortex
,
17
(
10
),
2407
2419
. ,
[PubMed]
Heinze
,
K.
,
Reniers
,
R. L.
,
Nelson
,
B.
,
Yung
,
A. R.
,
Lin
,
A.
,
Harrison
,
B. J.
, …
Wood
,
S. J.
(
2015
).
Discrete alterations of brain network structural covariance in individuals at ultra-high risk for psychosis
.
Biological Psychiatry
,
77
(
11
),
989
996
. ,
[PubMed]
Huang
,
X.
,
Cheng
,
Z.
,
Huang
,
Y.
,
Liang
,
C.
,
He
,
L.
,
Ma
,
Z.
, …
Liu
,
Z.
(
2018
).
CT-based radiomics signature to discriminate high-grade from low-grade colorectal adenocarcinoma
.
Academic Radiology
,
25
(
10
),
1285
1297
. ,
[PubMed]
Jin
,
D.
,
Wang
,
P.
,
Zalesky
,
A.
,
Liu
,
B.
,
Song
,
C.
,
Wang
,
D.
, …
Liu
,
Y.
(
2020a
).
Grab-AD: Generalizability and reproducibility of altered brain activity and diagnostic classification in Alzheimer’s disease
.
Human Brain Mapping
,
41
(
12
),
3379
3391
. ,
[PubMed]
Jin
,
D.
,
Zhou
,
B.
,
Han
,
Y.
,
Ren
,
J.
,
Han
,
T.
,
Liu
,
B.
, …
Liu
,
Y.
(
2020b
).
Generalizable, reproducible, and neuroscientifically interpretable imaging biomarkers for Alzheimer’s disease
.
Advanced Science (Weinh)
,
7
(
14
),
2000675
. ,
[PubMed]
Kassner
,
A.
, &
Thornhill
,
R. E.
(
2010
).
Texture analysis: A review of neurologic MR imaging applications
.
American Journal of Neuroradiology
,
31
(
5
),
809
816
. ,
[PubMed]
Kievit
,
R. A.
,
Davis
,
S. W.
,
Griffiths
,
J.
,
Correia
,
M. M.
,
Cam
,
C.
, &
Henson
,
R. N.
(
2016
).
A watershed model of individual differences in fluid intelligence
.
Neuropsychologia
,
91
,
186
198
. ,
[PubMed]
Kim
,
H. J.
,
Shin
,
J. H.
,
Han
,
C. E.
,
Kim
,
H. J.
,
Na
,
D. L.
,
Seo
,
S. W.
, …
the Alzheimer’s Disease Neuroimaging Initiative
. (
2016
).
Using individualized brain network for analyzing structural covariance of the cerebral cortex in Alzheimer’s patients
.
Frontiers in Neuroscience
,
10
,
394
. ,
[PubMed]
Kong
,
R.
,
Li
,
J.
,
Orban
,
C.
,
Sabuncu
,
M. R.
,
Liu
,
H.
,
Schaefer
,
A.
, …
Yeo
,
B. T. T.
(
2019
).
Spatial topography of individual-specific cortical networks predicts human cognition, personality, and emotion
.
Cerebral Cortex
,
29
(
6
),
2533
2551
. ,
[PubMed]
Kong
,
X. Z.
,
Wang
,
X.
,
Huang
,
L.
,
Pu
,
Y.
,
Yang
,
Z.
,
Dang
,
X.
, …
Liu
,
J.
(
2014
).
Measuring individual morphological relationship of cortical regions
.
Journal of Neuroscience Methods
,
237
,
103
107
. ,
[PubMed]
Lambin
,
P.
,
Leijenaar
,
R. T. H.
,
Deist
,
T. M.
,
Peerlings
,
J.
,
de Jong
,
E. E. C.
,
van Timmeren
,
J.
, …
Walsh
,
S.
(
2017
).
Radiomics: The bridge between medical imaging and personalized medicine
.
Nature Reviews Clinical Oncology
,
14
(
12
),
749
762
. ,
[PubMed]
Leijenaar
,
R. T.
,
Carvalho
,
S.
,
Velazquez
,
E. R.
,
van Elmpt
,
W. J.
,
Parmar
,
C.
,
Hoekstra
,
O. S.
, …
Lambin
,
P.
(
2013
).
Stability of FDG-PET radiomics features: An integrated analysis of test-retest and inter-observer variability
.
Acta Oncologica
,
52
(
7
),
1391
1397
. ,
[PubMed]
Li
,
J.
,
Seidlitz
,
J.
,
Suckling
,
J.
,
Fan
,
F.
,
Ji
,
G. J.
,
Meng
,
Y.
, …
Liao
,
W.
(
2021
).
Cortical structural differences in major depressive disorder correlate with cell type-specific transcriptional signatures
.
Nature Communications
,
12
(
1
),
1647
. ,
[PubMed]
Li
,
Y.
,
Liu
,
Y.
,
Li
,
J.
,
Qin
,
W.
,
Li
,
K.
,
Yu
,
C.
, &
Jiang
,
T.
(
2009
).
Brain anatomical network and intelligence
.
PLoS Computational Biology
,
5
(
5
),
e1000395
. ,
[PubMed]
Li
,
Z. R.
,
Duan
,
H. C.
,
Zhao
,
K.
, &
Ding
,
Y. H.
(
2019
).
Stability of MRI radiomics features of hippocampus: An integrated analysis of test-retest and inter-observer variability
.
IEEE Access
,
7
,
97106
97116
.
Maclaren
,
J.
,
Han
,
Z.
,
Vos
,
S. B.
,
Fischbein
,
N.
, &
Bammer
,
R.
(
2014
).
Reliability of brain volume measurements: A test-retest dataset
.
Scientific Data
,
1
,
140037
. ,
[PubMed]
Montembeault
,
M.
,
Joubert
,
S.
,
Doyon
,
J.
,
Carrier
,
J.
,
Gagnon
,
J. F.
,
Monchi
,
O.
, …
Brambati
,
S. M.
(
2012
).
The impact of aging on gray matter structural covariance networks
.
NeuroImage
,
63
(
2
),
754
759
. ,
[PubMed]
Montembeault
,
M.
,
Rouleau
,
I.
,
Provost
,
J. S.
,
Brambati
,
S. M.
, &
the Alzheimer’s Disease Neuroimaging Initiative
. (
2016
).
Altered gray matter structural covariance networks in early stages of Alzheimer’s disease
.
Cerebral Cortex
,
26
(
6
),
2650
2662
. ,
[PubMed]
Morey
,
R. A.
,
Selgrade
,
E. S.
,
Wagner
,
H. R.
, 2nd
,
Huettel
,
S. A.
,
Wang
,
L.
, &
McCarthy
,
G.
(
2010
).
Scan-rescan reliability of subcortical brain volumes derived from automated segmentation
.
Human Brain Mapping
,
31
(
11
),
1751
1762
. ,
[PubMed]
Noble
,
S.
,
Scheinost
,
D.
, &
Constable
,
R. T.
(
2019
).
A decade of test-retest reliability of functional connectivity: A systematic review and meta-analysis
.
NeuroImage
,
203
,
116157
. ,
[PubMed]
Parmar
,
C.
,
Leijenaar
,
R. T.
,
Grossmann
,
P.
,
Rios Velazquez
,
E.
,
Bussink
,
J.
,
Rietveld
,
D.
, …
Aerts
,
H. J.
(
2015
).
Radiomic feature clusters and prognostic signatures specific for lung and head and neck cancer
.
Scientific Reports
,
5
,
11044
. ,
[PubMed]
Pichet Binette
,
A.
,
Gonneaud
,
J.
,
Vogel
,
J. W.
,
La Joie
,
R.
,
Rosa-Neto
,
P.
,
Collins
,
D. L.
, …
the PREVENT-AD Research Group
. (
2020
).
Morphometric network differences in ageing versus Alzheimer’s disease dementia
.
Brain
,
143
(
2
),
635
649
. ,
[PubMed]
Preusse
,
F.
,
van der Meer
,
E.
,
Deshpande
,
G.
,
Krueger
,
F.
, &
Wartenburger
,
I.
(
2011
).
Fluid intelligence allows flexible recruitment of the parieto-frontal network in analogical reasoning
.
Frontiers in Human Neuroscience
,
5
,
22
. ,
[PubMed]
Rhein
,
C.
,
Muhle
,
C.
,
Richter-Schmidinger
,
T.
,
Alexopoulos
,
P.
,
Doerfler
,
A.
, &
Kornhuber
,
J.
(
2014
).
Neuroanatomical correlates of intelligence in healthy young adults: The role of basal ganglia volume
.
PLoS One
,
9
(
4
),
e93623
. ,
[PubMed]
Rubinov
,
M.
, &
Sporns
,
O.
(
2010
).
Complex network measures of brain connectivity: Uses and interpretations
.
NeuroImage
,
52
(
3
),
1059
1069
. ,
[PubMed]
Santarnecchi
,
E.
,
Tatti
,
E.
,
Rossi
,
S.
,
Serino
,
V.
, &
Rossi
,
A.
(
2015
).
Intelligence-related differences in the asymmetry of spontaneous cerebral activity
.
Human Brain Mapping
,
36
(
9
),
3586
3602
. ,
[PubMed]
Seidlitz
,
J.
,
Vasa
,
F.
,
Shinn
,
M.
,
Romero-Garcia
,
R.
,
Whitaker
,
K. J.
,
Vertes
,
P. E.
, …
Bullmore
,
E. T.
(
2018
).
Morphometric similarity networks detect microscale cortical organization and predict inter-individual cognitive variation
.
Neuron
,
97
(
1
),
231
247
. ,
[PubMed]
Shrout
,
P. E.
, &
Fleiss
,
J. L.
(
1979
).
Intraclass correlations: Uses in assessing rater reliability
.
Psychological Bulletin
,
86
(
2
),
420
428
. ,
[PubMed]
Solari
,
S. V.
, &
Stoner
,
R.
(
2011
).
Cognitive consilience: Primate non-primary neuroanatomical circuits underlying cognition
.
Frontiers in Neuroanatomy
,
5
,
65
. ,
[PubMed]
Sorensen
,
L.
,
Igel
,
C.
,
Liv Hansen
,
N.
,
Osler
,
M.
,
Lauritzen
,
M.
,
Rostrup
,
E.
, …
the Australian Imaging Biomarkers and Lifestyle Flagship Study of Ageing
. (
2016
).
Early detection of Alzheimer’s disease using MRI hippocampal texture
.
Human Brain Mapping
,
37
(
3
),
1148
1161
. ,
[PubMed]
Sorensen
,
L.
,
Igel
,
C.
,
Pai
,
A.
,
Balas
,
I.
,
Anker
,
C.
,
Lillholm
,
M.
, …
the Australian Imaging Biomarkers and Lifestyle Flagship Study of Ageing
. (
2017
).
Differential diagnosis of mild cognitive impairment and Alzheimer’s disease using structural MRI cortical thickness, hippocampal shape, hippocampal texture, and volumetry
.
NeuroImage Clinical
,
13
,
470
482
. ,
[PubMed]
Spreng
,
R. N.
,
DuPre
,
E.
,
Ji
,
J. L.
,
Yang
,
G.
,
Diehl
,
C.
,
Murray
,
J. D.
, …
Anticevic
,
A.
(
2019
).
Structural covariance reveals alterations in control and salience network integrity in chronic schizophrenia
.
Cerebral Cortex
,
29
(
12
),
5269
5284
. ,
[PubMed]
Stikov
,
N.
,
Trzasko
,
J. D.
, &
Bernstein
,
M. A.
(
2019
).
Reproducibility and the future of MRI research
.
Magnetic Resonance in Medicine
,
82
(
6
),
1981
1983
. ,
[PubMed]
Tijms
,
B. M.
,
Series
,
P.
,
Willshaw
,
D. J.
, &
Lawrie
,
S. M.
(
2012
).
Similarity-based extraction of individual networks from gray matter MRI scans
.
Cerebral Cortex
,
22
(
7
),
1530
1541
. ,
[PubMed]
Tzourio-Mazoyer
,
N.
,
Landeau
,
B.
,
Papathanassiou
,
D.
,
Crivello
,
F.
,
Etard
,
O.
,
Delcroix
,
N.
, …
Joliot
,
M.
(
2002
).
Automated anatomical labeling of activations in SPM using a macroscopic anatomical parcellation of the MNI MRI single-subject brain
.
NeuroImage
,
15
(
1
),
273
289
. ,
[PubMed]
van den Heuvel
,
M. P.
,
Stam
,
C. J.
,
Kahn
,
R. S.
, &
Hulshoff Pol
,
H. E.
(
2009
).
Efficiency of functional brain networks and intellectual performance
.
Journal of Neuroscience
,
29
(
23
),
7619
7624
. ,
[PubMed]
Van Essen
,
D. C.
,
Ugurbil
,
K.
,
Auerbach
,
E.
,
Barch
,
D.
,
Behrens
,
T. E.
,
Bucholz
,
R.
, …
WU-Minn HCP Consortium
. (
2012
).
The Human Connectome Project: A data acquisition perspective
.
NeuroImage
,
62
(
4
),
2222
2231
. ,
[PubMed]
Wang
,
T.
,
Gong
,
J.
,
Duan
,
H. H.
,
Wang
,
L. J.
,
Ye
,
X. D.
, &
Nie
,
S. D.
(
2019
).
Correlation between CT based radiomics features and gene expression data in non-small cell lung cancer
.
Journal of X-Ray Science and Technology
,
27
(
5
),
773
803
. ,
[PubMed]
Warden
,
A. S.
, &
Mayfield
,
R. D.
(
2017
).
Gene expression profiling in the human alcoholic brain
.
Neuropharmacology
,
122
,
161
174
. ,
[PubMed]
Wee
,
C. Y.
,
Yap
,
P. T.
,
Shen
,
D.
, &
the Alzheimer’s Disease Neuroimaging Initiative
. (
2013
).
Prediction of Alzheimer’s disease and mild cognitive impairment using cortical morphological patterns
.
Human Brain Mapping
,
34
(
12
),
3411
3425
. ,
[PubMed]
Woolgar
,
A.
,
Duncan
,
J.
,
Manes
,
F.
, &
Fedorenko
,
E.
(
2018
).
The multiple-demand system but not the language system supports fluid intelligence
.
Nature Human Behaviour
,
2
(
3
),
200
204
. ,
[PubMed]
Xie
,
S.
,
Chen
,
L.
,
Zuo
,
N.
, &
Jiang
,
T.
(
2016
).
DiffusionKit: A light one-stop solution for diffusion MRI data analysis
.
Journal of Neuroscience Methods
,
273
,
107
119
. ,
[PubMed]
Xu
,
K.
,
Liu
,
Y.
,
Zhan
,
Y.
,
Ren
,
J.
, &
Jiang
,
T.
(
2018
).
BRANT: A versatile and extendable resting-state fMRI toolkit
.
Frontiers in Neuroinformatics
,
12
(
52
),
18191
18196
. ,
[PubMed]
Yao
,
Z.
,
Zhang
,
Y.
,
Lin
,
L.
,
Zhou
,
Y.
,
Xu
,
C.
,
Jiang
,
T.
, &
the Alzheimer’s Disease Neuroimaging Initiative
. (
2010
).
Abnormal cortical networks in mild cognitive impairment and Alzheimer’s disease
.
PLoS Computational Biology
,
6
(
11
),
e1001006
. ,
[PubMed]
Yu
,
K.
,
Wang
,
X.
,
Li
,
Q.
,
Zhang
,
X.
,
Li
,
X.
, &
Li
,
S.
(
2018
).
Individual morphological brain network construction based on multivariate Euclidean distances between brain regions
.
Frontiers in Human Neuroscience
,
12
,
204
. ,
[PubMed]
Zeng
,
H.
,
Shen
,
E. H.
,
Hohmann
,
J. G.
,
Oh
,
S. W.
,
Bernard
,
A.
,
Royall
,
J. J.
, …
Jones
,
A. R.
(
2012
).
Large-scale cellular-resolution gene profiling in human neocortex reveals species-specific molecular signatures
.
Cell
,
149
(
2
),
483
496
. ,
[PubMed]
Zhang
,
L.
,
Qiao
,
L.
,
Chen
,
Q.
,
Yang
,
W.
,
Xu
,
M.
,
Yao
,
X.
, …
Yang
,
D.
(
2016
).
Gray matter volume of the lingual gyrus mediates the relationship between inhibition function and divergent thinking
.
Frontiers in Psychology
,
7
,
1532
. ,
[PubMed]
Zhang
,
Z.
,
Liao
,
W.
,
Xu
,
Q.
,
Wei
,
W.
,
Zhou
,
H. J.
,
Sun
,
K.
, …
Lu
,
G.
(
2017
).
Hippocampus-associated causal network of structural covariance measuring structural damage progression in temporal lobe epilepsy
.
Human Brain Mapping
,
38
(
2
),
753
766
. ,
[PubMed]
Zhao
,
K.
,
Ding
,
Y. H.
,
Han
,
Y.
,
Fan
,
Y.
,
Alexander-Bloch
,
A. F.
,
Han
,
T.
, …
Co
,
M.-C. A. D. I.
(
2020
).
Independent and reproducible hippocampal radiomic biomarkers for multisite Alzheimer’s disease: Diagnosis, longitudinal progress and biological basis
.
Science Bulletin
,
65
(
13
),
1103
1113
.
Zheng
,
W.
,
Yao
,
Z.
,
Hu
,
B.
,
Gao
,
X.
,
Cai
,
H.
, &
Moore
,
P.
(
2015
).
Novel cortical thickness pattern for accurate detection of Alzheimer’s disease
.
Journal of Alzheimer’s Disease
,
48
(
4
),
995
1008
. ,
[PubMed]
Zielinski
,
B. A.
,
Gennatas
,
E. D.
,
Zhou
,
J.
, &
Seeley
,
W. W.
(
2010
).
Network-level structural covariance in the developing brain
.
Proceedings of the National Academy of Sciences of the United States of America
,
107
(
42
),
18191
18196
. ,
[PubMed]

Author notes

Competing Interests: The authors have declared that no competing interests exist.

Handling Editor: Vince Calhoun

This is an open-access article distributed under the terms of the Creative Commons Attribution 4.0 International License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited. For a full description of the license, please visit https://creativecommons.org/licenses/by/4.0/legalcode.

Supplementary data