## Correlation Tensor Magnetic Resonance Imaging

Rafael Neto Henriques1, Sune N Jespersen2,3 and Noam Shemesh1*

1Champalimaud Research, Champalimaud Centre for the Unknown, Lisbon, Portugal 2Center of Functionally Integrative Neuroscience (CFIN) and MINDLab, Department of Clinical Medicine, Aarhus University, Aarhus, Denmark. 3Department of Physics and Astronomy, Aarhus University, Aarhus, Denmark

*Corresponding author:

Dr. Noam Shemesh, Champalimaud Research, Champalimaud Centre for the Unknown, Av.

Brasilia 1400-038, Lisbon, Portugal

E-mail: noam.shemesh@neuro.fchampalimaud.org;

Phone number: +351 210 480 000 ext. #4467

## Abstract

Diffusional  Kurtosis  Magnetic  Resonance  Imaging  (DKI)  quantifies  the  extent  of  nonGaussian water diffusion, which has been shown to be a sensitive biomarker for microstructure in health and disease. However, DKI is not specific to any microstructural property per se since kurtosis may emerge from several different sources. Q-space trajectory encoding schemes have been proposed for decoupling kurtosis arising from the variance of mean diffusivities (isotropic kurtosis)  from kurtosis  driven by  microscopic  anisotropy (anisotropic kurtosis). Still, these methods assume that the system is comprised of multiple Gaussian diffusion components with vanishing intra-compartmental kurtosis (associated with restricted diffusion). Here, we develop a more general framework for resolving the underlying kurtosis sources without relying on the multiple Gaussian diffusion approximation. We introduce Correlation Tensor MRI (CTI) -an approach harnessing the versatility of double diffusion encoding (DDE) and its sensitivity to displacement  correlation  tensors  capable  of  explicitly  decoupling  isotropic  and  anisotropic kurtosis  components  from  intra-compartmental  kurtosis  effects  arising  from  restricted  (and time-dependent)  diffusion.  Additionally,  we  show  that,  by  subtracting  these  isotropic  and anisotropic kurtosis components from the total diffusional kurtosis, CTI provides an index that is potentially sensitive to intra-compartmental kurtosis. The theoretical foundations of CTI, as well as the first proof-of-concept CTI ex vivo experiments in mouse brains at ultrahigh field of 16.4 T, are presented. We find that anisotropic and isotropic kurtosis can decouple microscopic anisotropy from substantial partial volume effects between tissue and free water. Our intracompartmental kurtosis index exhibited positive values in both white and grey matter tissues. Simulations  in  different  synthetic  microenvironments  show,  however,  that  our  current  CTI protocol for estimating intra-compartmental kurtosis is limited by higher order terms that were not taken into account in this study. CTI measurements were then extended to in vivo settings and used to map heathy rat brains at 9.4 T. These in vivo CTI results were found to be consistent with our ex vivo findings. Although future studies are still required to assess and mitigate the higher order effects on the intra-compartmental kurtosis index, our results show that CTI's more general estimates of anisotropic and isotropic kurtosis contributions are already ripe for future in vivo studies, which can have significant impact our understanding of the mechanisms underlying diffusion metrics extracted in health and disease.

## 1. Introduction

Sensing microstructural features of biological systems noninvasively is vital for understanding how large-scale biological systems evolve over time. Tissue microarchitecture is constantly remodelled,  whether  due  to  normal  processes  such  as  development,  learning  and  aging (Falangola et al., 2008; Moseley, 2002; Neil et al., 1998; Pfefferbaum et al., 2000), or due to abnormal processes such as disease progression or acute insults to the tissue (Cheung et al., 2012;  Fieremans  et  al.,  2013;  Moseley  et  al.,  1990a).  Often,  the  microstructural  changes precede functional outcomes: for example, spine density increases rapidly before learning has taken  place  (Xu  et  al.,  2009);  subtle  changes  in  cellular  density  and  structure  precede  the functional deficits incurred in neurological disorders (Hanisch and Kettenmann, 2007); and malignant  transformations  can  occur  well  before  tumours  can  be  detected  (Peinado  et  al., 2017). All these reinforce the need for accurate in vivo mapping of microstructural properties (Le Bihan and Johansen-Berg, 2012).

Diffusion MRI (dMRI (Le Bihan and Johansen-Berg, 2012)), mainly based on variants of  the  Single  Diffusion  Encoding  (SDE,  Shemesh  et  al.,  2016)  methodology  developed originally by Stejskal and Tanner (Stejskal and Tanner, 1965), has become the mainstay of contemporary non-invasive  microstructural imaging. Water molecules traverse  microscopic length scales on typical MR-relevant observation times at body temperature, and their diffusion properties are influenced by the presence of restricting boundaries, such as cell membranes and other subcellular structures. dMRI capitalizes on this 'endogenous sensor' by sensitizing the MRI  signal  towards  molecular  displacements  in  a  given  orientation,  thereby  enabling  the quantification  of  water  diffusion  properties  (Assaf  and  Cohen,  1998;  Jensen  et  al.,  2005; Moseley  et  al.,  1990a).  In  many  cases,  it  is  assumed  that  water  diffusion  can  be  fully characterized  by  a  single  apparent  diffusion  tensor  (Basser  et  al.,  1994).  Diffusion  Tensor

Imaging (DTI) can extract this apparent tensor from multiple diffusion-weighted measurements; the t ensor's trace has been shown to be sensitive, for example, towards early phases of acute stroke (Reith et al., 1995), while the tensor's orientation can be used to recover the absolute orientation of coherently aligned white matter tracts (Catani et al., 2002; Jones, 2008; Mori et al., 1999), which has been instrumental  to, e.g., surgical planning (Berman, 2009).

DTI  and  similar  methods  represent  the  dMRI  signal  as  being  sufficiently  well characterized by Gaussian diffusion (Basser, 1995; Dell'Acqua et al., 2007; Desco teaux et al., 2009; Tournier et al., 2007). Implicitly, this means that the (logarithm of the) diffusion signal is represented only up to first order in b-value, where the b-value represents the strength of diffusion weighting (Le Bihan et al., 1986; Le Bihan and Breton, 1985). However, it has been recognized  already  early  on  that  diffusion  in  biological  systems  is  generally  non-Gaussian (Assaf and Cohen, 1998; Mulkern et al., 1999; Sukstanskii and Yablonskiy, 2002; Yablonskiy et al., 2003), and that characterizing the non-Gaussian effects may provide deeper insights into tissue  microstructure.  To  capture  non-Gaussian  diffusion  effects  via  signal  representations, Jensen et al. developed diffusional kurtosis imaging (DKI) (Jensen et al., 2005). In DKI, the signal  is  expanded  using  cumulants up to second order in b-value, quantifying the leading deviation  from  Gaussian  diffusion  and  giving  rise  to  a  source  of  contrast  based  on  nonGaussian  properties  of  the  signal.  DKI  has  been  shown  to  be  more  sensitive  than  its  DTI counterpart towards quantifying microstructural changes related to, e.g. aging, development and disease (Cheung et al., 2012; Falangola et al., 2008; Fieremans et al., 2013; Gong et al., 2013; Helpern et al., 2011; Hui et al., 2012; Rudrapatna et al., 2014; Wang et al., 2011).

Despite  the  utility  of  DTI  and  DKI,  both  methods  conflate  mesoscopic  orientation dispersion and true microstructural properties (De Santis et al., 2014; Henriques et al., 2015; Jones et al.,  2013;  Szczepankiewicz et al., 2015).  For example, consider a system with  an

ensemble of microscopic components (or microenvironments) represented by diffusion tensors with identical trace and anisotropy, dispersed along a mesoscopic orientation distribution (Fig. 1A). The microscopic features in this system consist of microscopic anisotropy (the anisotropy of the microenvironment tensors in their own eigenframe), the distribution of diffusion tensor traces  (in  the  system  considered  above,  just  a  delta  function),  and  the  extent  of  restricted diffusion  (in  the  system  above -none).  The  main  mesoscopic  feature  is  the  orientation dispersion (the level of coherence between the different orientations of the otherwise identical tensors).  In  such  a  system,  the  total  diffusion  anisotropy  as  measured  by  DTI  is  highly dependent on the degree of the mesoscopic orientation dispersion and can span the full range between zero (completely randomly oriented tensors) and one (no orientation dispersion).  For any finite orientation dispersion, anisotropy will in general appear smaller in DTI compared with the anisotropy of the underlying microscopic components. The diffusion signal decay in this system will also exhibit non-Gaussian characteristics at higher b-values, with kurtosis emerging from both the microscopic properties and the orientation dispersion. However, no link between kurtosis values and microscopic features can be established without imposing priors about the underlying tissue (Fieremans et al., 2011; Henriques et al., 2019; Jensen et al., 2005; Jensen and Helpern, 2010).

More  generally,  kurtosis  derived  from  SDE  methods is  inherently limited  in  its specificity towards microstructure since it can emerge from multiple sources such as:

- (i) anisotropic diffusion within orientationally dispersed yet otherwise identical microenvironments (e.g., an ensemble of identical diffusion tensors dispersed along an orientation distribution)  (Henriques et al., 2019; Kaden et al., 2016; Kroenke et al., 2004; Szczepankiewicz et al., 2015; Yablonskiy and Sukstanskii, 2010) -Fig. 1A;
- (ii) diffusion arising from a distribution of Gaussian  components  (isotropic tensors) in different microenvironments  (e.g.,  an  ensemble  of  isotropic  tensors  that  are

polydisperse in tensor magnitude)  (Fieremans et al., 2011; Jensen et al., 2005; Jensen and Helpern, 2010; Sukstanskii and Yablonskiy, 2002) -Fig. 1B;

(iii) restricted, time-dependent non-Gaussian diffusion (e.g., diffusion within reflecting barriers) (Callaghan, 1995; Callaghan et al., 1991; Dhital et al., 2018; Henriques et al., 2019; Jespersen, 2018; Paulsen et al., 2015) -Fig. 1C.

Kurtosis from all the above sources can be modulated by other factors such as exchange (Jensen et al., 2005; Jensen and Helpern, 2010; Kärger, 1985; Ning et al., 2018) or e.g. surface relaxation. In realistic tissues diffusional kurtosis can arise from a combination of different sources  (Fig.  1D),  however,  current  state-of-the-art  SDE  methods  cannot  resolve  their differential contributions (De Santis et al., 2014; Henriques et al., 2019, 2015; Jones et al., 2013).

Attempting to increase the specificity of diffusion-driven information, several microstructural models have been proposed for directly relating diffusion-weighted signals to tissue  properties  (Jelescu  and  Budde,  2017;  Nilsson  et  al.,  2013;  Novikov  et  al.,  2018a; Yablonskiy  and  Sukstanskii,  2010). However,  due  to  flat fitting landscapes, several assumptions  and  constraints  are  necessary  for  stabilizing  model  fitting  (Assaf  et  al.,  2004; Fieremans et al., 2011; Jelescu et al., 2016; Jespersen et al., 2007; Stanisz et al., 1997; Zhang et al., 2012). Recent studies have shown that the underlying assumptions and constraints imposed  by  these  methods  can significantly  compromise  the  specificity  of  the extracted parameters (Henriques et al., 2019; Lampinen et  al., 2019, 2017; Novikov et al., 2018b).

As an alternative to microstructural models, more specific diffusion characterization can be obtained using advanced pulse sequences (Mitra, 1995; Shemesh et al., 2016; Wong et al.,  1995).  Double  diffusion  encoding  (DDE)  sequences  have  been  proposed  to  measure

microscopic anisotropy independently of mesoscopic orientation dispersion (Callaghan  and Komlosh, 2002; Cory et al., 1990; Hui and Jensen, 2015; Jespersen et al., 2013; Lawrenz and Finsterbusch,  2013;  Mitra,  1995;  Shemesh  and  Cohen,  2011)   Multi-dimensional  diffusion . encoding (MDE) sequences (de Almeida Martins and Topgaard, 2016; Eriksson et al., 2015, 2013; Topgaard, 2015; Valette et al., 2012; Wong et al., 1995) were also developed for similar purposes, with the explicit assumption that diffusion in tissues can be represented by a sum of non-exchanging  Gaussian  diffusion  components.  In  particular,  q-space  trajectory  encoding (QTE)  was  proposed  to  resolve  kurtosis  into  two  different  sources (Lasič  et  al.,  2014; Szczepankiewicz  et  al.,  2019,  2016,  2015;  Topgaard,  2019,  2017;  Westin  et  al.,  2016):  1) anisotropic kurtosis 𝐾𝑎𝑛𝑖𝑠𝑜 which is related with the non-Gaussian signal decay arising from microscopic anisotropy; and 2) isotropic kurtosis 𝐾𝑖𝑠𝑜 which underpins the non-Gaussian signal decay  arising  from  the  variance  of  mean  diffusivities  (polydispersity  in  diffusion  tensor magnitudes). Although recent studies showed that disentangling these two sources of kurtosis can be useful to contrast tissues with disparate microstructural features such as different tumour types (Nilsson et al., 2020; Szczepankiewicz et al., 2016, 2015) the validity of the extracted measures may be compromised by factors not considered by the underlying model, such as restricted diffusion and effects of diffusion time dependence (Jespersen et al., 2019, 2018).

In this study, we aim to more generally resolve the different kurtosis sources from the cumulant expansion of double diffusion encoding signals. Importantly, DDE signals can be used  to  extract  information  about  the  displacement  correlation  tensor 𝒁 (Jespersen,  2012; Jespersen et al., 2013; Jespersen and Buhl, 2011), which itself contains direct information about isotropic and anisotropic kurtosis contributions decoupled from intra-compartmental kurtosis effects. We additionally show that an index sensitive to intra-compartmental kurtosis arising from restricted diffusion can then be inferred by subtracting the total non-Gaussian information captured by the kurtosis tensor and the other inter-compartmental kurtosis source

contributions. Given the prominence of the correlation tensor in this methodology, we term this approach Correlation Tensor Imaging (CTI). We present the theory underpinning CTI and provide its first experimental contrasts in ex vivo mouse brains and in vivo rat brains. While ex vivo experiments were designed to  assess  the  full  potential  of  CTI  with  high  quality  data acquired  without  scanning  time  constraints,  the in  vivo experiments  were  performed  to demonstrate CTI's applicability and feasibility under in vivo conditions. Potential implications for future application of CTI are discussed.

.

## 2. Methods

## 2.1. Theory

## 2.1.1. Kurtosis Sources

Multiple Gaussian diffusion approximation: When the system consists of spins diffusing in non-exchanging microenvironments exclusively characterized by individual Gaussian diffusion tensors 𝑫 𝑐 (no restricted diffusion, no time dependence), the total signal measured from an SDE experiment can be described by the following equation (expressed with Einstein summation convention):

$$E _ { \Delta } ( q ) = \langle \exp [ - q _ { i } q _ { j } \Delta D _ { i j } ^ { c } ] \rangle$$

where 𝐸∆(𝒒) is  the  normalized diffusion-weighted signal decay for a given q-vector 𝒒 with magnitude 𝑞 = 𝛾𝛿𝑔 ( 𝛾, 𝛿, and 𝑔 are the gyromagnetic ratio constant, gradient pulse duration and gradient intensity, respectively), ∆ is  the  time  interval between the  gradient of a single diffusion encoding module (diffusion time), 𝐷𝑖𝑗 𝑐 are  the  elements of  an individual diffusion tensor 𝑫 𝑐 , and 〈∙〉 represents the average over different compartments (in biological tissues this will  include  any  type  of  intra-cellular  or  extra-cellular  contributions).  In  this  study,  the superscript 'c' is used to indicate properties of an individual component/compartment (n.b., in this study, we use 'component' to reflect gaussian diffusion, 'compartment' to imply restricted diffusion, and 'microenvironment' is used quite loosely, i.e., it can reflect either components or  compartments.  None  of  the  above  is  necessarily  assigned  to  a  particular  biological component (e.g., intra/extra-cellular spaces)).

According to  the  standard  DKI  framework  (Jensen  et  al.,  2005),  the  total  diffusion tensor 𝑫 and total kurtosis tensor 𝑾 of this system can be computed by decomposing equation 1 up to fifth order in 𝑞 (or up to the second order in 𝑏 = ∆𝑞 2 ):

$$\log E _ { \Delta } ( q ) \ = - q _ { i } q _ { j } \Delta D _ { i j } + \frac { 1 } { 6 } q _ { i } q _ { j } q _ { k } q _ { l } \Delta ^ { 2 } \bar { D } ^ { 2 } W _ { i j k l } + \mathcal { O } ( q ^ { 6 } )$$

