MRI and M/EEG studies of the White Matter Development in Human Fetuses and Infants: Review and Opinion

1. Post-mortem identification of the developing bundles

Because this section is mainly a prerequisite to understand the following in vivo findings (I.2.), we here summarize the developmental descriptions that have been reported on post-mortem specimen so far, whatever the investigation technique: histological examinations, anatomical and diffusion MRI (see [7–10] for further details). These approaches provide complementary information, each one being intrinsically limited to identify the tracts precisely (particularly MRI as detailed in I.3.). Consequently this global report strongly relies on the current state of technological art.

Before the formation of stable postnatal connectivity, transient fetal patterns of connectivity are observed [8–10]. During the early fetal period (9-15PCW), major fiber pathways are growing and pathfinding within the intermediate zone. Some afferent (e.g. thalamo-cortical fibers in the internal capsule) and efferent (cortico-subcortical) fibers are already observed. Before penetrating into the cortical plate, thalamo-cortical fibers make connections with neurons of the subplate zone, a crucial structure that develops between 13 and 15PCW [11]. Several limbic bundles (e.g. fornix, cingulum) are also identified [9], and their entire trajectories are visible at 17PCW with diffusion imaging andtractography [12].

During the midfetal period (15-23PCW), major efferent fibers penetrate their targets in the striatum, the pons and the spinal cord. Periventricular crossroads are developing, notably with the growth of cortico-cortical callosal and associative fibers within the intermediate zone [13]. The corpus callosum, the uncinate and inferior longitudinal fasciculi are detected with DTI [14]. The radial coherence of white matter observed at these ages with HARDI likely reflects a mixture of radial glial fibers, penetrating blood vessels, and radial axons [15]. In specimens younger than 20PCW, the radial migration pathway seems more developed posteriorly than anteriorly [16]. Recently, axonal tracts were also identified with digitalization and structure tensor analysis of histological sections stained with antibody against neural cell adhesion molecule (NCAM) [17]. At 20PCW, the intermediate zone contains different populations of fibers: bundles of fibers with a predominant tangential orientation, arching in a C-shape from the exit of the internal capsule around and parallel to the lateral ventricle toward the corpus callosum, and a more diffusely organized fraction of fibers with radial orientation [17].

During the early preterm period (24–28PCW), afferent fibers that were waiting in the subplate start making connections within the cortical plate in most brain regions, leading to the establishment of the permanent connectivity with neurons of the future cortical layer IV. Nevertheless a transient fetal circuitry still exists within the subplate [18], and this compartment reaches its maximum volume and thickness around 30PCW [19]. Limbic cortico-cortical connections are well developed in the cingulate, entorhinal and hippocampal cortices [20]. In the lateral neocortex, the fronto-occipital and uncinate fascicles are visible [21, 22].

During the late preterm period (29–34 PCW), long associative and commissural bundles are quickly developing, likely originating from pyramidal neurons of cortical layer III [23]. Short cortico-cortical fibers are growing and entering the cortex through the remnant subplate [21, 24].

After term birth, U-fibers develop [22], and cortico-cortical connectivity is reorganized by several processes including the development of dendritic arborization, the overproduction of synapses and dendritic spines and their later elimination [25–27]. Pruning of axonal fibers like callosal connections [28] also probably extends until the end of the first postnatal year [11, 22, 29].

2. In vivo identification of bundles with diffusion MRI

In vivo diffusion MRI of fetal and preterm brains has confirmed these post-mortem studies but with lower spatial resolution. In utero, the pyramidal tract, the splenium and genu of the corpus callosum are observed early on from 18 gestational weeks (GW) [30–32]. The uncinate and inferior fronto-occipital fascicles are depicted as early as 20GW, while the inferior longitudinal fasciculus, cingulum and fornix can be visualized in 3D in older fetuses during the third trimester [33]. In preterms, subcortical projection and association bundles are also identified[34, 35].

