Oscillator decomposition of infant fNIRS data.

The functional near-infrared spectroscopy (fNIRS) can detect hemodynamic responses in the brain and the data consist of bivariate time series of oxygenated hemoglobin (oxy-Hb) and deoxygenated hemoglobin (deoxy-Hb) on each channel. In this study, we investigate oscillatory changes in infant fNIRS signals by using the oscillator decompisition method (OSC-DECOMP), which is a statistical method for extracting oscillators from time series data based on Gaussian linear state space models. OSC-DECOMP provides a natural decomposition of fNIRS data into oscillation components in a data-driven manner and does not require the arbitrary selection of band-pass filters. We analyzed 18-ch fNIRS data (3 minutes) acquired from 21 sleeping 3-month-old infants. Five to seven oscillators were extracted on most channels, and their frequency distribution had three peaks in the vicinity of 0.01-0.1 Hz, 1.6-2.4 Hz and 3.6-4.4 Hz. The first peak was considered to reflect hemodynamic changes in response to the brain activity, and the phase difference between oxy-Hb and deoxy-Hb for the associated oscillators was at approximately 230 degrees. The second peak was attributed to cardiac pulse waves and mirroring noise. Although these oscillators have close frequencies, OSC-DECOMP can separate them through estimating their different projection patterns on oxy-Hb and deoxy-Hb. The third peak was regarded as the harmonic of the second peak. By comparing the Akaike Information Criterion (AIC) of two state space models, we determined that the time series of oxy-Hb and deoxy-Hb on each channel originate from common oscillatory activity. We also utilized the result of OSC-DECOMP to investigate the frequency-specific functional connectivity. Whereas the brain oscillator exhibited functional connectivity, the pulse waves and mirroring noise oscillators showed spatially homogeneous and independent changes. OSC-DECOMP is a promising tool for data-driven extraction of oscillation components from biological time series data.

Developmental changes in cortical sensory processing during wakefulness and sleep.

Infants are exposed to auditory and visual information during sleep as well as wakefulness. Little is known, however, about the differences in cortical processing of sensory input between these different behavioral states. In the present study, cortical hemodynamic responses to auditory and visual stimuli during wakefulness and sleep were measured in infants aged 2-10 months using functional near-infrared spectroscopy (fNIRS). While asynchronously presented auditory and visual stimuli during wakefulness induced focal responses in the corresponding sensory regions of the occipital and temporal cortices, the responses to the same stimuli during sleep were dramatically different. Auditory stimuli during sleep induced global responses over the frontal, temporal, and occipital regions, and the response pattern did not change between 2 and 10 months of age. In contrast, visual stimuli during sleep induced responses in the occipital cortex, and the response pattern exhibited developmental changes from a pattern of activation to one of deactivation around a half year of age. The functional connectivity among the cortical regions was generally higher during sleep than during wakefulness. The hemoglobin phase of oxygenation and deoxygenation (hPod) and the phase locking index of hPod (hPodL) showed general developmental changes and behavioral state dependent differences but no significant differences were seen between the stimulus types. The results suggest that the behavioral states have a fundamental impact on cortical sensory processing; (1) sensory processing during wakefulness is performed in more localized regions, (2) auditory processing is active during both wakefulness and sleep, (3) visual processing undergoes development of inhibitory mechanisms during sleep, and (4) these phenomena primarily reflect neural development rather than vascular and metabolic development.

Detection of resting state functional connectivity using partial correlation analysis: A study using multi-distance and whole-head probe near-infrared spectroscopy.

Multichannel near-infrared spectroscopy (NIRS) is a functional neuroimaging modality that enables easy-to-use and noninvasive measurement of changes in blood oxygenation levels. We developed a clinically-applicable method for estimating resting state functional connectivity (RSFC) with NIRS using a partial correlation analysis to reduce the influence of extraneural components. Using a multi-distance probe arrangement NIRS, we measured resting state brain activity for 8min in 17 healthy participants. Independent component analysis was used to extract shallow and deep signals from the original NIRS data. Pearson's correlation calculated from original signals was significantly higher than that calculated from deep signals, while partial correlation calculated from original signals was comparable to that calculated from deep (cerebral-tissue) signals alone. To further test the validity of our method, we also measured 8min of resting state brain activity using a whole-head NIRS arrangement consisting of 17 cortical regions in 80 healthy participants. Significant RSFC between neighboring, interhemispheric homologous, and some distant ipsilateral brain region pairs was revealed. Additionally, females exhibited higher RSFC between interhemispheric occipital region-pairs, in addition to higher connectivity between some ipsilateral pairs in the left hemisphere, when compared to males. The combined results of the two component experiments indicate that partial correlation analysis is effective in reducing the influence of extracerebral signals, and that NIRS is able to detect well-described resting state networks and sex-related differences in RSFC.

