Rongpin Wang
Molly Wilkinson3‡
Tara Kane
Emi Takahashi- 1Division of Newborn Medicine, Department of Medicine, Boston Children's Hospital, Harvard Medical School, Boston, MA, United States
- 2Department of Radiology, Guizhou Provincial People's Hospital, Guiyang, China
- 3Department of Behavioral Neuroscience, Northeastern University, Boston, MA, United States
There has been evidence that during brain development, emerging thalamocortical (TC) and corticothalamic (CT) pathways converge in some brain regions and follow each other's trajectories to their final destinations. Corpus callosal (CC) pathways also emerge at a similar developmental stage, and are known to converge with TC pathways in specific cortical regions in mature brains. Given the functional relationships between TC and CC pathways, anatomical convergence of the two pathways are likely important for their functional integration. However, it is unknown (1) where TC and CT subcortically converge in the human brain, and (2) where TC and CC converge in the cortex of the human brain, due to the limitations of non-invasive methods. The goals of this study were to describe the spatio-temporal relationships in the development of the TC/CT and CC pathways in the human brain, using high-angular resolution diffusion MR imaging (HARDI) tractography. Emerging cortical, TC and CC pathways were identified in postmortem fetal brains ranging from 17 gestational weeks (GW) to 30 GW, as well as in vivo 34–40 GW newborns. Some pathways from the thalami were found to be converged with pathways from the cerebral cortex as early as 17 GW. Such convergence was observed mainly in anterior and middle regions of the brain until 21 GW. At 22 GW and onwards, posterior pathways from the thalami also converged with cortical pathways. Many CC pathways reached the full length up to the cortical surface as early as 17 GW, while pathways linked to thalami (not only TC axons but also including pathways linked to thalamic neuronal migration) reached the cortical surface at and after 20 GW. These results suggest that CC pathways developed earlier than the TC pathways. The two pathways were widespread at early stages, but by 40 GW they condensed and formed groups of pathways that projected to specific regions of the cortex and overlapped in some brain regions. These results suggest that HARDI tractography has the potential to identify developing TC/CT and CC pathways with the timing and location of their convergence in fetal stages persisting in postnatal development.
Introduction
Thalamocortical/corticothalamic (TC/CT) and corpus callosal (CC) pathways are major white matter tracts in the brain which are conserved across species throughout evolution, and partake in fundamental brain functions. The TC/CT system is important for consciousness (Bagshaw et al., 2014), working memory (Sarnthein et al., 2005), sensory relay (Jones, 2007), executive control, and salience processing (Seeley et al., 2007). The CC pathway connects the left and right hemispheres and has multiple functional integrational roles in visual processing (Bocci et al., 2014), motor movement (Fabri and Polonara, 2013), and somatosensory processing (Fabri et al., 1999; Marcano-Reik and Blumberg, 2008; Magnuson et al., 2014).
During brain development, the TC axons go through a number of steps in order to reach their final destinations (Rakic, 1975, 1988), following a distinct pathway stemming from the thalamus and reaching into the cerebral cortex (López-Bendito and Molnar, 2003; Garel and Rubenstein, 2004). Around the second and third gestational weeks (GW) in humans, the neocortex and thalamus start to link with each other through specific and reciprocal connections (Leyva-Díaz and López-Bendito, 2013), and continue increasing the number of connections between them. It is still under debate how human TC pathways find a way to go through brain regions to finally project to the cortex.
The “handshake hypothesis” was proposed by Blakemore and Molnár (1990) to explain how TC pathways ascending through the internal capsule project to their cortical targets with assistance from reciprocal descending cortical pathways (Molnár et al., 1998). Métin and Godement (1996) used mice and hamsters and suggested that the ganglionic eminence (GE) may be an intermediate target in the pathfinding of axons between the cortex and the thalamus.
Anatomically, while the TC pathways project to a wide range of areas throughout the brain including visual, auditory, motor, and somatosensory areas, the CC pathways are known to project to only specific regions. For example, in the primary visual cortex, the CC pathways project specifically to the border of V1 and V2 areas, in mice (Wang et al., 2007) and monkeys (Dehay et al., 1986; Kennedy et al., 1986; Dehay and Kennedy, 1988), while the TC pathways project to more broad areas throughout the primary visual cortex, with most connecting to layer IV (Robertson, 1976; Killackey and Belford, 1979; Burton and Kopf, 1984; Afifi and Bergman, 1997; Van Essen, 2005; Viaene et al., 2011). Given the functional relationships between the TC and CC pathways, anatomical convergence of the two pathways at the regions where CC pathways are projecting should be important for their functional integration. However, it is still unknown where these two pathways converge in the human brain due to the limitation of spatial resolution of non-invasive methods. Moreover, their anatomic relationships in more anterior brain regions (other than the visual area) and their spatial relationships during fetal developmental stages are elusive, even in animal brains.
Previous diffusion tensor MR imaging (DTI) tractography studies have been done on the TC (e.g., Johansen-Berg et al., 2005; Yamada et al., 2007; Jaermann et al., 2008; Hong et al., 2010; Klein et al., 2010; Zhang et al., 2010) and CC tracts (e.g., Häberling et al., 2011; Luders et al., 2012; Whitehead et al., 2014), with regionally differential developmental courses (Gilmore et al., 2007). High-angular resolution diffusion MR imaging (HARDI) tractography enables identification of water diffusivity in many different directions in each imaging voxel, which theoretically provides an advantage over DTI (Frank, 2002; Tournier et al., 2004, 2007). HARDI theoretically allows reconstruction of complex crossing tissue coherence in the brain (Tuch et al., 2003) and even in immature fetal brains (e.g., Takahashi et al., 2012), which are typically more challenging to segment due to a surplus of unmyelinated pathways. Although, past diffusion tractography studies were able to detect overall patterns of trajectories of TC and CC pathways in adults and pediatric populations (Jaermann et al., 2008; Lebel et al., 2010; Vasung et al., 2011; Poh et al., 2015), and our past study successfully revealed potential neuronal migration pathways of thalamic neurons and TC development in fetal brains (Wilkinson et al., 2016), there has been no report that aimed to identify convergence of TC/CT pathways and convergence of TC/CC pathways in human fetuses.
The goals of this study were to describe the spatio-temporal relationships of the development of the TC/CT and TC/CC pathways in human fetal brains, using HARDI tractography. One of our hypotheses based on the literature is that TC projections contact with pathways from the cortex. This “handshaking” has been proposed to be seen within the GE, and we specifically analyzed such convergence of TC and cortical pathways within the GE. Another hypothesis in this study is that the convergence of TC and CC pathways in the cortex can be detected with HARDI. We expected that TC and CC pathways converge in specific cortical regions, based on the knowledge that TC pathways project to broad areas in the brain but CC pathways only project to the borders of visual areas.
Materials and Methods
Postmortem Fetal Brain Specimens and in Vivo Subjects
The data of nine postmortem fetal brains were used in this study. Six postmortem brains (one 17, two 18, and three 20 GW) were obtained from the Department of Pathology, Brigham and Women's Hospital (Boston, MA, USA), under a protocol approved by the hospital's institutional review board for human research, and three postmortem brains (19, 21, and 22 GW) were provided by the Allen Institute Brain Bank (Seattle, WA, USA). Postmortem specimens were obtained with full parental consent, whereby the primary cause of death was complications of prematurity. The brains were grossly normal, and standard autopsy examination of all brains revealed minimal or no pathologic abnormalities at the macroscopic level. All research protocols were approved by the Institutional Review Board (IRB) at each institute.
All living participants (one 30 GW, one 34 GW, and three 40 GW) had clinically-indicated brain MRI studies that were diagnostically interpreted to show no abnormalities. Indications for imaging included concern for hypoxic ischemic injury, apnea, and transient choreiform movements after an upper respiratory tract infection. None had clinical concerns for a congenital malformation or genetic disorder. This in vivo part of this paper is a retrospective study of clinical MRI data, and the Institutional Review Board at Boston Children's Hospital deemed this an exempt project because the research is retrospective and involved existing data with no risk to patient confidentiality.
Tissue Preparation for HARDI
At the time of autopsy, all brains were immersion fixed and the brains from the Brigham and Women's Hospital were stored in 4% paraformaldehyde, and the brains from the Allen Institute Brain Bank were stored in 4% periodate-lysine-paraformaldehyde (PLP). During MR image acquisition, the brains from the Brigham and Women's Hospital were placed in Fomblin solution (Ausimont, Thorofare, NJ) (e.g., Takahashi et al., 2012) and the brains from the Allen Institute Brain Bank were placed in 4% PLP. These different kinds of solutions in which the brains from different institutes were placed tend to change the background contrast (i.e., we see dark background outside of the brain using Fomblin, and bright background using PLP), but they do not specifically change diffusion properties (e.g., FA and ADC) within the brain.
Diffusion MRI Procedures
The postmortem brain specimens from the Brigham and Women's Hospital were imaged using a 4.7T Bruker Biospec MR system (specimens from 17 GW to 31 GW) and specimens form the Allen Institute Brain Bank were imaged using a 3T Siemens MR system [three fetal (19, 20, 22 GW) specimens] at A. A. Martinos Center, Massachusetts General Hospital, Boston, MA. We used the 3T system as the brains from the Brain Bank were in cranio and did not fit in the 4.7T bore. To improve the imaging quality and obtain the best signal to noise and high spatial resolution, we used a custom-made MR coil with the 3T system that just fit the specimens. Different scanner systems were used to accommodate the different brain sizes, and magnetic resonance (MR) coils that best fit each brain sample were used to ensure optimal imaging.
For the brains from the Brigham and Women's hospital, a 3-dimensional (3-D) diffusion-weighted spin-echo echo-planar imaging (SE-EPI) sequence was used with a repetition time/echo time (TR/TE) of 1,000/40 ms, with an imaging matrix of 112 × 112 × 112 pixels. Sixty diffusion-weighted measurements (with the strength of the diffusion weighting, b = 8,000 s/mm2) and one non-diffusion-weighted measurement (no diffusion weighting or b = 0 s/mm2) were acquired with the duration of the diffusion gradients, δ = 12.0 ms, and the time interval between the start of the two diffusion gradients, Δ = 24.2 ms. The total acquisition time was ~2 h for each imaging session. The spatial resolution was 415 × 500 × 550 μm for the specimens at 20 and 21 gestational weeks; 440 × 520 × 720 μm for the specimen at 3 postnatal months. We determined the highest spatial resolution for each brain specimen with an acceptable signal-to-noise ratio of more than 130, and within a reasonable scan time of 2 h. For the brains from the Brain Bank, diffusion-weighted data were acquired over two averages using a steady state free precession sequence (TR/TE = 24.82/18.76 ms, α = 60°; 400 μm isotropic resolution). Diffusion weighting was isotropically distributed along 44 directions (b = 730 s/mm2) with 4 b = 0 images.
The brains of living patients were imaged on a 3T Siemens MR system, Boston Children's Hospital, Boston, MA. The diffusion pulse sequence used for imaging live participants was a diffusion-weighted spin-echo echo-planar imaging (SE-EPI) sequence, TR/TE 8320/88 ms, with an imaging matrix of 128 × 128 × 64 pixels. The spatial resolution was 2 × 2 × 2 mm. Thirty diffusion-weighted measurements (b = 1,000 s/mm2) and 5 non-diffusion-weighted measurements (b = 0 s/mm2) were acquired with δ = 40 ms and Δ = 68 ms.
Diffusion Data Reconstruction for Tractography
DiffusionToolkit and TrackVis (trackvis.org) were used to reconstruct and visualize tractography pathways. Tractography pathways were reconstructed using a HARDI model with a streamline/FACT algorithm and a 45° angle threshold. No threshold of fractional anisotropy (FA) was used for the fiber reconstruction. Brain mask volumes were used to terminate tractography structures instead of the standard FA threshold (Takahash et al., 2010; Takahashi et al., 2011, 2012, 2014; Song et al., 2015; Das and Takahashi, 2017; Vasung et al., 2017), because progressive myelination and crossing fibers in the developing brain can result in low FA values that may potentially incorrectly terminate tractography tracing in brain regions with low FA values.
Tract Delineation
Cortical, TC/CT, and CC pathways were identified in each subject by placing regions of interest (ROIs) in the cerebral cortex, thalami, and corpus callosum regions, referring to DTI and tractography atlases (Mori, 2007; Catani and Thiebaut de Schotten, 2008; Oishi et al., 2010). The size of all the ROIs were carefully optimized not to include other white matter pathways as well as not to miss the cortical, TC/CT, and CC pathways by changing the size and location several times. We also occasionally used additional ROIs to exclude clearly different pathways from the pathways of interest.
The color-coding of tractography pathways in Figure 1 was determined for a visualization purpose: thalamocortical pathways in blue and cortical pathways in yellow. The color-coding of tractography pathways in Figure 2 was based on FA values (thalamocortical in red to yellow, and corpus callosal pathways in dark to light blue).
Figure 1. Sagittal overviews of all diffusion tractography pathways starting/ending in the thalami (blue pathways) and pathways within a sagittal plane starting/ending in the cortical plate (CP) in early developmental stages or within the gray matter cortex in later stages (yellow pathways). Ages range from 17 to 40 GW. Red arrows: pathways from the thalami converged with some pathways from the cerebral cortex as early as 17 GW, mainly in anterior and middle regions of the brain until 22 GW. Green arrows: posterior pathways from the thalami also converging with cortical pathways at 30 GW and onwards. Orange arrows: a few exceptional pathways turning perpendicularly, along with majority of short-range, straight cortical pathways. Pink arrows: several groups of fiber bundles at 30 GW and onwards, making u-shaped trajectories.
Figure 2. Figure 2. Overview of thalamocortical (TC, blue) and corpus callosal (CC, red-orange-yellow) pathways from 17 to 40 GW detected by diffusion tractography, color-coded by fractional anisotropy (FA) values. While it appears that some pathways lay outside of the brain surface, these are due to the angles of the images. Light-blue arrows: pathways linked to thalami that were greater in number compared to those in posterior brain regions (17–21 GW). Blue arrows: TC pathways in the middle and posterior brain regions increased and equivalent in number with pathways linked to thalami (22 GW-). Green arrows: increased anterior TC pathways reaching the cortex by 20 GW. Red arrows: CC pathways connecting bilateral medial walls of the brain observed fully penetrated into the upper edge of the cortical plate at 17 GW. By 30 GW, global TC and CC pathways reached the cortex (pink and purple arrows, respectively). Orange arrows: condensed groups of TC pathways identified by 30 GW. Gray arrows: TC and CC bands further being condensed and many CC and TC bands overlapping with each other as the two pathways converged at 40 GW.
Results
Convergence of Emerging Thalamic and Cortical Pathways
Figure 1 shows all pathways starting/ending in the thalami (blue pathways) and pathways within a sagittal plane starting/ending in the CP in early developmental stages or within the gray matter cortex in later stages (yellow pathways; see Method section). Throughout development, TC pathways grew toward the cortex and cortical pathways descended toward the center of the brain. Some of the pathways from the thalami were found converged with some pathways from the cerebral cortex as early as 17 GW (Figure 1, red arrows). Such convergence was observed mainly in anterior and middle regions of the brain until 22 GW. At 30 GW and onwards, posterior pathways from the thalami also converged with cortical pathways (Figure 1, green arrows). A majority of the pathways linked to the thalami were terminated around the lower edge of the cortical pathways located within the cortical plate (CP) until 22 GW, with only a few pathways penetrating into the CP, while such long penetrating pathways were found more in number at and after 30 GW, and by term, these pathways were found at the upper edge of the cerebral cortex. The majority of short-range cortical pathways were straight until 22 GW with a few exceptional pathways turning perpendicularly (Figure 1, orange arrows). At 30 GW and onwards, several groups of fiber bundles were observed making u-shaped trajectories (Figure 1, pink arrows).
Convergence of Thalamocortical and Corpus Callosal Pathways
In early fetal stages, pathways linked to thalami were greater in number compared to those in posterior brain regions (17–21 GW in Figure 2, light-blue arrows). Starting from 22 GW, TC pathways in the middle and posterior brain regions increased in number, equivalent with pathways linked to thalami (Figure 2, blue arrows). No TC pathways reached the upper edge of the CP in the middle of the brain at 22 GW. By 20 GW, more anterior TC pathways reached the cortex (Figure 2, green arrows).
At 17 GW, CC pathways connecting bilateral medial walls of the brain were observed already fully penetrated into the upper edge of the CP (Figure 2, red arrows). By 30 GW, these pathways had a higher FA and were thick and dense. Interestingly, when CC pathways still had relatively lower FA values (around 0.1 in red/orange at 18 GW), the center of the TC pathways already showed higher FA values (more than 0.2 in light blue).
By 30 GW, global TC and CC pathways reached the cortex (pink and purple arrows, respectively, in Figure 2). By 30 GW, condensed groups of TC pathways were identified and they continued to thicken until term (Figure 2, orange arrows). The same trend was observed in CC pathways after 30 GW, and at the same time, we observed that CC pathways became sparser than earlier stages. At 40 GW, TC, and CC bands further condensed and many CC and TC bands overlapped with each other as the two pathways converged (Figure 2, gray arrows).
Discussion
Emerging cortical, thalamocortical (TC) and corpus callosal (CC) pathways of human postmortem fetal brains and in vivo subjects' brains ranging from 17 to 40 GW were successfully detected and studied using HARDI tractography. Descending cortical and ascending TC pathways were observed touching and running along each other as early as 17 GW. Many CC pathways reached the full length up to the cortical surface as early as 17 GW, while pathways linked to thalami (not only TC axons but also including pathways linked to thalamic neuronal migration) reached the cortical surface at and after 20 GW. These results suggest that CC pathways developed earlier than the TC pathways. Both TC and CC pathways were widespread at early stages, but projected to specific regions of the cortex and overlapped in some regions by 40 GW. Given that TC and CC pathways closely interact each other not only during development but also after maturation, it is important to better understand structural basis of their relationships during early brain development.
We previously detected pathways likely linked to neuronal migration from the ventricular zone and GE to the thalami (Wilkinson et al., 2016). Between the ventricular zone and thalami, more tractography pathways were found in anterior compared to posterior regions, which was in agreement with postnatal observations that the anterior TC segment had a higher track count and volume than the posterior segment. The neuronal migration pathways from the ventricle in the anterior thalamic regions were replaced with white matter-like pathways that did not terminate in the ventricular zone, running through deep white matter regions between 22 and 30 GW. This time frame is also when our study saw the convergence of the TC and CC pathways, suggesting that this period is crucial to the co-development of these pathways. Regarding the GE pathways, three differential pathways likely linked to neuronal migration were detected from the anterior GE to the middle part of the thalami, from the dorsal GE to the dorsal thalami, and from the mid-posterior GE to the mid-posterior thalami (Wilkinson et al., 2016). It is possible that some of these pathways were included as the detected pathways that “handshake”; pathways between the thalamus and the GE in this study. Unfortunately, it is challenging to determine using HARDI whether the detected pathways are corresponding to neuronal migration pathways or if they have already been replaced by axonal pathways.
This study showed that the onset of the TC pathway development was more delayed than that of the CC pathway, and the anterior TC pathway that extended to the frontal cortex developed first. It is also shown that the middle and anterior parts of the TC pathway that connect to other areas like the somatosensory cortex and the visual area developed later. While CC pathways became more sparse and appeared to be pruning after 22 GW, the TC pathway seemed to be thickening and began extending at around 20 GW which continued to term. Both the TC and the CC pathways were found to form condensed groups of pathways that target specific areas of the cortex, and many of these groups followed similar trajectories. While one of our past research studies (Cohen et al., 2016) showed that the CC and TC pathways took similar postnatal growth patterns, the current results suggest that in earlier developmental stages, the CC structurally matures first.
Although the TC and CC projection pathways are closely connected throughout development and are involved in many of the same systems, their reactions to critical damage are seemingly different (Miller et al., 1991). Following the formation of lesions of the thalamus in hamsters, rerouting of pathways and altering of previously existing typical pathways were found. One alteration was the addition of CC pathways coming from the temporal cortex, which are not seen in the control group. Compared to TC pathways, CC pathways have the unique ability to reorganize after damage. Miller and colleagues have hypothesized that TC pathways may be unable to alter after damage because they have developed further at this point than CC pathways have, allowing CC projections to have more plasticity further into development. Other studies have shown that disturbances with either or both of these pathways can coincide with multiple disorders, including schizophrenia (Mitelman et al., 2009; Marenco et al., 2012; Collinson et al., 2014; Klingner et al., 2014), autism (Gozzi et al., 2012; Nair et al., 2013; Hanaie et al., 2014), and bipolar disorder (Anticevic et al., 2014; Yamada et al., 2015). Callosal dysfunction can also lead to epilepsy (Andrade et al., 2014), major depressive disorder (Cyprien et al., 2014), obsessive-compulsive disorder (Jose et al., 2015) and subcortical ischemic vascular dementia (Wu et al., 2015). Not only does damage to the pathways themselves result in a number of disorders and impairments, but environmental exposures, including visual stress, can result in negative effects on the TC pathway, leading to abnormal interactions of the CC pathway (Yinon and Hammer, 1990). Yinon and Hammer provided evidence that the TC pathways are necessary for CC pathways to accurately transfer visual information between hemispheres. Another study supporting the necessary functional interactions between TC and CC pathways was performed by Pallas et al. (1999). They rerouted TC pathways in deaf-induced ferrets, which resulted in the ferrets having altered CC connections, supporting the adaptation of one pathway in response to damage of the other, in order to decrease the injurious effects.
Development of TC pathways to term suggests implications on TC pathways in the birth of premature infants. Ball et al. (2013) used tractography to establish the TC connectome of both infants born at term and premature infants, born between 23 and 34 GW. Their study provides evidence that preterm birth has a detrimental impact on TC pathways, with a decrease in TC connections to multiple brain regions, including the supplementary motor areas, occipital lobe, temporal gyri, and frontal cortices. As premature infants are born before the TC pathways have fully developed, it is no surprise that preterm birth is associated with an increased risk of multiple developmental disorders, including problems with language, attention, memory, behavior, and learning.
The current study is evidence that HARDI tractography is able to determine the development patterns of TC and CC pathways from gestational development, continuing into the postnatal time frame. Not only can researching the development patterns of the TC and CC pathways individually help further our understanding of brain development, but studying the TC and CC pathways together will help develop an understanding of the networks these pathways are involved in. Understanding human brain development at both specific and broad ranges, will lead to appropriate pre- and post-natal care for early detection and treatment of developmental disorders.
The cortical pathways in Figure 1 were identified using ROIs for the cortex and the GE. For Figure 2, CC pathways were identified using ROIs in the corpus callosum. Therefore, we believe we did not show contaminated radial pathways in this paper. Potential limitations of this study include the small number of specimens across gestational ages, the use of multiple MRI acquisition systems, and the nature of diffusion tractography. The limited number of specimens available in this analysis (one per developmental stage) would impose substantial uncertainty to any reported quantitative findings and so we have relegated this topic to future work. We believe that the issue of different acquisition systems is not significant as tractography of both systems yielded similar trends. Given the nature of diffusion tractography, even with HARDI, it is possible that the pathways observed to terminate within the GE structure may be due to strong coherent tissue structures in the GE (e.g., Miyazaki et al., 2016). However, this would not be the case, because even at 17 GW, when the GE is still strongly coherent, several pathways from the thalamus reached at the bottom of the cortical plate running through the GE.
Partial volume effects with limited spatial resolution of diffusion MRI could cause tractography pathways erroneously terminating at the brain surface. However, in our data, there are many tractography pathways that terminate before reaching the upper layer of the cortical mantle (e.g., Vasung et al., 2017). Another explanation for tractography pathways that seem to terminate at the brain surface is that the location of the anatomical image is not always in-plane with the tractography pathways shown in the figures. Since tractography pathways are 3-dimensional, selecting one structural plane to overlay the pathways on it could cause misunderstanding for the spatial relationships between tractography pathways and anatomical image.
Since our results may indicate that CC pathways develop earlier than TC pathways, and TC pathways seem to follow the CC pathways and converge with them in late gestational ages, it is possible that CC pathways are influential in the TC development. However, it is also possible that the TC pathways find their destinations independently from CC pathways. It would be necessary to study this topic in patients with callosal agenesis to answer this question, a study which is currently underway in our group.
Although it was not always clear from the data whether the brain regions where “condensed” TC and CC pathways were observed are corresponding to the locations of gyral structures, the borders of visual areas where TC and CC pathways are supposed to converge have not been linked to gyral curvatures in the literature (e.g., Caspers et al., 2015). It would be interesting to study in more detail the locations of convergent regions of TC and CC pathways in human brains during development in the future.
Author Contributions
MW and ET designed the study. ET obtained the data. RW, MW, TK, and ET processed the data. MW, TK, and ET analyzed and wrote the paper.
Conflict of Interest Statement
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Acknowledgments
This work was supported by NICHD (R01HD078561, R21HD069001) (ET) and NINDS (R03NS091587) (ET). This research was carried out in part at the Athinoula A. Martinos Center for Biomedical Imaging at the Massachusetts General Hospital, using resources provided by the Center for Functional Neuroimaging Technologies, NIH P41RR14075, a P41 Regional Resource supported by the Biomedical Technology Program of the National Center for Research Resources (NCRR), National Institutes of Health. This work also involved the use of instrumentation supported by the NCRR Shared Instrumentation Grant Program (NIH S10RR023401, S10RR019307, and S10RR023043) and High-End Instrumentation Grant Program (NIH S10RR016811). This study was conducted partly using postmortem human brain specimens from the tissue collection at the Department of Neurobiology at Yale University School of Medicine (supported by grant NIH MH081896), which form a part of the BrainSpan Consortium collection (http://www.brainspan.org). The other brain specimens were kindly provided by Dr. Rebecca D. Folkerth, Brigham and Women's Hospital. We thank Ms. Ayano Hara for data processing support, Mr. Alec Nelson and Mr. Rodney Sadberry for editorial assistance, and Dr. Lana Vasung for valuable comments for a previous version of this manuscript.
References
Afifi, A. K., and Bergman, R. A. (1997). Functional Neuroanatomy. New York, NY: McGraw-Hill Professional.
Andrade, C. S., Leite, C. C., Otaduy, M. C., Lyra, K. P., Valente, K. D., Yasuda, C. L., et al. (2014). Diffusion abnormalities of the corpus callosum in patients with malformations of cortical development and epilepsy. Epilepsy Res. 108, 1533–1542. doi: 10.1016/j.eplepsyres.2014.08.023
Anticevic, A., Cole, M. W., Repovs, G., Murray, J. D., Brumbaugh, M. S., Winkler, A. M., et al. (2014). Characterizing thalamo-coritcal disturbances in schizophrenia and bipolar illness. Cereb. Cortex 24, 3116–3130. doi: 10.1093/cercor/bht165
Bagshaw, A. P., Rollings, D. T., Khalsa, S., and Cavanna, A. E. (2014). Multimodal neuroimaging investigations of alterations to consciousness: the relationship between absence epilepsy and sleep. Epilepsy Behav. 30, 33–37. doi: 10.1016/j.yebeh.2013.09.027
Ball, G., Boardman, J. P., Aljabar, P., Pandit, A., Arichi, T., Merchant, N., et al. (2013). The influence of preterm birth on the developing thalamocortical connectome. Cortex 49, 1711–1721. doi: 10.1016/j.cortex.2012.07.006
Blakemore, C., and Molnár, Z. (1990). Factors involved in the establishment of specific interconnections between thalamus and cerebral cortex. Cold Spring Harb. Symp. Quant. Biol. 55, 491–504. doi: 10.1101/SQB.1990.055.01.048
Bocci, T., Pietrasanta, M., Cerri, C., Restani, L., Caleo, M., and Sartucci, F. (2014). Visual callosal connections: role in visual processing in health and disease. Rev. Neurosci. 25, 113–127. doi: 10.1515/revneuro-2013-0025
Burton, H., and Kopf, E. M. (1984). Connections between the thalamus and the somatosensory areas of the anterior ectosylvian gyrus in the cat. J. Comp. Neurol. 224, 173–205. doi: 10.1002/cne.902240203
Caspers, S., Axer, M., Caspers, J., Jockwitz, C., Jütten, K., Reckfort, J., et al. (2015). Target sites for transcallosal fibers in human visual cortex - A combined diffusion and polarized light imaging study. Cortex 72, 40–53. doi: 10.1016/j.cortex.2015.01.009
Catani, M., and Thiebaut de Schotten, M. (2008). A diffusion tensor imaging tractography atlas for virtual in vivo dissections. Cortex 44, 1105–1132. doi: 10.1016/j.cortex.2008.05.004
Cohen, A. H., Wang, R., Wilkinson, M., MacDonald, P., Lim, A. R., and Takahashi, E. (2016). Development of human white matter fiber pathways: from newborn to adult ages. Int. J. Dev. Neurosci. 50, 26–38. doi: 10.1016/j.ijdevneu.2016.02.002
Collinson, S. L., Gan, S. C., Woon, P. S., Kuswanto, C., Sum, M. Y., Yang, G. L., et al. (2014). Corpus callosum morphology in first-episode and chronic schizophrenia: combined magnetic resonance and diffusion tensor imaging study of Chinese Singaporean patients. Br. J. Psychiatry 204, 55–60. doi: 10.1192/bjp.bp.113.127886
Cyprien, F., Courtet, P., Poulain, V., Maller, J., Meslin, C., Bonafé, A., et al. (2014). Corpus callosum size may predict late-life depression in women: a 10-year follow-up study. J. Affect. Disord. 165, 16–23. doi: 10.1016/j.jad.2014.04.040
Das, A., and Takahashi, E. (2017). Neuronal migration and axonal pathways linked to human fetal insular development revealed by diffusion MR tractography. Cereb. Cortex. doi: 10.1093/cercor/bhx224. [Epub ahead of print].
Dehay, C., and Kennedy, H. (1988). The maturational status of thalamocortical and callosal connections of visual areas V1 and V2 in the newborn monkey. Behav. Brain Res. 29, 237–244. doi: 10.1016/0166-4328(88)90028-9
Dehay, C., Kennedy, H., and Bullier, J. (1986). Callosal connectivity of areas V1 and V2 in the newborn monkey. J. Comp. Neurol. 254, 20–33. doi: 10.1002/cne.902540103
Fabri, M., and Polonara, G. (2013). Functional topography of human corpus callosum: an FMRI mapping study. Neural Plast. 2013:251308. doi: 10.1155/2013/251308
Fabri, M., Polonara, G., Quattrini, A., Salvolini, U., Del Pesce, M., and Manzoni, T. (1999). Role of the corpus callosum in the somatosensory activation of the ipsilateral cerebral cortex: an fMRI study of callosotomized patients. Eur. J. Neurosci. 11, 3983–3994. doi: 10.1046/j.1460-9568.1999.00829.x
Frank, L. (2002). Characterization of anisotropy in high angular resolution diffusion-weighted mri. Magn. Reson. Med. 47, 1083–1099. doi: 10.1002/mrm.10156
Garel, S., and Rubenstein, J. L. (2004). Intermediate targets in formation of topographic projections: inputs from the thalamocortical system. Trends Neurosci. 27, 533–539. doi: 10.1016/j.tins.2004.06.014
Gilmore, J. H., Lin, W., Corouge, I., Vetsa, Y. S., Smith, J. K., Kang, C., et al. (2007). Early postnatal development of corpus callosum and corticospinal white matter assessed with quantitative tractography. AJNR Am. J. Neuroradiol. 28, 1789–1795. doi: 10.3174/ajnr.A0751
Gozzi, M., Nielson, D. M., Lenroot, R. K., Ostuni, J. L., Luckenbaugh, D. A., Thurm, A. E., et al. (2012). A magnetization transfer imaging study of corpus callosum myelination in young children with autism. Biol. Psychiatry 72, 215–220. doi: 10.1016/j.biopsych.2012.01.026
Häberling, I. S., Badzakova-Trajkov, G., and Corballis, M. C. (2011). Callosal tracts and patterns of hemispheric dominance: a combined fMRI and DTI study. Neuroimage 54, 779–786. doi: 10.1016/j.neuroimage.2010.09.072
Hanaie, R., Mohri, I., Kagitani-Shimono, K., Tachibana, M., Matsuzaki, J., Watanabe, Y., et al. (2014). Abnormal corpus callosum connectivity, socio-communicative deficits, and motor deficits in children with autism spectrum disorder: a diffusion tensor imaging study. J. Autism Dev. Disord. 44, 2209–2220. doi: 10.1007/s10803-014-2096-8
Hong, J. H., Son, S. M., and Jang, S. H. (2010). Identification of spinothalamic tract and its related thalamocortical fibers in human brain. Neurosci. Lett. 468, 102–105. doi: 10.1016/j.neulet.2009.10.075
Jaermann, T., De Zanche, N., Staempfli, P., Pruessmann, K. P., Valavanis, A., Boesiger, P., et al. (2008). Preliminary experience with visualization of intracortical fibers by focused high-resolution diffusion tensor imaging. AJNR Am. J. Neuroradiol. 29, 146–150. doi: 10.3174/ajnr.A0742
Johansen-Berg, H., Behrens, T. E., Sillery, E., Ciccarelli, O., Thompson, A. J., Smith, S. M., et al. (2005). Functional-anatomical validation and individual variation of diffusion tractography-based segmentation of the human thalamus. Cereb. Cortex 15, 31–39. doi: 10.1093/cercor/bhh105
Jose, D., Narayanaswamy, J. C., Agarwal, S. M., Kalmady, S. V., Venkatasubramanian, G., and Reddy, Y. C. (2015). Corpus callosum abnormalities in medication-naive adult patients with obsessive compulsive disorder. Psychiatry Res. 231, 341–345. doi: 10.1016/j.pscychresns.2015.01.019
Kennedy, H., Dehay, C., and Bullier, J. (1986). Organization of the callosal connections of visual areas V1 and V2 in the macaque monkey. J. Comp. Neurol. 247, 398–415. doi: 10.1002/cne.902470309
Killackey, H. P., and Belford, G. R. (1979). The formation of afferent patterns in the somatosensory cortex of the neonatal rat. J. Comp. Neurol. 183, 285–303. doi: 10.1002/cne.901830206
Klein, J. C., Rushworth, M. F., Beherens, T. E., Mackay, C. E., Crespigny, A. J., D'Arceuil, H., et al. (2010). Topography of connections between human prefrontal cortex and mediodorsal thalamus studied with diffusion tractography. Neuroimage 51, 555–564. doi: 10.1016/j.neuroimage.2010.02.062
Klingner, C. M., Langbein, K., Dietzek, M., Smesny, S., Witte, O. W., Sauer, H., et al. (2014). Thalamocortical connectivity during resting state in schizophrenia. Eur. Arch. Psychiatry Clin. Neurosci. 264, 111–119. doi: 10.1007/s00406-013-0417-0
Lebel, C., Caverhill-Godkewitsch, S., and Baeulieu, C. (2010). Age-related regional variations of the corpus callosum identified by diffusion tensor tractography. Neuroimage 52, 20–31. doi: 10.1016/j.neuroimage.2010.03.072
Leyva-Díaz, E., and López-Bendito, G. (2013). In and out from the cortex: development of major forebrain connections. Neuroscience 254, 26–44. doi: 10.1016/j.neuroscience.2013.08.070
López-Bendito, G., and Molnar, Z. (2003). Thalamocortical development: how are we going to get there? Nat. Rev. Neurosci. 4, 276–289. doi: 10.1038/nrn1075
Luders, E., Phillips, O. R., Clark, K., Kurth, F., Toga, A. W., and Narr, K. L. (2012). Bridging the hemispheres in meditation: thicker callosal regions and enhanced fractional anisotropy (FA) in long-term practitioners. Neuroimage 61, 181–187. doi: 10.1016/j.neuroimage.2012.02.026
Magnuson, M. E., Thompson, G. J., Pan, W. J., and Keilholz, S. D. (2014). Effects of severing the corpus callosum on electrical and BOLD functional connectivity and spontaneous dynamic activity in the rat brain. Brain Connect. 4, 15–29. doi: 10.1089/brain.2013.0167
Marcano-Reik, A. J., and Blumberg, M. S. (2008). The corpus callosum modulates spindle-burst activity within homotopic regions of somatosensory cortex in newborn rats. Eur. J. Neurosci. 28, 1457–1466. doi: 10.1111/j.1460-9568.2008.06461.x
Marenco, S., Stein, J. L., Savostyanova, A. A., Sambataro, F., Tan, H. Y., Goldman, A. L., et al. (2012). Investigation of anatomical thalamo-cortical connectivity and FMRI activation in schizophrenia. Neuropsychopharmacology 37, 499–507. doi: 10.1038/npp.2011.215
Métin, C., and Godement, P. (1996). The ganglionic eminence may be an intermediate target for corticofugal and thalamocortical axons. J. Neurosci. 16, 3219–3235.
Miller, B., Windrem, M. S., and Finlay, B. L. (1991). Thalamic ablations and neocortical development: alterations in thalamic and callosal connectivity. Cereb. Cortex 1, 241–261. doi: 10.1093/cercor/1.3.241
Mitelman, S. A., Nikiforova, Y. K., Canfield, E. L., Hazlett, E. A., Brickman, A. M., Shihabuddin, L., et al. (2009). A longitudinal study of the corpus callosum in chronic schizophrenia. Schizophr. Res. 114, 144–153. doi: 10.1016/j.schres.2009.07.021
Miyazaki, Y., Song, J. W., and Takahashi, E. (2016). Asymmetry of radial and symmetry of tangential neuronal migration pathways in developing human fetal brains. Front. Neuroanat. 10:2. doi: 10.3389/fnana.2016.00002
Molnár, Z., Adams, R., and Blakemore, C. (1998). Mechanisms underlying the early establishment of thalamocortical connections in the rat. J. Neurosci. 18, 5723–5745.
Nair, A., Treiber, J. M., Shukla, D. K., Shih, P., and Müller, R. A. (2013). Impaired thalamocortical connectivity in autism spectrum disorder: a study of functional and anatomical connectivity. Brain 136, 1942–1955. doi: 10.1093/brain/awt079
Oishi, K., Faria, A. V., van Zijl, P. C. M., and Mori, S. (2010). MRI Atlas of Human White Matter, 2nd Edn. Amsterdam: Elsevier.
Pallas, S. L., Littman, T., and Moore, D. R. (1999). Cross-modal reorganization of callosal connectivity without altering thalamocortical projections. Proc. Natl. Acad. Sci. U.S.A. 96, 8751–8756. doi: 10.1073/pnas.96.15.8751
Poh, J. S., Li, Y., Ratnarajah, N., Fortier, M. V., Chong, Y. S., Kwek, K., et al. (2015). Developmental synchrony of thalamocortical circuits in the neonatal brain. Neuroimage 1, 168–176. doi: 10.1016/j.neuroimage.2015.03.039
Rakic, P. (1975). “Timing of major ontogenetic events in the visual cortex of the monkey,” in Brain Mechanisms in Mental Retardation, eds N. A. Buckwald, and M. Brazier (New York, NY: Academic Press), 3–40.
Rakic, P. (1988). Defects of neuronal migration and the pathogenesis of cortical malformations. Prog. Brain Res. 73, 15–37. doi: 10.1016/S0079-6123(08)60494-X
Robertson, R. T. (1976). Thalamic projections to visually responsive regions of parietal cortex. Brain Res. Bull. 1, 459–469. doi: 10.1016/0361-9230(76)90116-7
Sarnthein, J., Morel, A., von Stein, A., and Jeanmonod, D. (2005). Thalamocortical theta coherence in neurological patients at rest and during a working memory task. Int. J. Psychophysiol. 57, 87–96. doi: 10.1016/j.ijpsycho.2005.03.015
Seeley, W. W., Menon, V., Schatzberg, A. F., Keller, J., Glover, G. H., Kenna, H., et al. (2007). Dissociable intrinsic connectivity networks for salience processing and executive control. J. Neurosci. 27, 2349–2356. doi: 10.1523/JNEUROSCI.5587-06.2007
Song, J. W., Mitchell, P. D., Kolasinski, J., Ellen Grant, P., Galaburda, A. M., and Takahashi, E. (2015). Asymmetry of white matter pathways in developing human brains. Cereb. Cortex 25, 2883–2893. doi: 10.1093/cercor/bhu084
Takahashi, E., Dai, G., Rosen, G. D., Wang, R., Ohki, K., Folkerth, R. D., et al. (2011). Developing neocortex organization and connectivity in cats revealed by direct correlation of diffusion tractography and histology. Cereb. Cortex 21, 200–211. doi: 10.1093/cercor/bhq084
Takahashi, E., Dai, G., Wang, R., Ohki, K., Rosen, G. D., Galaburda, A. M., et al. (2010). Development of cerebral fiber pathways in cats revealed by diffusion spectrum imaging. Neuroimage 49, 1231–1240. doi: 10.1016/j.neuroimage.2009.09.002
Takahashi, E., Folkerth, R. D., Galaburda, A. M., and Grant, P. E. (2012). Emerging cerebral connectivity in the human fetal brain: an MR tractography study. Cereb. Cortex 22, 455–464. doi: 10.1093/cercor/bhr126
Takahashi, E., Hayashi, E., Schmahmann, J. D., and Grant, P. E. (2014). Development of cerebellar connectivity in human fetal brains revealed by high angular resolution diffusion tractography. Neuroimage 96, 326–333. doi: 10.1016/j.neuroimage.2014.03.022
Tournier, J. D., Calamante, F., and Connelly, A. (2007). Robust determination of the fibre orientation distribution in diffusion MRI: non-negativity constrained super-resolved spherical deconvolution. Neuroimage 35, 1459–1472. doi: 10.1016/j.neuroimage.2007.02.016
Tournier, J. D., Calamante, F., Gadian, D. G., and Connelly, A. (2004). Direct estimation of the fiber orientation density function from diffusion-weighted MRI data using spherical deconvolution. Neuroimage 23, 1176–1185. doi: 10.1016/j.neuroimage.2004.07.037
Tuch, D. S., Reese, T. G., Wiegell, M. R., and Wedeen, V. J. (2003). Diffusion MRI of complex neural architecture. Neuron 40, 885–895. doi: 10.1016/S0896-6273(03)00758-X
Van Essen, D. C. (2005). Corticocortical and thalamocortical information flow in the primate visual system. Prog. Brain Res. 149, 173–185. doi: 10.1016/S0079-6123(05)49013-5
Vasung, L., Jovanov-Milošević, N., Pletikos, M., Mori, S., Judaš, M., and Kostović, I. (2011). Prominent periventricular fiber system related to ganglionic eminence and striatum in the human fetal cerebellum. Brain Struct. Funct. 215, 237–253. doi: 10.1007/s00429-010-0279-4
Vasung, L., Raguz, M., Kostovic, I., and Takahashi, E. (2017). Spatiotemporal relationship of brain pathways during human fetal development using high-angular resolution diffusion MR imaging and histology. Front. Neurosci. 11:348. doi: 10.3389/fnins.2017.00348
Viaene, A. N., Petrof, I., and Sherman, S. M. (2011). Synaptic properties of thalamic input to the subgranular layers of primary somatosensory and auditory cortices in the mouse. J. Neurosci. 31, 12738–12747. doi: 10.1523/JNEUROSCI.1565-11.2011
Wang, Q., Gao, E., and Burkhalter, A. (2007). In vivo transcranial imaging of connections in mouse visual cortex. J. Neurosci. Methods 159, 268–276. doi: 10.1016/j.jneumeth.2006.07.024
Whitehead, M. T., Raju, A., and Choudhri, A. F. (2014). Normal centrolineal myelination of the callosal splenium reflects the development of the cortical origin and size of its commissural fibers. Neuroradiology 56, 333–338. doi: 10.1007/s00234-014-1323-5
Wilkinson, M., Kane, T., Wang, R., and Takahashi, E. (2016). Migration pathways of thalamic neurons and development of thalamocortical connections in humans revealed by diffusion MR tractography. Cereb. Cortex. doi: 10.1093/cercor/bhw339. [Epub ahead of print].
Wu, X. P., Gao, Y. J., Yang, J. L., Xu, M., and Sun, D. H. (2015). Quantitative measurement to evaluate morphological changes of the corpus callosum in patients with subcortical ischemic vascular dementia. Acta Radiol. 56, 214–218. doi: 10.1177/0284185114520863
Yamada, K., Nagakane, Y., Yoshikawa, K., Kizu, O., Ito, H., Kubota, T., et al. (2007). Somatotopic organization of thalamocortical projection fibers as assessed with MR tractography. Radiology 242, 840–845. doi: 10.1148/radiol.2423060297
Yamada, S., Takahashi, S., Ukai, S., Tsuji, T., Iwatani, J., Tsuda, K., et al. (2015). Microstructural abnormalities in anterior callosal fibers and their relationship with cognitive function in major depressive disorder and bipolar disorder: a tract-specific analysis study. J. Affect. Disord. 174, 542–548. doi: 10.1016/j.jad.2014.12.022
Yinon, U., and Hammer, A. (1990). Post critical period plasticity of callosal transfer to visual cortex cells of cats following early conditioning of monocular deprevation and late optic chiasm transection. Brain Res. 516, 84–90. doi: 10.1016/0006-8993(90)90900-V
Keywords: thalamocortical, corticothalamic, corpus callosal, development, fetus, human, tractography
Citation: Wang R, Wilkinson M, Kane T and Takahashi E (2017) Convergence of Cortical, Thalamocortical, and Callosal Pathways during Human Fetal Development Revealed by Diffusion MRI Tractography. Front. Neurosci. 11:576. doi: 10.3389/fnins.2017.00576
Received: 08 February 2017; Accepted: 02 October 2017;
Published: 01 November 2017.
Edited by:
Ahmet O. Caglayan, Istanbul Bilim University, TurkeyReviewed by:
Zeljka Krsnik, School of Medicine, University of Zagreb, CroatiaAndras Jakab, University Children's Hospital Zurich, Switzerland
Copyright © 2017 Wang, Wilkinson, Kane and Takahashi. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) or licensor are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Emi Takahashi, emi@nmr.mgh.harvard.edu
†Present Address: Rongpin Wang, Department of Radiology, Guizhou Provincial People's Hospital, Guiyang, China
‡These authors have contributed equally to this work.