In the recent years, the network architecture of the preterm brain has also been detailed based on connectivity matrices, by measuring for all pairs of brain regions, the degree of anatomical connections (i.e. the number of reconstructed streamlines, eventually weighted by DTI anisotropy as a local measure of the connections strength and integrity). Already at 30 post-menstrual weeks (PMW), such a structural connectome demonstrates a small-world modular organization like in the adult brain [36]. Cortical hubs are highly connected to form a “rich club” that may underlie early network communication [37]. Until term-equivalent age, the small-world topology further increases [36], and the network structure gets more clustered between 27 and 45PMW, with the apparent development of many connections [38]. While these approaches of structural connectivity can provide whole-brain descriptions of networks without any a priori anatomical assumption on specific pathways, their biological interpretation throughout development should remain cautious relatively to plausible neural connectivity and myelination (see [39] for a more general discussion on the anatomical reality of small-world architecture reported in brain imaging studies).

During early infancy, although weakly myelinated, all major white matter bundles are identified with diffusion imaging and tractography: commissural bundles of the corpus callosum (genu, body and splenium), projection bundles (cortico-spinal and spino-thalamic tracts, optic radiations, anterior limb of the internal capsule), limbic bundles (fornix and cingulum) and associative bundles (external capsule, uncinate, arcuate, superior and inferior longitudinal fascicles) [40, 41] (Figs. 1, 2, 4a-d). The trajectory and morphology of these bundles remain stable, as revealed by a longitudinal study, between birth and 2 years of age [42]. Short-range connections, such as the ones projecting above and below the superior temporal sulcus, are also depicted in the infant brain by 3T-MR imaging at high spatial resolution [43] (Fig. 4c). At the whole brain level, an increase in global efficiency and integration, and a decrease in segregation are observed on the structural connectome during the first two post-natal years [44].

3. Caveats

Several issues should be kept in mind when diffusion imaging and tractography are used to decipher the developing organization of white matter fibers. First, estimating the real axonal trajectories is not so simple, and notably during the fetal life to disentangle them from the radial glia fibers. The trajectory of reconstructed fibers is highly dependent on the acquisition parameters (e.g. the number of directions), the diffusion model and the tractography approach, and results may not be comparable across acquisition protocols, across studies, or even across ages due to the growing brain size, the myelination and density of axons, and the complex crossings of bundles with other structures (e.g. other bundles, glial fibers). DTI and deterministic tractography can highlight relatively simple trajectories, but only more elaborated models (e.g. HARDI) and algorithms (e.g. deterministic with regularization, probabilistic) can highlight accurate and realistic tracts (see Fig. 2a,b showing reconstructions of the cortico-spinal tract). Furthermore, a bundle may not be dissected, particularly in the developing brain, because of an insufficient spatial resolution compared to its small size, or because of its low maturation compared with surrounding fibers. When bundles are crossing, the situation is even worse, and exact trajectories may hardly be delineated. Besides, diffusion imaging is not sensitive to the functional direction of fibers that is from the neural soma to the periphery or from region A to B. Thus, to which direction the fiber is growing cannot be directly inferred at the early stages of brain development. One can only try to assess whether connections exist or not between regions A and B.

Second, there is no simple way so far to measure fibers overproduction and pruning, yet a major developmental process. Throughout infancy, macroscopic changes in bundles volume are mostly driven by the increase in axonal diameter and fiber myelination, rather than by fiber loss. In the infrequent event that fibers are identified at an early stage but not at a more mature stage, one may wonder whether it reflects a genuine pruning phenomenon or artefactual differences in tract reconstructions (see Figs. 1 and 2 for comparative reconstructions in infants and adults). Similarly as for the tract trajectories, the number of reconstructed streamlines strongly depends on acquisition and post-processing parameters. It is impacted by effective spatial resolution relatively to brain size, by white matter maturation as detailed below, etc. Thus this measure should be considered with caution notably in connectivity matrices, and it may provide little information on the real axonalcount.

