Global and regional cortical connectivity maturation index (CCMI) of developmental human brain with quantification of short-range association tracts
Minhui Ouyang1, Tina Jeon1, Jennifer Muller1, Virendra Mishra2, Haixiao Du3, Yu Wang3, Yun Peng4, Bo Hong5, and Hao Huang1,6

1Department of Radiology, Children's Hospital of Philadelphia, Philadelphia, PA, United States, 2Cleveland Clinic Lou Ruvo Center for Brain Health, Las Vegas, NV, United States, 3Department of Electronic Engineering, Tsinghua University, Beijing, China, People's Republic of, 4Department of Radiology, Beijing Children's Hospital, Capital Medical University, Beijing, China, People's Republic of, 5Department of Biomedical Engineering, School of Medicine, Tsinghua University, Beijing, China, People's Republic of, 6Department of Radiology, Perelman School of Medicine, University of Pennsylvania, Philadelphia, PA, United States

Synopsis

Disturbance of precisely balanced strengthening of certain axons and pruning of others in developmental human brains is associated with mental disorders such as autism and schizophrenia. To characterize this balance, we defined a cortical connectivity maturation index (CCMI) derived from short-range association tracts traced with diffusion MRI tractography. The brain CCMI values were measured with diffusion MRI and T1-weighted datasets of 21 healthy subjects with age of 2-25 years. CCMI in all cortical regions decreased in early developmental stage and increased later, yet with distinctive trajectories. The observed CCMI dynamics may be underlaid by heterogeneous pruning among cortical regions.

Purpose

From early childhood to adulthood, synaptogenesis and synaptic pruning reshape the architecture of neuronal connections in developmental human brains. Disturbance of precisely balanced strengthening of certain axons and pruning of others may cause mental disorders such as autism and schizophrenia [1-2]. To characterize this balance, we defined and measured a cortical connectivity maturation index (CCMI) derived from short-range association tracts with diffusion MRI (dMRI) tractoghraphy traced from each parcellated cortical gyrus with T1-weighted images. The goal is to develop a novel maturational index sensitive to connectivity of specific cortical regions in normal and pathological human brain development.

Methods

Participants: 21 healthy subjects with age of 2-25 years were included in present study. Acquisition of dMRI and T1-weighted image: All MR scans were performed on a 3T Philips Achieva MR system. dMRI were acquired using single-shot EPI with SENSE=2.3 and other parameters were: TR/TE=7960/83 ms, FOV=224x224mm2, imaging resolution=2x2x2mm3, 65 slices, 30 independent diffusion-weighted directions, b-value=1000 sec/mm2. T1-weighted images were acquired using MPRAGE sequence with imaging resolution=1x1x1mm3, 160 slices. Fiber tracing from a parcellated cortical gyrus: With T1-weighted image, the brain cortical surface was rendered and parcellated into 68 gyral labels [3] using Freesurfer (http://surfer.nmr.mgh.harvard.edu). Fiber assignment of continuous tractography (FACT) [4] was used to trace the whole brain fibers for all subjects in DiffusionToolkit (http://www.trackvis.org/dtk/) with FA threshold of 0.2 and angular threshold of 50o. Using inferior parietal gyrus (IPG) as an example, the parcellated cortical ribbon transformed from T1-weighted image space to dMRI space (Figure 1a) was then dilated by 8 mm (shown in green overlaid on b0 image in Figure 1b) with in-house program to get through the dense white matter (WM) zone for initiating fiber tracking [5]. All association fibers traced from IPG as seed ROI were shown in Figure 1c. Categorization of long- and short-range fibers based on termination location of the other end of fibers: The adjacent and non-adjacent gyral labels of each cortical gyrus were identified. Using IPG (shown in green in the 3D reconstructed brain) as an example (Figure 1d), its adjacent gyri are superior parietal gyrus(yellow), lateral occipital gyrus(red) and supra marginal gyrus(blue) and all other gyri are non-adjacent gyri for IPG. Then association fibers initiated from IPG can be categorized into short- and long-range based on the other end of fibers terminating in adjacent and non-adjacent gyri to IPG, respectively (Figure 1e). Regional CCMI of IPG was calculated as the ratio between number of short-range association fiber and total number of association fibers initiated from IPG. Global and regional CCMI developmental curve analysis: Global CCMI was calculated as the ratio between the total number of short-range association fibers and total number associate fibers from all 68 gyri. To investigate the developmental curve of CCMI, the following equation was used for fitting a quadratic or cubic model between y (global or regional CCMI) and age t, y=β01t+β2t23t3 using R software (https://www.r-project.org/), where β0, β1, β2 and β3 were parameters to be estimated (β3 is for cubic fitting only) and ε was an error term. The age at which lowest CCMI value of whole-brain and three representative cortical regions: primary somatosensory cortex (S1), visual cortex (V1) and prefrontal cortex were calculated.

Results

Figure 2 demonstrates the developmental curve of global and regional CCMI (fitted with quadratic and cubic model with a 95% confidence interval in shadow, respectively). Global CCMI decreases initially from 2 to 16 years of age, and then followed by a rise from 16 to 25 years of age (Figure 2a). In addition, we examined the developmental curve of regional CCMI in three representative functional regions of human brain: S1, V1 and prefrontal cortex (Figure 2b-2d). Similar as global CCMI, the developmental curves of regional CCMI from all three functional regions show initial decrease followed by a later rise. The developmental curves also suggested spatiotemperally heterogeneous CCMI dynamics. Specifically, the ages of minimum CCMI vary across these functional regions. The age of minimum CCMI for prefrontal cortex (~17 years) is later than those for S1 (~9 years) and V1 (~14 years), suggesting later maturation of higher-order functional regions than primary sensorimotor regions.

Discussion and conclusion

The CCMI is spatiotemporally heterogeneous, decreasing in early childhood and increasing later with minimum CCMI observed at various ages among cortical gyri. The proposed CCMI reflecting the balance between strengthening of certain axons and pruning of others may be used to characterize normal brain development and detect atypical development under pathological conditions. Acquisition and analysis of dMRI data from more subjects are under way.

Acknowledgements

This study is funded by NIH MH092535 and MH092535-S1.

References

[1]Paus et al (2008) Nature Rev Neurosci 9:947. [2]Courchesne and Pierce (2005) Curr Opin Neurobiol 15:225. [3]Desikan et al (2006) Neuroimage 31:968. [4]Mori et al (1999) Ann Neurol 45:265. [5] Reveley et al (2015) Proc Natl Acad Sci USA 112:E2802

Figures

Figure 1: The schematic pipeline of cortical parcellation (a), fiber tracing (b,c) and categorization of long- and short- range fibers (d, e) from a certain cortical gyrus.

Figure 2: Developmental curve of global CCMI fitted with quadratic model (a) and regional CCMI fitted with cubic model from three representative regions: primary somatosensory cortex (S1) (b), visual cortex (V1) (c) and prefrontal cortex (d), respectively. Each black circle in the scatter plot represents the CCMI from one subject.



Proc. Intl. Soc. Mag. Reson. Med. 24 (2016)
0560