Measuring the dynamic longitudinal cortex development in infants by reconstruction of temporally consistent cortical surfaces

nBEST: Deep-learning-based non-human primates Brain Extraction and Segmentation Toolbox across ages, sites and species

Accurate processing and analysis of non-human primate (NHP) brain magnetic resonance imaging (MRI) serves an indispensable role in understanding brain evolution, development, aging, and diseases. Despite the accumulation of diverse NHP brain MRI datasets at various developmental stages and from various imaging sites/scanners, existing computational tools designed for human MRI typically perform poor on NHP data, due to huge differences in brain sizes, morphologies, and imaging appearances across species, sites, and ages, highlighting the imperative for NHP-specialized MRI processing tools. To address this issue, in this paper, we present a robust, generic, and fully automated computational pipeline, called non-human primates Brain Extraction and Segmentation Toolbox (nBEST), whose main functionality includes brain extraction, non-cerebrum removal, and tissue segmentation. Building on cutting-edge deep learning techniques by employing lifelong learning to flexibly integrate data from diverse NHP populations and innovatively constructing 3D U-NeXt architecture, nBEST can well handle structural NHP brain MR images from multi-species, multi-site, and multi-developmental-stage (from neonates to the elderly). We extensively validated nBEST based on, to our knowledge, the largest assemblage dataset in NHP brain studies, encompassing 1,469 scans with 11 species (e.g., rhesus macaques, cynomolgus macaques, chimpanzees, marmosets, squirrel monkeys, etc.) from 23 independent datasets. Compared to alternative tools, nBEST outperforms in precision, applicability, robustness, comprehensiveness, and generalizability, greatly benefiting downstream longitudinal, cross-sectional, and cross-species quantitative analyses. We have made nBEST an open-source toolbox (https://github.com/TaoZhong11/nBEST) and we are committed to its continual refinement through lifelong learning with incoming data to greatly contribute to the research field.

Longitudinally consistent registration and parcellation of cortical surfaces using semi-supervised learning

Temporally consistent and accurate registration and parcellation of longitudinal cortical surfaces is of great importance in studying longitudinal morphological and functional changes of human brains. However, most existing methods are developed for registration or parcellation of a single cortical surface. When applying to longitudinal studies, these methods independently register/parcellate each surface from longitudinal scans, thus often generating longitudinally inconsistent and inaccurate results, especially in small or ambiguous cortical regions. Essentially, longitudinal cortical surface registration and parcellation are highly correlated tasks with inherently shared constraints on both spatial and temporal feature representations, which are unfortunately ignored in existing methods. To this end, we unprecedentedly propose a novel semi-supervised learning framework to exploit these inherent relationships from limited labeled data and extensive unlabeled data for more robust and consistent registration and parcellation of longitudinal cortical surfaces. Our method utilizes the spherical topology characteristic of cortical surfaces. It employs a spherical network to function as an encoder, which extracts high-level cortical features. Subsequently, we build two specialized decoders dedicated to the tasks of registration and parcellation, respectively. To extract more meaningful spatial features, we design a novel parcellation map similarity loss to utilize the relationship between registration and parcellation tasks, i.e., the parcellation map warped by the deformation field in registration should match the atlas parcellation map, thereby providing extra supervision for the registration task and augmented data for parcellation task by warping the atlas parcellation map to unlabeled surfaces. To enable temporally more consistent feature representation, we additionally enforce longitudinal consistency among longitudinal surfaces after registering them together using their concatenated features. Experiments on two longitudinal datasets of infants and adults have shown that our method achieves significant improvements on both registration/parcellation accuracy and longitudinal consistency compared to existing methods, especially in small and challenging cortical regions.

Severe Bronchopulmonary Dysplasia Adversely Affects Brain Growth in Preterm Infants

INTRODUCTION

Bronchopulmonary dysplasia (BPD) is associated with neurodevelopmental outcomes of preterm infants, but its effect on brain growth in preterm infants after the neonatal period is unknown. This study aimed to evaluate the effect of severe BPD on brain growth of preterm infants from term to 18 months of corrected age (CA).

METHODS

Sixty-three preterm infants (42 with severe BPD and 21 without severe BPD) who underwent magnetic resonance imaging at term equivalent age (TEA) and 18 months of CA were studied by using the Infant Brain Extraction and Analysis Toolbox (iBEAT). We measured segmented brain volumes and compared brain volume and brain growth velocity between the severe BPD group and the non-severe BPD group.

RESULTS

There was no significant difference in brain volumes at TEA between the groups. However, the brain volumes of the total brain and cerebral white matter in the severe BPD group were significantly smaller than those in the non-severe BPD group at 18 months of CA. The brain growth velocities from TEA to 18 months of CA in the total brain, cerebral cortex, and cerebral white matter in the severe BPD group were lower than those in the non-severe BPD group.

CONCLUSION

Brain growth in preterm infants with severe BPD from TEA age to 18 months of CA is less than that in preterm infants without severe BPD.

Cortical growth from infancy to adolescence in preterm and term-born children

Early life experiences can exert a significant influence on cortical and cognitive development. Very preterm birth exposes infants to several adverse environmental factors during hospital admission, which affect cortical architecture. However, the subsequent consequence of very preterm birth on cortical growth from infancy to adolescence has never been defined; despite knowledge of critical periods during childhood for establishment of cortical networks. Our aims were to: chart typical longitudinal cortical development and sex differences in cortical development from birth to adolescence in healthy term-born children; estimate differences in cortical development between children born at term and very preterm; and estimate differences in cortical development between children with normal and impaired cognition in adolescence. This longitudinal cohort study included children born at term (≥37 weeks' gestation) and very preterm (<30 weeks' gestation) with MRI scans at ages 0, 7 and 13 years (n = 66 term-born participants comprising 34 with one scan, 18 with two scans and 14 with three scans; n = 201 very preterm participants comprising 56 with one scan, 88 with two scans and 57 with three scans). Cognitive assessments were performed at age 13 years. Cortical surface reconstruction and parcellation were performed with state-of-the-art, equivalent MRI analysis pipelines for all time points, resulting in longitudinal cortical volume, surface area and thickness measurements for 62 cortical regions. Developmental trajectories for each region were modelled in term-born children, contrasted between children born at term and very preterm, and contrasted between all children with normal and impaired cognition. In typically developing term-born children, we documented anticipated patterns of rapidly increasing cortical volume, area and thickness in early childhood, followed by more subtle changes in later childhood, with smaller cortical size in females than males. In contrast, children born very preterm exhibited increasingly reduced cortical volumes, relative to term-born children, particularly during ages 0-7 years in temporal cortical regions. This reduction in cortical volume in children born very preterm was largely driven by increasingly reduced cortical thickness rather than area. This resulted in amplified cortical volume and thickness reductions by age 13 years in individuals born very preterm. Alterations in cortical thickness development were found in children with impaired language and memory. This study shows that the neurobiological impact of very preterm birth on cortical growth is amplified from infancy to adolescence. These data further inform the long-lasting impact on cortical development from very preterm birth, providing broader insights into neurodevelopmental consequences of early life experiences.

Parenting Influences on Frontal Lobe Gray Matter and Preterm Toddlers' Problem-Solving Skills

Children born preterm often face challenges with self-regulation during toddlerhood. This study examined the relationship between prematurity, supportive parent behaviors, frontal lobe gray matter volume (GMV), and emotion regulation (ER) among toddlers during a parent-assisted, increasingly complex problem-solving task, validated for this age range. Data were collected from preterm toddlers (n = 57) ages 15-30 months corrected for prematurity and their primary caregivers. MRI data were collected during toddlers' natural sleep. The sample contained three gestational groups: 22-27 weeks (extremely preterm; EPT), 28-33 weeks (very preterm; VPT), and 34-36 weeks (late preterm; LPT). Older toddlers became more compliant as the Tool Task increased in difficulty, but this pattern varied by gestational group. Engagement was highest for LPT toddlers, for older toddlers, and for the easiest task condition. Parents did not differentiate their support depending on task difficulty or their child's age or gestational group. Older children had greater frontal lobe GMV, and for EPT toddlers only, more parent support was related to larger right frontal lobe GMV. We found that parent support had the greatest impact on high birth risk (≤27 gestational weeks) toddler brain development, thus early parent interventions may normalize preterm child neurodevelopment and have lasting impacts.

A critical period plasticity framework for the sensorimotor-association axis of cortical neurodevelopment

To understand human brain development it is necessary to describe not only the spatiotemporal patterns of neurodevelopment but also the neurobiological mechanisms that underlie them. Human neuroimaging studies have provided evidence for a hierarchical sensorimotor-to-association (S-A) axis of cortical neurodevelopment. Understanding the biological mechanisms that underlie this program of development using traditional neuroimaging approaches has been challenging. Animal models have been used to identify periods of enhanced experience-dependent plasticity - 'critical periods' - that progress along cortical hierarchies and are governed by a conserved set of neurobiological mechanisms that promote and then restrict plasticity. In this review we hypothesize that the S-A axis of cortical development in humans is partly driven by the cascading maturation of critical period plasticity mechanisms. We then describe how recent advances in in vivo neuroimaging approaches provide a promising path toward testing this hypothesis by linking signals derived from non-invasive imaging to critical period mechanisms.

Mapping Genetic Topography of Cortical Thickness and Surface Area in Neonatal Brains

Adult twin neuroimaging studies have revealed that cortical thickness (CT) and surface area (SA) are differentially influenced by genetic information, leading to their spatially distinct genetic patterning and topography. However, the postnatal origins of the genetic topography of CT and SA remain unclear, given the dramatic cortical development from neonates to adults. To fill this critical gap, this study unprecedentedly explored how genetic information differentially regulates the spatial topography of CT and SA in the neonatal brain by leveraging brain magnetic resonance (MR) images from 202 twin neonates with minimal influence by the complicated postnatal environmental factors. We capitalized on infant-dedicated computational tools and a data-driven spectral clustering method to parcellate the cerebral cortex into a set of distinct regions purely according to the genetic correlation of cortical vertices in terms of CT and SA, respectively, and accordingly created the first genetically informed cortical parcellation maps of neonatal brains. Both genetic parcellation maps exhibit bilaterally symmetric and hierarchical patterns, but distinct spatial layouts. For CT, regions with closer genetic relationships demonstrate an anterior-posterior (A-P) division, while for SA, regions with greater genetic proximity are typically within the same lobe. Certain genetically informed regions exhibit strong similarities between neonates and adults, with the most striking similarities in the medial surface in terms of SA, despite their overall substantial differences in genetic parcellation maps. These results greatly advance our understanding of the development of genetic influences on the spatial patterning of cortical morphology.SIGNIFICANCE STATEMENT Genetic influences on cortical thickness (CT) and surface area (SA) are complex and could evolve throughout the lifespan. However, studies revealing distinct genetic topography of CT and SA have been limited to adults. Using brain structural magnetic resonance (MR) images of twins, we unprecedentedly discovered the distinct genetically-informed parcellation maps of CT and SA in neonatal brains, respectively. Each genetic parcellation map comprises a distinct spatial layout of cortical regions, where vertices within the same region share high genetic correlation. These genetic parcellation maps of CT and SA of neonates largely differ from those of adults, despite their highly remarkable similarities in the medial cortex of SA. These discoveries provide important insights into the genetic organization of the early cerebral cortex development.

Fine-grained functional parcellation maps of the infant cerebral cortex

Resting-state functional MRI (rs-fMRI) is widely used to examine the dynamic brain functional development of infants, but these studies typically require precise cortical parcellation maps, which cannot be directly borrowed from adult-based functional parcellation maps due to the substantial differences in functional brain organization between infants and adults. Creating infant-specific cortical parcellation maps is thus highly desired but remains challenging due to difficulties in acquiring and processing infant brain MRIs. In this study, we leveraged 1064 high-resolution longitudinal rs-fMRIs from 197 typically developing infants and toddlers from birth to 24 months who participated in the Baby Connectome Project to develop the first set of infant-specific, fine-grained, surface-based cortical functional parcellation maps. To establish meaningful cortical functional correspondence across individuals, we performed cortical co-registration using both the cortical folding geometric features and the local gradient of functional connectivity (FC). Then we generated both age-related and age-independent cortical parcellation maps with over 800 fine-grained parcels during infancy based on aligned and averaged local gradient maps of FC across individuals. These parcellation maps reveal complex functional developmental patterns, such as changes in local gradient, network size, and local efficiency, especially during the first 9 postnatal months. Our generated fine-grained infant cortical functional parcellation maps are publicly available at https://www.nitrc.org/projects/infantsurfatlas/ for advancing the pediatric neuroimaging field.

Four-dimensional mapping of dynamic longitudinal brain subcortical development and early learning functions in infants

Brain subcortical structures are paramount in many cognitive functions and their aberrations during infancy are predisposed to various neurodevelopmental and neuropsychiatric disorders, making it highly essential to characterize the early subcortical normative growth patterns. This study investigates the volumetric development and surface area expansion of six subcortical structures and their associations with Mullen scales of early learning by leveraging 513 high-resolution longitudinal MRI scans within the first two postnatal years. Results show that (1) each subcortical structure (except for the amygdala with an approximately linear increase) undergoes rapid nonlinear volumetric growth after birth, which slows down at a structure-specific age with bilaterally similar developmental patterns; (2) Subcortical local area expansion reveals structure-specific and spatiotemporally heterogeneous patterns; (3) Positive associations between thalamus and both receptive and expressive languages and between caudate and putamen and fine motor are revealed. This study advances our understanding of the dynamic early subcortical developmental patterns.

Modeling individual differences in the timing of change onset and offset

Individual differences in the timing of developmental processes are often of interest in longitudinal studies, yet common statistical approaches to modeling change cannot directly estimate the timing of when change occurs. The time-to-criterion framework was recently developed to incorporate the timing of a prespecified criterion value; however, this framework has difficulty accommodating contexts where the criterion value differs across people or when the criterion value is not known a priori, such as when the interest is in individual differences in when change starts or stops. This article combines aspects of reparameterized quadratic models and multiphase models to provide information on the timing of change. We first consider the more common situation of modeling decelerating change to an offset point, defined as the point in time at which change ceases. For increasing trajectories, the offset occurs when the criterion attains its maximum ("inverted J-shaped" trajectories). For decreasing trajectories, offset instead occurs at the minimum. Our model allows for individual differences in both the timing of offset and ultimate level of the outcome. The same model, reparameterized slightly, captures accelerating change from a point of onset ("J-shaped" trajectories). We then extend the framework to accommodate "S-shaped" curves where both the onset and offset of change are within the observation window. We provide demonstrations that span neuroscience, educational psychology, developmental psychology, and cognitive science, illustrating the applicability of the modeling framework to a variety of research questions about individual differences in the timing of change. (PsycInfo Database Record (c) 2023 APA, all rights reserved).

Altered topological organization of resting-state functional networks in children with infantile spasms

Covering neuroimaging evidence has demonstrated that epileptic symptoms are associated with the disrupted topological architecture of the brain network. Infantile spasms (IS) as an age-specific epileptic encephalopathy also showed abnormal structural or functional connectivity in specific brain regions or specific networks. However, little is known about the topological alterations of whole-brain functional networks in patients with IS. To fill this gap, we used the graph theoretical analysis to investigate the topological properties (whole-brain small-world property and modular interaction) in 17 patients with IS and 34 age- and gender-matched healthy controls. The functional networks in both groups showed efficient small-world architecture over the sparsity range from 0.05 to 0.4. While patients with IS showed abnormal global properties characterized by significantly decreased normalized clustering coefficient, normalized path length, small-worldness, local efficiency, and significantly increased global efficiency, implying a shift toward a randomized network. Modular analysis revealed decreased intra-modular connectivity within the default mode network (DMN) and fronto-parietal network but increased inter-modular connectivity between the cingulo-opercular network and occipital network. Moreover, the decreased intra-modular connectivity in DMN was significantly negatively correlated with seizure frequency. The inter-modular connectivity between the cingulo-opercular and occipital network also showed a significant correlation with epilepsy frequency. Together, the current study revealed the disrupted topological organization of the whole-brain functional network, which greatly advances our understanding of neuronal architecture in IS and may contribute to predict the prognosis of IS as disease biomarkers.

Mother’s physical activity during pregnancy and newborn’s brain cortical development

Background

Physical activity is known to improve mental health, and is regarded as safe and desirable for uncomplicated pregnancy. In this novel study, we aim to evaluate whether there are associations between maternal physical activity during pregnancy and neonatal brain cortical development.

Methods

Forty-four mother/newborn dyads were included in this longitudinal study. Healthy pregnant women were recruited and their physical activity throughout pregnancy were documented using accelerometers worn for 3–7 days for each of the 6 time points at 4–10, ∼12, ∼18, ∼24, ∼30, and ∼36 weeks of pregnancy. Average daily total steps and daily total activity count as well as daily minutes spent in sedentary/light/moderate/vigorous activity modes were extracted from the accelerometers for each time point. At ∼2 weeks of postnatal age, their newborns underwent an MRI examination of the brain without sedation, and 3D T1-weighted brain structural images were post-processed by the iBEAT2.0 software utilizing advanced deep learning approaches. Cortical surface maps were reconstructed from the segmented brain images and parcellated to 34 regions in each brain hemisphere, and mean cortical thickness for each region was computed for partial correlation analyses with physical activity measures, with appropriate multiple comparison corrections and potential confounders controlled.

Results

At 4–10 weeks of pregnancy, mother’s daily total activity count positively correlated (FDR corrected P ≤ 0.05) with newborn’s cortical thickness in the left caudal middle frontal gyrus (rho = 0.48, P = 0.04), right medial orbital frontal gyrus (rho = 0.48, P = 0.04), and right transverse temporal gyrus (rho = 0.48, P = 0.04); mother’s daily time in moderate activity mode positively correlated with newborn’s cortical thickness in the right transverse temporal gyrus (rho = 0.53, P = 0.03). At ∼24 weeks of pregnancy, mother’s daily total activity count positively correlated (FDR corrected P ≤ 0.05) with newborn’s cortical thickness in the left (rho = 0.56, P = 0.02) and right isthmus cingulate gyrus (rho = 0.50, P = 0.05).

Conclusion

We identified significant relationships between physical activity in healthy pregnant women during the 1st and 2nd trimester and brain cortical development in newborns. Higher maternal physical activity level is associated with greater neonatal brain cortical thickness, presumably indicating better cortical development.

Spherical Transformer on Cortical Surfaces

Motivated by the recent great success of attention modeling in computer vision, it is highly desired to extend the Transformer architecture from the conventional Euclidean space to non-Euclidean spaces. Given the intrinsic spherical topology of brain cortical surfaces in neuroimaging, in this study, we propose a novel Spherical Transformer, an effective general-purpose backbone using the self-attention mechanism for analysis of cortical surface data represented by triangular meshes. By mapping the cortical surface onto a sphere and splitting it uniformly into overlapping spherical surface patches, we encode the long-range dependency within each patch by the self-attention operation and formulate the cross-patch feature transmission via overlapping regions. By limiting the self-attention computation to local patches, our proposed Spherical Transformer preserves detailed contextual information and enjoys great efficiency with linear computational complexity with respect to the patch size. Moreover, to better process longitudinal cortical surfaces, which are increasingly popular in neuroimaging studies, we unprecedentedly propose the spatiotemporal self-attention operation to jointly extract the spatial context and dynamic developmental patterns within a single layer, thus further enlarging the expressive power of the generated representation. To comprehensively evaluate the performance of our Spherical Transformer, we validate it on a surface-level prediction task and a vertex-level dense prediction task, respectively, i.e., the cognition prediction and cortical thickness map development prediction, which are important in early brain development mapping. Both applications demonstrate the competitive performance of our Spherical Transformer in comparison with the state-of-the-art methods.

Fast Spherical Mapping of Cortical Surface Meshes Using Deep Unsupervised Learning

Spherical mapping of cortical surface meshes provides a more convenient and accurate space for cortical surface registration and analysis and thus has been widely adopted in neuroimaging field. Conventional approaches typically first inflate and project the original cortical surface mesh onto a sphere to generate an initial spherical mesh which contains large distortions. Then they iteratively reshape the spherical mesh to minimize the metric (distance), area or angle distortions. However, these methods suffer from two major issues: 1) the iterative optimization process is computationally expensive, making them not suitable for large-scale data processing; 2) when metric distortion cannot be further minimized, either area or angle distortion is minimized at the expense of the other, which is not flexible to generate application-specific meshes based on both of them. To address these issues, for the first time, we propose a deep learning-based algorithm to learn the mapping between the original cortical surface and spherical surface meshes. Specifically, we take advantage of the Spherical U-Net model to learn the spherical diffeomorphic deformation field for minimizing the distortions between the icosahedron-reparameterized original surface and spherical surface meshes. The end-to-end unsupervised learning scheme is very flexible to incorporate various optimization objectives. We further integrate it into a coarse-to-fine multi-resolution framework for better correcting fine-scaled distortions. We have validated our method on 800+ cortical surfaces, demonstrating reduced distortions than FreeSurfer (the most popularly used tool), while speeding up the process from 20 min to 5 s.

