Skip to main content

ORIGINAL RESEARCH article

Front. Neuroimaging, 25 May 2022
Sec. Brain Imaging Methods
This article is part of the Research Topic Advances in Brain Functional and Structural Networks Modeling via Graph Theory View all 7 articles

A Riemannian Revisiting of Structure–Function Mapping Based on Eigenmodes

  • Centre Inria d'Université Côte d'Azur, Valbonne, France

Understanding the link between brain structure and function may not only improve our knowledge of brain organization, but also lead to better quantification of pathology. To quantify this link, recent studies have attempted to predict the brain's functional connectivity from its structural connectivity. However, functional connectivity matrices live in the Riemannian manifold of the symmetric positive definite space and a specific attention must be paid to operate on this appropriate space. In this work we investigated the implications of using a distance based on an affine invariant Riemannian metric in the context of structure–function mapping. Specifically, we revisit previously proposed structure–function mappings based on eigendecomposition and test them on 100 healthy subjects from the Human Connectome Project using this adapted notion of distance. First, we show that using this Riemannian distance significantly alters the notion of similarity between subjects from a functional point of view. We also show that using this distance improves the correlation between the structural and functional similarity of different subjects. Finally, by using a distance appropriate to this manifold, we demonstrate the importance of mapping function from structure under the Riemannian manifold and show in particular that it is possible to outperform the group average and the so–called glass ceiling on the performance of mappings based on eigenmodes.

1. Introduction

Current imaging technology allows us to observe the brain's connectivity both from a functional and structural point of view, for example using resting state functional magnetic resonance imaging (MRI) and diffusion MRI tractography, respectively. These two views are not decoupled and provide complementary information of the connectivity, and more general architecture, of the brain. Indeed, a current representation of the brain is that of a network, where the white matter substrates allows communication between distant brain regions (Sporns et al., 2005). Given this perspective, it is natural to try to understand and quantify the link between brain structure and function. One strategy to quantify this link is to predict, at least partially, the brain's function from its structure, so called structure–function mapping. In addition to the overarching goal of understanding the brain's organization, the ability to predict the brain's function from its structure has the potential to improve our understanding and quantification of brain pathologies (Wang et al., 2017).

Recently, many different strategies have been proposed to predict functional brain connectivity from structural connectivity. These range from sophisticated biophysical models (Galán, 2008; Honey et al., 2009; Deco et al., 2011; Deco and Jirsa, 2012; Messé et al., 2015a,b) to simpler approaches that do not attempt direct physical modeling (Deligianni et al., 2013; Abdelnour et al., 2014, 2018; Meier et al., 2016; Saggio et al., 2016; Liang and Wang, 2017; Becker et al., 2018; Nozari et al., 2020; Benkarim et al., 2021)'. A subset of these models specifically rely on the eigendecomposition of the structural connectivity matrix or its Laplacian (Robinson et al., 2016; Preti and Van De Ville, 2019) with promising applications to pathology (Raj et al., 2012; Abdelnour et al., 2021). Focusing on this type of approach, our group recently showed that these mappings can be unified under a general formulation (Deslauriers-Gauthier et al., 2020), highlighting the relation between series expansion and eigenmode approaches (Tewarie et al., 2020). This allowed us to highlight the link between the different models and directly compare them. While some models achieved high prediction accuracy, we found that all current model based on the eigendecomposition were outperformed by a simple group average. A possible explanation for this result is that the variability of the structural and functional connectivity is not sufficient in healthy subjects to be captured by the models (Deslauriers-Gauthier et al., 2020). A second possibility, the focus of this work, is that the tools used to train the models and evaluate their performance, i.e., the Frobenius norm and Pearson correlation, are not adapted to the space of functional connectivity matrices. By adapting the notion of distance to the manifold of functional connectivity matrices, it may be possible to improve the performance of existing structure–function mappings.

In this work, we revisit existing structure–function mapping models based on eigenmodes from a Riemannian perspective. Instead of computing distances between functional matrices in the vector space of symmetric matrices, we instead propose to use an affine invariant metric inducing a Riemannian distance. In addition to altering the objective function used to identify the parameters of the different mappings, this also changes the benchmark used to evaluate them. In our previous work (Deslauriers-Gauthier et al., 2020), we showed that a simple group average, which was used as a reference model, outperformed all eigenmode models. However, the use of a Riemannian distance for positive define matrices leads to a unique and well defined Fréchet mean (You and Park, 2021) which we propose as a more appropriate reference over the direct mean used previously. With this change of perspective, we demonstrate the importance of mapping function from structure under the Riemannian manifold, we improve the prediction of functional connectivity from structural connectivity, and we overcome the performance of the group average and the so–called glass ceiling on the performance of mappings based on eigenmodes.

2. Materials and Methods

2.1. Structure–Function Mapping

The objective of structure–function mapping is to identify a function f that minimizes the functional

L(f)=1Kk=1Kd2(f(Sk),Fk)    (1)