Development of a neurofeedback protocol targeting the frontal pole using near-infrared spectroscopy.

AIM: Neurofeedback has been studied with the aim of controlling cerebral activity. Near-infrared spectroscopy is a non-invasive neuroimaging technique used for measuring hemoglobin concentration changes in cortical surface areas with high temporal resolution. Thus, near-infrared spectroscopy may be useful for neurofeedback, which requires real-time feedback of repeated brain activation measurements. However, no study has specifically targeted neurofeedback, using near-infrared spectroscopy, in the frontal pole cortex. METHODS: We developed an original near-infrared spectroscopy neurofeedback system targeting the frontal pole cortex. Over a single day of testing, each healthy participant (n = 24) received either correct or incorrect (Sham) feedback from near-infrared spectroscopy signals, based on a crossover design. RESULTS: Under correct feedback conditions, significant activation was observed in the frontal pole cortex (P = 0.000073). Additionally, self-evaluation of control and metacognitive beliefs were associated with near-infrared spectroscopy signals (P = 0.006). CONCLUSION: The neurofeedback system developed in this study might be useful for developing control of frontal pole cortex activation.

A brain of two halves: insights into interhemispheric organization provided by near-infrared spectroscopy.

The discovery of functional lateralization and localization of the brain marked the beginning of a new era in neuroscience. While the past 150 years of research have provided a great deal of knowledge of hemispheric differences and functional relationships, the precise organization of functional laterality remains a topic of intense debate. Here I will shed light on the functional organization of the two hemispheres by reviewing some of the most recent functional near-infrared spectroscopy (NIRS) studies that have reported hemispheric differences in activation patterns. Most NIRS studies using visual stimuli, which revealed functional differentiation between the hemispheres, have reported unilateral activation, i.e., significant levels of activation in only one hemisphere. Auditory stimuli, including speech sounds, elicited bilateral activation, while the limited number of studies on young infants revealed primarily unilateral activation. The stimulus modality and the age of the participants therefore determine whether the resulting cortical activation is unilateral or bilateral. By combining a review of the existing literature with NIRS results regarding homologous connectivity across hemispheres, I hypothesized that the origin of functional lateralization changes from the independence of each hemispheric region, to mutual inhibition between homologous regions during development. Future studies applying multi-modal measurements along with NIRS and spatiotemporal analyses will further deepen our understanding of the interhemispheric organization of brain function.

Effect of auditory input on activations in infant diverse cortical regions during audiovisual processing.

A fundamental question with regard to perceptual development is how multisensory information is processed in the brain during the early stages of development. Although a growing body of evidence has shown the early emergence of modality-specific functional differentiation of the cortical regions, the interplay between sensory inputs from different modalities in the developing brain is not well understood. To study the effects of auditory input during audio-visual processing in 3-month-old infants, we evaluated the spatiotemporal cortical hemodynamic responses of 50 infants while they perceived visual objects with or without accompanying sounds. The responses were measured using 94-channel near-infrared spectroscopy over the occipital, temporal, and frontal cortices. The effects of sound manipulation were pervasive throughout the diverse cortical regions and were specific to each cortical region. Visual stimuli co-occurring with sound induced the early-onset activation of the early auditory region, followed by activation of the other regions. Removal of the sound stimulus resulted in focal deactivation in the auditory regions and reduced activation in the early visual region, the association region of the temporal and parietal cortices, and the anterior prefrontal regions, suggesting multisensory interplay. In contrast, equivalent activations were observed in the lateral occipital and lateral prefrontal regions, regardless of sound manipulation. Our findings indicate that auditory input did not generally enhance overall activation in relation to visual perception, but rather induced specific changes in each cortical region. The present study implies that 3-month-old infants may perceive audio-visual multisensory inputs by using the global network of functionally differentiated cortical regions.

A NIRS-fMRI study of resting state network.