Harmonization of Multi-site Cortical Data Across the Human Lifespan

Neuroimaging data harmonization has become a prerequisite in integrative data analytics for standardizing a wide variety of data collected from multiple studies and enabling interdisciplinary research. The lack of standardized image acquisition and computational procedures introduces non-biological variability and inconsistency in multi-site data, complicating downstream statistical analyses. Here, we propose a novel statistical technique to retrospectively harmonize multi-site cortical data collected longitudinally and cross-sectionally between birth and 100 years. We demonstrate that our method can effectively eliminate non-biological disparities from cortical thickness and myelination measurements, while preserving biological variation across the entire lifespan. Our harmonization method will foster large-scale population studies by providing comparable data required for investigating developmental and aging processes.

Influence of gonadal steroids on cortical surface area in infancy

Sex differences in the human brain emerge as early as mid-gestation and have been linked to sex hormones, particularly testosterone. Here, we analyzed the influence of markers of early sex hormone exposure (polygenic risk score (PRS) for testosterone, salivary testosterone, number of CAG repeats, digit ratios, and PRS for estradiol) on the growth pattern of cortical surface area in a longitudinal cohort of 722 infants. We found PRS for testosterone and right-hand digit ratio to be significantly associated with surface area, but only in females. PRS for testosterone at the most stringent P value threshold was positively associated with surface area development over time. Higher right-hand digit ratio, which is indicative of low prenatal testosterone levels, was negatively related to surface area in females. The current work suggests that variation in testosterone levels during both the prenatal and postnatal period may contribute to cortical surface area development in female infants.