where SkM×M and FkM×M are the structural and functional matrices of the kth subject and M is the number of considered brain regions. The function f:ℝM × M → ℝM × M is the so called structure–function mapping whose parameters are optimized by minimizing Equation (1). Note that f could be defined indendently for each subject, but we focus here on the situation where a single mapping is optimized across all subjects. The function d(A, B) measures the distance between matrices A and B. If the Frobenius norm is used, Equation (1) reduces to the usual formulation of the problem which is to minimize

L(f)=1Kk=1K||f(Sk)-Fk||F2.    (2)

Once the parameters of the mapping f minimizing Equation (1) have been identified, its performance can be assessed on a new set of subjects by evaluating

E(f)=1K-Kk=K+1Ndt2(f(Sk),Fk)    (3)

where dt is the testing function. Interestingly, while the Frobenius norm or Euclidean distance is sometimes used, a more common choice is to use the Pearson correlation (Becker et al., 2018). In this case, the parameters of the mapping identified may be suboptimal as the optimization criterion differs for the evaluation criterion. However, optimizing the correlation directly is challenging as it is scale invariant and thus does not have a unique solution. In addition, the correlation is not a metric as it does not satisfy the triangle inequality which may lead to unintuitive comparison results. A metric can be calculated from the correlation via its relation to the cosine similarity, but there is no intrinsic justification for the use of angular distance between functional matrices. When the Frobenius norm is used for both optimization and evaluation, the model parameters are optimal, but the notion of distance is not adapted to the manifold of functional connectivity matrices as we discuss in the following section.

2.2. Riemannian Structure–Function Mapping

In the context of structure–function mapping, the functional connectivity matrices are typically computed as the Pearson correlation between the resting state BOLD fMRI signal of different brain areas. Let D ∈ ℝM × T be the data matrix where each row contains the demeaned and normalized time series of a brain area. The functional matrix can be written as F = DDT whose singular value decomposition is UΣ2UT. Functional connectivity matrices are therefore symmetric positive semidefinite. In practice, because the time series of a brain area is unlikely to be expressed exactly as a linear combination of the other areas, they are in fact symmetric positive definite (SPD). The space of SPD matrices is a convex smooth manifold of the Euclidean space, a convex half–cone in the vector space of symmetric matrices. Because of this conical structure, it is natural to consider Riemannian metrics on the space of SPD matrices. In particular, affine–invariant metrics have been suggested and induce a Riemannian distance

d2(A,B)=Tr(log2(A-1/2BA-1/2))    (4)

where log is the matrix logarithm and Tr is the trace. These notions have been used in the context of brain connectivity to classify connectivity matrices (Dodero et al., 2015) and regression (Wong et al., 2018). The Riemannian nature of the problem has also been observed by Benkarim et al. (2021), although they rely on the Euclidean distance for optimization and the correlation for evaluation. You and Park (2021) also re–visited Riemannian geometry for functional connectivity matrices, but they did not consider structure–function mapping specifically. Using a Riemannian distance over the Euclidean distance or correlation has several advantages. First, the distance is not computed elementwise but rather considers the general structure of the matrix (You and Park, 2021). This is consistent with the notion that the structure–function relationship is not one to one (Mišić et al., 2016) and instead operates on a network level (Atasoy et al., 2016; Xie et al., 2021). Second, matrices with non-positive or infinite eigenvalues are both infinitely far from any SPD matrices (Lenglet et al., 2006; Pennec et al., 2006), capturing the notion that they do not correspond to valid functional connectivity matrices. Finally, this space is smooth and allows the computation of a unique geodesic, and thus a unique minimal distance, between any pair of functional connectivity matrices. We thus propose to revisit the structure–function mappings based on eigenmodes using the Riemannian distance of Equation (4).

2.3. Reference Model

In our previous work (Deslauriers-Gauthier et al., 2020), we argued that the performance of structure–function mappings should be compared to a suitable reference model. A natural benchmark is the Fréchet mean F̄ which minimizes kd2(F̄,Fk) and can thus be seen as a zeroth order mapping, one that does not make use of structural information. When the Euclidean norm is used to estimate model parameters, it can be computed as the average of the functional matrices of the training set F̄E=1Kk=1KFk. However, when the distance is given by Equation (4), the Fréchet mean cannot be computed directly. It must instead be estimated using the following iterative scheme

F̄t+1=F̄t1/2exp(1Kk=1Klog(F̄t-1/2FkF̄t-1/2))F̄t1/2    (5)

which coverages to the unique value F̄R minimizing kd2(F̄,Fk) (Pennec et al., 2006). We use the definition above as a the reference model for our proposed methodology.

2.4. Optimization Procedure

As was shown in our previous work (Deslauriers-Gauthier et al., 2020), all existing mappings based on eigenmodes can be written as

f(S)=n=0N-1g(λn)h(un)+C    (6)

where C ∈ ℝN × N is a constant symmetric matrix and where λn ∈ ℝ and unN are the nth eigenvalue and eigenvector of S, respectively. For example, to reproduce the spectral mapping of Becker et al. (2018), we can define

g(λ)=m=0Papλpand h(un)=RununTRT    (7)

where am ∈ ℝ and R is a rotation matrix. The formulation in Equation (6) allows us to devise a general optimization scheme for all existing mappings. The objective is to minimize Equation (1) where the distance is given by Equation (4), that is