where 𝐷𝑖𝑗 are  the  elements of the total diffusion tensor, 𝑊𝑖𝑗𝑘𝑙 are  the  elements of  the  total kurtosis tensor, and 𝐷 ̅ is the mean diffusivity (i.e., 𝐷 ̅ = trace(𝑫)/3 ).

To factor out mesoscopic orientation dispersion of microenvironments (Callaghan et al., 1979; Jespersen et al., 2013; Kaden et al., 2016; Lasič et al., 2014) , it is useful to compute the  powder-averaged  decays 𝐸 ̅ ∆ from  the  average  of  diffusion-weighted  decays  measured across q-vector samples 𝒒𝑖 with different evenly distributed directions and constant magnitude 𝑞 :

$$\bar { E } _ { \Delta } ( q ) \, = \frac { 1 } { N _ { g } } \Sigma _ { i = 1 } ^ { N _ { g } } E _ { \Delta } ( q _ { i } \, )$$

where 𝑁𝑔 is  the  number  of  gradient  directions.  As  for  standard  DKI  (Jensen  et  al.,  2005), equation 3 can be expanded in cumulants to extract the total excess-kurtosis 𝐾𝑇 of the powderaveraged signal decays (Henriques, 2018; Henriques et al., 2019; Westin et al., 2016):

$$\log \bar { E } _ { \Delta } ( q ) \ = - \bar { D } \Delta q ^ { 2 } + \frac { 1 } { 6 } K _ { T } \bar { D } ^ { 2 } \Delta ^ { 2 } q ^ { 4 } + \mathcal { O } ( q ^ { 6 } )$$

with

$$K _ { T } = \frac { 6 } { 5 } \frac { \langle V _ { \lambda } ( D ^ { c } ) \rangle } { \bar { D } ^ { 2 } } + 3 \, \frac { V ( \bar { D } ^ { c } ) } { \bar { D } ^ { 2 } }$$

where 𝐷 ̅ c is  the  mean  diffusivity  of  individual  tensors  (i.e., 𝐷 ̅ c = trace(𝑫 𝑐 ) 3 ), 𝑉 𝜆 (𝑫 𝑐 ) is  the eigenvalue variance of an individual diffusion tensor 𝑫 𝑐 , and 𝑉(𝐷 c ) = 〈𝐷 ̅ ̅ c 2 〉 - 〈𝐷 ̅ c 〉 2 is the variance of mean diffusivities across microenvironments.

Equation 5 shows that, when microenvironments can be fully characterized by nonexchanging Gaussian diffusion components, 𝐾𝑇 can be fully decomposed by the anisotropic and isotropic kurtosis sources 𝐾𝑎𝑛𝑖𝑠𝑜 and 𝐾𝑖𝑠𝑜 :

$$K _ { T } \, = K _ { a n i s o } + K _ { i s o }$$

with and

$$K _ { i s o } \, = 3 \, \frac { V ( \bar { D } ^ { c } ) } { \bar { D } ^ { 2 } }$$

It  is  important  to  note  that  when  applying  equation  2  to  equations  4  and  5,  the relationship between the total average excess-kurtosis 𝐾𝑇 and the kurtosis tensor 𝑾 is given by the following expression:

$$K _ { T } \, = \, \bar { W } + \Psi$$

where 𝑊 ̅ is the mean kurtosis tensor defined by (Hansen et al., 2016, 2013), i.e.:

$$\bar { W } = \frac { 1 } { 5 } ( W _ { 1 1 1 1 } + W _ { 2 2 2 2 } + W _ { 3 3 3 3 } + 2 W _ { 1 1 2 2 } + 2 W _ { 1 1 3 3 } + 2 W _ { 2 2 3 3 } ) \quad \quad ( 1 0 )$$

and Ψ is a factor dependent on the mesoscopic orientation dispersion, that can be computed from the elements of the total diffusion tensor 𝑫 yielding:

$$\Psi = \frac { 2 } { 5 } \frac { D _ { 1 1 } ^ { \ 2 } + D _ { 2 2 } ^ { \ 2 } + D _ { 3 3 } ^ { \ 2 } + 2 D _ { 1 2 } ^ { \ 2 } + 2 D _ { 1 3 } ^ { \ 2 } + 2 D _ { 2 3 } ^ { \ 2 } } { \bar { D } ^ { 2 } } - \frac { 6 } { 5 } }$$

Intra-compartmental kurtosis effects: Equation 6 does not consider non-Gaussian (restricted, time-dependent) diffusion arising from water molecules reflections on restricting boundaries (e.g., as can occur in intra- and extra-cellular compartments in tissues). This means that every component of the system is characterized only by its individual Gaussian diffusion (c.f. Eq.1).

$$K _ { a n i s o } \, = \frac { 6 } { 5 } \frac { \langle V _ { \lambda } ( D ^ { c } ) \rangle } { \bar { D } ^ { 2 } }$$

To consider the restricted diffusion effects up to fourth order in 𝑞 , it becomes necessary to add individual intra-compartmental kurtosis tensors 𝑾 𝑐 in equation 1 (Jespersen et al., 2019):

$$E _ { \Delta } ( q ) = \langle \exp \left [ - q _ { i } q _ { j } \Delta D _ { i j } ^ { c } + \frac { 1 } { 6 } q _ { i } q _ { j } q _ { k } q _ { l } \Delta ^ { 2 } \bar { D } ^ { c } ^ { 2 } W _ { i j k l } ^ { c } \, \right ] \rangle$$

where 𝑊𝑖𝑗𝑘𝑙 𝑐 are the elements of the individual intra-compartmental kurtosis tensor 𝑾 𝑐 arising from  restricted  or  time-dependent  diffusion.  We  also  note  that  for  compartments  with  a restricted diffusion-driven 𝑾 𝑐 , 𝐷𝑖𝑗 𝑐 &lt; 𝐷0 also reflects the time-dependent restricted diffusion effects  on  the  diffusion  tensor  which  are  omitted  for  clarity  in  equation  12, 𝐷0 is  the  free diffusion coefficient.

Applying  this  new  representation  to  equations  3-4  and  computing  the  cumulant expansion up to the fourth order in 𝑞 , the total excess-kurtosis 𝐾𝑇 is now given by (Jespersen et al., 2019):

$$K _ { T } = \frac { 6 } { 5 } \frac { \langle V _ { \lambda } ( \bar { D } ^ { c } ) \rangle } { \bar { D } ^ { 2 } } + 3 \, \frac { V ( \bar { D } ^ { c } ) } { \bar { D } ^ { 2 } } + \frac { \langle ( \bar { D } ^ { c } ) ^ { 2 } K _ { T } ^ { c } \rangle } { \bar { D } ^ { 2 } }$$

where 𝐾𝑇 𝑐 is the excess-kurtosis of powder averaged signals for an individual compartment 'c' , which can be computed as (𝑊1111 𝑐 +𝑊2222 + 𝑊3333 + 2𝑊 1122 𝑐 𝑐 𝑐 +2𝑊1133 + 2𝑊 2233 )/5 𝑐 𝑐 .  We note in passing that the individual 𝐾𝑇 𝑐 arises solely from restriction, while 𝐾𝑇 emerges from all sources as explicitly shown in equation 13.

Equation  13  shows  that,  when  non-Gaussian  diffusion  arises  from  reflection  at boundaries, the apparent total excess-kurtosis 𝐾𝑇 encompasses contributions from the apparent isotropic  and  anisotropic  kurtosis  sources  ( 𝐾𝑎𝑛𝑖𝑠𝑜 and 𝐾𝑖𝑠𝑜 ,  c.f.  equations  7  and  8)  and contributions from individual excess-kurtosis 𝐾𝑇 𝑐 averaged across all compartments -that will be referred to as the intra-compartmental kurtosis ( 𝐾𝑖𝑛𝑡𝑟𝑎 ) - i.e.:

$$K _ { T } \, = \, K _ { a n i s o } + K _ { i s o } + K _ { i n t r a }$$

with

$$K _ { i n t r a } \, = \frac { \langle ( \bar { D } ^ { c } ) ^ { 2 } K _ { T } ^ { c } \rangle } { \bar { D } ^ { 2 } },$$

Note that 𝐾𝑖𝑛𝑡𝑟𝑎 is not based on a specific microstructural model and thus it generally captures information from all restricting compartments, e.g., in white or gray matter tissues, whether arising,  e.g.,  from  intra-axonal,  intra-cellular, or  extra-cellular  spaces.  However, contributions  from  elements  with  higher  mean  diffusivity 𝐷 ̅ c will  have  contribute  more  to 𝐾𝑖𝑛𝑡𝑟𝑎 as apparent from the numerator in equation 15.

## 2.1.2. Correlation tensor imaging

The correlation  tensor imaging  (CTI) framework  is based on  the  correlation tensors  in the cumulant expansion of double diffusion encoding (DDE) signals. Fig. 2A shows an illustration of the DDE gradient waveform and its parameters. DDE comprises two diffusion encoding modules characterized by different q-vectors ( 𝒒1 and 𝒒2 ) and diffusion times ( ∆1 and ∆2 ). The time interval between the two diffusion encoding modules is termed the mixing time ( 𝜏𝑚 ). In this study, both diffusion times are equal, ∆1= ∆1≡ ∆ .

Although, up to the fourth order in 𝑞 , DDE signals can be related to a single 6th order kurtosis tensor (Hui and Jensen, 2015; Jensen et al., 2014), our correlation tensor imaging (CTI) approach is based on DDE's cumulant expansion formulated by Jespersen (2012) and the signal is expressed in terms of five unique second- and fourth-order tensors:

$$\log E _ { \Delta } ( q _ { 1 }, q _ { 2 } ) \, = - \left ( q _ { 1 i } q _ { 1 j } + q _ { 2 i } q _ { 2 j } \right ) \Delta D _ { i j } + q _ { 1 i } q _ { 2 j } Q _ { i j }$$

$$+ \frac { 1 } { 6 } ( q _ { 1 i } q _ { 1 j } q _ { 1 k } q _ { 1 l } + q _ { 2 i } q _ { 2 j } q _ { 2 k } q _ { 2 l } ) \Delta ^ { 2 } \overline { D } ^ { 2 } W _ { i j k l }$$

$$+ \frac { 1 } { 4 } q _ { 1 i } q _ { 1 j } q _ { 2 k } q _ { 2 l } Z _ { i j k l }$$

$$+ \frac { 1 } { 6 } ( q _ { 1 i } q _ { 1 j } q _ { 1 k } q _ { 2 l } + q _ { 2 i } q _ { 2 j } q _ { 2 k } q _ { 1 l } ) S _ { i j k l } + \mathcal { O } ( q ^ { 6 } )$$

(16)

where, in addition to elements of the diffusion and kurtosis tensors ( 𝐷𝑖𝑗 and 𝑊𝑖𝑗𝑘𝑙 ), 𝑄𝑖𝑗 are the elements of a 2nd order correlation tensor 𝑸 which encodes information on the time dependence of 𝑫 , and 𝑍𝑖𝑗𝑘𝑙 and 𝑆𝑖𝑗𝑘𝑙 are the elements of the 4th order displacement correlation tensors 𝒁 and 𝑺 (Jespersen, 2012).

## 2.1.3. Kurtosis separation using CTI

Anisotropic and isotropic kurtosis sources:  At the long mixing time regime (and in the absence of  flow  and  exchange  (Jespersen,  2012)), 𝒁 becomes  proportional  to  the  diffusion  tensor covariance 𝑪 (Jespersen  et  al.,  2013;  Topgaard,  2017;  Valiullin,  2017).  Explicitly, 𝑍 𝑖𝑗𝑘𝑙 → 4∆ 𝟐 𝐶𝑖𝑗𝑘𝑙 - which can then be used to generally estimate the ensemble averaged microscopic anisotropy 〈𝑉 𝜆 (𝑫 𝑐 )〉 , and variance of mean diffusivities 𝑉(𝐷 c ̅ ) using the following equations (Jespersen et al., 2013; Topgaard, 2017; Valiullin, 2017):

$$\langle V _ { \lambda } ( D ^ { c } ) \rangle = & \, \frac { 2 } { 9 } \, \left [ C _ { 1 1 1 } + D _ { 1 1 } \, ^ { 2 } + C _ { 2 2 2 } + D _ { 2 2 } \, ^ { 2 } + C _ { 3 3 3 3 } + D _ { 3 3 } \, ^ { 2 } - C _ { 1 1 2 2 } - D _ { 1 1 } D _ { 2 2 } - C _ { 1 1 3 3 } \\ & - D _ { 1 1 } D _ { 3 3 } - C _ { 2 2 3 3 } - D _ { 2 2 } D _ { 3 3 } \\ & + 3 ( C _ { 1 2 1 2 } + D _ { 1 2 } \, ^ { 2 } + C _ { 1 3 1 3 } + D _ { 1 3 } \, ^ { 2 } + C _ { 2 3 2 3 } + D _ { 2 3 } \, ^ { 2 } ) \right ]$$

(17)

and

$$V ( \bar { D } ^ { \text{c} } ) = \frac { 1 } { 9 } ( C _ { 1 1 1 1 } + C _ { 2 2 2 2 } + C _ { 3 3 3 3 } + 2 C _ { 1 1 2 2 } + 2 C _ { 1 1 3 3 } + 2 C _ { 2 2 3 3 } )$$

It is important to note that equations 17 and 18 do not rely on the multiple Gaussian diffusion assumption  and  thus,  in  the  absence  of  flow  and  exchange, 〈𝑉 𝜆 (𝑫 𝑐 )〉 and 𝑉(𝐷 c ̅ ) can  be generally converted to the anisotropic and isotropic kurtosis sources using equations 7 and 8.

Intra-compartmental kurtosis index: Given that the total excess-kurtosis 𝐾𝑇 can be obtained from 𝑾 and 𝑫 (eq. 9) and that the two inter-compartmental kurtosis sources 𝐾𝑎𝑛𝑖𝑠𝑜 and 𝐾𝑖𝑠𝑜

can be derived from 𝒁 and 𝑫 (eqs. 7-8, 17-18), an index sensitive to the intra-compartmental kurtosis can be computed from the following simple subtraction:

$$K _ { i n t r a } = K _ { T } - K _ { a n i s o } - K _ { i s o }$$

since 𝐾𝑇 contains the information on all three kurtosis sources.

Higher order terms: Note that, since CTI is based on the truncated cumulant expansion of DDE signals, the accuracy of all indices, and especially the 𝐾𝑖𝑛𝑡𝑟𝑎 index, may be compromised by the higher order effects neglected in the above expressions (c.f. 𝒪(𝑞 6 ) in Eq. 16), particularly if  the  data  are not acquired at sufficiently low b-values (Chuhutin et al., 2017; Ianuş et al., 2018). In this study, we only explored these higher order effects empirically (i.e., by probing b-value dependencies and using simulations to provide intuition). Exploring ways to minimize these effects was beyond the scope of this study and will be reported in due course.

## 2.1.4. Acquisition requirements for CTI

Suppressing 𝑄𝑖𝑗 and 𝑆𝑖𝑗𝑘𝑙 :  As shown above, only 𝑫 𝑾 , and 𝒁 are  necessary to resolve the contributions of the different kurtosis sources. DDE has an appealing 'built in' suppressor of the 𝑸 and 𝑺 tensors via the mixing time: at long 𝜏𝑚 , both tensors vanish (Jespersen, 2012). Alternatively, these tensors can be cancelled out by combining DDE acquisitions with inverted 𝒒2 vectors:

$$\frac { \log E _ { \Delta } ( q _ { 1 }, q _ { 2 } ) } { 2 } + \frac { \log E _ { \Delta } ( q _ { 1 }, - q _ { 2 } ) } { 2 } =$$

$$- ( q _ { 1 i } q _ { 1 j } + q _ { 2 i } q _ { 2 j } ) \Delta D _ { i j }$$

$$+ \frac { 1 } { 6 } ( q _ { 1 i } q _ { 1 j } q _ { 1 k } q _ { 1 l } + q _ { 2 i } q _ { 2 j } q _ { 2 k } q _ { 2 l } ) \Delta ^ { 2 } \overline { D } ^ { 2 } W _ { i j k l }$$

$$+ \frac { 1 } { 4 } q _ { 1 i } q _ { 1 j } q _ { 2 k } q _ { 2 l } Z _ { i j k l } + \mathcal { O } ( q ^ { 6 } )$$

To reduce the number of parameters to be fitted, the contributions of tensors Q and 𝑺 were here suppressed in this way.

