key: cord-0729520-bi07wvub
authors: Liu, Siman; Wang, Yin-Shan; Zhang, Qing; Zhou, Quan; Cao, Li-Zhi; Jiang, Chao; Zhang, Zhe; Yang, Ning; Dong, Qi; Zuo, Xi-Nian
title: Chinese Color Nest Project : An accelerated longitudinal brain-mind cohort
date: 2021-10-11
journal: Dev Cogn Neurosci
DOI: 10.1016/j.dcn.2021.101020
sha: 75130c25182b41fc2cb12962133a6d2623f2cf71
doc_id: 729520
cord_uid: bi07wvub

The ongoing Chinese Color Nest Project (CCNP) was established to create normative charts for brain structure and function across the human lifespan, and link age-related changes in brain imaging measures to psychological assessments of behavior, cognition, and emotion using an accelerated longitudinal design. In the initial stage, CCNP aims to recruit 1520 healthy individuals (6–90 years), which comprises three phases: developing (devCCNP: 6–18 years, N = 480), maturing (matCCNP: 20–60 years, N = 560) and aging (ageCCNP: 60–84 years, N = 480). In this paper, we present an overview of the devCCNP, including study design, participants, data collection and preliminary findings. The devCCNP has acquired data with three repeated measurements from 2013 to 2017 in Southwest University, Chongqing, China (CCNP-SWU, N = 201). It has been accumulating baseline data since July 2018 and the second wave data since September 2020 in Chinese Academy of Sciences, Beijing, China (CCNP-CAS, N = 168). Each participant in devCCNP was followed up for 2.5 years at 1.25-year intervals. The devCCNP obtained longitudinal neuroimaging, biophysical, social, behavioral and cognitive data via MRI, parent- and self-reported questionnaires, behavioral assessments, and computer tasks. Additionally, data were collected on children’s learning, daily life and emotional states during the COVID-19 pandemic in 2020. We address data harmonization across the two sites and demonstrated its promise of characterizing the growth curves for the overall brain morphometry using multi-center longitudinal data. CCNP data will be shared via the National Science Data Bank and requests for further information on collaboration and data sharing are encouraged.

Construction of normative charts for brain structure and function across the human lifespan in typically developing individuals is essential for characterizing brain morphology and connectivity changes and monitoring neurodevelopmental status (Harms et al., 2018; Peterson et al., 2018) . The age-related trajectory of the brain in the life cycle can be used to elucidate the normative pattern of brain development, maturation and aging, and further reveal similarities and differences at different developmental stages (Betzel et al., 2014; Courchesne et al., 2000; Zuo et al., 2010 Zuo et al., , 2017 . Growth charts for height, weight and head circumference provide important information about children's physical health, which enables pediatricians and parents to identify abnormal physical growth (Cole, 2012; Kløvgaard et al., 2018; Ouyang et al., 2018; Yang et al., 2015) . Similarly, inadequate or excessive brain growth as estimated by deviations from the normalized distribution in growth curves may represent abnormal neurodevelopmental status that needs closer monitoring or interventions (Peterson et al., 2018) .

With respect to the use of brain growth charts in clinical practices, many neuropsychiatric disorders, such as autism spectrum disorder (ASD), psychotic disorder, bipolar disorder and attention deficit hyperactivity disorder (ADHD), are theorized to have origins in abnormal brain development (Ecker et al., 2015; Hazlett et al., 2017; Jalbrzikowski et al., 2019; Norman et al., 2016; Skåtun et al., 2016) . The presence and characterization of deviations from normative brain trajectories related to neuropsychiatric conditions may help reveal associated pathological and physiological changes, allowing for the construction of neuroimaging biomarkers and shedding new light on clinical diagnosis and treatment Shaw, 2016; Zabihi et al., 2019; Ziegler et al., 2014) . From an intervention perspective, brain growth charts also provide optimal periods for intervention and normative age-expected values that are required for further assessment of the effectiveness of an intervention program (Kulkarni et al., 2017) . If participants with abnormal brain development return to follow the normative pattern after receiving a specific intervention, this intervention could be considered effective in attenuating neurocognitive impairment or promoting brain development and successful aging.

Given the implications of brain growth charts in basic research and clinical application, researchers have sought to establish brain growth charts and use them to predict cognitive deficits and neuropsychiatric conditions (Erus et al., 2015; Kessler et al., 2016) . However, reliable and accurate charts for the brain are still not available (see a recent seminal work from Bethlehem et al., 2021) . The underlying reason is that most studies rely on cross-sectional design and focus on a narrow range of ages and sample sizes, which limit their ability to accurately reflect the continuity of individual development over the life course. In cross-sectional studies, the age-related changes shown in developmental trajectories often confound age and cohort effects (Spuling et al., 2015; Sutin et al., 2013; Twenge et al., 2017) . Standard (i.e., pure or single-cohort) longitudinal design overcomes the shortcomings of cross-sectional design by recruiting a single cohort of participants of the same age and following them over a period of time. However, this design inevitably raises issues of time demands and sample attrition over time, and also elicits practice effects due to repeated assessments, which reduces its feasibility and efficiency in life-span studies involving a wide range of age groups (Thompson et al., 2011; Hayward and Krause, 2016; Twenge et al., 2017; Zuo et al., 2017) .