L(f)=1Kk=1KTr(log2(Fk-1/2f(Sk)Fk-1/2))    (8)

By substituting the general form in Equation (6), we obtain

L(f)=1Kk=1KTr(log2(Fk-1/2(n=0N-1g(λn)h(un)+C)Fk-1/2))        =1Kk=1KTr(log2(n=0N-1g(λn)hk(un)+Ck))

where

hk=Fk-1/2h(un)Fk-1/2and Ck=Fk-1/2CFk-1/2.    (9)

To minimize this functional, we iteratively minimize over the parameters of each term g, h′, and C′ while keeping the others fixed in a manner similar to the coordinate descent algorithm. For the eigenmode weights g, the parameters were optimized using the Broyden–Fletcher–Goldfarb–Shanno (BFGS) algorithm via scipy (Virtanen et al., 2020). For h′ and C′, the parameters were optimized using Pymanopt1 to ensure the parameters remained on the manifold of rotation and symmetric matrices, respectively. In practice, the performance of all models converged after approximately three iterations.

2.5. Evaluation

We evaluated the performance of the mappings proposed by Abdelnour et al. (2014, 2018) and spectral mappings (Meier et al., 2016; Becker et al., 2018; Deslauriers-Gauthier et al., 2020) which were all implemented using the general model in Equation (6). The specific form of the mappings is provided in Table 1. For the specific cases of the diffusion mappings, the Laplacian of the structural matrix, defined as L = ID−1/2SD−1/2 where D is the diagonal degree matrix, was used as input. For models relying on spectral mapping, the maximum order was set to M = 6. For each mapping, the parameters were learned using a leave-one-out cross-validation strategy, meaning K′ = 100 and K = 99 in Equation (3). The performance of each mapping was then evaluated by measuring the distance between the predicted and observed functional connectivity matrices using Equation (4) on the left out subject. Finally, the performance of Fréchet mean in Equation (5) was evaluated using the same strategy and used as a benchmark.

TABLE 1
www.frontiersin.org

Table 1. Definition of the structure–function mappings based on eigenmodes and their number of degrees of freedom.

2.6. Data and Data Processing

The data was provided by the Human Connectome Project2 (HCP) and our processing started from the minimally processed data (Glasser et al., 2013) for both diffusion and functional MRI data. The rest of pipeline is very similar to the one used in our previous work (Deslauriers-Gauthier et al., 2020), which we describe here again for completeness. For 100 unrelated subjects provided by the HCP, the brain was extracted using FSL bet (Smith, 2002) and the white matter, gray matter, and cerebrospinal fluid were segmented using FSL fast (Zhang et al., 2001). Using the FreeSurfer3 segmentation (Fischl, 2012) and the diffusion weighted images, a single fiber response function was computed for the white matter, gray matter, and cerebrospinal fluid using MRtrix34 (Tournier et al., 2019). These response function were then used to compute a fiber orientation distribution function for each voxel using contrained spherical deconvolution (Tournier et al., 2007). Five million streamlines were generated using anatomically constrained probabilistic tractography using a step of 0.3 mm, a maximum length of 400 mm, and backtracking (Smith et al., 2012). In an effort to make connectomes more quantitative and to reduce the impact of false positive connections (Maier-Hein et al., 2017) by re-establishing the biological interpretability of streamline-based structural connections, the tractograms were filtered using SIFT2 (Smith et al., 2015) assigning a weight to each streamline representing its cross-sectional area. The cortical surface extracted with FreeSurfer was parcelated using the Schaefer atlas (Schaefer et al., 2018) into N = 200 regions. Given the parcellation and the streamlines, the structural connectome was built by summing the SIFT2 weights of streamlines connecting two cortical regions. The connectomes were symmetrized by summing the (i, j) and (j, i) entries of the connectome. Finally, the structural connectomes were normalized by dividing by the sum of the off diagonal entries and the diagonal set to zeros. While some authors removed weak connections we followed the recommendations of Civier et al. (2019) and omitted pruning thus keeping weak connections. The rational is that weak connections have little impact on values derived from connectivity matrices. It is therefore preferable to simplify the processing pipeline and reduce the number of arbitrary parameters by omitting this step.

The resting state functional MRI data provided by the HCP offers two sessions each containing two acquisitions, one with left to right encoding and one with right to left encoding. Here, we made use of the first (REST1) session and preprocessed the two acquisitions as follows. The time series were filtered using a butterworth bandpass filter with critical frequencies 0.01 and 0.1 Hz (van den Heuvel and Hulshoff Pol, 2010). The first 20 volumes of each acquisition were dropped and the remaining volumes linearly detrended (Caballero-Gaudes and Reynolds, 2017). The movement parameters and their derivatives were regressed out and the data was motion scrubbed (Power et al., 2012). After this preprocessing, we concatenated both acquisitions to produce a single dataset per subject. The mean signal for each parcel was computed and used to build the functional connectivity matrix using Pearson correlation. To preserve the SPD nature of the functional connectomes, the negative entries of the matrices were preserved, in contrast to our previous work.

3. Results and Discussion

3.1. Riemannian vs. Euclidean Distances

