Asymmetry of Radial and Symmetry of Tangential Neuronal Migration Pathways in Developing Human Fetal Brains
Yuta Miyazaki
Jae W. Song
Emi Takahashi- 1Department of Medicine, Chiba University School of Medicine, Chiba, Japan
- 2 Department of Radiology and Biomedical Imaging, Yale University School of Medicine, New Haven, CT, USA
- 3Division of Newborn Medicine, Department of Medicine, Boston Children’s Hospital, Harvard Medical School, Boston, MA, USA
- 4Athinoula A. Martinos Center for Biomedical Imaging, Massachusetts General Hospital, Harvard Medical School, Charlestown, MA, USA
- 5Fetal-Neonatal Neuroimaging and Developmental Science Center, Boston Children’s Hospital, Harvard Medical School, Boston, MA, USA
The radial and tangential neural migration pathways are two major neuronal migration streams in humans that are critical during corticogenesis. Corticogenesis is a complex process of neuronal proliferation that is followed by neuronal migration and the formation of axonal connections. Existing histological assessments of these two neuronal migration pathways have limitations inherent to microscopic studies and are confined to small anatomic regions of interest (ROIs). Thus, little evidence is available about their three-dimensional (3-D) fiber pathways and development throughout the entire brain. In this study, we imaged and analyzed radial and tangential migration pathways in the whole human brain using high-angular resolution diffusion MR imaging (HARDI) tractography. We imaged ten fixed, postmortem fetal (17 gestational weeks (GW), 18 GW, 19 GW, three 20 GW, three 21 GW and 22 GW) and eight in vivo newborn (two 30 GW, 34 GW, 35 GW and four 40 GW) brains with no neurological/pathological conditions. We statistically compared the volume of the left and right radial and tangential migration pathways, and the volume of the radial migration pathways of the anterior and posterior regions of the brain. In specimens 22 GW or younger, the volume of radial migration pathways of the left hemisphere was significantly larger than that of the right hemisphere. The volume of posterior radial migration pathways was also larger when compared to the anterior pathways in specimens 22 GW or younger. In contrast, no significant differences were observed in the radial migration pathways of brains older than 22 GW. Moreover, our study did not identify any significant differences in volumetric laterality in the tangential migration pathways. These results suggest that these two neuronal migration pathways develop and regress differently, and radial neuronal migration varies regionally based on hemispheric and anterior-posterior laterality, potentially explaining regional differences in the amount of excitatory neurons that migrate along the radial scaffold.
Introduction
Corticogenesis is a complex process that occurs during a critical period of brain development. This process begins with subcortical neuroproliferation, followed by neuronal migration and the formation and refinement of axonal connections (Caviness et al., 1996). Neuronal subtypes of different origins migrate differentially along transient pathways to ultimately establish connections that result in a highly conserved, six-layered neocortical structure. In developing human brains, the radial and tangential migration pathways are two primary sources of neural progenitor cells.
Radial glia are recognized as a primary source of progenitor cells and also serve as a scaffold for migrating neurons in the developing brain. During corticogenesis, a complex process of both symmetric and asymmetric divisions of radial glial cells from the ventricular zone (VZ) has been described, with symmetric divisions attributed to increasing the progenitor pool and asymmetric divisions resulting in neuronal generation (Cai et al., 2002). Neuron subtypes that originate from the radial migration pathway include excitatory glutamatergic projection neurons and GABAergic neurons (Haubensak et al., 2004), which migrate along radial glial fascicles toward the cortical plate (CP; Rakic, 1972, 1988; Bystron et al., 2008; Petanjek et al., 2009a,b). One histologic study showed vimentin immunoreactive radial glia in the developing visual cortex in the occipital lobe of the human brain at specific pre- (14 gestational weeks (GW), 19 GW, 36 GW) and post-term (79 GW) ages, describing the fate and destination of radial glia towards this region during development (Honig et al., 1996). Other studies investigated radial pathways with immunohistochemistry at 10–19 GW (Bayatti et al., 2008) and with lipophilic tracing methods at 14–15 GW (Hansen et al., 2010). However, detailed information regarding radial glial pathways remains unknown as these studies lacked specimens from a critical period of brain development (19–36 GW).
Another major neuronal migratory pathway is the tangential migration pathway, which runs tangential to the CP and from which inhibitory GABAergic interneuron progenitors arise (Wonders and Anderson, 2006). The tangential migration pathway originates from the ganglionic eminence (GE), and although the VZ/subventricular zone (SVZ) are also suggested to be important sources of GABAergic interneurons (Radonjić et al., 2014; Al-Jaberi et al., 2015; Arshad et al., 2015), the GE has been suggested to be a major pool of GABAergic interneurons (Hansen et al., 2013; Ma et al., 2013; Keverne, 2014). Both pathways have been extensively studied in rodents (e.g., Price and Thurlow, 1988; Walsh and Cepko, 1988; Anderson et al., 1997), primates (Rakic, 1972; Radonjić et al., 2014) and to a lesser degree, humans (e.g., Rakić and Sidman, 1969; Letinic and Rakic, 2001; Al-Jaberi et al., 2015). Taken together with the literature on rodents and primates, recent studies on human fetal brains showed that certain aspects of neuronal migration are unique to human corticogenesis, emphasizing the need to examine trajectories of neuronal migration pathways directly in the developing human brain (Letinic and Rakic, 2001; Al-Jaberi et al., 2015). However, due to limitations inherent to microscopic studies, only small regions of neuronal migratory streams have been histologically investigated, and a systematic understanding of the neuronal migratory streams during the human fetal stage is still elusive.
Initial efforts to assess migratory streams include whole brain analysis by magnetic resonance imaging (MRI; Zhang et al., 2013), and analyses of hemispheric asymmetries (Snyder et al., 1995; Kovalev et al., 2003). Several groups have imaged diffusion coherence of migration pathways by conventional diffusion tensor imaging (DTI) tractography (McKinstry et al., 2002; Huang et al., 2009), but many questions still remain regarding the time-course of the development and regression of the migratory streams as well as the need for comprehensive analyses utilizing whole brain specimens and wider age ranges. Moreover, conventional DTI methods are challenged by the reconstruction of crossing neural fibers in the brain (Tuch et al., 2003).
High-angular resolution diffusion imaging (HARDI) tractography is a technique that enables reconstruction of complicated crossing neural fibers in the brain (Tuch et al., 2003), even in immature brains (Takahashi et al., 2011, 2012, 2013, 2014; Song et al., 2014). Given the technical advantage of HARDI tractography, the present study aims to analyze radial and tangential migration pathways in developing human brains using this technique to understand patterns of neuronal migration during a critical period of brain development.
Disorganized neuronal migration from the radial and tangential pathways has been implicated in many neurological and psychiatric disorders (Marín and Rubenstein, 2001; Gressens, 2006; Volpe, 2009) such as intractable epilepsy (Cusmai et al., 2014), developmental dyslexia (Adler et al., 2013; Platt et al., 2013), autism (Wegiel et al., 2010), and schizophrenia (Benes and Berretta, 2001; McIntosh et al., 2008; Hori et al., 2014; Muraki and Tanigaki, 2015). Thus, a better understanding of this complex developmental process is critical to gain further insight into such disease processes.
Materials and Methods
Fetal Brain Specimens
Ten fetal brain specimens were imaged for this study. IRB committees at Partners and Boston Children’s Hospital (BCH) approved the use of postmortem specimens for MRI studies. The brains were acquired after obtaining informed consent by the parents. Table 1 reports the details of GW, the sources of specimens, and whether the brain was imaged ex- or in vivo. Seven postmortem brain specimens were provided by the Department of Pathology, Brigham and Women’s Hospital (BWH; Boston, MA, USA), and three postmortem brain specimens were obtained from the Allen Institute Brain Bank (AIBB; Seattle, WA, USA). These brains were grossly normal, and standard autopsy examination of all brains undergoing postmortem HARDI revealed minimal or no pathologic abnormalities at the macroscopic level. Any cases with suspected abnormality, malformations and disruption were excluded from this study. Ex vivo fetal and postnatal brains and living postnatal infant brains were included in the study, which allowed us to reveal brain structure in unforeseen detail. Although ex vivo imaging is superior, we have previously shown that ex- and in vivo imaging produces comparable tracking results in brains at the fetal, neonatal, and toddler ages (Xu et al., 2014).
Table 1. Characteristics of brain specimens.
In Vivo Subjects
Eight living participants underwent clinically-indicated brain MRI studies that were interpreted to show no abnormalities. The IRB at BCH approved the retrospective use of in vivo imaging data. 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. All MR image acquisitions were performed under protocols approved by each hospital’s institutional review board for human research.
Tissue Preparation for HARDI
At the time of autopsy, all brains were immersion fixed. The brains from BWH were stored in 4% paraformaldehyde, and the brains from AIBB were stored in 4% periodate-lysine-paraformaldehyde (PLP). During MR image acquisition, BWH brains were placed in Fomblin solution (Ausimont, Thorofare, NJ, USA; Takahashi et al., 2012), and AIBB brains were placed in 4% PLP. While these different solutions tend to change the background contrast (i.e., a dark background outside of the brain is often visualized using Fomblin, and a bright background using PLP), these solutions do not specifically change diffusion properties (e.g., FA and ADC) within the brain parenchyma (Kolasinski et al., 2013; Song et al., 2014).
Diffusion MRI Procedures
Six postmortem brain specimens from BWH were imaged with a 4.7T Bruker Biospec MR system, and three specimens from the AIBB were imaged with a 3T Siemens MR system at the A. A. Martinos Center, Massachusetts General Hospital, Boston, MA, USA. The 3T system was used to accommodate the AIBB brains that 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 custom-made radio-frequency (RF) coils with one channel on the 4.7T and 3T systems (Takahashi et al., 2012; Xu et al., 2014). We used multiple scanner systems to utilize the best-fit RF coils to ensure optimal imaging with good signal to noise.
For the BWH brains, a three-dimensional (3-D) diffusion-weighted spin-echo echo-planar imaging (SE-EPI) sequence was used with a repetition time/echo time (TR/TE) of 1000/40 ms, with an imaging matrix of 112 × 112 × 112 pixels. Sixty diffusion-weighted measurements (with the strength of the diffusion weighting, b = 8000 s/mm2) and one non-diffusion-weighted measurement (no diffusion weighting or b = 0 s/mm2) were acquired with δ = 12.0 ms and Δ = 24.2 ms. The spatial resolution was 440 × 500 × 500 μm. For the brains from the AIBB, diffusion-weighted data were acquired over two averages using a steady-state free-precession sequence with TR/TE = 24.82/18.76 ms, α = 60°, and the spatial resolution was 400 × 400 × 400 μm. Diffusion weighting was isotropically distributed along 44 directions (b = 730 s/mm2) with 4 b = 0 images. We determined the highest spatial resolution for each brain specimen with an acceptable signal-to-noise ratio of more than 130 for tractography.
The brains of living participants were imaged on a 3T Siemens MR system, BCH, Boston, MA, USA. The diffusion pulse sequence used was a diffusion-weighted 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 = 1000 s/mm2) and five non-diffusion-weighted measurements (b = 0 s/mm2) were acquired with δ = 40 ms and Δ = 68 ms.
Reconstruction and Identification of Tractography Pathways
We used Diffusion Toolkit and TrackVis1 to reconstruct and visualize neural migration pathways. A streamline algorithm for diffusion tractography was used (Mori et al., 1999), described in previous publications (Schmahmann et al., 2007; D’Arceuil et al., 2008; Takahashi et al., 2010, 2011, 2012). The term “streamline” refers to connecting tractography pathways using a local maximum or maxima. This method is true for both DTI and HARDI. The streamline technique is limited in its ability to resolve crossing pathways when used with the traditional DTI technique because one simply connects the direction of the principal eigenvector on a tensor to produce the DTI tractography pathways. This feature is a recognized limitation of DTI (Mori et al., 1999). Hence, in the current study, we used HARDI, which can theoretically detect multiple local maxima on an orientation distribution function (ODF). Using each local maxima on an ODF, we applied the streamline algorithm to initiate and continue tractography (Tuch et al., 2003), thus enabling us to identify crossing pathways within a voxel.
Trajectories were propagated by consistently pursuing the orientation vector of least curvature. Tracking was terminated when the angle between two consecutive orientation vectors was greater than the given threshold (40°) or when the fibers extended outside of the brain surface using a brain mask. The brain mask volumes were used to terminate tractography structures instead of the FA threshold (Schmahmann et al., 2007; Wedeen et al., 2008; Takahashi et al., 2010, 2012; Vishwas et al., 2010), because progressive myelination and crossing fibers in the developing brain can result in low FA values and the use of an FA threshold may potentially incorrectly terminate tractography tracking in regions with low FA values.
Regions of Interest (ROIs)
Using previously reported methodology (Takahashi et al., 2012; Kolasinski et al., 2013), a regions of interest (ROIs) approach was used to identify radial and tangential tractogrpahy pathways. ROIs for the CP, VZ/SVZ, and GE were identified (Figures 1A,B). We used ROIs of the region of VZ/SVZ and CP for identification of the radial migration pathway (Figure 1A). Tangential pathways were identified using ROIs of the GE (Figure 1B). Previous studies also have shown the GE structure by manually identifying its location (Huang et al., 2009; Kostovic and Vasung, 2009; Takahashi et al., 2012; Song et al., 2014).
Figure 1. Regions of Interest (ROIs) used to identify radial and tangential tractograpy pathways. (A) ROIs placed in cortical plate (CP, blue) and ventricular/subventricular zone (VZ/SVZ; left: red, right: orange), (B) ROIs placed in a middle part of ganglionic eminence (GE; left: green, right: pink).
Statistical Analyses
Volume measurements are automatically derived by TrackVis for each migration pathway. The volume of pathways was calculated by counting the number of voxels that touched or passed through the detected pathways, and was compared at different developmental time points.
We defined left and right radial pathways separately between the VZ/SVZ and CP in each hemisphere. We also defined the anterior and posterior radial pathways as the front or posterior half of the radial pathway. Due to challenges in defining the “middle” coronal plane in early fetal stages from only anatomic information without gyral or other mature subcortical structures, we instead measured the length between the anterior and posterior edges of the brain and defined a middle coronal plane to separate the anterior and posterior radial pathways. Left/right hemispheric asymmetry was tested for both the radial and tangential migration pathways, and anterior-posterior differences were tested for the radial pathway.
We statistically analyzed asymmetry with a laterality index (LI), calculated by the following formula:
LIs range from −2 to 2, with positive and negative LI values corresponding to left- and rightward, and anterior and posterior asymmetries, respectively (Caviness et al., 1996). The asymmetries of the radial and tangential migration pathways using LIs were statistically analyzed with SPSS version 19.0.0 (IBM SPSS, Chicago, IL, USA). Statistical significance was set at P < 0.01.
Results
Qualitative Description of the Development and Regression Time-Lines of Radial and Tangential Migration Pathways
Development and regression of the radial and tangential migration pathways are shown in Figures 2, 3. Radial and tangential migration pathways were both identified in the early age group (17–22 GW), while only radial pathways were identified in the late age group (30–40 GW). The tangential migration pathways in the late age group (30 and 34 GW) became sparse, and after 34 GW they were not identified by tractography.
Figure 2. Radial pathways. Tractography at representative ages of the radial pathways are shown in sagittal views. White arrowheads show brain regions where the radial pathways became progressively sparse. Scale bar = 5 cm for 17–21 GW, and 10 cm for 22–40 GW.
Figure 3. Tangential migration pathways. Tractography at representative ages of the tangential migration pathways are shown in sagittal views. White arrowheads show examples of the tangential migration tract branching into the cerebral walls. Scale bar = 5 cm for 17–21 GW, and 10 cm for 22–40 GW.
Radial Migration Pathways
Figure 2 shows reconstructed radial tractography pathways (perpendicular to the cortical mantle) running through the cerebral wall in the sagittal views (Figure 2). Radial pathways gradually decreased density by 21 GW (Figure 2, white arrowheads). Residual fiber pathways were still observed at 40 GW in a few brains suggestive of individual variability. Gradual disappearance of the radial pathways was not evenly observed in the whole brain. The disappearance of radial pathways by tractography started in the dorsal parietal and superior frontal regions at 21 GW and progressed in an anterior, inferior direction towards the inferior frontal lobe after 21 GW (Figure 2, white arrowheads). At 30 GW, a large segment of the radial pathways persisted in the occipital lobe. However, by 40 GW, the density of the radial pathways appeared substantially decreased, globally.
Tangential Migration Pathways
Reconstructed GE tractography pathways are shown in Figure 3. Pathways likely corresponding to the GE migration stream were very thick at 17 GW and revealed some radial pathways emanating from the stream at the most distal anterior and posterior regions as it penetrated the cortical mantle (Figure 3, white arrowheads). These radial pathways were contiguous with tangentially oriented pathways, as previously described by Kolasinski et al. (2013). Until 22 GW, tangential migration pathways appeared tightly packed with gradual thinning in the sagittal planes (Figure 3). After 30 GW, the tangential migration pathways became sparser, and after 34 GW, the tangential migration pathways were not identifiable by tractography (Figure 3, white arrowheads).
Quantitative Results
Hemispheric Asymmetry
Next, we studied the volumetric laterality indices (LIs) for the radial and tangential migration pathways in both the early and late ages (Figures 4A,B). Significant leftward asymmetry in the radial pathways in the early age group was observed (p = 0.0045). In contrast, asymmetry in the late age group was not observed (p > 0.01). When combining all ages, significant leftward asymmetry of the radial pathways emerged (p = 0.0019). By contrast, statistically significant hemispheric asymmetry was not observed in any age groups for the GE pathway (p > 0.01).
Figure 4. Laterality indices (LIs) for hemispheric asymmetry (A,B) and anterior-posterior asymmetry (C) of radial pathways (A,C) and tangential pathways (B). Radial pathways in the early age group revealed significant hemispheric asymmetry (p = 0.0045) (A). In contrast, significant hemispheric asymmetry was not observed in the tangential pathways (p > 0.01) (B). In the early age group, the posterior region was significantly larger than the anterior region (p < 0.0001, [1.75 × 10−5]) (C). In contrast, no significant difference was observed in the late age group (C). Diamond shape: early age group (17–22 GW); square: late age group (30–40 GW).
Anterior-Posterior Differences
Results of comparisons between the volume of the anterior and posterior regions of the identified radial pathways are shown in Figure 4C. In the early age group, the posterior region was significantly larger than the anterior region (p < 0.01). In contrast, no significant difference was observed in the late age group of the radial pathway (p > 0.01). Combining all ages, no significant volumetric difference was observed between the anterior and posterior regions of the radial pathways (p = 0.55).
Discussion
We investigated two major neural migration pathways, the radial and tangential migration pathways, in human fetal and newborn brains (17–40 GW) using HARDI tractography. Due to limitations of microscopy, previous studies have performed only limited regional analyses of these migration pathways in developing human brains. Our new findings are: (1) significant left-right hemispheric asymmetry and anterior-posterior asymmetry of radial migration pathways between 17 GW and 22 GW by analyzing whole human fetal brain with MR imaging; (2) there were no significant asymmetry of tangential migration pathways between 17 GW and 22 GW, and no significant asymmetry of radial migration pathways between 30 GW and 40 GW. Our previous study (Kolasinski et al., 2013) reported on the radial and tangential migration pathways, but the study just used three specimens between 19 GW and 21 GW and qualitatively described trajectories of pathways without investigating asymmetry. These results suggest that these two neuronal migration pathways develop and regress differently. This study showed, for the first time, hemispheric and anterior-posterior asymmetry of neuronal migration pathways in human fetal brains.
Radial Migration Pathway—Development and Regression Time Line
Many lines of evidence suggest that the majority of radial pathways contribute excitatory glutamatergic projection neurons (Rakic, 1972, 1988; Bystron et al., 2008; Petanjek et al., 2009a,b). Therefore, the volume of the radial pathway may, in part, reflect the relative degree of density of this subtype of glutamatergic neurons originating from this pathway. Our study revealed that the disappearance of radial pathways began at 21 GW. This thinning of radial tractography pathways was most conspicuous in the dorsal parietal and frontal regions and progressed to the inferior frontal lobe around 22 GW. This observation suggests a determined spatiotemporal direction of migration of excitatory neurons that arise from the radial pathway.
Radial Migration Pathway—Hemispheric Asymmetry
While the radial tractography pathways showed structural asymmetries, the tangential migration pathways revealed no significant hemispheric asymmetry in volume in this study. These results suggest that, at least in the pathways that supply inhibitory interneurons to the cortex, there is no demonstrable laterality by tractography.
The volume of the radial pathway in the left hemisphere was significantly larger than that of the right hemisphere in the early age group (17–22 GW) in our study. In contrast, no significant difference between the left- and right hemispheres was observed in the late age group (30–40 GW). Previous studies reported that the left hemisphere is volumetrically larger than the right hemisphere in right-handed adults (Büchel et al., 2004), and it has been suggested that the leftward asymmetry in Broca’s area contributes to human language acquisition (Schenker et al., 2010). Such normal hemispheric asymmetry in adults is often not observed in neurological and psychiatric disorders (e.g., Mackay et al., 2003). However, to date, very few studies have explored the developmental origins of this asymmetry. Although handedness was not confirmed, a larger left temporal area was already observed in the fetal stage by 30 GW (Chi et al., 1977). Even earlier, leftward asymmetry of gyral folding in specific temporal areas started to be observed from 23 GW (Habas et al., 2012). By term (or term-equivalent ages), leftward asymmetries were also detected in some tractography pathways likely corresponding to the superior longitudinal fasciculus, thalamocortical pathways, and the corticospinal tract (Liu et al., 2010). Given the link between neuronal migration and gyral formation in cobblestone lissencephaly (Dobyns et al., 1996), classical lissencephaly (Glesson, 2000; Lammens, 2000; Leventer et al., 2000), and polymicrogyria (Chang et al., 2004; Guerrini and Filippi, 2005; Jansen and Andermann, 2005), one can speculate that the leftward volumetric asymmetry of radial neuronal migration pathways may contribute to hemispheric asymmetry of the number of neurons in the cortex, followed by hemispheric asymmetry of axonal pathways from the neurons. In fact, the earliest gestational age where Habas et al. (2012) found hemispheric asymmetry of gyral structures (23 GW) coincides with the early age group of the current study (17–22 GW) when we observed leftward asymmetry. Future studies are necessary to directly link the transient leftward asymmetry of the radial pathway to the hemispheric asymmetry in mature brains.
Radial Migration Pathway—Anterior-Posterior Asymmetry
In the current study, the volume of the radial pathway in the posterior region was significantly larger compared to the anterior region when examining ages up to 22 GW in our study. This quantitative analysis was in agreement with the visual tractography of the radial pathways, which revealed a gradual thinning of this pathway in the parieto-occipital regions at 21 GW. This thinning of pathways appeared to progress towards the frontal region by 22 GW. Previous studies reported a higher neural density in the posterior than in the anterior brain regions in primates (Cahalane et al., 2012; Charvet et al., 2015), presumably in adult mature brains, although the studies did not specify the ages of the brains. Our result suggests that a higher density of the radial tractography pathways in the posterior region contributes to the anterior-posterior gradient in neuron density.
While the regression of the radial pathways in the dorsal parietal and frontal regions continued through 21 GW, the regression of the radial pathways in the inferior frontal regions seemed to begin after 22 GW. Although radial pathways in the occipitoparietal regions also decreased in these ages, the late-onset regression in the inferior frontal regions seemed greater than in the occipitoparietal regions, which likely contributed to the relative increase in volume of the posterior radial pathways compared to those in the anterior region. At 30 GW, almost no anterior radial pathways were visible, while there were some residual radial pathways in the posterior regions. By 34 GW and 40 GW, radial pathways almost all disappeared throughout the brain. One reason why no significant difference between the anterior and posterior regions was observed in the late age group may be due to the uneven number of specimens in each age group within the late age group. For instance, in the late age group, only 30 GW (n = 2 of 8) brains revealed radial pathways remaining in the posterior regions, yet 40 GW brains (n = 4 of 8), which revealed sparse to no pathways in the posterior regions, comprised half of the specimens in this late age group. For future studies, increasing the number of brains and evening out the distribution of ages, specifically of specimens in the 23–29 GW, would be important to gain a more comprehensive understanding of the regression of radial migration pathways.
Tangential Migration Pathway
The tangential migration pathways were observed in the early age group (17–22 GW) and at the beginning of the late age group (30–34 GW), but were not visible at 40 GW. This observation is consistent with the literature describing the gradual disappearance of the GE by 34 GW (Rakić and Sidman, 1969). Many lines of evidence in primate and human studies report that a substantial percentage of cortical GABAergic neurons arise from the GE (Letinic et al., 2002; Rakic and Zecevic, 2003; Yu and Zecevic, 2011) with the peak of GABAergic neurogenesis occurring at the end of the second trimester of gestation (Petanjek et al., 2008; 15–24 GW). Similarly, in our study, the tangential migration pathways were not distinctly identifiable by tractography in brains developed past the second trimester. This finding is consistent with the literature reporting this pathway as a transient anatomic structure. However, during the early ages of development when the GE was identifiable by tractography, no significant asymmetry was observed.
While Arshad et al. (2015) claimed that interneuron generation in the GE extends beyond 35 GW (not PCW, however), we only showed tangential tractography pathways until 34 GW in the current study. A limitation of our study is a sampling gap between 34–40 GW. Thus it is possible that tractography may have been detectable past 34 GW, if samples were available. Another potential reason for these controversial results is the age grouping in the study by Arshad et al. (2015). They grouped 29–35 GW and 35–40 GW together. The only figure that they showed for the results from the 35–40 GW group Figure 6, Arshad et al. (2015) is not really convincing regarding the existence of interneuron generation at this age range. Other figures showing results only until 35 GW showed some evidence for interneuron generation, but because 35 GW was grouped with 29–34 GW, it is possible that a younger population primarily contributed to their results.
The coherent diffusivity within the GE visualized by tractography during the early fetal ages is likely attributable to the tangential migratory routes on scaffolds of corticofugal fibers. This directionality of migration has been shown in investigational studies using tracer-labeling (Denaxa et al., 2001) and fate mapping studies (Wichterle et al., 2001). These studies have shown that during the period of corticogenesis, migrating neurons follow multiple tangential routes to their destination in the developing cortex. More specifically, cells from the lateral GE migrate ventrally and anteriorly to give rise to neurons in the striatum. Cells from the medial GE have been shown to migrate dorsally and posteriorly (Wichterle et al., 2001). The temporospatial migratory pattern likely contributes to the diffusivity seen within the GE in this study. Although many studies confirming this tangential migratory pattern have been performed in mice, similar findings have also been described in humans (Hansen et al., 2013; Ma et al., 2013) and further support our tractography findings in human fetal brains.
Conclusion
In this study, we examined whether volumetric asymmetries by tractography were present in two critical neuronal migratory pathways, the radial and tangential migration pathways, from which progenitor neurons have been reported to migrate in the process of corticogenesis in developing human fetal brains. Several conclusions may be drawn from this study. First, our results confirm a postero-dorsal to antero-ventral direction of migration of the radial pathway based on the gradual disappearance of tractography pathways. By tractography, the volume of radial pathways was greater in the posterior region than in the anterior brain region at early stages of development (22 GW and younger). Second, volumetric asymmetry for only the radial pathways was statistically significant in the early age group (22 GW and younger). Third, no volumetric hemispheric asymmetry of the tangential pathway emerged, contrasting with the radial pathway. Taken together, these findings suggest hemispheric volumetric asymmetry related to neural progenitors from the radial glial pathway is established early in development, which contrasts with the GE. Describing the anatomic development and regression of these neuronal migratory pathways is critical as many neurologic and psychological disorders are described to have underlying disorganized corticogenesis. Future investigations aim to include a more comprehensive range of developmental stages to enable a more thorough description of the developmental changes, as well as comparisons of normal and diseased brains to understand what may contribute to the underlying problems in the formation of brain connectivity.
Author Contributions
ET performed the experiments. YM, JWS, and ET analyzed the data 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 Boston Children’s Hospital (BCH), NICHD (R01HD078561, R21HD069001, 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, Boston Children’s Hospital and Harvard Medical School.
Footnotes
References
Adler, W. T., Platt, M. P., Mehlhorn, A. J., Haight, J. L., Currier, T. A., Etchegaray, M. A., et al. (2013). Position of neocortical neurons transfected at different gestational ages with shRNA targeted against candidate dyslexia susceptibility genes. PLoS One 8:e65179. doi: 10.1371/journal.pone.0065179
Al-Jaberi, N., Lindsay, S., Sarma, S., Bayatti, N., and Clowry, G. J. (2015). The early fetal development of human neocortical GABAergic interneurons. Cereb. Cortex 25, 631–645. doi: 10.1093/cercor/bht254
Anderson, S. A., Eisenstat, D. D., Shi, L., and Rubenstein, J. L. R. (1997). Interneuron migration from basal forebrain to neocortex: dependence on dlx genes. Science 278, 474–476. doi: 10.1126/science.278.5337.474
Arshad, A., Vose, L. R., Vinukonda, G., Hu, F., Yoshikawa, K., Csiszar, A., et al. (2015). Extended production of cortical interneurons into the third trimester of human gestation. Cereb. Cortex doi: 10.1093/cercor/bhv074 [Epub ahead of print].
Bayatti, N., Moss, J. A., Sun, L., Ambrose, P., Ward, J. F. H., Lindsay, S., et al. (2008). A molecular neuroanatomical study of the developing human neocortex from 8 to 17 postconceptional weeks revealing the early differentiation of the subplate and subventricular zone. Cereb. Cortex 18, 1536–1548. doi: 10.1093/cercor/bhm184
Benes, F. M., and Berretta, S. (2001). GABAergic interneurons: implications for understanding schizophrenia and bipolar disorder. Neuropsychopharmacology 25, 1–27. doi: 10.1016/s0893-133x(01)00225-1
Büchel, C., Raedler, T., Sommer, M., Sach, M., Weiller, C., and Koch, M. A. (2004). White matter asymmetry in the human brain: a diffusion tensor MRI study. Cereb. Cortex 14, 945–951. doi: 10.1093/cercor/bhh055
Bystron, I., Blakemore, C., and Rakic, P. (2008). Development of the human cerebral cortex: Boulder Committee revisited. Nat. Rev. Neurosci. 9, 110–122. doi: 10.1038/nrn2252
Cahalane, D. J., Charvet, C. J., and Finlay, B. L. (2012). Systematic, balancing gradients in neuron density and number across the primate isocortex. Front. Neuroanat. 6:28. doi: 10.3389/fnana.2012.00028
Cai, L., Hayes, N. L., Takahashi, T., Caviness, V. S. Jr., and Nowakowski, R. S. (2002). Size distribution of retrovirally marked lineages matches prediction from population measurements of cell cycle behavior. J. Neurosci. Res. 69, 731–744. doi: 10.1002/jnr.10398
Caviness, V. S., Meyer, J., Makris, N., and Kennedy, D. N. (1996). MRI-based topographic parcellation of human neocortex: an anatomically specified method with estimate of reliability. J. Cogn. Neurosci. 8, 566–587. doi: 10.1162/jocn.1996.8.6.566
Chang, B. S., Piao, X., Giannini, C., Cascino, G. D., Scheffer, I., Woods, C. G., et al. (2004). Bilateral generalized polymicrogyria (BGP): a distinct syndrome of cortical malformation. Neurology 62, 1722–1728. doi: 10.1212/01.wnl.0000125187.52952.e9
Charvet, C. J., Cahalane, D. J., and Finlay, B. L. (2015). Systematic, cross-cortex variation in neuron numbers in rodents and primates. Cereb. Cortex 25, 147–160. doi: 10.1093/cercor/bht214
Chi, J. G., Dooling, E. C., and Gilles, F. H. (1977). Left-right asymmetries of the temporal speech areas of the human fetus. Arch. Neurol. 34, 346–348. doi: 10.1001/archneur.1977.00500180040008
Cusmai, R., Verrotti, A., Moavero, R., Curatolo, P., Battaglia, D., Matricardi, S., et al. (2014). Rufinamide for the treatment of refractory epilepsy secondary to neuronal migration disorders. Epilepsy Res. 108, 542–546. doi: 10.1016/j.eplepsyres.2014.01.013
D’Arceuil, H., Liu, C., Levitt, P., Thompson, B., Kosofsky, B., and de Crespigny, A. (2008). Three-dimensional high-resolution diffusion tensor imaging and tractography of the developing rabbit brain. Dev. Neurosci. 30, 262–275. doi: 10.1159/000110503
Denaxa, M., Chan, C. H., Schachner, M., Pamavelas, J. G., and Karagogeos, D. (2001). The adhesion moleceule TAG-1 mediates the migration of cortical interneurons from the ganglionic eminence along the corticofugal fiber system. Development 128, 4635–4644.
Dobyns, W. B., Andermann, E., Andermann, F., Czapanksy-Beilman, D., Dubeau, F., Dulac, O., et al. (1996). X-linked malformations of neuronal migration. Neurology 47, 331–339. doi: 10.1212/wnl.47.2.331
Glesson, J. C. (2000). Classical lissencephaly and double cortex (subcortical band heterotopia): L1S1 and doublecortin. Curr. Opin. Neurol. 13, 121–125. doi: 10.1097/00019052-200004000-00002
Gressens, P. (2006). Pathogenesis of migration disorders. Curr. Opin. Neurol. 19, 135–140. doi: 10.1097/01.wco.0000218228.73678.e1
Guerrini, R., and Filippi, T. (2005). Neuronal migration disorders, genetics and epileptogenesis. J. Child Neurol. 20, 287–299. doi: 10.1177/08830738050200040401
Habas, P. A., Scott, J. A., Roosta, A., Rajagopalan, V., Kim, K., Rousseau, F., et al. (2012). Early folding patterns and asymmetries of the normal human brain detected from in utero MRI. Cereb. Cortex 22, 13–25. doi: 10.1093/cercor/bhr053
Hansen, D. V., Lui, J. H., Flandin, P., Yoshikawa, K., Rubenstein, J. L., Alvarez-Buylla, A., et al. (2013). Non-epithelial stem cells and cortical interneuron production in the human ganglionic eminences. Nat. Neurosci. 16, 1576–1587. doi: 10.1038/nn.3541
Hansen, D. V., Lui, J. H., Parker, P. R. L., and Kriegstein, A. R. (2010). Neurogenic radial glia in the outer subventricular zone of human neocortex. Nature 464, 554–561. doi: 10.1038/nature08845
Haubensak, W., Attardo, A., Denk, W., and Huttner, W. B. (2004). Neurons arise in the basal neuroepithelium of the early mammalian telencephalon: a major site of neurogenesis. Proc. Natl. Acad. Sci. U S A 101, 3196–3201. doi: 10.1073/pnas.0308600100
Honig, L. S., Herrmann, K., and Shatz, C. J. (1996). Developmental changes revealed by immunohistochemical markers in human cerebral cortex. Cereb. Cortex 6, 794–806. doi: 10.1093/cercor/6.6.794
Hori, K., Nagai, T., Shan, W., Sakamoto, A., Taya, S., Hashimoto, R., et al. (2014). Cytoskeletal regulation by AUTS2 in neuronal migration and neuritogenesis. Cell Rep. 9, 2166–2179. doi: 10.1016/j.celrep.2014.11.045
Huang, H., Xue, R., Zhang, J., Ren, T., Richards, L. J., Yarowsky, P., et al. (2009). Anatomical characterization of human fetal brain development with diffusion tensor magnetic resonance imaging. J. Neurosci. 29, 4263–4273. doi: 10.1523/jneurosci.2769-08.2009
Jansen, A., and Andermann, E. (2005). Genetics of the polymicrogyria syndrome. J. Med. Genet. 42, 369–378. doi: 10.1136/jmg.2004.023952
Keverne, E. B. (2014). Significance of epigenetics for understanding brain development, brain evolution and behaviour. Neuroscience 264, 207–217. doi: 10.1016/j.neuroscience.2012.11.030
Kolasinski, J., Takahashi, E., Stevens, A. A., Benner, T., Fischi, B., Zöllei, L., et al. (2013). Raidal and tangential neuronal migration pathways in the human fetal brain: anatomically distinct patterns of diffusion MRI coherence. Neuroimage 79, 412–422. doi: 10.1016/j.neuroimage.2013.04.125
Kostovic, I., and Vasung, L. (2009). Insights from in vitro fetal magnetic resonance imaging of cerebral development. Semin. Perinatol. 33, 220–233. doi: 10.1053/j.semperi.2009.04.003
Kovalev, V. A., Kruggel, F., and von Cramon, D. Y. (2003). Gender and age effects in structural brain asymmetry as measured by MRI texture analysis. Neuroimage 19, 895–905. doi: 10.1016/s1053-8119(03)00140-x
Lammens, M. (2000). Neuronal migration disorders in man. Eur. J. Morphol. 38, 327–333. doi: 10.1076/ejom.38.5.0327
Letinic, K., and Rakic, P. (2001). Telencephalic origin of human thalamic GABAergic neurons. Nat. Neurosci. 4, 931–936. doi: 10.1038/nn0901-931
Letinic, K., Zoncu, R., and Rakic, P. (2002). Origin of GABAergic neurons in the human neocortex. Nature 417, 645–649. doi: 10.1038/nature00779
Leventer, R. J., Mills, P. L., and Dobyns, W. B. (2000). X-linked malformations of cortical development. Am. J. Med. Genet. 97, 213–220. doi: 10.1002/1096-8628(200023)97:3<;213::aid-ajmg1039>3.0.co;2-w
Liu, Y., Balériaux, D., Kavec, M., Metens, T., Absil, J., Denolin, V., et al. (2010). Structural asymmetries in motor and language networks in a population of healthy preterm neonates at term equivalent age: a diffusion tensor imaging and probabilistic tractography study. Neuroimage 51, 783–788. doi: 10.1016/j.neuroimage.2010.02.066
Ma, T., Wang, C., Wang, L., Zhou, X., Tian, M., Zhang, Q., et al. (2013). Subcortical origins of human and monkey neocortical interneurons. Nat. Neurosci. 16, 1588–1597. doi: 10.1038/nn.3536
Mackay, C. E., Barrick, T. R., Roberts, N., DeLisi, L. E., Maes, F., Vandermeulen, D., et al. (2003). Application of a new image analysis technique to study brain asymmetry in schizophrenia. Psychiatry Res. 124, 25–35. doi: 10.1016/s0925-4927(03)00088-x
Marín, O., and Rubenstein, J. L. (2001). A long, remarkable journey: tangential migration in the telencephalon. Nat. Rev. Neurosci. 2, 780–790. doi: 10.1038/35097509
McIntosh, A. M., Muñoz Maniega, S., Lymer, G. K. S., McKirdy, J., Hall, J., Sussmann, J. E. D., et al. (2008). White matter tractography in bipolar disorder and schizophrenia. Biol. Psychiatry 64, 1088–1092. doi: 10.1016/j.biopsych.2008.07.026
McKinstry, R. C., Mathur, A., Miller, J. H., Ozcan, A., Snyder, A. Z., Schefft, G. L., et al. (2002). Radial organization of developing preterm human cerebral cortex revealed by non-invasive water diffusion anisotropy MRI. Cereb. Cortex 12, 1237–1243. doi: 10.1093/cercor/12.12.1237
Mori, S., Crain, B. J., Chacko, V. P., and van Zijl, P. C. (1999). Three-dimensional tracking of axonal projections in the brain by magnetic resonance imaging. Ann. Neurol. 45, 265–269. doi: 10.1002/1531-8249(199902)45:2<265::aid-ana21>3.0.co;2-3
Muraki, K., and Tanigaki, K. (2015). Neuronal migration abnormalities and its possible implications for schizophrenia. Front. Neurosci. 9:74. doi: 10.3389/fnins.2015.00074
Petanjek, Z., Berger, B., and Esclapez, M. (2009a). Origins of cortical GABAergic neurons in the cynomolgus monkey. Cereb. Cortex 19, 249–262. doi: 10.1093/cercor/bhn078
Petanjek, Z., Kostović, I., and Esclapez, M. (2009b). Primate-specific origins and migration of cortical GABAergic neurons. Front. Neuroanat. 3:26. doi: 10.3389/neuro.05.026.2009
Petanjek, Z., Dujmović, A., Kostović, I., and Esclapez, M. (2008). Distinct origin of GABA-ergic neurons in forebrain of man, nonhuman primates and lower mammals. Coll. Antropol. 32, 9–17.
Platt, M. P., Adler, W. T., Mehlhorn, A. J., Johnson, G. C., Wright, K. A., Choi, R. T., et al. (2013). Embryonic disruption of the candidate dyslexia susceptibility gene homolog Kiaa0319-like results in neuronal migration disorders. Neuroscience 248, 585–593. doi: 10.1016/j.neuroscience.2013.06.056
Price, J., and Thurlow, L. (1988). Cell lineage in the rat cerebral cortex: a study using retroviral-mediated gene transfer. Development 104, 473–482.
Radonjić, N. V., Ayoub, A. E., Memi, F., Yu, X., Maroof, A., Jakovcevski, I., et al. (2014). Diversity of cortical interneurons in primates: the role of the dorsal proliferative niche. Cell Rep. 9, 2139–2151. doi: 10.1016/j.celrep.2014.11.026
Rakic, P. (1972). Mode of cell migration to the superficial layers of fetal monkey neocortex. J. Comp. Neurol. 145, 61–83. doi: 10.1002/cne.901450105
Rakic, P. (1988). Specification of cerebral cortical areas. Science 241, 170–176. doi: 10.1126/science.3291116
Rakić, P., and Sidman, R. L. (1969). Telancephalic origin of pulvinar neurons in the fetal human brain. Z. Anat. Entwicklungsgesch. 129, 53–82. doi: 10.1007/bf00521955
Rakic, S., and Zecevic, N. (2003). Emerging complexity of layer I in human cerebral cortex. Cereb. Cortex 13, 1072–1083. doi: 10.1093/cercor/13.10.1072
Schenker, N. M., Hopkins, W. D., Spocter, M. A., Garrison, A. R., Stimpson, C. D., Erwin, J. M., et al. (2010). Broca’s area homologue in chimpanzees (Pan troglodytes): probabilistic mapping, asymmetry and comparison to humans. Cereb. Cortex 20, 730–742. doi: 10.1093/cercor/bhp138
Schmahmann, J. D., Pandya, D. N., Wang, R., Dai, G., D’Arceuil, H. E., de Crespigny, A. J., et al. (2007). Association fibre pathways of the brain: parallel observations from diffusion spectrum imaging and autoradiography. Brain 130, 630–653. doi: 10.1093/brain/awl359
Snyder, P. J., Bilder, R. M., Wu, H., Bogerts, B., and Lieberman, J. A. (1995). Cerebellar volume asymmetries are related to handedness: a quantitative MRI study. Neuropsychologia 33, 407–419. doi: 10.1016/0028-3932(94)00125-9
Song, J. W., Mitchell, P. D., Kolasinski, J., Ellen Grant, P., Galaburda, A. M., and Takahashi, E. (2014). 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
Takahashi, E., Song, J. W., Folkerth, R. D., Grant, P. E., and Schmahmann, J. D. (2013). Detection of postmortem human cerebellar cortex and white matter pathways using high angular resolution diffusion tractography: a feasibility study. Neuroimage 68, 105–111. doi: 10.1016/j.neuroimage.2012.11.042
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
Vishwas, M. S., Chitnis, T., Pienaar, R., Healy, B. C., and Grant, P. E. (2010). Tract-based analysis of callosal, projection and association pathways in pediatric patients with multiple sclerosis: a preliminary study. AJNR Am. J. Neuroradiol. 31, 121–128. doi: 10.3174/ajnr.A1776
Volpe, J. J. (2009). Brain injury in premature infants: a complex amalgam of destructive and developmental disturbances. Lancet Neurol. 8, 110–124. doi: 10.1016/s1474-4422(08)70294-1
Walsh, C., and Cepko, C. L. (1988). Clonally related cortical cells show several migration patterns. Science 241, 1342–1345. doi: 10.1126/science.3137660
Wedeen, V. J., Wang, R. P., Schmahmann, J. D., Benner, T., Tseng, W. Y. I., Dai, G., et al. (2008). Diffusion spectrum magnetic resonance imaging (DSI) tractography of crossing fibers. Neuroimage 41, 1267–1277. doi: 10.1016/j.neuroimage.2008.03.036
Wegiel, J., Kuchna, I., Nowicki, K., Imaki, H., Wegiel, J., Marchi, E., et al. (2010). The neuropathology of autism: defects of neurogenesis and neuronal migration and dysplastic changes. Acta Neuropathol. 119, 755–770. doi: 10.1007/s00401-010-0655-4
Wichterle, H., Turnbull, D. H., Nery, S., Fishell, G., and Alvarez-Buylla, A. (2001). In utero fate mapping reveals distinct migratory pathways and fates of neurons born in the mammalian basal forebrain. Development 128, 3759–3771.
Wonders, C. P., and Anderson, S. A. (2006). The origin and specification of cortical interneurons. Nat. Rev. Neurosci. 7, 687–696. doi: 10.1038/nrn1954
Xu, G., Takahashi, E., Folkerth, R. D., Haynes, R. L., Volpe, J. J., Grant, P. E., et al. (2014). Radial coherence of diffusion tractography in the cerebral white matter of the human fetus: neuroanatomic insights. Cereb. Cortex 24, 579–592. doi: 10.1093/cercor/bhs330
Yu, X., and Zecevic, N. (2011). Dorsal radial glial cells have the potential to generate cortical interneurons in human but not in mouse brain. J. Neurosci. 31, 2413–2420. doi: 10.1523/JNEUROSCI.5249-10.2011
Keywords: development, radial migration, ganglionic eminence, human, diffusion imaging, tractography
Citation: Miyazaki Y, Song JW 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
Received: 14 September 2015; Accepted: 02 January 2016;
Published: 25 January 2016.
Edited by:
Hao Huang, University of Pennsylvania, USAReviewed by:
Marina Bentivoglio, University of Verona, ItalyGavin John Clowry, Newcastle University, UK
Copyright © 2016 Miyazaki, Song and Takahashi. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution and 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