Computational neuroanatomy of baby brains: A review

Intrinsic Infant Hippocampal Function Supports Inhibitory Processing

Impaired cerebral inhibition is commonly observed in neurodevelopmental disorders and may represent a vulnerability factor for their development. The hippocampus plays a key role in inhibition among adults and undergoes significant and rapid changes during early brain development. Therefore, the structure represents an important candidate region for early identification of pathology that is relevant to inhibitory dysfunction. To determine whether hippocampal function corresponds to inhibition in the early postnatal period, the present study evaluated relationships between hippocampal activity and sensory gating in infants 4-20 weeks of age (N = 18). Resting-state functional magnetic resonance imaging was used to measure hippocampal activity, including the amplitude of low-frequency fluctuations (ALFFs) and fractional ALFF. Electroencephalography during a paired-stimulus paradigm was used to measure sensory gating (P50). Higher activity of the right hippocampus was associated with better sensory gating (P50 ratio), driven by a reduction in response to the second stimulus. These findings suggest that meaningful effects of hippocampal function can be detected early in infancy. Specifically, higher intrinsic hippocampal activity in the early postnatal period may support effective inhibitory processing. Future work will benefit from longitudinal analysis to clarify the trajectory of hippocampal function, alterations of which may contribute to the risk of neurodevelopmental disorders and represent an intervention target.

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.

Longitudinal trajectories of brain development from infancy to school age and their relationship to literacy development

Reading is one of the most complex skills that we utilize daily, and it involves the early development and interaction of various lower-level subskills, including phonological processing and oral language. These subskills recruit brain structures, which begin to develop long before the skill manifests and exhibit rapid development during infancy. However, how longitudinal trajectories of early brain development in these structures supports long-term acquisition of literacy subskills and subsequent reading is unclear. Children underwent structural and diffusion MRI scanning at multiple timepoints between infancy and second grade and were tested for literacy subskills in preschool and decoding and word reading in early elementary school. We developed and implemented a reproducible pipeline to generate longitudinal trajectories of early brain development to examine associations between these trajectories and literacy (sub)skills. Furthermore, we examined whether familial risk of reading difficulty and a child's home literacy environment, two common literacy-related covariates, influenced those trajectories. Results showed that individual differences in curve features (e.g., intercepts and slopes) for longitudinal trajectories of volumetric, surface-based, and white matter organization measures in left-hemispheric reading-related regions and tracts were linked directly to phonological processing and indirectly to second-grade decoding and word reading skills via phonological processing. Altogether, these findings suggest that the brain bases of phonological processing, previously identified as the strongest behavioral predictor of reading and decoding skills, may already begin to develop early in infancy but undergo further refinement between birth and preschool. The present study underscores the importance of considering academic skill acquisition from the very beginning of life.

Constructing fine-grained spatiotemporal neonatal functional atlases with spectral functional network learning