Third, the bundle maturational stage highly affects the tracts delineation. On the one hand, the tight organization of fibers in compact and highly coherent bundles leads to intrinsic anisotropy even in the absence of myelin [45], for instance in the corpus callosum of preterm newborns [46]. Studies in rat pups have shown that the first evidence of anisotropy precedes initial myelin [47], and that this early anisotropy may be related to sodium-channel activity [48]. But on the other hand, myelination highly increases the fibers anisotropy (see part II), and thus the contrast between the bundle, surrounding tissue and crossing fibers. This is particularly an issue when one aims to assess whether connections exist or not in the developing brain. It should be kept in mind that reconstructed fibers definitely remain a rough representation of axonal fibers. Recently, Li and colleagues [49] described with DTI an increase in the density of fibers connecting the cortex (projecting to or originating from it) during the first post-natal year. Rather than a recognized biological phenomenon, their observation probably relies on the maturation of fibers that already exist at birth but could not be tracked. In the same way, weakly mature bundles may appear shorter on tractography reconstructions in infants than in adults, just because of the low anisotropy at the distant extremity related to the myelination progression from the cell soma to the axonal end. This phenomenon may lead to erroneous interpretation on the presence/absence of a pathway. For instance, the arcuate fasciculus may wrongly seem not to fully connect the posterior temporal and inferior frontal regions in newborns [50, 51], as a result from an artefact due to the low myelination of its frontal extremity [43] (Fig. 2c,d). The reconstruction becomes even more complex at places where bundles maturing over different time-periods do cross. It is for example the case for the crossing of the arcuate fasciculus and cortico-spinal tract at the level of the corona radiata (Fig. 2e). The same phenomena may explain why some late-maturing associative bundles (e.g. the superior longitudinal fasciculus) display changes in fiber orientations during the first post-natal months [52] though they have grown prenatally. Finally, these caveats (e.g. missing or shorter connections) also highly impact connectome approaches. In regards to these different issues, section II.5 summarizes the bottlenecks to be solved in order to accurately track bundles of the immature brain despite fiber pruning, bundles crossing and asynchronous maturation.

1. Current knowledge based on histological examinations

Concurrently and subsequently to the organization of general neural connectivity, the white matter bundles progressively mature and become functionally more efficient through myelination (i.e. myelin formation around axons). This process includes several steps, and proceeds from the neuron body to the periphery [53]. It occurs in the human brain from the second part of pregnancy to the end of adolescence, with a peak during the first post-natal year, and it develops at different ages and rates depending on the bundles and the networks [54–56]. Thus the myelination progression varies across cerebral regions [56–62]: it follows a caudo-rostral gradient, progresses from the center to the periphery, and it occurs earlier and faster in sensory pathways (somatosensory, vision, audition) than in motor ones, and in projection fibers than in associative ones. Left-right asynchronies are also observed [43, 63], which are interesting to consider in relation to the functional asymmetries of the motor, linguistic and social cognition networks.

To classify white matter regions according to their myelination patterns, eight sub-groups have been proposed. Based on post-mortem examinations, they integrate whether myelin is present (A) or not (B) at birth, and the age periods at which mature myelin is observed (stages 1–4) [57]. For instance, the cortico-spinal tract in the midbrain and pons, the posterior limb of the internal capsule and the central corona radiata belong to sub-group A1; optic and auditory radiations, the corpus callosum body and splenium belong to sub-group B1; the anterior limb of the internal capsule, the external capsule and cingulum belong to sub-group B2, etc.

2. In vivo quantification based on MRI parameters

Since a few years, several MRI parameters (Fig. 3) have been used to quantify the myelination process in the developing white matter [7]. As detailed below, intense changes are measured after term-birth and during the first post-natal months (Fig. 3 right row), then the rate of modifications slows down during toddlerhood and more drastically thereafter until young adulthood.

