Macroscopic Structural and Connectome Mapping of the Mouse Brain Using Diffusion Magnetic Resonance Imaging    

Materials and Reagents

1. 5 ml syringe (Sigma-Aldrich, catalog number: Z683582-100EA)

2. Vacutainer safety-lock blood collection set (25 G × 3/4” × 12”, 0.5 × 19 × 305 mm, Becton Dickinson, catalog number: 367283)

3. Nylon Zip ties (4” and 8” in length, LECO plastics, part# L-4-18, L-8-50)

4. 50 ml conical tubes (Corning, catalog number: 352070)

5. BALB/cByJ mice (Jackson Laboratories, catalog number: 001026, 8-10 weeks old, males and females)

6. Chloral hydrate (Sigma, catalog number: 102425)

7. Heparin (Sigma, catalog number: H3393-50KU)

8. PBS (Corning, catalog number: 21-031)

9. Gadodiamide (Omniscan, CAS# 131410-48-5)

10. 10% Formalin solution (PolyScience, catalog number: 08279-20)

11. Perfluoropolyether (Fomblin^®, PerkinElmer LLC, CAS# 69991067-9, Sigma-Aldrich, catalog number: 317926)

Equipment

1. Tools for routine transcardiac perfusion in mice (peristaltic pump and tubing, sharp small scissors, blunt tweezers, 21 G infusion butterfly (Becton Dickinson, catalog number: 367281), large container for blood collection, top of an insulated foam box to pin the mouse, and 23 G needles.

2. Horizontal 7 Tesla (T) Magnetic Resonance (MR) system (Bruker Biospin, Billercia, MA, USA) or other high-field (7T or greater) MRI system

3. 4-channel receive only cryogenic probe (Bruker Biospin, Billerica, MA, USA)

4. 72 mm inner diameter volume transmit coil (Bruker Biospin, Billerica, MA, USA)

5. Animal holder for the cryogenic probe (Bruker Biospin, Billerica, MA, USA)

6. Vacuum and vacuum chamber (e.g., 1-gal)

Software

1. Paravision (PV 6.0.1 or later)

2. Matlab R2019b or later (www.mathworks.com)

3. DTIStudio (www.mristudio.org)

4. AMIRA (thermofisher.com, version 5.0 or later)

5. DiffeoMap (www.mristudio.org) or ANTs (http://stnava.github.io/ANTs/)

6. Mrtrix (www.mrtrix.org)

7. Graph theoretical network analysis toolbox (GRETNA) (www.nitrc.org/projects/gretna)

8. GraphPad Prism (Version 8.4.3 for Windows, GraphPad Software, La Jolla California USA) (www.graphpad.com)

Procedure

1. Ex-vivo brain sample preparation

   1. Anesthetize the mouse with chloral hydrate (intraperitoneal injection 100 mg/kg in sterile PBS).

   2. Transcardially perfuse the mouse with 35 ml of cold PBS/heparin (50 units/ml) solution followed by 35 ml of 10% formalin. The perfusion rate is approximately 12 ml/min for PBS and formalin, with good perfusion assessed by the liver changing color from dark red to brownish/grey and the animal carcass becoming stiff.

   3. Decapitate the mouse at the mid-cervical line (around C3-C4), making sure not to damage the spinal cord, and place the head in 50 ml 10% Formalin solution at 4°C for 24 h in a 50 ml conical tube.

   4. After 24 h post-fixation, replace the formalin with PBS. Samples can be stored at 4°C at this point until ready to be scanned.

   5. Replace the PBS solution with 50 ml of 2 mM gadodiamide solution in PBS.

   6. Store the sample at 4°C for one week for the gadodiamide to diffuse into the tissue.

   7. Trim the skin and muscle tissues but keep the skull and eyeballs intact (Figure 1A). Remove the mandible bone and the tongue.

   8. Place one brain in the barrel of a 5 ml syringe, with the nose facing the hub of the syringe, then place 2-3 small pieces of bent zip-tie at the bottom and back of the brain to properly fix the specimen within the syringe barrel (Figure 1A).

   9. Replacing the cap of the syringe with a loosely tied vacutainer, fill the syringe with perfluoropolyether (Fomblin^®), insert the plunger, flip the syringe so that the hub points upward, and remove the cap.

   10. Place the syringe with its hub pointing upward in a vacuum chamber for 30 min to remove air bubbles (Figure 1B).

   11. Remove the syringe from the vacuum chamber, push out the remaining air and PBS in the barrel, and seal the cap of the syringe by tightening the vacutainer (Figure 1C).

       
       

2. MR data acquisition

   1. Place the syringe horizontally in the animal holder for the cryogenic probe and adjust the sample position so that the dorsal part of the brain is as close to the cryogenic coil as possible to maximize sensitivity. Use tape to fix the syringe to the animal holder (Figure 1D).

   2. Insert the animal holder into the magnet under the cryogenic probe (Figure 1E).

   3. Acquire a pilot scan using the Bruker Localizer protocol. For any MRI studies, Localizer is the very first scan that acquires reference images of the subject in three orthogonal planes. The images of the resulting scan appear in the ‘geometry editor,’ where the first three viewports show the reference brain slices in axial, sagittal, and coronal orientations. Therefore, the Localizer provides a quick view of the specimen in the magnet (Figure 1F). Check whether the sample is in the most sensitive region of the cryogenic probe with no apparent tilt toward the left or right sides. Adjust the position of the subject and re-run the Localizer protocol prior to proceeding to the next step.

   4. Adjust the tuning and match of the cryogenic probe and acquire a map of the main magnetic field (B₀ field) over the entire sample.

   5. Use the Bruker MapShim procedure to adjust shimming currents to achieve a relatively homogeneous B₀ field. In brief, select the specific scan to calculate the shim. Then choose Map_shim from the setup tab and define the target volume of interest in cubic, cylinder, or ellipsoid shapes. Shift, resize or rotate the target volume in the geometry editor such that the volume covers the entire specimen. Run the scan to compute the optimum shim values in the target volume based on the B₀ map measured in the previous step.

   6. Acquire high angular resolution diffusion weighted imaging (HARDI) of the whole mouse brain using a modified 3D gradient and spin echo (GRASE) sequence (Wu et al., 2013) (an alternative is the 3D multi-shot diffusion weighted echo planar imaging (EPI) sequence provided by Bruker) and with the following imaging parameters:

      1. Echo time (TE)/repetition time (TR): 33/400 min.

      2. Number of non-diffusion weighted images (b₀s): 2.

      3. Number of diffusion weighted images (DWIs): 60, auto-generated by the sequence.

      4. b-value: 5,000 s/mm².

      5. Resolution: 100 µm isotropic.

   
   
   Figure 1. Preparation for MRI.

   A. Sample preparation: Remove the tissues outside of the skull carefully without damaging the eyeballs (top panel). Place the brain in a 5ml syringe with small pieces of zip-ties to fix its position (bottom panel). B-C. Remove air from the syringe: Connect the syringe to a loosely tied vacutainer filled with Fomblin^® (shown in C) and place in the vacuum chamber (shown in B). Remove the vacutainer and turn on the vacuum for 30 min to remove air bubbles. Push out the remaining air after vacuum and seal the top by tightening the vacutainer. D. Place the sample in a manufacturer-made sample holder designed for the cryogenic probe. D’. A zoom-in view of the sample. E. Insert the sample holder into the magnet bore of the magnet. E’. A zoom-in view of the holder (indicated by the white arrow). F. Three orthogonal plane images acquired using the Localizer protocol on a 7 Tesla Bruker preclinical MRI system.

Data analysis

1. Data pre-processing

   For each dataset, perform the following steps accordingly:

   1. Motion correction: Using DTIStudio (Jiang et al., 2006), align all DWIs to the average of b₀s to remove small sample displacements due to vibrations and B₀ field drift during the long scan (Figure 2A).

   2. Skull-stripping: Use AMIRA segmentation editor to remove non-brain tissues and define the subject specific whole brain mask (Figure 2B).

   3. Estimation of diffusion tensor: From the raw data, compute the tensor model (Mori and Zhang, 2006) within the respective brain mask using weighted linear least squares estimations as implemented in MRtrix (command: dwi2tensor) (Figure 2C) (Tournier et al., 2012).

   4. Computation of average DWIs and fractional anisotropy (FA): Compute the average DWI (aDWI) from 60 DWIs using Matlab, and calculate the FA map from the tensor using MRtrix (command: tensor2metric) (Basser et al., 1994; Tournier et al., 2012).

   5. Image registration and transfer of atlas labels into subject’s native space: Using DiffeoMap, normalize the aDWI and FA maps to an MRI-based atlas (Chuang et al., 2011; Arefin et al., 2019) via multi-channel (aDWI + FA) large deformation diffeomorphic metric mapping (LDDMM) (Ceritoglu et al., 2009) (Figure 3A). Next, transfer the structural labels (i.e., brain regions or nodes) to the subject’s-native space using the inverse mapping from LDDMM (Figure 3B, also see Note 1). If DiffeoMap is not available, ANTs (http://stnava.github.io/ANTs/) can be used instead (Figure 2D).

      
      
      Figure 2. Illustration of the data pre-processing steps.

      A. Motion correction using DTIStudio. Run Automatic Image Registration (circled) to align all diffusion weighted images (DWIs) to the non-diffusion-weighted image (b₀). B. Use the AMIRA segmentation editor to generate a binary mask (purple) for the brain. C. Schematic diagram of the steps to compute the tensor and FA from the rawdata using Mrtrix. D. Use DiffeoMap for image registration and transformation of atlas labels into subject’s native image space.

      
      

      
      Figure 3. Image registration pipeline.

      A. Co-registration of the MR data (subject) into group averaged mouse brain atlas template using multi-channel LDDMM. B. Transformation of structural labels from MRI-based atlas to subject’s native space.

      
      

2. Data post-processing

   Assessment of brain microstructural changes

   At first, compute the Jacobian determinant value for each voxel from the mapping between atlas and subject images generated by LDDMM, and conduct whole-brain voxel-based morphometric analysis in Matlab to identify local volumetric changes affected by rearing, sex, and their interaction (2 × 2 ANOVA, FDR corrected, α = 0.1, P < 0.0105, cluster size > 25 voxels). See the Matlab codes (source data 1) used to conduct 2 × 2 ANOVA (White et al., 2020). Then, similarly perform 2 × 2 ANOVA (FDR corrected, α = 0.1, P < 0.007, cluster size > 25 voxels) to examine the voxel-wise changes in FA (White et al., 2020). These analyses will provide unbiased overviews of morphometric changes due to rearing, sex, and rearing by sex interaction.

   
   Selection of brain regions (nodes) for structural connectivity assessment

   Identify nodes that show rearing-mediated volumetric and FA changes to investigate structural connectivity alterations between nodes, as well as modifications in the brain global and regional network properties (also see Note 2). These nodes will be identical for both left and right hemispheres.

   
   Assessment of brain structural connectivity using fiber tractography

   Upon pre-processing the data and selection of potential brain nodes, execute the following steps accordingly for each individual subject to map axonal projections between nodes using probabilistic fiber tractography in MRtrix:

   Step 1: From the pre-processed raw data, estimate the response function for spherical deconvolution (command: dwi2response) (Tournier et al., 2012, Tournier et al., 2013). Specify the algorithm name ‘tournier’ (other options: dhollander, manual, fa, msmt_5tt, tax), gradient table, brain mask, and the maximum harmonic degree (l_max = 6).

   Step 2: Estimate the whole brain fiber orientation distribution (FOD) map from the pre-processed raw data and respective response function (command: dwi2fod) (Tournier et al., 2007). Define the algorithm name ‘CSD,’ gradient table, and brain mask.

   Step 3. Generate the whole brain fiber tractogram from the FOD map (command: tckgen) (Tournier et al., 2009). Use the whole brain mask as the ‘seed region’ to enable tracking fibers throughout the brain for whole brain tractography (whole brain tractogram) (Figure 4A). Set the tractography method to probabilistic, the FOD amplitude cut-off to 0.05, the minimum length of the fiber to 3 mm, and the target number of the streamlines to be counted to 5 million.

   Step 4: For node-to-node tractography, the whole brain tractography in step 3 may not generate enough streamlines for small nodes (e.g., amygdala). Further increasing the total number of streamlines (> 5 million) may not resolve this issue but requires significant computational resources. In this case, extract the regions of interest (ROIs) from the atlas co-registered into the subject’s native space using Matlab. Next, define a specific node as ‘seed region’ to initiate the fiber tracking from and another node as ‘target’ to define the fiber termination point. Then use these two nodes to extract the streamlines connecting two nodes (seed and target) using the tckedit command (Figure 4B). Consider two nodes as ‘connected’ if there is at least one streamline terminating at the target node; otherwise, they are ‘not connected.’

   
   
   Figure 4. Fiber tractography pipeline.

   A. Estimation of mouse whole brain fiber tractogram from the fiber orientation distribution (FOD) map. Red, green, and blue colors represent the fiber projections in x, y, and z-axis, respectively. Five million fibers were generated from each subject; 100 K streamlines were extracted for better visualization of the brain structures. B. Extraction of fibers connecting two specific nodes (seed = amygdala and target = PFC).

   
   

   Generating brain structural connectome matrix

   Repeat step 4 to estimate the structural connections between all possible pairs of nodes (ignore intra-regional connectivity) for both hemispheres (Figure 5A). For example, for 14 nodes in one hemisphere, the total number of tractograms would be the number of nodes N = 14 multiplied by N-1, or 14 × 13 = 182. Finally, for M number of seed regions and N number of target regions, generate an M × N matrix individually for the left and the right hemispheres. Assign the seed and target regions in horizontal and vertical axis, respectively, so that each cell represents the number of streamlines connecting the corresponding seed and target nodes (Figure 5B). Consider the number of streamlines between nodes as a measure of the connection strength. Generate the connectome matrix for all subjects and name them according to the subject IDs.

   
   
   Figure 5. Generation of the mouse brain structural connectome.

   A. Extraction of fibers connecting seed and target nodes. B. Generation of structural connectome from the tractograms estimated from selected seed and target nodes. Blue cells correspond to the tractograms shown in A, and white cells indicate intra-regional connectivity (not counted). C. Use the GRETNA software to compute global and regional brain network properties. Panels on the left list all possible properties available for computation. Select the properties based on the study design and transfer them to the pipeline option on the right panel using the respective arrows. Load the connectome matrix for all subjects belonging to one group with specific group ID and then load for the next group with different ID. Specify the output folder to store the results and define the network configuration. Finally, hit the ‘Run’ button to start computation.

   
   

   Brain network properties analysis

   Use the Matlab based Graph theoretical network analysis toolbox (GRETNA) to compute the brain global and regional network properties (Wang et al., 2015). Perform the following steps accordingly for brain network-based analysis (Figure 5C):

   1. Create an individual data folder containing two sub-folders for left and right hemispheres for each group.

   2. Save the connectome matrices as ‘.mat’ files in the respective folders.

   3. Open GRETNA in Matlab and select ‘Network Analysis’ (GRETNA >> Network Analysis).

   4. In the ‘Brain Connectivity Matrix’ tab, load all connectivity matrices of one hemisphere from one group and assign the group ID. Do the same for the other group.

   5. Locate a directory for saving the results in the ‘output directory’ tab.

   6. Next, select network properties to be computed from the Global Network Metrics and Nodal and Modular Network Metrics tabs.

   7. For global brain network analysis, select ‘Global – Small-World (SW)’ and ‘Global – Efficiency (Geff).’ For regional network properties, select ‘Nodal – Clustering Coefficient (NCp),’ ‘Nodal – Efficiency (Neff),’ and ‘Nodal – Degree Centrality (Dcent).’ Other properties can be selected as per the study design or requirements.

   8. Configure the brain network in the ‘Network Analysis’ tab as follows:

      Parameters	Value
      Sign of matrix	Absolute
      Thresholding method	Network sparsity
      Threshold sequence	0.05, 0.1, 0.15 (or as per the study design)
      Network type	Weighted
      Random network number	1,000

      
      

   9. Recheck the loaded data and the network configuration. Hit the ‘Run’ button if everything looks alright. Computation time depends on the number of subjects, size of the connectome matrices, random network number, and the threshold sequence.

   10. Once the computation is done, results can be retrieved from the output directory. For further assistance, please refer to the following manual from Neuroimaging Tools and Resources Collaboratory (NITRC): https://www.nitrc.org/docman/view.php/668/2262/manual_v2.0.0.pdf.

   
   Statistical analysis of the estimated structural connectivity and brain network properties

   To investigate the effect of rearing and sex on brain structural connectivity and brain network properties, perform a two-way ANOVA with rearing condition (CTL or UPS) and sex as fixed factors, followed by post-hoc comparisons using Tukey’s HSD or Sidak’s test using GraphPad Prism.

Notes

1. It is very important to check whether structural labels were correctly transferred and show good agreement with the corresponding structures. We recommend refining the segmentation manually, slice by slice, along the axial orientation, forfeiting attention to the other two orientations as well as to the slices preceding and following if necessary.

2. Selection of nodes for brain network analysis is crucial. Using unbiased voxel-based analyses, identify only those nodes which show UPS-mediated volumetric and FA alterations and are highly connected based on the Allen Mouse Brain Connectivity Atlas (Oh et al., 2014). Furthermore, selected nodes should be non-overlapping, having a unique set of connections to other nodes, and well delineated using a standard parcellation scheme that is comparable across species (Kaiser, 2011).

Acknowledgments

This work was supported by NIMH R01MH119164 (AK and JZ), NIMH R01MH118332 (AK and JZ), and R01NS102904 (JZ). For the original research paper where this protocol has been used, see White et al. (2020).

Competing interests

The authors declare no conflict of interest.

Ethics

Animal experimentation: All studies were approved by the Institutional Animal Care and Use Committee (IACUC) at Yale University, protocol #2020-10981, and were conducted in accordance with the recommendations of the NIH Guide for the Care and the Use of Laboratory Animals.

References

1. Aggarwal, M., Mori, S., Shimogori, T., Blackshaw, S. and Zhang, J. (2010). Three-dimensional diffusion tensor microimaging for anatomical characterization of the mouse brain. Magn Reson Med 64(1): 249-261.
2. Dorr, A. E., Lerch, J. P., Spring, S., Kabani, N. and Henkelman, R. M. (2008). High resolution three-dimensional brain atlas using an average magnetic resonance image of 40 adult C57Bl/6J mice. Neuroimage 42(1): 60-9.
3. Badea, A., Ng, K. L., Anderson, R. J., Zhang, J., Miller, M. I. and O'Brien, R. J. (2019). Magnetic resonance imaging of mouse brain networks plasticity following motor learning. PLoS One 14(5): e0216596.
4. Basser, P. J., Mattiello, J. and LeBihan, D. (1994). MR diffusion tensor spectroscopy and imaging. Biophys J 66(1): 259-267.
5. Calabrese, E., Badea, A., Cofer, G., Qi, Y. and Johnson, G. A. (2015). A Diffusion MRI Tractography Connectome of the Mouse Brain and Comparison with Neuronal Tracer Data. Cereb Cortex 25(11): 4628-4637.
6. Ceritoglu, C., Oishi, K., Li, X., Chou, M. C., Younes, L., Albert, M., Lyketsos, C., van Zijl, P. C., Miller, M. I. and Mori, S. (2009). Multi-contrast large deformation diffeomorphic metric mapping for diffusion tensor imaging. Neuroimage 47(2): 618-627.
7. Chuang, N., Mori, S., Yamamoto, A., Jiang, H., Ye, X., Xu, X., Richards, L. J., Nathans, J., Miller, M. I., Toga, A. W., Sidman, R. L. and Zhang, J. (2011). An MRI-based atlas and database of the developing mouse brain. Neuroimage 54(1): 80-89.
8. Edwards, T. J., Fenlon, L. R., Dean, R. J., Bunt, J., Consortium, I. R. C., Sherr, E. H. and Richards, L. J. (2020). Altered structural connectivity networks in a mouse model of complete and partial dysgenesis of the corpus callosum. Neuroimage 217116868.
9. Feo, R. and Giove, F. (2019). Towards an efficient segmentation of small rodents brain: A short critical review. J Neurosci Methods 32382-89.
10. Helmstaedter, M., Briggman, K. L., Turaga, S. C., Jain, V., Seung, H. S. and Denk, W. (2013). Connectomic reconstruction of the inner plexiform layer in the mouse retina. Nature 500(7461): 168-174.
11. Jiang, H., van Zijl, P. C., Kim, J., Pearlson, G. D. and Mori, S. (2006). DtiStudio: resource program for diffusion tensor computation and fiber bundle tracking. Comput Methods Programs Biomed 81(2): 106-116.
12. Kaffman, A., White, J. D., Wei, L., Johnson, F. K. and Krystal, J. H. (2019). Enhancing the Utility of Preclinical Research in Neuropsychiatry Drug Development. Methods Mol Biol 20113-22.
13. Kaiser, M., (2011). A tutorial in connectome analysis: topological and spatial features of brain networks. Neuroimage 57(3): 892-907.
14. Le Bihan, D., (2003). Looking into the functional architecture of the brain with diffusion MRI. Nat Rev Neurosci 4(6): 469-80.
15. Lerch, J. P., van der Kouwe, A. J., Raznahan, A., Paus, T., Johansen-Berg, H., Miller, K. L., Smith, S. M., Fischl, B. and Sotiropoulos, S. N. (2017). Studying neuroanatomy using MRI. Nat Neurosci. 20(3): 314-326..
16. Moldrich, R. X., Pannek, K., Hoch, R., Rubenstein, J. L., Kurniawan, N. D. and Richards, L. J. (2010). Comparative mouse brain tractography of diffusion magnetic resonance imaging. Neuroimage 51(3): 1027-1036.
17. Mori, S. and van Zijl, P. C. (2002). Fiber tracking: principles and strategies- a technical review. NMR Biomed 15(7-8): 468-480.
18. Mori, S. and Zhang, J. (2006). Principles of diffusion tensor imaging and its applications to basic neuroscience research. Neuron 51(5): 527-539.
19. Muller, H. P., Roselli, F., Rasche, V. and Kassubek, J. (2020). Diffusion Tensor Imaging-Based Studies at the Group-Level Applied to Animal Models of Neurodegenerative Diseases. Front Neurosci 14734.
20. Novikov, D. S., (2021). The present and the future of microstructure MRI: From a paradigm shift to normal science. J Neurosci Methods 351108947.
21. Oh, S. W., Harris, J. A., Ng, L., Winslow, B., Cain, N., Mihalas, S., Wang, Q., Lau, C., Kuan, L. and Henry, A. M. (2014). A mesoscale connectome of the mouse brain. Nature 508(7495): 207-214.
22. Pallast, N., Wieters, F., Nill, M., Fink, G. R. and Aswendt, M. (2020). Graph theoretical quantification of white matter reorganization after cortical stroke in mice. Neuroimage 217116873.
23. Saleeba, C., Dempsey, B., Le, S., Goodchild, A. and McMullan, S. (2019). A Student's Guide to Neural Circuit Tracing. Front Neurosci 13897.
24. Takemura, S. Y., Bharioke, A., Lu, Z., Nern, A., Vitaladevuni, S., Rivlin, P. K., Katz, W. T., Olbris, D. J., Plaza, S. M. and Winston, P. (2013). A visual motion detection circuit suggested by Drosophila connectomics. Nature 500(7461): 175-181.
25. Toga, A. W., Clark, K. A., Thompson, P. M., Shattuck, D. W. and Van Horn, J. D. (2012). Mapping the human connectome. Neurosurgery 71(1): 1-5.
26. Tournier, J. D., Calamante, F. and Connelly, A. (2009). Improved probabilistic streamlines tractography by 2^nd order integration over fibre orientation distributions. Proc Intl Soc Mag Reson Med 181670.
27. Tournier, J. D., Calamante, F. and Connelly, A. (2012). MRtrix: Diffusion tractography in crossing fiber regions. Int J Imaging Syst Technol 22(1): 53-66.
28. Tournier, J. D., Calamante, F. and Connelly, A. (2007). Robust determination of the fibre orientation distribution in diffusion MRI: non-negativity constrained super-resolved spherical deconvolution. Neuroimage 35(4): 1459-1472.
29. Tournier, J. D., Calamante, F. and Connelly, A. (2007). Robust determination of the fibre orientation distribution in diffusion MRI: non-negativity constrained super-resolved spherical deconvolution. Neuroimage 35(4): 1459-1472.
30. Tournier, J. D., Calamante, F. and Connelly, A. (2013). Determination of the appropriate b value and number of gradient directions for high-angular-resolution diffusion-weighted imaging. NMR Biomed 26(12): 1775-1786.
31. Tuch, D. S., Reese, T. G., Wiegell, M. R., Makris, N., Belliveau, J. W. and Wedeen, V. J. (2002). High angular resolution diffusion imaging reveals intravoxel white matter fiber heterogeneity. Magn Reson Med 48(4): 577-582.
32. Van Essen, D. C., Smith, S. M., Barch, D. M., Behrens, T. E., Yacoub, E., Ugurbil, K. and Consortium, W. U. (2013). The WU-Minn Human Connectome Project: an overview. Neuroimage 8062-79.
33. Wang, J., Wang, X., Xia, M., Liao, X., Evans, A. and He, Y. (2015). GRETNA: a graph theoretical network analysis toolbox for imaging connectomics. Front Hum Neurosci 9386.
34. Wedeen, V. J., Wang, R. P., Schmahmann, J. D., Benner, T., Tseng, W. Y., Dai, G., Pandya, D. N., Hagmann, P., D'Arceuil, H. and de Crespigny, A. J. (2008). Diffusion spectrum magnetic resonance imaging(DSI) tractography of crossing fibers. Neuroimage 41(4): 1267-1277.
35. White, J. D., Arefin, T. M., Pugliese, A., Lee, C. H., Gassen, J., Zhang, J. and Kaffman, A. (2020). Early life stress causes sex-specific changes in adult fronto-limbic connectivity that differentially drive learning. Elife 9: e58301.
36. Wu, D., Reisinger, D., Xu, J., Fatemi, S. A., van Zijl, P. C., Mori, S. and Zhang, J. (2014). Localized diffusion magnetic resonance micro-imaging of the live mouse brain. Neuroimage 9112-20.
37. Wu, D., Xu, J., McMahon, M. T., van Zijl, P. C., Mori, S., Northington, F. J. and Zhang, J. (2013). In vivo high-resolution diffusion tensor imaging of the mouse brain. Neuroimage 8318-26.
38. Wu, D. and Zhang, J. (2016). In vivo mapping of macroscopic neuronal projections in the mouse hippocampus using high-resolution diffusion MRI. Neuroimage 12584-93.
39. Xiong, J., Ren, J., Luo, L. and Horowitz, M. (2018). Mapping Histological Slice Sequences to the Allen Mouse Brain Atlas Without 3D Reconstruction. Front Neuroinform 1293.

1. Arefin, T. M., Lee, C. H., White, J. D., Zhang, J. and Kaffman, A. (2021). Macroscopic Structural and Connectome Mapping of the Mouse Brain Using Diffusion Magnetic Resonance Imaging. Bio-protocol 11(22): e4221. DOI: 10.21769/BioProtoc.4221.
   
2. White, J. D., Arefin, T. M., Pugliese, A., Lee, C. H., Gassen, J., Zhang, J. and Kaffman, A. (2020). Early life stress causes sex-specific changes in adult fronto-limbic connectivity that differentially drive learning. Elife 9: e58301.