The dilemma between cross-sectional and standard longitudinal designs highlights the need for using accelerated longitudinal designs (ALD), which recruit multiple cohorts of individuals of different ages and follow them over a shorter period of time (Galbraith et al., 2017; Joiner et al., 2018) . The multicohort approach is less time-consuming and more practical for creating growth curves over a long interval of the life course. Indeed, previous studies have mapped developmental trajectories over a time span of decades using ALD, the most famous of which is The Seattle Longitudinal Study, investigating longitudinally cognitive development from young adulthood to old age (Schaie et al., 2004; Hülür et al., 2019) . In the field of developmental cognitive neuroscience, researchers have used brain imaging methods to create accelerated longitudinal cohorts (see Table 1 for a non-exhaustive list). These studies have characterized the developmental trajectories of brain function and cognitive processes from childhood to adulthood, and further examine the neural correlates underlying the developmental changes of cognitive performances (Darki & Klingberg, 2015; Montez et al., 2017) . For example, Simmonds et al. (2017) conducted an accelerated longitudinal study over 10 years in 129 healthy participants, aged 8-30 years, who were asked to return for annual visits (mean = 2.8 visits/participant), and found that increases in encoding/retrieval activity in visual cortex accompanied increases in working memory accuracy while decreases in maintenance activity in prefrontal/subcortical regions accompanied decreases in working memory latency. Despite these advances, there are very few accelerated longitudinal cohorts (e. g., LCBC) available for characterizing brain and psycho-behavioral development across the postnatal lifespan. Meanwhile, previous studies revealed cultural or ethnic effects on development of brain and mind (Fan et al., 2015; Dong et al., 2020; Qiu, 2020) . Many large-scale neuroimaging samples include almost universally European or American (see Fig. 1 in Rosenberg et al., 2018) , missing cultural diversity.