Diffusion-weighting  requirements:  Analogously  to  Diffusional  Kurtosis  Imaging (DKI),  fitting  tensors  associated  with  the 𝑞 4 cumulant  (i.e. 𝑾 and 𝒁 ),  CTI  requires  data acquired  with  at  least  three  different  diffusion  gradient  intensities  (shells).  In  addition, asymmetric DDE gradients are required for decoupling of elements 𝑊𝑖𝑗𝑘𝑙 and 𝑍𝑖𝑗𝑘𝑙 .  In  this study, a minimal protocol for CTI was designed based on the following eight gradient intensity combinations of 𝑞 ∈ {0, 𝑞 𝑎 , 𝑞 𝑏 } (Fig. 2B):

- 1) |𝒒1 |, |𝒒 2 | = 𝑞 𝑎 , 𝑞 𝑎 ;
- 2) |𝒒1 |, |𝒒 2 | = 𝑞 𝑎 , 𝑞 𝑎 , with inverted 𝒒2 direction;
- 3) |𝒒1 |, |𝒒 2 | = 𝑞 𝑎 , 0 ;
- 4) |𝒒1 |, |𝒒 2 | = 0, 𝑞 𝑎 ;
- 5) |𝒒1 |, |𝒒 2 | = 𝑞 𝑏 , 𝑞 𝑏 ;
- 6) |𝒒1 |, |𝒒 2 | = 𝑞 𝑏 , 𝑞 𝑏 , with inverted 𝒒2 direction;
- 7) |𝒒1 |, |𝒒 2 | = 𝑞 𝑏 , 0 ; and
- 8) |𝒒1 |, |𝒒 2 | = 0, 𝑞 𝑏 .

The magnitudes 𝑞𝑎 and 𝑞𝑏 were defined for a given 𝑏𝑚𝑎𝑥 value ( 𝑞𝑎 = √𝑏𝑚𝑎𝑥/2(∆ - 𝛿/3) and 𝑞𝑏 = √𝑏𝑚𝑎𝑥/(∆ - 𝛿/3) ).

Gradient direction requirements: To resolve the anisotropic information of 𝒁 , different pairs of gradient directions are required (Jespersen et al., 2013). The orientations defined in the DDE 5-design (Jespersen et al., 2013) are adequate for this purpose, namely, twelve pairs of parallel 𝒒1 -𝒒2 directions and sixty pairs of perpendicular 𝒒1 -𝒒2 directions. In this study, to decrease the difference between the number of parallel and perpendicular 𝒒1 -𝒒2 directions, 45

extra DDE experiments with parallel 𝒒1 -𝒒2 directions were acquired, yielding a total of 117 (57 parallel + 60 perpendicular) 𝒒1 -𝒒2 combination of directions. The parallel directions of these latter 𝒒 1 -𝒒2 parallel pairs were evenly sampled on a spherical 3-dimensional grid. All 117 direction combinations were repeated for the 𝒒1 -𝒒2 magnitude combinations 1-8 described above.

## 2.2. MRI experiments

All  animal  experiments  were  preapproved  by  the  competent  institutional  and  national authorities, and carried out according to European Directive 2010/63.

Ex vivo CTI experiments. Brain specimens were extracted via transcardial perfusion with 4% Paraformaldehyde (PFA) from two adult mice (N=2 C57BL/6J males 13 weeks old, weights 23/24g, respectively, grown with a 12h/12h light/dark cycle with ad libitum access to food  and  water).  After  extraction  from  the  skull,  both  brains  were  immersed  in  4%  PFA solution for 24 h, and then washed in Phosphate-Buffered Saline (PBS) solution for at least 48 h.  The  specimens  were  then  placed  in  a  10-mm  NMR  tube  filled  with  Fluorinert  (Sigma Aldrich, Lisbon, PT), secured with a stopper from above to prevent floating, and the NMR tube was sealed using paraffin film. MRI scans were prefomed on a 16.4 T Aeon Ascend Bruker scanner (Karlsruhe, Germany) equipped with an AVANCE IIIHD console, and a Micro5 probe with gradient coils capable of producing up to 3000 mT/m in all directions. Using the probe's variable temperature capability, we maintained the samples at 37oC. The samples were allowed to equilibrate with the surroundings for at least 3 h prior to commencement of diffusion MRI experiments.

Double diffusion encoding data was acquired for five coronal slices using an in house written EPI-based DDE pulse sequence. The diffusion encoding gradient pulse separation 

and mixing time 𝜏𝑚 were set to 13 ms, while the pulsed gradient duration  was set to 1.5 ms (Fig.  2A).  Data  were  acquired  for  the  minimal  CTI  protocol  using  the  117  DDE  pairs  of directions and repeated for all eight 𝒒1 -𝒒2 magnitude combinations (parameters described in section  2.1.4), in addition to sixty acquisitions without any diffusion-weighted sensitization (zero b-value). To assess the effects of higher order terms not considered by CTI, the entire set of minimal protocol acquisitions were repeated for seven evenly sampled 𝑏𝑚𝑎𝑥 values (1.00, 1.25,  1.50,  1.75,  2.00,  2.25,  and  2.5  ms/  m2).  For  all  experiments,  the  following  common parameters were used: TR/TE = 2200/52 ms, Field of View = 10.4  10.4 mm2, matrix size 80  80, leading to an in-plane voxel resolution of 130  130  m2, slice thickness = 0.9 mm, slice gap = 0.6 mm, number of segments = 2, number of averages = 8, partial fourier  effective acceleration = 1.42. For each 𝑏𝑚𝑎𝑥 value, the total acquisition time was about 9.75 h.

In  addition  to  the  diffusion-weighted  data,  coronal  T2-weighted  images  with  high resolution and high SNR were acquired for anatomical reference and assessment of partial volume effects. This data was acquired using a RARE sequence with the following parameters: TR = 4250 ms, effective TE = 22 ms, RARE factor = 8, Field of View = 10  10 mm2, matrix size 126  126, in-plane voxel resolution = 79.4  79.4  m2, slice thickness = 79.4  m (40 slices), number of averages = 230, partial fourier effective acceleration = 1.05.

In vivo CTI experiments. To assess CTI's ability to characterize tissues in vivo within feasible scan times, data was acquired on two rats in vivo (N=2 Long Evans females, 14/15 weeks old, weights 264/254 g, respectively, grown in a 12h/12h light/dark cycle with ad libitum access to food and water) under anesthesia (Isoflurane ~2.5% in 28% oxygen). All in vivo data were acquired on a 9.4 T Bruker Biospec MRI scanner equipped with an 86 mm quadrature coil  for  transmission  and  a  4-element  array  cryocoil  for  reception.  Due  to  acquisition  time constraints, in vivo double diffusion encoding data were acquired only for one CTI minimal

protocol with 𝑏𝑚𝑎𝑥 = 2 ms/  m2  (  = 12 ms,  = 12 ms,  = 3 ms, 117 pairs of direction for 8 gradient  intensity  combinations,  c.f.  Fig.  2B)  for  three  coronal  images  using  per  slice respiratory gating. A large number of 180 b-values = 0 data were also acquired to ensure a high ratio between the number of non-diffusion and diffusion-weighted acquisitions (Alexander and Barker, 2005; Jones et al., 1999)). Other acquisition parameters included: TR/TE = 3000 / 48.5 ms, Field of View = 20  20 mm2, matrix size 100  100, in-plane resolution of 200  200  m2, slice thickness = 1 mm, slice gap = 1.8 mm, number of EPI segments = 1, number of averages = 2, partial fourier effective acceleration = 1.40. For each animal, the acquisition time of all of the diffusion data was about 2 hours.

Data processing: All diffusion-weighted datasets were first preprocessed by realigning the data using a sub-pixel registration technique (Guizar-Sicairos et al., 2008). CTI was then directly fitted to the data using a weighted-linear-least squares fitting procedure, implemented in-house, to equation 20. No terms associated with cumulants higher than the fifth order in 𝑞 were considered in our fitting procedure to ensure the precision of kurtosis estimates (Kiselev, 2017). Two different analyses were performed for  each ex  vivo mouse brain dataset: (i) to assess the contrasts of different kurtosis sources with maximum precision, CTI estimates were obtained in a single fit incorporating all diffusion-weighted datasets of all 𝑏𝑚𝑎𝑥 values, making a  total  of  56  combination  of 𝒒1 -𝒒2 magnitudes and 420 b-value = 0 acquisitions; (ii) CTI estimates  were  additionally  produced  for  individual 𝑏𝑚𝑎𝑥 sets  in  the  minimal  protocols acquired, to test the robustness of kurtosis estimates with different b-values and to assess the effects  of  higher  order  terms.  For  the in  vivo acquisitions,  CTI  estimates  of  each  rat  brain dataset  was  processed  for  the  single  acquired  minimal  protocol  defined  with  a 𝑏𝑚𝑎𝑥 = 2 ms/  m 2 . In addition to visual inspection of the CTI derived maps, kurtosis estimates were also extracted from regions of interest (ROIs) which were bilateratly drawn on all data slices.

## 2.3. Simulations

To support the interpretation of the results from the MRI experiments, the CTI approach was also subjected to extensive numerical simulations using noise-free synthetic signals in which ground truth kurtosis sources are known a priori (simulations with added synthetic noise are reported in supplementary material). The simulations were performed using the minimal CTI protocol defined by the eight gradient intensity combinations decribed in Fig. 2B, which were repeated for the gradient directions of Jespersen's 5-design  in  addition  to  45  extra  parallel directions. It is important to note that the eight gradient intensity combinations decribed in Fig. 2B are defined for a given 𝑏𝑚𝑎𝑥 value, and thus apparent kurtosis estimates can be obtained as a  function  of 𝑏𝑚𝑎𝑥 .  To  assess  the  effects  of  higher  order  terms,  kurtosis  estimates  were independently produced for seven different 𝑏𝑚𝑎𝑥 values (i.e., 1.00, 1.25, 1.50, 1.75, 2.00, 2.25, and 2.5 ms/  m2). The diffusion gradient separation times, pulse durations and mixing time for these simulations were set to  = 13 ms,  = 1.5 ms, and 𝜏𝑚 = 13 ms, respectively (the same parameters used in the ex vivo mouse brain experiments). The synthetic signals were produced according to two different ground-truth scenarios comprising a mix of Gaussian components and a mix of Gaussian and/or restricted compartments:

1) Gaussian microenvironements according to a two-component model: DDE signals were first produced for two well-aligned axially-symmetric Gaussian diffusion components (no restricted diffusion). The axial and radial diffusivities for the first component were set to 2 and 0  m2/ms, while the axial and radial diffusivities for the second component were set to 1.5 and 0.5  m2/ms. Volume fractions for both components were set to 0.5. Based on these values, ground truth kurtosis was computed using equations 7 and 8. Note that for these  simulations 𝐾𝑖𝑛𝑡𝑟𝑎 = 0 .  To  assess  the  dependence  of  the  signals  on  mesoscopic orientations, simulations were also repeated for different levels of orientation dispersion.

For this, simulations of 10000 replicas of the symmetric Gaussian diffusion components were  produced.  The  directions  of  these  replicas  were  sampled  based  on  a  Watson distribution (Watson, 1965) which can be produced with arbitrary dispersion levels. For this  study,  different  dispersion  levels  were  tested  by  changing  the  Watson  distribution concentration  parameter 𝜅 from  0  to  16.58,  where 𝜅 =  0  corresponds  to  completely randomly oriented microenvironments, while 𝜅 = 16.58 corresponds to a low dispersion of  10 degrees (according to  the  dispersion  angle definition proposed by  (Riffert  et  al., 2014)).

2) Gaussian  component  and  restricted  microenvironements:  synthetic  DDE  signals  for non-zero 𝐾𝑖𝑛𝑡𝑟𝑎 were  produced  by  incorporating  a  spherical  compartment  (restricted diffusion) to the two Gaussian components described above. The signals for the restricted compartment  were  produced  using  the  MISST  package  (Drobnjak  et  al.,  2011). Simulations were repeated for three different sphere diameters ( 𝑑𝑖 = 5, 7.5, and 10  m, approximately mimicking typical cell soma sizes in the nervous system) and the volume fractions for both Gaussian components and for the restricted compartment were set to have equal contributions (i.e. 𝑓 1 = 𝑓 2 = 𝑓 3 = 1/3 ). Other simulation parameters were as follows: intrisic diffusivity = 2  m2/ms; simulations sampling time = 0.015 ms; number of spherical Bessel orders = 70 (Callaghan, 1997; Drobnjak et al., 2011). The ground truth individual mean diffusivity 𝐷 ̅ c and individual excess-kurtosis 𝐾 ̅ c of the spherical restricted compartment  were  determined  using  extra  SDE  simulations.  To  achieve  the  apparent values  for 𝑞 → 0 ,  these  ground  truth  values  were  computed  by  fitting  the  standard diffusional kurtosis imaging (DKI) equation (Jensen et al., 2005) to the synthetic signals simulated  for  a  maximum  gradient  intensity  set  to 𝑞𝑚𝑎𝑥 = 0.1/ 𝑑𝑖 . To  avoid  fitting instabilities,  SDE  synthetic  signals  were  evenly  sampled  for  2500  b-values  from  0  to 𝑏𝑚𝑎𝑥 = (2𝜋𝑞𝑚𝑎𝑥 ) 2 (∆ - 𝛿/3) .  The  ground  truth 𝐾𝑎𝑛𝑖𝑠𝑜 , 𝐾𝑖𝑠𝑜 and 𝐾𝑖𝑛𝑡𝑟𝑎 for  the  total

synthetic signals were computed using equations 7, 8, and 15. These simulations were also repeated for different levels of orientation dispersion of the axially-symmetric compnents. For this, 10000 replicas of the Gaussian and spherical compartments were sampled based on  a  Watson  distribution  with  a  varying  dispersion  level.  For  the  sake  of  simplicity, simulations of restricted spherical compartments and dispersed Gaussian compartments were only produced for the larger sphere diameter of 10  m, which corresponds to the scenario with larger 𝐾𝑖𝑛𝑡𝑟𝑎 .

## 3. Results

## 3.1. MRI experiments

Raw data of ex vivo DDE-MRI experiments for different b-values are presented in Fig. 3. ROIs placed in white matter (WM), grey matter (GM) and cerebral ventricles (CV) at b = 0 (Fig. 3A) exhibited signal-to-noise ratios (SNRs) of 70  20, 110  10, and 192  2 for all WM, GM, and CV ROIs, respectively. For a representative slice, Fig. 3B-D shows the powder-averaged DDE signal  decays  at  total  b-value  of  1,  3  and  5  ms/  m2,  respectively,  and  for  parallel  DDE experiments ( 𝐸 ̅ ∆ (𝒒1 ↑↑ 𝒒2 ) , Fig. 3 B1, C1, D1), anti-parallel DDE experiments ( 𝐸 ̅ ∆ (𝒒1 ↑↓ 𝒒2 ) , Fig. 3 B2, C2, D2), and perpendicular DDE experiments ( 𝐸 ̅ ∆ (𝒒1 ⊥ 𝒒2), Fig. 3 B3, C3, D3), respectively. At these b-values, individual diffusion-weighted images were characterized by SNRs of 41  9, 19  1, and 10  1 for all WM ROIs, and SNRs of 62  4, 20  1, and 8  1 for all GM ROIs (the SNR values are from powder-averaged parallel DDE experiments, Fig. 3 B1, C1, D1).

An important assumption of CTI, as presented here, is that the long mixing time regime has been achieved. This can be empirically tested by comparing data acquired with parallel and anti-parallel diffusion pairs. Maps corresponding to the signal ratios between parallel and antiparallel DDE measurements are shows in panels B4, C4, and D4 of Fig. 3. The ratio maps show values near unity, indicating that the long mixing time regime assumption ( lim  𝐸 (𝒒 1 𝜏→∞ ̅ ∆ ↑↑ 𝒒2 )/ 𝐸 (𝒒1 ̅ ∆ ↑↓ 𝒒2 ) → 1 , (Jespersen and Buhl, 2011; Koch and Finsterbusch, 2008; Ozarslan, 2009)) is practically fulfilled for most voxels. For comparison, the maps of the ratio between parallel and perpendicular DDE signals ( 𝐸 ̅ ∆ (𝒒1 ↑↑ 𝒒2 )/ 𝐸 ∆ ̅ (𝒒1 ⊥ 𝒒2 ) ) are shown in panels B5, C5, and D5  of  Fig.  3.  These  reveal  the  expected  higher  values  for  white  matter  regions  where

microscopic anisotropy is known to be higher. As predicted by (Ianuş et al., 2018; Jespersen et al., 2013), this contrast increased with higher b-value (Fig. 3B5, Fig. 3C5, and Fig. 3D5).