Gyral peaks: Novel gyral landmarks in developing macaque brains

Cerebral cortex development undergoes a variety of processes, which provide valuable information for the study of the developmental mechanism of cortical folding as well as its relationship to brain structural architectures and brain functions. Despite the variability in the anatomy–function relationship on the higher‐order cortex, recent studies have succeeded in identifying typical cortical landmarks, such as sulcal pits, that bestow specific functional and cognitive patterns and remain invariant across subjects and ages with their invariance being related to a gene‐mediated proto‐map. Inspired by the success of these studies, we aim in this study at defining and identifying novel cortical landmarks, termed gyral peaks, which are the local highest foci on gyri. By analyzing data from 156 MRI scans of 32 macaque monkeys with the age spanned from 0 to 36 months, we identified 39 and 37 gyral peaks on the left and right hemispheres, respectively. Our investigation suggests that these gyral peaks are spatially consistent across individuals and relatively stable within the age range of this dataset. Moreover, compared with other gyri, gyral peaks have a thicker cortex, higher mean curvature, more pronounced hub‐like features in structural connective networks, and are closer to the borders of structural connectivity‐based cortical parcellations. The spatial distribution of gyral peaks was shown to correlate with that of other cortical landmarks, including sulcal pits. These results provide insights into the spatial arrangement and temporal development of gyral peaks as well as their relation to brain structure and function.

Quantitative Structural Brain Magnetic Resonance Imaging Analyses: Methodological Overview and Application to Rett Syndrome

Congenital genetic disorders often present with neurological manifestations such as neurodevelopmental disorders, motor developmental retardation, epilepsy, and involuntary movement. Through qualitative morphometric evaluation of neuroimaging studies, remarkable structural abnormalities, such as lissencephaly, polymicrogyria, white matter lesions, and cortical tubers, have been identified in these disorders, while no structural abnormalities were identified in clinical settings in a large population. Recent advances in data analysis programs have led to significant progress in the quantitative analysis of anatomical structural magnetic resonance imaging (MRI) and diffusion-weighted MRI tractography, and these approaches have been used to investigate psychological and congenital genetic disorders. Evaluation of morphometric brain characteristics may contribute to the identification of neuroimaging biomarkers for early diagnosis and response evaluation in patients with congenital genetic diseases. This mini-review focuses on the methodologies and attempts employed to study Rett syndrome using quantitative structural brain MRI analyses, including voxel- and surface-based morphometry and diffusion-weighted MRI tractography. The mini-review aims to deepen our understanding of how neuroimaging studies are used to examine congenital genetic disorders.