The early stages of human development are increasingly acknowledged as pivotal in laying the groundwork for subsequent behavioral and cognitive development. Spatiotemporal (4D) brain functional atlases are important in elucidating the development of human brain functions. However, the scarcity of such atlases for early life stages stems from two primary challenges: (1) the significant noise in functional magnetic resonance imaging (fMRI) that complicates the generation of high-quality atlases for each age group, and (2) the rapid and complex changes in the early human brain that hinder the maintenance of temporal consistency in 4D atlases. This study tackles these challenges by integrating low-rank tensor learning with spectral embedding, thereby proposing a novel, data-driven 4D functional atlas generation framework based on spectral functional network learning (SFNL). This method utilizes low-rank tensor learning to capture common functional connectivity (FC) patterns across different ages, thus optimizing FCs for each age group to improve the temporal consistency of functional networks. Incorporating spectral embedding aids in mitigating potential noise in FC networks derived from fMRI data by reconstructing networks in the spectral space. Utilizing SFNL-generated functional networks enables the creation of consistent and highly qualified spatiotemporal functional atlases. The framework was applied to the developing Human Connectome Project (dHCP) dataset, generating the first neonatal 4D functional atlases with fine-grained temporal and spatial resolutions. Experimental evaluations focusing on functional homogeneity, reliability, and temporal consistency demonstrated the superiority of our framework compared to existing methods for constructing 4D atlases. Additionally, network analysis experiments, including individual identification, functional systems development, and local efficiency assessments, further corroborate the efficacy and robustness of the generated atlases. The 4D atlases and related codes will be made publicly accessible (https://github.com/zhaoyunxi/neonate-atlases).

Early autism diagnosis based on path signature and Siamese unsupervised feature compressor

Autism spectrum disorder has been emerging as a growing public health threat. Early diagnosis of autism spectrum disorder is crucial for timely, effective intervention and treatment. However, conventional diagnosis methods based on communications and behavioral patterns are unreliable for children younger than 2 years of age. Given evidences of neurodevelopmental abnormalities in autism spectrum disorder infants, we resort to a novel deep learning-based method to extract key features from the inherently scarce, class-imbalanced, and heterogeneous structural MR images for early autism diagnosis. Specifically, we propose a Siamese verification framework to extend the scarce data, and an unsupervised compressor to alleviate data imbalance by extracting key features. We also proposed weight constraints to cope with sample heterogeneity by giving different samples different voting weights during validation, and used Path Signature to unravel meaningful developmental features from the two-time point data longitudinally. We further extracted machine learning focused brain regions for autism diagnosis. Extensive experiments have shown that our method performed well under practical scenarios, transcending existing machine learning methods and providing anatomical insights for autism early diagnosis.

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.

Comprehensive analysis of synthetic learning applied to neonatal brain MRI segmentation

Brain segmentation from neonatal MRI images is a very challenging task due to large changes in the shape of cerebral structures and variations in signal intensities reflecting the gestational process. In this context, there is a clear need for segmentation techniques that are robust to variations in image contrast and to the spatial configuration of anatomical structures. In this work, we evaluate the potential of synthetic learning, a contrast-independent model trained using synthetic images generated from the ground truth labels of very few subjects. We base our experiments on the dataset released by the developmental Human Connectome Project, for which high-quality images are available for more than 700 babies aged between 26 and 45 weeks postconception. First, we confirm the impressive performance of a standard UNet trained on a few volumes, but also confirm that such models learn intensity-related features specific to the training domain. We then confirm the robustness of the synthetic learning approach to variations in image contrast. However, we observe a clear influence of the age of the baby on the predictions. We improve the performance of this model by enriching the synthetic training set with realistic motion artifacts and over-segmentation of the white matter. Based on extensive visual assessment, we argue that the better performance of the model trained on real T2w data may be due to systematic errors in the ground truth. We propose an original experiment allowing us to show that learning from real data will reproduce any systematic bias affecting the training set, while synthetic models can avoid this limitation. Overall, our experiments confirm that synthetic learning is an effective solution for segmenting neonatal brain MRI. Our adapted synthetic learning approach combines key features that will be instrumental for large multisite studies and clinical applications.

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.

Deep learning-assisted preclinical MR fingerprinting for sub-millimeter T1 and T2 mapping of entire macaque brain

PURPOSE

Preclinical MR fingerprinting (MRF) suffers from long acquisition time for organ-level coverage due to demanding image resolution and limited undersampling capacity. This study aims to develop a deep learning-assisted fast MRF framework for sub-millimeter T1 and T2 mapping of entire macaque brain on a preclinical 9.4 T MR system.

METHODS

Three dimensional MRF images were reconstructed by singular value decomposition (SVD) compressed reconstruction. T1 and T2 mapping for each axial slice exploited a self-attention assisted residual U-Net to suppress aliasing-induced quantification errors, and the transmit-field (B1 + ) measurements for robustness against B1 + inhomogeneity. Supervised network training used MRF images simulated via virtual parametric maps and a desired undersampling scheme. This strategy bypassed the difficulties of acquiring fully sampled preclinical MRF data to guide network training. The proposed fast MRF framework was tested on experimental data acquired from ex vivo and in vivo macaque brains.

RESULTS

The trained network showed reasonable adaptability to experimental MRF images, enabling robust delineation of various T1 and T2 distributions in the brain tissues. Further, the proposed MRF framework outperformed several existing fast MRF methods in handling the aliasing artifacts and capturing detailed cerebral structures in the mapping results. Parametric mapping of entire macaque brain at nominal resolution of 0.35× $$ \times $$ 0.35× $$ \times $$ 1 mm3 can be realized via a 20-min 3D MRF scan, which was sixfold faster than the baseline protocol.

CONCLUSION

Introducing deep learning to MRF framework paves the way for efficient organ-level high-resolution quantitative MRI in preclinical applications.

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.

Fully automated measurement of intracranial CSF and brain parenchyma volumes in pediatric hydrocephalus by segmentation of clinical MRI studies

BACKGROUND

Brain parenchyma (BP) and intracranial cerebrospinal fluid (iCSF) volumes measured by fully automated segmentation of clinical brain MRI studies may be useful for the diagnosis and follow-up of pediatric hydrocephalus. However, previously published segmentation techniques either rely on dedicated sequences, not routinely used in clinical practice, or on spatial normalization, which has limited accuracy when severe brain distortions, such as in hydrocephalic patients, are present.

PURPOSE

We developed a fully automated method to measure BP and iCSF volumes from clinical brain MRI studies of pediatric hydrocephalus patients, exploiting the complementary information contained in T2- and T1-weighted images commonly used in clinical practice.

METHODS

The proposed procedure, following skull-stripping of the combined volumes, performed using a multiparametric method to obtain a reliable definition of the inner skull profile, maximizes the CSF-to-parenchyma contrast by dividing the T2w- by the T1w- volume after full-scale dynamic rescaling, thus allowing separation of iCSF and BP through a simple thresholding routine.

RESULTS

Validation against manual tracing on 23 studies (four controls and 19 hydrocephalic patients) showed excellent concordance (ICC > 0.98) and spatial overlap (Dice coefficients ranging from 77.2% for iCSF to 96.8% for intracranial volume). Accuracy was comparable to the intra-operator reproducibility of manual segmentation, as measured in 14 studies processed twice by the same experienced neuroradiologist. Results of the application of the algorithm to a dataset of 63 controls and 57 hydrocephalic patients (19 with parenchymal damage), measuring volumes' changes with normal development and in hydrocephalic patients, are also reported for demonstration purposes.

CONCLUSIONS

The proposed approach allows fully automated segmentation of BP and iCSF in clinical studies, also in severely distorted brains, enabling to assess age- and disease-related changes in intracranial tissue volume with an accuracy comparable to expert manual segmentation.

Early development of the functional brain network in newborns

During the prenatal period and the first postnatal years, the human brain undergoes rapid growth, which establishes a preliminary infrastructure for the subsequent development of cognition and behavior. To understand the underlying processes of brain functioning and identify potential sources of developmental disorders, it is essential to uncover the developmental rules that govern this critical period. In this study, graph theory modeling and network science analysis were employed to investigate the impact of age, gender, weight, and typical and atypical development on brain development. Local and global topologies of functional connectomes obtained from rs-fMRI data were collected from 421 neonates aged between 31 and 45 postmenstrual weeks who were in natural sleep without any sedation. The results showed that global efficiency, local efficiency, clustering coefficient, and small-worldness increased with age, while modularity and characteristic path length decreased with age. The normalized rich-club coefficient displayed a U-shaped pattern during development. The study also examined the global and local impacts of gender, weight, and group differences between typical and atypical cases. The findings presented some new insights into the maturation of functional brain networks and their relationship with cognitive development and neurodevelopmental disorders.

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.

Statistical Learning Methods for Neuroimaging Data Analysis with Applications

The aim of this review is to provide a comprehensive survey of statistical challenges in neuroimaging data analysis, from neuroimaging techniques to large-scale neuroimaging studies and statistical learning methods. We briefly review eight popular neuroimaging techniques and their potential applications in neuroscience research and clinical translation. We delineate four themes of neuroimaging data and review major image processing analysis methods for processing neuroimaging data at the individual level. We briefly review four large-scale neuroimaging-related studies and a consortium on imaging genomics and discuss four themes of neuroimaging data analysis at the population level. We review nine major population-based statistical analysis methods and their associated statistical challenges and present recent progress in statistical methodology to address these challenges.

Longitudinal mapping of the development of cortical thickness and surface area in rhesus macaques during the first three years

Studying dynamic spatiotemporal patterns of early brain development in macaque monkeys is critical for understanding the cortical organization and evolution in humans, given the phylogenetic closeness between humans and macaques. However, due to huge challenges in the analysis of early brain Magnetic Resonance Imaging (MRI) data typically with extremely low contrast and dynamic imaging appearances, our knowledge of the early macaque cortical development remains scarce. To fill this critical gap, this paper characterizes the early developmental patterns of cortical thickness and surface area in rhesus macaques by leveraging advanced computing tools tailored for early developing brains based on a densely sampled longitudinal dataset with 140 rhesus macaque MRI scans seamlessly covering from birth to 36 mo of age. The average cortical thickness exhibits an inverted U-shaped trajectory with peak thickness at around 4.3 mo of age, which is remarkably in line with the age of peak thickness at 14 mo in humans, considering the around 3:1 age ratio of human to macaque. The total cortical surface area in macaques increases monotonically but with relatively lower expansions than in humans. The spatial distributions of thicker and thinner regions are quite consistent during development, with gyri having a thicker cortex than sulci. By 4 mo of age, over 81% of cortical vertices have reached their peaks in thickness, except for the insula and medial temporal cortices, while most cortical vertices keep expanding in surface area, except for the occipital cortex. These findings provide important insights into early brain development and evolution in primates.

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.

BIBSNet: A Deep Learning Baby Image Brain Segmentation Network for MRI Scans

Objectives

Brain segmentation of infant magnetic resonance (MR) images is vitally important in studying developmental mental health and disease. The infant brain undergoes many changes throughout the first years of postnatal life, making tissue segmentation difficult for most existing algorithms. Here, we introduce a deep neural network BIBSNet (Baby and Infant Brain Segmentation Neural Network), an open-source, community-driven model that relies on data augmentation and a large sample size of manually annotated images to facilitate the production of robust and generalizable brain segmentations.

Experimental Design

Included in model training and testing were MR brain images on 84 participants with an age range of 0-8 months (median postmenstrual ages of 13.57 months). Using manually annotated real and synthetic segmentation images, the model was trained using a 10-fold cross-validation procedure. Testing occurred on MRI data processed with the DCAN labs infant-ABCD-BIDS processing pipeline using segmentations produced from gold standard manual annotation, joint-label fusion (JLF), and BIBSNet to assess model performance.

Principal Observations

Using group analyses, results suggest that cortical metrics produced using BIBSNet segmentations outperforms JLF segmentations. Additionally, when analyzing individual differences, BIBSNet segmentations perform even better.

Conclusions

BIBSNet segmentation shows marked improvement over JLF segmentations across all age groups analyzed. The BIBSNet model is 600x faster compared to JLF and can be easily included in other processing pipelines.

Development of cortical folds in the human brain: An attempt to review biological hypotheses, early neuroimaging investigations and functional correlates

The folding of the human brain mostly takes place in utero, making it challenging to study. After a few pioneer studies looking into it in post-mortem foetal specimen, modern approaches based on neuroimaging have allowed the community to investigate the folding process in vivo, its normal progression, its early disturbances, and its relationship to later functional outcomes. In this review article, we aimed to first give an overview of the current hypotheses on the mechanisms governing cortical folding. After describing the methodological difficulties raised by its study in fetuses, neonates and infants with magnetic resonance imaging (MRI), we reported our current understanding of sulcal pattern emergence in the developing brain. We then highlighted the functional relevance of early sulcal development, through recent insights about hemispheric asymmetries and early factors influencing this dynamic such as prematurity. Finally, we outlined how longitudinal studies have started to relate early folding markers and the child's sensorimotor and cognitive outcome. Through this review, we hope to raise awareness on the potential of studying early sulcal patterns both from a fundamental and clinical perspective, as a window into early neurodevelopment and plasticity in relation to growth in utero and postnatal environment of the child.

An attention-based context-informed deep framework for infant brain subcortical segmentation

Precise segmentation of subcortical structures from infant brain magnetic resonance (MR) images plays an essential role in studying early subcortical structural and functional developmental patterns and diagnosis of related brain disorders. However, due to the dynamic appearance changes, low tissue contrast, and tiny subcortical size in infant brain MR images, infant subcortical segmentation is a challenging task. In this paper, we propose a context-guided, attention-based, coarse-to-fine deep framework to precisely segment the infant subcortical structures. At the coarse stage, we aim to directly predict the signed distance maps (SDMs) from multi-modal intensity images, including T1w, T2w, and the ratio of T1w and T2w images, with an SDM-Unet, which can leverage the spatial context information, including the structural position information and the shape information of the target structure, to generate high-quality SDMs. At the fine stage, the predicted SDMs, which encode spatial-context information of each subcortical structure, are integrated with the multi-modal intensity images as the input to a multi-source and multi-path attention Unet (M2A-Unet) for achieving refined segmentation. Both the 3D spatial and channel attention blocks are added to guide the M2A-Unet to focus more on the important subregions and channels. We additionally incorporate the inner and outer subcortical boundaries as extra labels to help precisely estimate the ambiguous boundaries. We validate our method on an infant MR image dataset and on an unrelated neonatal MR image dataset. Compared to eleven state-of-the-art methods, the proposed framework consistently achieves higher segmentation accuracy in both qualitative and quantitative evaluations of infant MR images and also exhibits good generalizability in the neonatal dataset.

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).