Given the robustness of the raw data and the fulfillment of the long mixing time regime, the  correlation  tensor  metrics  were  first  extracted  from  the  extensively  sampled  b-value protocol (i.e. data were fit using all 56 acquired 𝒒1 -𝒒2 magnitude combinations together). Fig. 4  presents  the  kurtosis  source  separation  maps  for  each  of  the  five  slices  acquired  for  a representative  mouse  brain. 𝐾𝑇 , 𝐾𝑎𝑛𝑖𝑠𝑜 , 𝐾𝑖𝑠𝑜 ,  and 𝐾𝑖𝑛𝑡𝑟𝑎 estimates  (shown  in  panels  A-D, respectively) evidenced drastically different contrasts. Notably, 𝐾𝑇 is, as expected, higher than any of its sources (Fig. 4A). 𝐾𝑎𝑛𝑖𝑠𝑜 is revealed to be the largest source contributing to the total kurtosis in white matter (e.g. regions pointed by white arrows, Fig. 4B). On the other hand, 𝐾𝑖𝑠𝑜 shows relatively low intensities for both white and grey matter, with the exception of areas where partial volume effects arising from free water in cerebral ventricles are dominant (e.g. regions  pointed  by  grey  arrows,  Fig.  4C).  Large  partial  volume  effects  in  these  areas  are supported  by  high-resolution  anatomical  mapping  (supplementary  Fig.  S1),  that  directly showed areas with free water overlapping with other tissues. 𝐾𝑖𝑛𝑡𝑟𝑎 maps mainly exhibited positive values (25th and 75th percentiles for all data voxels were 0.286 and 0.430, respectively). Relative to 𝐾𝑎𝑛𝑖𝑠𝑜 , both 𝐾𝑖𝑠𝑜 and 𝐾𝑖𝑛𝑡𝑟𝑎 maps seem to be more sensitive to image artefacts (e.g. dorsal-ventral oscillations pointed by the red arrows, which are likely to be a combined effect of Gibbs Ringing and partial Fourier). All these results were found to be consistent between the N=2 ex vivo mouse brains scanned, as shown in supplementary Fig. S2. For instance, 𝐾𝑖𝑛𝑡𝑟𝑎 maps for mouse brain #2 show consistent positive values (25th and 75th percentiles for all data voxels are 0.320 and 0.431, respectively).

Kurtosis estimates extracted for individual 𝑏𝑚𝑎𝑥 minimal protocols are shown in Fig. 5, namely, 𝐾𝑇 , 𝐾𝑎𝑛𝑖𝑠𝑜 , 𝐾𝑖𝑠𝑜 , and 𝐾𝑖𝑛𝑡𝑟𝑎 maps are displayed for both mice. The kurtosis mean values and standard deviation across the animals are plotted as a function of the protocol values

𝑏𝑚𝑎𝑥 in Fig. 5A3-D3. For reference, the mean kurtosis estimates of the extensively sampled bvalue protocol are plotted  (dotted lines) in these latter panels. Across the two mouse brain specimens, 𝐾𝑇 maps  are  consistent  only  for  the  higher  b-values  and  evidence  implausible negative values at low b-values (Fig. 5A1-2). 𝐾𝑎𝑛𝑖𝑠𝑜 estimates decreased as 𝑏𝑚𝑎𝑥 increased (Fig. 5B) -this b-value dependence can be particularly appreciated by observing the 𝐾𝑎𝑛𝑖𝑠𝑜 estimates extracted from the white matter ROIs (Fig. 5B3). Similarly to the 𝐾𝑇 maps, 𝐾𝑖𝑠𝑜 maps were visually nosier at lower b-values (Fig. 5C1 and Fig. 5C2). The b-value dependence for mean 𝐾𝑖𝑠𝑜 is  visually less obvious than the other kurtosis estimates (Fig. 5C3). 𝐾𝑖𝑛𝑡𝑟𝑎 maps consistently exhibited positive values at the higher b-values; however, negative values were present in the noisier 𝐾𝑖𝑛𝑡𝑟𝑎 maps at lower 𝑏𝑚𝑎𝑥 values (Fig. 5D).

Results from the in vivo rat experiments are shown in Fig. 6. White matter and grey matter ROIs placed in  the superior regions of the rat brain (green and red ROIs in Fig. 6) exhibited, respectively, SNRs of 41  5 and 41  4 for the non-diffusion weighted data, SNRs of 12  2 and 10  2 for total b-value = 2 ms/  m2, and SNRs of 5.3  0.9 and 3.1  0.4 for total b-value = 4 ms/  m2. Inferior areas of the brain in our in vivo data (regions associated with a lower cryocoil sensitivity) showed lower signal intensities and thus lower SNRs, particularly, a grey matter ROIs placed in the inferior region of the brain (blue ROI in Fig. 6A) evidenced SNRs of 26  4, 6.2  0.8, and 2.4  0.4 for the total b-values of 0, 2 and 4 ms/  m2, respectively. 𝐾𝑇 , 𝐾𝑎𝑛𝑖𝑠𝑜 , 𝐾𝑖𝑠𝑜 , and 𝐾𝑖𝑛𝑡𝑟𝑎 maps for all slices and for both animals are shown in Fig. 6 B-E. Note  that  different  kurtosis  types  are  displayed  with  different  color  bar  ranges  for  better contrast visualization. Consistent with the ex vivo CTI results, 𝐾𝑇 was higher than any of its sources  (Fig.  6B)  and 𝐾𝑎𝑛𝑖𝑠𝑜 dominated  in  white  matter  (Fig.  6C). In  vivo 𝐾𝑖𝑠𝑜 maps  also demonstrated the abovementioned sensitivity towards partial volume effects between cerebral tissue and cerebral ventricle (free water pointed by black arrows), but also appeared higher in WM (white arrows) (Fig. 6D). High 𝐾𝑖𝑠𝑜 estimates were also present  in the grey matter of

inferior brain regions (pointed by red arrows). For both in vivo rats, 𝐾𝑖𝑛𝑡𝑟𝑎 was consistently positive (Fig. 6E, 25th-75th percentile range of all slice voxels is 0.21-0.45 for Rat 1 and 0.180.43 for Rat 2).

## 3.2. Numerical simulations

To further validate CTI and investigate how the different sources of kurtosis would vary with b-value,  numerical  simulations  for  several  plausible  diffusion  conditions  were  performed . When the system consists of Gaussian components, namely, perfectly aligned 'sticks' (with zero radial diffusivity) and a tensor (Fig. 7A), 𝐾𝑇 and 𝐾𝑎𝑛𝑖𝑠𝑜 extracted from CTI shows a weak dependence on b-value; in particular, the biases introduced by higher order effects are lower than 1.2% (Fig. 7A1-2) . On the other hand , the apparent 𝐾𝑖𝑠𝑜 extracted from CTI appears lower with increasing b-value due to higher order terms, reaching negative biases larger than 50% at 𝑏𝑚𝑎𝑥 = 2.5 ms/  m2 (Fig.  7A3).  In  this  system, 𝐾𝑖𝑛𝑡𝑟𝑎 ground  truth  is  identically  zero; however, higher order terms induce a positive apparent 𝐾𝑖𝑛𝑡𝑟𝑎 (maximum bias of ~0.01 for 𝑏𝑚𝑎𝑥 = 2.5 ms/  m2, Fig. 7A4). When orientation dispersion is added to the same system (Fig. 7B), the extracted parameters have different dependencies on b-value. 𝐾𝑇 and 𝐾𝑎𝑛𝑖𝑠𝑜 extracted from  CTI  approach  their  nominal  ground  truth  only  at  b  =  0,  while  they  are  increasingly underestimated  at  higher  b-values  (Fig.  7B1-2). 𝐾𝑖𝑠𝑜 and 𝐾𝑖𝑛𝑡𝑟𝑎 extracted  from  CTI  now become increasingly overestimated with higher b-values (Fig. 7B3-4). Particularly, 𝐾𝑖𝑛𝑡𝑟𝑎 bias becomes  positive,  higher  than  0.2  for  the  larger 𝑏𝑚𝑎𝑥 , when  Gaussian  components  are completely randomly oriented (yellow curve of Fig. 7B4)

To  assess  how  intra-compartmental  restricted  diffusion  may  affect  the  b-value dependence  of  the  extracted  parameters,  simulations  incorporating  an  impermeable  sphere along with the stick and tensor model above, were performed (Fig.  8). Without orientation

dispersion (Fig. 8A), 𝐾𝑇 , 𝐾𝑎𝑛𝑖𝑠𝑜 and 𝐾𝑖𝑠𝑜 values approach their ground-truth at low b-values but  are  all  underestimated  at  higher  b-values  (Fig.  8A1-3).  Although  large  positive  bias  is observed at higher b-values, apparent 𝐾𝑖𝑛𝑡𝑟𝑎 values are negative for low b-values as expected for the negative concave signal decays profiles of restricted diffusion (Fig. 8A4). 𝐾𝑖𝑛𝑡𝑟𝑎 bias trends depend on the size of the sphere, yet apparent 𝐾𝑖𝑛𝑡𝑟𝑎 values approached their respective negative ground truth values at low b-values (Fig. 8A4). Fig. 8B shows the simulations when a restricted sphere of 10  m is added to orientationally-dispersed replicas of the two-Gaussian compartments; the trends are similar to those reported above, although the underestimation of 𝐾𝑖𝑛𝑡𝑟𝑎 becomes higher in this scenario (Fig. 8B4).

## 4. Discussion

Since its inception, DKI has played an important role in microstructural characterizations. In many cases, kurtosis measurements appeared more sensitive to disease or other processes, such as development and ageing, compared with their diffusion tensor counterparts (Cheung et al., 2012; Falangola et al., 2008; Henriques, 2018; Lin et al., 2018; Rudrapatna et al., 2014; Sun et al., 2015). Clinical applications of diffusional kurtosis MRI abound, and deeper investigations into kurtosis features, such as its time dependence, are being vigorously studied (Grussu et al., 2019; Jensen and Helpern, 2010; Jespersen et al., 2018; Lee et al., 2019; Pyatigorskaya et al., 2014). Nearly invariably, these measurements are performed using single diffusion encoding pulse  sequences;  however,  SDE  methods  cannot  separate  the  different  sources  of  kurtosis, which  reduces  the  specificity  of  DKI.  QTE  approaches  have  recently  been  gaining  much interest for their ability to portray anisotropic and isotropic kurtosis sources separately (Lasič et al., 2014; Szczepankiewicz et al., 2019, 2016, 2015; Topgaard, 2019, 2017; Westin et al.,

2016); however, the strong model assumptions underlying QTE may limit the confidence in the specificity and the interpretation of the extracted kurtosis sources (Jespersen et al., 2019). Specifically, time dependent and intra-compartmental diffusional kurtosis effects may corrupt QTE estimates in a non-linear way. Other recent studies have attempted to decouple sources of inter- and intra-compartmental kurtosis using a frequency modulation of specific symmetrized DDE  experiments  (Ji  et  al.,  2019;  Paulsen  et  al.,  2015).  Nevertheless,  that  approach  is confounded by orientation dispersion (Paulsen et al., 2015), i.e., the frequency modulation in these  symmetrized  DDE  experiments  conflates  the  different  sources  of  kurtosis:  intracompartmental kurtosis, anisotropic kurtosis, and orientation dispersion.

## 4.1. CTI theoretical potentials and experimental requirements