Resting state functional connectivity, which is defined as temporal correlation of spontaneous activity between diverse brain regions, has been reported to form resting state networks (RSNs), consisting of a specific set of brain regions, based on functional magnetic resonance imaging (fMRI). Recently, studies using near-infrared spectroscopy (NIRS) reported that NIRS signals also show temporal correlation between different brain regions. The local relationship between NIRS and fMRI signals has been examined by simultaneously recording these signals when participants perform tasks or respond to stimuli. However, the NIRS-fMRI signal relationship during the resting state has been reported only between NIRS signals obtained within limited regions and whole brain fMRI signals. Therefore, it remains unclear whether NIRS signals obtained at diverse regions correlate with regional fMRI signals close to the NIRS measurement channels, especially in relation to the RSNs. In this study, we tested whether the signals measured by these different modalities during the resting state have the consistent characteristics of the RSNs. Specifically, NIRS signals during the resting state were acquired over the frontal, temporal, and occipital cortices while whole brain fMRI data was simultaneously recorded. First, by projecting the NIRS channel positions over the cerebral cortical surface, we identified the most likely anatomical locations of all NIRS channels used in the study. Next, to investigate the regional signal relationship between NIRS and fMRI, we calculated the cross-correlation between NIRS signals and fMRI signals in the brain regions adjacent to each NIRS channel. For each NIRS channel, we observed the local maxima of correlation coefficients between NIRS and fMRI signals within a radius of 2 voxels from the projection point. Furthermore, we also found that highly correlated voxels with the NIRS signal were mainly localized within brain tissues for all NIRS channels, with the exception of 2 frontal channels. Finally, by calculating the correlation between NIRS signals at a channel and whole brain fMRI signals, we observed that NIRS signals correlate with fMRI signals not only within brain regions adjacent to NIRS channels but also within distant brain regions constituting RSNs, such as the dorsal attention, fronto-parietal control, and default mode networks. These results support the idea that NIRS signals obtained at several cortical regions during the resting state mainly reflect regional spontaneous hemodynamic fluctuations that originate from spontaneous cortical activity, and include information that characterizes the RSNs. Because NIRS is relatively easy to use and a less physically demanding neuroimaging technique, our findings should facilitate a broad application of this technique to examine RSNs, especially for clinical populations and conditions unsuitable for fMRI.

Functional development in the infant brain for auditory pitch processing.

Understanding how the developing brain processes auditory information is a critical step toward the clarification of infants' perception of speech and music. We have reported that the infant brain perceives pitch information in speech sounds. Here, we used multichannel near-infrared spectroscopy to examine whether the infant brain is sensitive to information of pitch changes in auditory sequences. Three types of auditory sequences with distinct temporal structures of pitch changes were presented to 3- and 6-month-old infants: a long condition of 12 successive tones constructing a chromatic scale (600 ms), a short condition of four successive tones constructing a chromatic scale (200 ms), and a random condition of random tone sequences (50 ms per tone). The difference among the conditions was only in the sequential order of the tones, which causes pitch changes between the successive tones. We found that the bilateral temporal regions of both ages of infants showed significant activation under the three conditions. The stimulus-dependent activation was observed in the right temporoparietal region of the both infant groups; the 3- and 6-month-old infants showed the most prominent activation under the random and short conditions, respectively. Our findings indicate that the infant brain, which shows functional differentiation and lateralization in auditory-related areas, is capable of responding to more than single tones of pitch information. These results suggest that the right temporoparietal region of the infants increases sensitivity to auditory sequences, which have temporal structures similar to those of syllables in speech sounds, in the course of development.

Activation and deactivation in response to visual stimulation in the occipital cortex of 6-month-old human infants.

In an infant's developing cortex, the explanation for the mechanisms underlying the activations and deactivations in response to visual stimuli remains controversial. While previous near-infrared spectroscopy (NIRS) studies in awake infants have demonstrated cortical activations in response to meaningful/attractive visual stimuli, functional magnetic resonance imaging (fMRI) studies performed on sleeping infants showed negative blood oxygenation level-dependent (BOLD) responses to high-luminance unpatterned stimulations, such as a photic stimulation. To examine the effect of the characteristics of visual stimuli on cortical processing in awake infants, we measured cortical hemodynamic responses in 6-month-old infants during the presentation of a high-luminance unpatterned stimulus by using a NIRS system with 94 measurement channels. Results from 35 infants showed dissociated cortical responses between the occipital region and the other parts of the cortex, including the temporal and prefrontal regions. Although the visual stimulus produced sustained increases in oxygenated hemoglobin (oxy-Hb) signals in the temporal and prefrontal regions, it produced a transient increase in oxy-Hb signals followed by a salient decrease in oxy-Hb signals during a trial in a focal region of the occipital visual region. This suggests that the deactivation of the occipital visual region in response to visual stimulation is not a phenomenon that occurs only in the sleeping state, but that a high-luminance unpatterned stimulus can induce deactivation even in the awake infants.