The principal work of the ongoing Chinese Color Nest Project (CCNP) is to collect high-quality data on brain magnetic resonance imaging (MRI) scans and socio-psycho-behavioral factors using an ALD encompassing three repeated measurements on healthy individuals aged 6-84 years at enrollment. The overall aim of CCNP is to build normative charts for brain structure and function across the human lifespan, and to further examine how age-related changes in brain anatomical and function are related to those in psycho-behavioral functions, such as cognitive and behavioral abilities. As illustrated in Fig. 1A , the general CCNP design comprises three phases: developing (devCCNP: 6-18 years, N = 24 n), maturing (matCCNP: 18-60 years, N = 28 n) and aging (ageCCNP: 60-84 years, N = 24 n) and n is the number of sex-specific participants (female or male) in each age cohort. At its pilot stage, CCNP targets an accelerated longitudinal cohort of 1520 healthy participants (i.e., n = 20) implemented at two cities (Chongqing and (Russell et al., 2021; Giedd et al., 2015; Evans, 2006) Beijing). Specifically, the devCCNP started in 2013, which has successfully acquired CCNP-Southwest University (CCNP-SWU) samples with three repeated measurements at the trial stage (2013) (2014) (2015) (2016) (2017) . At the second stage, devCCNP is collecting CCNP-Chinese Academy of Sciences (CCNP-CAS) baseline data since July 2018 and the second-wave data since September 2020. Meanwhile, the matCCNP was initiated in 2015 and has collected more than 1000 repeated measures in Beijing. This dataset will be published soon as part of an open resource called 3R-BRAIN, for reliable, reproducible and replicable brain imaging research (Zuo et al., , 2019 . The ageCCNP will be initiated in 2022 at Beijing and Chongqing, respectively. In this paper, we present an overview of devCCNP, including the study design, participants, data collection and preliminary findings.

Before we started this project, we were finishing a high-resolution (4 mm) network map of the human brain connectome using functional imaging data from more than 1000 healthy volunteers across different stages of the postnatal lifespan (6-85 years) based upon the 1000 Functional Connectomes Project (Biswal et al., 2010; Zuo et al., 2012) . This map was shaped like a colored nest and inspired us to name the project (Fig. 1B ). In addition, color nest represents protection, vitality and development, and implies the concept of respecting individual differences of human beings and promoting healthy development and well-being. CCNP is devoted to collecting data on brain structure and function across different stages of the human lifespan. While CCNP was design comprises three phases: developing (devCCNP: 6-18 years, N=480, 12 age cohorts, interval = 15 months), maturing (matCCNP: 18-60 years, N=560, 14 age cohorts, interval = 39 months) and aging (ageCCNP: 60-84 years, N=480, 12 age cohorts, interval = 27 months). It has three waves of measurements including baseline at Wave 1 (purple), followup1 at Wave 2 (blue) and followup2 at Wave 3 (green). There are 40 participants (20 males and 20 females) in each age cohort. (B) A high-resolution network map of the human brain connectome, which was shaped like a colored nest, and inspired us to name the project. currently at its pilot stage and designed to cover the postnatal development (6-90 years), we plan to expand the CCNP design to include participants much younger in future, and even fetuses.

The CCNP uses the accelerated longitudinal design (ALD), which allows controlling for cohort and time of measurement effects with the additional advantage of collecting data on a wider age range in less time (Thompson et al., 2011; Rioux and Little, 2020) . This structured multi-cohort longitudinal design is particularly advantageous for lifespan trajectory studies, and optimal when accounting for missing data. As illustrated in Fig. 1A , at the pilot stage, CCNP distributes the ALD (baseline age = 6-18 years, 12 age cohorts, 3 waves, interval = 15 months) to devCCNP, and the ALD (baseline age = 60-84 years, 12 age cohorts, 3 waves, interval = 27 months) to ageALD while matALD receives the ALD (baseline age = 18-60 years, 14 age cohorts, 3 waves, interval = 39 months). While attrition is an important factor for us to decide the specific sampling intervals, we also consider other scientific factors for developmental brain-mind association studies. The CCNP design ensures an overlap of two measurements between successive age cohorts. The specific sampling intervals indicate that the duration between two neighboring waves is 1.25 years, 2.25 years or 3.25 years, i.e., integral years and one season, preventing multiple visits from the same participant during a season (Spring, Summer, Fall, Winter). As a result of the season effect balanced, for a certain participant, three repeated measurement occasions (i.e., waves) were distributed in various seasons. We provided additional details on the ALD design and participants below.

The trial stage implementation of devCCNP took place at Beibei District, Chongqing, China from March 2013 to January 2017. The Institute of Psychology, Chinese Academy of Sciences (CAS) and the Faculty of Psychology, Southwest University (SWU) are responsible for this pilot study called CCNP-SWU. The ongoing CCNP-CAS began in September 2017 in Beijing, China. Ethical approval was provided by the Institutional Review Board of the CAS Institute of Psychology. To provide better representation of the Chinese population, we aimed to include participants from cities differing in economy and region of the country. CCNP-SWU recruited participants in Chongqing, a second-tier city located in the southwest of China. By contrast, CCNP-CAS recruited participants in Beijing, a first-tier city located in the north of China. Going forward, CCNP will recruit additional participants from more cities (48 cities in next ten years). Furthermore, to ensure demographic consistency across the three phases of CCNP (i.e., devCCNP, matCCNP and ageCCNP), each sub-sample from one of the three phases (e.g., devCCNP) would be recruited to match to other sub-sample from different phases (e.g., matCCNP or age CCNP) across a broad range of social and demographic features. We noted that the measures in the devCCNP would be compared to those that are being used in the matCCNP and ageCCNP if the measures are used for all age groups (e.g., biophysical measures, and MRI scans), or a set of tests designed and normed for different age groups, such as Eysenck Personality Questionnaire and Wechsler Intelligence Scale.

The pilot devCCNP aimed to recruit 480 age-sex stratified children distributed in 12 age cohorts at Chongqing ( Fig. 2A , CCNP-SWU: 8 boys and 8 girls in each cohort) and Beijing (Fig. 2B , CCNP-CAS: 12 boys and 12 girls in each cohort). A priori power analysis conducted in G*Power 3.1 (Faul et al., 2009 ) suggested that a sample size of 123 was necessary to achieve the specified statistical power of 0.95 (two tails, a medium effect size f2 = 0.15, alpha = 0.01, and three predictors) for a single regression coefficient in linear multiple regression. It is challenging to estimate sample sizes for building growth curves, because it is potentially more complicate and thus need much larger sample than a common regression. During the pilot stage, the consortium aims to obtain 4560 samples from 1520 participants, which is about 37 times of the common regression derived sample size.

The measures used in devCCNP across all waves to assess biophysical health and socio-psycho-behavioral factors are documented in Table 2 . Objective biophysical measures were collected at each wave, including height, weight, head circumference, and biomarkers of cardiovascular health (i.e., blood pressure and pulse). Regarding psychological measures, the study primarily focused on cognition, personality, and issues pertaining to social-emotional functioning. Widely used instruments with high reliability and validity were selected. The questionnaires were primarily built to collect socio-psycho-behavioral data on sociodemographic characteristics, child behavioral problems (Achenbach and Edelbrock, 1983) , life events (Liu et al., 1991) , self-concept (Piers, 1984) , emotions and affects such as stress, anxiety, depression, loneliness, and positive and negative affect (Asher et al., 1984; Cohen et al., 1983; Kovacs, 1992; La Greca et al., 1988; March et al., 1997; Spielberger, 1983; Watson et al., 1988) , emotional quotient (Bar-On and Parker, 2000), reading (Xue et al., 2013) , personality characteristics (Eysenck and Eysenck, 1975) , creativity (Hwang et al., 2007; Torrance, 1984) , as well as video game use. Questionnaires on sociodemographic information and CBCL were filled out by the parents, while the rest of (B) devCCNP-CAS contains 288 participants (12 boys and 12 girls in each cohort, 12 age cohorts). Three measuring waves are represented by purple (Wave 1), blue (Wave 2) and green (Wave 3) boxes. The duration between two neighboring waves is designed as 1.25 years for the purpose of avoiding a kid always visited at a same season. Each cohort covers one integer age-interval at its baseline (C01: 6-7 years; C02: 7-8 years; …; C12: 17-18 years). the questionnaires were completed by the children.

The Chinese Version of WISC-IV was used to measure cognitive ability (Baron, 2005; Watkins and Smith, 2013; Wechsler, 2003; Zhang, 2009 ). The scale consists of 10 core subtests and 4 supplementary subtests, which provides a full-scale Intelligence Quotient (FSIQ) reflecting the general cognitive ability and four indices: Verbal Comprehension Index (VCI), Perceptual Reasoning Index (PRI), Working Memory Index (WMI) and Processing Speed Index (PSI). Three E-Prime computer tasks, including the Attention Network Test, Task-Switch Paradigm and Digit N-back Task, were used to assess children's alerting, orienting, executive attention, cognitive flexibility and working memory (Fan et al., 2002; Schuch and Koch, 2003; Zhu et al., 2006) .

Imaging data at CCNP-SWU were collected across all waves on a Siemens Trio 3 T scanner at Southwest University, while imaging data at CCNP-CAS were collected on a GE Discovery MR750 3 T scanner at the CAS Institute of Psychology. Two resting-state functional MRI (rfMRI) scans were acquired and separated by a T1-weighted MRI scan. The main acquisition parameters for T1-weighted MRI and rfMRI are listed in Supplementary Table S2 . A T2-weighted MRI scan was also performed to evaluate brain lesions and improve cross-registration. All the Child Questionnaires and Cognitive Battery Tests were completed after the imaging scans, in order to avoid the influence of physiological fatigue and affective state. Structural MRI images were first anonymized to remove all the personal information from the raw MRI data. We obscured the facial information by using the face-masking tool (Milchenko and Marcus, 2013) customized with the Chinese pediatric templates developed by our lab (Dong et al., 2020) . The anonymized images were then denoised by spatially adaptive non-local means and corrected for intensity normalization in the Connectome Computation System (Xu et al., 2015) . To extract individual brains, we trained a deep learning method ) using a small set of semi-automatically extracted brains in the CCNP-SWU samples, and then applied to all the devCCNP samples. The preprocessed brain volumes were all in the native space and fed into FreeSurfer (version 6.0) pipeline (Fischl, 2012) to obtain general morphological measurements of different brain tissues from aseg.stats and lh/rh.aparc.stats including estimated intracranial volume (eTIV), total gray matter volume (GMV), subcortical gray matter volume (sGMV), cerebral white matter volume (cWMV), mean cortical thickness (CT) and white matter surface area (SA).