In this study, we sought to develop a methodology capable of resolving anisotropic kurtosis ( 𝐾𝑎𝑛𝑖𝑠𝑜 ) , isotropic kurtosis ( 𝐾𝑖𝑠𝑜 ), and intra-compartmental kurtosis ( 𝐾𝑖𝑛𝑡𝑟𝑎 ) for any number of non-exchanging tissue compartments (in the absence of flow) and without relying on the assumption of Gaussian (or multiple Gaussian) diffusion. Our correlation tensor imaging (CTI) approach allows the estimation of the correlation tensor Z , which is expressed in DDE but not SDE signals. We have demonstrated, as previously predicted, that the Z tensor can provide the sought-after information, provided  that the  long mixing time regime  is reached (Jespersen, 2012; Jespersen et al., 2013) so that the contributions of the Q and S tensors are approximately zero (mixing time approximation #1) and the Z tensor is approximately equal to 4∆ 𝟐 𝑪 (mixing time  approximation  #2).  Although  the  signals  arising  from  the Q and S tensors  can  be eliminated by acquiring DDE experiments with an inverted gradient direction (c.f. Eq. 20), reaching the long mixing time is still necessary to fulfil the mixing time approximation #2 such that the anisotropic and isotropic kurtosis sources can be accurately extracted from the Z tensor

(c.f. Eqs. 17 and 18). The long mixing time regime need not be an implicit assumption of CTI since this regime can be empirically identified by comparing DDE signals with parallel and antiparallel experiments  (Jespersen  and  Buhl,  2011;  Lawrenz  and  Finsterbusch,  2013; Ozarslan,  2009),  as  indeed  done  in  this  study.  Since  both  mixing  time  approximations  are associated  with  the  same characteristic length, observing that the signals from parallel and antiparallel experiments are identical (i.e., their ratios are near unity) is likely to be a good indication that both Q S , →0 and that 𝒁 ≈ 4∆ 𝟐 𝑪 . Our finding that the ratio between parallel and antiparallel experiments is close to unity in the entire brain (Fig. 3), already at a mixing time of 13 ms, is consistent with previous studies suggesting the long mixing time is reached in neural tissues rather rapidly (Henriques et al., 2019; Ianuş et al., 2018; Shemesh et al., 2012, 2011;  Shemesh  and  Cohen,  2011)  and  consistent  with  our  simulations  which  showed  that 𝐾𝑎𝑛𝑖𝑠𝑜 and 𝐾𝑖𝑠𝑜 estimates are close to their ground truth values (at least for 𝑏𝑚 → 0 ).

To map the correlation tensor directly, typical DDE experiments with |𝒒1 | = |𝒒 2 | are insufficient. Although such experiments can resolve 𝐾𝑎𝑛𝑖𝑠𝑜 and other microscopic anisotropy measures,  e.g. (Ianuş  et  al.,  2018;  Jespersen  et  al.,  2013) , asymmetric  DDE  wavevector magnitudes  are   necessary for decoupling the  full  set  of  elements in the Z tensor  from  the W tensor (c.f.  Eq.  (24))  and for  measuring 𝐾𝑖𝑠𝑜 and 𝐾𝑖𝑛𝑡𝑟𝑎 . An  efficient  way  to  construct  an asymmetric  DDE  protocol  is  to  incorporate  measurements  whereby  one  of  the  DDE wavevectors is set to zero, effectively making the acquisition scheme a combination of SDE and DDE measurements (linear and planar encoding). However, we note that setting one of the wavevectors to zero is not necassary and in other applications, finite yet unequal magnitudes may be beneficial. In future studies, the CTI acquisition protocol could be further optimized by exploring  other  b-value  and  gradient  direction combinations potentially  better  suited for extraction of different kurtosis sources or for accelerating acquisition times.

We note that although the mathematical framework of CTI was derived based on the assumption that tissue can be represented by a number of non-exchanging signal contributions, exchange across components can only decrease the kurtosis values measured here (Fieremans et al., 2010; Jensen et al., 2005; Jensen and Helpern, 2010; Ning et al., 2018).

## 4.2. Sources of non-Gaussian diffusion in brain tissues

Fitting  the ex  vivo dataset  with  heavily  sampled  b-values  yielded  robust  maps  of  kurtosis sources in the mouse brain (Fig. 4). The usefulness of CTI can be gleaned by noticing that 𝐾𝑇 is similarly high in white matter (white arrows) and in areas near the ventricles (grey arrow). Our CTI approach resolved different sources of kurtosis in these two areas: in white matter, 𝐾𝑎𝑛𝑖𝑠𝑜 dominates  the  high  excess-kurtosis,  while  in  the  areas  close  to  the  ventricles, 𝐾𝑖𝑠𝑜 underlies the high excess-kurtosis, due to a wide distribution of diffusivities arising from the large partial volume effects between tissue water and freely diffusing water in the cerebral ventricles. In deep white matter and grey matter regions, 𝐾𝑖𝑠𝑜 is the smaller kurtosis source, which is consistent with recent studies showing that isotropic diffusion encoding signals exhibit smaller  deviations  from  mono-exponential  decay  in  similar  areas  (Dhital  et  al.,  2018; Szczepankiewicz et al., 2015).

Theoretically, negative 𝐾𝑖𝑛𝑡𝑟𝑎 measurements would be expected for the concave signal decay arising from fully restricted diffusion, that in some circumstances can lead to diffusion diffraction  effects  at  high  b-values  (Avram  et  al.,  2004;  Callaghan,  1995;  Callaghan  et  al., 1991;  Cohen  and  Assaf,  2002).  On  the  other  hand,  positive 𝐾𝑖𝑛𝑡𝑟𝑎 (scalar)  values  may  be expected for convex signal decays arising from diffusion in compartments with variable crosssection  dimensions,  or  in  the  presence  of  sub-structures  (Dhital  et  al.,  2018;  Novikov  and Kiselev,  2010).  Vanishing 𝐾𝑖𝑛𝑡𝑟𝑎 could  be  expected  at  very  long  diffusion  time  regimes

(Jespersen et al., 2007; Novikov et al., 2018a). For systems comprising an ensemble of different compartment types,  the observed 𝐾𝑖𝑛𝑡𝑟𝑎 index  corresponds to a weighted sum of  all  intracompartmental kurtosis sources (c.f. Eq. 15), where compartments with higher diffusivities have higher weights (e.g., perhaps diffusion in extra-cellular spaces). In this study, CTI's 𝐾𝑖𝑛𝑡𝑟𝑎 measurements from both ex vivo and in vivo experiments consistently evidenced positive values ranging from ~0.3 to ~0.6. Although these positive values seem to suggest that positive intracompartmental  kurtosis  sources  prevail  in  brain  tissues,  it  is  important  to  stress  that  our simulations suggest that 𝐾𝑖𝑛𝑡𝑟𝑎 estimates measured with the current protocol are likely to be compromised by high order terms not accounted for by CTI ( vide infra ). Although our results in general revealed consistent positive values across animals (both ex vivo and in vivo ), it is currently difficult to infer more subtle details on how 𝐾𝑖𝑛𝑡𝑟𝑎 varies along different brain areas due to the small dimensions of the mice/rat brain structures and due to small differences in slice positioning and orientation across the animals, in addition to the bias from higher order terms. Further studies are thus required for asserting the nature of intra-compartmental kurtosis.

## 4.3. Trade-off between accuracy and precision

Like all kurtosis measurements, higher order terms can play an important role in the accuracy and precision of the estimated metrics (Chuhutin et al., 2017; Ianuş et al., 2018) . Since kurtosis is  formally  defined  at  b-value  =  0,  it  (and  any  other  metric  defined  from  the  cumulant expansion) is inherently biased when measured at any finite b-value. Similarly, a biased metric will be measured in the CTI framework. We thus sought to explore the impacts of higher order effects by examining CTI 's b-value dependence in a controlled way -namely, by varying 𝑏𝑚𝑎𝑥 in the minimal protocol (c.f. Fig. 5). The higher order terms are expected to have a larger impact on parameter accuracy at higher b-values; however, lower b-values might not provide sufficient

diffusion-weighting  for  precise  kurtosis  estimation.  Indeed,  our  results  show  that  kurtosis estimates are not qualitatively consistent across e x vivo mouse brain specimens for 𝑏𝑚𝑎𝑥 values lower than 1 ms/  m2, and thus, we suggest that a 𝑏𝑚𝑎𝑥 value between 1.5 and 2 ms/  m2 may provide, in practice, an optimal trade-off between precision and accuracy. Kurtosis estimates for higher 𝑏𝑚𝑎𝑥 vales (2.25 and 2.5 ms/  m2, c.f. Fig. 5) indicate that higher order terms might introduce negative biases in 𝐾𝑎𝑛𝑖𝑠𝑜 and positive biases in 𝐾𝑖𝑛𝑡𝑟𝑎 . 𝐾𝑖𝑠𝑜 estimates exhibited both positive and negative bias trends depending on the specific brain regions -a  feature which should be considered in future studies.

Our simulations confirmed the expected trends between CTI's accuracy and precision (c.f.,  Fig.  7,  Fig  8,  Fig.  S3,  and  Fig.  S4):  firstly,  negative  bias  trends  for  synthetic 𝐾𝑎𝑛𝑖𝑠𝑜 estimates and positive bias trends for synthetic 𝐾𝑖𝑛𝑡𝑟𝑎 estimates were observed for orientation dispersion and/or non-Gaussian effects from restricted compartments (Fig. 7B, Fig. 8A, and Fig. 8B); secondly, our simulations confirmed that the signs of 𝐾𝑖𝑠𝑜 bias can vary with different microstructural scenarios. Particularly, pronounced higher order effects in 𝐾𝑖𝑠𝑜 are expected in regions containing partial volume effects with free water, since free water components have a strong signal attenuation due to high diffusivity. Indeed, these pronounced higher order effects were confirmed by supplementary simulations in which a free water component was added (supplementary Fig. S4). The supplementary simulations incorporating synthetic noise not only confirmed CTI 's low precision for 𝑏𝑚𝑎𝑥 values lower than 1.5 ms/  m2 (for plausible diffusion scenarios in neural tissues, at least) but also showed that CTI's accuracy can also be affected by Rician noise biases (supplementary Fig. S3). Although Rician biases may not be relevant for our ex vivo results due the high SNR, these biases may explain the higher 𝐾𝑖𝑠𝑜 estimates in inferior  grey  matter  regions  of  the in  vivo data  which  presented  lower  SNRs  (Fig.6D,  red arrows).

Given their observed b-value dependence, and similarly to all other metrics derived from the cumulant expansion, CTI measures cannot be considered to be completely accurate in  comparison  to  their  ground  truth  values.  Nevertheless,  our  results  show  that  CTI  still provides a robust characterization of the isotropic and anisotropic kurtosis sources. Specifically, despite the imperfect accuracy, 𝐾𝑎𝑛𝑖𝑠𝑜 would still be sensitive to differences in microscopic anisotropy and clearly highlights regions of high anisotropy, while 𝐾𝑖𝑠𝑜 would still  be  sensitive  tensor  magnitude  variance  across  different  compartments . Note  that  this observation is expected to be correct even for voxels containing partial volume effects (e.g., with free water) in which higher order effects are expected to be highly pronounced - according to  our  simulations, 𝐾𝑖𝑠𝑜 in  voxels  containing  free  water  (Fig.  S4)  is  higher  than  the 𝐾𝑖𝑠𝑜 estimates of voxels with no added free water (Fig. 7). Thus, CTI is able to robustly characterize the correct trends in the investigated system, also in terms of 𝐾𝑖𝑠𝑜 estimates.

Regarding the 𝐾𝑖𝑛𝑡𝑟𝑎 estimates, our simulations confirm that this metric is sensitive to restricted diffusion. Particularly, Fig. 7A4 shows that, for a given 𝑏𝑚𝑎𝑥 and a fixed degree of orientation dispersion, lower 𝐾𝑖𝑛𝑡𝑟𝑎 estimates are still associated with bigger compartments. However,  our  simulations  also  revealed  that  the  specific  protocol  used  here  for  CTI  is susceptible to higher order effects, and may introduce significant biases in 𝐾𝑖𝑛𝑡𝑟𝑎 estimates (in some cases, higher than the ground truth by orders of magnitude). Since these biases depend on confounding effects such as degree of mesoscopic orientation dispersion (Fig. 7B4), higher order terms can compromise the specificity of 𝐾𝑖𝑛𝑡𝑟𝑎 extracted using the current CTI protocol. Despite this issue, it is important to stress that 𝐾𝑖𝑛𝑡𝑟𝑎 estimates reported in this study represent the first attempt of measuring intra-compartmental kurtosis separately from the anisotropic and isotropic  kurtosis  sources.  We  expect  that  our  study  will  motivate  further  developments  in CTI's framework to minimize confounding effects on 𝐾𝑖𝑛𝑡𝑟𝑎 estimates ( vide infra ).

## 4.4. Future CTI vistas and clinical relevance

Recent studies based on isotropic diffusion encoding strategies (e.g. QTE) are showing that resolving  different  kurtosis  sources  can  be  clinically  useful,  because  they  can  distinguish tissues with different underlying microstructural features such as different tumor types (Nilsson et  al.,  2020;  Szczepankiewicz et al.,  2016,  2015).  Since the  CTI  approach provides  a more general and less assumption-driven strategy to decouple these different kurtosis sources, it can potentially  provide a  more  robust  characterization  of such tumours  or of  other  disorders  in which microscopic features are important. For instance, CTI-driven 𝐾𝑎𝑛𝑖𝑠𝑜 may more robustly detect the higher microscopic anisotropy of meningioma cells in comparison to the expected low anisotropy of glioblastoma cells. Previous studies have shown that diffusional kurtosis is a highly sensitive marker in stroke patients with stroke (Cheung et al., 2012; Rudrapatna et al., 2014). The different kurtosis sources afforded by CTI can likely provide important insights on, for example, the different microscopic/mesoscopic mechanisms underlying stroke, which are currently still under debate (Budde and Frank, 2010; Moseley et al., 1990b; van der Toorn et al.,  1996). For basic research studies, resolving different kurtosis sources can be a valuable technique for decoupling microstructural features from mesoscopic confounding factors. For example, CTI can potentially be applied in studies of brain aging to decouple the effects of white and grey matter degeneration (likely detected by decreases of 𝐾𝑎𝑛𝑖𝑠𝑜 ) from the known increases of free water partial volume effects due to the gross morphological atrophy (likely detected by increases of 𝐾𝑖𝑠𝑜 ) (Henriques, 2018; Metzler-Baddeley et al., 2012). CTI-driven 𝐾𝑖𝑠𝑜 could also be sensitive to edema in neurodegenerative diseases, such as Parkinson diseases and Multiple Sclerosis, or in traumatic injury (Donkin and Vink, 2010; Gelfand et al., 2012; Lee et al., 2006). In addition to 𝐾𝑎𝑛𝑖𝑠𝑜 and 𝐾𝑖𝑠𝑜 estimates, intra-compartmental kurtosis may be relevant in the future for investigating pathologies or processes in which microstructural modulations occur (e.g. axonal injury). Moreover, measuring intra-compartmental kurtosis via

CTI could be of interest for validation of other diffusion MRI methods that rely on assumptions such  as    vanishing  intra-compartmental  kurtosis  or  time-independent  diffusion  in  tissues (Henriques et al., 2015; Jespersen et al., 2007; Novikov et al., 2019, 2018a; Szczepankiewicz et  al.,  2016,  2015;  Westin  et  al.,  2016;  Zhang  et  al.,  2012).  We  once  more  reiterate  the importance of attenuating the higher order effects for resolving more accurate 𝐾𝑖𝑛𝑡𝑟𝑎 values.

All  these  potential  clinical  and  basic  research  applications  motivate  the  future application of CTI to investigate porous systems at large and neural tissues in particular, as well as translation of CTI to clinical scanners. Although this study focused on CTI's proof -ofconcept, our first in vivo contrasts of the rat brain show that consistent kurtosis source maps can be obtained from living animals (c.f. Fig. 6). As done for previous microscopic anisotropy measurements  of  DDE  (Kerkelä  et  al.,  2019;  Yang  et  al.,  2018)  and  multi-dimensional diffusion encoding (Sjölund et al., 2015), the clinical feasibility of CTI can be further promoted in future studies by refining its acquisition parameters and finding more optimal trends between precision, accuracy, and acquisition time. Although in this study, DDE data was processed with minimal pre-processing steps (only sub-pixel realigning was performed to correct data from motion and signals drifts), the CTI methodology can also be further improved by incorporating state-of-the-art denoising and artefact suppression algorithms. Even if future studies suggest that the CTI methodology is not yet compatible with clinical scanning acquisition times, the translation of the more general CTI approach to clinical scanners could still be fundamental to validate and calibrate the faster QTE acquisitions under different experimental conditions.

## 4.5. Methodological Limitations and further developments

As any other technique based on the cumulant expansion of diffusion-weighted signal decays (Chuhutin et al., 2017; Ianuş et al., 2018) ,  CTI measures are biased by higher order terms.

Though their accuracy may be affected, CTI's 𝐾𝑎𝑛𝑖𝑠𝑜 and 𝐾𝑖𝑠𝑜 may still be very useful if they provide a precise characterization of the different diffusional kurtosis sources, as mentioned above. By contrast, 𝐾𝑖𝑛𝑡𝑟𝑎 estimates obtained using the current CTI protocol are more severely affected by higher order effects as mentioned repeatedly above and shown explicitly in (Fig. 7A4 and Fig.8A4), leading to potentially more confounds in its interpretation. Still, we expect that these higher order effects can be mitigated significantly. This can be achieved by, e.g.: 1) designing DDE b-value combinations that decrease the high order effects differences between the kurtosis tensor W and the covariance tensor C - since 𝐾𝑖𝑛𝑡𝑟𝑎 information is captured by the subtraction of these tensors, minimizing the high order effects differences of these tensors can effectively suppress the biases in 𝐾𝑖𝑛𝑡𝑟𝑎 estimates; 2) applying CTI after powder-averaging DDE signals, thereby effectively making them independent to tissue dispersion (Eriksson et al., 2015; Jespersen et al., 2013; Kaden et al., 2016; Lasič et al., 20 14) - although this strategy may not eliminate the biases of higher order terms, it can remove the confounding dependence on orientation dispersion on all CTI metrics;  3) incorporating parameters in the CTI equation to fit and remove the impact of higher order terms, such as previously done for microscopic anisotropy estimates from DDE powder-averaged signals (Ianuş et al., 2018) . Note that future studies should carefully assess the robustness of such bias minimization strategies since they might  also  affect  CTI's  precision.  For  instance,  previous  studies  DKI studies,  showed  that incorporating higher order terms in the cumulant expansion equation significantly compromises the precision of the lower order terms (Chuhutin et al., 2017; Kiselev, 2017). These mitigation strategies are currently being implemented and will be reported in due course.

## 5. Conclusion

We provide a first general framework for measuring the diffusion correlation tensor directly from  double  diffusion  encoded  signals  and  without  a-priori  assumptions.  The  ensuing Correlation  Tensor  Imaging  (CTI)  approach  was  here  theoretically  derived,  explored  via simulations and applied to characterize the different sources of diffusional kurtosis in rodent brains both ex vivo and in vivo . Our theory shows that CTI resolves kurtosis sources, and can provide quantitative indices related to the anisotropic and isotropic diffusion variances ( 𝐾𝑎𝑛𝑖𝑠𝑜 and 𝐾𝑖𝑠𝑜 )  without relying on the Gaussian diffusion assumption; in addition, CTI offers an index  sensitive  to  intra-compartmental  kurtosis  ( 𝐾𝑖𝑛𝑡𝑟𝑎 ) arising  from  restricted  diffusion. Although future CTI protocol optimization is still required for mitigating higher order term effects on 𝐾𝑖𝑛𝑡𝑟𝑎 estimates, our results suggest that separating kurtosis sources is a promising vista for imparting  specificity  on  kurtosis measurements,  without  relying  on  specific microstructural assumption or constraints. All these features augur well for future development and implementation of CTI measures for basic research and biomedical applications.

## Acknowledgments

This study was funded by the European Research Council (ERC) (agreement No. 679058). The authors acknowledge the vivarium of the Champalimaud Centre for the Unknow, a facility of CONGENTO which is a research infrastructure co-financed by Lisboa Regional Operational Programme (Lisboa 2020), under the PORTUGAL 2020 Partnership Agreement through the European Regional Development Fund (ERDF) and Fundacao para a cienecia  e tecnologia (Portugal), project LISBOA-01-0145-FEDER-022170. The authors also want to thank Prof Dr Valerij  G.  Kiselev  (Freiburg  University)  and  Mr.  Leevi  Kerkela  (UCL)  for  insightful discussions and suggestions, Ms. Teresa Serradas Duarte and Dr. Daniel Nunes for assistance in the preparation of the ex vivo mouse brain specimens, and Andrada Ianus ̧for sharing her previous DDE sequence implementations.

## References

Alexander, D.C., Barker, G.J., 2005. Optimal imaging parameters for fiber-orientation estimation in diffusion MRI. Neuroimage 27, 357-367. https://doi.org/10.1016/j.neuroimage.2005.04.008

Assaf, Y., Cohen, Y., 1998. Non-mono-exponential attenuation of water and N-acetyl aspartate signals due to diffusion in brain tissue. J. Magn. Reson. 131, 69-85. https://doi.org/10.1006/jmre.1997.1313

Assaf, Y., Freidlin, R.Z., Rohde, G.K., Basser, P.J., 2004. New modeling and experimental framework to characterize hindered and restricted water diffusion in brain white matter. Magn. Reson. Med. 52, 965-978. https://doi.org/10.1002/mrm.20274

Avram, L., Assaf, Y., Cohen, Y., 2004. The effect of rotational angle and experimental parameters on the diffraction patterns and micro-structural information obtained from qspace diffusion NMR: implication for diffusion in white matter fibers. J. Magn. Reson. 169, 30-38. https://doi.org/10.1016/J.JMR.2004.03.020

Basser, P.J., 1995. Inferring microstructural features and the physiological state of tissues from diffusion-weighted images. NMR Biomed. 8, 333-344. https://doi.org/10.1002/nbm.1940080707

Basser, P.J., Mattiello, J., LeBihan, D., 1994. MR diffusion tensor spectroscopy and imaging. Biophys. J. 66, 259-267. https://doi.org/10.1016/S0006-3495(94)80775-1

Berman, J., 2009. Diffusion MR Tractography As a Tool for Surgical Planning. Magn. Reson. Imaging Clin. N. Am. 17, 205-214. https://doi.org/10.1016/j.mric.2009.02.002 Budde, M.D., Frank, J.A., 2010. Neurite beading is sufficient to decrease the apparent diffusion coefficient after ischemic stroke. Proc. Natl. Acad. Sci. U. S. A. 107, 14472-7.

https://doi.org/10.1073/pnas.1004841107

Callaghan, P.T., 1997. A Simple Matrix Formalism for Spin Echo Analysis of Restricted Diffusion under Generalized Gradient Waveforms. J. Magn. Reson. 129, 74-84. https://doi.org/10.1006/jmre.1997.1233

Callaghan, P.T., 1995. Pulsed-Gradient Spin-Echo NMR for Planar, Cylindrical, and Spherical Pores under Conditions of Wall Relaxation. J. Magn. Reson. Ser. A 113, 5359. https://doi.org/10.1006/JMRA.1995.1055

Callaghan, P.T., Coy, A., MacGowan, D., Packer, K.J., Zelaya, F.O., 1991. Diffraction-like effects in NMR diffusion studies of fluids in porous solids. Nature 351, 467-469. https://doi.org/10.1038/351467a0

Callaghan, P.T., Jolley, K.W., Lelievre, J., 1979. Diffusion of water in the endosperm tissue of wheat grains as studied by pulsed field gradient nuclear magnetic resonance. Biophys. J. 28, 133-41. https://doi.org/10.1016/S0006-3495(79)85164-4

Callaghan, P.T., Komlosh, M.E., 2002. Locally anisotropic motion in a macroscopically isotropic system: displacement correlations measured using double pulsed gradient spinecho NMR. Magn. Reson. Chem. 40, S15-S19. https://doi.org/10.1002/mrc.1122

Catani, M., Howard, R.J., Pajevic, S., Jones, D.K., 2002. Virtual in vivo interactive dissection of white matter fasciculi in the human brain. Neuroimage 17, 77-94.

Cheung, J.S., Wang, E., Lo, E.H., Sun, P.Z., 2012. Stratification of heterogeneous diffusion MRI ischemic lesion with kurtosis imaging: evaluation of mean diffusion and kurtosis MRI mismatch in an animal model of transient focal ischemia. Stroke 43, 2252-4. https://doi.org/10.1161/STROKEAHA.112.661926

Chuhutin, A., Hansen, B., Jespersen, S.N., 2017. Precision and accuracy of diffusion kurtosis estimation and the influence of b-value selection. NMR Biomed. 30.

https://doi.org/10.1002/NBM.3777

Cohen, Y., Assaf, Y., 2002. Highb-value q-space analyzed diffusion-weighted MRS and MRI in neuronal tissues - a technical review. NMR Biomed. 15, 516-542. https://doi.org/10.1002/nbm.778

Cory, D.G., Garroway, A.N., Miller, J.B., 1990. Applications of spin transport as a probe of local geometry. Polym Prepr 31, 149.

de Almeida Martins, J.P., Topgaard, D., 2016. Two-Dimensional Correlation of Isotropic and Directional Diffusion Using NMR. Phys. Rev. Lett. 116. https://doi.org/10.1103/PhysRevLett.116.087601

De Santis, S., Drakesmith, M., Bells, S., Assaf, Y., Jones, D.K., 2014. Why diffusion tensor MRI does well only some of the time: Variance and covariance of white matter tissue microstructure attributes in the living human brain. Neuroimage 89, 35-44. https://doi.org/10.1016/j.neuroimage.2013.12.003

Dell'Acqua, F., Rizzo, G., Scifo, P., Clarke, R.A., Scotti, G., Fazio, F., 2007. A Model-Based Deconvolution Approach to Solve Fiber Crossing in Diffusion-Weighted MR Imaging. IEEE Trans. Biomed. Eng. 54, 462-472. https://doi.org/10.1109/TBME.2006.888830

Descoteaux, M., Deriche, R., Knosche, T.R., Anwander, A., 2009. Deterministic and Probabilistic Tractography Based on Complex Fibre Orientation Distributions. IEEE Trans. Med. Imaging 28, 269-286. https://doi.org/10.1109/TMI.2008.2004424

Dhital, B., Kellner, E., Kiselev, V.G., Reisert, M., 2018. The absence of restricted water pool in brain white matter. Neuroimage 182, 398-406. https://doi.org/10.1016/j.neuroimage.2017.10.051

Donkin, J.J., Vink, R., 2010. Mechanisms of cerebral edema in traumatic brain injury: therapeutic developments. Curr. Opin. Neurol. 23, 293-299.

https://doi.org/10.1097/WCO.0b013e328337f451

Drobnjak, I., Zhang, H., Hall, M.G., Alexander, D.C., 2011. The matrix formalism for generalised gradients with time-varying orientation in diffusion NMR. J. Magn. Reson. 210, 151-157. https://doi.org/10.1016/j.jmr.2011.02.022

Eriksson, S., Lasič, S., Nilsson, M., Westin, C.-F., Topgaard, D., 2015. NMR diffusionencoding with axial symmetry and variable anisotropy: Distinguishing between prolate and oblate microscopic diffusion tensors with unknown orientation distribution. J.

Chem. Phys. 142, 104201. https://doi.org/10.1063/1.4913502

Eriksson, S., Lasic, S., Topgaard, D., 2013. Isotropic diffusion weighting in PGSE NMR by magic-angle spinning of the q-vector. J. Magn. Reson. 226, 13-18. https://doi.org/10.1016/j.jmr.2012.10.015

Falangola, M.F., Jensen, J.H., Babb, J.S., Hu, C., Castellanos, F.X., Di Martino, A., Ferris, S.H., Helpern, J.A., 2008. Age-related non-Gaussian diffusion patterns in the prefrontal brain. J. Magn. Reson. Imaging 28, 1345-1350. https://doi.org/10.1002/jmri.21604

Fieremans, E., Benitez, A., Jensen, J.H., Falangola, M.F., Tabesh, A., Deardorff, R.L., Spampinato, M.V.S., Babb, J.S., Novikov, D.S., Ferris, S.H., Helpern, J.A., 2013. Novel White Matter Tract Integrity Metrics Sensitive to Alzheimer Disease Progression. Am.

- J. Neuroradiol. 34, 2105-2112. https://doi.org/10.3174/ajnr.A3553

Fieremans, E., Jensen, J.H., Helpern, J.A., 2011. White matter characterization with diffusional kurtosis imaging. Neuroimage 58, 177-188.

https://doi.org/10.1016/j.neuroimage.2011.06.006

Fieremans, E., Novikov, D.S., Jensen, J.H., Helpern, J.A., 2010. Monte Carlo study of a twocompartment exchange model of diffusion. NMR Biomed. 23, 711-24. https://doi.org/10.1002/nbm.1577

Gelfand, J.M., Nolan, R., Schwartz, D.M., Graves, J., Green, A.J., 2012. Microcystic macular oedema in multiple sclerosis is associated with disease severity. Brain 135, 1786-1793. https://doi.org/10.1093/brain/aws098

Gong, N.-J., Wong, C.-S., Chan, C.-C., Leung, L.-M., Chu, Y.-C., 2013. Correlations between microstructural alterations and severity of cognitive deficiency in Alzheimer's disease and mild cognitive impairment: a diffusional kurtosis imaging study. Magn. Reson. Imaging 31, 688-694. https://doi.org/10.1016/j.mri.2012.10.027

Grussu, F., Ianuş, A., Tur, C., Prados, F., Schneider, T., Kaden, E., Ourselin, S., Drobnjak, I., Zhang, H., Alexander, D.C., Gandini Wheeler Kingshott, C.A.M., 2019. Relevance of -time dependence for clinically viable diffusion imaging of the spinal cord. Magn. -Reson. Med. 81, 1247-1264. https://doi.org/10.1002/mrm.27463

Guizar-Sicairos, M., Thurman, S.T., Fienup, J.R., 2008. Efficient subpixel image registration algorithms. Opt. Lett. 33, 156-8.

Hanisch, U.-K., Kettenmann, H., 2007. Microglia: active sensor and versatile effector cells in the normal and pathologic brain. Nat. Neurosci. 10, 1387-1394. https://doi.org/10.1038/nn1997

Hansen, B., Lund, T.E., Sangill, R., Jespersen, S.N., 2013. Experimentally and computationally fast method for estimation of a mean kurtosis. Magn. Reson. Med. 69, 1754-1760. https://doi.org/10.1002/mrm.24743

Hansen, B., Lund, T.E., Sangill, R., Stubbe, E., Finsterbusch, J., Jespersen, S.N., 2016. Experimental considerations for fast kurtosis imaging. Magn. Reson. Med. 76, 14551468. https://doi.org/10.1002/mrm.26055

Helpern, J.A., Adisetiyo, V., Falangola, M.F., Hu, C., Di Martino, A., Williams, K., Castellanos, F.X., Jensen, J.H., 2011. Preliminary evidence of altered gray and white

matter microstructural development in the frontal lobe of adolescents with attentiondeficit hyperactivity disorder: a diffusional kurtosis imaging study. J. Magn. Reson. Imaging 33, 17-23. https://doi.org/10.1002/jmri.22397

Henriques, R.N., 2018. Advanced Methods for Diffusion MRI Data Analysis and their Application to the Healthy Ageing Brain. University of Cambridge. https://doi.org/https://doi.org/10.17863/CAM.29356

Henriques, R.N., Correia, M.M., Nunes, R.G., Ferreira, H.A., 2015. Exploring the 3D geometry of the diffusion kurtosis tensor-Impact on the development of robust tractography procedures and novel biomarkers. Neuroimage 111. https://doi.org/10.1016/j.neuroimage.2015.02.004

Henriques, R.N., Jespersen, S.N., Shemesh, N., 2019. Microscopic anisotropy misestimation in spherical-mean single diffusion encoding MRI. Magn. Reson. Med. https://doi.org/10.1002/mrm.27606

Hui, E.S., Fieremans, E., Jensen, J.H., Tabesh, A., Feng, W., Bonilha, L., Spampinato, M. V., Adams, R., Helpern, J.A., 2012. Stroke Assessment With Diffusional Kurtosis Imaging. Stroke 43, 2968-2973. https://doi.org/10.1161/STROKEAHA.112.657742

Hui, E.S., Jensen, J.H., 2015. Double-pulsed diffusional kurtosis imaging for the in vivo assessment of human brain microstructure. Neuroimage 120, 371-381. https://doi.org/10.1016/j.neuroimage.2015.07.013

Ianuş, A., Jespersen, S.N., Serradas Duarte, T., Alexander, D.C., Drobnjak, I., Shemesh, N., 2018. Accurate estimation of microscopic diffusion anisotropy and its time dependence in the mouse brain. Neuroimage 183, 934-949.

https://doi.org/10.1016/j.neuroimage.2018.08.034

Jelescu, I.O., Budde, M.D., 2017. Design and Validation of Diffusion MRI Models of White

Matter. Front. Phys. 5, 61. https://doi.org/10.3389/fphy.2017.00061 parameter estimation for multi-compartmental diffusion in neuronal tissue. NMR

Jelescu, I.O., Veraart, J., Fieremans, E., Novikov, D.S., 2016. Degeneracy in model Biomed. 29, 33-47. https://doi.org/10.1002/nbm.3450

Jensen, J.H., Helpern, J.A., 2010. MRI quantification of non-Gaussian water diffusion by kurtosis analysis. NMR Biomed. 23, 698-710. https://doi.org/10.1002/nbm.1518 Jensen, J.H., Helpern, J.A., Ramani, A., Lu, H., Kaczynski, K., 2005. Diffusional kurtosis imaging: the quantification of non-gaussian water diffusion by means of magnetic resonance imaging. Magn. Reson. Imaging 53. https://doi.org/10.1002/mrm.20508 Jensen, J.H., Hui, E.S., Helpern, J.A., 2014. Double-pulsed diffusional kurtosis imaging.

NMR Biomed. 27, 363-370. https://doi.org/10.1002/nbm.3094

Jespersen, S.N., 2018. White matter biomarkers from diffusion MRI. J. Magn. Reson. 291, 127-140. https://doi.org/10.1016/j.jmr.2018.03.001

Jespersen, S.N., 2012. Equivalence of double and single wave vector diffusion contrast at low diffusion weighting. NMR Biomed. 25, 813-818. https://doi.org/10.1002/nbm.1808 Jespersen, S.N., Buhl, N., 2011. The displacement correlation tensor: Microstructure, ensemble anisotropy and curving fibers. J. Magn. Reson. 208, 34-43. https://doi.org/10.1016/j.jmr.2010.10.003

Jespersen, S.N., Kroenke, C.D., Østergaard, L., Ackerman, J.J.H., Yablonskiy, D.A., 2007.

Modeling dendrite density from magnetic resonance diffusion measurements. Neuroimage 34, 1473-1486. https://doi.org/10.1016/j.neuroimage.2006.10.037 Jespersen, S.N., Lundell, H., Sønderby, C.K., Dyrby, T.B., 2013. Orientationally invariant metrics of apparent compartment eccentricity from double pulsed field gradient

diffusion experiments. NMR Biomed. 26, 1647-62. https://doi.org/10.1002/nbm.2999 Jespersen, S.N., Olesen, J.L., Hansen, B., Shemesh, N., 2018. Diffusion time dependence of microstructural parameters in fixed spinal cord. Neuroimage 182, 329-342. https://doi.org/10.1016/J.NEUROIMAGE.2017.08.039

- Jespersen, S.N., Olesen, J.L., Ianuş, A., Shemesh, N., 2019. Effects of nongaussian diffusion on 'isotropic diffusion' measurements: An ex-vivo microimaging and simulation study.
- J. Magn. Reson. 300, 84-94. https://doi.org/10.1016/J.JMR.2019.01.007
- Ji, Y., Paulsen, J., Zhou, I.Y., Lu, D., Machado, P., Qiu, B., Song, Y.-Q., Sun, P.Z., 2019. In vivo microscopic diffusional kurtosis imaging with symmetrized double diffusion encoding EPI. Magn. Reson. Med. 81, 533-541. https://doi.org/10.1002/mrm.27419 Jones, D.K., 2008. Studying connections in the living human brain with diffusion MRI. Cortex 44, 936-952. https://doi.org/10.1016/j.cortex.2008.05.002

Jones, D.K., Horsfield, M.A., Simmons, A., 1999. Optimal strategies for measuring diffusion in anisotropic systems by magnetic resonance imaging. Magn. Reson. Med. 42, 515525. https://doi.org/10.1002/(SICI)1522-2594(199909)42:3&lt;515::AID-

MRM14&gt;3.0.CO;2-Q

Jones, D.K., Knösche, T.R., Turner, R., 2013. White matter integrity, fiber count, and other fallacies: The do's and don'ts of diffusion MRI. Neuroimage 73, 239-254. https://doi.org/10.1016/j.neuroimage.2012.06.081

Kaden, E., Kruggel, F., Alexander, D.C., 2016. Quantitative mapping of the per-axon diffusion coefficients in brain white matter. Magn. Reson. Med. 75, 1752-1763. https://doi.org/10.1002/mrm.25734

Kärger, J., 1985. NMR self-diffusion studies in heterogeneous systems. Adv. Colloid Interface Sci. 23, 129-148. https://doi.org/10.1016/0001-8686(85)80018-X

Kerkelä, L., Henriques, R.N., Hall, M.G., Clark, C.A., Shemesh, N., 2019. Validation and noise robustness assessment of microscopic anisotropy estimation with clinically feasible double diffusion encoding MRI. Magn. Reson. Med. mrm.28048. https://doi.org/10.1002/mrm.28048

Kiselev, V.G., 2017. Fundamentals of diffusion MRI physics. NMR Biomed. 30, e3602. https://doi.org/10.1002/nbm.3602

Koch, M.A., Finsterbusch, J., 2008. Compartment size estimation with double wave vector diffusion-weighted imaging. Magn. Reson. Med. 60, 90-101. https://doi.org/10.1002/mrm.21514

Kroenke, C.D., Ackerman, J.J.H., Yablonskiy, D.A., 2004. On the nature of the NAA diffusion attenuated MR signal in the central nervous system. Magn. Reson. Med. 52, 1052-1059. https://doi.org/10.1002/mrm.20260

Lampinen, B., Szczepankiewicz, F., Mårtensson, J., van Westen, D., Sundgren, P.C., Nilsson, M., 2017. Neurite density imaging versus imaging of microscopic anisotropy in diffusion MRI: A model comparison using spherical tensor encoding. Neuroimage 147, 517-531. https://doi.org/10.1016/j.neuroimage.2016.11.053

Lampinen, B., Szczepankiewicz, F., Novén, M., Westen, D., Hansson, O., Englund, E., Mårtensson, J., Westin, C., Nilsson, M., 2019. Searching for the neurite density with diffusion MRI: Challenges for biophysical modeling. Hum. Brain Mapp. 40, 2529-2545. https://doi.org/10.1002/hbm.24542

Lasič, S., Szczepankiewicz, F., Eriksson, S., Nilsson, M., Topgaard, D., 2014. Microanisotropy imaging: quantification of microscopic diffusion anisotropy and orientational order parameter by diffusion MRI with magic-angle spinning of the qvector. Front. Phys. 2, 11. https://doi.org/10.3389/fphy.2014.00011

Lawrenz, M., Finsterbusch, J., 2013. Double-wave-vector diffusion-weighted imaging reveals microscopic diffusion anisotropy in the living human brain. Magn. Reson. Med. 69, 1072-1082. https://doi.org/10.1002/mrm.24347

Le Bihan, D., Breton, E., 1985. Imagerie de diffusion in-vivo par résonance magnétique nucléaire. Comptes-Rendus l'Academie des Sci. 93, 27-34.

Le Bihan, D., Breton, E., Lallemand, D., Grenier, P., Cabanis, E., Laval-Jeantet, M., 1986. MR imaging of intravoxel incoherent motions: application to diffusion and perfusion in neurologic disorders. Radiology 161, 401-407.

https://doi.org/10.1148/radiology.161.2.3763909

Le Bihan, D., Johansen-Berg, H., 2012. Diffusion MRI at 25: exploring brain tissue structure and function. Neuroimage 61, 324-41.

https://doi.org/10.1016/j.neuroimage.2011.11.006

Lee, H.H., Fieremans, E., Novikov, D., 2019. Exploring the effect of varying axonal shape on the transverse diffusion inside EM-reconstructed axons using 3d Monte Carlo simulations, in: ISMRM. Montreal, p. 0650.

Lee, P.H., Shin, D.H., Kim, J.W., Song, Y.S., Kim, H.S., 2006. Parkinsonism with basal ganglia lesions in a patient with uremia: Evidence of vasogenic edema. Parkinsonism Relat. Disord. 12, 93-96. https://doi.org/10.1016/j.parkreldis.2005.07.009

Lin, L., Bhawana, R., Xue, Y., Duan, Q., Jiang, R., Chen, H., Chen, X., Sun, B., Lin, H., 2018. Comparative Analysis of Diffusional Kurtosis Imaging, Diffusion Tensor Imaging, and Diffusion-Weighted Imaging in Grading and Assessing Cellular Proliferation of Meningiomas. Am. J. Neuroradiol. https://doi.org/10.3174/AJNR.A5662

Metzler-Baddeley, C., O'Sullivan, M.J., Bells, S., Pasternak, O., Jones, D.K., 2012. How and how not to correct for CSF-contamination in diffusion MRI. Neuroimage 59, 1394-

1403. https://doi.org/10.1016/j.neuroimage.2011.08.043

Mitra, P.P., 1995. Multiple wave-vector extensions of the NMR pulsed-field-gradient spinecho diffusion measurement. Phys. Rev. B 51, 15074-15078. https://doi.org/10.1103/PhysRevB.51.15074

Mori, S., Crain, B.J., Chacko, V.P., Van Zijl, P.C.M., 1999. Three-dimensional tracking of axonal projections in the brain by magnetic resonance imaging. Ann. Neurol. 45, 265269. https://doi.org/10.1002/1531-8249(199902)45:2&lt;265::AID-ANA21&gt;3.0.CO;2-3 Moseley, M.E., 2002. Diffusion tensor imaging and aging - a review. NMR Biomed. 15, 553-560. https://doi.org/10.1002/nbm.785

Moseley, M.E., Cohen, Y., Kucharczyk, J., Mintorovitch, J., Asgari, H.S., Wendland, M.F., Tsuruda, J., Norman, D., 1990a. Diffusion-weighted MR imaging of anisotropic water diffusion in cat central nervous system. Radiology 176, 439-445. https://doi.org/10.1148/radiology.176.2.2367658

Moseley, M.E., Cohen, Y., Mintorovitch, J., Chileuitt, L., Shimizu, H., Kucharczyk, J., Wendland, M.F., Weinstein, P.R., 1990b. Early detection of regional cerebral ischemia in cats: Comparison of diffusion- and T2-weighted MRI and spectroscopy. Magn. Reson. Med. 14, 330-346. https://doi.org/10.1002/mrm.1910140218

Mulkern, R. V., Gudbjartsson, H., Westin, C.-F., Zengingonul, H.P., Gartner, W., Guttmann, C.R.G., Robertson, R.L., Kyriakos, W., Schwartz, R., Holtzman, D., Jolesz, F.A., Maier, S.E., 1999. Multi-component apparent diffusion coefficients in human brain. NMR Biomed. 12, 51-62. https://doi.org/10.1002/(SICI)1099-1492(199902)12:1&lt;51::AIDNBM546&gt;3.0.CO;2-E

Neil, J.J., Shiran, S.I., McKinstry, R.C., Schefft, G.L., Snyder, A.Z., Almli, C.R., Akbudak, E., Aronovitz, J.A., Miller, J.P., Lee, B.C., Conturo, T.E., 1998. Normal brain in human

newborns: apparent diffusion coefficient and diffusion anisotropy measured by using diffusion tensor MR imaging. Radiology 209, 57-66. https://doi.org/10.1148/radiology.209.1.9769812

Nilsson, M., Szczepankiewicz, F., Brabec, J., Taylor, M., Westin, C.-F., Golby, A., van Westen, D., Sundgren, P.C., 2020. Tensor-valued diffusion MRI in under 3 minutes: An initial survey of microscopic anisotropy and tissue heterogeneity in intracranial tumors. Magn. Reson. Imaging 83, 608-620.

Nilsson, M., van Westen, D., Ståhlberg, F., Sundgren, P.C., Lätt, J., 2013. The role of tissue microstructure and water exchange in biophysical modelling of diffusion in white matter. Magn. Reson. Mater. Physics, Biol. Med. 26, 345-370.

https://doi.org/10.1007/s10334-013-0371-x

Ning, L., Nilsson, M., Lasič, S., Westin, C.-F., Rathi, Y., 2018. Cumulant expansions for measuring water exchange using diffusion MRI. J. Chem. Phys. 148, 074109. https://doi.org/10.1063/1.5014044

Novikov, D.S., Fieremans, E., Jespersen, S.N., Kiselev, V.G., 2019. Quantifying brain microstructure with diffusion MRI: Theory and parameter estimation. NMR Biomed. 32, e3998. https://doi.org/10.1002/nbm.3998

Novikov, D.S., Kiselev, V.G., 2010. Effective medium theory of a diffusion-weighted signal.

NMR Biomed. 23, 682-697. https://doi.org/10.1002/nbm.1584

Novikov, D.S., Kiselev, V.G., Jespersen, S.N., 2018a. On modeling. Magn. Reson. Med. 79, 3172-3193. https://doi.org/10.1002/mrm.27101

Novikov, D.S., Veraart, J., Jelescu, I.O., Fieremans, E., 2018b. Rotationally-invariant mapping of scalar and orientational metrics of neuronal microstructure with diffusion MRI. Neuroimage 174, 518-538. https://doi.org/10.1016/j.neuroimage.2018.03.006

Ozarslan, E., 2009. Compartment shape anisotropy (CSA) revealed by double pulsed field gradient MR. J. Magn. Reson. 199, 56-67. https://doi.org/10.1016/j.jmr.2009.04.002

Paulsen, J.L., Özarslan, E., Komlosh, M.E., Basser, P.J., Song, Y.-Q., 2015. Detecting compartmental non-Gaussian diffusion with symmetrized double-PFG MRI. NMR Biomed. 28, 1550-6. https://doi.org/10.1002/nbm.3363

Peinado, H., Zhang, H., Matei, I.R., Costa-Silva, B., Hoshino, A., Rodrigues, G., Psaila, B., Kaplan, R.N., Bromberg, J.F., Kang, Y., Bissell, M.J., Cox, T.R., Giaccia, A.J., Erler, J.T., Hiratsuka, S., Ghajar, C.M., Lyden, D., 2017. Pre-metastatic niches: organ-specific homes for metastases. Nat. Rev. Cancer 17, 302-317. https://doi.org/10.1038/nrc.2017.6 Pfefferbaum, A., Sullivan, E. V, Hedehus, M., Lim, K.O., Adalsteinsson, E., Moseley, M.E., 2000. Age-related decline in brain white matter anisotropy measured with spatially corrected echo-planar diffusion tensor imaging. Magn. Reson. Med. 44, 259-68.

Pyatigorskaya, N., Le Bihan, D., Reynaud, O., Ciobanu, L., 2014. Relationship between the diffusion time and the diffusion MRI signal observed at 17.2 tesla in the healthy rat brain cortex. Magn. Reson. Med. 72, 492-500. https://doi.org/10.1002/mrm.24921 Reith, W., Hasegawa, Y., Latour, L.L., Dardzinski, B.J., Sotak, C.H., Fisher, M., 1995. Multislice diffusion mapping for 3-D evolution of cerebral ischemia in a rat stroke model. Neurology 45, 172-7. https://doi.org/10.1212/wnl.45.1.172

Riffert, T.W., Schreiber, J., Anwander, A., Knösche, T.R., 2014. Beyond fractional Neuroimage 100, 176-191. https://doi.org/10.1016/j.neuroimage.2014.06.015

Rudrapatna, S.U., Wieloch, T., Beirup, K., Ruscher, K., Mol, W., Yanev, P., Leemans, A., sensitivity and specificity of detecting microstructural alterations in brain tissue anisotropy: Extraction of bundle-specific structural metrics from crossing fiber models. van der Toorn, A., Dijkhuizen, R.M., 2014. Can diffusion kurtosis imaging improve the

chronically after experimental stroke? Comparisons with diffusion tensor imaging and histology. Neuroimage 97, 363-373. https://doi.org/10.1016/j.neuroimage.2014.04.013

Shemesh, N., Adiri, T., Cohen, Y., 2011. Probing Microscopic Architecture of Opaque Heterogeneous Systems Using Double-Pulsed-Field-Gradient NMR. J. Am. Chem. Soc. 133, 6028-6035. https://doi.org/10.1021/ja200303h

Shemesh, N., Barazany, D., Sadan, O., Bar, L., Zur, Y., Barhum, Y., Sochen, N., Offen, D., Assaf, Y., Cohen, Y., 2012. Mapping apparent eccentricity and residual ensemble anisotropy in the gray matter using angular double-pulsed-field-gradient MRI. Magn. Reson. Med. 68, 794-806. https://doi.org/10.1002/mrm.23300

Shemesh, N., Cohen, Y., 2011. Microscopic and compartment shape anisotropies in gray and white matter revealed by angular bipolar double-PFG MR. Magn. Reson. Med. 65, 1216-1227. https://doi.org/10.1002/mrm.22738

Shemesh, N., Jespersen, S.N., Alexander, D.C., Cohen, Y., Drobnjak, I., Dyrby, T.B., Finsterbusch, J., Koch, M.A., Kuder, T., Laun, F., Lawrenz, M., Lundell, H., Mitra, P.P., Nilsson, M., Özarslan, E., Topgaard, D., Westin, C.-F., 2016. Conventions and nomenclature for double diffusion encoding NMR and MRI. Magn. Reson. Med. 75, 82-87. https://doi.org/10.1002/mrm.25901

Sjölund, J., Szczepankiewicz, F., Nilsson, M., Topgaard, D., Westin, C.-F., Knutsson, H., 2015. Constrained optimization of gradient waveforms for generalized diffusion encoding. J. Magn. Reson. 261. https://doi.org/10.1016/j.jmr.2015.10.012

Stanisz, G.J., Szafer, A., Wright, G.A., Henkelman, R.M., Szafer, A., 1997. An analytical model of restricted diffusion in bovine optic nerve. Magn. Reson. Med. 37, 103-111. https://doi.org/10.1002/mrm.1910370115

Stejskal, E.O., Tanner, J.E., 1965. Spin Diffusion Measurements: Spin Echoes in the

Presence of a Time Dependent Field Gradient, The Journal of Chemical Physics. -American Institute of Physics. https://doi.org/10.1063/1.1695690

Sukstanskii, A.L., Yablonskiy, D.A., 2002. Effects of Restricted Diffusion on MR Signal Formation. J. Magn. Reson. 157, 92-105. https://doi.org/10.1006/JMRE.2002.2582 Sun, K., Chen, X., Chai, W., Fei, X., Fu, C., Yan, X., Zhan, Y., Chen, K., Shen, K., Yan, F., 2015. Breast Cancer: Diffusion Kurtosis MR Imaging-Diagnostic Accuracy and Correlation with Clinical-Pathologic Factors. Radiology 277, 46-55.

https://doi.org/10.1148/radiol.15141625

Szczepankiewicz, F., Lasič, S., van Westen, D., Sundgren, P.C., Englund, E., Westin, C.-F., Ståhlberg, F., Lätt, J., Topgaard, D., Nilsson, M., 2015. Quantification of microscopic diffusion anisotropy disentangles effects of orientation dispersion from microstructure: applications in healthy volunteers and in brain tumors. Neuroimage 104, 241-52. https://doi.org/10.1016/j.neuroimage.2014.09.057

Szczepankiewicz, F., Sjö Lund, J., Ståhlberg, F., Lä Tt, J., Nilsson, M., 2019. Tensor-valued diffusion encoding for diffusional variance decomposition (DIVIDE): Technical feasibility in clinical MRI systems. PLoS One 14, e0214238. https://doi.org/10.1371/journal.pone.0214238

Szczepankiewicz, F., van Westen, D., Englund, E., Westin, C.-F., Ståhlberg, F., Lätt, J.,

Sundgren, P.C., Nilsson, M., 2016. The link between diffusion MRI and tumor heterogeneity: Mapping cell eccentricity and density by diffusional variance decomposition (DIVIDE). Neuroimage 142, 522-532. https://doi.org/10.1016/J.NEUROIMAGE.2016.07.038

Topgaard, D., 2019. Diffusion tensor distribution imaging. NMR Biomed. 32, e4066.

https://doi.org/10.1002/nbm.4066

Topgaard, D., 2017. Multidimensional diffusion MRI. J. Magn. Reson. 275, 98-113. https://doi.org/10.1016/J.JMR.2016.12.007

Topgaard, D., 2015. Isotropic diffusion weighting using a triple-stimulated echo pulse sequence with bipolar gradient pulse pairs. Microporous Mesoporous Mater. 205, 4851. https://doi.org/10.1016/J.MICROMESO.2014.08.023

Tournier, J.-D., Calamante, F., Connelly, A., 2007. Robust determination of the fibre spherical deconvolution. Neuroimage 35, 1459-1472.

orientation distribution in diffusion MRI: Non-negativity constrained super-resolved https://doi.org/10.1016/J.NEUROIMAGE.2007.02.016

Valette, J., Giraudeau, C., Marchadour, C., Djemai, B., Geffroy, F., Ghaly, M.A., Le Bihan, D., Hantraye, P., Lebon, V., Lethimonnier, F., 2012. A new sequence for single-shot diffusion-weighted NMR spectroscopy by the trace of the diffusion tensor. Magn. Reson. Med. 68, 1705-1712. https://doi.org/10.1002/mrm.24193

Valiullin, R., 2017. Diffusion NMR of confined systems : fluid transport in porous solids and heterogeneous materials. Royal Society of Chemistry, Cambridge.

van der Toorn, A., Syková, E., Dijkhuizen, R.M., Vořišek, I., Vargová, L., Škobisová, E., van Lookeren Campagne, M., Reese, T., Nicolay, K., 1996. Dynamic changes in water ADC, energy metabolism, extracellular space volume, and tortuosity in neonatal rat brain during global ischemia. Magn. Reson. Med. 36, 52-60.

https://doi.org/10.1002/mrm.1910360110

Wang, J.-J., Lin, W.-Y., Lu, C.-S., Weng, Y.-H., Ng, S.-H., Wang, C.-H., Liu, H.-L., Hsieh, R.-H., Wan, Y.-L., Wai, Y.-Y., 2011. Parkinson Disease: Diagnostic Utility of Diffusion Kurtosis Imaging. Radiology 261, 210-217. https://doi.org/10.1148/radiol.11102277 Watson, G.S., 1965. Equatorial Distributions on a Sphere. Biometrika 52, 193.

## https://doi.org/10.2307/2333824

Westin, C.-F., Knutsson, H., Pasternak, O., Szczepankiewicz, F., Özarslan, E., Van Westen,

- D., Mattisson, C., O'donnell, L.J., Kubicki, M., Topgaard, D., Nilsson, M., 2016. Qspace trajectory imaging for multidimensional diffusion MRI of the human brain.

Neuroimage 135, 345-362. https://doi.org/10.1016/j.neuroimage.2016.02.039

Wong, E.C., Cox, R.W., Song, A.W., 1995. Optimized isotropic diffusion weighting. Magn. Reson. Med. 34, 139-143. https://doi.org/10.1002/mrm.1910340202

Xu, T., Yu, X., Perlik, A.J., Tobin, W.F., Zweig, J.A., Tennant, K., Jones, T., Zuo, Y., 2009. Rapid formation and selective stabilization of synapses for enduring motor memories. Nature 462, 915-919. https://doi.org/10.1038/nature08389

Yablonskiy, D.A., Bretthorst, G.L., Ackerman, J.J.H., 2003. Statistical model for diffusion attenuated MR signal. Magn. Reson. Med. 50, 664-669. https://doi.org/10.1002/mrm.10578

Yablonskiy, D.A., Sukstanskii, A.L., 2010. Theoretical models of the diffusion weighted MR signal. NMR Biomed. 23, 661-681. https://doi.org/10.1002/nbm.1520

Yang, G., Tian, Q., Leuze, C., Wintermark, M., McNab, J.A., 2018. Double diffusion encoding MRI for the clinic. Magn. Reson. Med. 80, 507-520.

https://doi.org/10.1002/mrm.27043

Zhang, H., Schneider, T., Wheeler-Kingshott, C.A., Alexander, D.C., 2012. NODDI: Practical in vivo neurite orientation dispersion and density imaging of the human brain. Neuroimage 61, 1000-1016. https://doi.org/10.1016/j.neuroimage.2012.03.072

## Figures

Figure  1 -Illustration  of  inter-compartmental  and  intra-compartmental  kurtosis  sources . A) Kurtosis can emerge in the dMRI signal decay if replicas of a single diffusion tensor ( 𝑫 𝑐 ) are mesoscopically dispersed along multiple orientations. In this scenario, kurtosis emerges from the variance across the different eigenvalues of the individual diffusion tensor 𝑫 𝑐 , and is hence termed here "anisotropic kurtosis". B) Kurtosis can also arise from polydisperse tensor traces, without orientation dispersion. In this scenario, kurtosis can be fully determined by the diffusion variance across the different mean  diffusivities 𝐷 ̅ 𝑐 ; hence, it is termed isotropic kurtosis. C) Finally, kurtosis can be a result of non-Gaussian, restricted diffusion within individual compartments with reflecting barriers. This kurtosis source is termed here intra-compartmental kurtosis. Given that (A) and (B) represent ensemble properties (n.b., if the gaussian diffusion tensors comprising the ensemble were isolated, each individual component would not exhibit any kurtosis), they can be considered as "inter-compartmental" kurtosis sources; (C) emerges even in a single compartment. D) In realistic tissues, diffusional kurtosis can be expected to emerge from a combination of all the  different  sources  of  inter-  and  intra-compartmental  kurtosis  -  the  total kurtosis is the sum of the three above kurtosis sources.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00122v2-figure-80da8f3f-4adc-4abf-ac0d-0d1e3efb6704.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00122v2-figure-318d5c2d-d92d-45a8-844b-9e79a7f04d5a.jpg)