A 4D infant brain volumetric atlas based on the UNC/UMN baby connectome project (BCP) cohort

Spatiotemporal (four-dimensional) infant-dedicated brain atlases are essential for neuroimaging analysis of early dynamic brain development. However, due to the substantial technical challenges in the acquisition and processing of infant brain MR images, 4D atlases densely covering the dynamic brain development during infancy are still scarce. Few existing ones generally have fuzzy tissue contrast and low spatiotemporal resolution, leading to degraded accuracy of atlas-based normalization and subsequent analyses. To address this issue, in this paper, we construct a 4D structural MRI atlas for infant brains based on the UNC/UMN Baby Connectome Project (BCP) dataset, which features a high spatial resolution, extensive age-range coverage, and densely sampled time points. Specifically, 542 longitudinal T1w and T2w scans from 240 typically developing infants up to 26-month of age were utilized for our atlas construction. To improve the co-registration accuracy of the infant brain images, which typically exhibit dynamic appearance with low tissue contrast, we employed the state-of-the-art registration method and leveraged our generated reliable brain tissue probability maps in addition to the intensity images to improve the alignment of individual images. To achieve consistent region labeling on both infant and adult brain images for facilitating region-based analysis across ages, we mapped the widely used Desikan cortical parcellation onto our atlas by following an age-decreasing mapping manner. Meanwhile, the typical subcortical structures were manually delineated to facilitate the studies related to the subcortex. Compared with the existing infant brain atlases, our 4D atlas has much higher spatiotemporal resolution and preserves more structural details, and thus can boost accuracy in neurodevelopmental analysis during infancy.

SPHERICAL TRANSFORMER FOR QUALITY ASSESSMENT OF PEDIATRIC CORTICAL SURFACES

Brain cortical surfaces, which have an intrinsic spherical topology, are typically represented by triangular meshes and mapped onto a spherical manifold in neuroimaging analysis. Inspired by the strong capability of feature learning in Convolutional Neural Networks (CNNs), spherical CNNs have been developed accordingly and achieved many successes in cortical surface analysis. Motivated by the recent success of the transformer, in this paper, for the first of time, we extend the transformer into the spherical space and propose the spherical transformer, which can better learn contextual and structural features than spherical CNNs. We applied the spherical transformer in the important task of automatic quality assessment of infant cortical surfaces, which is a necessary procedure to identify problematic cases due to extremely low tissue contrast and strong motion effects in pediatric brain MRI studies. Experiments on 1,860 infant cortical surfaces validated its superior effectiveness and efficiency in comparison with spherical CNNs.

White matter myelination during early infancy is linked to spatial gradients and myelin content at birth

Development of myelin, a fatty sheath that insulates nerve fibers, is critical for brain function. Myelination during infancy has been studied with histology, but postmortem data cannot evaluate the longitudinal trajectory of white matter development. Here, we obtained longitudinal diffusion MRI and quantitative MRI measures of longitudinal relaxation rate (R1) of white matter in 0, 3 and 6 months-old human infants, and developed an automated method to identify white matter bundles and quantify their properties in each infant's brain. We find that R1 increases from newborns to 6-months-olds in all bundles. R1 development is nonuniform: there is faster development in white matter that is less mature in newborns, and development rate increases along inferior-to-superior as well as anterior-to-posterior spatial gradients. As R1 is linearly related to myelin fraction in white matter bundles, these findings open new avenues to elucidate typical and atypical white matter myelination in early infancy.

Existence of Functional Connectome Fingerprint during Infancy and Its Stability over Months

The functional connectome fingerprint is a cluster of individualized brain functional connectivity patterns that are capable of distinguishing one individual from others. Although its existence has been demonstrated in adolescents and adults, whether such individualized patterns exist during infancy is barely investigated despite its importance in identifying the origin of the intrinsic connectome patterns that potentially mirror distinct behavioral phenotypes. To fill this knowledge gap, capitalizing on a longitudinal high-resolution structural and resting-state functional MRI dataset with 104 human infants (53 females) with 806 longitudinal scans (age, 16-876 d) and infant-specific functional parcellation maps, we observe that the brain functional connectome fingerprint may exist since infancy and keeps stable over months during early brain development. Specifically, we achieve an ∼78% individual identification rate by using ∼5% selected functional connections, compared with the best identification rate of 60% without connection selection. The frontoparietal networks recognized as the most contributive networks in adult functional connectome fingerprinting retain their superiority in infants despite being widely acknowledged as rapidly developing systems during childhood. The existence and stability of the functional connectome fingerprint are further validated on adjacent age groups. Moreover, we show that the infant frontoparietal networks can reach similar accuracy in predicting individual early learning composite scores as the whole-brain connectome, again resembling the observations in adults and highlighting the relevance of functional connectome fingerprint to cognitive performance. For the first time, these results suggest that each individual may retain a unique and stable marker of functional connectome during early brain development.SIGNIFICANCE STATEMENT Functional connectome fingerprinting during infancy featuring rapid brain development remains almost uninvestigated even though it is essential for understanding the early individual-level intrinsic pattern of functional organization and its relationship with distinct behavioral phenotypes. With an infant-tailored functional connection selection and validation strategy, we strive to provide the delineation of the infant functional connectome fingerprint by examining its existence, stability, and relationship with early cognitive performance. We observe that the brain functional connectome fingerprint may exist since early infancy and remains stable over months during the first 2 years. The identified key contributive functional connections and networks for fingerprinting are also verified to be highly predictive for cognitive score prediction, which reveals the association between infant connectome fingerprint and cognitive performance.

Developmental abnormalities of structural covariance networks of cortical thickness and surface area in autistic infants within the first 2 years

Converging evidence supports that a collection of brain regions is functionally or anatomically abnormal in autistic subjects. Structural covariance networks (SCNs) representing patterns of coordinated regional maturation are widely used to study abnormalities associated with neurodisorders. However, the possible developmental changes of SCNs in autistic individuals during the first 2 postnatal years, which features dynamic development and can potentially serve as biomarkers, remain unexplored. To fill this gap, for the first time, SCNs of cortical thickness and surface area were constructed and investigated in infants at high familial risk for autism and typically developing infants in this study. Group differences of SCNs emerge at 12 months of age in surface area. By 24 months of age, the autism group shows significantly increased integration, decreased segregation, and decreased small-worldness, compared with controls. The SCNs of surface area are deteriorated and shifted toward randomness in autistic infants. The abnormal brain regions changed during development, and the group differences of the left lateral occipital cortex become more prominent with age. These results indicate that autism has more significant influences on coordinated development of surface area than that of cortical thickness and the occipital cortex maybe an important biomarker of autism during infancy.

Genetic Influences on Longitudinal Trajectories of Cortical Thickness and Surface Area during the First 2 Years of Life

Genetic influences on cortical thickness (CT) and surface area (SA) are known to vary across the life span. Little is known about the extent to which genetic factors influence CT and SA in infancy and toddlerhood. We performed the first longitudinal assessment of genetic influences on variation in CT and SA in 501 twins who were aged 0-2 years. We observed substantial additive genetic influences on both average CT (0.48 in neonates, 0.37 in 1-year-olds, and 0.44 in 2-year-olds) and total SA (0.59 in neonates, 0.74 in 1-year-olds, and 0.73 in 2-year-olds). In addition, we found strong heritability of the change in average CT (0.49) from neonates to 1-year-olds, but not from 1- to 2-year-olds. Moreover, we found strong genetic correlations for average CT (rG = 0.92) between 1- and 2-year-olds and strong genetic correlations for total SA across all timepoints (rG = 0.96 between neonates and 1-year-olds, rG = 1 between 1- and 2-year-olds). In addition, we found CT and SA are strongly genetic correlated at birth, but weaken over time. Overall, results suggest a dynamic genetic relationship between CT and SA during first 2 years of life and provide novel insights into how genetic influences shape the cortical structure during early brain development.

Cortical development coupling between surface area and sulcal depth on macaque brains