Individual variability in the nonlinear development of the corpus callosum during infancy and toddlerhood: a longitudinal MRI analysis.

The human brain spends several years bootstrapping itself through intrinsic and extrinsic modulation, thus gradually developing both spatial organization and functions. Based on previous studies on developmental patterns and inter-individual variability of the corpus callosum (CC), we hypothesized that inherent variations of CC shape among infants emerge, depending on the position within the CC, along the developmental timeline. Here we used longitudinal magnetic resonance imaging data from infancy to toddlerhood and investigated the area, thickness, and shape of the midsagittal plane of the CC by applying multilevel modeling. The shape characteristics were extracted using the Procrustes method. We found nonlinearity, region-dependency, and inter-individual variability, as well as intra-individual consistencies, in CC development. Overall, the growth rate is faster in the first year than in the second year, and the trajectory differs between infants; the direction of CC formation in individual infants was determined within six months and maintained to two years. The anterior and posterior subregions increase in area and thickness faster than other subregions. Moreover, we clarified that the growth rate of the middle part of the CC is faster in the second year than in the first year in some individuals. Since the division of regions exhibiting different tendencies coincides with previously reported divisions based on the diameter of axons that make up the region, our results suggest that subregion-dependent individual variability occurs due to the increase in the diameter of the axon caliber, myelination partly due to experience and axon elimination during the early developmental period.

Spontaneous variability predicts compensative motor response in vocal pitch control.

Our motor system uses sensory feedback to keep desired performance. From this view, motor fluctuation is not simply 'noise' inevitably caused in the nervous system but would play a role in generating variations to explore better outcomes via sensory feedback. Vocalization system offers a good model for studying such sensory-motor interactions since we regulate vocalization by hearing our own voice. This behavior is typically observed as compensatory responses in vocalized pitch, or fundamental frequency (fo), when artificial fo shifts were induced in the auditory feedback. However, the relationship between adaptive regulation and motor exploration in vocalization has remained unclear. Here we investigated behavioral variability in spontaneous vocal fo and compensatory responses against fo shifts in the feedback, and demonstrated that larger spontaneous fluctuation correlates with greater compensation in vocal fo. This correlation was found in slow components (≤ 5 Hz) of the spontaneous fluctuation but not in fast components (between 6 and 30 Hz), and the slow one was amplified during the compensatory responses. Furthermore, the compensatory ratio was reduced when large fo shifts were applied to the auditory feedback, as if reflecting the range of motor exploration. All these findings consistently suggest the functional role of motor variability in the exploration of better vocal outcomes.

Development of global cortical networks in early infancy.

Human cognition and behaviors are subserved by global networks of neural mechanisms. Although the organization of the brain is a subject of interest, the process of development of global cortical networks in early infancy has not yet been clarified. In the present study, we explored developmental changes in these networks from several days to 6 months after birth by examining spontaneous fluctuations in brain activity, using multichannel near-infrared spectroscopy. We set up 94 measurement channels over the frontal, temporal, parietal, and occipital regions of the infant brain. The obtained signals showed complex time-series properties, which were characterized as 1/f fluctuations. To reveal the functional connectivity of the cortical networks, we calculated the temporal correlations of continuous signals between all the pairs of measurement channels. We found that the cortical network organization showed regional dependency and dynamic changes in the course of development. In the temporal, parietal, and occipital regions, connectivity increased between homologous regions in the two hemispheres and within hemispheres; in the frontal regions, it decreased progressively. Frontoposterior connectivity changed to a "U-shaped" pattern within 6 months: it decreases from the neonatal period to the age of 3 months and increases from the age of 3 months to the age of 6 months. We applied cluster analyses to the correlation coefficients and showed that the bilateral organization of the networks begins to emerge during the first 3 months of life. Our findings suggest that these developing networks, which form multiple clusters, are precursors of the functional cerebral architecture.

Frequency-specific functional connectivity in the brain during resting state revealed by NIRS.