Figure 2 -Acquisition requirements for a minimal protocol designed for CTI: A) Parameters and waveform of  a  standard  DDE  pulse  sequence  ( ∆1 and ∆2 are  the  diffusion  gradient  separation  times, 𝛿 is  the  diffusion gradient pulse duration, and 𝜏𝑚 is the mixing time between the two diffusion encoding modules marked by the blue and red lines); B) The eight gradient intensity combinations used for the minimal protocol of CTI. These gradient intensity combinations are acquired for 117 directions pairs for 𝒒1 -𝒒2 directions ( Jespersen's 5 design + 45 parallel DDE experiments) and can be acquired for different 𝑏𝑚𝑎𝑥 values.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00122v2-figure-7eb43d6c-16b0-44fc-8d6c-667ae20c271c.jpg)

Figure 3 -Raw double diffusion encoding (DDE) data: A) Non-diffusion weighted images of the five coronal slices acquired  from both  mouse  brains;  regions  of  interest  (ROIs) were  manually  defined.  White  matter ROIs were drawn in corpus callosum white matter ( WM1r and WM1l ), internal capsule white matter ( WM2r and WM2l ), and external capsule white matter ( WM3r and WM3l ), while GM ROIs comprised grey matter of motor cortex ( GM1r and GM1l ), somatosensory cortex ( GM2r and GM2l ), visual cortex, and auditory cortex ( GM3r and WM3l ).  For reference,  ROIs were also drawn in the cerebral ventricles (large contribution of  free water diffusion, CVr and CVl ); B) Powder average DDE data decays at a total b-value = 1 ms/μm 2 for slice #4 of mouse specimen #1 and for the following experiments: B1) parallel DDE experiments; B2) anti-parallel DDE experiments; B3) perpendicular DDE experiments; B4) ratio between parallel and anti-parallel experiments; and B5) ratio between parallel and perpendicular experiments; C) Powder averaged DDE signal decays at a total bvalue = 3 ms/μm 2 for slice #4 of mouse specimen #1 and for: C1) parallel DDE experiments; C2) anti-parallel DDE experiments; C3) perpendicular DDE experiments; C4) ratio between parallel and anti-parallel experiments; and C5) ratio between parallel and perpendicular experiments; D) Powder averaged DDE signal decays at total bvalue = 5 ms/μm 2 for slice #4 of  mouse specimen #1 and for: D1) parallel DDE experiments; D2) anti-parallel DDE experiments; D3) perpendicular DDE experiments; D4) ratio between parallel and anti-parallel experiments; and D5) ratio between parallel and perpendicular experiments.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00122v2-figure-a4a52daa-79f0-4497-a129-2694f469f7b5.jpg)