It is important to consider MRI instrument and data processing factors when designing multi-site longitudinal studies. We implemented the same ALD for multiple subprojects (i.e., devCCNP-SWU and devCCNP-CAS). Different age groups (cohorts) were tested on same scanners for same subprojects while not on the same scanners for different subprojects. To reduce variability or bias introduced by different scanning platforms, participants from different subprojects would be tested with the same field strength and acquisition parameters. Data and its analytic protocol are harmonized and optimized across different scanners or subprojects for achieving highly reliable multimodal MRI measurements. To ensure growth charts minimally affected by scanning platform, we applied the longitudinal ComBat method (Beer et al., 2020) to harmonize the extracted morphological measurements prior to the growth curve modelling. Generalized Additive Mixed Model (GAMM) was used to model the brain morphological growth for each sex with the harmonized samples from both devCCNP-SWU and devCCNP-CAS (van Duijvenvoorde et al., 2019). The GAMMs took participant as a random intercept by mgcv package in R 4.1.

In devCCNP-SWU, both children and their parents/guardians were approached through lectures that explained this large-scale longitudinal cohort and science popularization activities regarding brain development in primary and middle/high schools for students in grades 1-11 (i. e., grade 2 in senior high school). By contrast, devCCNP-CAS comprises community samples recruited via the public uses social media wechat and science popularization activities. The exclusion criteria included the following: (a) children's birth weight less than 1.5 kg or greater than 4.2 kg; (b) children with a gestational or perinatal history of pathologies or risk factors; (c) children's height, weight, or head circumference less than the third percentile of the growth curve of Chinese children and adolescents; (d) children with a history of head injury; (e) children or one of their family members diagnosed with a neuropsychiatric disorder, such as schizophrenia, ASD, ADHD, bipolar disorder, alcohol or other drug use disorder; (f) Wechsler Intelligence Scale for Children-IV (WISC-IV) full-scale Intelligence Quotient (FSIQ) standard score lower than 80. Both children and their parents/guardians volunteered to participate in this study and signed the informed consent before participation.

The devCCNP-SWU recruitment took place between March and December in 2013. A total of 198 children volunteered to take part in this study at baseline, of which 6 were excluded due to brain cyst (n = 2), depression disorder (n = 1), claustrophobia (n = 1), or WISC-IV FSIQ standard score below 80 (n = 2). Therefore, 192 children were eligible for further participation, which formed the sample for the baseline assessment of CCNP-SWU. Because of attrition, of the 192 original samples at baseline (Wave 1), 152 participated in the 15-month follow-up (Wave 2), representing a 79.17% re-participation rate from Wave 1 to Wave 2. Seven children were further included at Wave 2; of the 159 children participating at Wave 2, 106 retained in the 30-month follow-up (Wave 3), representing a 66.67% re-participation rate from Wave 2 to Wave 3. Also, 2 children were further included at Wave 3. Eventually, devCCNP-SWU sample consisted of 201 children who participated at least one wave, of which 100 (49.75%) participated in all three waves, 58 (28.86%) participated in two waves, and 43 (21.39%) participated in only one wave. Data collection of devCCNP-SWU for Wave 1 was conducted from December 2013 to July 2014; data collection for Wave 2 was conducted from April to August in 2015; data collection for Wave 3 was conducted from September 2016 to January 2017. The mean time lag between Wave 1 and Wave 2 was 1.27 years (SD = 0.06; range = 1.16-1.50 years), and between Wave 2 and Wave 3 was 1.37 years (SD = 0.07; range = 1.22-1.56 years). Attrition analyses were conducted to compare the devCCNP-SWU families who participated in all three measurement waves (n = 100) with these who participated in Wave 1 and then dropped out at either Wave 2 or Wave 3 (n = 92) on sociodemographic characteristics (see Supplementary Table S1 ). The attrition analyses indicated that the retained families did not differ from the drop-out families in the distribution of child sex, family size, urban areas or rural areas, and parental education status. However, children who attended measurement 3 times were younger than those who participated in Wave 1 and subsequently dropped out. We also conducted attrition analyses on main psychological outcomes (i.e., problem total scores, personality, social anxiety, depression, and WISC-IV scores) and brain morphological indexes. The results showed that there were no significant differences in the main psychological and brain morphological measures between retained families and drop-out families after FDR correction (FDR corrected ps > 0.075).