With myelination, water content in voxels decreases, thus the proton density decreases, while the macromolecule tissue volume [64] increases [65]. This mechanism, together with changes in water molecules compartmentalization [66] and increase of protein and lipid contents [67, 68], lead to decreases in T1 and T2 relaxation times with age [69–71] (Fig. 3a,b), with different time courses between T1 and T2 [72]. Indeed, T1 shortening starts first during the “pre-myelinating” state, while T2 shortening occurs later on with the chemical maturation of the myelin sheath [56, 67, 73]. Myelinated white matter regions may already be found on T1w and T2w images during the preterm period, mainly from 36w GA on [74], but changes are mainly seen after term and vary across the brain because of the asynchronous myelination process [40, 65].

DTI parameters capture some aspects of this maturational pattern [75, 76]. During the preterm period, diffusivities decrease, while anisotropy increases in most white matter regions [77] except at cross-roads locations [78]. During infancy and childhood, transverse diffusivity decreases more than longitudinal diffusivity [42, 79, 80] (Fig. 3d,e), and anisotropy keeps increasing (Fig. 3c) except where the crossing bundles are at different maturational stages creating an ambiguous signal [7]. These parameters might be sensitive to distinct microstructural mechanisms such as the proliferation of glial cell bodies and extension of oligodendroglial processes during the first stage of myelination (“pre-myelination”) [81, 82], and the ensheathment of oligodendroglial processes around the axons during the “true” myelination process [79]. Among DTI parameters, transverse diffusivity is probably the best marker of myelination [83–85]. Aside from the myelination effects, anisotropy and longitudinal diffusivity are more dependent on the tract intrinsic architecture (e.g. fibers compactness, crossings) than transverse diffusivity. Indeed, the variability in anisotropy and longitudinal diffusivity observed across bundles seem relatively stable during infancy and adulthood, despite important age-related changes [7, 43, 86]. Like for T1 and T2, DTI parameters demonstrate an asynchronous evolution across bundles [79], even within a functional network [43] (Fig. 4c). A central-to-peripheral and posterior-to-anterior gradient of maturation is observed during the preterm period [77], as well as a posterior-to-anterior gradient duringchildhood [80].

The progression of myelination along tracts can be exquisitely studied with tractography-based quantification of DTI parameters, as recently reported for main bundles during toddlerhood [87], childhood and adolescence [88]. For instance, during the preterm period, the cortico-spinal tract starts maturing at the level of the internal capsule [89]. During infancy, the optic radiations show an early central-to-peripheral wave of maturation, followed by a peripheral-to-central wave, which might correspond to thalamo-cortical and cortico-thalamic fibers respectively [86]. Because DTI parameters are also affected by the tract shape, compactness, and by the crossings and partial volume effects with neighboring structures, such analyses require a normalization of the values by the same measures done at the mature stage to disentangle maturation per se from these other factors [7, 43, 79, 86].

Several other MRI parameters can also be measured to quantify the myelin amount in white matter. Magnetic susceptibility decreases during childhood with maturation [90]. Magnetization transfer imaging (MTI) is sensitive to hydrogen atoms bound to macromolecules, such as proteins and lipids found in myelin [91], and magnetization transfer ratio (MTR) increases with myelination [68], showing asynchrony across white matter regions [92]. However, during the preterm period, MTR is unexpectedly high in the unmyelinated callosal fibers [81], and it decreases in the frontal intermediate zone where cross-roads develop [78]. These observations suggest that MTR is not only sensitive to myelin-associated macromolecules, but also to the macromolecular density of axonal cytoskeleton components such as microtubules and neurofilaments. A recent study in a cuprizone mouse model of demyelination has shown that among MRI parameters (T1/T2, DTI, MTI), the bound pool fraction (f) from MTI is the best indicator of the myelin sheath fraction, while T1 relates to fraction of myelinated axons, and longitudinal diffusivity relates to the fraction of non-myelinated cells; these two latter parameters thus characterize the tissue micro-architecture rather thanmyelination [93].

3. Refined quantification based on multi-compartment models

