Large scale brain network dynamics in very preterm children and relationship with socio-emotional outcomes ========================================================================================================== * Vanessa Siffredi * Maria Chiara Liverani * Lorena G. A. Freitas * D. Tadros * Y. Farouj * Cristina Borradori Tolsa * Dimitri Van De Ville * Petra Susan Hüppi * Russia Hà-Vinh Leuchter ## Abstract Children born very preterm (VPT; < 32 completed weeks of gestation) are at high risk of neurodevelopmental and neurobehavioral difficulties associated with atypical brain maturation. The analysis of large-scale brain network dynamics during rest allows to investigate brain functional connectivity and their association with behavioural outcomes. Of interest, prematurity has been associated with atypical socio-emotional development with significant implications for the forming of peer relationships, adaptive functioning, academic achievement and mental health. In this study, we extracted dynamic functional connectivity by using the innovation-driven co-activation patterns (iCAPs) framework in preterm and full-term children aged 6 to 9 to explore changes in spatial organisation, laterality and temporal dynamics of spontaneous large-scale brain activity. Multivariate pattern analysis was used to explore potential biomarkers for socio-emotional difficulties in preterm children. Results show a spatial organisation of 13 networks retrieved comparable to full-term controls. Dynamic features and lateralisation of network brain activity were also comparable across groups for all large-scale brain networks. Despite apparent similarities in terms of dynamical functional connectivity parameters, multivariate pattern analysis unveiled group differences in their associations with socio-emotional abilities. While a pattern of decreased engagement in certain brain networks were associated with better socio-emotional abilities in full-term controls; in the VPT group, better socio-emotional abilities were associated with coordination of activity across different networks, i.e., coupling duration between different pairs of networks. It is possible that group differences reflect reduced degree of maturation of functional architecture in the VPT group for socio-emotional abilities. Keywords * preterm * dynamic functional connectivity * socio-emotional abilities * resting-state fMRI * children ## Introduction Preterm birth occurs during key phases of interrelated neurobiological processes underlying brain development. This seem particularly relevant in infants born very preterm (VPT; < 32 completed weeks of gestation). As a consequence, VPT children are at high risk of neurodevelopmental and neurobehavioral difficulties associated with atypical brain structural and functional maturation (Volpe, 2009). Magnetic resonance imaging (MRI) techniques provide powerful, non-invasive tools to delineate aberrant brain development related to prematurity. In particular, resting state-functional MRI (rs-fMRI) has emerged as a promising tool for studying neural networks underlying typical and atypical brain development (Ellard et al., 2021). Functional connectivity (FC) derived from rs-fMRI is defined as the temporal dependence of neuronal activity patterns of anatomically separated brain regions (Fox and Raichle, 2007; Van Den Heuvel and Pol, 2010). FC between regions in the absence of goal-directed activity and stimulation are used to identify networks with synchronous, spontaneous neuronal activity, termed resting state networks (Fox and Raichle, 2007; Van Den Heuvel and Pol, 2010). The functional relevance of resting state networks has been inferred not only from their spatial overlap with regions known to underpin sensorimotor and cognitive functions but also from their overlap with known white-matter anatomical pathways (Biswal et al., 1995; Fox et al., 2005; Greicius et al., 2009; Skudlarski et al., 2008). To date, most studies on rs-fMRI in VPT individuals of different ages from newborn to adolescent have used static FC, i.e., the correlation between the activation in different brain regions over the whole scanning time (Preti et al., 2017). Using a seed based correlation approach - in which few region of interest (ROI) time series are selected a priori and voxelwise cross-correlation is computed across the whole brain (Fox et al., 2005) - studies in VPT individuals at different ages have found altered FC in various brain networks using ROIs such as superior temporal, amygdala, posteromedial and lateral parietal, prefrontal or lateral and superior sensorimotor regions (Carter et al., 2004; Johns et al., 2019; Papini et al., 2016; Wehrle et al., 2018; Wilke et al., 2014). Similarly, alterations of different resting state networks have been found in VPT individuals when using whole-brain data-driven approaches, including FC of the salience, default mode or frontal networks (Ball et al., 2016; Damaraju et al., 2010; Lordier et al., 2019; White et al., 2014). However, recent studies suggest that stationary measures might be too simplistic to capture the full extent of resting state activity as that they ignore the inherently dynamic nature of brain FC (Christoff et al., 2016; Karahanoğlu and Van De Ville, 2017; Preti et al., 2017). In this perspective, dynamic approaches have been developed with the potential to identify meaningful variations over time in FC between different brain regions. Among various methods to identify dynamic FC, the innovation-driven co-activation patterns (iCAPs) framework detects moments of significantly changing brain activity to extract large-scale brain networks and their dynamic properties (Karahanoğlu and Van De Ville, 2015; Karahanoğlu and Van De Ville, 2017; Zöller et al., 2018). Using a whole-brain data-driven approach, this framework not only allows to extract spatial and temporal characteristics of large scale brain networks but also temporal overlaps between these different networks. Using the iCAPs framework, recent studies conducted in children, adolescents and adults have found significant association between large scale brain dynamics and psychiatric symptoms including depression, anxiety and psychotic symptoms (Piguet et al., 2021; Zöller et al., 2019). In VPT children, studies show that approximately 25% of them experience behavioural problems (Arpi and Ferrari, 2013; Johnson and Marlow, 2011). Of particular interest, socio-emotional abilities have been found to have important implications for the forming of peer relationships, adaptive functioning, academic achievement and mental health (Izard, 2011; Montagna and Nosarti, 2016a; Reyes et al., 2021). Prematurity has been associated with atypical socio-emotional development as early as the first year of life and extend from difficulties in emotional information processing, social understanding, emotion regulation, socialising, peer relationship as well as internalising problems (Bhutta et al., 2002; Healy et al., 2013; Hille et al., 2001; Johnson et al., 2015; Johnson et al., 2010; Jones et al., 2013; Landry et al., 1990; Langerock et al., 2013; Montagna and Nosarti, 2016b; Treyvaud et al., 2013; Witt et al., 2014). Importantly, deficits in socio-emotional processing and regulation in early life are considered precursors of later psychiatric and mental health problem (Briggs-Gowan and Carter, 2008; Carter et al., 2004). In VPT individuals, these difficulties have indeed been found to be long-lasting, with consequences observed in social, occupational and family functioning through adolescence into adulthood (Mendonça et al., 2019; Montagna and Nosarti, 2016b; Saigal et al., 2016); as well as increased risk of developing psychiatric disorders as adults, including depression, bipolar affective disorder, anxiety disorder and schizophrenia (Nosarti et al., 2012; Räikkönen et al., 2008; Walshe et al., 2008). In line with these studies, alteration in structural brain architecture has been associated with atypical socio-emotional development in VPT children (Fischi-Gomez et al., 2015; Kanel et al., 2021). Consequently, understanding brain functional underpinning of socio-emotional abilities in young VPT children might therefore inform on potential biomarker of socio-emotional difficulties and vulnerability of future psychiatric disorders. In this study, we extracted dynamic FC by using the iCAPs framework in young preterm and full-term children aged 6 to 9 to explore changes in spatial organisation, laterality and temporal dynamics of spontaneous large-scale brain activity. Multivariate pattern analysis was used to explore potential biomarkers for socio-emotional difficulties in preterm children. ## Methods and Materials ### Participants 227 VPT children born before 32 gestational weeks between 01.01.2008 and 01.05.2013, in the Neonatal Unit at the Geneva University Hospital (Switzerland) and followed up at the Division of Child Development and Growth, were invited to participate in the “Via-à-vis interventional study” between January 2017 and July 2019. VPT children were excluded if they had an intelligence quotient below 70, sensory or physical disabilities (cerebral palsy, blindness, hearing loss), or an insufficient understanding of French. A total of 45 VPT participants aged between 6 to 9-year-old were enrolled. Moreover, 17 term-born controls aged between 6 to 9-year-old were recruited through the community. Described in details below, of the 62 enrolled participants, 9 participants were excluded as they did not complete both the brain MRI scan and the neuropsychological assessment (VPT, n=7; full-term controls, n=2); and 13 participants were excluded due to high level of motion artefacts in rs-fMRI sequence described in the previous section (VPT, n=10; full-term controls, n=3). The final sample included 40 participants between 6 and 9 years of age: 28 VPT and 12 full-term participants. This study was approved by the Swiss Ethics Committees on research involving humans, ID: 2015-00175. Written informed consent was obtained from the principal caregiver and from the participant. ### Neuropsychological measures The Kaufman Assessment Battery for Children – 2nd Edition (K-ABC-II; (Kaufman and Kaufman, 2013)) was used to evaluate the Fluid-Crystallized Index (FCI) as a measure of general intellectual functioning. For children younger than 7 years of age, the FCI is derived from a linear combination of 11 core subtests; and for older children, it is derived from 10 core subtests. The FCI have a mean of 100 and a standard deviation of 15. Participants’ socio-emotional outcomes was assessed using four different measures. First, subtests of the Developmental Neuropsychological Assessment - 2nd Edition (NEPSY-II (Korkman et al., 2007)) including: a) the Affect Recognition subtest giving a total score assessing facial emotional recognition; and b) the Theory of Mind subtest giving a total score measuring the ability to understand mental contents, such as belief, intention or deception. As the Theory of Mind subtest does not provide a standard score, raw scores for both measures of the NEPSY-II were regressed on age at testing and standardised residuals was used as a score, called affect recognition and theory of mind. Second, the Internalised Score of the Strength and Difficulties Questionnaire – parent version (SDQ) was used to assess emotional and peer problems in daily life (Goodman, 1997, 2001). It rates participant’s internalised difficulties over the previous 6 months. The Internalised Score of the SDQ is scored on a Likert scale and is the sum of the emotional and peer problems scales. As standardised scores are not available for this measure, raw scores were regressed on age at testing and standardised residuals was used as a score, called internalised problems. Higher internalised scores reflect increased internalised difficulties in daily life. Finally, the Emotional Control Scale of the Behaviour Rating Inventory of Executive Function, parent version (BRIEF (Gioia et al. (2000)) was used to measure the extent to which the child is able to mediate emotional responses in daily life. As standardised scores are available for this measure, standardised scores were used (mean = 50, SD = 10), called emotional control. Higher emotional control scores reflect increased difficulties in emotional control. ### Magnetic Resonance Imaging (MRI) acquisition and preprocessing #### Magnetic Resonance Imaging acquisition MRI data were acquired at the Campus Biotech in Geneva, Switzerland, using a Siemens 3T Magnetom Prisma scanner. All participants completed a simulated “mock” MRI session prior to their MRI scan. This preparation process was conducted by trained research staff and allowed participants to familiarise themselves with the scanner and the scanning process, eventually raising any concerns they might have had prior to the MRI scan. Furthermore, this process is known to facilitated acquisition of good quality MRI images in children (de Bie et al., 2010; Tamnes et al., 2018). Structural T1-weighted MP-RAGE (magnetization-prepared rapid gradient-echo) sequences was acquired using the following parameters: voxel size = 0.9 × 0.9 × 0.9 mm; repetition time (TR) = 2,300 ms; echo time (TE) = 2.32 ms; inversion time (TI) = 900 ms; flip angle (FA) = 8°; and field of view (Fov) = 240 mm. Resting state functional images were T2*-weighted with a multislice gradient-echo-planar imaging (EPI) sequence of 64 slices; voxel size = 2 × 2 × 2 mm; TR = 720 ms; TE = 33 ms; Fov = 208 mm. During the rs-fMRI sequence, children were asked to keep their eyes closed and engage in mind wandering. Finally, a fieldmap was acquired each time a participant entered the scanner, with TR = 627 ms; TE1 = 5.19 ms; TE2 = 7.65 ms; and FA = 60°. #### Resting-State functional MRI data preprocessing Our data were preprocessed using SPM12 (Wellcome Department of Imaging Neuroscience, UCL, UK) in MATLAB R2016a (The MathWorks, Inc., Natick, Massachusetts, United States). Rs-fMRI data were converted from the native DICOM to NIFTI format and the preprocessing pipeline described by Freitas, Liverani and colleagues was used (Freitas et al., 2021). Rs-fMRI images were spatially realigned and unwarped, respectively, to correct for motion artefacts and potential geometric distortions. The unwarping step allows to improve the co-registration between structural and functional images and to reduce the distortion variability across subjects during spatial normalisation to a common space (Hutton et al., 2002). Functional images were then coregistered to their corresponding structural images in subject space and smoothed with a Gaussian filter of full width at half maximum (FWHM) = 6 mm. Structural images were segmented with the SPM12 segmentation algorithm to automatically identify different tissue types within the images, i.e., grey matter, white matter, cerebrospinal fluid (Ashburner and Friston, 2005); and a study-specific template was generated using Diffeomorphic Anatomical Registration using Exponential Lie algebra (DARTEL (Ashburner, 2007)) that will be used in the Innovation-Driven Co-Activation Patterns (iCAPs) framework described below. Finally, the first five rs-fMRI images were excluded, and average white matter and cerebrospinal fluid signals were regressed out from the BOLD time series. #### Head motion For rs-fMRI data, the mean framewise displacement for each frame was computed to quantify the extent of head motion from volume to volume for each participant (Power et al., 2012; Power et al., 2014). Following Power and colleagues’ recommendations, volume censoring (“scrubbing”) for motion correction were applied to frames with a mean framewise displacement above 0.5 mm, as well as one frame before and two after. Moreover, if more than 20% of the frames had to be scrubbed, the participant was removed from further analyses. Based on these criteria, 10 VPT and 3 full-term control participants were excluded from further analyses. Innovation-Driven Co-Activation Patterns (iCAPs) and extraction of iCAPs activation measures #### Innovation-Driven Co-Activation Patterns (iCAPs) The innovation-driven-co-activation-pattern (iCAP) analysis is a novel state-of-the-art rs-fMRI analysis tool that allows to derive a set of whole-brain spatial patterns of regions whose activity simultaneously increases or decreases, thus characterised by similar functional dynamic behaviour (Karahanoğlu and Van De Ville, 2015). For a comprehensive explanation of the methodology and procedure, please refer to Karahanoglu and colleagues (2015) and to Zöller and colleague (2019). We tailored the openly available MATLAB code ([https://c4science.ch/source/iCAPs/](https://c4science.ch/source/iCAPs/)) MATLAB vR2016a (The MathWorks, Inc., Natick, MA) to apply the iCAPs framework in VPT and full-term participants. The overall routine is composed of 4 steps: * Total activation (TA): In this first step, we employed TA which applies a voxel-wise hemodynamically-informed deconvolution (Farouj et al., 2017; Karahanoğlu et al., 2013) to the fMRI timeseries in a way that promotes the rareness of activity transients and spatially coherent activations. TA provides three types of information: (1) activity-related signals that are denoised fMRI signals; (2) sustained, or block-type, activity-inducing signals that are deconvolved signals; (3) innovation signals that are the derivative of the activity-inducing signals and encode transient brain activity episodes as spikes. * Significant transients detection: *Innovation* signals are computed as the temporal derivative of the deconvolved signals. The obtained signals can be seen as a representation in terms of transients in neural activity, where large amplitude transients implicitly identify change-points. Significant transients were determined using a two-step thresholding procedure. A temporal threshold estimated from a surrogate distribution that keeps only transient larger than 95% or lower than 5%. Then, a spatial thresholding procedure was applied, in which a frame was considered significant if at least 5% of the grey matter voxels were active. The frames showing significant transients are called *innovation* frames, and allow to identify time-points when a given region in the brain undergoes an increase or decrease in activity. * Aggregation: The significant, i.e. *innovation*, frames were warped into MNI (Montreal Neurologic Institute) space via a study-specific DARTEL template previously created (see preprocessing). All frames were then aggregated for clustering. * Clustering: The retained *innovation* frames, underwent K-means clustering to obtain simultaneously transitioning brain patterns, i.e., the iCAPs. The optimum number of 13 clusters was determined by consensus clustering and following recommendation of previous studies (Karahanoğlu et al., 2013; Monti et al., 2003; Zöller et al., 2019) (Supplementary Figure S1 & S2). * Time courses extraction: Time courses were obtained for all iCAPs using spatiotemporal transient-informed regression (Zöller et al., 2018). #### Extraction of temporal properties For computation of temporal properties, iCAPs’ time series were recovered by backprojecting each iCAP into subjects’ activity-inducing signals; i.e., block-type activity representations recovered by TA. For each iCAP, we then computed two measures representing the temporal characterisation of iCAPs: 1) occurrence, i.e., the number of activation blocks, 2) total duration, i.e., the total duration of overall activation as percentage of the total nonmotion scanning time. To explore the dynamic interactions during resting state, coupling and anticoupling duration of each pair of iCAPs were calculated as time points of same-signed or oppositely signed coactivation measured as percentage of the total nonmotion scanning time or as Jaccard score; i.e., percent joint activation time of the two respective iCAPs. #### Extraction of laterality measure To explore the laterality of brain networks, the iCAPs maps were co-registered MNI symmetrical template, available at [http://www.bic.mni.mcgill.ca/ServicesAtlases/ICBM152NLin2009](http://www.bic.mni.mcgill.ca/ServicesAtlases/ICBM152NLin2009). Based on previous studies, the laterality of activity maps aimed at exploring possible asymmetry effect between the two hemispheres by comparing lateralised amplitude maps of the iCAPs (Karolis et al., 2019). Therefore, the amplitudes of these patterns reflect the mean activity amplitude for each voxel when contributing to a certain network. In order to obtain an Amplitude Laterality Index (A-LI) for each voxel, we flipped the left hemisphere maps and subtracted them from unflipped right hemisphere maps (Siffredi et al., 2021). Positive and negative values in these A-LI maps reflect, respectively, right and left lateralisation. These maps were then averaged for each iCAPs in order to obtain one A-LI for each iCAPs and each participant. ### Statistical analyses #### Group comparisons of iCAPs activation measures Measures of occurrence and total duration of each iCAPs as well as coupling and anticoupling between each pair of iCAPs were compared between the VPT and full-term groups two-sample t-tests. Similarly, A-LI laterality measure was compared between the VPT and full-term groups using two-sample t-tests. The p values were corrected for multiple comparisons with the false discovery rate (FDR, (Benjamini and Hochberg, 1995)). Analyses were performed using the R software version 4.0.3, and R studio version 1.3.1093 (Team, 2020; Team, 2013). #### Multivariate correlation between iCAPs’ time courses and socio-emotional measures To evaluate multivariate patterns of correlation between iCAPs temporal characteristic and socio-emotional measures, we used partial least squares correlation (PLSC). A publicly available PLSC implementation in MATLAB was used: [https://github.com/danizoeller/myPLS](https://github.com/danizoeller/myPLS) (Kebets et al., 2019; Zoller et al., 2019). PLSC is a data-driven multivariate technique that maximizes the covariance between two matrices by identifying latent components which are linear combinations of the two matrices, i.e., socio-emotional and iCAPs temporal characteristics measures (McIntosh and Lobaugh, 2004). For the socio-emotional scores, the 4 scores of affect recognition, theory of mind, internalised problems and emotional control were considered. For affect recognition, theory of mind and internalised problems, raw scores were regressed on age at testing and socio-economic status (i.e., Largo score) confounds, and the standardised residuals were used for each score. For emotional control, standard scores were regressed on socio-economic status (i.e., Largo score) and the standardised residuals were used. Socio-emotional scores were stored in a 40 × 4 matrix denoted X, in which each row represents one participant and the matrix’s 4 columns consist of the 4 socio-emotional scores considered. For each iCAPs temporal characteristics, PLSC was repeated four times, each time considering only one temporal characteristics, i.e., occurrence, total duration, coupling, anticoupling. First, occurrence measures were gathered in a 40 × 13 matrix denoted Y, with each row matching one participant and each column representing the occurrence of each of the 13 iCAPs. Second, total duration measures were gathered in a 40 × 13 matrix denoted Y, with each row matching one participant and each column representing the total duration of each of the 13 iCAPs. Third, coupling measures were stored in a 40 × 78 matrix denoted Y, in which each row represents one participant and the matrix’s 78 columns consists of the coupling duration for each pair of iCAPs. Finally, anticoupling measures were stored in a 40 × 78 matrix denoted Y, in which each row represents one participant and the matrix’s 78 columns consists of the anticoupling duration for each pair of iCAPs. A cross-covariance matrix was then computed between X (participants x socio-emotional scores) and Y (participants x iCAPs temporal characteristics). Singular value decomposition was then applied to this cross-covariance matrix, resulting in latent components. Statistical significance of multivariate correlation patterns, i.e., latent component, was assessed with permutation testing (1000 permutations) and considered robust at p<0.01 following guidelines of previous studies (McIntosh and Lobaugh, 2004; Zöller et al., 2019). Stability of saliences were estimated using bootstrapping (500 bootstrap samples with replacement). Bootstrap ratio z-scores for each socio-emotional measure and iCAPs temporal characteristics were obtained by dividing each socio-emotional and iCAPs temporal characteristics weight by its bootstrap-estimated standard deviation, and a p-value was obtained for each bootstrap ratio z-score. The contribution of socio-emotional and iCAPs temporal characteristics weights for a given latent component was considered robust at p < 0.01 (i.e., absolute bootstrap ratio z-scores above 3 or below −3). ## Results ### Participant characteristics Of the 62 participants enrolled, 9 participants were excluded as they did not complete both the brain MRI scan and the neuropsychological assessment (VPT, n=7; FT, n=2); and 13 participants were excluded due to high level of motion artefacts in rs-fMRI sequence described in the previous section (VPT, n=10; FT, n=3). The final sample included 28 VPT and 12 full-term participants between 6 and 9 years of age (Table 1). Baseline characteristics were comparable between VPT and full-term participants for sex, age at assessment and the Fluid-Crystallized Index. Socio-economic status, as measured by the Largo score (Largo et al., 1989), showed significant group difference, with lower socio-economic status (i.e., higher Largo score) in the VPT group compared to the full-term group. Socio-emotional outcomes were comparable between the VPT and full-term groups across all scores, i.e., affect recognition, theory of mind, internalised problems and emotional control scores. View this table: [Table 1.](http://medrxiv.org/content/early/2022/03/21/2022.03.18.22272505/T1) Table 1. Neonatal and demographic characteristics, as well as socio-emotional outcomes of the VPT and full-term participants ### Extracted innovation-driven co-activation patterns (iCAPs) The iCAPs framework was applied to rs-fMRI scans of both VPT and full-term participants. The 13 extracted spatial maps represented in each iCAP correspond to well-known resting state networks and were reminiscent of common task-related and cognitive networks typically observed in fMRI studies (Van Den Heuvel and Pol, 2010). They also correspond to iCAP networks identified in previous studies (Bommarito et al., 2021; Piguet et al., 2021; Siffredi et al., 2021; Zöller et al., 2019). Specifically, the obtained networks included sensory-related networks, i.e., sensori-motor/auditory, primary and secondary visual; higher level cognitive network, i.e., the fronto-parietal, frontal and fronto-temporal left and right networks; and default-mode network (DMN) decomposed into precuneus/posterior DMN and posterior DMN. The remaining iCAPs comprised anterior insula/amygdala and cerebellar network decomposed into anterior cerebellum/vermis, posterior cerebellum and cerebellum/visual networks (Figure 1 and Supplementary Table S1). ![Figure 1.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2022/03/21/2022.03.18.22272505/F1.medium.gif) [Figure 1.](http://medrxiv.org/content/early/2022/03/21/2022.03.18.22272505/F1) Figure 1. Spatial patterns of the 13 innovation-driven coactivation patterns (iCAPs) retrieved from all participants include: precuneus/posterior default mode network (DMN), sensori-motor/auditory (AUD), anterior cerebellum/vermis, fronto-parietal, posterior default mode network (pDMN), posterior cerebellum, fronto-temporal right, anterior insula/amygdala, fronto-temporal left, secondary visual, primary visual, frontal, cerebellum/visual. Locations denote displayed slices in Montreal Neurological Institute coordinates. ![Figure 1](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2022/03/21/2022.03.18.22272505/F2.medium.gif) [Figure 1](http://medrxiv.org/content/early/2022/03/21/2022.03.18.22272505/F2) Figure 1 Associations between socio-emotional abilities and occurrence of the 13 iCAPs in the very preterm (VPT) and full-term (FT) groups based on the PLSC analysis. The diverging graph shows bootstrap ratio z-scores (x-axis) for each socio-emotional measures (y-axis) for the VPT and the FT groups in blue and yellow, respectively. Socio-emotional measures with an absolute bootstrap ratio z-score ≥ 3 or ≤ −3 (i.e., indicated by a back dash-dotted line on the graph) yield a robust contribution to the component (i.e., indicated by a black star). Only robust bootstrap ratio z-scores (absolute bootstrap ratio z-scores above or equal to 3) for iCAPs’ occurrences are shown. ### Group comparison of temporal properties of networks We compared the VPT and the full-term groups for each temporal property considered, i.e., the occurrence and the total duration for each iCAPs, as well as the percent joint activation time of coupling and anticoupling for each pair of iCAPs. There was no significant group difference for any of the temporal measures after multiple comparison correction (Supplementary Tables S2, S3, S4, S5). ### Group comparison of laterality of networks We compared the VPT and the full-term groups for the A-LI measure. There was no significant group difference for the laterality measure after multiple comparison correction (Supplementary Table S6). ### Association between iCAPs’ dynamics and socio-emotional measures First, the PLSC analysis applied on socio-emotional scores and occurrence of the 13 iCAPs in the preterm and the full-term controls groups identified one statistically significant latent components: latent component 1 (p = 0.006). Latent component 1 revealed a pattern of significant association in the full-term group only (Figure 1). In full-term children, increased socio-emotional abilities (i.e., affect recognition and theory of mind) were associated with increased occurrence of iCAP3-anterior cerebellum/vermis, iCAP6-posterior cerebellum, iCAP9-fronto-temporal left, iCAP10-visual secondary. In the VPT group, there was no significant association between socio-emotional abilities and occurrences of the 13 iCAPs network. Second, the PLSC analysis applied on socio-emotional measures and the total duration of the 13 iCAPs in the VPT and the full-term groups show no significant latent component. Third, the PLSC analysis applied on socio-emotional measures and coupling duration of the 13 iCAPs in the VPT and the full-term groups identified one statistically significant latent components: latent component 1 (p = 0.002). Latent component 1 revealed a pattern of significant association in the VPT group only (Figure 2). In VPT children, increased socio-emotional abilities (i.e., decreased difficulties in emotional control and internalised problems) were associated with patterns of increase and decrease duration of coupling between the different iCAP networks. More specifically, in VPT, better socio-emotional abilities were associated with increased coupling duration between iCAPS 8 to 12 (anterior insula/amygdala to frontal network), iCAPS 9 to 10 (fronto-temporal left to secondary visual), iCAPS 4 to 6 (fronto-parietal to posterior cerebellum) and iCAPS 1 to 3 (precuneus, pDMN to anterior cerebellum/vermis); and with decreased coupling duration between iCAPS 5 to 8 (posterior DMN to anterior insula/amygdala), iCAPS 1 to 6 (precuneus, pDMN to posterior cerebellum), iCAPS 6 to 9 (posterior cerebellum to fronto-temporal left), iCAPS 2 to 9 (sensori-motor/auditory to fronto-temporal left), iCAPS 9 to 13 (fronto-temporal left to cerebellum/visual), iCAPS 11 to 13 (primary visual to cerebellum/visual), iCAPS 5 to 10 (posterior DMN to secondary visual), iCAPS 4 to 10 (fronto-parietal to secondary visual), iCAPS 6 to 13 (posterior cerebellum to cerebellum/visual), iCAPS 9 to 12 (fronto-temporal left to frontal), iCAPS 2 to 13 (sensori-motor/auditory to cerebellum/visual), iCAPS 2 to 3 (sensori-motor/auditory to 15 anterior cerebellum/vermis), iCAPS 8 to 10 (anterior insula/amygdala to secondary visual), iCAPS 2 to 12 (sensori-motor/auditory to frontal). ![Figure 2.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2022/03/21/2022.03.18.22272505/F3.medium.gif) [Figure 2.](http://medrxiv.org/content/early/2022/03/21/2022.03.18.22272505/F3) Figure 2. Associations between socio-emotional abilities and coupling of the 13 pairs of iCAPs in the very preterm (VPT) and full-term (FT) groups based on the PLSC analysis. The diverging graph shows bootstrap ratio z-scores (x-axis) for each socio-emotion measure (y-axis) for the VPT and the FT groups in blue and yellow, respectively. Socio-emotional measures with an absolute bootstrap ratio z-score ≥ 3 or ≤ −3 (i.e., indicated by a back dash-dotted line on the graph) yield a robust contribution to the component (i.e., indicated by a black star). Only robust bootstrap ratio z-scores (absolute bootstrap ratio z-scores ≥ 3 or ≤ −3) for iCAPs’ couplings are shown. Finally, the PLSC analysis applied on socio-emotional measures and anticoupling duration of the 13 iCAPs in the VPT and the full-term groups show no significant latent component. Original saliences as well as their bootstrap-estimated standard deviations and bootstrap ratio z-scores for the PLSC analyses showing a significant latent component are reported in Supplementary Table S8 and S9. ## Discussion Using a novel whole-brain data-driven approach tailored to disentangle transient activity of coordinated brain networks from spontaneous rs-fMRI measurements, we unravelled large-scale resting state dynamics in VPT and full-term children aged 6 to 9 years. In this cohort of well-functioning VPT children, spatial organisation of the 13 networks retrieved was comparable to full-term controls. Dynamic features and lateralisation of network brain activity were also comparable across groups for the 13 large-scale brain networks. Despite apparent similarities in terms of dynamical FC parameters, multivariate pattern analysis unveiled group differences in their associations with socio-emotional abilities. While a pattern of decreased engagement in certain brain networks were associated with better socio-emotional abilities in full-term controls; in the VPT group, better socio-emotional abilities were associated with coordination of activity across different networks, i.e., coupling duration between different pairs of networks. In this cohort of well-functioning VPT children, overall dynamic FC characteristics were comparable to the full-term control group. First, spatial organisation of spontaneous brain activity was similar across the two groups. Second, temporal dynamics of iCAPs brain network were comparable in terms of occurrence and duration of engagement in a given network as well as in terms of coupling and anticoupling between pairs of networks. Moreover, the lateralisation of the large-scale brain networks retrieved was also comparable across groups. This pattern of dynamic FC similarities between VPT and full-term children are not consistent with previous studies using static FC (Degnan et al., 2015; Mossad et al., 2021; Rowlands et al., 2016; Wehrle et al., 2018). To the best of our knowledge, there is currently no study in school-age preterm children exploring dynamic FC. Nevertheless, Stoecklein and colleagues (2020) recently used FC variability or moment-to-moment variations of BOLD signals that represent neural dynamic range or flexibility completed in preterm and full-term infants (Stoecklein et al., 2020). Whole-brain FC variability was highly similar between 50 preterm and 25 full-term infants at term equivalent age suggesting that prematurity does not influence significantly FC variability at term equivalent age. In the current study, children with an IQ below 70 as well as children with high rates of motion artefact in rs-fMRI time series were excluded. It is possible that only well-functioning children with good self-regulation abilities were included in the final analyses. A recent study from Bolton and colleagues (2020) specifically examined motion characteristics during rs-fMRI acquisition in healthy adults (Bolton et al., 2020). Results show a broad array of behavioural and clinical characteristics related to motion including, among other, self-regulation, anxiety and depression. During childhood, it is possible that the association between motion and behavioural characteristics might be even more exacerbate resulting in the exclusion of children with increased behavioural difficulties in daily life and therefore increased atypicality in functional brain organisation. The iCAPs framework has previously being found to provide relevant makers of psychiatric symptomatology in adults and children, such as depression, anxiety and psychotic symptoms (Piguet et al., 2021; Zöller et al., 2019). In the current study, despite comparable dynamic features of resting state networks between VPT and full-term controls, different patterns of association with socio-emotional outcomes were observed across the two groups. In full-term controls, increased socio-emotional abilities were associated with reduced occurrence of different resting state networks, including the anterior cerebellum/vermis, posterior cerebellum, fronto-temporal left, and visual secondary networks. At the opposite, in VPT children, increased socio-emotional abilities were associated with a large pattern of increased and decreased coupling between networks. This is in line with previous studies conducted in the VPT population, suggesting that cognitive abilities rely more on coupling between networks in comparison to full-term controls for whom cognitive abilities is more related to the involvement of a given network (Barnes-Davis et al., 2018; Bäuml et al., 2015; Degnan et al., 2015; Mossad et al., 2021). Altogether, the findings suggest that, while in the full-term control group socio-emotional abilities was associated with different degree of integration of functional brain networks, i.e., recruitment of intra-network FC; in the VPsT group, socio-emotional abilities were associated with different degree of segregation between networks, i.e. increase and decrease in recruitment of inter-network FC. According to previous studies, typical development of functional architecture appears to show ongoing changes with an increase in integration along with a decrease in segregation of functional networks as children get older (Dosenbach et al., 2010; Menon, 2013; Rohr et al., 2018). Structurally this increased segregation and loss of increased integration has been shown at different ages in preterm compare to full-term infants (de Almeida et al., 2021; Fischi-Gomez et al., 2016). It is therefore possible that the group difference observed in the pattern of association between dynamic FC features and socio-emotional abilities reflect reduced degree of maturation of functional architecture in the VPT group. While providing new insight into dynamic functional organisaton in children born VPT, the current study has a number of limitations which need to be considered. First, the sample size of the study is rather small and warrants validation in larger samples. The small sample size may have limited the identification of subtle group differences in FC dynamics. Moreover, study participants comprised only well-functioning VPT children. The exclusion of 13 participants due to significant head motion during rs-fMRI acquisition might have accentuate this point even more through the exclusion of children with increased behavioural difficulties. Taking these limitations into account, dynamic features of FC might represent a relevant neuroimaging markers and inform on potential mechanisms through which preterm birth leads to neurodevelopmental deficits. ## Conclusions Leveraging recent advances in analysis of dynamic features of FC, we explored precise moments of brain network activation and interaction in VPT compared to full-term children from 6 to 9 years of age. In the current study, brain networks were comparable across groups not only in terms of spatial organisation but also in terms of laterality and dynamic characteristics. Despite group similarities in dynamical FC parameters, multivariate pattern analysis revealed different pattern of association with socio-emotional abilities in the VPT and full-term groups, that might reflect reduced degree of maturation of functional architecture in the VPT group. ## Supporting information Supplementary Materials [[supplements/272505_file03.pdf]](pending:yes) ## Data Availability Deidentified individual participant data (including data dictionaries) will be made available, in addition to study protocols, the statistical analysis plan, and the informed consent form. The data will be made available upon publication to researchers who provide a methodologically sound proposal for use in achieving the goals of the approved proposal. Proposals should be submitted to Russia.HaVinhLeuchter{at}unige.ch. ## Funding This work was supported by the Swiss National Science Foundation, No. 324730_163084 [PI: P.S. Hüppi]. ## Acknowledgements We thank and acknowledge all participating young adolescents and families who made this research possible. We also thank the Fondation Campus Biotech Geneva (FCBG), a foundation of the Swiss Federal Institute of Technology Lausanne (EPFL), the University of Geneva (UniGe), and the University Hospitals of Geneva (HUG); the Research Platform of the University Hospitals of Geneva (HUG) for their practical help. * Received March 18, 2022. * Revision received March 18, 2022. * Accepted March 21, 2022. * © 2022, Posted by Cold Spring Harbor Laboratory The copyright holder for this pre-print is the author. All rights reserved. The material may not be redistributed, re-used or adapted without the author's permission. ## References 1. Arpi, E., Ferrari, F., 2013. Preterm birth and behaviour problems in infants and preschool-age children: A review of the recent literature. Developmental Medicine sand Child Neurology 55, 788–796. 2. Ashburner, J., 2007. A fast diffeomorphic image registration algorithm. Neuroimage 38, 95–113. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.neuroimage.2007.07.007&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=17761438&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F03%2F21%2F2022.03.18.22272505.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000249992900011&link_type=ISI) 3. Ashburner, J., Friston, K.J., 2005. Unified segmentation. Neuroimage 26, 839–851. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.neuroimage.2005.02.018&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=15955494&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F03%2F21%2F2022.03.18.22272505.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000230211100020&link_type=ISI) 4. Ball, G., Aljabar, P., Arichi, T., Tusor, N., Cox, D., Merchant, N., Nongena, P., Hajnal, J.V., Edwards, A.D., Counsell, S., 2016. Machine-learning to characterise neonatal functional connectivity in the preterm brain. Neuroimage 124, 267–275. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.neuroimage.2015.08.055&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=26341027&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F03%2F21%2F2022.03.18.22272505.atom) 5. Barnes-Davis, M.E., Merhar, S.L., Holland, S.K., Kadis, D.S., 2018. Extremely preterm children exhibit increased interhemispheric connectivity for language: findings from fMRI-constrained MEG analysis. Developmental science 21, e12669. 6. Bäuml, J.G., Daamen, M., Meng, C., Neitzel, J., Scheef, L., Jaekel, J., Busch, B., Baumann, N., Bartmann, P., Wolke, D.J., 2015. Correspondence between aberrant intrinsic network connectivity and gray-matter volume in the ventral brain of preterm born adults. Cerebral cortex 25, 4135–4145. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/cercor/bhu133&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=24935776&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F03%2F21%2F2022.03.18.22272505.atom) 7. Benjamini, Y., Hochberg, Y., 1995. Controlling the false discovery rate: a practical and powerful approach to multiple testing. Journal of the Royal statistical society: series B 57, 289–300. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/S0166-4328(01)00297-2&link_type=DOI) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=A1995QE45300017&link_type=ISI) 8. Bhutta, A.T., Cleves, M.A., Casey, P.H., Cradock, M.M., Anand, K.J., 2002. Cognitive and behavioral outcomes of school-aged children who were born preterm: a meta-analysis. Jama 288, 728–737. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1001/jama.288.6.728&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=12169077&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F03%2F21%2F2022.03.18.22272505.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000177408100034&link_type=ISI) 9. Biswal, B., Zerrin Yetkin, F., Haughton, V.M., Hyde, J.S., 1995. Functional connectivity in the motor cortex of resting human brain using echo-planar MRI. Magnetic Resonance in Medicine 34, 537–541. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1002/mrm.1910340409&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=8524021&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F03%2F21%2F2022.03.18.22272505.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=A1995RX75100008&link_type=ISI) 10. Bolton, T.A., Kebets, V., Glerean, E., Zöller, D., Li, J., Yeo, B.T., Caballero-Gaudes, C., Van De Ville, D., 2020. Agito ergo sum: Correlates of spatio-temporal motion characteristics during fMRI. Neuroimage 209, 116433. 11. Bommarito, G., Tarun, A., Farouj, Y., Preti, M.G., Petracca, M., Droby, A., El Mendili, M.M., Inglese, M., Van De Ville, D., 2021. Altered anterior default mode network dynamics in progressive multiple sclerosis. Multiple Sclerosis Journal, 13524585211018116. 12. Briggs-Gowan, M.J., Carter, A.S., 2008. Social-emotional screening status in early childhood predicts elementary school outcomes. Pediatrics 121, 957–962. [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6MTA6InBlZGlhdHJpY3MiO3M6NToicmVzaWQiO3M6OToiMTIxLzUvOTU3IjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjIvMDMvMjEvMjAyMi4wMy4xOC4yMjI3MjUwNS5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 13. Carter, A.S., Briggs-Gowan, M.J., Davis, N.O., 2004. Assessment of young children’s social-emotional development and psychopathology: Recent advances and recommendations for practice. Journal of Child Psychology 45, 109–134. 14. Christoff, K., Irving, Z.C., Fox, K.C., Spreng, R.N., Andrews-Hanna, J.R., 2016. Mind-wandering as spontaneous thought: a dynamic framework. Nature reviews neuroscience 17, 718–731. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/nrn.2016.113&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=27654862&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F03%2F21%2F2022.03.18.22272505.atom) 15. Damaraju, E., Phillips, J., Lowe, J.R., Ohls, R., Calhoun, V.D., Caprihan, A., 2010. Resting-state functional connectivity differences in premature children. Frontiers in systems neuroscience 4, 23. 16. de Almeida, J.S., Meskaldji, D.-E., Loukas, S., Lordier, L., Gui, L., Lazeyras, F., Hüppi, P.S., 2021. Preterm birth leads to impaired rich-club organization and fronto-paralimbic/limbic structural connectivity in newborns. Neuroimage 225, 117440. 17. de Bie, H.M., Boersma, M., Wattjes, M.P., Adriaanse, S., Vermeulen, R.J., Oostrom, K.J., Huisman, J., Veltman, D.J., Delemarre-Van de Waal, H.A.J.E.j.o.p., 2010. Preparing children with a mock scanner training protocol results in high quality structural and functional MRI scans. 169, 1079–1085. 18. Degnan, A.J., Wisnowski, J.L., Choi, S., Ceschin, R., Bhushan, C., Leahy, R.M., Corby, P., Schmithorst, V.J., Panigrahy, A., 2015. Altered structural and functional connectivity in late preterm preadolescence: An anatomic seed-based study of resting state networks related to the posteromedial and lateral parietal cortex. PLoS One 10, e0130686. 19. Dosenbach, N.U., Nardos, B., Cohen, A.L., Fair, D.A., Power, J.D., Church, J.A., Nelson, S.M., Wig, G.S., Vogel, A.C., Lessov-Schlaggar, C.N., 2010. Prediction of individual brain maturity using fMRI. Science 329, 1358–1361. [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6Mzoic2NpIjtzOjU6InJlc2lkIjtzOjEzOiIzMjkvNTk5Ny8xMzU4IjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjIvMDMvMjEvMjAyMi4wMy4xOC4yMjI3MjUwNS5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 20. Ellard, K.K., Uribe, S., Funes, C.J., 2021. Resting-state functional magnetic resonance imaging (rsfMRI) in bipolar and unipolar depression. The Neuroscience of Depression. Elsevier, pp. 277–299. 21. Farouj, Y., Karahanoğlu, F.I., Van De Ville, D., 2017. Regularized spatiotemporal deconvolution of fMRI data using gray-matter constrained total variation. 2017 IEEE 14th International Symposium On Biomedical Imaging (Isbi 2017). Ieee, pp. 472–475. 22. Fischi-Gomez, E., Muñoz-Moreno, E., Vasung, L., Griffa, A., Borradori-Tolsa, C., Monnier, M., Lazeyras, F., Thiran, J.-P., Hüppi, P.S., 2016. Brain network characterization of high-risk preterm-born school-age children. NeuroImage: Clinical 11, 195–209. 23. Fischi-Gomez, E., Vasung, L., Meskaldji, D.E., Lazeyras, F., Borradori-Tolsa, C., Hagmann, P., Barisnikov, K., Thiran, J.P., Huppi, P.S., 2015. Structural Brain Connectivity in School-Age Preterm Infants Provides Evidence for Impaired Networks Relevant for Higher Order Cognitive Skills and Social Cognition. Cereb Cortex 25, 2793–2805. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/cercor/bhu073&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=24794920&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F03%2F21%2F2022.03.18.22272505.atom) 24. Fox, M.D., Raichle, M.E., 2007. Spontaneous fluctuations in brain activity observed with functional magnetic resonance imaging. Nature reviews neuroscience 8, 700–711. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/nrn2201&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=17704812&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F03%2F21%2F2022.03.18.22272505.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000248846200016&link_type=ISI) 25. Fox, M.D., Snyder, A.Z., Vincent, J.L., Corbetta, M., Van Essen, D.C., Raichle, M.E., 2005. The human brain is intrinsically organized into dynamic, anticorrelated functional networks. Proceedings of the National Academy of Sciences 102, 9673–9678. [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NDoicG5hcyI7czo1OiJyZXNpZCI7czoxMToiMTAyLzI3Lzk2NzMiO3M6NDoiYXRvbSI7czo1MDoiL21lZHJ4aXYvZWFybHkvMjAyMi8wMy8yMS8yMDIyLjAzLjE4LjIyMjcyNTA1LmF0b20iO31zOjg6ImZyYWdtZW50IjtzOjA6IiI7fQ==) 26. Freitas, L.G., Liverani, M.C., Siffredi, V., Schnider, A., Tolsa, C.B., Leuchter, R.H.-V., Van De Ville, D., Hüppi, P.S., 2021. Altered orbitofrontal activation in preterm-born young adolescents during performance of a reality filtering task. NeuroImage: Clinical 30, 102668. 27. Gioia, G., Isquith, P., Guy, S., Kenworthy, L., 2000. BRIEF – Behavior Rating Inventory of Executive Function. Professional manual. Psychological Assessment Resources Inc, 28. Odessa, FL. Goodman, R., 1997. The Strengths and Difficulties Questionnaire: a research note. J Child Psychol Psychiatry 38, 581–586. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1111/j.1469-7610.1997.tb01545.x&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=9255702&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F03%2F21%2F2022.03.18.22272505.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=A1997XN15800011&link_type=ISI) 29. Goodman, R., 2001. Psychometric properties of the strengths and difficulties questionnaire. J Am Acad Child Adolesc Psychiatry 40, 1337–1345. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1097/00004583-200111000-00015&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=11699809&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F03%2F21%2F2022.03.18.22272505.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000171752900015&link_type=ISI) 30. Greicius, M.D., Supekar, K., Menon, V., Dougherty, R.F., 2009. Resting-state functional connectivity reflects structural connectivity in the default mode network. Cerebral cortex 19, 72–78. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/cercor/bhn059&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=18403396&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F03%2F21%2F2022.03.18.22272505.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000261679400007&link_type=ISI) 31. Healy, E., Reichenberg, A., Nam, K.W., Allin, M.P., Walshe, M., Rifkin, L., Murray, R.M., Nosarti, C., 2013. Preterm birth and adolescent social functioning–alterations in emotion-processing brain areas. The Journal of pediatrics 163, 1596–1604. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.jpeds.2013.08.011&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=24070828&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F03%2F21%2F2022.03.18.22272505.atom) 32. Hille, E.T., den Ouden, A.L., Saigal, S., Wolke, D., Lambert, M., Whitaker, A., Pinto-Martin, J.A., Hoult, L., Meyer, R., Feldman, J.F., 2001. Behavioural problems in children who weigh 1000 g or less at birth in four countries. The Lancet 357, 1641–1643. 33. Hutton, C., Bork, A., Josephs, O., Deichmann, R., Ashburner, J., Turner, R., 2002. Image distortion correction in fMRI: a quantitative evaluation. Neuroimage 16, 217–240. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1006/nimg.2001.1054&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=11969330&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F03%2F21%2F2022.03.18.22272505.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000176624700019&link_type=ISI) 34. Izard, C.E., 2011. Forms and functions of emotions: Matters of emotion–cognition interactions. Emotion Review 3, 371–378. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1177/1754073911410737&link_type=DOI) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000306274900003&link_type=ISI) 35. Johns, C.B., Lacadie, C., Vohr, B., Ment, L.R., Scheinost, D., 2019. Amygdala functional connectivity is associated with social impairments in preterm born young adults. NeuroImage: Clinical 21, 101626. 36. Johnson, S., Gilmore, C., Gallimore, I., Jaekel, J., Wolke, D., 2015. The long-term consequences of preterm birth: what do teachers know? Developmental Medicine and Child Neurology 57, 571–577. 37. Johnson, S., Hollis, C., Kochhar, P., Hennessy, E., Wolke, D., Marlow, N., 2010. Psychiatric disorders in extremely preterm children: longitudinal finding at age 11 years in the EPICure study. Journal of the American Academy of Child Adolescent Psychiatry 49, 453–463. e451. 38. Johnson, S., Marlow, N., 2011. Preterm birth and childhood psychiatric disorders. Pediatric research 69, 11–18. 39. Jones, K.M., Champion, P.R., Woodward, L.J., 2013. Social competence of preschool children born very preterm. Early Hum Dev 89, 795–802. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.earlhumdev.2013.06.008&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F03%2F21%2F2022.03.18.22272505.atom) 40. Kanel, D., Vanes, L.D., Pecheva, D., Hadaya, L., Falconer, S., Counsell, S.J., Edwards, D.A., Nosarti, C., 2021. Neonatal White Matter Microstructure and Emotional Development during the Preschool Years in Children Who Were Born Very Preterm. eNeuro 8. 41. Karahanoğlu, F.I., Caballero-Gaudes, C., Lazeyras, F., Van De Ville, D., 2013. Total activation: fMRI deconvolution through spatio-temporal regularization. Neuroimage 73, 121–134. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.neuroimage.2013.01.067&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=23384519&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F03%2F21%2F2022.03.18.22272505.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000317084500012&link_type=ISI) 42. Karahanoğlu, F.I., Van De Ville, D., 2015. Transient brain activity disentangles fMRI resting-state dynamics in terms of spatially and temporally overlapping networks. Nature communications 6, 1–10. 43. Karahanoğlu, F.I., Van De Ville, D., 2017. Dynamics of large-scale fMRI networks: Deconstruct brain activity to build better models of brain function. Current Opinion in Biomedical Engineering 3, 28–36. 44. Karolis, V.R., Corbetta, M., Thiebaut de Schotten, M., 2019. The architecture of functional lateralisation and its relationship to callosal connectivity in the human brain. Nature communications 10, 1–9. 45. Kaufman, A.S., Kaufman, N.L., 2013. Kaufman assessment battery for children. Encyclopedia of Special Education: A Reference for the Education of Children, Adolescents, Adults with Disabilities Other Exceptional Individuals. 46. Kebets, V., Holmes, A.J., Orban, C., Tang, S., Li, J., Sun, N., Kong, R., Poldrack, R. A., Yeo, B.T.T., 2019. Somatosensory-Motor Dysconnectivity Spans Multiple Transdiagnostic Dimensions of Psychopathology. Biological Psychiatry . 47. Korkman, M., Kirk, U., Kemp, S., 2007. A developmental neuropsychological assessment 2nd Edition: NEPSY-II. Pearson, San Antonio, Texas. 48. Landry, S.H., Chapieski, M.L., Richardson, M.A., Palmer, J., Hall, S., 1990. The social competence of children born prematurely: Effects of medical complications and parent behaviors. Child Dev 61, 1605–1616. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.2307/1130768&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=2245750&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F03%2F21%2F2022.03.18.22272505.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=A1990ED70000026&link_type=ISI) 49. Langerock, N., van Hanswijck de Jonge, L., Graz, M.B., Hüppi, P.S., Tolsa, C.B., Barisnikov, K., 2013. Emotional reactivity at 12 months in very preterm infants born at< 29 weeks of gestation. Infant Behavior and Development 36, 289–297. 50. Largo, R.H., Pfister, D., Molinari, L., Kundu, S., Lipp, A., Duc, G., 1989. Significance of prenatal, perinatal and postnatal factors in the development of AGA preterm infants at five to seven years. Dev Med Child Neurol 31, 440–456. [PubMed](http://medrxiv.org/lookup/external-ref?access_num=2680687&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F03%2F21%2F2022.03.18.22272505.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=A1989AL55400004&link_type=ISI) 51. Lordier, L., Meskaldji, D.-E., Grouiller, F., Pittet, M.P., Vollenweider, A., Vasung, L., Borradori-Tolsa, C., Lazeyras, F., Grandjean, D., Van De Ville, D., Hüppi, P., 2019. Music in premature infants enhances high-level cognitive brain networks. Proceedings of the National Academy of Sciences 116, 12103–12108. [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NDoicG5hcyI7czo1OiJyZXNpZCI7czoxMjoiMTE2LzI0LzEyMTAzIjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjIvMDMvMjEvMjAyMi4wMy4xOC4yMjI3MjUwNS5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 52. McIntosh, A.R., Lobaugh, N.J.J.N., 2004. Partial least squares analysis of neuroimaging data: applications and advances. 23, S250–S263. 53. Mendonça, M., Bilgin, A., Wolke, D., 2019. Association of preterm birth and low birth weight with romantic partnership, sexual intercourse, and parenthood in adulthood: a systematic review and meta-analysis. JAMA Network Open 2, e196961–e196961. 54. Menon, V., 2013. Developmental pathways to functional brain networks: emerging principles. Trends in cognitive sciences 17, 627–640. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.tics.2013.09.015&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=24183779&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F03%2F21%2F2022.03.18.22272505.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000328719600008&link_type=ISI) 55. Montagna, A., Nosarti, C., 2016a. Socio-Emotional Development Following Very Preterm Birth: Pathways to Psychopathology. Front Psychol 7, 80. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.3389/fpsyg.2016.00080&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=26903895&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F03%2F21%2F2022.03.18.22272505.atom) 56. Montagna, A., Nosarti, C., 2016b. Socio-emotional development following very preterm birth: pathways to psychopathology. Frontiers in psychology 7, 80. 57. Monti, S., Tamayo, P., Mesirov, J., Golub, T., 2003. Consensus clustering: a resampling-based method for class discovery and visualization of gene expression microarray data. Machine learning 52, 91–118. 58. Mossad, S.I., Vandewouw, M.M., Smith, M.L., Taylor, M.J., 2021. The preterm social brain: Altered functional networks for Theory of Mind in very preterm children. Brain communications 3, fcaa237. 59. Nosarti, C., Reichenberg, A., Murray, R.M., Cnattingius, S., Lambe, M.P., Yin, L., MacCabe, J., Rifkin, L., Hultman, C.M., 2012. Preterm birth and psychiatric disorders in young adult life. Archives of general psychiatry 69, 610–617. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1001/archgenpsychiatry.2011.1374&link_type=DOI) 60. Papini, C., White, T., Montagna, A., Brittain, P., Froudist-Walsh, S., Kroll, J., Karolis, V., Simonelli, A., Williams, S., Murray, R., Nosarti, C., 2016. Altered resting-state functional connectivity in emotion-processing brain regions in adults who were born very preterm. Psychological Medicine 46, 3025–3039. [PubMed](http://medrxiv.org/lookup/external-ref?access_num=27523311&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F03%2F21%2F2022.03.18.22272505.atom) 61. Piguet, C., Karahanoğlu, F.I., Saccaro, L.F., Van De Ville, D., Vuilleumier, P., 2021. Mood disorders disrupt the functional dynamics, not spatial organization of brain resting state networks. NeuroImage: Clinical 32, 102833. 62. Power, J.D., Barnes, K.A., Snyder, A.Z., Schlaggar, B.L., Petersen, S.E., 2012. Spurious but systematic correlations in functional connectivity MRI networks arise from subject motion. Neuroimage 59, 2142–2154. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.neuroimage.2011.10.018&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=22019881&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F03%2F21%2F2022.03.18.22272505.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000299494000017&link_type=ISI) 63. Power, J.D., Mitra, A., Laumann, T.O., Snyder, A.Z., Schlaggar, B.L., Petersen, S.E., 2014. Methods to detect, characterize, and remove motion artifact in resting state fMRI. NeuroImage 84, 320–341. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.neuroimage.2013.08.048&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=23994314&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F03%2F21%2F2022.03.18.22272505.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000328868600030&link_type=ISI) 64. Preti, M.G., Bolton, T.A., Van De Ville, D., 2017. The dynamic functional connectome: State-of-the-art and perspectives. Neuroimage 160, 41–54. 65. Räikkönen, K., Pesonen, A.-K., Heinonen, K., Kajantie, E., Hovi, P., Järvenpää, A.-L., Eriksson, J.G., Andersson, S., 2008. Depression in young adults with very low birth weight: the Helsinki study of very low-birth-weight adults. Archives of general psychiatry 65, 290–296. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1001/archgenpsychiatry.2007.40&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=18316675&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F03%2F21%2F2022.03.18.22272505.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000253672200006&link_type=ISI) 66. Reyes, L.M., Jaekel, J., Bartmann, P., Wolke, D., 2021. Peer relationship trajectories in very preterm and term individuals from childhood to early adulthood. Journal of Developmental Behavioral Pediatrics 42, 621–630. 67. Rohr, C.S., Arora, A., Cho, I.Y., Katlariwala, P., Dimond, D., Dewey, D., Bray, S., 2018. Functional network integration and attention skills in young children. Developmental cognitive neuroscience 30, 200–211. 68. Rowlands, M.A., Scheinost, D., Lacadie, C., Vohr, B., Li, F., Schneider, K.C., Constable, R.T., Ment, L.R., 2016. Language at rest: A longitudinal study of intrinsic functional connectivity in preterm children. NeuroImage: Clinical 11, 149–157. 69. Saigal, S., Day, K.L., Van Lieshout, R.J., Schmidt, L.A., Morrison, K.M., Boyle, M.H., 2016. Health, wealth, social integration, and sexuality of extremely low-birth-weight prematurely born adults in the fourth decade of life. JAMA pediatrics 170, 678–686. 70. Siffredi, V., Farouj, Y., Tarun, A., Anderson, V., Wood, A.G., McIlroy, A., Leventer, R.J., Spencer-Smith, M.M., Van De Ville, D., 2021. Large-scale functional network dynamics in human callosal agenesis: Increased subcortical involvement and preserved laterality. Neuroimage 243, 118471. 71. Skudlarski, P., Jagannathan, K., Calhoun, V.D., Hampson, M., Skudlarska, B.A., Pearlson, G., 2008. Measuring brain connectivity: diffusion tensor imaging validates resting state temporal correlations. Neuroimage 43, 554–561. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.neuroimage.2008.07.063&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=18771736&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F03%2F21%2F2022.03.18.22272505.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000262300200014&link_type=ISI) 72. Stoecklein, S., Hilgendorff, A., Li, M., Förster, K., Flemmer, A.W., Galiè, F., Wunderlich, S., Wang, D., Stein, S., Ehrhardt, H., 2020. Variable functional connectivity architecture of the preterm human brain: impact of developmental cortical expansion and maturation. Proceedings of the National Academy of Sciences 117, 1201–1206. [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NDoicG5hcyI7czo1OiJyZXNpZCI7czoxMDoiMTE3LzIvMTIwMSI7czo0OiJhdG9tIjtzOjUwOiIvbWVkcnhpdi9lYXJseS8yMDIyLzAzLzIxLzIwMjIuMDMuMTguMjIyNzI1MDUuYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9) 73. Tamnes, C.K., Roalf, D.R., Goddings, A.-L., Lebel, C., 2018. Diffusion MRI of white matter microstructure development in childhood and adolescence: Methods, challenges and progress. Developmental cognitive neuroscience 33, 161–175. 74. 1. RStudio, I. Team, R., 2020. RStudio: Integrated Development for R. In: RStudio, I. (Ed.), Boston, MA, p. [http://www.rstudio.com/](http://www.rstudio.com/). 75. Team, R.C., 2013. R: A language and environment for statistical computing. 76. Treyvaud, K., Ure, A., Doyle, L.W., Lee, K.J., Rogers, C.E., Kidokoro, H., Inder, T.E., Anderson, P.J., 2013. Psychiatric outcomes at age seven for very preterm children: rates and predictors. Journal of Child Psychology 54, 772–779. 77. Van Den Heuvel, M.P., Pol, H.E.H., 2010. Exploring the brain network: a review on resting-state fMRI functional connectivity. European neuropsychopharmacology 20, 519–534. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.euroneuro.2010.03.008&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=20471808&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F03%2F21%2F2022.03.18.22272505.atom) 78. Volpe, J.J., 2009. Brain injury in premature infants: a complex amalgam of destructive and developmental disturbances. The Lancet Neurology 8, 110–124. 79. Walshe, M., Rifkin, L., Rooney, M., Healy, E., Nosarti, C., Wyatt, J., Stahl, D., Murray, R., Allin, M., 2008. Psychiatric disorder in young adults born very preterm: role of family history. European Psychiatry 23, 527–531. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.eurpsy.2008.06.004&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=18774275&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F03%2F21%2F2022.03.18.22272505.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000260366600012&link_type=ISI) 80. Wehrle, F.M., Michels, L., Guggenberger, R., Huber, R., Latal, B., O’Gorman, R.L., Hagmann, C.F., 2018. Altered resting-state functional connectivity in children and adolescents born very preterm. NeuroImage: Clinical 20, 1148–1156. 81. White, T.P., Symington, I., Castellanos, N.P., Brittain, P.J., Walsh, S.F., Nam, K.-W., Sato, J.R., Allin, M.P., Shergill, S.S., Murray, R.M., Williams, J.M.G., Nosarti, C., 2014. Dysconnectivity of neurocognitive networks at rest in very-preterm born adults. NeuroImage: Clinical 4, 352–365. 82. Wilke, M., Hauser, T.K., Krägeloh-Mann, I., Lidzba, K., 2014. Specific impairment of functional connectivity between language regions in former early preterms. Human Brain Mapping 35, 3372–3384. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1002/hbm.22408&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=24243552&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F03%2F21%2F2022.03.18.22272505.atom) 83. Witt, A., Theurel, A., Tolsa, C.B., Lejeune, F., Fernandes, L., van Hanswijck de Jonge, L., Monnier, M., Graz, M.B., Barisnikov, K., Gentaz, E., 2014. Emotional and effortful control abilities in 42-month-old very preterm and full-term children. Early Hum Dev 90, 565–569. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.earlhumdev.2014.07.008&link_type=DOI) 84. Zoller, D., Sandini, C., Karahanoglu, F.I., Padula, M.C., Schaer, M., Eliez, S., Van De Ville, D., 2019. Large-Scale Brain Network Dynamics Provide a Measure of Psychosis and Anxiety in 22q11.2 Deletion Syndrome. Biological Psychiatry: Cognitive Neuroscience and Neuroimaging 4, 881–892 85. Zöller, D., Sandini, C., Karahanoğlu, F.I., Padula, M.C., Schaer, M., Eliez, S., Van De Ville, D., 2019. Large-scale brain network dynamics provide a measure of psychosis and anxiety in 22q11. 2 deletion syndrome. Biological Psychiatry: Cognitive Neuroscience Neuroimaging 4, 881–892. 86. Zöller, D.M., Bolton, T.A., Karahanoğlu, F.I., Eliez, S., Schaer, M., Van De Ville, D., 2018. Robust recovery of temporal overlap between network activity using transient-informed spatio-temporal regression. IEEE Transactions on Medical Imaging 38, 291–302. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1109/TMI.2018.2863944&link_type=DOI)