The devCCNP-CAS has been recruiting participants since July 2018 and has enrolled 338 children thus far. One child was excluded due to WISC-IV FSIQ standard score below 80. The baseline data have been obtained from 168 eligible children since July 2018, of which 57 have already completed data collection at Wave 2 since September 2020. The rest of the children who have been enrolled but not assessed would be invited to complete the assessment as soon as possible. Fig. 3 presents age and sex distributions of devCCNP-SWU and devCCNP-CAS samples. Sociodemographic characteristics for devCCNP-SWU and devCCNP-CAS samples are summarized in Table 3 . Pearson's Chi-squared tests were used to compare sample characteristics between devCCNP-SWU and devCCNP-CAS samples. The results indicated that children in devCCNP-CAS were more likely to live in urban areas and had highly educated mothers and fathers than children in devCCNP-SWU.

After assessment at each wave, every family received a personal development report for participation. The development report provided feedback on five aspects of children's development, including physiological characteristics (i.e., height, weight, head circumference, blood pressure, and pulse), cognitive ability, social-emotional development (e. g., social anxiety, depression, stress perception and behavioral problems), personality, and brain development. The brain development report included the global tissue morphology (i.e., eTIV, GMV, sGMV, cWMV) and the network-level volumetric features (Yeo et al., 2011) . The parents were shown a percentile score for their child (see Fig. 4A and B for a sample report of a girl). We note that this is an objective report of the test results and blanket recommendations instead of interventional strategies that may alter the behaviors of the participants or their parents. In addition, we emphasized that the participant's performance was always affected by the physical and emotional state while their developmental processes are dynamic and more informative and valid based on the final reports derived from the 5-year following assessments.

Main variables of boys and girls across two age groups (younger group less than 12 years old vs. older group greater than 12 years old) at devCCNP-SWU baseline are presented in Table 4 . The raw scores of extraversion, neuroticism and psychoticism assessed by Eysenck Personality Questionnaire (EPQ) were used. Independent samples t-tests were conducted to compare the means of main variables between boys and girls across two age groups. There were sex differences in body mass index (BMI), CBCL problem total scores, depression, neuroticism, psychoticism, eTIV, GMV, sGMV, cWMV, CT and SA in younger group. As for the older group, there were significant sex differences in height, social anxiety measured by Social Anxiety Scale for Children (SASC), psychoticism, FSIQ, verbal comprehension, processing speed, eTIV, GMV, sGMV, cWMV, and SA. Table 5 presents main variables of boys and girls across two age groups at devCCNP-CAS baseline. There were detectable sex differences in processing speed, eTIV, GMV, sGMV, cWMV, and SA for younger group, and in cWMV and SA for older group.

During the outbreak of Coronavirus Disease-2019 (COVID-19), we did not collect devCCNP data from January to August in 2020. Data collection has returned to normal since September 2020. However, we designed questionnaires to collect data on children's learning and daily life during COVID-19 pandemic (see Table 6 ). A total of 66 children (mean age = 11.38 years, SD = 3.33; 32 boys, 34 girls) from devCCNP-SWU (n = 11) and devCCNP-CAS (n = 55) completed the questionnaires online from March to October in 2020. About three in five children (59.09%) felt a bit nervous or panic when COVID-19 spread rapidly. The COVID-19 pandemic exerted slight impact on daily routine for more than 60% of the households (62.12%), and exerted moderate or strong impact on daily routine for nearly 30% of the households (28.79%). As illustrated in Table 6 , more than half of the children reported that their learning efficiency, learning consciousness and exercise time were decreased. Meanwhile, screen time was dramatically increased for more than 60% of the children. Children's emotional states during COVID-19 pandemic were also measured by four scales: Depression Anxiety Stress Scale (DASS; Lovibond and Lovibond, 1995) , Children's Depression Inventory (CDI; Kovacs, 1992) , Depression Self-rating Scale for Children (DSRSC; Birleson, 1981) , and Childhood Anxiety Sensitivity Index (CASI; Silverman et al., 1999) . The descriptive statistics were presented in Supplementary Table S3 .

Growth curves are generated with both site-harmonized and raw data. While the overall shapes the growth curves are similar, the longitudinal ComBat improves the growth curve modeling, showing more powered longitudinal changes. In the main text, we visualized the growth curves using the ComBat derived data in Fig. 4 (see Supplementary Figure 1 for those using the raw data). Specifically, eTIV, sGMV and WMV exhibited a growth pattern with increasing volumes when growing up (Fig. 4C ). In contrast, CT was growing as a linearly thinning process with development while the growth curves of GMV and SA demonstrated somehow nonlinearity, indicating an inverted-U growing shape (Fig. 4D) .