Recently, several modelling approaches based on relaxometry or diffusion data have been proposed to characterize distinct pools of water molecules inside each voxel. These pools can be distinguished on the basis of their different relaxometry or diffusion characteristics. Such decomposition is supposed to provide more accurate information on the tissue microstructure than considering univariate parameters at the voxel level. Nevertheless, it requires the acquisition of many more data to estimate the pools characteristics, limiting its actual use in fetuses and infants.

On the one hand, analyses based on multi-component relaxation (MCR) assume that at least two or three pools of water contribute to the MR signal sensitized to T1 and T2 relaxometries [94]: water located within the myelin sheath (with relatively short T1 and T2 relaxation times), intra-axonal, intra-cellular and interstitial water (i.e. water outside of the myelin sheath, with intermediate T1 and T2), and free water (with long T1 and T2, close to values of cerebro-spinal fluid). Different methods have been used so far in the developing brain, such as multi-component T2 relaxometry [95] or multi-component driven equilibrium single pulse measures [96]. The main potential of these approaches is to estimate the volume fraction of water related to myelin (often named the “myelin water fraction” MWF). This fraction drastically increases with age in the white matter (Fig. 3f), shows a spatio-temporal progression that is coherent with histological studies of myelination [72, 97], and follows a nonlinear growth pattern during infancy and toddlerhood [98]. To reduce the MRI acquisition time required for the estimation of MWF maps in infants, we recently proposed an approach based on the assumption that T1 and T2 characteristics of water pools remain stable throughout development and can be fixed to adult estimates, and that only the fraction changes of the different pools account for the known sequence of T1 and T2 changes during development [99]. Such an hypothesis may be a good approximation given previous description of stable myelin chemical properties during maturation [73].

On the other hand, different multi-component approaches have been proposed to analyze the MR diffusion signal acquired with multiple shells (i.e. sensitized to multiple b-values) and with multiple gradient directions. Designed for the adult brain, these complex models (e.g. CHARMED [100] or NODDI [101]) enable to characterize the tissue microstructural properties and eventually the spatial distribution of fibers within white matter voxels, by modeling distinct water pools (e.g. water with hindered diffusion in the extra-cellular space, water with restricted diffusion in the intra-cellular compartment). Recent studies have tested these approaches in the developing brain by fixing the pools characteristics (diffusivities) to adult estimates. This is a reasonable hypothesis given that intra- and extra-cellular axial diffusivities, estimated with diffusion kurtosis imaging (DKI), appear stable throughout development [102]. In newborns, the intra-neurite volume fraction from NODDI may inform on the maturation of white matter fibers, while the index of neurite orientation dispersion seems to reflect the presence of fiber crossings and fanning [103]. Thus these parameters may help distinguish between bundles with similar cellular structure but different myelination (e.g posterior vs anterior limb of the internal capsule), or reciprocally with similar maturation but different fiber microstructural organization (e.g. external capsule vs periventricular crossroads) [103]. During the first three post-natal years, intra-axonal water fraction and tortuosity of the extra-axonal space increase, with asynchronous patterns across distinct white matter regions such as the corpus callosum genu and splenium, and the posterior limb of the internal capsule [102]. Actually, during childhood and adolescence, the increase in DTI anisotropy observed in most bundles seems to relate to the increase in NODDI neurite density rather than on changes in orientation dispersion [104]. Despite these approaches look very promising to precisely characterize the tissue microstructure, diffusion imaging only provides indirect information on the myelin amount: indeed no pool of water related to myelin can be taken into account in diffusion models because of too short T2 characteristics.

4. Refined quantification based on multi-parametric approaches

Despite some correlations between specific pairs of parameters and over specific developmental periods [7], MRI parameters (T1, T2, DTI, MTR, MWF, NODDI, etc.) do capture different tissue properties, and provide complementary information on maturation mechanisms. Thus they show different age-related trajectories (for recent example over the lifespan [65]), which makes it possible to differentiate white matter structures based on their maturation patterns (for recent example during the preterm period [78, 81]). A clever combination or integration of these parameters may thus provide a unique opportunity to infer the bundles microstructural organization and myelination, and to put in relation networks maturation and infants’ learning.