First, we investigated the impact of measuring distances between functional matrices using the Euclidean or Riemannian metrics. To this end, we measured the distance between the functional matrices for every pair of subjects using both metrics. In addition, for each subject, we identified the other subject having the closest functional matrix, i.e., the nearest neighbor, again for both distances. The resulting distance and nearest neighbor matrices are illustrated in Figure 1. The Riemannian distances have a mean of 31.5∓2.3 with a range of 26.0 to 40.4 and the Euclidean distribution has a mean of 51.4∓9.3 with a range of 30.3 to 90.1. This change in distribution indicates that using the Riemannian distance significantly changes the relationship between functional matrices. To further confirm this observation, we note that the nearest neighbors for each distances do not agree. Indeed, only nine subjects had the same nearest neighbors using both distances, while 91 subjects differed. From these results, we conclude that the notion of similarity between subjects is significantly altered by the change in metric. Consequently, the relationship between structure and function will also be affect and reflected in the structure–function mappings.

FIGURE 1
www.frontiersin.org

Figure 1. Histogram of the distances between the functional matrices of every subject using the Euclidean and Riemannian metrics. The nearest neighbor, minimizing each distance, are also presented. Teal and blue elements correspond to Euclidean and Riemannian nearest neighbors respectively and yellow indicates agreement between the two metrics.

To investigate whether these changes improved the structure–function relation, we additionally computed the distance between structural matrices of every subject. Figure 2 illustrates the functional distance as a function of the structural distance for every subject. Both distances where converted to z-scores before plotting. We can observe that subjects that have a short structural distance are more likely to have a short functional distance when measured with the Riemannian distance rather than the Euclidean distance. To quantify this observation, we computed the coefficient of determination R2 for each distance and obtained 0.016 and 0.135 for the Euclidean and Riemannian distances, respectively. In other words, subjects that have a similar structural matrix also have a similar functional matrix, but only when distances are measured using the Riemannian metric. Note that these results are independent of any mapping strategy and instead provide evidence in favor of Riemannian metric on a more fundamental level. It is indeed expected that subjects with a similar structure should also have a similar function, but the R2 score near zero obtained by the Euclidean metrics would instead indicate a decorrelation between these two views of the brain.

FIGURE 2
www.frontiersin.org

Figure 2. Distances between the structural matrices as a function of the distance between functional matrices for every subject. Both distances where converted to z-scores before plotting. On the (Left) plot, the distance between functional matrices is measured using the Euclidean metric and on the (Right) with the Riemannian metric. The coefficients of determination R2 are 0.016 and 0.135 for the Euclidean and Riemannian distances, respectively.

To leverage this correlation between the structural and functional distances, we devised an additional mapping and included it in the following experiments. In this new mapping, the predicted functional connectivity matrix is given the Riemannian mean of the P training subjects whose structural connectivity matrices are closest to that the of the new subject, i.e., the P nearest structural neighbors. The rationale is that subjects having a similar structure will also have a similar function and averaging only the structural nearest neighbors will lead to an improve functional estimation. In the following section we refer to this mapping as the Riemannian nearest neighbors mean and empirically choose P = 50.

3.2. Structure–Function Mapping

Figure 3 illustrates the performance of the different mappings optimized using the procedure in Section 2.4. Notably, the mapping devised in Deslauriers-Gauthier et al. (2020) by adding a constant to the spectral mapping of Becker et al. (2018) obtains a mean squared distance (MSD) of 449.9 and outperforms the Riemannian mean which obtained a distance of 450.5. While the performance gain is admittedly small, we found the differences in errors to be statistically significant (p < 0.001). Another notable change is the inversion of the performance between the diffusion model of Abdelnour et al. (2018) and the spectral model of Meier et al. (2016) with respect to our previous study (Deslauriers-Gauthier et al., 2020). Using the Euclidean metric, the diffusion model was outperformed by the spectral model, a result that was justified by the parsimonious nature of the diffusion model. Using the Riemannian metric, the diffusion model obtained an MSD of 558.7 whereas the spectral model obtained 578.4. Except for this change, the performance relationship between previously proposed mapping strategies followed the same trends as in our previous work. In general, mappings with more degrees of freedom (see Table 1) provided better predictions, but at the cost of increased complexity. Finally, the Riemannian nearest neighbors mean mapping obtained the best performance overall with a MSD of 448.9. This improvement in performance illustrates that the correlation between structural and functional distances, only observed using the Riemannian metric, can be leveraged to improve our prediction of function from structure.

FIGURE 3
www.frontiersin.org

Figure 3. Mean squared Riemannian distance of five previously proposed mapping based on eigenmodes, the reference mapping (Riemannian mean), and the Riemannian nearest neighbors mean (P = 50).

3.3. Limitations

One of the fundamental assumptions in structure–function mapping is that changes in the structural connectivity of subjects affects their functional connectivity. However, this study considered healthy young adults, a very homogeneous cohort with similar structural connectivity. As such, the structural differences between subjects are subtle and may contribute to the high prediction accuracy of the mean. A second limitation is the use of existing mappings that were not explicitly designed with the tools of Riemannian geometry. This may also contribute to the high performance of the mean, which is optimal in a Riemannian sense. An interesting avenue of research is the design of mappings that leverage the topology of the space of SPD matrices to improve the predictions performance. A simple proof-of-concept for this approach is the Riemannian nearest neighbors mean which obtained the lowest MSD, and thus the best functional prediction.