A previous review was published in Chinese by our group, which comprehensively described the devCCNP-SWU protocol for experimental design, sample selection and data collection (Yang et al., 2017) . Using the devCCNP-SWU baseline brain imaging data from 84 participants, we previously reported that children exhibited similar region-specific asymmetry of dorsal anterior cingulate cortex (dACC) as in adults, and further revealed that dACC functional connectivity with default, frontoparietal and visual networks showed a region-specific asymmetry . This baseline data has been released as part of the Consortium for Reliability and Reproducibility (CoRR; Zuo et al., 2014) , the IPCAS 7 site (https://doi.org/10.15387/fcp_indi.corr. ipcas7), which has been listed as one of the existing, ongoing large-scale developmental dataset (Rosenberg et al., 2018) . Head motion data during mock-scanning from devCCNP-CAS are recently demonstrated with frequency-specific evidence to support motion as a developmental trait across children and adolescents by the development of a neuroinformatic tool, namely DREAM (Gong et al., 2021) .

3 Fig. 4 . Brain morphological growth curves in schoolchildren. Individual scores are plotted onto the normative growth charts of the global tissue morphological measurements (A) and volumetric measures of the seven common brain networks (see Yeo et al., 2011 for more details on the networks) (B). Growth curves are generated with their confidence intervals for the estimated intracranial volume (eTIV), subcortical gray matter volume (sGMV), cerebral white matter volume (WMV) (C), and total gray matter volume (GMV), mean cortical thickness (CT) and white matter surface area (SA) (D) males and females, respectively. These curves base the growth charts for the individual assessments.

The full set of devCCNP data is increasingly appreciated by collaborative studies on school-aged children and adolescents. It supports the building of age-specific cranio-cortical correspondences for school-aged children and adolescents . We recently observed that social anxiety was positively correlated with the GMV in an area of the orbital-frontal cortex, and its functional connectivity with the amygdala (Mao et al., 2020) . Using the longitudinal data, we also charted the growth curves of human amygdala across school ages (Zhou et al., 2021) . A standardized protocol on charting brain development during school age has been developed to generate the corresponding brain templates and model growth charts, revealing the differences in brain morphological growth between Chinese and American population, in particularly around puberty (Dong et al., 2020) . Using resting-state fMRI data from devCCNP-SWU, we revealed age-dependent changes in the macroscale organization of cortex. The findings suggest that the scheduled maturation of functional connectivity gradient shifts may be critically important for understanding how cognitive and behavioral capabilities are refined across development, marking puberty-related changes . As part of an international consortium for human lifespan brain chart recently initiated (Bethlehem et al., 2021) , CCNP contributes to the largest world-wide MRI samples (N > 120,000) for building the normative brain charts for the human lifespan (0-100 years).

CCNP is a cohort to combine longitudinal brain imaging, biophysical, social, behavioral and cognitive measures across nearly the full lifespan, enhancing the culture diversity among the existing accelerated longitudinal datasets (Table 1 ). The pilot CCNP demonstrated that the accelerated longitudinal and freely available data on high-resolution structural and functional images beginning in childhood and extending into early adulthood with the same imaging protocols have been practical to acquire. Such data are invaluable for the construction of normative growth charts to describe how brain morphological metrics and functional connectivity change as a function of child age across sexs Dong et al., 2020; Holla et al., 2020a Holla et al., , 2020b Qiu, 2020) . In addition, longitudinal psycho-behavioral data will enable the scientific community to investigate how age-related changes in brain structure and function are associated with age-related changes in cognitive function and behavioral performances. By harnessing big data, the field of developmental cognitive neuroscience is rapidly transitioning into a population science, namely developmental population cognitive neuroscience (Zuo et al., 2018) . This raises lots of challenges and opportunities (Becht and Mills, 2020; Fair et al., 2021; Marek et al., 2020; Mills and Tamnes, 2014) , such as the how to combine the existing data resources and delineating the individual differences in brain and mind development. One promising solution would be to plan more comprehensive and sophisticated large-scale national-level cohort based upon the rich experiences from the existing cohort. ABCD (Garavan et al., 2018) and Lifebrain (Walhovd et al., 2018) are two national projects using longitudinal designs in the United States and Europe, respectively. In China, the National Brain Project (Poo et al., 2016) has initiated a large longitudinal cohort on school-aged brain and mind development. The devCCNP has shared its experiences in building the national cohort, and in future, will distribute its initial practice into other 48 sites across the country for a final national representative sample including 50 sites (i.e., 38,000 participants, 114,000 samples). This will build normative charts for translating the basic research on individual differences in brain and mind development into educational and clinical conditions across the lifespan (Bethlehem et al., 2021) . We also hope that the final CCNP will offer a unique resource for examining cultural or ethnic effects on development of brain and mind (Fan et al., 2015; Qiu, 2020) .