In complement to myelin-sensitive and diffusion imaging techniques, a composite framework has been proposed recently to estimate the g-ratio [105], i.e. the ratio between the inner and outer diameters of the myelin sheath, which is supposed to be a fundamental characteristic of white matter fibers. The approach relies on MTI to estimate the myelin content and on NODDI to assess the intra-axonal fraction. It provides accurate estimations of g-ratios in the corpus callosum compared with histological examinations in the monkey brain, and looks promising to detect demyelinating lesions in adult patients with multiple sclerosis [105]. Recently, similar approaches using MWF instead of MTR, have been tested to quantify the white matter myelination in the developing brain of newborns [106], infants and children [107]. Because g-ratio depends on at least two different MRI techniques, this parameter may be even more sensitive to acquisition protocols than previous parameters, making comparisons across studies difficult. Besides, this approach only provides a rough estimate of the complex processes involved in the development of white matter microstructure and myelination. Among unrevealed aspects, g-ratio is a voxel-wise aggregate index that, like other MRI parameters, merges the properties of distinct fibers: e.g. fibers to be pruned do not myelinate but probably may have similar properties of diameter and trajectory than myelinating fibers; crossing fibers may show different g-ratios and age-related changes according to asynchronous myelination, etc.

In the last years, a few other multi-parametric approaches have shown potential to quantify the progression of white matter maturation throughout infancy. Let us remember that MRI parameters vary across both bundles and white matter regions also in the adult brain, in relation to fibers geometry, compactness, crossings, etc. Highlighting maturational effects in the developing brain thus requires either considering the developmental trajectories toward adulthood to evaluate the asymptotes of maturation, or normalizing infant measurements by adult references. Recently, modelling the evolution of proton density, T1w and T2w signals, or DTI diffusivities has provided quantitative absolute rates of change and relative time shifts across regions [108]. Combining the time trajectories of DTI parameters has also accurately highlighted their evolution across bundles [109]. Besides, clustering brain voxels over a group of infants on the basis of T1, T2 and DTI parameters, has uncovered four types of white matter regions with different compactness and maturation [110]: the spatial distribution of these clusters was anatomically relevant while the approach was free from any anatomical hypothesis. Finally, a maturational distance to the adult stage has been proposed to summarize changes in T1, T2 and DTI diffusivities during early infancy while taking into account their possible correlations [40]. In addition to outperforming univariate approaches, this approach enables the quantitative timing of maturation asynchrony (in weeks) across projection, association, limbic and callosal bundles, showing that most intense changes occur during the first post-natal year(Fig. 4a,b).

5. Some reflections for the accurate exploration of bundles architecture and maturation

To sum up, several MRI techniques have been recently proposed to explore the organization and maturation of white matter networks in the fetal and infant brain. Nevertheless, several concerns have been raised on the correspondence between these measures and the biological mechanisms occurring during development, mainly because these mechanisms were not considered when the techniques were first imagined and implemented for the adult brain. So far major developmental processes (e.g. axonal overproduction-pruning) remain unexplored, and complex situations (e.g. asynchronous myelination of crossing fibers) are insufficiently considered.

Regarding the development of white matter architecture, a first bottleneck relies on the accuracy of bundles reconstruction. As detailed previously, current tractography approaches are directly affected by the maturational stage. They provide limited information on fibers growth, pruning and crossings, and often mix these mechanisms at the voxel level. Yet, the known patterns of maturation might help to disentangle these complex situations: myelination progresses from neural soma to periphery, pruning and stabilized fibers have distinct myelination, fibers belonging to different bundles mature at different times, etc. Considering these maturational characteristics and asynchronies might facilitate the reconstruction of bundles and the differentiation between separate pools of fibers. A reverse bottleneck relies on the specificity of maturation measures within each bundle. Although promising, most univariate, multi-compartment or multi-parametric approaches provide aggregate measures at the voxel level, whereas a tract-level resolution should be achieved even in the complex places of bundles crossings.