It is also important to note that the impact of the different processing steps on our ability to predict function from structure in a Riemannian setting was not evaluated. Indeed, both the functional and structural connectivity pipelines require the selection of parameters which will modify the estimated connectivity. For example, structural matrices can be constructed by counting streamlines or by using a weighting strategy such as COMMIT (Schiavi et al., 2020) or SIFT2 as was done here. To investigate the importance of this specific choice, we reproduced the experiments using streamline counting structural connectomes. The results obtained (available in Supplementary Materials), are close to those produced by SIFT2 connectomes and do not alter the conclusions. However, this may not be the case for every step of the pipelines and an exhaustive search may lead to further improvement of the predictions.

Finally, the structural connectivity matrices of this study only considered white matter cortico-cortical connections. Other connections, for example those including subcortical nodes and superficial white matter fibers (Reveley et al., 2015), are not encoded in the matrices. As such, it is unlikely that functional connectivity can be completely predicted by the structural matrices. Adding these missing connections may very well improve the performance of structure-function mappings, in addition to providing a more complete view of structural connectivity.

4. Conclusions

In this work, we investigated the implications of using a distance based on an affine invariant Riemannian metric in the context of structure–function mapping. We argued that functional matrices, if they are not thresholded, are symmetric positive definite and therefore live in a convex half–cone in the vector space of symmetric matrices. By using a distance appropriate to this manifold, we showed that a previously proposed mapping is able to outperform the group average, a result which was not obtained with an Euclidean metric. While the improvement in performance with respect to the mean was moderate, it should be observed that the mappings were not specifically designed with this manifold in mind and could thus be improved. In addition, while we chose to focus on mappings based on eigenmodes, our findings may also have interesting extensions to deep learning approaches for structure-function mapping (Ji et al., 2021). Indeed, the already promising performance of these approaches may be further enhanced by the use of Riemannian geometry. In conclusion, our results show that using the correct distance for SPD matrices improved the performance of existing structure–function mappings and may lead to better prediction with more specific models.

Data Availability Statement

Publicly available datasets were analyzed in this study. This data can be found here: Human Connectome Project (HCP) Young Adult, https://db.humanconnectome.org.

Author Contributions

SD-G: conceptualization, methodology, software, investigation, writing—original draft, and visualization. MZ: conceptualization, methodology, and writing—review and editing. HL: methodology and writing—review and editing. RD: conceptualization, methodology, writing—review and editing, supervision, and funding acquisition. All authors contributed to the article and approved the submitted version.

Funding

This work has received funding from the European Research Council (ERC) under the European Union's Horizon 2020 research and innovation program (ERC Advanced Grant agreement no. 694665: CoBCoM—Computational Brain Connectivity Mapping). This work has been supported by the French government, through the 3IA Côte d'Azur Investments in the Future project managed by the National Research Agency (ANR) with the reference number ANR-19-P3IA-0002. The authors are grateful to the OPAL infrastructure from Université Côte d'Azur for providing resources and support.

Conflict of Interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Publisher's Note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fnimg.2022.850266/full#supplementary-material

Footnotes

References

Abdelnour, F., Dayan, M., Devinsky, O., Thesen, T., and Raj, A. (2018). Functional brain connectivity is predictable from anatomic network's Laplacian eigen-structure. NeuroImage 172, 728–739. doi: 10.1016/j.neuroimage.2018.02.016

PubMed Abstract | CrossRef Full Text | Google Scholar

Abdelnour, F., Dayan, M., Devinsky, O., Thesen, T., and Raj, A. (2021). Algebraic relationship between the strutural network's laplacian and functional network's adjacency matrix is preserved in temporal lobe epilepsy subjects. NeuroImage 228:117705. doi: 10.1016/j.neuroimage.2020.117705

PubMed Abstract | CrossRef Full Text | Google Scholar

Abdelnour, F., Voss, H. U., and Raj, A. (2014). Network diffusion accurately models the relationship between structural and functional brain connectivity networks. NeuroImage 90, 335–347. doi: 10.1016/j.neuroimage.2013.12.039

PubMed Abstract | CrossRef Full Text | Google Scholar

Atasoy, S., Donnelly, I., and Pearson, J. (2016). Human brain networks function in connectome specific hamonic waves. Nat. Commun. 7:10340. doi: 10.1038/ncomms10340

PubMed Abstract | CrossRef Full Text | Google Scholar

Becker, C. O., Pequito, S., Pappas, G. J., Miller, M. B., Grafton, S. T., Bassett, D. S., et al. (2018). Spectral mapping of brain functional connectivity from diffusion imaging. Sci. Rep. 8:1411. doi: 10.1038/s41598-017-18769-x

PubMed Abstract | CrossRef Full Text | Google Scholar

Benkarim, O., Paquola, C., Park, B.-Y., Royer, J., Rodríguez-Cruces, R., de Wael, R. V., et al. (2021). A riemannian approach to predicting brain function from the structural connectome. bioRxiv [Preprint]. doi: 10.1101/2021.10.27.465906