Figure  4 -CTI  kurtosis  measures  for  all  five  slices  extracted  from  all  data  acquired  in  mouse  brain specimen #1. A) total kurtosis of powder-average signals; B) anisotropic kurtosis; C) isotropic kurtosis; and D) intra-compartmental kurtosis.  White arrows highlight white matter regions with high anisotropy; grey arrows point to regions contaminated by free water partial volume effects; and red arrows point towards regions corrupted by Gibbs ringing artefacts.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00122v2-figure-a83ada15-3aae-4095-91b8-2d57e3548f98.jpg)

Figure  5 -CTI  kurtosis  measures  extracted  from  minimal  protocols  for  different bmax values.  A) total kurtosis; B) anisotropic kurtosis; C) isotropic kurtosis; and D) intra-compartmental kurtosis. Parametric maps in the left of each panel are plotted for sub-protocols with bmax = 1, 1.5, 2, 2.5 ms/  m2 and for the mouse brain specimen #1 (A1, B1, C1, and D1) and mouse brain specimen #2 (A2, B2, C2, and D2); while the kurtosis mean and standard deviation across animals for seven ROIs are shown in panels A3, B3, C3, and D3. For a reference, the mean kurtosis estimates of the extensive sampled b-value protocol are plotted by the dotted lines in panels A3, B3, C3, and D3.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00122v2-figure-a9b4d121-ed7f-4159-ad07-12ae984cc69a.jpg)