iBEAT V2.0: a multisite-applicable, deep learning-based pipeline for infant cerebral cortical surface reconstruction

The human cerebral cortex undergoes dramatic and critical development during early postnatal stages. Benefiting from advances in neuroimaging, many infant brain magnetic resonance imaging (MRI) datasets have been collected from multiple imaging sites with different scanners and imaging protocols for the investigation of normal and abnormal early brain development. However, it is extremely challenging to precisely process and quantify infant brain development with these multisite imaging data because infant brain MRI scans exhibit (a) extremely low and dynamic tissue contrast caused by ongoing myelination and maturation and (b) inter-site data heterogeneity resulting from the use of diverse imaging protocols/scanners. Consequently, existing computational tools and pipelines typically perform poorly on infant MRI data. To address these challenges, we propose a robust, multisite-applicable, infant-tailored computational pipeline that leverages powerful deep learning techniques. The main functionality of the proposed pipeline includes preprocessing, brain skull stripping, tissue segmentation, topology correction, cortical surface reconstruction and measurement. Our pipeline can handle both T1w and T2w structural infant brain MR images well in a wide age range (from birth to 6 years of age) and is effective for different imaging protocols/scanners, despite being trained only on the data from the Baby Connectome Project. Extensive comparisons with existing methods on multisite, multimodal and multi-age datasets demonstrate superior effectiveness, accuracy and robustness of our pipeline. We have maintained a website, iBEAT Cloud, for users to process their images with our pipeline ( http://www.ibeat.cloud ), which has successfully processed over 16,000 infant MRI scans from more than 100 institutions with various imaging protocols/scanners.

3D-MASNet: 3D mixed-scale asymmetric convolutional segmentation network for 6-month-old infant brain MR images

Precise segmentation of infant brain magnetic resonance (MR) images into gray matter (GM), white matter (WM), and cerebrospinal fluid (CSF) are essential for studying neuroanatomical hallmarks of early brain development. However, for 6-month-old infants, the extremely low-intensity contrast caused by inherent myelination hinders accurate tissue segmentation. Existing convolutional neural networks (CNNs) based segmentation models for this task generally employ single-scale symmetric convolutions, which are inefficient for encoding the isointense tissue boundaries in baby brain images. Here, we propose a 3D mixed-scale asymmetric convolutional segmentation network (3D-MASNet) framework for brain MR images of 6-month-old infants. We replaced the traditional convolutional layer of an existing to-be-trained network with a 3D mixed-scale convolution block consisting of asymmetric kernels (MixACB) during the training phase and then equivalently converted it into the original network. Five canonical CNN segmentation models were evaluated using both T1- and T2-weighted images of 23 6-month-old infants from iSeg-2019 datasets, which contained manual labels as ground truth. MixACB significantly enhanced the average accuracy of all five models and obtained the most considerable improvement in the fully convolutional network model (CC-3D-FCN) and the highest performance in the Dense U-Net model. This approach further obtained Dice coefficient accuracies of 0.931, 0.912, and 0.961 in GM, WM, and CSF, respectively, ranking first among 30 teams on the validation dataset of the iSeg-2019 Grand Challenge. Thus, the proposed 3D-MASNet can improve the accuracy of existing CNNs-based segmentation models as a plug-and-play solution that offers a promising technique for future infant brain MRI studies.

Computational methods for the analysis of early-pregnancy brain ultrasonography: a systematic review

BACKGROUND

Early screening of the brain is becoming routine clinical practice. Currently, this screening is performed by manual measurements and visual analysis, which is time-consuming and prone to errors. Computational methods may support this screening. Hence, the aim of this systematic review is to gain insight into future research directions needed to bring automated early-pregnancy ultrasound analysis of the human brain to clinical practice.

METHODS

We searched PubMed (Medline ALL Ovid), EMBASE, Web of Science Core Collection, Cochrane Central Register of Controlled Trials, and Google Scholar, from inception until June 2022. This study is registered in PROSPERO at CRD42020189888. Studies about computational methods for the analysis of human brain ultrasonography acquired before the 20th week of pregnancy were included. The key reported attributes were: level of automation, learning-based or not, the usage of clinical routine data depicting normal and abnormal brain development, public sharing of program source code and data, and analysis of the confounding factors.

FINDINGS

Our search identified 2575 studies, of which 55 were included. 76% used an automatic method, 62% a learning-based method, 45% used clinical routine data and in addition, for 13% the data depicted abnormal development. None of the studies shared publicly the program source code and only two studies shared the data. Finally, 35% did not analyse the influence of confounding factors.

INTERPRETATION

Our review showed an interest in automatic, learning-based methods. To bring these methods to clinical practice we recommend that studies: use routine clinical data depicting both normal and abnormal development, make their dataset and program source code publicly available, and be attentive to the influence of confounding factors. Introduction of automated computational methods for early-pregnancy brain ultrasonography will save valuable time during screening, and ultimately lead to better detection, treatment and prevention of neuro-developmental disorders.

FUNDING

The Erasmus MC Medical Research Advisor Committee (grant number: FB 379283).

Associations between cortical thickness and anxious/depressive symptoms differ by the quality of early care

A variety of childhood experiences can lead to anxious/depressed (A/D) symptoms. The aim of the present study was to explore the brain morphological (cortical thickness and surface area) correlates of A/D symptoms and the extent to which these phenotypes vary depending on the quality of the parenting context in which children develop. Structural magnetic resonance imaging (MRI) scans were acquired on 45 children with Child Protective Services (CPS) involvement due to risk of not receiving adequate care (high-risk group) and 25 children without CPS involvement (low-risk group) (range_age = 8.08-12.14; M_age = 10.05) to assess cortical thickness (CT) and cortical surface area (SA). A/D symptoms were measured using the Child Behavioral Checklist. The association between A/D symptoms and CT, but not SA, differed by risk status such that high-risk children showed decreasing CT as A/D scores increased, whereas low-risk children showed increasing CT as A/D scores increased. This interaction was specific to CT in prefrontal, frontal, temporal, and parietal cortical regions. The groups had marginally different A/D scores, in the direction of higher risk being associated with lower A/D scores. Results suggest that CT correlates of A/D symptoms are differentially shaped by the quality of early caregiving experiences and should be distinguished between high- and low-risk children.

Template-based graph registration network for boosting the diagnosis of brain connectivity disorders

Brain graphs are powerful representations to explore the biological roadmaps of the human brain in its healthy and disordered states. Recently, a few graph neural networks (GNNs) have been designed for brain connectivity synthesis and diagnosis. However, such non-Euclidean deep learning architectures might fail to capture the neural interactions between different brain regions as they are trained without guidance from any prior biological template-i.e., template-free learning. Here we assume that using a population-driven brain connectional template (CBT) that captures well the connectivity patterns fingerprinting a given brain state (e.g., healthy) can better guide the GNN training in its downstream learning task such as classification or regression. To this aim we design a plug-in graph registration network (GRN) that can be coupled with any conventional graph neural network (GNN) so as to boost its learning accuracy and generalizability to unseen samples. Our GRN is a graph generative adversarial network (gGAN), which registers brain graphs to a prior CBT. Next, the registered brain graphs are used to train typical GNN models. Our GRN can be integrated into any GNN working in an end-to-end fashion to boost its prediction accuracy. Our experiments showed that GRN remarkably boosted the prediction accuracy of four conventional GNN models across four neurological datasets.

[Neuroradiological and pathohistological markers of the main epileptogenic substrates in children.Cortical malformations]

High-resolution MRI is an important tool in the diagnosis of structural epilepsy in determining the seizure initiation zones, identification of the mechanisms of epileptogenesis in predicting outcomes and preventing postoperative complications in patients. In this article we demonstrate the neuroradiological and pathohistological characteristics of the main epileptogenic substrates in children using modern classification. The first part of the article is devoted to cortical malformations as the most common epileptogenic cerebral disorders.

Predicting the evolution trajectory of population-driven connectional brain templates using recurrent multigraph neural networks

The mapping of the time-dependent evolution of the human brain connectivity using longitudinal and multimodal neuroimaging datasets provides insights into the development of neurological disorders and the way they alter the brain morphology, structure and function over time. Recently, the connectional brain template (CBT) was introduced as a compact representation integrating a population of brain multigraphs, where two brain regions can have multiple connections, into a single graph. Given a population of brain multigraphs observed at a baseline timepoint t₁, we aim to learn how to predict the evolution of the population CBT at follow-up timepoints t>t₁. Such model will allow us to foresee the evolution of the connectivity patterns of healthy and disordered individuals at the population level. Here we present recurrent multigraph integrator network (ReMI-Net^⋆) to forecast population templates at consecutive timepoints from a given single timepoint. In particular, we unprecedentedly design a graph neural network architecture to model the changes in the brain multigraph and identify the biomarkers that differentiate between the typical and atypical populations. Addressing such issues is of paramount importance in diagnosing neurodegenerative disorders at early stages and promoting new clinical studies based on the pinned-down biomarker brain regions or connectivities. In this paper, we demonstrate the design and use of the ReMI-Net^⋆ model, which learns both the multigraph node level and time level dependencies concurrently. Thanks to its novel graph convolutional design and normalization layers, ReMI-Net^⋆ predicts well-centered, discriminative, and topologically sound connectional templates over time. Additionally, the results show that our model outperforms all benchmarks and state-of-the-art methods by comparing and discovering the atypical connectivity alterations over time. Our ReMI-Net^⋆ code is available on GitHub at https://github.com/basiralab/ReMI-Net-Star.

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.

Mapping developmental regionalization and patterns of cortical surface area from 29 post-menstrual weeks to 2 years of age

Surface area of the human cerebral cortex expands extremely dynamically and regionally heterogeneously from the third trimester of pregnancy to 2 y of age, reflecting the spatial heterogeneity of the underlying microstructural and functional development of the cerebral cortex. However, little is known about the developmental patterns and regionalization of cortical surface area during this critical stage, due to the lack of high-quality imaging data and accurate computational tools for pediatric brain MRI data. To fill this critical knowledge gap, by leveraging 1,037 high-quality MRI scans with the age between 29 post-menstrual weeks and 24 mo from 735 pediatric subjects in two complementary datasets, i.e., the Baby Connectome Project (BCP) and the developing Human Connectome Project (dHCP), and state-of-the-art dedicated image-processing tools, we unprecedentedly parcellate the cerebral cortex into a set of distinct subdivisions purely according to the developmental patterns of the cortical surface. Our discovered developmentally distinct subdivisions correspond well to structurally and functionally meaningful regions and reveal spatially contiguous, hierarchical, and bilaterally symmetric patterns of early cortical surface expansion. We also show that high-order association subdivisions, where cortical folds emerge later during prenatal stages, undergo more dramatic cortical surface expansion during infancy, compared with the central regions, especially the sensorimotor and insula cortices, thus forming a distinct central-pole division in early cortical surface expansion. These results provide an important reference for exploring and understanding dynamic early brain development in health and disease.

Path Signature Neural Network of Cortical Features for Prediction of Infant Cognitive Scores

Studies have shown that there is a tight connection between cognition skills and brain morphology during infancy. Nonetheless, it is still a great challenge to predict individual cognitive scores using their brain morphological features, considering issues like the excessive feature dimension, small sample size and missing data. Due to the limited data, a compact but expressive feature set is desirable as it can reduce the dimension and avoid the potential overfitting issue. Therefore, we pioneer the path signature method to further explore the essential hidden dynamic patterns of longitudinal cortical features. To form a hierarchical and more informative temporal representation, in this work, a novel cortical feature based path signature neural network (CF-PSNet) is proposed with stacked differentiable temporal path signature layers for prediction of individual cognitive scores. By introducing the existence embedding in path generation, we can improve the robustness against the missing data. Benefiting from the global temporal receptive field of CF-PSNet, characteristics consisted in the existing data can be fully leveraged. Further, as there is no need for the whole brain to work for a certain cognitive ability, a top K selection module is used to select the most influential brain regions, decreasing the model size and the risk of overfitting. Extensive experiments are conducted on an in-house longitudinal infant dataset within 9 time points. By comparing with several recent algorithms, we illustrate the state-of-the-art performance of our CF-PSNet (i.e., root mean square error of 0.027 with the time latency of 518 milliseconds for each sample).

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 cascaded nested network for 3T brain MR image segmentation guided by 7T labeling

Accurate segmentation of the brain into gray matter, white matter, and cerebrospinal fluid using magnetic resonance (MR) imaging is critical for visualization and quantification of brain anatomy. Compared to 3T MR images, 7T MR images exhibit higher tissue contrast that is contributive to accurate tissue delineation for training segmentation models. In this paper, we propose a cascaded nested network (CaNes-Net) for segmentation of 3T brain MR images, trained by tissue labels delineated from the corresponding 7T images. We first train a nested network (Nes-Net) for a rough segmentation. The second Nes-Net uses tissue-specific geodesic distance maps as contextual information to refine the segmentation. This process is iterated to build CaNes-Net with a cascade of Nes-Net modules to gradually refine the segmentation. To alleviate the misalignment between 3T and corresponding 7T MR images, we incorporate a correlation coefficient map to allow well-aligned voxels to play a more important role in supervising the training process. We compared CaNes-Net with SPM and FSL tools, as well as four deep learning models on 18 adult subjects and the ADNI dataset. Our results indicate that CaNes-Net reduces segmentation errors caused by the misalignment and improves segmentation accuracy substantially over the competing methods.

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.

(Un)common space in infant neuroimaging studies: A systematic review of infant templates

In neuroimaging, spatial normalization is an important step that maps an individual's brain onto a template brain permitting downstream statistical analyses. Yet, in infant neuroimaging, there remain several technical challenges that have prevented the establishment of a standardized template for spatial normalization. Thus, many different approaches are used in the literature. To quantify the popularity and variability of these approaches in infant neuroimaging studies, we performed a systematic review of infant magnetic resonance imaging (MRI) studies from 2000 to 2020. Here, we present results from 834 studies meeting inclusion criteria. Studies were classified into (a) processing data in single subject space, (b) using an off the shelf, or "off the shelf," template, (c) creating a study specific template, or (d) using a hybrid of these methods. We found that across the studies in the systematic review, single subject space was the most used (no common space). This was the most used common space for diffusion-weighted imaging and structural MRI studies while functional MRI studies preferred off the shelf atlases. We found a pattern such that more recently published studies are more commonly using off the shelf atlases. When considering special populations, preterm studies most used single subject space while, when no special populations were being analyzed, an off the shelf template was most common. The most used off the shelf templates were the UNC Infant Atlases (24%). Using a systematic review of infant neuroimaging studies, we highlight a lack of an established "standard" template brain in these studies.

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.

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.

A brain extraction algorithm for infant T2 weighted magnetic resonance images based on fuzzy c-means thresholding

It is challenging to extract the brain region from T2-weighted magnetic resonance infant brain images because conventional brain segmentation algorithms are generally optimized for adult brain images, which have different spatial resolution, dynamic changes of imaging intensity, brain size and shape from infant brain images. In this study, we propose a brain extraction algorithm for infant T2-weighted images. The proposed method utilizes histogram partitioning to separate brain regions from the background image. Then, fuzzy c-means thresholding is performed to obtain a rough brain mask for each image slice, followed by refinement steps. For slices that contain eye regions, an additional eye removal algorithm is proposed to eliminate eyes from the brain mask. By using the proposed method, accurate masks for infant T2-weighted brain images can be generated. For validation, we applied the proposed algorithm and conventional methods to T2 infant images (0–24 months of age) acquired with 2D and 3D sequences at 3T MRI. The Dice coefficients and Precision scores, which were calculated as quantitative measures, showed the highest values for the proposed method as follows: For images acquired with a 2D imaging sequence, the average Dice coefficients were 0.9650 ± 0.006 for the proposed method, 0.9262 ± 0.006 for iBEAT, and 0.9490 ± 0.006 for BET. For the data acquired with a 3D imaging sequence, the average Dice coefficient was 0.9746 ± 0.008 for the proposed method, 0.9448 ± 0.004 for iBEAT, and 0.9622 ± 0.01 for BET. The average Precision was 0.9638 ± 0.009 and 0.9565 ± 0.016 for the proposed method, 0.8981 ± 0.01 and 0.8968 ± 0.008 for iBEAT, and 0.9346 ± 0.014 and 0.9282 ± 0.019 for BET for images acquired with 2D and 3D imaging sequences, respectively, demonstrating that the proposed method could be efficiently used for brain extraction in T2-weighted infant images.

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.

Reference-Relation Guided Autoencoder with Deep CCA Restriction for Awake-to-Sleep Brain Functional Connectome Prediction

The difficulty of acquiring resting-state fMRI of early developing children under the same condition leads to a dedicated protocol, i.e., scanning younger infants during sleep and older children during being awake, respectively. However, the obviously different brain activities of sleep and awake states arouse a new challenge of awake-to-sleep connectome prediction/translation, which remains unexplored despite its importance in the longitudinally-consistent delineation of brain functional development. Due to the data scarcity and huge differences between natural images and geometric data (e.g., brain connectome), existing methods tailored for image translation generally fail in predicting functional connectome from awake to sleep. To fill this critical gap, we unprecedentedly propose a novel reference-relation guided autoencoder with deep CCA restriction (R²AE-dCCA) for awake-to-sleep connectome prediction. Specifically, 1) A reference-autoencoder (RAE) is proposed to realize a guided generation from the source domain to the target domain. The limited paired data are thus greatly augmented by including the combinations of all the age-restricted neighboring subjects as the references, while the target-specific pattern is fully learned; 2) A relation network is then designed and embedded into RAE, which utilizes the similarity in the source domain to determine the belief-strength of the reference during prediction; 3) To ensure that the learned relation in the source domain can effectively guide the generation in the target domain, a deep CCA restriction is further employed to maintain the neighboring relation during translation; 4) New validation metrics dedicated for connectome prediction are also proposed. Experimental results showed that our proposed R²AE-dCCA produces better prediction accuracy and well maintains the modular structure of brain functional connectome in comparison with state-of-the-art methods.

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.