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 T
1-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=224x224mm
2,
imaging resolution=2x2x2mm
3, 65 slices, 30 independent diffusion-weighted
directions, b-value=1000 sec/mm
2. T
1-weighted
images
were acquired using MPRAGE sequence with imaging resolution=1x1x1mm
3,
160 slices.
Fiber tracing from a parcellated cortical gyrus: With T
1-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 50
o. Using
inferior parietal gyrus (IPG) as an example, the parcellated cortical ribbon
transformed from T
1-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=β0+β1t+β2t2+β3t3+ε
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