Main variables of boys and girls across two age groups at devCCNP-SWU baseline (N = 192 This study has several limitations. First, we are aware of the inevitable sample attrition across three waves in CCNP-SWU. Although attrition analyses suggested that child sex and family background factors did not contribute to the attrition, children who completed all three waves were younger than those who participated in Wave 1 and then dropped out of the study at either Wave 2 or Wave 3, which may introduce bias, especially when constructing growth curves and describing age-related changes. Second, CCNP focuses exclusively on Chinese Han population, and therefore the findings may not be generalizable to other ethnic groups. There is a physiological variation between different ethnic groups, and cognitive neuroscience research on ethnic minority groups in China is also an important direction for future studies. Finally, the sample is limited to a few regions, which may limit its ability to represent Chinese general population. Due to the complexity and high cost of MRI brain scanning, it is difficult to carry out large-scale and coordinated nationwide research, which highlights the need for conducting collaborative research and promoting MRI data sharing.

Baseline CCNP-SWU data on brain imaging are available to researchers via the CoRR, which is committed to open science by aggregating and sharing MRI data from multiple sources to establish testretest reliability and reproducibility in functional connectomics . The rest of the data will be publicly shared via the National Science Data Bank and fully available to the research community when acquisition is completed for the pilot CCNP. At this stage, data are only available to researchers and collaborators of CCNP. More information about CCNP can be found at: http://deepneuro.bnu.edu.cn/?p=163 or https://github.com/zuoxinian/CCNP. Requests for further information and collaboration are encouraged and considered by principal investigator Xi-Nian Zuo [xinian.zuo@bnu.edu.cn] .

The Chinese Color Nest Consortium members are at http://deepne uro.bnu.edu.cn/?p=163.

This work was supported by the Start-up Funds for Leading Talents at Beijing Normal University, the Key-Area Research and Development Program of Guangdong Province (2019B030335001), the National Basic Science Data Center "Chinese Data-sharing Warehouse for In-vivo Imaging Brain" (NBSDC-DB-15), the Beijing Municipal Science and Technology Commission (Z161100002616023, Z181100001518003), the Major Project of National Social Science Foundation of China (20&ZD296), the CAS-NWO Programme (153111KYSB20160020), the Guangxi BaGui Scholarship (201621) 

The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper.

Baseline CCNP-SWU data on brain imaging are available to researchers via the Consortium for Reliability and Reproducibility (https://doi.org/10.15387/fcp_indi.corr.ipcas7), which is committed to open science by aggregating and sharing MRI data from multiple sources to establish test-retest reliability and reproducibility in functional connectomics. The rest of the data will be publicly shared and fully available to the research community when acquisition is completed. At this stage, data are only available to researchers and collaborators of CCNP. More information about CCNP can be found at https://github.com/zuoxinian/CCNP. Requests for further information and collaboration are encouraged and considered by the principal investigator (PI) Xi-Nian Zuo [xinian.zuo@bnu.edu.cn or zuoxn@psych.ac.cn].

Manual for the Child Behavior Checklist and Revised Child Behavior Profile

Loneliness in children

Test review: Wechsler Intelligence Scale for Children-Fourth Edition (WISC-IV)

BarOn Emotional Quotient Inventory: Youth Version. Multi-Health system

Modeling individual differences in brain development

Alzheimer's Disease Neuroimaging, I., 2020. Longitudinal ComBat: a method for harmonizing longitudinal multi-scanner imaging data

Brain charts for the human lifespan

Dispersion of functional gradients across the adult lifespan

Changes in structural and functional connectivity among resting-state networks across the human lifespan

The validity of depressive disorder in childhood and the development of a self-rating scale: a research report

Toward discovery science of human brain function

Development of hippocampal-prefrontal cortex interactions through adolescence

Longitudinal development of hippocampal subregions from early-to mid-childhood

Modeling longitudinal changes in hippocampal subfields and relations with memory from early-to mid-childhood

Brain mechanisms underlying neuropsychiatric symptoms in Alzheimer's disease: a systematic review of symptom-general and -specific lesion patterns

A global measure of perceived stress

The development of growth references and growth charts

Normal brain development and aging: quantitative analysis at in vivo MR imaging in healthy volunteers

The americleft psychosocial outcomes project: a multicenter approach to advancing psychosocial outcomes for youth with cleft lip and palate

The role of fronto-parietal and fronto-striatal networks in the development of working memory: a longitudinal study

Charting brain growth in tandem with brain templates for schoolchildren

Shifting gradients of macroscale cortical organization mark the transition from childhood to adolescence

Neuroimaging in autism spectrum disorder: brain structure and function across the lifespan

Imaging patterns of brain development and their relationship to cognition

Brain development cooperative group. The NIH MRI study of normal brain development

Manual of Eysenck Personality Questionnaire (Adult and Junior)

Neuroscience in the Era of Networks and Big Data: Strengths, Weaknesses, Opportunities, and Threats

Modeling the 3D geometry of the cortical surface with genetic ancestry

Development of the default-mode network during childhood and adolescence: a longitudinal resting-state fMRI study

Testing the efficiency and independence of attentional networks

Statistical power analyses using G*Power 3.1: tests for correlation and regression analyses

FreeSurfer

Development and aging of cortical thickness correspond to genetic organization patterns

The genetic organization of longitudinal subcortical

Accelerated longitudinal designs: an overview of modelling, power, costs and handling missing data

White matter microstructural change contributes to worse cognitive function in patients with type 2 diabetes

Recruiting the ABCD sample: design considerations and procedures

Child psychiatry branch of the National Institute of Mental Health longitudinal structural magnetic resonance imaging study of human brain development

DREAM: a toolbox to decode rhythms of the brain system

Mapping domain-and age-specific functional brain activity for children's cognitive and affective development

Forms of attrition in a longitudinal study of religion and health in older adults and implications for sample bias

Early brain development in infants at high risk for autism spectrum disorder

Age, cohort, and period effects on metamemory beliefs

A series of five population-specific Indian brain templates and atlases spanning ages 6-60 years

Prediction complements explanation in understanding the developing brain

Cohort differences in cognitive aging: the role of perceived work environment

Multiple representation skills and creativity effects on mathematical problem solving using a multimedia whiteboard system

Population normative models of human brain growth across development

Structural brain alterations in youth with psychosis and bipolar spectrum symptoms

Affective experience across the adult lifespan: an accelerated longitudinal design

Growth charting of brain connectivity networks and the identification of attention impairment in youth

Growth of children in Greenland exceeds the World Health Organization growth charts

Children's Depression inventory: Manual. Multi-Health Systems

Brain developmental trajectories in children and young adults with isolated cleft lip and/or cleft palate

Endoscopic treatment versus shunting for infant hydrocephalus in Uganda

Development of the social anxiety scale for children: reliability and concurrent validity

Brain imaging research in autism spectrum disorders: in search of neuropathology and health across the lifespan

Longitudinal volumetric brain changes in autism spectrum disorder ages 6-35 years

Maturation of the human striatal dopamine system revealed by PET and quantitative MRI

Establishment of adolescent life event scale and reliability and validity test (in Chinese)

AIMS Longitudinal European Autism Project (LEAP): design and methodologies to identify and validate stratification biomarkers for autism spectrum disorders

The structure of negative emotional states: comparison of the Depression Anxiety Stress Scales (DASS) with the Beck Depression and Anxiety Inventories

The Multidimensional Anxiety Scale for Children (MASC): factor structure, reliability, and validity

Towards reproducible brainwide association studies

Conceptualizing mental disorders as deviations from normative functioning

Neuroimaging markers of adolescent depression in the National Consortium on Alcohol and Neurodevelopment in Adolescence (NCANDA) study

Obscuring surface anatomy in volumetric imaging data

Methods and considerations for longitudinal structural brain imaging analysis across development

The expression of established cognitive brain states stabilizes with working memory development

Neuropsychological development in adolescents: longitudinal associations with white matter microstructure

Structural and functional brain abnormalities in attention-deficit/hyperactivity disorder and obsessive-compulsive disorder: a comparative meta-analysis

The YOUth study: rationale, design, and study procedures

Growth patterns from birth to 24 months in Chinese children: a birth cohorts study across China

Dopamine-related striatal neurophysiology is associated with specialization of frontostriatal reward circuitry through adolescence

The Netherlands Study of Depression and Anxiety (NESDA): rationale, objectives and methods

Normative human brain volume growth

Piers-Harris Childrens Self-Concept Scale: Revised Manual

China brain project: basic neuroscience, brain diseases, and brain-inspired computing

A 16-year study of longitudinal volumetric brain development in males with autism

Habitual sleep is associated with both source memory and hippocampal subfield volume during early

Developmental variation in amygdala volumes: modeling differences across time, age, and puberty

Underused methods in developmental science to inform policy and practice

Child brain growth standard: age and ethnicity dependent

The Seattle longitudinal study: relationship between personality and cognition

The role of response selection for inhibition of task sets in task shifting

Global Imaging Genetics in Adolescents, C., 2019. Precision medicine and global mental health

The Cambridge Centre for Ageing and Neuroscience (Cam-CAN) study protocol: a cross-sectional, lifespan, multidisciplinary examination of healthy cognitive ageing

Maps of the development of the brain's functional architecture: could they provide growth charts for psychiatry?

Developing functional network connectivity of the dorsal anterior cingulate cortex mediates externalizing psychopathology in adolescents with child neglect

Factor structure of the childhood anxiety sensitivity index

Protracted development of executive and mnemonic brain systems underlying working memory in adolescence: a longitudinal fMRI study

Global brain connectivity alterations in patients with schizophrenia and bipolar spectrum disorders

Manual for the State-Trait Anxiety Inventory STAI (Form Y)

Changing predictors of self-rated health: disentangling age and cohort effects

The effect of birth cohort on well-being: the legacy of economic hard times

Intelligence development and school adjustment of school-age children and adolescents: a follow-up cohort study

Alzheimer's Disease Neuroimaging Initiative. Design considerations for characterizing psychiatric trajectories across the lifespan: application to effects of APOE-epsilon4 on cerebral cortical thickness in Alzheimer's disease

The role of creativity in identification of the gifted and talented

Open science datasets from PREVENT-AD

Age, time period, and birth cohort differences in self-esteem: reexamining a cohort-sequential longitudinal study

A three-wave longitudinal study of subcortical-cortical resting-state connectivity in adolescence: Testing age-and puberty-related changes

Longitudinal changes in amygdala, hippocampus and cortisol development following early caregiving adversity

Neurodevelopmental origins of lifespan changes in brain and cognition

Healthy minds 0-100 years: optimising the use of European brain imaging cohorts

Dorsal anterior cingulate cortex in typically developing children: laterality analysis

U-net model for brain extraction: trained on humans for transfer to non-human primates

Long-term stability of the wechsler intelligence scale for children-fourth edition

Development and validation of brief measures of positive and negative affect: the PANAS scales

Wechsler Intelligence Scale for Children-Fourth Edition. Psychological Corporation

A connectome computation system for discovery science of brain

The stability of literacy-related cognitive contributions to Chinese character naming and reading fluency

Comparison of the China growth charts with the WHO growth standards in assessing malnutrition of children

Chinese Color Nest Project (CCNP): growing up in China

Early prevention of cognitive impairment in the community population: the Beijing Aging Brain Rejuvenation Initiative

The organization of the human cerebral cortex estimated by intrinsic functional connectivity

Fractionating autism based on neuroanatomical normative modeling

Dissecting the heterogeneous cortical anatomy of autism spectrum disorder using normative models

The revision of WISC-IV Chinese version (in Chinese)

The Consortium on Vulnerability to Externalizing Disorders and Addictions

Transcranial brain atlas for school-aged children and adolescents

Tracing human amygdala across school age. bioRxiv

fMRI revealed neural substrate for reversible working memory dysfunction in subclinical hypothyroidism

Alzheimer's disease neuroimaging initiative. Individualized Gaussian process-based prediction and detection of local and global gray matter abnormalities in elderly subjects

Network centrality in the human functional connectome

Human connectomics across the life span

Growing together and growing apart: regional and sex differences in the lifespan developmental trajectories of functional homotopy

Developmental population neuroscience: emerging from ICHBD

Harnessing reliability for neuroscience research

Statistician Ting Xu is acknowledged for advice with data extraction and statistical analyses. Daniel Margulies (Université de Paris), Avram Holmes (Yale University) and Irene Droney (Child Mind Institute) are acknowledged for assistance on proofreading the manuscript. Research staff members Hao-Ming Dong, Zhu-Qing Gong, Lei Zhang and Xue-Ru Fan deserve gratitude for data extraction. Education coordinators Xu Chen, Jiang Qiu, Antao Chen, Ting-Yong Feng and Xun Liu are acknowledged for help with practical arrangements. We are also grateful to all families participating in this cohort project and all the support from schools and the communities.

Supplementary data associated with this article can be found in the online version at doi:10.1016/j.dcn.2021.101020.