Understanding the development of white matter networks might require studying both bundles architecture and maturation jointly, by modelling different pools of fibers with different spatial architectures, microstructural and myelination characteristics. Some aspects might be dealt with at the acquisition step. Recently, complex MR sequences have been implemented for the adult brain to explore fibers with different orientations based on diffusion MRI while measuring their myelin properties with T1 relaxometry [111]. Nevertheless, technological challenges remain to be addressed for infant imaging, for instance in terms of acquisition time, gradient and radiofrequency exposure. Alternatively, some implementations might be achieved in the post-processing steps, like constraining the quantification procedure by the bundles spatial trajectory (i.e. considering a certain degree of regularization of MRI parameters quantified along fibers).

1. Visual system

The visual modality develops mainly after birth, and infants’ visual capacities drastically improve in a few months. Early visual cortical responses in newborns have been mostly studied using EEG recordings during pattern reversal stimulation. In adults, visual evoked potentials (VEPs) recorded in occipital regions at short latencies consist of a positive peak at around 100 ms (P100), surrounded by two negative peaks at around 70 and 145 ms (N70 and N145) [118]. At term birth, the small positive component P1 (∼P100 in adults) is detected at a latency that decreases strongly and quickly with age [119–121], from around 260 ms in neonates to around 110–120 ms at 12–14 weeks of age, depending on the patterns size [122] (Fig. 5a). The earlier negative component N1 (∼N70) develops in general after 2 weeks of age, and reaches a latency below 100 ms by 5 weeks of age [123, 124]. The second negative component N2 (∼N145) is observed from 2 months of age [118]. P1 and N2 latencies reach their adult values around 7–12 months, whereas N1 latency decreases until 24–66 months [118].

These observations on VEPs responses suggest progressive increase in conduction velocity in afferent visual pathways related to fiber myelination [125], together with other maturational processes (e.g. retina and cortical development). A few recent studies have compared these functional changes to the microstructural evolution of visual bundles throughout development, focusing on the latency of P1 and its transfer from contralateral to ipsilateral hemisphere in case of lateral stimuli. During infancy, the increase in P1 conduction velocity has been related to the maturation of optic radiations (transverse diffusivity measured with DTI) [86]. This relation was specific (i.e. not observed for other white matter bundles) and not explained by intra-individual differences in infants’ age. This study has been recently confirmed for visual lateral stimuli, and extended to cortico-cortical connections: the conduction velocity of P1 inter-hemispheric transfer was related to the maturation of the visual fibers of the corpus callosum [126]. In another study in 6–12-year-old children, the P1 latency from visual evoked fields (VEFs) was inversely related to the white matter integrity (FA measured with DTI) in visual and motor association regions [127]. Again this effect was independent of age differences across children, and it was not observed for N1 and N2 latencies.

Besides developmental studies, this issue of structure-function relationships has been also addressed in adults to characterize visual processing variability across subjects. A study combining VEFs and DTI has shown an inverse relationship between P1 latency and the microstructure (FA) of bilateral parietal and right lateral frontal white matter [128], so with a different spatial dependence than for children [127]. Based on VEPs and DTI, the inter-hemispheric transfer time of visual information for P1 (but not N1) has been inversely related to FA in the visual callosal fibers [129]. Nevertheless, these latter studies have not taken into account the inter-individual variations in the bundles length related to brain size and shape, which probably impact latency differences in a non-negligible way. Using diffusion imaging (AxCaliber approach), it has been proposed that the conduction velocity for N2 inter-hemispheric transfer was related to the axon diameter of visual callosal fibers [130]. Whereas such a demonstration in the human brain is very exciting, this study has been criticized in several aspects [131]: notably the estimated transfer times were very short (∼5 ms), N2 is a relatively late visual component, and earlier components would have been preferable in order to target activity in the primary visual cortex.

2. Somatosensory and auditory systems