Postnatal development of cerebral cortex is associated with a variety of neuronal processes and is thus critical to development of brain function and cognition. Longitudinal changes of cortical morphology and topology, such as postnatal cortical thinning and flattening have been widely studied. However, thorough and systematic investigation of such cortical change, including how to quantify it from multiple spatial directions and how to relate it to surface topology, is rarely found. In this work, based on a longitudinal macaque neuroimaging dataset, we quantified local changes in gyral white matter's surface area and sulcal depth during early development. We also investigated how these two metrics are coupled and how this coupling is linked to cortical surface topology, underlying white matter, and positions of functional areas. Semi-parametric generalized additive models were adopted to quantify the longitudinal changes of surface area (A) and sulcal depth (D), and the coupling patterns between them. This resulted in four classes of regions, according to how they change compared with global change throughout early development: slower surface area change and slower sulcal depth change (slowA_slowD), slower surface area change and faster sulcal depth change (slowA_fastD), faster surface area change and slower sulcal depth change (fastA_slowD), and faster surface area change and faster sulcal depth change (fastA_fastD). We found that cortex-related metrics, including folding pattern and cortical thickness, vary along slowA_fastD-fastA_slowD axis, and structural connection-related metrics vary along fastA_fastD-slowA_slowD axis, with which brain functional sites align better. It is also found that cortical landmarks, including sulcal pits and gyral hinges, spatially reside on the borders of the four patterns. These findings shed new lights on the relationship between cortex development, surface topology, axonal wiring pattern and brain functions.

Longitudinal brain atlases of early developing cynomolgus macaques from birth to 48 months of age

Longitudinal brain imaging atlases with densely sampled time-points and ancillary anatomical information are of fundamental importance in studying early developmental characteristics of human and non-human primate brains during infancy, which feature extremely dynamic imaging appearance, brain shape and size. However, for non-human primates, which are highly valuable animal models for understanding human brains, the existing brain atlases are mainly developed based on adults or adolescents, denoting a notable lack of temporally densely-sampled atlases covering the dynamic early brain development. To fill this critical gap, in this paper, we construct a comprehensive set of longitudinal brain atlases and associated tissue probability maps (gray matter, white matter, and cerebrospinal fluid) with totally 12 time-points from birth to 4 years of age (i.e., 1, 2, 3, 4, 5, 6, 9, 12, 18, 24, 36, and 48 months of age) based on 175 longitudinal structural MRI scans from 39 typically-developing cynomolgus macaques, by leveraging state-of-the-art computational techniques tailored for early developing brains. Furthermore, to facilitate region-based analysis using our atlases, we also provide two popular hierarchy parcellations, i.e., cortical hierarchy maps (6 levels) and subcortical hierarchy maps (6 levels), on our longitudinal macaque brain atlases. These early developing atlases, which have the densest time-points during infancy (to the best of our knowledge), will greatly facilitate the studies of macaque brain development.

Infants' cortex undergoes microstructural growth coupled with myelination during development

Development of cortical tissue during infancy is critical for the emergence of typical brain functions in cortex. However, how cortical microstructure develops during infancy remains unknown. We measured the longitudinal development of cortex from birth  to six months of age  using multimodal quantitative imaging of cortical microstructure. Here we show that infants' cortex undergoes profound microstructural tissue growth during the first six months of human life. Comparison of postnatal to prenatal transcriptomic gene expression data demonstrates that myelination and synaptic processes are dominant contributors to this postnatal microstructural tissue growth. Using visual cortex as a model system, we find hierarchical microstructural growth: higher-level visual areas have less mature tissue at birth than earlier visual areas but grow at faster rates. This overturns the prominent view that visual areas that are most mature at birth develop fastest. Together, in vivo, longitudinal, and quantitative measurements, which we validated with ex vivo transcriptomic data, shed light on the rate, sequence, and biological mechanisms of developing cortical systems during early infancy. Importantly, our findings propose a hypothesis that cortical myelination is a key factor in cortical development during early infancy, which has important implications for diagnosis of neurodevelopmental disorders and delays in infants.

Maternal Obesity during Pregnancy is Associated with Lower Cortical Thickness in the Neonate Brain

BACKGROUND AND PURPOSE

Recent studies have suggested that maternal obesity during pregnancy is associated with differences in neurodevelopmental outcomes in children. In this study, we aimed to investigate the relationships between maternal obesity during pregnancy and neonatal brain cortical development.

MATERIALS AND METHODS

Forty-four healthy women (28 normal-weight, 16 obese) were prospectively recruited at <10 weeks' gestation, and their healthy full-term neonates (23 boys, 21 girls) underwent brain MR imaging. All pregnant women had their body composition (fat mass percentage) measured at ∼12 weeks of pregnancy. All neonates were scanned at ∼2 weeks of age during natural sleep without sedation, and their 3D T1-weighted images were postprocessed by the new iBEAT2.0 software. Brain MR imaging segmentation and cortical surface reconstruction and parcellation were completed using age-appropriate templates. Mean cortical thickness for 34 regions in each brain hemisphere defined by the UNC Neonatal Cortical Surface Atlas was measured, compared between groups, and correlated with maternal body fat mass percentage, controlled for neonate sex and race, postmenstrual age at MR imaging, maternal age at pregnancy, and the maternal intelligence quotient and education.

RESULTS

Neonates born to obese mothers showed significantly lower (P ≤ .05, false discovery rate-corrected) cortical thickness in the left pars opercularis gyrus, left pars triangularis gyrus, and left rostral middle frontal gyrus. Mean cortical thickness in these frontal lobe regions negatively correlated (R = -0.34, P = .04; R = -0.50, P = .001; and R = -0.42, P = .01; respectively) with the maternal body fat mass percentage measured at early pregnancy.

CONCLUSIONS

Maternal obesity during pregnancy is associated with lower neonate brain cortical thickness in several frontal lobe regions important for language and executive functions.

A Deep Network for Joint Registration and Parcellation of Cortical Surfaces

Cortical surface registration and parcellation are two essential steps in neuroimaging analysis. Conventionally, they are performed independently as two tasks, ignoring the inherent connections of these two closely-related tasks. Essentially, both tasks rely on meaningful cortical feature representations, so they can be jointly optimized by learning shared useful cortical features. To this end, we propose a deep learning framework for joint cortical surface registration and parcellation. Specifically, our approach leverages the spherical topology of cortical surfaces and uses a spherical network as the shared encoder to first learn shared features for both tasks. Then we train two task-specific decoders for registration and parcellation, respectively. We further exploit the more explicit connection between them by incorporating the novel parcellation map similarity loss to enforce the boundary consistency of regions, thereby providing extra supervision for the registration task. Conversely, parcellation network training also benefits from the registration, which provides a large amount of augmented data by warping one surface with manual parcellation map to another surface, especially when only few manually-labeled surfaces are available. Experiments on a dataset with more than 600 cortical surfaces show that our approach achieves large improvements on both parcellation and registration accuracy (over separately trained networks) and enables training high-quality parcellation and registration models using much fewer labeled data.

Learning 4D Infant Cortical Surface Atlas with Unsupervised Spherical Networks

Spatiotemporal (4D) cortical surface atlas during infancy plays an important role for surface-based visualization, normalization and analysis of the dynamic early brain development. Conventional atlas construction methods typically rely on classical group-wise registration on sub-populations and ignore longitudinal constraints, thus having three main issues: 1) constructing templates at discrete time points; 2) resulting in longitudinal inconsistency among different age's atlases; and 3) taking extremely long runtime. To address these issues, in this paper, we propose a fast unsupervised learning-based surface atlas construction framework incorporating longitudinal constraints to enforce the within-subject temporal correspondence in the atlas space. To well handle the difficulties of learning large deformations, we propose a multi-level multimodal spherical registration network to perform cortical surface registration in a coarse-to-fine manner. Thus, only small deformations need to be estimated at each resolution level using the registration network, which further improves registration accuracy and atlas quality. Our constructed 4D infant cortical surface atlas based on 625 longitudinal scans from 291 infants is temporally continuous, in contrast to the state-of-the-art UNC 4D Infant Surface Atlas, which only provides the atlases at a few discrete sparse time points. By evaluating the intra- and inter-subject spatial normalization accuracy after alignment onto the atlas, our atlas demonstrates more detailed and fine-grained cortical patterns, thus leading to higher accuracy in surface registration.

S3Reg: Superfast Spherical Surface Registration Based on Deep Learning