CrossRef Full Text | Google Scholar

Caballero-Gaudes, C., Reynolds, R. C. (2017). Methods for cleaning the BOLD fMRI signal. NeuroImage 154, 128–149. doi: 10.1016/j.neuroimage.2016.12.018

PubMed Abstract | CrossRef Full Text | Google Scholar

Civier, O., Smith, R. E., Yeh, C. H., Connelly, A., and Calamante, F. (2019). Is removal of weak connections necessary for graph-theoretical analysis of dense weighted structural connectomes from diffusion MRI? NeuroImage 194, 68–81. doi: 10.1016/j.neuroimage.2019.02.039

PubMed Abstract | CrossRef Full Text | Google Scholar

Deco, G., Jirsa, V. K. (2012). Ongoing cortical activity at rest: criticality, multistability, and ghost attractors. J. Neurosci. 32, 3366–3375. doi: 10.1523/JNEUROSCI.2523-11.2012

PubMed Abstract | CrossRef Full Text | Google Scholar

Deco, G., Jirsa, V. K., McIntosh, A. R. (2011). Emerging concepts for the dynamical organization of resting-state activity in the brain. Nat. Rev. Neurosci. 12, 43–56. doi: 10.1038/nrn2961

PubMed Abstract | CrossRef Full Text | Google Scholar

Deligianni, F., Varoquaux, G., Thirion, B., Sharp, D. J., Ledig, C., Leech, R., et al. (2013). A framework for inter-subject prediction of functional connectivity from structural networks. IEEE Trans. Med. Imaging 32, 2200–2214. doi: 10.1109/TMI.2013.2276916

PubMed Abstract | CrossRef Full Text | Google Scholar

Deslauriers-Gauthier, S., Zucchelli, M., Frigo, M., and Deriche, R. (2020). A unified framework for multimodal structure-function mapping based on eigenmodes. Med. Image Anal. 66:101799. doi: 10.1016/j.media.2020.101799

PubMed Abstract | CrossRef Full Text | Google Scholar

Dodero, L., Minh, H. Q., Biagio, M. S. V. M., and Sona, D. (2015). “Kernel-based classification for brain connectivity graphs on the Riemannian manifold of positive definite matrices,” in 2015 IEEE 12th International Symposium on Biomedical Imaging (ISBI) (Brooklyn, NY), 42–45. doi: 10.1109/ISBI.2015.7163812

CrossRef Full Text | Google Scholar

Fischl, B.. (2012). FreeSurfer. NeuroImage 62, 774–781. doi: 10.1016/j.neuroimage.2012.01.021

PubMed Abstract | CrossRef Full Text | Google Scholar

Galán, R. F.. (2008). On how network architecture determines the dominant patterns of spontaneous neural activity. PLoS ONE 3:e2148. doi: 10.1371/annotation/2c9bfbcb-6b96-4d77-bfe3-10c5988150b8

PubMed Abstract | CrossRef Full Text | Google Scholar

Glasser, M. F., Sotiropoulos, S. N., Wilson, J. A., Coalson, T. S., Fischl, B., Andersson, J. L., et al. (2013). The minimal preprocessing pipelines for the Human Connectome Project. NeuroImage 80, 105–124. doi: 10.1016/j.neuroimage.2013.04.127

PubMed Abstract | CrossRef Full Text | Google Scholar

Honey, C. J., Sporns, O., Cammoun, L., Gigandet, X., Thiran, J. P., Meuli, R., et al. (2009). Predicting human resting-state functional connectivity from structural connectivity. Proc. Natl. Acad. Sci. U.S.A. 106, 2035–2040. doi: 10.1073/pnas.0811168106

PubMed Abstract | CrossRef Full Text | Google Scholar

Ji, Y., Deslauriers-Gauthier, S., and Deriche, R. (2021). “Structure-function mapping via graph neural networks,” in Machine Learning in Clinical Neuroimaging. MLCN 2021. Lecture Notes in Computer Science, Vol. 13001, eds A. Abdulkadir, S. M. Kia, M. Habes, V. Kumar, J. M. Rondina, C. Tax, and T. Wolfers (Cham: Springer). doi: 10.1007/978-3-030-87586-2_14

PubMed Abstract | CrossRef Full Text | Google Scholar

Lenglet, C., Rousson, M., Deriche, R., and Faugeras, O. (2006). Statistics on the manifold of multivariate normal distributions: THEORY and applications to diffusion tensor MRI processing. J. Math. Imag. Vis. 25, 423–444. doi: 10.1007/s10851-006-6897-z

CrossRef Full Text | Google Scholar

Liang, H., and Wang, H. (2017). Structure-function network mapping and its assessment via persistent homology. PLoS Comput. Biol. 13:e1005325. doi: 10.1371/journal.pcbi.1005325

PubMed Abstract | CrossRef Full Text | Google Scholar

Maier-Hein, K. H., Neher, P. F., Houde, C., Côté, M.-A., Garyfallidis, E., Zhong, J., et al. (2017). The challenge of mapping the human connectome based on diffusion tractography. Nat. Commun. 8:1349. doi: 10.1038/s41467-017-01285-x

