Anatomically optimised structural connectome derivation for the neonatal brain

Structural connectomics provides a means of characterising white matter (WM) network organisation using diffusion MRI (dMRI), but reconstructing neonatal structural connectomes remains challenging due to the immature and heterogeneous nature of the developing brain. Although recent neonatal studies have increasingly adopted approaches beyond conventional diffusion tensor imaging (DTI), there remains scope to incorporate more advanced tractography and streamline post-processing approaches within neonatal pipelines. Here, we propose a neonatal structural connectome pipeline integrating surface-based anatomical constraints, ensemble tractography, hierarchical anatomical refinement of pathways (HARP) filtering, and spherical-deconvolution informed filtering (SIFT2) reweighting. We acquired 3T MRI data from 300 neonates (185 preterm, 115 term-born) to assess the impact of pipeline configuration on network organisation. Mean gestational age (GA) at birth was 29.31 (23.43-32.86) weeks in preterm infants and 39.52 (37.00-42.43) weeks in term infants. At scan, mean GA was 40.85 (36.29-49.00) weeks in preterm infants and 42.06 (38.29-48.00) weeks in term infants. Four connectome variants were generated: raw surfaced-based tractograms, HARP-filtered, SIFT2-weighted, and combined HARP+SIFT2. Pipeline performance was primarily evaluated using an existing brain age prediction model: connectome-based predictive model (CPM) and our derived connectomes to estimate GA at scan in term and preterm infants. In addition, global graph theoretical metric analyses were conducted to provide a complementary assessment of pipeline configuration differences in network organisation. In total, seven graph theoretical metrics were calculated and evaluated: global efficiency, local efficiency, transitivity, modularity, assortativity, degree variance and statistical complexity. Across all configurations, connectomes supported accurate brain age prediction, with mean absolute errors (MAE) of 0.76-0.79 weeks in term infants and 1.00-1.06 weeks in preterm infants. Prediction error and brain age gap variability were greater in preterm than term infants across configurations, with HARP+SIFT2 achieving the lowest preterm MAE (1.00 weeks) and brain age gap (BAG) variability (SD = 1.38 weeks). However, differences in prediction performance across configurations were minimal, suggesting limited influence of pipeline choice on brain age prediction. However, analysis of global network metrics demonstrated that pipeline choice significantly influenced downstream measures of network organisation, with HARP+SIFT2 being the only configuration to identify significant preterm vs. term differences across all seven metrics examined. These findings demonstrate the importance of considering tractography and streamline-processing choices when constructing neonatal structural connectomes and provide a framework for more anatomically informed investigation of early brain development.