Cortical surface registration is an essential step and prerequisite for surface-based neuroimaging analysis. It aligns cortical surfaces across individuals and time points to establish cross-sectional and longitudinal cortical correspondences to facilitate neuroimaging studies. Though achieving good performance, available methods are either time consuming or not flexible to extend to multiple or high dimensional features. Considering the explosive availability of large-scale and multimodal brain MRI data, fast surface registration methods that can flexibly handle multimodal features are desired. In this study, we develop a Superfast Spherical Surface Registration (S3Reg) framework for the cerebral cortex. Leveraging an end-to-end unsupervised learning strategy, S3Reg offers great flexibility in the choice of input feature sets and output similarity measures for registration, and meanwhile reduces the registration time significantly. Specifically, we exploit the powerful learning capability of spherical Convolutional Neural Network (CNN) to directly learn the deformation fields in spherical space and implement diffeomorphic design with "scaling and squaring" layers to guarantee topology-preserving deformations. To handle the polar-distortion issue, we construct a novel spherical CNN model using three orthogonal Spherical U-Nets. Experiments are performed on two different datasets to align both adult and infant multimodal cortical features. Results demonstrate that our S3Reg shows superior or comparable performance with state-of-the-art methods, while improving the registration time from 1 min to 10 sec.

Retinotopic organization of visual cortex in human infants

Vision develops rapidly during infancy, yet how visual cortex is organized during this period is unclear. In particular, it is unknown whether functional maps that organize the mature adult visual cortex are present in the infant striate and extrastriate cortex. Here, we test the functional maturity of infant visual cortex by performing retinotopic mapping with functional magnetic resonance imaging (fMRI). Infants aged 5-23 months had retinotopic maps, with alternating preferences for vertical and horizontal meridians indicating the boundaries of visual areas V1 to V4 and an orthogonal gradient of preferences from high to low spatial frequencies. The presence of multiple visual maps throughout visual cortex in infants indicates a greater maturity of extrastriate cortex than previously appreciated. The areas showed subtle age-related fine-tuning, suggesting that early maturation undergoes continued refinement. This early maturation of area boundaries and tuning may scaffold subsequent developmental changes.

Developmental changes of the central sulcus morphology in young children

The human brain grows rapidly in early childhood, reaching 95% of its final volume by age 6. Understanding brain growth in childhood is important both to answer neuroscience questions about anatomical changes in development, and as a comparison metric for neurological disorders. Metrics for neuroanatomical development including cortical measures pertaining to the sulci can be instrumental in early diagnosis, monitoring, and intervention for neurological diseases. In this paper, we examine the development of the central sulcus in children aged 12-60 months from structural magnetic resonance images. The central sulcus is one of the earliest sulci to develop at the fetal stage and is implicated in diseases such as Attention Deficit Hyperactive Disorder and Williams syndrome. We investigate the relationship between the changes in the depth of the central sulcus with respect to age. In our results, we observed a pattern of depth present early on, that had been previously observed in adults. Results also reveal the presence of a rightward depth asymmetry at 12 months of age at a location related to orofacial movements. That asymmetry disappears gradually, mostly between 12 and 24 months, and we suggest that it is related to the development of language skills.

Cortical Structure and Cognition in Infants and Toddlers

Cortical structure has been consistently related to cognitive abilities in children and adults, yet we know little about how the cortex develops to support emergent cognition in infancy and toddlerhood when cortical thickness (CT) and surface area (SA) are maturing rapidly. In this report, we assessed how regional and global measures of CT and SA in a sample (N = 487) of healthy neonates, 1-year-olds, and 2-year-olds related to motor, language, visual reception, and general cognitive ability. We report novel findings that thicker cortices at ages 1 and 2 and larger SA at birth, age 1, and age 2 confer a cognitive advantage in infancy and toddlerhood. While several expected brain-cognition relationships were observed, overlapping cortical regions were also implicated across cognitive domains, suggesting that infancy marks a period of plasticity and refinement in cortical structure to support burgeoning motor, language, and cognitive abilities. CT may be a particularly important morphological indicator of ability, but its impact on cognition is relatively weak when compared with gestational age and maternal education. Findings suggest that prenatal and early postnatal cortical developments are important for cognition in infants and toddlers but should be considered in relation to other child and demographic factors.

The maturation and cognitive relevance of structural brain network organization from early infancy to childhood

The interactions of brain regions with other regions at the network level likely provide the infrastructure necessary for cognitive processes to develop. Specifically, it has been theorized that in infancy brain networks become more modular, or segregated, to support early cognitive specialization, before integration across networks increases to support the emergence of higher-order cognition. The present study examined the maturation of structural covariance networks (SCNs) derived from longitudinal cortical thickness data collected between infancy and childhood (0–6 years). We assessed modularity as a measure of network segregation and global efficiency as a measure of network integration. At the group level, we observed trajectories of increasing modularity and decreasing global efficiency between early infancy and six years. We further examined subject-based maturational coupling networks (sbMCNs) in a subset of this cohort with cognitive outcome data at 8–10 years, which allowed us to relate the network organization of longitudinal cortical thickness maturation to cognitive outcomes in middle childhood. We found that lower global efficiency of sbMCNs throughout early development (across the first year) related to greater motor learning at 8–10 years. Together, these results provide novel evidence characterizing the maturation of brain network segregation and integration across the first six years of life, and suggest that specific trajectories of brain network maturation contribute to later cognitive outcomes.

Multi-Site Infant Brain Segmentation Algorithms: The iSeg-2019 Challenge