Analyses of spontaneous hemodynamic fluctuations observed on functional magnetic resonance imaging (fMRI) have revealed the existence of temporal correlations in signal changes between widely separated brain regions during the resting state, termed "resting state functional connectivity." Recent studies have demonstrated that these correlations are also present in the hemodynamic signals measured by near infrared spectroscopy (NIRS). However, it is still uncertain whether frequency-specific characteristics exist in these signals. In the present study, we used multichannel NIRS to investigate the frequency dependency of functional connectivity between diverse regions in the cerebral cortex by decomposing fluctuations of oxygenated hemoglobin (oxy-Hb) and deoxygenated hemoglobin (deoxy-Hb) signals into various frequency bands. First, within a wide frequency range (0.009-0.1Hz), we observed that both oxy-Hb and deoxy-Hb signals showed functional connectivity within local regions and between contralateral hemispheric regions of the cortex. Next, by decomposing measured fluctuations into narrower frequency components, we determined that only oxy-Hb signals showed frequency-specific functional connectivity between the frontal and occipital regions, emerging in a narrow frequency range (0.04-0.1Hz). To clarify the coherency of functional connectivity, we calculated the average coherence values between selected channel pairs. This approach demonstrated that functional connectivity based on the oxy-Hb signals between homologous cortical regions of the contralateral hemisphere (homologous connectivity) showed high coherence over a wide frequency range (0.009-0.1Hz), whereas connectivity between the prefrontal and occipital regions (fronto-posterior connectivity) showed high coherence only within a specific narrow frequency range (0.04-0.1Hz). Our findings suggest that homologous connectivity may reflect synchronization of neural activation over a wide frequency range through direct neuroanatomical connections, whereas fronto-posterior connectivity as revealed by high coherence only within a specific narrow frequency range corresponding to the time scale of typical hemodynamic response to a single event may reflect synchronization of transient neural activation among distant cortical regions. The present study demonstrated that NIRS provides a powerful tool to elucidate network properties of the cortex during resting state.

Best practices for fNIRS publications.

The application of functional near-infrared spectroscopy (fNIRS) in the neurosciences has been expanding over the last 40 years. Today, it is addressing a wide range of applications within different populations and utilizes a great variety of experimental paradigms. With the rapid growth and the diversification of research methods, some inconsistencies are appearing in the way in which methods are presented, which can make the interpretation and replication of studies unnecessarily challenging. The Society for Functional Near-Infrared Spectroscopy has thus been motivated to organize a representative (but not exhaustive) group of leaders in the field to build a consensus on the best practices for describing the methods utilized in fNIRS studies. Our paper has been designed to provide guidelines to help enhance the reliability, repeatability, and traceability of reported fNIRS studies and encourage best practices throughout the community. A checklist is provided to guide authors in the preparation of their manuscripts and to assist reviewers when evaluating fNIRS papers.

Errata: Best practices for fNIRS publications.

[This corrects the article DOI: 10.1117/1.NPh.8.1.012101.].

General to specific development of functional activation in the cerebral cortexes of 2- to 3-month-old infants.

A critical issue in the functional development of the cerebral cortex is whether cortical regions are functionally differentiated in early infancy. Although a growing number of neuroimaging studies have revealed that functional differentiation between early sensory and association regions of the cortex is already present at 3 months of age, it is unclear how functional regions per se emerge in the earlier developmental period. Here, we present 3 possible hypotheses regarding the functional development of the cerebral cortex as follows: (1) functionally differentiated regions are prespecified in the early developmental period; (2) functional activations appear in a hierarchical order from early sensory regions to the association regions; and (3) functional activation patterns change in a general-to-specific manner, thereby increasing the localization of regions activated by a particular stimulus and increasing the exclusivity of the response to specific stimuli within a particular cortical region. In the present study, we used multichannel near-infrared spectroscopy (NIRS) to measure cortical hemodynamic responses to 2 different video images of colorful mobile objects and black-and-white checkerboard pattern reversals over the occipital and prefrontal regions in awake 2-month-old infants. Both visual stimuli produced comparative activations over broad regions of the cortex including the early sensory and association regions, supporting the general-to-specific development (Hypothesis 3). This result suggests that functional cortical regions emerge between 2 and 3 months of age for visual perception.

Prefrontal cortical involvement in young infants' analysis of novelty.

