Abstract
Introduction:
Gyrification is the intricate process through which the mammalian cerebral cortex develops its characteristic pattern of sulci and gyri. Monitoring gyrification provides valuable insights into brain development and identifies potential abnormalities at an early stage. This study analyzes the cortical structure in neurotypical and pathological (spina bifida) fetuses using various shape descriptors to shed light on the gyrification process during pregnancy.
Methods:
We compare morphometric properties encoded by commonly used scalar point-wise curvature-based signatures—such as mean curvature (H), Gaussian curvature (K), shape index (SI), and curvedness (C)—with multidimensional point-wise shape signatures, including spectral geometry processing methods like the Heat Kernel Signature (HKS) and Wave Kernel Signature (WKS), as well as the Signature of Histograms of Orientations (SHOT), which combines histogram and signature techniques. These latter signatures originate from computer graphics techniques and are rarely applied in the medical field. We propose a novel technique to derive a global descriptor from a given point-wise signature, obtaining GHKS, GWKS, and GSHOT. The extracted signatures are then evaluated using Support Vector Regression (SVR)-based algorithms to predict fetal gestational age (GA).
Results:
GSHOT better encodes the GA to other global multidimensional point-wise shape signatures (GHKS, GWKS) and commonly used scalar point-wise curvature-based signatures (C, H, K, SI, FI), achieving a prediction R2 of 0.89 and a mean absolute error of 6 days in neurotypical fetuses, and a R2 of 0.64 and a mean absolute error of 10 days in pathological fetuses.
Conclusion:
GSHOT provides researchers with an advanced tool to capture more nuanced aspects of fetal brain development and, specifically, of the gyrification process.
1 Introduction
Magnetic resonance imaging (MRI) has become a pivotal tool for the study of brain morphology and understanding structural alterations associated with various pathological conditions. Various geometric quantities can be exploited to summarize the morphometric information, providing valuable insights into population-based studies, contributing to our understanding of brain-related disorders, and paving the way for more personalized approaches to diagnosis and treatment. Advanced shape analysis techniques allow us to explore new dimensions of brain morphology beyond traditional measures (e.g., brain volumes and surface areas). Some of these techniques utilize spectral () and local extrinsic geometric () properties to gain deeper insights into brain shape characteristics. Spectral shape analysis techniques involve analyzing the shape of an object based on its spectral properties, and this is typically accomplished by encoding the shape through a differential operator and computing its eigendecomposition. On the other hand, local geometric properties focus on analyzing the shape at a more localized and detailed level, typically involving specific quantities such as curvature, surface normals, or deformations at specific points or regions of the shape. In this fashion, shapes can be compared by measuring similarities between these features. Thus, the efficacy of shape descriptors can be assessed in terms of discriminativeness and robustness against shape variations due to noise or deformations ().
In this work, we analyze the fetal brain cortical structure using different shape descriptors to enhance our comprehension of the fetal brain gyrification process, i.e., the formation of gyri and sulci. The gyrification process plays a crucial role in brain development, contributing significantly to overall growth, organization, and functionality. Just like fetal ultrasound provides an estimate of gestational age (GA) by measuring basic morphometric features (such as skull size or femur length), tools linking brain morphological MR images to the central nervous system development can be a valuable resource for monitoring pregnancy and detecting fetal diseases in their initial stages. Here, we estimate the fetus GA in weeks, comparing the cortical structure morphometric properties encoded with the commonly used scalar point-wise curvature-based descriptors to those derived via multidimensional point-wise shape signatures which are widely used in computer graphics analysis. We examine several scalar point-wise signatures based on curvature: mean curvature (H), which measures extrinsic curvature or folding; Gaussian curvature (K), which measures intrinsic curvature or distortion; shape index (SI) and folding index (FI), indicators of shape and folding patterns (). Moreover, we consider the curvedness (C), a signature incorporating information from both H and K, a valuable measure of the gyrification process (, ). On the other hand, we examine three different multidimensional point-wise shape signatures that are rarely applied in the medical field: the Heat Kernel Signature [HKS, ()], the Wave Kernel Signature [WKS, ()], and the Signature of Histograms of OrienTations [SHOT, ()]. HKS is derived from the heat equation, a partial differential equation that describes the heat diffusion across the surface over time. Similarly, WKS is derived from the solution of another partial differential equation, the wave equation, which describes the evolution of waves across the surface over time. SHOT is computed by dividing the neighborhood around each point into multiple cells and calculating histograms of relative orientations of the normals in each cell. These histograms are then concatenated to obtain the final signature, which results in a compact representation of the local geometric properties of the shape. In the last years, several studies have been proposed to analyze fetal brain gyrification by extracting cortical surface morphometric properties. In this context, scalar point-wise curvature-based measures are the gold standard for assessment of neurodevelopment (, , –). Other novel techniques based on sulcal pattern analysis can be employed to observe geometric and topological patterning of early sulcal folds, including 3D positions, sulcal basin surface area, and depth (, –).
Furthermore, we define a novel procedure to extract a global encoding framework from these multidimensional point-wise shape signatures, leading to their global version. We namely refer to the global descriptors produced through our pipeline as Global plus the name of the input pointwise descriptor we use, Global Heat Kernel Signature (GHKS), Global Wave Kernel Signature (GWKS), and Global Signature of Histograms of OrienTations (GSHOT). The global descriptor has several excellent advantages such as it allows for shape comparisons using minimal shape preprocessing, it is robust to noise since it implicitly employs surface smoothing by neglecting higher frequencies of the shape, and finally, it encodes isometric invariance properties of the shape, which are crucial to deal with shape deformations.
We tested our descriptors in the context of the fetal brain gyrification process. A linear Support Vector Regression (SVR)-based approach () was employed to predict the fetus GA from the cortical structure morphometric properties encoded by descriptors. Experiments on a public dataset of 80 fetuses (n = 31 neurotypical and n = 49 pathological (), and two public atlases of 18 and 16 fetuses (, ) showed promising prediction results in distinguishing the fetal brain gyrification process.
2 Methods
The proposed approach comprises five main steps: data gathering, cortical structure reconstruction, computation of shape descriptors, and GA prediction.
2.1 Data
We included data from different sources in this study. In particular, we used two publicly available fetal brain atlases (, ) and one publicly available fetal brain dataset (). For each source, we used all the provided data without assessing the quality of the fetal brain high-resolution reconstruction and tissue segmentation. A brain atlas is a digital representation of the human brain population, which highlights common structural features and provides a reference point for researchers and clinicians to compare and analyze specific brain regions. On the other hand, a brain dataset refers to a collection of brain images of real fetuses, thus characterized by unique variations.
The fetal brain atlas introduced by Gholipour et al. () (hereafter, “CRL atlas”) is defined at the GA range of 21–38 weeks. It consists of an age-specific T2-weighted (T2w) template and label images of 124 brain tissues, including gray matter (GM) and white matter (WM). The fetal brain atlas introduced by Uus et al. () (hereafter, “dHCP fetal atlas”) is defined at 21–36 weeks. It includes age-specific T2w templates and 19 brain tissue labels, separate for each hemisphere.
On the other hand, the fetal brain dataset introduced initially by Payette et al. (22) and later updated () (hereafter, “FeTA dataset”) consists of MRI-reconstructed images of 80 fetuses (n = 49 pathological and n = 31 neurotypical) defined in the GA range of 20–35 weeks. Each subject was released with a T2w template brain reconstruction (reconstructed with either NiftyMIC1, MIALSRTK2, or Simple IRTK3) with the corresponding seven brain tissue label images. Pathological subjects included fetuses with spina bifida either before or after fetal spinal lesion repair surgery, as these were the only publicly available pathological datasets (23). A summary of the available cohort of fetuses is reported in Table 1.
Table 1
| Data | MRI contrast | Tissue labels | GA range | Cohort | Public link |
|---|---|---|---|---|---|
| CRL atlas | T2w | 124 | 21–38 weeks | 18 | http://crl.med.harvard.edu/research/fetal_brain_atlas |
| dHCP fetal atlas | T2w | 19 | 21–36 weeks | 16 | https://gin.g-node.org/kcl_cdb/fetal_brain_mri_atlas |
| FeTA dataset | T2w | 7 | 20–35 weeks | 80 (31 neurotypical 49 pathological) | https://www.synapse.org/#!Synapse:syn25649159/wiki/610007 |
Summary of the publicly available fetal brain MRI atlases and dataset used in our study.
2.2 Cortical structure reconstruction
Our study focused on the cortical structure of the fetal brain, which is defined as the external layer of the parenchymal tissue and will become the cortical GM in the mature brain (Figure 1). During embryonic development, the brain is surrounded by a thin layer (darker than other tissues in the T2w images) called the cortical plate (CP). In the beginning, the CP is a flat and smoothed structure; as the brain grows and enlarges, it thickens and differentiates into different cortical layers. Visually, this results in folds, or gyri, and grooves, or sulci, that give the brain its characteristic wrinkled appearance. By the end of fetal development, the differentiated CP becomes the outermost layer of the brain, known as the cortical GM.
Figure 1
The CP folding process can be monitored using its boundary surfaces, i.e., the external and internal surfaces. The external surface separates the parenchyma from the cerebrospinal fluid. On the other hand, the inner surface divides the CP and the WM structure. We decided to focus on the inner cortical surface since the interface between WM and cortex is more stable and less prone to segmentation errors due to partial volume effects than the cortex-cerebrospinal fluid interface (
For each data set, we generated the inner cortical volume by merging the already validated tissue segmentation labels, which encompass the WM to the inner structures of the brain. Although the different datasets were generated using different segmentation protocols (
Figure 2

Inner cortical surface of the fetal brain from 28 to 38 gestational weeks. The surfaces depicted in the figure are generated from the previously quoted Gholipour et al. (
2.3 Shape descriptors
To capture the most informative intrinsic geometric properties of the inner cortical surface shape, we computed both scalar point-wise curvature-based signatures (C, H, K, SI, FI) and multidimensional point-wise shape signatures (HKS, WKS, SHOT).
Scalar point-wise curvature-based signatures are computed for each vertex of the surface mesh using the FreeSurfer function mris_curvature_stats (
Multidimensional point-wise shape signatures are computed accordingly. In detail, the HKS and WKS descriptors are implemented with an in-house MATLAB code. Here, we used k = 100 eigenvalues and scaled the temporal domain logarithmically in n = 10 time values, as suggested by Sun et al. (
Figure 3

An example of a global descriptor construction for the HKS signature (GHKS). (a) Each point of the inner cortical shape of the brain is colored according to the heat kernel (HK) value at time ti. (b) These values are then gathered into histograms for each scale ti. (c) The histograms are concatenated, leading to the global signature. The brain's surfaces shown in the figure are generated from a 33-week fetus of the previously quoted Gholipour et al. (
To aid in the data visualization of the signatures extracted by each descriptor, we performed a Principal Component Analysis (PCA). We derived a 2D scatter plot representing each descriptor's first and second principal components. These components contain the most relevant variations shown in the dataset, proving its capability to encode the changes in shape.
2.4 Gestational age prediction
We implemented a ML experiment to investigate whether and which derived signatures include the information associated with CP development, and if they can be used to predict a subject's GA.
We used a SVR algorithm to predict the GA in weeks for each image sample from the features extracted with the shape descriptor from the inner cortical surface mesh. SVR is one of the most powerful supervised machine learning approaches used for regression tasks (25), which aims to find a hyperplane maximizing the margin while minimizing errors in a high-dimensional feature space. It is an extension of the support vector machine classification algorithm but it predicts continuous output values instead of class labels. In the present study, only linear kernels were employed since nonlinear methods may require sample sizes that are too large to generalize well (26). Furthermore, we decided to use the z-score method to normalize the data, estimating the mean and the standard deviation on the training set and applying them to normalize the test set. We trained each SVR model on the signatures extracted from the atlases (
3 Results
The proposed GA prediction method was employed for the characterization of the fetal brain gyrification process occuring during pregnancy. A dataset of 114 images (n = 31 neurotypical and n = 49 pathological from public dataset, and n = 34 neurotypical from online available atlases) has been evaluated. After the surface construction, the scalar point-wise curvature-based signatures (C, H, K, SI, FI) and the global multidimensional point-wise shape signatures (GHKS, GWKS, and GSHOT) were computed and normalized by the z-score technique. The GA prediction procedure is employed as described in Section 2.4. Notably, each linear-based kernel SVR algorithm was trained on the shape signatures extracted from the inner cortical surface of fetuses included in atlases (
Table 2
| Performance metric [weeks] | Scalar point-wise curvature-based signatures | Global multidimensional point-wise shape signatures | ||||||
|---|---|---|---|---|---|---|---|---|
| C | H | K | SI | FI | GHKS | GWKS | GSHOT | |
| RMSE | 2.81 | 4.81 | 4.11 | 3.87 | 1.76 | 2.04 | 2.36 | 1.18 |
| MAE | 2.40 | 4.32 | 3.52 | 3.53 | 1.41 | 1.63 | 1.65 | 0.91 |
| R2 | 0.36 | −0.86 | −0.36 | −0.20 | 0.75 | 0.67 | 0.55 | 0.89 |
Gestational age prediction in neurotypical fetuses, expressed in weeks.
The goodness-of-fit of the individual linear-SVR models is evaluated by measuring the mean absolute error (MAE), the root mean square error (RMSE), and the coefficient of determination (R2). The scalar point-wise curvature-based signatures (C, H, K, SI, FI) and the global multidimensional point-wise shape signatures (GHKS, GWKS, GSHOT) are compared.
The best-performing metrics are shown in bold.
Table 3 shows the GA prediction performance of the considered individual SVR models tested on the pathological subset of the FeTA dataset. The results highlight a larger prediction error compared to the neurotypical subset of the FeTA dataset. Similarly, GSHOT outperforms other descriptors, achieving a prediction R2 of 0.64 and a corresponding MAE of 10.1 days. However, the agreement between predictions and ground truths for this model was poor based on Lin's concordance correlation coefficient (ρc = 0.80, 95%CI = 0.70–0.87).
Table 3
| Performance metric [weeks] | Scalar point-wise curvature-based signatures | Global multidimensional point-wise shape signatures | ||||||
|---|---|---|---|---|---|---|---|---|
| C | H | K | SI | FI | GHKS | GWKS | GSHOT | |
| RMSE | 2.96 | 8.50 | 6.25 | 6.04 | 3.27 | 3.31 | 3.53 | 1.87 |
| MAE | 2.62 | 8.10 | 5.75 | 5.55 | 2.47 | 2.84 | 2.71 | 1.44 |
| R2 | 0.10 | −6.46 | −3.03 | −2.77 | −0.10 | −0.13 | −0.28 | 0.64 |
Gestational age prediction in pathological (spina bifida) fetuses, expressed in weeks.
The goodness-of-fit of the individual linear-SVR models is evaluated by measuring the mean absolute error (MAE), the root mean square error (RMSE), and the coefficient of determination (R2). The scalar point-wise curvature-based signatures (C, H, K, SI, FI) and the global multidimensional point-wise shape signatures (GHKS, GWKS, GSHOT) are compared.
The best-performing metrics are shown in bold.
In Figure 4, the GA prediction obtained with the best global multidimensional point-wise shape signature (GSHOT) and with the best scalar point-wise curvature-based signature (FI) is visualized in the true vs. predicted response plot (top row), and the prediction model is evaluated using the residual plot (bottom row). All the GSHOT points are close to the diagonal line, suggesting an excellent estimation of the SVR model. Moreover, the points of FI are more dispersed than those of GSHOT. Finally, the residuals obtained in the GA prediction from the pathological population are much larger than those obtained from the neurotypical population. An exhaustive visualization of the results obtained from each descriptor is reported in the Supplementary Figures S1, S2, which shows larger prediction errors and a similar trend for both neurotypical and pathological populations.
Figure 4

A visualization of the results in GA prediction using GSHOT and FI. The figure displays the true vs. predicted response plots at the top and the residual plots at the bottom. Black points represent neurotypical fetuses, while pathological (spina bifida) fetuses are described by blue stars. The FeTA dataset was used for GA prediction (
PCA analysis shows that shape descriptors codify the largest part of the relevant information about the inner surface of the CP in the first few components. Among all the shape descriptors, GSHOT and FI provide the best performances. Their combined first two components explain 98.8% of the variability in the FeTA neurotypical fetuses and more than 95% in the FeTA pathological fetuses (Supplementary Table S1). Figure 5 shows that GSHOT provides a clear graphical representation of the fetal evolution. The different GA samples are distributed in a distinct “U” shape, with an increase from right to left on the first component (x-axis). The GA ranges behave symmetrically on the second component (y-axis), growing downwards on the negative axis values and upwards on the positive axis values. This geometric behavior interpretation cannot be inferred from FI results as it revealed lower discretization ability.
Figure 5

GSHOT and FI first and second principal components. Neurotypical fetuses are represented by points, while pathological (spina bifida) fetuses are represented by stars. The colorbar displays the color code used for the fetus GA ranges. These results are derived from Payette et al. (
The other descriptors have a graphical representation worse than FI (see Supplementary Figure S3).
4 Discussion
Gyrification in the human fetus occurs from the 10th gestational week and continues hierarchically almost up until the last weeks of pregnancy (29). During gyrification, the smooth surface of the fetal brain develops folds and wrinkles, increasing the surface area of the cerebral cortex. This folding is essential for accommodating the large number of cortical neurons and connections within the limited space of the skull. However, disruptions or abnormalities during this process can lead to cortical malformations such as lissencephaly [smooth brain, (30)], or polymicrogyria [excessive folding, (31)]. These malformations are associated with various neurological disorders and cognitive impairments. Therefore, understanding its construction mechanism is crucial for studying brain development, function, and disorders (32).
Here, we introduced novel multidimensional point-wise shape signatures (HKS, WKS, SHOT) to analyze the inner cortical surface development in the fetal brain by an innovative procedure. These signatures are a well-established method in computer graphics to analyze the geometry of an object based on its spectral properties. However, they have rarely been applied in the medical field due to the intricate nature of medical data, the heavy workload and professional expertise required, and the need for thorough validation and approval processes to comply with regulatory standards (
This study presents some limitations. First, other measures related to cortical folding, such as the gyrification index (34) and sulcal depth, can be studied. Furthermore, given cortical folding alterations associated with several pathologies (e.g., ventriculomegaly), cortical thickness is another measure worth investigating, considering partial volume effects. Second, other regression models (e.g., relevance vector regression RVR and Gaussian process regression GPR) can be tested by applying different kernel functions (e.g., polynomial, radial basis) to achieve the highest prediction performances. Third, we identified a larger prediction error in the pathological subset of the FeTA dataset used as test set. Unfortunately, we were not able to establish if this was caused by the clinical condition or by a model error. On the other hand, spina bifida is characterized by several malformative aspects both in the neonatal and in fetal central nervous system, encompassing not only altered gyrifications, but also altered brain and infratentorial structures size, and corpus callosum hypoplasia or partial dysgenesis (35–37). It is therefore reasonable to assume that the implemented shape descriptors highlighted a different developmental pattern in the gyrification process, leading to an error in the estimated GA that can be used as an indicator of pathology. Future works will investigate the brain's structure surface development across its different regions by using GSHOT to uncover new insight into the neurodevelopment process. Moreover, the quality of the T2w brain reconstructions and relative tissue label maps were not investigated as it is out of the scope of this study and has been previously addressed (
5 Conclusion
In this work, global multidimensional point-wise shape signatures (GHKS, GWKS, and GSHOT) are exploited to improve the prediction of GA in neurotypical and pathological fetuses.
GSHOT outperforms other global multidimensional point-wise signatures and scalar point-wise curvature-based signatures (C, H, K, SI, FI), providing researchers with a more sophisticated tool to capture more nuanced aspects of shapes. This approach enhances the accuracy and effectiveness of shape analysis tasks such as classification, segmentation, or matching, potentially leading to new methods for early detection of fetal diseases. In addition, a novel exploration of the fetal brain based on this approach can potentially uncover new insight into the structures development of the brain.
Finally, this innovative procedure for extracting multidimensional global descriptors from a given point-wise signature can also be applied in different scenarios of shape analysis within computer graphics.
Statements
Data availability statement
Publicly available datasets were analyzed in this study. This data can be found here: http://crl.med.harvard.edu/research/fetal_brain_atlas; https://gin.g-node.org/kcl_cdb/fetal_brain_mri_atlas; https://www.synapse.org/#!Synapse:syn25649159/wiki/610007.
Author contributions
TC: Data curation, Formal Analysis, Methodology, Resources, Software, Validation, Visualization, Writing – original draft, Writing – review & editing. LS: Resources, Writing – review & editing. AB: Supervision, Writing – review & editing. PB: Funding acquisition, Project administration, Writing – review & editing. SM: Conceptualization, Funding acquisition, Resources, Supervision, Writing – review & editing. DP: Conceptualization, Funding acquisition, Resources, Supervision, Writing – review & editing.
Funding
The author(s) declare financial support was received for the research, authorship, and/or publication of this article. This work was partially supported by grants from the Italian Ministry of Health to Paolo Brambilla (RF-2019-12371349, “Ricerca corrente 2023” and “Ricerca corrente 2024”) and Denis Peruzzo (“Ricerca corrente 2023” and “Ricerca corrente 2024”), and by a grant from the Italian Ministry of University and Research to Simone Melzi (“Dipartimenti di Eccellenza 2023-2027”, Department of Informatics, Systems and Communication of the University of Milano-Bicocca).
Acknowledgments
We gratefully acknowledge the support of NVIDIA Corporation with the RTX A5000 GPUs granted to the University of Milano-Bicocca through the Academic Hardware Grant Program for the project “Learned representations for implicit binary operations on real-world 2D-3D data”.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fped.2024.1471080/full#supplementary-material
Footnotes
1.^https://github.com/gift-surg/NiftyMIC
2.^https://github.com/Medical-Image-Analysis-Laboratory/mialsuperresolutiontoolkit
References
1.
LévyBZhangH. Spectral mesh processing. ACM SIGGRAPH 2010 Courses (2010). p. 1–312. 10.1145/1837101.1837109
2.
SaltiSTombariFDi StefanoL. SHOT: unique signatures of histograms for surface and texture description. Comput Vis Image Underst. (2014) 125:251–64. 10.1016/j.cviu.2014.04.011
3.
CastellaniUMirtuonoPMurinoVBellaniMRambaldelliGTansellaMet alA new shape diffusion descriptor for brain classification. Med Image Comput Comput Assist Interv. (2011) 14(Pt 2):426–33. 10.1007/978-3-642-23629-7_52
4.
DemirciNHollandMA. Cortical thickness systematically varies with curvature and depth in healthy human brains. Hum Brain Mapp. (2022) 43(6):2064–84. 10.1002/hbm.25776
5.
HuHHChenHYHungCIGuoWYWuYT. Shape and curvedness analysis of brain morphology using human fetal magnetic resonance images in utero. Brain Struct Funct. (2013) 218(6):1451–62. 10.1007/s00429-012-0469-3
6.
WuJAwateSPLichtDJClouchouxCdu PlessisAJAvantsBBet alAssessment of MRI-based automated fetal cerebral cortical folding measures in prediction of gestational age in the third trimester. AJNR Am J Neuroradiol. (2015) 36(7):1369–74. 10.3174/ajnr.A4357
7.
SunJOvsjanikovMGuibasL. A concise and provably informative multi-scale signature based on heat diffusion. Comput Graph Forum. (2009) 28(5):1383–92. 10.1111/j.1467-8659.2009.01515.x
8.
AubryMSchlickeweiUCremersD. The wave kernel signature: a quantum mechanical approach to shape analysis. 2011 IEEE Iinternational Cconference on Ccomputer Vvision Wworkshops (ICCV Wworkshops). IEEE (2011). p. 1626–33.
9.
BatchelorPGCastellano SmithADHillDLHawkesDJCoxTCDeanAF. Measures of folding applied to the development of the human fetal brain. IEEE Trans Med Imaging. (2002) 21(8):953–65. 10.1109/TMI.2002.803108
10.
PienaarRFischlBCavinessVMakrisNGrantPE. A methodology for analyzing curvature in the developing brain from preterm to adult. Int J Imaging Syst Technol. (2008) 18(1):42–68. 10.1002/ima.20138
11.
ClouchouxCKudelskiDGholipourAWarfieldSKViseurSBouyssi-KobarMet alQuantitative in vivo MRI measurement of cortical development in the fetus. Brain Struct Funct. (2012) 217(1):127–39. 10.1007/s00429-011-0325-x
12.
HabasPAScottJARoostaARajagopalanVKimKRousseauFet alEarly folding patterns and asymmetries of the normal human brain detected from in utero MRI. Cerebral Cortex. (2012) 22(1):13–25. 10.1093/cercor/bhr053
13.
ShimonyJSSmyserCDWidemanGAlexopoulosDHillJHarwellJet alComparison of cortical folding measures for evaluation of developing human brain. NeuroImage. (2016) 125:780–90. 10.1016/j.neuroimage.2015.11.001
14.
TaruiTImKMadanNMadankumarRSkotkoBGSchwartzAet alQuantitative MRI analyses of regional brain growth in living fetuses with down syndrome. Cerebral Cortex. (2020) 30(1):382–90. 10.1093/cercor/bhz094
15.
BenkarimOMHahnerNPiellaGGratacosEGonzález BallesterMAEixarchEet alCortical folding alterations in fetuses with isolated non-severe ventriculomegaly. NeuroImage Clin. (2018) 18:103–14. 10.1016/j.nicl.2018.01.006
16.
TaruiTMadanNFarhatNKitanoRCeren TanritanirAGrahamGet alDisorganized patterns of sulcal position in fetal brains with agenesis of corpus Callosum. Cerebral Cortex. (2018) 28(9):3192–203. 10.1093/cercor/bhx191
17.
TaruiTMadanNGrahamGKitanoRAkiyamaSTakeokaEet alComprehensive quantitative analyses of fetal magnetic resonance imaging in isolated cerebral ventriculomegaly. NeuroImage Clin. (2023) 37:103357. 10.1016/j.nicl.2023.103357
18.
DruckerHBurgesCJKaufmanLSmolaAVapnikV. Support vector regression machines. Advances in neural information processing systems 9 (1996).
19.
PayetteKLiHBde DumastPLicandroRJiHSiddiqueeMMRet alFetal brain tissue annotation and segmentation challenge results. Med Image Anal. (2023) 88:102833. Advance online publication. 10.1016/j.media.2023.102833
20.
GholipourARollinsCKVelasco-AnnisCOuaalamAAkhondi-AslAAfacanOet alA normative spatiotemporal MRI atlas of the fetal brain for automatic segmentation and analysis of early brain growth. Sci Rep. (2017) 7(1):476. 10.1038/s41598-017-00525-w
21.
UusAKyriakopoulouVCordero GrandeLChristiaensDPietschMPriceAet alMulti-channel spatio-temporal MRI atlas of the normal fetal brain development from the developing human connectome project. G-Node. (2023). 10.12751/g-node.ysgsy1
22.
PayetteKde DumastPKebiriHEzhovIPaetzoldJCShitSet alAn automatic multi-tissue human fetal brain segmentation benchmark using the fetal tissue annotation dataset. Sci Data. (2021) 8(1):167. 10.1038/s41597-021-00946-3
23.
CiceriTCasartelliLMontanoFConteSSquarcinaLBertoldoAet alFetal brain MRI atlases and datasets: a review. NeuroImage. (2024) 292:120603. 10.1016/j.neuroimage.2024.120603
24.
Rodriguez-CarranzaCEMukherjeePVigneronDBarkovichJStudholmeC. A framework for in vivo quantification of regional brain folding in premature neonates. NeuroImage. (2008) 41(2):462–78. 10.1016/j.neuroimage.2008.01.008
25.
ZhouWYanZZhangL. A comparative study of 11 non-linear regression models highlighting autoencoder, DBN, and SVR, enhanced by SHAP importance analysis in soybean branching prediction. Sci Rep. (2024) 14(1):5905. 10.1038/s41598-024-55243-x
26.
RasmussenPMMadsenKHLundTEHansenLK. Visualization of nonlinear kernel models in neuroimaging by sensitivity maps. NeuroImage. (2011) 55(3):1120–31. 10.1016/j.neuroimage.2010.12.035
27.
LawrenceILinK. A concordance correlation coefficient to evaluate reproducibility. Biometrics. (1989) 45(1):255–68. 10.2307/2532051
28.
McBrideGB. (2005). NIWA client report: HAM2005-062, 45, 307-310. A proposal for strength-of-agreement criteria for Lin’s concordance correlation coefficient.
29.
RajagopalanVScottJHabasPAKimKCorbett-DetigJRousseauFet alLocal tissue growth patterns underlying normal fetal human brain gyrification quantified in utero. J Neurosci. (2011) 31(8):2878–87. 10.1523/JNEUROSCI.5458-10.2011
30.
MochidaGH. Genetics and biology of microcephaly and lissencephaly. Semin Pediatr Neurol. (2009) 16(3):120–6. 10.1016/j.spen.2009.07.001
31.
SquierWJansenA. Polymicrogyria: pathology, fetal origins and mechanisms. Acta Neuropathol Commun. (2014) 2:80. 10.1186/s40478-014-0080-3
32.
FernándezVLlinares-BenaderoCBorrellV. Cerebral cortex expansion and folding: what have we learned?EMBO J. (2016) 35(10):1021–44. 10.15252/embj.201593701
33.
DahdouhSLimperopoulosC. Unsupervised fetal cortical surface parcellation. Proc SPIE Int Soc Opt Eng. (2016) 9784:97840J. 10.1117/12.2212805
34.
ZillesKArmstrongESchleicherAKretschmannHJ. The human pattern of gyrification in the cerebral cortex. Anat Embryol. (1988) 179(2):173–9. 10.1007/BF00304699
35.
MuftiNAertsenMEbnerMFidonLPatelPRahmanMBAet alCortical spectral matching and shape and volume analysis of the fetal brain pre- and post-fetal surgery for spina bifida: a retrospective study. Neuroradiology. (2021) 63(10):1721–34. 10.1007/s00234-021-02725-8
36.
MuftiNSaccoAAertsenMUshakovFOurselinSThomsonDet alWhat brain abnormalities can magnetic resonance imaging detect in foetal and early neonatal spina bifida: a systematic review. Neuroradiology. (2022) 64(2):233–45. 10.1007/s00234-021-02853-1
37.
MuftiNChappellJAertsenMEbnerMFidonLDeprestJet alAssessment of longitudinal brain development using super-resolution magnetic resonance imaging following fetal surgery for open spina bifida. Ultrasound Obstet Gynecol. (2023) 62(5):707–20. 10.1002/uog.26244
Summary
Keywords
fetal brain, gestational age prediction, shape descriptors, cortical surface, MRI
Citation
Ciceri T, Squarcina L, Bertoldo A, Brambilla P, Melzi S and Peruzzo D (2024) Fetal gestational age prediction via shape descriptors of cortical development. Front. Pediatr. 12:1471080. doi: 10.3389/fped.2024.1471080
Received
26 July 2024
Accepted
30 October 2024
Published
20 November 2024
Volume
12 - 2024
Edited by
Sahar Ahmad, University of North Carolina at Chapel Hill, United States
Reviewed by
Alena Uus, King’s College London, United Kingdom
Samson Nivins, Karolinska Institutet (KI), Sweden
Updates

Check for updates
Copyright
© 2024 Ciceri, Squarcina, Bertoldo, Brambilla, Melzi and Peruzzo.
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Paolo Brambilla paolo.brambilla1@unimi.it
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.