To better understand early brain development in health and disorder, it is critical to accurately segment infant brain magnetic resonance (MR) images into white matter (WM), gray matter (GM), and cerebrospinal fluid (CSF). Deep learning-based methods have achieved state-of-the-art performance; h owever, one of the major limitations is that the learning-based methods may suffer from the multi-site issue, that is, the models trained on a dataset from one site may not be applicable to the datasets acquired from other sites with different imaging protocols/scanners. To promote methodological development in the community, the iSeg-2019 challenge (http://iseg2019.web.unc.edu) provides a set of 6-month infant subjects from multiple sites with different protocols/scanners for the participating methods. T raining/validation subjects are from UNC (MAP) and testing subjects are from UNC/UMN (BCP), Stanford University, and Emory University. By the time of writing, there are 30 automatic segmentation methods participated in the iSeg-2019. In this article, 8 top-ranked methods were reviewed by detailing their pipelines/implementations, presenting experimental results, and evaluating performance across different sites in terms of whole brain, regions of interest, and gyral landmark curves. We further pointed out their limitations and possible directions for addressing the multi-site issue. We find that multi-site consistency is still an open issue. We hope that the multi-site dataset in the iSeg-2019 and this review article will attract more researchers to address the challenging and critical multi-site issue in practice.

Spherical Deformable U-Net: Application to Cortical Surface Parcellation and Development Prediction

Convolutional Neural Networks (CNNs) have achieved overwhelming success in learning-related problems for 2D/3D images in the Euclidean space. However, unlike in the Euclidean space, the shapes of many structures in medical imaging have an inherent spherical topology in a manifold space, e.g., the convoluted brain cortical surfaces represented by triangular meshes. There is no consistent neighborhood definition and thus no straightforward convolution/pooling operations for such cortical surface data. In this paper, leveraging the regular and hierarchical geometric structure of the resampled spherical cortical surfaces, we create the 1-ring filter on spherical cortical triangular meshes and accordingly develop convolution/pooling operations for constructing Spherical U-Net for cortical surface data. However, the regular nature of the 1-ring filter makes it inherently limited to model fixed geometric transformations. To further enhance the transformation modeling capability of Spherical U-Net, we introduce the deformable convolution and deformable pooling to cortical surface data and accordingly propose the Spherical Deformable U-Net (SDU-Net). Specifically, spherical offsets are learned to freely deform the 1-ring filter on the sphere to adaptively localize cortical structures with different sizes and shapes. We then apply the SDU-Net to two challenging and scientifically important tasks in neuroimaging: cortical surface parcellation and cortical attribute map prediction. Both applications validate the competitive performance of our approach in accuracy and computational efficiency in comparison with state-of-the-art methods.

Robust Cortical Thickness Morphometry of Neonatal Brain and Systematic Evaluation Using Multi-Site MRI Datasets

The human brain grows the most dramatically during the perinatal and early post-natal periods, during which pre-term birth or perinatal injury that may alter brain structure and lead to developmental anomalies. Thus, characterizing cortical thickness of developing brains remains an important goal. However, this task is often complicated by inaccurate cortical surface extraction due to small-size brains. Here, we propose a novel complex framework for the reconstruction of neonatal WM and pial surfaces, accounting for large partial volumes due to small-size brains. The proposed approach relies only on T1-weighted images unlike previous T2-weighted image-based approaches while only T1-weighted images are sometimes available under the different clinical/research setting. Deep neural networks are first introduced to the neonatal magnetic resonance imaging (MRI) pipeline to address the mis-segmentation of brain tissues. Furthermore, this pipeline enhances cortical boundary delineation using combined models of the cerebrospinal fluid (CSF)/GM boundary detection with edge gradient information and a new skeletonization of sulcal folding where no CSF voxels are seen due to the limited resolution. We also proposed a systematic evaluation using three independent datasets comprising 736 pre-term and 97 term neonates. Qualitative assessment for reconstructed cortical surfaces shows that 86.9% are rated as accurate across the three site datasets. In addition, our landmark-based evaluation shows that the mean displacement of the cortical surfaces from the true boundaries was less than a voxel size (0.532 ± 0.035 mm). Evaluating the proposed pipeline (namely NEOCIVET 2.0) shows the robustness and reproducibility across different sites and different age-groups. The mean cortical thickness measured positively correlated with post-menstrual age (PMA) at scan (p < 0.0001); Cingulate cortical areas grew the most rapidly whereas the inferior temporal cortex grew the least rapidly. The range of the cortical thickness measured was biologically congruent (1.3 mm at 28 weeks of PMA to 1.8 mm at term equivalent). Cortical thickness measured on T1 MRI using NEOCIVET 2.0 was compared with that on T2 using the established dHCP pipeline. It was difficult to conclude that either T1 or T2 imaging is more ideal to construct cortical surfaces. NEOCIVET 2.0 has been open to the public through CBRAIN (https://mcin-cnim.ca/technology/cbrain/), a web-based platform for processing brain imaging data.

DIKA-Nets: Domain-invariant knowledge-guided attention networks for brain skull stripping of early developing macaques

As non-human primates, macaques have a close phylogenetic relationship to human beings and have been proven to be a valuable and widely used animal model in human neuroscience research. Accurate skull stripping (aka. brain extraction) of brain magnetic resonance imaging (MRI) is a crucial prerequisite in neuroimaging analysis of macaques. Most of the current skull stripping methods can achieve satisfactory results for human brains, but when applied to macaque brains, especially during early brain development, the results are often unsatisfactory. In fact, the early dynamic, regionally-heterogeneous development of macaque brains, accompanied by poor and age-related contrast between different anatomical structures, poses significant challenges for accurate skull stripping. To overcome these challenges, we propose a fully-automated framework to effectively fuse the age-specific intensity information and domain-invariant prior knowledge as important guiding information for robust skull stripping of developing macaques from 0 to 36 months of age. Specifically, we generate Signed Distance Map (SDM) and Center of Gravity Distance Map (CGDM) based on the intermediate segmentation results as guidance. Instead of using local convolution, we fuse all information using the Dual Self-Attention Module (DSAM), which can capture global spatial and channel-dependent information of feature maps. To extensively evaluate the performance, we adopt two relatively-large challenging MRI datasets from rhesus macaques and cynomolgus macaques, respectively, with a total of 361 scans from two different scanners with different imaging protocols. We perform cross-validation by using one dataset for training and the other one for testing. Our method outperforms five popular brain extraction tools and three deep-learning-based methods on cross-source MRI datasets without any transfer learning.

Multi-Regression based supervised sample selection for predicting baby connectome evolution trajectory from neonatal timepoint

The connectional map of the baby brain undergoes dramatic changes over the first year of postnatal development, which makes its mapping a challenging task, let alone learning how to predict its evolution. Currently, learning models for predicting brain connectomic developmental trajectories remain broadly absent despite their great potential in spotting atypical neurodevelopmental disorders early. This is most likely due to the scarcity and often incompleteness of longitudinal infant neuroimaging studies for training such models. In this paper, we propose the first approach for progressively predicting longitudinal development of brain networks during the postnatal period solely from a baseline connectome around birth. To this end, a supervised multi-regression sample selection strategy is designed to learn how to identify the best set of neighbors of a testing baseline connectome to eventually predict its evolution trajectory at follow-up timepoints. However, given that the training dataset may have missing samples (connectomes) at certain timepoints, this may affect the training of the predictive model. To overcome this problem, we perform a low-rank tensor completion based on a robust principal component analysis to impute the missing training connectomes by linearly approximating similar complete training networks. In the prediction step, our sample selection strategy aims to preserve spatiotemporal relationships between consecutive timepoints. Therefore, the proposed method learns how to identify the set of the local closest neighbors to a target network by training an ensemble of bidirectional regressors leveraging temporal dependency between consecutive timepoints with a recall to the baseline observations to progressively predict the evolution of a testing network over time. Our method achieves the best prediction results and better captures the dynamic changes of each brain connectome over time in comparison to its ablated versions using leave-one-out cross-validation strategy.

Longitudinal Mapping of Cortical Thickness Measurements: An Alzheimer's Disease Neuroimaging Initiative-Based Evaluation Study

Longitudinal studies of development and disease in the human brain have motivated the acquisition of large neuroimaging data sets and the concomitant development of robust methodological and statistical tools for quantifying neurostructural changes. Longitudinal-specific strategies for acquisition and processing have potentially significant benefits including more consistent estimates of intra-subject measurements while retaining predictive power. Using the first phase of the Alzheimer's Disease Neuroimaging Initiative (ADNI-1) data, comprising over 600 subjects with multiple time points from baseline to 36 months, we evaluate the utility of longitudinal FreeSurfer and Advanced Normalization Tools (ANTs) surrogate thickness values in the context of a linear mixed-effects (LME) modeling strategy. Specifically, we estimate the residual variability and between-subject variability associated with each processing stream as it is known from the statistical literature that minimizing the former while simultaneously maximizing the latter leads to greater scientific interpretability in terms of tighter confidence intervals in calculated mean trends, smaller prediction intervals, and narrower confidence intervals for determining cross-sectional effects. This strategy is evaluated over the entire cortex, as defined by the Desikan-Killiany-Tourville labeling protocol, where comparisons are made with the cross-sectional and longitudinal FreeSurfer processing streams. Subsequent linear mixed effects modeling for identifying diagnostic groupings within the ADNI cohort is provided as supporting evidence for the utility of the proposed ANTs longitudinal framework which provides unbiased structural neuroimage processing and competitive to superior power for longitudinal structural change detection.

Neurocranium thickness mapping in early childhood

The neurocranium changes rapidly in early childhood to accommodate the growing brain. Developmental disorders and environmental factors such as sleep position may lead to abnormal neurocranial maturation. Therefore, it is important to understand how this structure develops, in order to provide a baseline for early detection of anomalies. However, its anatomy has not yet been well studied in early childhood due to the lack of available imaging databases. In hospitals, CT is typically used to image the neurocranium when a pathology is suspected, but the presence of ionizing radiation makes it harder to construct databases of healthy subjects. In this study, instead, we use a dataset of MRI data from healthy normal children in the age range of 6 months to 36 months to study the development of the neurocranium. After extracting its outline from the MRI data, we used a conformal geometry-based analysis pipeline to detect local thickness growth throughout this age span. These changes will help us understand cranial bone development with respect to the brain, as well as detect abnormal variations, which will in turn inform better treatment strategies for implicated disorders.

Semi-supervised Transfer Learning for Infant Cerebellum Tissue Segmentation

To characterize early cerebellum development, accurate segmentation of the cerebellum into white matter (WM), gray matter (GM), and cerebrospinal fluid (CSF) tissues is one of the most pivotal steps. However, due to the weak tissue contrast, extremely folded tiny structures, and severe partial volume effect, infant cerebellum tissue segmentation is especially challenging, and the manual labels are hard to obtain and correct for learning-based methods. To the best of our knowledge, there is no work on the cerebellum segmentation for infant subjects less than 24 months of age. In this work, we develop a semi-supervised transfer learning framework guided by a confidence map for tissue segmentation of cerebellum MR images from 24-month-old to 6-month-old infants. Note that only 24-month-old subjects have reliable manual labels for training, due to their high tissue contrast. Through the proposed semi-supervised transfer learning, the labels from 24-month-old subjects are gradually propagated to the 18-, 12-, and 6-month-old subjects, which have a low tissue contrast. Comparison with the state-of-the-art methods demonstrates the superior performance of the proposed method, especially for 6-month-old subjects.

Unsupervised Learning for Spherical Surface Registration

Current spherical surface registration methods achieve good performance on alignment and spatial normalization of cortical surfaces across individuals in neuroimaging analysis. However, they are computationally intensive, since they have to optimize an objective function independently for each pair of surfaces. In this paper, we present a fast learning-based algorithm that makes use of the recent development in spherical Convolutional Neural Networks (CNNs) for spherical cortical surface registration. Given a set of surface pairs without supervised information such as ground truth deformation fields or anatomical landmarks, we formulate the registration as a parametric function and learn its parameters by enforcing the feature similarity between one surface and the other one warped by the estimated deformation field using the function. Then, given a new pair of surfaces, we can quickly infer the spherical deformation field registering one surface to the other one. We model this parametric function using three orthogonal Spherical U-Nets and use spherical transform layers to warp the spherical surfaces, while imposing smoothness constraints on the deformation field. All the layers in the network are well-defined and differentiable, thus the parameters can be effectively learned. We show that our method achieves accurate cortical alignment results on 102 subjects, comparable to two state-of-the-art methods: Spherical Demons and MSM, while runs much faster.

Disentangled Intensive Triplet Autoencoder for Infant Functional Connectome Fingerprinting

Functional connectome "fingerprint" is a highly characterized brain pattern that distinguishes one individual from others. Although its existence has been demonstrated in adults, an unanswered but fundamental question is whether such individualized pattern emerges since infancy. This problem is barely investigated despites its importance in identifying the origin of the intrinsic connectome patterns that mirror distinct behavioral phenotypes. However, addressing this knowledge gap is challenging because the conventional methods are only applicable to developed brains with subtle longitudinal changes and typically fail on the dramatically developing infant brains. To tackle this challenge, we invent a novel model, namely, disentangled intensive triplet autoencoder (DI-TAE). First, we introduce the triplet autoencoder to embed the original connectivity into a latent space with higher discriminative capability among infant individuals. Then, a disentanglement strategy is proposed to separate the latent variables into identity-code, age-code, and noise-code, which not only restrains the interference from age-related developmental variance, but also captures the identity-related invariance. Next, a cross-reconstruction loss and an intensive triplet loss are designed to guarantee the effectiveness of the disentanglement and enhance the inter-subject dissimilarity for better discrimination. Finally, a variance-guided bootstrap aggregating is developed for DI-TAE to further improve the performance of identification. DI-TAE is validated on three longitudinal resting-state fMRI datasets with 394 infant scans aged 16 to 874 days. Our proposed model outperforms other state-of-the-art methods by increasing the identification rate by more than 50%, and for the first time suggests the plausible existence of brain functional connectome "fingerprint" since early infancy.

Infant FreeSurfer: An automated segmentation and surface extraction pipeline for T1-weighted neuroimaging data of infants 0-2 years

The development of automated tools for brain morphometric analysis in infants has lagged significantly behind analogous tools for adults. This gap reflects the greater challenges in this domain due to: 1) a smaller-scaled region of interest, 2) increased motion corruption, 3) regional changes in geometry due to heterochronous growth, and 4) regional variations in contrast properties corresponding to ongoing myelination and other maturation processes. Nevertheless, there is a great need for automated image-processing tools to quantify differences between infant groups and other individuals, because aberrant cortical morphologic measurements (including volume, thickness, surface area, and curvature) have been associated with neuropsychiatric, neurologic, and developmental disorders in children. In this paper we present an automated segmentation and surface extraction pipeline designed to accommodate clinical MRI studies of infant brains in a population 0-2 year-olds. The algorithm relies on a single channel of T1-weighted MR images to achieve automated segmentation of cortical and subcortical brain areas, producing volumes of subcortical structures and surface models of the cerebral cortex. We evaluated the algorithm both qualitatively and quantitatively using manually labeled datasets, relevant comparator software solutions cited in the literature, and expert evaluations. The computational tools and atlases described in this paper will be distributed to the research community as part of the FreeSurfer image analysis package.