Our knowledge of infant perception and cognition is primarily based on habituation and dishabituation, but the underlying neural mechanisms for these processes per se remain unclear. It has been argued that habituation is related to building internal representations of repeated stimuli in the central nervous system, whereas dishabituation is related to an increased attention to novel items and events. This leads to a hypothesis that a distributed network including the sensory, association and prefrontal cortical regions of young infants is involved in those processes, in contrast with the classical developmental view that onset of the functioning of the prefrontal cortex is delayed. Here we examined the time evolution of spatio-temporal hemodynamic responses related to the auditory habituation and dishabituation in the temporal and prefrontal regions of 3-month-old infants by using multichannel near-infrared spectroscopy. We found that the temporal regions remained activated by repetitive auditory stimuli; however, the prefrontal regions exhibited phasic activation in relation to novel stimuli. The dissociated activation pattern between the temporal and prefrontal regions suggests that distinct cortical regions play distinct functional roles in auditory habituation and dishabituation, and that the prefrontal cortex is involved in perceiving invariance or novelty of the immediate environment in early infancy.

Unconscious and Distinctive Control of Vocal Pitch and Timbre During Altered Auditory Feedback.

Vocal control plays a critical role in smooth social communication. Speakers constantly monitor auditory feedback (AF) and make adjustments when their voices deviate from their intentions. Previous studies have shown that when certain acoustic features of the AF are artificially altered, speakers compensate for this alteration in the opposite direction. However, little is known about how the vocal control system implements compensations for alterations of different acoustic features, and associates them with subjective consciousness. The present study investigated whether compensations for the fundamental frequency (F0), which corresponds to perceived pitch, and formants, which contribute to perceived timbre, can be performed unconsciously and independently. Forty native Japanese speakers received two types of altered AF during vowel production that involved shifts of either only the formant frequencies (formant modification; Fm) or both the pitch and formant frequencies (pitch + formant modification; PFm). For each type, three levels of shift (slight, medium, and severe) in both directions (increase or decrease) were used. After the experiment, participants were tested for whether they had perceived a change in the F0 and/or formants. The results showed that (i) only formants were compensated for in the Fm condition, while both the F0 and formants were compensated for in the PFm condition; (ii) the F0 compensation exhibited greater precision than the formant compensation in PFm; and (iii) compensation occurred even when participants misperceived or could not explicitly perceive the alteration in AF. These findings indicate that non-experts can compensate for both formant and F0 modifications in the AF during vocal production, even when the modifications are not explicitly or correctly perceived, which provides further evidence for a dissociation between conscious perception and action in vocal control. We propose that such unconscious control of voice production may enhance rapid adaptation to changing speech environments and facilitate mutual communication.

Functional activation in diverse regions of the developing brain of human infants.

To understand the functional organization of the human cortex during the early postnatal period, we performed neuroimaging studies for visual perception in awake 3-month-old infants. Cortical hemodynamic responses to 2 different video images, moving mobile objects and black-and-white checkerboard pattern reversals, were observed using multichannel near-infrared spectroscopy (NIRS). Although a focal region of the occipital cortex was equally activated by both stimuli, the occipitotemporal region was activated only by the mobile objects. A possible explanation of the result is that the former and the latter regions are involved in the primary processing of visual stimuli and perception of objects with complex visual features, respectively. Furthermore, the prefrontal region was distinctly activated by the mobile objects. These results suggest that the early sensory region and the higher sensory/association and prefrontal regions are functionally differentiated by 3 months of age and that diverse regions of the cortex including the prefrontal region function during perception of visual events.

Spatial variation in the hemoglobin phase of oxygenation and deoxygenation in the developing cortex of infants.

Spontaneous low-frequency oscillatory changes in oxygenated hemoglobin (oxy-Hb) and deoxygenated hemoglobin (deoxy-Hb) are observed using functional near-infrared spectroscopy (fNIRS). A previous study showed that the time-averaged phase difference between oxy-Hb and deoxy-Hb changes, referred to as hemoglobin phase of oxygenation and deoxygenation (hPod), is sensitive to the development of the cortex. We examined phase-locking index of hPod, referred to as [Formula: see text], in addition to hPod, in neonates and 3- and 6-month-old infants using the 94-channel fNIRS data, which covered large lateral regions of the cortex. The results showed that (1) developmental changes in hPod exhibited spatial dependency; (2) [Formula: see text] increased between the neonate group and 3-month-old infant group over the posterior, but not anterior, regions of the cortex; and (3) the cortical regions of each age group were clustered in several domains with specific characteristics of hPod and [Formula: see text]. This study indicates that the neonatal cortex is composed of regions with specific characteristics of hPod and [Formula: see text], and drastic changes occur between the neonatal period and 3 months of age. This study suggests that hPod and [Formula: see text] are sensitive to the cortical region-specific development of the circulatory, blood flow, metabolic, and neurovascular functions in young infants.