So far, very few studies have related anatomical and functional maturation in the somatosensory and auditory modalities, despite intense changes in electrophysiological responses throughout infancy. Cortical somatosensory evoked potentials (SEPs) have been measured with EEG following electrical stimulation of the median nerve, at an intensity just above the motor threshold. In the adult brain, a negative peak (N20) and a positive peak (P25) are measured over centro-parietal electrodes at around 20 and 25 ms [132], and are supposed to reflect activity in the primary somatosensory cortex (S1) [133]. The SEPs development is intense in infants. At term birth, N20 is observed at around 25 ms [134], and its latency strongly decreases as a function of age during infancy [134, 135]. P25 becomes the major component by 2-3 weeks of age [136]. More recently, somatosensory evoked magnetic fields (SEFs) have been measured with MEG. In adults, two early responses with opposite polarity are observed following median nerve or tactile stimulation (responses noted N20 m and P30 m [137], or M30 and M50 [138]). In full-term newborns, two responses of the same polarity (M30 and M60) are detected with generators plausibly located in S1 [137] and latencies that decrease with age [139] (Fig. 5b). The maturation of SEPs/SEFs latencies have not been related to the myelination of corresponding white matter bundles so far. A single study in adults has related the inter-individual variability in SEPs-N140 inter-hemispheric transfer time to the axon diameter for callosal fibers supposed to connect S1 cortices [130], but these observations raised similar criticisms than for the visual system [131].

The auditory modality is already functional in utero, but its development is more protracted throughout infancy and toddlerhood than the visual and somatosensory modalities. In adults, auditory evoked potentials (AEPs) recorded at the vertex display several peaks: P1 (∼50–100 ms), N1 (∼100–150 ms), P2 (∼175–200 ms), N2 (∼270 ms), etc. [140, 141]. During the first months of life, the auditory response expands over almost all scalp electrodes and has large amplitude, making the comparison with adult components difficult. Two positive peaks are recorded over frontal electrodes, simultaneously to negative peaks over temporal and occipital electrodes at around 110–270 ms and 250–410 ms respectively. These latencies decrease exponentially with age [140, 142] (Fig. 5c). Regarding auditory evoked fields, a decrease in the latency of P2 m is observed during infancy and toddlerhood [143]. This decrease has been related during childhood to the maturation of acoustic radiations based on DTI anisotropy [144]. Nevertheless, this correlation was not significant when taking into account the dependance of both measures on age. This casted doubt on the structure-function specificity reported in this study.

3. Perspectives

Reliable correlations between the latency or conduction velocity of evoked responses and the myelination of underlying white matter bundles remain scarce so far. Further studies in the somatosensory and auditory modalities are required to confirm the results obtained in the visual modality. The comparison of anatomo-functional changes across modalities within the same subjects are needed to characterize the asynchronous development of brain networks and explore their sensitivity to distinct critical periods. For instance, the acceleration of early M/EEG components seems more prolonged in the auditory modality than in the visual modality (Fig. 5d). It might be related to the different maturational tempos of the optic and acoustic radiations. If it is the case, M/EEG/DTI combination will provide in vivo assessments of the bundles functional integrity, a useful information in a lot of pathologies (premature births, autism, demyelinating diseases to cite a few). Aside from the latency of evoked responses, several other functional parameters might be compared to the brain anatomical changes throughout development. The evolution in morphology and amplitude of peaks, notably during the preterm period, might be related to microstructural changes in the cortex, and to the cortical expansion going with the folding process. Complexity measures of EEG signals (e.g. multiscale entropy which strongly increases throughout development [145]) might also be compared more robustly across functional modalities than evoked responses. While these measures in adults do not differ between the visual and auditory modalities, auditory responses show lower complexity than visual ones in children below 5 years of age [145]. Such discrepancies between networks need to be modelled to understand how integration is realized during development for example in the case of audio-visual integration during language acquisition. We advocate for a more systematic combination of DTI and M/EEG measures in the same children to explore these questions and open new understanding of human brain maturation and of its pathologies.