Figure 6 -CTI kurtosis measures extracted from in vivo rat brain data. A) Non-diffusion weighted images for all three acquired slices of both rats (white and grey matter ROIs are manually defined in slice #2 of rat brain #1 for SNR estimation); B) maps of the total kurtosis for all three acquired slices in both rats; C) maps of the anisotropic kurtosis for all three acquired slices in both rats; D) maps of the isotropic kurtosis for all three acquired slices of both rats (magenta arrows point to areas where partial volume effects between tissue and free water of cerebral ventricles is high; white arrows point to white matter); and E) maps of the intra-compartmental kurtosis for all three acquired slices in both rats.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00122v2-figure-2f84f09c-b8a3-4b7d-aa25-bd5cf77719a3.jpg)

Figure 7 -CTI kurtosis measures for synthetic signals of environments containing two types of Gaussian components. A) Simulations performed based on two aligned Gaussian components (axial and radial diffusivities for the first component are 2 and 0  m2/ms, respectively, while the axial and radial diffusivities for the second compartment are 1.5 and 0.5  m2/ms, respectively) - total kurtosis, anisotropic kurtosis, isotropic kurtosis, and intra-compartmental kurtosis estimates are plotted as a function of bmax from panels A1 to A4, respectively. B) Simulations  performed  based  on  replicas  of  the  two  Gaussian  components  dispersing  at  different  degrees (dispersion angles of 10o and 20o are plotted with purple and green lines, while completely powder-averaged (p.a.) replicas  are  plotted  with  the  yellow  line)  -  total  kurtosis,  anisotropic  kurtosis,  isotropic  kurtosis,  and  intracompartmental kurtosis estimates are plotted as a function of bmax in panels B1 to B4, respectively. Ground truth values are marked by the black dashed lines. Note that individual kurtosis estimates in these curves were obtained using 8 q 1 and q 2 gradient intensity combinations which are fully defined given a single bmax value according to the minimal CTI protocol.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00122v2-figure-8cbb67f7-9557-4d03-99a0-3478a8f7f88f.jpg)

Figure 8 -CTI kurtosis metrics for synthetic signals of environments containing two types of Gaussian components and a spherical compartment in which restricted diffusion occurs. A) Simulations based on two coherently aligned Gaussian components and a restricted spherical compartment for different diameters (spherical diameters  of  5,  7.5,  and  10  m  are  plotted  in  purple,  green  and  yellow  lines,  respectively)  -  total  kurtosis, anisotropic kurtosis, isotropic kurtosis, and intra-compartmental kurtosis estimates are plotted as a function of bmax  from  panels  A1  to  A4. B) Simulations  performed  based  on  replicas  of  two  Gaussian  components  with different levels of dispersion and a restricting spherical compartment with diameter of 10  m (dispersion angles of 10o and 20o are plotted with purple and green lines, while completely powder-averaged (p.a.) replicas are plotted with the yellow line)- total kurtosis, anisotropic kurtosis, isotropic kurtosis, and intra-compartmental kurtosis estimates are plotted as a function of bmax in panels B1 to B4, respectively. Ground truth values are marked by dashed  lines.  For  these  simulations,  the  different  underlying  components  were  set  to  have  similar  signal contributions. Note that individual kurtosis estimates  in these curves are obtained using 8 q 1  and q 2  gradient intensity combinations which are fully defined given a single bmax value according to the minimal CTI protocol.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00122v2-figure-3fb10420-95fb-4749-b521-1a66060b6201.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00122v2-figure-e7982a51-d244-4327-b970-b29a90eba0ae.jpg)

Supplementary Figure S1 -Partial volume effects between tissue and cerebral ventricle as estimated from high-resolution data. A) All high-resolution structural images (slice thickness = 79.4  m) overlapping with the thicker diffusion-weighted image slice #3 (with thickness = 900  m) -in all 11 structural images the cerebral ventricles (CV) were manually outlined (blue outlines on each panel). B) Tissue volume fraction map overlaid in averaged high-resolution structural image -tissue volume fraction is estimated by first projecting the 11 CV ROIs to the averaged structural imaging and then dividing the number of non-overlaying ROIs at each voxel position by the total number of CV ROIs (i.e. 11) -these volume fraction estimates are only performed on voxels overlaid by at least one ROI. C) For comparison, isotropic kurtosis estimated from diffusion-weighted data of slice #3 is displayed. The volume fraction profile in panel B is consistent with the profiles of the high isotropic kurtosis in panel C. This suggests that high isotropic kurtosis values close to the cerebral ventricles are a consequence of the partial volume effects between tissue and free water tissue. This figure was produced for mouse specimen #2.

Supplementary Figure S2 -CTI kurtosis measures for all five slices extracted from all data acquired in mouse brain specimen #2. A) total  kurtosis  of  powder-average  signals; B) anisotropic  kurtosis; C) isotropic kurtosis; and D) intra-compartmental kurtosis.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00122v2-figure-0a3db7bc-f2bf-4754-ac80-c4400567bb0a.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00122v2-figure-2ce3ca85-3c87-4e5c-b2ac-cce5139d4beb.jpg)

bmax (mslum )

max

Supplementary Figure S3 - CTI kurtosis measures for simulated signals corrupted with synthetic Rician distributed  noise  (SNR=40)  for  two  tissue  dispersion  levels:  A) Simulations  of  two  types  of  Gaussian components dispersing at a dispersion angles of 10o- the median, 25th and 75th percentile values for total kurtosis, anisotropic kurtosis, isotropic kurtosis, and intra-compartmental kurtosis are plotted as a function of bmax in panels A1 to A4, respectively - solid colored lines represents the median values, the 25th and 75th percentiles are plotted by the grey lines, and the ground truth values are marked by the black dashed lines. B) Simulations of two types of randomly oriented Gaussian components (i.e. powder-averaged compartments) -the median values of total kurtosis,  anisotropic  kurtosis,  isotropic  kurtosis,  and  intra-compartmental  kurtosis  estimates  are  plotted  as  a function of bmax in panels B1 to B4, respectively. Note that: 1) the SNR of simulations was set according to the estimates obtained from the superior ROIs of the in vivo rat data; 2) Median and percentile values are computed for 10000 noise initializations; 3) individual kurtosis estimates for noise iteration are obtained using 8 q 1 and q 2 gradient intensity combinations which are fully defined given a single bmax value according to the minimal CTI protocol.

Supplementary Figure S4 -CTI kurtosis measures for synthetic signals containing a free water component:

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00122v2-figure-abb46fa6-f9c2-431f-bc2d-c5deaf9a76e9.jpg)

A) total  kurtosis  estimates; B) anisotropic  kurtosis  estimates; C) isotropic  kurtosis  estimates;  and D) intracompartmental  kurtosis  estimates. These  simulations  were  performed  based  on  replicas  of  two  Gaussian components (components 1 and 2) dispersing at different degrees in addition to a free water Gaussian component. For  these  simulations,  the  different  compartments  were  set  to  have  similar  signal  contributions  (i.e.  volume fractions  of  1/3),  axial  and  radial  diffusivities  for  the  first  component  are  2  and  0  m2/ms,  axial  and  radial diffusivities for the second compartment are 1.5 and 0.5  m2/ms, while the isotropic diffusivity for compartment 3 was set to 3  m2/ms. In each panel, results for the dispersion angles of 10o and 20o are plotted in purple and green lines, while completely powder-averaged (p.a.) replicas are plotted in the yellow line. Ground truth values are marked in black dashed lines. Note that individual kurtosis estimates on these curves are obtained using 8 q 1 and q 2 gradient intensity combinations which are fully defined given a single bmax value according to the minimal CTI protocol.
Rafael Neto Henriques Sune Nørhøj Jespersen Noam Shemesh
Correlation Tensor Magnetic Resonance Imaging
Diffusional Kurtosis Imaging (DKI) is a sensitive biomarker for
microstructure in health and disease. However, DKI is not specific to any
microstructural property since it may emerge from several different sources.
Q-space trajectory encoding has been proposed for decoupling isotropic from
anisotropic kurtosis. Still, this method assume that the system is comprised of
multiple Gaussian diffusion components. Here, we develop a more general
framework for resolving the underlying kurtosis sources. We introduce
Correlation Tensor MRI (CTI) - an approach harnessing the versatility of double
diffusion encoding (DDE) and capable of explicitly decoupling isotropic and
anisotropic kurtosis components from intra-compartmental kurtosis effects
arising from restricted diffusion. Additionally, CTI provides an index that is
potentially sensitive to intra-compartmental kurtosis. The theoretical
foundations of CTI, as well as the first proof-of-concept CTI ex vivo
experiments in mouse brains at a field of 16.4T, are presented. We find that
anisotropic and isotropic kurtosis can decouple microscopic anisotropy from
substantial partial volume effects between tissue and free water. Our
intra-compartmental kurtosis index exhibited positive values in both white and
grey matter tissues. Simulations in different microenvironments show, however,
that our current CTI protocol for estimating intra-compartmental kurtosis is
limited by higher-order terms that were not considered in this study. CTI
measurements were then extended to in vivo settings and used to map healthy rat
brains at 9.4T. These in vivo CTI results were consistent with our ex vivo
findings. Although future studies are still required to mitigate the
higher-order effects on the intra-compartmental kurtosis index, results show
that CTI's more general estimates of anisotropic and isotropic kurtosis
contributions are already ripe for future studies.