PubMed Abstract | CrossRef Full Text | Google Scholar

Meier, J., Tewarie, P., Hillebrand, A., Douw, L., van Dijsk, B. W., Stufflebeam, S. M., et al. (2016). A mapping between structural and functional brain networks. Brain Connect. 6, 298–311. doi: 10.1089/brain.2015.0408

PubMed Abstract | CrossRef Full Text | Google Scholar

Messé, A., Benali, H., and Marrelec, G. (2015a). Relating structural and functional connectivity in MRI: a simple model for a complex brain. IEEE Trans. Med. Imaging 34, 27–37. doi: 10.1109/TMI.2014.2341732

PubMed Abstract | CrossRef Full Text | Google Scholar

Messé, A., Hutt, M. T., Konig, P., Hilgetag, C. C. (2015b). A closer look at the apparent correlation of structural and functional connectivity in excitable neural networks. Sci. Rep. 5:7870. doi: 10.1038/srep07870

PubMed Abstract | CrossRef Full Text | Google Scholar

Mišić, B., Betzel, R. F., de Reus, M. A., van den Heuvel, M. P., Berman, M. G., McIntosh, A. R., et al. (2016). Network-level structure-function relationships in human neocortex. Cereb. Cortex 26, 3285–3296. doi: 10.1093/cercor/bhw089

PubMed Abstract | CrossRef Full Text | Google Scholar

Nozari, E., Bertolero, M. A., Stiso, J., Caciagli, L., Cornblath, E. J., He, X., et al. (2020). Is the brain macroscopically linear? A system identification of resting state dynamics. arXiv. doi: 10.1101/2020.12.21.423856

CrossRef Full Text | Google Scholar

Pennec, X., Fillard, P., and Ayache, N. (2006). A Riemannian framework for tensor computing. Int. J. Comput. Vis. 66, 41–66. doi: 10.1007/s11263-005-3222-z

CrossRef Full Text | Google Scholar

Power, J. D., Barnes, K. A., Snyder, A. Z., Schlaggar, B. L., Petersen, S. E. (2012). Spurious but systematic correlations in functional connectivity MRI networks arise from subject motion. NeuroImage 59, 2142–2154. doi: 10.1016/j.neuroimage.2011.10.018

PubMed Abstract | CrossRef Full Text | Google Scholar

Preti, M. G., Van De Ville, D. (2019). Decoupling of brain function from structure reveals regional behavioral specialization in humans. Nat. Commun. 10:4747. doi: 10.1038/s41467-019-12765-7

PubMed Abstract | CrossRef Full Text | Google Scholar

Raj, A., Kuceyeski, A., and Weiner, M. (2012). A network diffusion model of disease progression in dementia. Neuron 73, 1204–1215. doi: 10.1016/j.neuron.2011.12.040

PubMed Abstract | CrossRef Full Text | Google Scholar

Reveley, C., Seth, A. K., Pierpaoli, C., Silva, A. C., Yu, D., Saunders, R. C., et al. (2015). Superficial white matter fiber systems impede detection of long-range cortical connections in diffusion MR tractography. Proc. Natl. Acad. Sci. U.S.A. 112, E2820–E2828. doi: 10.1073/pnas.1418198112

PubMed Abstract | CrossRef Full Text | Google Scholar

Robinson, P. A., Aquino, K. M., Griffiths, J., Sarkar, S., Mehta-Pandejee, G. (2016). Eigenmodes of brain activity: neural field theory predictions and comparison with experiment. NeuroImage 142, 79–98. doi: 10.1016/j.neuroimage.2016.04.050

PubMed Abstract | CrossRef Full Text | Google Scholar

Saggio, M. L., Ritter, P., Jirsa, V. K. (2016). Analytical operations relate strutural and functional connectivity in the brain. PLoS ONE 11:e157292. doi: 10.1371/journal.pone.0157292

PubMed Abstract | CrossRef Full Text | Google Scholar

Schaefer, A., Kong, R., Gordon, E. M., Laumann, T. O., Zuo, X.-N., Holmes, A. J., et al. (2018). Local-global parcellation of the human cerebral cortex from intrinsic functional connectivity MRI. Cereb. Cortex 28, 3095–3114. doi: 10.1093/cercor/bhx179

PubMed Abstract | CrossRef Full Text | Google Scholar

Schiavi, S., Ocampo-Pineda, M., Barakovic, M., Petit, L., Descoteaux, M., Thiran, J.-P., et al. (2020). A new method for accurate in vivo mapping of human brain connections using microstructural and anatomical information. Sci. Adv. 6:eaba8245. doi: 10.1126/sciadv.aba8245

PubMed Abstract | CrossRef Full Text | Google Scholar

Smith, R. E., Tournier, J. D., Calamante, F., and Connelly, A. (2012). Anatomically-constrained tractography: improved diffusion MRI streamlines tractography through effective use of anatomical information. NeuroImage 62, 1924–1938. doi: 10.1016/j.neuroimage.2012.06.005

PubMed Abstract | CrossRef Full Text | Google Scholar

