# a Chinese multi-modal neuroimaging data release for increasing diversity of human brain mapping

Peng Gao, Hao-Ming Dong, Si-Man Liu, Xue-Ru Fan, Chao Jiang, Yin-Shan

Wang, Daniel Margulies, Hai-Fang Li, Xi-Nian Zuo

► To cite this version:

Peng Gao, Hao-Ming Dong, Si-Man Liu, Xue-Ru Fan, Chao Jiang, et al., a Chinese multi-modal neuroimaging data release for increasing diversity of human brain mapping. Scientific Data, 2022, 9 (1), pp.286. 10.1038/s41597-022-01413-3. http://hal-04307121

HAL Id: hal-04307121

https://hal.science/hal-04307121v1

Submitted on 25 Nov 2023

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>
scientific **data**

<figure><img src="image_2.png" /></figure>

Check for updates

OPEN

DATA DESCRIPTOR

# A Chinese multi-modal neuroimaging data release for increasing diversity of human brain mapping

Peng Gao<sup>1,9</sup>, Hao-Ming Dong<sup>2,3,9</sup>, Si-Man Liu<sup>4</sup>, Xue-Ru Fan<sup>2</sup>, Chao Jiang<sup>1,5</sup>, Yin-Shan Wang<sup>2</sup>, Daniel Margulies<sup>1,6</sup>, Hai-Fang Li<sup>1,7</sup> & Xi-Nian Zuo<sup>1,2,3,4,7,8</sup>