Neural and psychophysiological correlates of social communication development: Evidence from sensory processing, motor, cognitive, language and emotional behavioral milestones across infancy

This article presents a literature review focusing on the neural and psychophysiological correlates associated with social communication development in infancy. Studies presenting evidence on infants' brain activity and developments in infant sensory processing, motor, cognitive, language, and emotional abilities are described in regard to the neuropsychophysiological processes underlying the emergence of these specific behavioral milestones and their associations with social communication development. Studies that consider specific age-related characteristics across the infancy period are presented. Evidence suggests that specific neural and physiological signatures accompany age-related social communication development during the first 18 months of life.

Environmental Influences on Infant Cortical Thickness and Surface Area

Cortical thickness (CT) and surface area (SA) vary widely between individuals and are associated with intellectual ability and risk for various psychiatric and neurodevelopmental conditions. Factors influencing this variability remain poorly understood, but the radial unit hypothesis, as well as the more recent supragranular cortex expansion hypothesis, suggests that prenatal and perinatal influences may be particularly important. In this report, we examine the impact of 17 major demographic and obstetric history variables on interindividual variation in CT and SA in a unique sample of 805 neonates who received MRI scans of the brain around 2 weeks of age. Birth weight, postnatal age at MRI, gestational age at birth, and sex emerged as important predictors of SA. Postnatal age at MRI, paternal education, and maternal ethnicity emerged as important predictors of CT. These findings suggest that individual variation in infant CT and SA is explained by different sets of environmental factors with neonatal SA more strongly influenced by sex and obstetric history and CT more strongly influenced by socioeconomic and ethnic disparities. Findings raise the possibility that interventions aimed at reducing disparities and improving obstetric outcomes may alter prenatal/perinatal cortical development.

Individual identification and individual variability analysis based on cortical folding features in developing infant singletons and twins

Studying the early dynamic development of cortical folding with remarkable individual variability is critical for understanding normal early brain development and related neurodevelopmental disorders. This study focuses on the fingerprinting capability and the individual variability of cortical folding during early brain development. Specifically, we aim to explore (a) whether the developing neonatal cortical folding is unique enough to be considered as a "fingerprint" that can reliably identify an individual within a cohort of infants; (b) which cortical regions manifest more individual variability and thus contribute more for infant identification; (c) whether the infant twins can be distinguished by cortical folding. Hence, for the first time, we conduct infant individual identification and individual variability analysis involving twins based on the developing cortical folding features (mean curvature, average convexity, and sulcal depth) in 472 neonates with 1,141 longitudinal MRI scans. Experimental results show that the infant individual identification achieves 100% accuracy when using the neonatal cortical folding features to predict the identities of 1- and 2-year-olds. Besides, we observe high identification capability in the high-order association cortices (i.e., prefrontal, lateral temporal, and inferior parietal regions) and two unimodal cortices (i.e., precentral gyrus and lateral occipital cortex), which largely overlap with the regions encoding remarkable individual variability in cortical folding during the first 2 years. For twins study, we show that even for monozygotic twins with identical genes and similar developmental environments, their cortical folding features are unique enough for accurate individual identification; and in some high-order association cortices, the differences between monozygotic twin pairs are significantly lower than those between dizygotic twins. This study thus provides important insights into individual identification and individual variability based on cortical folding during infancy.

Structural neuroimaging

Characterizing the neuroanatomical correlates of brain development is essential in understanding brain-behavior relationships and neurodevelopmental disorders. Advances in brain MRI acquisition protocols and image processing techniques have made it possible to detect and track with great precision anatomical brain development and pediatric neurologic disorders. In this chapter, we provide a brief overview of the modern neuroimaging techniques for pediatric brain development and review key normal brain development studies. Characteristic disorders affecting neurodevelopment in childhood, such as prematurity, attention deficit hyperactivity disorder (ADHD), autism spectrum disorder (ASD), epilepsy, and brain cancer, and key neuroanatomical findings are described and then reviewed. Large datasets of typically developing children and children with various neurodevelopmental conditions are now being acquired to help provide the biomarkers of such impairments. While there are still several challenges in imaging brain structures specific to the pediatric populations, such as subject cooperation and tissues contrast variability, considerable imaging research is now being devoted to solving these problems and improving pediatric data analysis.