Smith, R. E., Tournier, J. D., Calamante, F., and Connelly, A. (2015). SIFT2: Enabling dense quantitative assessment of brain white matter connectivity using streamlines tractography. NeuroImage 119, 338–351. doi: 10.1016/j.neuroimage.2015.06.092

PubMed Abstract | CrossRef Full Text | Google Scholar

Smith, S. M.. (2002). Fast robust automated brain extraction. Hum. Brain Mapp. 17, 143–155. doi: 10.1002/hbm.10062

PubMed Abstract | CrossRef Full Text | Google Scholar

Sporns, O., Tonomi, G., Kötter, R. (2005). The human connectome: a structural description of the human brain. PLoS Comput. Biol. 1:e10042. doi: 10.1371/journal.pcbi.0010042

PubMed Abstract | CrossRef Full Text | Google Scholar

Tewarie, P., Prasse, B., Meier, J. M., Santos, F. A., Douw, L., Shoonheim, M. M. (2020). Mapping functional brain networks from the structural connectome: relating the series expansion and eigenmode approaches. NeuroImage 216:116805. doi: 10.1016/j.neuroimage.2020.116805

PubMed Abstract | CrossRef Full Text | Google Scholar

Tournier, J. D., Calamante, F., and Connelly, A. (2007). Robust determination of the fibre orientation distribution in diffusion MRI: Non-negativity constrained super-resolved spherical deconvolution. NeuroImage 35, 1459–1472. doi: 10.1016/j.neuroimage.2007.02.016

PubMed Abstract | CrossRef Full Text | Google Scholar

Tournier, J. D., Smith, R., Raffelt, D., Tabbara, R., Dhollander, T., Pietsch, M., et al. (2019). MRtrix3: A fast, flexible and open software framework for medical image processing and visualisation. NeuroImage 202:116137. doi: 10.1016/j.neuroimage.2019.116137

PubMed Abstract | CrossRef Full Text | Google Scholar

van den Heuvel, M. P., Hulshoff Pol, H. E. (2010). Exploring the brain network: a review on resting-state fMRI functional connectivity. Eur. Neuropsychopharmacol. 20, 519–534. doi: 10.1016/j.euroneuro.2010.03.008

PubMed Abstract | CrossRef Full Text | Google Scholar

Virtanen, P., Gommers, R., Oliphant, T. E., Haberland, M., Reddy, T., Cournapeau, D., et al. (2020). SciPy 1.0: fundamental algorithms for scientific computing in Python. Nat. Methods 17, 261–272. doi: 10.1038/s41592-020-0772-5

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, M. B., Owen, J. P., Mukherjee, P., and Raj, A. (2017). Brain network eigenmodes provide a robust and compact representation of the structural connectome in health and disease. PLOS Comput. Biol. 13:e1005550. doi: 10.1371/journal.pcbi.1005550

PubMed Abstract | CrossRef Full Text | Google Scholar

Wong, E., Anderson, J. S., Zielinski, B. A., Fletcher, P. T. (2018). “Riemannian regression and classification models of brain networks applied to autism,” in Connectomics in NeuroImaging, eds G. Wu, I. Rekik, M. Schirmer, A. Chung, and B. Munsell (Cham: Springer), 78–87. doi: 10.1007/978-3-030-00755-3_9

PubMed Abstract | CrossRef Full Text | Google Scholar

Xie, X., Cai, C., Damasceno, P. F., Nagarajan, S. S., and Raj, A. (2021). Emergence of canonical functional networks form the structral connectome. NeuroImage 237:118190. doi: 10.1016/j.neuroimage.2021.118190

PubMed Abstract | CrossRef Full Text | Google Scholar

You, K., Park, H.-J. (2021). Re-visiting Riemannian geometry of symmetric positive definite matrices for the analysis of functional connectivity. NeuroImage 225:117464. doi: 10.1016/j.neuroimage.2020.117464

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhang, Y., Brady, M., and Smith, S. (2001). Segmentation of brain MR images through a hidden Markov random field model and the expectation-maximization algorithm. IEEE Trans. Med. Imag. 20, 45–57. doi: 10.1109/42.906424

PubMed Abstract | CrossRef Full Text | Google Scholar

Keywords: brain structure-function mapping, functional connectivity, structural connectivity, eigenvalue decomposition, Riemannian distance

Citation: Deslauriers-Gauthier S, Zucchelli M, Laghrissi H and Deriche R (2022) A Riemannian Revisiting of Structure–Function Mapping Based on Eigenmodes. Front. Neuroimaging 1:850266. doi: 10.3389/fnimg.2022.850266

Received: 07 January 2022; Accepted: 23 March 2022;
Published: 25 May 2022.

Edited by:

David Wack, University at Buffalo, United States

Reviewed by:

Alonso Ramirez-Manzanares, Centro de Investigación en Matemáticas, Mexico
Tuo Zhang, Northwestern Polytechnical University, China

Copyright © 2022 Deslauriers-Gauthier, Zucchelli, Laghrissi and Deriche. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Samuel Deslauriers-Gauthier, samuel.deslauriers-gauthier@inria.fr

Disclaimer: All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.