The big-data use is becoming a standard practice in the neuroimaging field through data-sharing initiatives. It is important for the community to realize that such open science effort must protect personal, especially facial information when raw neuroimaging data are shared. An ideal tool for the face anonymization should not disturb subsequent brain tissue extraction and further morphological measurements. Using the high-resolution head images from magnetic resonance imaging (MRI) of 215 healthy Chinese, we discovered and validated a template effect on the face anonymization. Improved facial anonymization was achieved when the Chinese head templates but not the Western templates were applied to obscure the faces of Chinese brain images. This finding has critical implications for international brain imaging data-sharing. To facilitate the further investigation of potential culture-related impacts on and increase diversity of data-sharing for the human brain mapping, we released the 215 Chinese multi-modal MRI data into a database for imaging Chinese young brains, namely 'I See your Brains (ISYB)', to the public via the Science Data Bank (https://doi.org/10.11922/sciencedb.00740).

## Background & Summary

With the worldwide increase in the acquisition of neuroimaging data using magnetic resonance imaging (MRI), the use of big data is becoming a standard practice in the field through international data-sharing initiatives<fnref n="1" />⁻⁴. Data sharing can bring many benefits to neuroscience research and increase sample sizes, which allow for greater precision and improve the statistical power to detect smaller effects, although smaller effects may be associated with either smaller biological effects or potentially confounding factors⁵. This has accelerated many fields into their new stages or even entirely new field, such as the developmental population neuroscience⁶,⁷ or statistical connectomics⁸. On the one hand, it can facilitate researchers to conduct test-retest studies when more neuroimaging public databases are available and improve the reproducibility of research results⁹,¹⁰. Data reuse can also make peer review better⁹,¹¹. On the other hand, it encourages cooperation and communication between institutions and organizations, and thus creates a win-win situation for both researchers giving and receiving data¹².

As the concept of sharing data continues to grow, more and more public databases are emerging to support researchers such as the Human Connectome Project (HCP)¹³ and its derivatives across human lifespan¹⁴ (e.g., the Adolescent Brain Cognition Development, ABCD¹⁵), UK Biobank¹⁶, the Autism Brain Imaging Data Exchange (ABIDE)¹⁷,¹⁸, the 1000 Functional Connectomes Project¹⁹. These efforts on big-data sharing has propelled the cognitive neurology field forward (see more data projects from our recent reviews²⁰,²¹). However, most of the existing databases are based on Caucasian populations. This preventing us from identifying ethnic or culture-related differences in neuroscience associations. Moreover, China has the largest population in the world and there is still a need to expand and improve the public database of the Chinese brain mapping

<footnote n="1">¹College of Information and Computer, Taiyuan University of Technology, Taiyuan, 030024, China. ²State Key Laboratory of Cognitive Neuroscience and Learning, Beijing Normal University, Beijing, 100875, China. ³National Basic Science Data Center, Beijing, 100109, China. ⁴Institute of Psychology, Chinese Academy of Sciences, Beijing, 100101, China. ⁵School of Psychology, Capital Normal University, Beijing, 100048, China. ⁶Centre National de la Recherche Scientifique, Frontlab, Brain and Spinal Cord Institute, Paris, UMR 7225, France. ⁷Developmental Population Neuroscience Research Center, IDG/McGovern Institute for Brain Research, Beijing Normal University, Beijing, 100875, China. ⁸Key Laboratory of Brain and Education, School of Education Science, Nanning Normal University, Nanning, 530001, China. ⁹These authors contributed equally: Peng Gao, Hao-Ming Dong. ℮e-mail: lihaifang@tyut.edu.cn; xinian.zuo@bnu.edu.cn</footnote>
(few exceptions, e.g., the Chinese Color Nest Project²⁰, the REST-meta-MDD project²², the Chinese Human Connectome Project²³ and the Consortium for Reliability and Reproducibility⁹).

Much like the genetics field²⁴, data-sharing faces many challenges in the neuroscience field. One of the challenges is to ensure that personal information, especially facial information, remains protected when raw MRI data are shared. This is especially challenging for high-resolution T1-weighted anatomical images, in which facial information has been fully embedded. The higher the image resolution is, the clearer the facial information becomes. It has been possible to restore subject's facial information through the combination of existing facial recognition and reconstruction technologies. As a result, personal information can easily be disclosed. Therefore, we must first remove or obscure the facial features in brain images, especially in high-resolution structural MRI images. In response to this, the Health Insurance Portability and Accountability Act (HIPAA) of the United States specifically stipulates that the subject's private photos or equivalent images must be removed from the medical images to protect their privacy²⁵. The European Union has general data protection, which requires research institutions to strengthen the privacy protection of personal data. Therefore, data protection should be considered and implemented during the design phase of scientific studies²⁴,²⁶.

In order to ensure participant's private security, all the aforementioned face masking algorithms are involved in taking a standard brain template as references, so as to blur the facial information in the template space. However, to date, previous studies have not dealt with whether a different facial mask will be generated due to a different ethnic template used during anonymization. This may reflect the situation that the existing public brain imaging databases rarely contain non-western samples. In this data descriptor, we release the'I See Your Brains' database for imaging Chinese young brains and demonstrate the template effects on face masking using head templates of different races. Specifically, ethnically diverse brain templates are employed to anonymize Chinese multimodal MRI data independently. The discrepancy between Chinese brain images anonymized using different templates of race is quantitatively examined. We hypothesized that the use of different (Chinese versus Western) templates will cause significant differences in the face masking performance. Use of the corresponding ethnic template will achieve a better performance in face anonymization.

## Methods

**Participant recruitment.** The database, Imaging Chinese Young Brains (namely, “I See Your Brains”, ISYB), is generated by the Institute of Psychology, Chinese Academy of Sciences (IPCAS), which is one contributing site of the Chinese imaging genetics cohort (CHIMGEN)²⁷ to enhance cross-ethnic and cross-geographic brain research. It is the largest prospective neuroimaging genetic cohort for Chinese Han adults with natural and socioeconomic measurements obtained from remote sensing. In accordance with CHIMGEN, all participants in ISYB were recruited by advertisements posted in colleges and communities. Participants were excluded if they met any of the following criteria: regular smoker, pregnancy, abnormal color discrimination, a history of alcohol or drug abuse, currently any medication, MRI contraindications, neuropsychiatric or severe somatic disorder and sedative-hypnotic medication within a month or any medication for major neuropsychiatric disorders. The ISYB study was approved by the IPCAS ethical committee and written informed consent was obtained from each participant.

**Image acquisition.** ISYB contains multi-modal neuroimaging data from 241 right-handed Chinese healthy volunteers. Each participant received MRI scans of the brain images at a 3.0 Tesla MRI scanner (GE MR750) at the IPCAS MRI Research Center including: (1) the high-resolution T1-weighted structural MRI (sMRI, matrix = 256 × 256, number of slices = 188, field of view = 256 × 256 mm, repetition time = 8.16 ms, echo time = 3.18 ms, flip angle = 12°, inversion time = 450 ms), (2) resting-state functional MRI (rfMRI, matrix = 64 × 64, number of slices = 36, field of view = 220 × 220 mm, repetition time = 2000ms, echo time = 30 ms, flip angle = 90°), (3) diffusion tensor MRI (dMRI, matrix = 128 × 128, number of slices = 50, field of view = 256 × 256 mm, repetition time = 6000 ms, echo time = 65 ms, flip angle = 90°, 64 diffusion-sensitisation directions at b = 1000 with b = 0), and (4) arterial spin labeling MRI (aMRI, matrix = 128 × 128, number of slices = 50, field of view = 240 × 240 mm, repetition time = 5046 ms, echo time = 11.09 ms, flip angle = 111°, inversion time = 2025 ms).

**Quality assurance.** All participants were included in the first-step quality assessments of the neuroimaging data. Their anatomical T1 sMRI images were visually inspected to ensure no substantial head motion and structural abnormalities. Two participants were excluded for excessive head motion. In the second-step quality assessments, the T1 sMRI and rfMRI image quality of ISYB datasets were examined using the **MRIQC** toolkit (https://github.com/poldracklab/mriqc), which incorporated a series of quantitative metrics as in the following list. More details of these metrics can be found in the previous work²⁸,²⁹.

* Spatial Metrics (sMRI, rfMRI):

1. Signal-to-Noise Ratio (SNR)³⁰: The mean value of gray matter divided by the standard deviation (SD) of the air.

2. Foreground to Background Energy Ratio (FBER): Mean energy of image values (i.e., mean of squares) within the head relative to outside the head.

3. Contrast-to-Noise Ratio (CNR) (only sMRI)³⁰: Calculated as the mean of the gray matter values minus the mean of the white matter values, divided by the standard deviation of the air values.

4. Entropy Focus Criteria (EFC)³¹: Shannon's entropy is used to summarize the principal directions distribution.

5. Smoothness of Voxels³²: The full-width half maximum (FWHM) of the spatial distribution of image intensities.
<figure><img src="image_3.png" /></figure>

**Fig. 1** Quality assessments on the ISYB sMRI data using both bar (the left column) and scatter (the right column) plots of the following quality metrics: (a) Signal-to-Noise Ratio (SNR) (b) Foreground to Background Energy Ratio (FBER) (c) Contrast-to-Noise Ratio (CNR) (d) Entropy Focus Criteria (EFC) (e) Full-Width Half Maximum (FWHM).

<figure><img src="image_4.png" /></figure>

**Fig. 2** Quality assessments on the ISYB rfMRI data using both bar (the left column) and scatter (the right column) plots of the following quality metrics: (a) Signal-to-Noise Ratio (SNR) (b) Foreground to Background Energy Ratio (FBER) (c) Entropy Focus Criteria (EFC) (d) Full-Width Half Maximum (FWHM), (e) Mean Framewise Displacement (FD_MEAN), (f) Derivative of time series and root-mean-square VARIance over voxels (DVARS), (g) Global Correlation (GCOR).

* Temporal Metrics (rfMRI):

1. Mean Framewise Displacement (FD)³³: A measure of subject head motion, which compares the motion between the current and previous volumes. This is calculated by summing the absolute value of displacement changes in the x, y and z directions and rotational changes about those three axes. The rotational changes are given distance values based on the changes across the surface of a 50 mm radius sphere.

2. Standardized DVARS³³: The spatial standard deviation of the temporal derivative of the data (D refer-ring to temporal derivative of time series, VARS referring to root-mean-square variance over voxels), normalized by the temporal standard deviation and temporal auto-correlation.

3. Global Correlation (GCOR)³⁴: The average of the entire brain correlation matrix, which is computed as the brain-wide average time series correlation over all possible combinations of voxels.

The mean and SD were calculated for each of the eight (five spatial metrics and three temporal metrics) to remove extreme values or outliers which is 3-SD above or below the mean. Taking the effects of head move-ment in rfMRI into consideration, the data with average FD greater than 0.20 mm³³ were further excluded. We note that there were 26 participants excluded by the quality assessment procedure, remaining 215 participants included for the subsequent analyses (156 females; 18–30 years old, mean = 22.55 yrs, SD = 2.68 yrs).

The assessments on the results of the quality assurance of structural MRI images and functional MRI images are summarized in Figs. 1 and 2, respectively. Details of individual quality reports are presented in the supple-mentary figures (SFig. 1 for sMRI and SFig. 2 for rfMRI).

**Standard brain templates.** Taking into consideration on the global scope of the open science commu-nity and data-sharing initiatives, it is critical to realize that the face masking procedure initially uses Western head template as reference. Extensive research has shown that the race factor played a major role when com-paring the group differences in brain morphology metrics. For example, previous studies have indicated signifi-cant morphological differences between Chinese and Western brains in size³⁵, shape³⁶, and anatomical volume³⁷.
Caucasian brains are significantly longer but present a smaller width and a decreased height compared with Chinese brains²³,³⁸. Thus, taking templates of different races as reference during registration exerts a significant impact on the results of the registration, which may lead to systematic differences in subsequent morphological measurements³⁹. Selecting ethnicity-matched templates to analyze structural MRI data may be preferable, which could improve the accuracy of the registration and obtain robust morphological results⁴⁰. This is attributed to the differences in the brain structure between China and the West population, which cause some systematic bias in the registration between individual images and template images and reduce the accuracy in morphological measurements. Such bias has been reproduced in diverse samples across imaging protocols and lifespan studies⁴¹⁻⁴⁴.

Ethnicity-specific brain templates are discrepant regarding the diversity in the brain structures between the East and the West. To test the template effects on face masking, we chose two Chinese brain templates (CN200 and Chinese2020): (1) CN200 is a head template assembled by MRI T1 images using the 3 T Siemens Prisma scanner from 200 healthy Chinese Han participants from the Chinese Human Connectome Project²³; (2) Chinese2020 is a head template generated based on the MRI T1 images under multiple 1.5 T magnets collected from more than 2,000 healthy participants across multiple centers and regions⁴⁵. The MNI152 head template was used as the Western template, which is the default head template as in the face masking toolkit.

**Face anonymization.** In order to effectively blur the part of the face in images and protect privacy as much as possible, researchers have developed several toolkits to remove facial information. For example, Quickshear segments brain images into facial features and brain, and configures the voxels of the facial features to 0, thus resulting in identifiable facial features cut effectively²⁵. FreeSurfer determines the location of facial information and creates a brain template after manually separating brain tissue from non-brain tissue. In the subsequent images, the voxel position of non-brain tissue is obtained through the established linear transformation relationship and the voxel value of the facial information position is eliminated⁴⁶. Although these two methods effectively remove facial privacy information, they also lose some areas of the original image at the same time (e.g., intracranial volume and cerebrospinal fluid volume), which may cause the problem of lacking important anatomical features when calculating certain structural metrics. To address this privacy concern, Milchenko and colleagues have developed a deface toolkit, namely face masking, which extracts the anatomical surface of the face in the volume space and then obscures it⁴⁷. The recognition rate of the image is reduced by changing the image resolution, leading to ideal anonymization of the data without disturbing reliability and validity of subsequent tissue extraction and measurements of brain morphology⁴⁸⁻⁵¹. This procedure has become a built-in function of the Extensible Neuroimaging Archive Toolkit⁵² for facial anonymization and been incorporated into the Human Connectome Project⁵³. For these reasons mentioned above, when data in the ISYB database were anonymized for facial information, we adopted the CN200 Chinese template. The purpose of which is to remove more privacy and obtain a better anonymization effect for facial information in Chinese brain data. The results and discussion of data anonymization based on different ethnic brain templates are included in Usage Notes session. The ISYB data are processed by the face masking toolkit.

## Data Records

As supports of the open sciences, we have released the multimodal neuroimaging data of the 215 participants in the ISYB database for both research and education in human brain mapping as well as culture association studies. Specifically, as demonstrated in the previous sections, we first applied the face masking pipeline based on CN200 brain template to anonymize the MRI data by blurring facial information while all the participants agreed to share their anonymized data to the public. We then organized the data according to the Brain Imaging Data Structure (BIDS: https://bids.neuroimaging.io)⁵⁴ and employed the **MRIQC** toolkit for monitoring the quality of all structural and functional images. This toolkit generated comprehensive reports on the data quality, which were presented as in the Supplementary Fig. 1 and Supplementary Fig. 2 at group level. All the BIDS-formatted data and related files on the QC reports have been publicly shared through the Science Data Bank (SciDB) at the National Basic Science Data Center (https://doi.org/10.11922/sciencedb.00740)⁵⁵.

The **ISYB_DATA** folder contains 215 folders labeled with participant ID numbers (e.g., sub-0002). Each participant's folder contains subfolders labeled *func* containing the rfMRI images (`.nii.gz`), *anat* containing defaced T1-weighted sMRI images (`.nii.gz`), *perf* containing aMRI images (`.nii.gz`), *dwi* containing dMRI images (`.nii.gz`, `.bval`, `.bvec`). The **ISYB_QC** folder contains QC reports generated by **MRIQC** toolkit for each participant. All QC results were integrated into webpages, and linked in an overview page on group analysis. Of note, the current ISYB release contained the demographic information (age and sex) but did not include any behavioral and clinical measurements.

## Technical Validation

As in the Consortium for Reliability and Reproducibility (CoRR)⁹, we demonstrate the utility of the released ISYB multi-modal MRI data. All the anonymized and quality-controlled data were preprocessed using the Connectome Computation System (CCS)⁵⁶, https://github.com/zuoxinian/CCS) with the most recent updates⁵⁷. This pipeline integrates multiple analytical software packages to achieve imaging processing of multi-modal MRI data. Considering the advantages of surface-based functional brain mapping⁵⁸, we reconstructed cortical surface models to generate NIFTI file for structural and functional metrics on the cortical surface. The structural image went through the following preprocessing steps: (1) spatially adaptive non-local means denoising, (2) rough inhomogeneity correction, (3) spatial normalization into the MNI standard brain space, (4) inhomogeneity correction, (5) intensity normalization, (6) brain extraction by non-local intracranial cavity extraction (NICE), and (7) gray and white matter segmentation, surface reconstruction. The rfMRI image preprocessing included (1) dropping off the first 5 EPI volumes, (2) removing and interpolating temporal spikes, (3) correcting acquisition timing among image slices and head motion among image volumes, (4) normalizing the 4D global mean
intensity to 10,000, (5) regressing out head motion artifacts and other spurious noise by using ICA-AROMA⁵⁹,⁶⁰, and (6) removing linear and quadratic trends from the rfMRI signals to mitigate the scanner-related influences⁵⁶.

The preprocessed NIFTI images were converted into the GIFTI format using the Ciftify toolbox⁶¹. The sMRI images were further preprocessed using the custom pipeline to generate highly refined cortical surface meshes which were spatially normalized to the MNI space and re-sampled to have 32k vertices, i.e., the fsaverage_LR32 cortical surface. Individual dMRI and aMRI images were preprocessed in individual volume space and then transformed onto the fsaverage_LR32 space. Individual surface-mapped rfMRI signals were then brought into the register across participants using a multi-modal surface matching algorithm to the and vectorized into the CIFTI format. This maps each surface vertex to an index in a vector by the Ciftify toolbox. A group-level surface mask was finally established by including every vertex where all the 215 participants obtained the full rfMRI signals. Based on the surface-based data, we derived two sets of brain imaging metrics using CCS for both structural and functional characteristics, respectively. Specifically, following derivatives were calculated:

* Cortical Thickness (CT)⁶²: The average distance between the gray/white boundary and the pial surface within each ROI.

* Surface Area (SA)⁶²: The sum of the areas of each tessellation falling within each ROI.

* Gray Matter Volume (GMV): The measure of the density of brain cell body in a particular region.

* Local Gyrification Index (LGI)⁶³: The ratio of a regional surface area for the pial surface to a smooth surface contour estimated to wraps around the pial surface.

* Fractional Anisotropy (FA)⁶⁴: An index for the amount of diffusion asymmetry within a voxel, defined in terms of its eigenvalues.

* Mean Diffusivity (MD)⁶⁴: The mean amount of diffusion in each of the principal directions calculated in the tensor.

* Regional Homogeneity (ReHo)⁵⁸: The synchronicity of a voxel's time series and that of its nearest neighbors based on Kendall's coefficient of concordance to measure the local brain functional homogeneity.

* Voxel-Mirrored Homotopic Connectivity (VMHC)⁶⁵: The functional connectivity between a pair of geometrically symmetric, inter-hemispheric voxels.

* Amplitude of Low-Frequency Fluctuations (ALFF)⁶⁶: The total power in the low frequency range (0.01–0.1 Hz) of an fMRI image, normalized by the total power across all frequencies measured in that same image.

* Cerebral Blood Flow (CBF)⁶⁷: The rate of delivery of arterial blood to the capillary bed in brain tissue and is quantified in milliliters of blood per 100 g of brain tissue per minute.

We note that each metric was calculated within the group-level surface mask. The CT metric was generated and extracted from the structural outcomes of **ciftify** in the fsaverage_LR32k space. The LGI metric was produced by the FreeSurfer scripts (https://surfer.nmr.mgh.harvard.edu/fwiki/LGI). The SA, ReHo, ALFF and VMHC metrics were calculated by the CCS scripts from the sMRI and rfMRI outcomes of **ciftify**. CCS also provided a computational workflow for dMRI data to derive both FA and MD for gray matter and white matter voxels in individual volume spaces. The results were named as FA-WM, MD-WM, FA-GM and MD-GM, respectively. The **bbregister** command in FreeSurfer was then applied to the native individual maps of CBF, SA, FA and MD to register these individual volume maps onto the fsaverage surface template. The **wb_command -metric-resample** command in the HCP workbench was used to map the native individual CBF, SA, FA and MD metrics to fs_LR and resample them to the 32k space. The vertex-wise group-level maps across all qualified participants are shown on the fsLR_32k surface as in Fig. 3 (mean) and Fig. 4 (SD).

## Usage Notes

We encourage other labs to use this dataset in publication under the requirement of citing this article and contact us for additional data sharing and cooperation. The user of the ISYB database should acknowledge the contributions of the original authors and research lab, and properly cite this article. We note that the impact of the brain template on the facial anonymization procedure is not negligible. In this section, to evaluate the anonymization of different racial brain template, we used high-resolution sMRI images from the database. The ISYB database contains the 215 Chinese multi-modal MRI images⁵⁵ and is part of the CHIMGEN consortium (see more details in the consortium paper²⁷). We split the ISYB dataset into two subsets for replicability validation (Discovery sample: ISYB-1, n = 89; Validation sample: ISYB-2, n = 126).

**Implementation of analytic framework.** As illustrated in Fig. 5, we implemented an analytical framework for evaluation of the template effects on face masking. We noted that the Western brain template (i.e., MNI152 template) was employed as the default brain template in the face masking toolkit. We thus substituted the MNI152 template with the CN200 and Chinese2020 templates to establish the corresponding relationship (i.e., registration) between the two Chinese templates and the individual data space. The corresponding relationship and different ethnic templates were paired and integrated into the face masking process. The three different pipelines were then used to implement the anonymous processing of the MRI data with the different templates. As for the evaluation of the face anonymization performance, we first subtracted the anonymized image from the original image to get the facial information, which was removed from the original face. We then identified the facial areas blurred by the face masking toolkit and created a mask for each blurred surface based on the individual MRI image.

Normalized mutual information (NMI) was used to quantify the efficiency of performing the face-masking process with different templates. NMI is widely used to evaluate the similarity between two sets of information. It serves as an important indicator in information theory and can be used to evaluate the spatial similarity or
<figure><img src="image_5.png" /><figcaption>**Fig. 3** The group mean maps of the ISYB multi-modal neuroimaging data on the fsLR_32k cortical surfaces. Given each vertex on the surface, the mean values were calculated across all the 215 participants for the 12 metrics based on the following three MRI modalities: (a) sMRI (CT, SA, GMV, LGI), (b) dMRI (FA-WM, MD-WM, FA-GM, MD-GM), (c) rfMRI and aMRI metrics (CBF, ALFF, ReHo and VMHC).</figcaption></figure>

<figure><img src="image_6.png" /><figcaption>**Fig. 4** The group variability maps of the ISYB multi-modal neuroimaging data on the fsLR_32k cortical surfaces. Given each vertex on the surface, the standard deviation values were calculated across all the 215 participants for the 12 metrics based on the following three MRI modalities: (a) sMRI (CT, SA, GMV, LGI), (b) dMRI (FA-WM, MD-WM, FA-GM, MD-GM), (c) rfMRI and aMRI metrics (CBF, ALFF, ReHo and VMHC).</figcaption></figure>

information shared between images. Its value ranges from 0 to 1. A higher NMI value indicates more spatial information shared between a pair of brain images. In the present study, a higher NMI value reveals more facial information picked up by the face-masking pipeline, and thus more efficient in removing private facial information from the original MRI head image. It means that the privacy information of the participant is blurred more efficiently. Therefore, we hypothesize that the face-masking pipelines with a race-appropriate head template would be more effective and recommendable for the anonymization of the facial information.

When the original brain images from Chinese (ISYB-1 and ISYB-2) passed through the face masking procedure with three different head templates (CN200, Chinese2020 and MNI152), six sets of facial masked images under the
<figure><img src="image_7.png" /><figcaption>**Fig. 5** The analytic framework of the evaluation of template effects on face anonymization. (a) The face masking toolkit takes raw individual unmasked MRI data with the registration from the unmasked head image to the head template. (b) The head templates include the default MNI152 and other two Chinese (CN200 and Chinese2020) templates. (c) During the face masking, a face mask is generated and used for restoring the face removed by the face masking pipeline. (d) The normalized mutual information between the removed face and the original face is calculated.</figcaption></figure>

<table><caption>**Table 1.** Statistics of the post-hoc tests on the face anonymization with the different templates.</caption><thead><tr><td>Dataset</td><td>Template_1</td><td>Template_2</td><td>Difference</td><td>Standard Error</td><td>p-value</td></tr></thead><tbody><tr><td>ISYB</td><td>CN200</td><td>Chinese2020</td><td>1.48 × 10<sup>-3</sup></td><td>5.03 × 10<sup>-4</sup></td><td>9.11 × 10<sup>-3</sup></td></tr><tr><td>ISYB</td><td>CN200</td><td>MNI152</td><td>3.41 × 10<sup>-3</sup></td><td>3.59 × 10<sup>-4</sup></td><td>9.56 × 10<sup>-10</sup></td></tr><tr><td>ISYB</td><td>Chinese2020</td><td>MNI152</td><td>1.93 × 10<sup>-3</sup></td><td>4.67 × 10<sup>-4</sup></td><td>1.12 × 10<sup>-4</sup></td></tr></tbody></table>

2 × 3 different conditions were generated. The NMI values between the facial information removed under different conditions and the corresponding original data were calculated for each participant. A two-way analysis of variance with repeated measures (database versus template) was performed on the NMI values with the head template as the within-participant factor for testing the interaction effect between template and database on the face masking pipeline. The statistical model revealed a highly significant template effect (DoF = 4, F = 7.08, p = 1.42 × 10⁻⁵). We further carried out a set of post-hoc paired two-sample tests to demonstrate the directions of the effect detected. Multiple comparisons of the estimated marginal means based on the three templates in the repeated measures model were also applied. The details of the results are presented in the Table 1. Specifically, for ISYB dataset, both CN200 and Chinese2020 templates exhibited better performances of face masking than the MNI152 template while the 3.0 T MRI derived CN200 template was better than the 1.5 T MRI derived Chinese2020 template. As illustrated in Fig. 6, such effects were highly replicable across the two data subsets (ISYB-1 versus ISYB-2).

**Guidance and implication.** With the continuous development of neuroimaging, the distinction in the brain structure across different races can become more refined. The comparison of templates of different races provides a basis for more in-depth research on the distinctions of separate races. However, many challenges have also emerged, especially when the private information is easily leaked through images. At the same time, due to the continuous development and maturity of deep learning technology, combined with big data, the entire face can be reconstructed and refabricated using a small portion of facial information. This forces us to pay special attention to protection of subjects' private information, especially facial information. Masking the face of the data by using the face masking toolkit effectively blur the subjects' faces, thereby protecting privacy. Meanwhile, it protects the integrity of the brain structure, which makes the structure MRI easier to interface with mainstream brain image analysis software. Follow-up analysis is more convenient without degenerating image quality. Since skull shape differs by racial and ethnic groups in length and width⁶⁸, systematic biases may be introduced when images of different races are registered to the same template in the face masking toolkit. This kind of registration difference based on data and templates which are not in the same race has been emphasized in preceding articles, which is inevitable due to genetic factors⁴¹. This study assumed that using different human brain templates to anonymize facial data would yield different results⁶⁹. We selected 215 data from Chinese database and three brain templates from two races, and used ethnicity-matched and ethnicity-mismatched templates to blur the facial information of the data, explored the impact of different ethnic templates on data anonymization.
<figure><img src="image_8.png" /><figcaption>**Fig. 6** The replicable template effects on face anonymization. Normalized mutual information between the face image from original data and the face image removed by the anonymization are visualized as violin plots with paired lines under the three different template (CN200, MNI152, Chinese2020) conditions. This template effect is replicated across three datasets: (a) ISYB-1, (b) ISYB-2, (c) ISYB. \*p\*<0.05, \*\*\*p\*<0.0001.</figcaption></figure>

After face masking was performed on Chinese data respectively using brain templates from two races, matching the data according to the ethnic template, quantifying the results according to the index of normalized mutual information, it can be concluded that: (1) There existed a significant difference in the anonymization results when the results based on the two Chinese templates were compared with that based on the Western template. The different ethnicity of the templates lead to the differences in the results when using different ethnic templates to blur face information in structural brain images by the face masking toolbox. Therefore, we believe that such variability is attributed to the race factor. (2) Choosing a template that matches the race of the data will achieve better results in comparison to with a mismatched template when using different templates to anonymize data. In other words, the template selected during anonymization should correspond to the race of the data so that better private information anonymity results can be achieved. In the process of masking the face, the surface mask algorithm proposed by Budin and colleagues was used and performed<sup>70</sup>. The face area was identified and positioned, and then smoothed by using a blur algorithm to obtain a contaminated face image. However, in the face mask algorithm, the data need to be registered on the specific template, and then the brain tissue position and skull boundary information could be found through the binary image, thereby blurring the part between the brain and skull. The registration effect varies from race to race, resulting in different positions of facial features and causing disparate phenomena based on different templates when using templates of different races as a reference to registration.

The effect of facial anonymization varies by brain templates to a large extent, the results demonstrate that there exist a difference in the effect of facial masking based on two Chinese templates, the difference may attributed to the fact that the two Chinese templates are constructed using data obtained with different field strengths (1.5T and 3.0T). What we want to emphasis is that the anonymization effect is determined by the ethnicity of the templates. The brain structure of individuals from different ethnic group is quite different. Although different templates can be used to anonymize structure MRI data, the use of ethnicity-matched templates will lead to the best anonymization effect. This also demonstrated that the accuracy of brain image data processing could be enhanced by using the ethnicity-appropriate templates. Improved facial anonymization, conserving less information, has been demonstrated when the Chinese head template is applied to obscure the faces of data acquired in a Chinese population. The customized pipeline using a native Chinese head template is thus recommended for the anonymization of sharing Chinese neuroimaging data. International data-sharing and the relevant large-scale team work on brain imaging such as building brain charts for the human lifespan (see a recent review<sup>21</sup> and progress<sup>71</sup>) must consider this effect carefully. Its impact on the generalizability of brain-mind association studies<sup>10,72</sup> is also an interesting topic in the future.

## Code availability

All the scripts and brain templates involved for processing is publicly available as part of the recent CCS updates<sup>57</sup> and can be visited at GitHub (https://github.com/zuoxinian/CCS/tree/master/projects/isybdemo).

Received: 4 November 2021; Accepted: 24 May 2022;

Published online: 09 June 2022
## References

1. Zuo, X. N. Mapping the miswired connectome in autism spectrum disorder. *Journal of the American Academy of Child & Adolescent Psychiatry* **59**, 348–349 (2020).

2. Milham, M. P. et al. Assessment of the impact of shared brain imaging data on the scientific literature. *Nature Communications* **9**, 2818 (2018).

3. Nichols, T. E. et al. Best practices in data analysis and sharing in neuroimaging using MRI. *Nature Neuroscience* **20**, 299–303 (2017).

4. Qu, H., Lei, H. & Fang, X. Big data and the brain: Peeking at the future. *Genomics, Proteomics & Bioinformatics* **17**, 333–336 (2019).

5. Smith, S. M. & Nichols, T. E. Statistical challenges in “big data” human neuroimaging. *Neuron* **97**, 263–268 (2018).

6. Zuo, X.-N. et al. Developmental population neuroscience: emerging from ICHBD. *Science Bulletin* **63**, 331–332 (2018).

7. Fair, D. A., Dosenbach, N. U., Moore, A. H., Satterthwaite, T. D. & Milham, M. P. Developmental cognitive neuroscience in the era of networks and big data: Strengths, weaknesses, opportunities, and threats. *Annual Review of Developmental Psychology* **3**, 249–275 (2021).

8. Chung, J. et al. Statistical connectomics. *Annual Review of Statistics and Its Application* **8**, 463–492 (2021).

9. Zuo, X. N. et al. An open science resource for establishing reliability and reproducibility in functional connectomics. *Scientific Data* **1**, 1–13 (2014).

10. Marek, S. et al. Reproducible brain-wide association studies require thousands of individuals. *Nature* **603**, 654–660 (2022).

11. Poldrack, R. A. & Gorgolewski, K. J. Making big data open: data sharing in neuroimaging. *Nature Neuroscience* **17**, 1510–1517 (2014).

12. Ascoli, G. A., Maraver, P., Nanda, S., Polavaram, S. & Armañanzas, R. Win-win data sharing in neuroscience. *Nature Methods* **14**, 112–116 (2017).

13. Van Essen, D. C. et al. The WU-Minn human connectome project: an overview. *Neuroimage* **80**, 62–79 (2013).

14. Harms, M. P. et al. Extending the human connectome project across ages: Imaging protocols for the lifespan development and aging projects. *Neuroimage* **183**, 972–984 (2018).

15. Casey, B. et al. The Adolescent Brain Cognitive Development (ABCD) study: Imaging acquisition across 21 sites. *Developmental Cognitive Neuroscience* **32**, 43–54 (2018).

16. Littlejohns, T. J. et al. The UK Biobank imaging enhancement of 100,000 participants: rationale, data collection, management and future directions. *Nature Communications* **11**, 2624 (2020).

17. Di Martino, A. et al. The autism brain imaging data exchange: towards a large-scale evaluation of the intrinsic brain architecture in autism. *Molecular Psychiatry* **19**, 659–667 (2014).

18. Martino, A. D. et al. Enhancing studies of the connectome in autism using the autism brain imaging data exchange ii. *Scientific Data* **4**, 170010 (2017).

19. Biswal, B. B. et al. Toward discovery science of human brain function. *Proceedings of the National Academy of Sciences* **107**, 4734–4739 (2010).

20. Liu, S. et al. Chinese Color Nest Project: An accelerated longitudinal brain-mind cohort. *Developmental Cognitive Neuroscience* **52**, 101020 (2021).

21. Chen, L. Z., Holmes, A. J., Zuo, X. N. & Dong, Q. Neuroimaging brain growth charts: A road to mental health. *Psychoradiology* **1**, 272–286 (2021).

22. Yan, C. G. et al. Reduced default mode network functional connectivity in patients with recurrent major depressive disorder. *Proceedings of the National Academy of Sciences* **116**, 9078–9083 (2019).

23. Yang, G. Y. et al. Sample sizes and population differences in brain template construction. *Neuroimage* **206**, 116318 (2020).

24. Byrd, J. B., Greene, A. C., Prasad, D. V., Jiang, X. & Greene, C. S. Responsible, practical genomic data sharing that accelerates research. *Nature Reviews Genetics* **21**, 615–629 (2020).

25. Schimke, N. & Hale, J. Quickshear defacing for neuroimages. *2nd USENIX Conference on Health Security and Privacy* (2011).

26. White, T., Blok, E. & Calhoun, V. D. Data sharing and privacy issues in neuroimaging research: Opportunities, obstacles, challenges, and monsters under the bed. *Human Brain Mapping* **43**, 278–291 (2022).

27. Xu, Q. et al. CHIMGEN: a chinese imaging genetics cohort to enhance cross-ethnic and cross-geographic brain research. *Molecular Psychiatry* **25**, 517–529 (2020).

28. Esteban, O. et al. MRIQC: Advancing the automatic prediction of image quality in MRI from unseen sites. *PLoS One* **12**, e0184661 (2017).

29. Esteban, O. et al. Crowdsourced MRI quality metrics and expert quality annotations for training of humans and machines. *Scientific Data* **6**, 30 (2019).

30. Magnotta, V. A., Friedman, L. & Birn, F. Measurement of signal-to-noise and contrast-to-noise in the fBIRN multicenter imaging study. *Journal of Digital Imaging* **19**, 140–147 (2006).

31. Atkinson, D., Hill, D. L., Stoyle, P. N., Summers, P. E. & Keevil, S. F. Automatic correction of motion artifacts in magnetic resonance images using an entropy focus criterion. *IEEE Transactions on Medical Imaging* **16**, 903–910 (1997).

32. Friedman, L. et al. Test–retest and between-site reliability in a multicenter fMRI study. *Human Brain Mapping* **29**, 958–972 (2008).

33. Power, J. D., Barnes, K. A., Snyder, A. Z., Schlaggar, B. L. & Petersen, S. E. Spurious but systematic correlations in functional connectivity MRI networks arise from subject motion. *Neuroimage* **59**, 2142–2154 (2012).

34. Saad, Z. S. et al. Correcting brain-wide correlation differences in resting-state fMRI. *Brain Connectivity* **3**, 339–352 (2013).

35. Bai, J. et al. Population differences in brain morphology and microstructure among Chinese, Malay, and Indian neonates. *PLoS One* **7**, e47816 (2012).

36. Lee, J. S. et al. Development of Korean standard brain templates. *Journal of Korean Medical Science* **20**, 483–488 (2005).

37. Xie, W. Z., Richards, J. E., Lei, D., Lee, K. & Gong, Q. Comparison of the brain development trajectory between Chinese and US children and adolescents. *Frontiers in Systems Neuroscience* **8**, 249 (2015).

38. Rao, N. P. et al. Population differences in brain morphology: Need for population specific brain template. *Psychiatry Research: Neuroimaging* **265**, 1–8 (2017).

39. Tang, Y. C. et al. The construction of a Chinese MRI brain atlas: a morphometric comparison study between chinese and caucasian cohorts. *Neuroimage* **51**, 33–41 (2010).

40. Xie, W. Z. et al. The construction of MRI brain/head templates for Chinese children from 7 to 16 years of age. *Developmental Cognitive Neuroscience* **15**, 94–105 (2015).

41. Tang, Y. C. et al. Brain structure differences between Chinese and Caucasian cohorts: A comprehensive morphometry study. *Human Brain Mapping* **39**, 2147–2155 (2018).

42. Richards, J. E. & Xie, W. Brains for all the ages: structural neurodevelopment in infants and children from a life-span perspective. In *Advances in child development and behavior*, **48**, 1–52 (Elsevier, 2015).

43. Zhao, T. D. et al. Unbiased age-specific structural brain atlases for Chinese pediatric population. *Neuroimage* **189**, 55–70 (2019).

44. Dong, H. M. et al. Charting brain growth in tandem with brain templates at school age. *Science Bulletin* **65**, 1924–1934 (2020).

45. Liang, P. P. et al. Construction of brain atlases based on a multi-center MRI dataset of 2020 Chinese adults. *Scientific Reports* **5**, 1–7 (2015).

46. Bischoff, G. A. et al. A technique for the deidentification of structural brain MR images. *Human Brain Mapping* **28**, 892–903 (2007).

47. Milchenko, M. & Marcus, D. Obscuring surface anatomy in volumetric imaging data. *Neuroinformatics* **11**, 65–75 (2013).

48. Zuo, X. N., Xu, T. & Milham, M. P. Harnessing reliability for neuroscience research. *Nature Human Behaviour* **3**, 768–771 (2019).

49. Xing, X. X. et al. The anatomy of reliability: a must read for future human brain mapping. *Science Bulletin* **63**, 1606–1607 (2018).

50. Marcus, D. S. et al. Open Access Series of Imaging Studies (OASIS): cross-sectional mri data in young, middle aged, nondemented, and demented older adults. *Journal of Cognitive Neuroscience* **19**, 1498–1507 (2007).

51. Lerch, J. P. et al. Studying neuroanatomy using MRI. *Nature Neuroscience* **20**, 314–326 (2017).
52. Marcus, D. S., Olsen, T. R., Ramaratnam, M. & Buckner, R. L. The extensible neuroimaging archive toolkit. *Neuroinformatics* 5, 11–33 (2007).

53. Glasser, M. F. et al. The human connectome project’s neuroimaging approach. *Nature Neuroscience* **19**, 1175–1187 (2016).

54. Gorgolewski, K. et al. The brain imaging data structure, a format for organizing and describing outputs of neuroimaging experiments. *Scientific Data* **3**, 160044 (2016).

55. Gao, P., Dong, H. M., Wang, Y. S., Yu, C. S. & Zuo, X. N. *Imaging Chinese Young Brains (I See Your Brain)* https://doi.org/10.11922/sciencedb.00740 (2021).

56. Xu, T., Yang, Z., Jiang, L. L., Xing, X. X. & Zuo, X. N. A Connectome Computation System for discovery science of brain. *Science Bulletin* **60**, 86–95 (2015).

57. Xing, X. X., Xu, T., Jiang, C., Wang, Y. S. & Zuo, X. N. Connectome Computation System: 2015–2021 updates. *Science Bulletin* **67**, 448–451 (2022).

58. Zuo, X. N. et al. Toward reliable characterization of functional homogeneity in the human brain: Preprocessing, scan duration, imaging resolution and computational space. *Neuroimage* **65**, 374–386 (2013).

59. Pruim, R. H. R. et al. ICA-AROMA: A robust ICA-based strategy for removing motion artifacts from fMRI data. *Neuroimage* **112**, 267–277 (2015).

60. Pruim, R. H. R., Mennes, M., Buitelaar, J. K. & Beckmann, C. F. Evaluation of ICA-AROMA and alternative strategies for motion artifact removal in resting state fMRI. *Neuroimage* **112**, 278–287 (2015).

61. Dickie, E. W. et al. Ciftify: A framework for surface-based analysis of legacy MR acquisitions. *Neuroimage* **197**, 818–826 (2019).

62. Panizzon, M. S. et al. Distinct genetic influences on cortical surface area and cortical thickness. *Cerebral Cortex* **19**, 2728–2735 (2009).

63. Lyu, I., Kim, S. H., Girault, J. B., Gilmore, J. H. & Styner, M. A. A cortical shape-adaptive approach to local gyration index. *Medical Image Analysis* **48**, 244–258 (2018).

64. Alexander, A. L., Lee, J. E., Lazar, M. & Field, A. S. Diffusion Tensor Imaging of the Brain. *Neurotherapeutics* **4**, 316–329 (2007).

65. Zuo, X. N. et al. Growing together and growing apart: Regional and sex differences in the lifespan developmental trajectories of functional homotopy. *Journal of Neuroscience* **30**, 15034–15043 (2010).

66. Zuo, X. N. et al. The oscillating brain: Complex and reliable. *Neuroimage* **49**, 1432–1445 (2010).

67. Borogovac, A. & Asllani, I. Arterial Spin Labeling (ASL) fMRI: Advantages, theoretical constraints and experimental challenges in neurosciences. *International Journal of Biomedical Imaging* **2012**, 818456 (2012).

68. Zulkifli, N. S. A. & Kamal, N. F. Does common standard brain template standardize for all brains regardless the different of age, gender and culture? *Journal of Life Science and Technologies* 150–153 (2013).

69. Bhalerao, G. V. et al. Construction of population-specific Indian MRI brain template: Morphometric comparison with Chinese and Caucasian templates. *Asian Journal of Psychiatry* **35**, 93–100 (2018).

70. Budin, F., Zeng, D., Ghosh, A. & Bullitt, E. Preventing facial recognition when rendering MR images of the head in three dimensions. *Medical Image Analysis* **12**, 229–239 (2008).

71. Bethlehem, R. A. I. et al. Brain charts for the human lifespan. *Nature* **604**, 525–533 (2022).

72. Yarkoni, T. The generalizability crisis. *Behavioral and Brain Sciences* **45**, e1 (2022).

## Acknowledgements

This work was supported in part by the National Natural Science Foundation of China (61976150), the Major Project of National Social Science Foundation of China (20&ZD296), the National Basic Science Data Center “Chinese Data-sharing Warehouse for In-vivo Imaging Brain” (NBSDC-DB-15), the Child Brain-Mind Development Cohort Study in China Brain Initiative (SQ2021AAA010024) and the Start-up Funds for Leading Talents at Beijing Normal University.

## Author contributions

X.N. Zuo, H.F. Li: study design and conceptualisation; data acquisition and formal data analysis; manuscript preparation; data curation; management of online repositories. P. Gao, H.M. Dong: data acquisition and formal data analysis; manuscript preparation; data curation; management of online repositories. All authors edited and revised the manuscript.

## Competing interests

The authors declare no competing interests.

## Additional information

**Supplementary information** The online version contains supplementary material available at https://doi.org/10.1038/s41597-022-01413-3.

**Correspondence** and requests for materials should be addressed to H.-FL. or X.-N.Z.

**Reprints and permissions information** is available at www.nature.com/reprints.

**Publisher’s note** Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

<figure><img src="image_9.png" /></figure>

**Open Access** This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.

© The Author(s) 2022