Anomalous cerebral morphology of pregnant women with cleft fetuses
Zhen Li
Chunlin Li2,3,4 
Xu Zhang
Qingqing Wu- 1Department of Ultrasound, Beijing Obstetrics and Gynecology Hospital, Capital Medical University, Beijing Maternal and Child Health Care Hospital, Beijing, China
- 2Beijing Advanced Innovation Center for Big Data-Based Precision Medicine, Capital Medical University, Beijing, China
- 3School of Biomedical Engineering, Capital Medical University, Beijing, China
- 4Beijing Key Laboratory of Fundamental Research on Biomechanics in Clinical Application, Capital Medical University, Beijing, China
- 5Department of Radiology, Beijing Obstetrics and Gynecology Hospital, Capital Medical University, Beijing Maternal and Child Health Care Hospital, Beijing, China
Objective: Pregnancy leads to long-lasting changes in brain structure for healthy women; however, little is known regarding alterations in the cortical features of pregnant women with malformed fetuses. Isolated clefts of the lip and/or palate (ICL/P) are the most common congenital anomaly in the craniofacial region, which is highly gene-associated. We speculated that pregnant women carrying fetuses with ICL/P may have associated risk genes and specific brain changes during pregnancy.
Methods: In this study, we investigated T1-weighted brain magnetic resonance imaging data from 48 pregnant women: 24 women carrying fetuses with ICL/P (ICL/P group) and 24 women carrying normal fetuses (normal controls), then explored intergroup differences in gray matter volume (GMV), cortical thickness (CT) and cortical complexity (gyrification).
Results: Compared with controls, the ICL/P group had decreased total intracranial volume (TIV) than normal controls; besides, they exhibited increased GMV in the left cuneus, decreased GMV in the right superior temporal gyrus; increased CT in the left precuneus and left superior parietal gyrus, decreased CT involving parsopercularis, fusiform, middle temporal in the left hemisphere and supramarginal, precentral gyrus (PreCG) in the right hemisphere; increased gyrification in the left insula and PreCG, the left middle temporal, and the right supratemporal gyrus.
Conclusion: Pregnant women with ICL/P fetuses had brain morphology changes involving language, auditory, vision, and sensory cortex, which may be their special brain changes compared to normal pregnant women. This study may provide clues for the early detection of fetuses with ICL/P, and be vital for preconception and prenatal counseling with non-invasive methods.
Introduction
Pregnancy leads to long-lasting structural and functional adaptations in the mother’s brain for non-human animals and humans, which are necessary for the onset, maintenance, and regulation of maternal behavior (Nelander et al., 2018; Feldman et al., 2019; Pawluski et al., 2022). Structural magnetic resonance imaging (MRI) on healthy primiparous women (first-time mothers) explored an increase of the pituitary gland (Gonzalez et al., 1988), transient reductions in overall brain size in pregnant women (Oatridge et al., 2002), and widespread reduced gray matter volume (GMV) covering the right temporal lobe, precuneus, prefrontal cortex compared with nulliparous, which was associated with decreased cognitive function (Hoekzema et al., 2017). Carmona et al. (2019) found a monthly rate of volumetric reductions of 0.09 mm3 for primiparous women, accompanied by decreases in cortical thickness (CT), surface area, local gyrification index, sulcal depth, and sulcal length, as well as increases in sulcal width. Furthermore, these motherhood-related brain structural changes predicted measures of postpartum maternal attachment (Hoekzema et al., 2020). There were also studies showing no significant changes in GMV, white matter volume (WMV), or brain volume for primiparous women during this period (Zheng et al., 2018).
At present, very little is known concerning maternal brain structural changes in pregnant women with malformed fetuses. This kind of research may provide important clues for preconception and prenatal counseling with non-invasive methods. Isolated clefts of the lip and/or palate (ICL/P) are the most common congenital anomaly in the craniofacial region, which is highly gene-associated. Numerous studies explored widespread brain structural changes in the ICL/P cohorts (Adamson et al., 2014; Conrad et al., 2015, 2021; Kuhlmann et al., 2021). We speculated that pregnant women carrying fetuses with ICL/P (pregnancies with ICL/P) may have associated risk genes and specific expression, for example, specific brain changes during pregnancy. Exploring these changes will add to the early detection of fetuses with ICL/P, and will be vital for preconception and prenatal counseling.
In our previous functional MRI (fMRI) study, pregnancies with ICL/P were shown to have altered functional connectivity and topological indices within neural networks of speech and language (Li et al., 2019). Brain morphological analysis is a stable and reliable method to investigate brain structures (Schmaal et al., 2017; Zhao et al., 2019). In this study, we analyzed their brain structural data and compare GMV, CT, and cortical complexity (gyrification) between pregnancies with ICL/P and pregnant women carrying healthy fetuses (normal controls). We speculated that pregnancies with ICL/P may have specific brain cortical indexes, which may reveal the structural basis for the altered functional indices shown in our previous fMRI study. Besides, these results may provide clues for the early detection of fetuses with ICL/P, and be vital for preconception and prenatal counseling. Furthermore, they might provide clues for interaction between fetuses and their mothers during pregnancy for future study.
Subjects
This study recruited 24 pregnant women carrying fetuses with ICL/P (ICL/P group) and 24 pregnant women carrying healthy fetuses [normal controls (NC) group] in Beijing Gynecology and Obstetrics Hospital affiliated with Capital Medical University from January 2018 to December 2019. Age, educational levels, and gestation weeks (GW) of both groups were matched (Table 1). The enrollment criteria mainly included: singleton pregnancy, carrying ICL/P fetuses, or normal fetuses. All fetuses were without intrauterine growth restriction or chromosome abnormality. The key exclusion criteria for mothers were brain structural abnormalities, neurological or psychiatric disorders, complications of pregnancy, and MRI contraindications. The research protocol was approved by the Medical Research Ethics Committee of Beijing Gynecology and Obstetrics Hospital. The methods were carried out in accordance with the Declaration of Helsinki. All participants provided written informed consent after being informed of the study details.
Table 1. Demographic data of the isolated clefts of the lip and/or palate (ICL/P) and normal controls (NC).
Magnetic resonance imaging scan protocol
Brain MRI scans were performed using a 3.0-T MR scanner (Discovery MR750, GE, Milwaukee, WI, United States), 32-channel head coil. Participants were placed on their backs so that they remained relaxed during the scan. A cushion was used to limit head movement. The T1-weighted brain structure scan parameters were as follows: repetition time = 8.2 ms; echo time = 3.2 ms; flip angle = 12; acquisition matrix = 256 × 256; voxel size = 1 mm × 1 mm × 1 mm; 164 contiguous axial slices.
Structure magnetic resonance imaging preprocessing
Brain MRI data were processed with the Computational Anatomy Toolbox 12 (CAT121) and Statistical Parametric Mapping12 (SPM122). Images were processed with the following steps: converted to Neuroimaging Informatics Technology Initiative (NIFTI) files, corrected for bias–field inhomogeneities, spatially normalized (using the Diffeomorphic Anatomical Registration using Exponentiated Lie algebra (DARTEL) algorithm (Ashburner, 2007)), segmented [into gray matter (GM), white matter (WM), and cerebrospinal fluid (CSF)] (Ashburner and Friston, 2005). Then GM maps were modulated and smoothed with a Gaussian kernel of 8 mm (full-width at half maximum). CT was computed with topological correction (based on spherical harmonics) (Yotter et al., 2011) and the established novel algorithm (for extracting the cortical surface) (Dahnke et al., 2013). The local curvature-based gyrification index was calculated based on absolute mean curvature (AMC) (Luders et al., 2006). Central cortical surfaces were created for both hemispheres separately. Surface-based images were resampled and smoothed with a Gaussian kernel of 20 mm (full-width at half maximum).
Statistical analysis
Demographics data were compared using SPSS 22.0 software and p < 0.05 was considered statistically significant. Normally distributed data were expressed as mean ± SD, and tested by a two-sample t-test; non-normal distributed data were expressed as median (interquartile range, IQR), and tested by Mann–Whitney U-test non-parametric test.
The GMV maps of the two groups were compared using Matlab2013b and SPM12 with a two-sample t-test. Intergroup differences in CT and gyrification were performed using CAT12 and SPM12 with a two-sample t-test, across each hemisphere. Total intracranial volume (TIV, for GMV comparison), age and GW, and educational years were added as covariates. Statistical significance was defined as p < 0.001, uncorrected with a minimum cluster extent of 10 voxels (for GMV); p < 0.001, uncorrected or p < 0.05, peak-level Family Wise Error (FWE)-corrected with a minimum cluster extent of 20 vertices (for CT and gyrification).
Results
Demographic data
There were no statistical differences in age, GW, and educational background between the ICL/P and NC groups (Table 1).
Intergroup gray matter volume difference between the isolated clefts of the lip and/or palate and normal controls groups
Compared with controls, the ICL/P group had decreased TIV than the NC group (Table 1). Besides, they exhibited increased GMV in the left cuneus and decreased GMV in the right superior temporal gyrus (p < 0.001, uncorrected, Supplementary Figure 1 and Supplementary Table 1).
Intergroup cortical thickness difference between the isolated clefts of the lip and/or palate and normal controls groups
Two-sample t-test showed that the ICL/P group exhibited increased CT in the left precuneus and left superior parietal gyrus, decreased CT involving parsopercularis, fusiform, middle temporal in the left hemisphere and the right precentral gyrus (PreCG) (p < 0.001, uncorrected, Figure 1 and Table 2). Furthermore, they had decreased CT in the right supramarginal gyrus (SMG) (p < 0.05, peak-level FWE-corrected, Figure 2 and Table 2). CT of the intergroup differential regions of both groups was displayed in Supplementary Figure 2.
Figure 1. Altered cortical thickness (CT) and gyrification between the isolated clefts of the lip and/or palate (ICL/P) group and normal controls (NC). A two-sample t-test was used; statistical significance was p < 0.001, uncorrected. Color bars indicate T-values. The 3D inflated brain maps show the spatial location of intergroup different regions, shown on lateral or medial views.
Table 2. Intergroup cortical thickness (CT) difference between the isolated clefts of the lip and/or palate (ICL/P) and normal controls (NC).
Figure 2. Altered cortical thickness (CT) of the right hemisphere between the isolated clefts of the lip and/or palate (ICL/P) group and normal controls (NC). Statistical significance was p < 0.001, uncorrected (left) and p < 0.05, peak-level FWE-corrected (right). Color bars indicate T-values. The 3D brain maps [supplied by SPM12 (Statistical Parametric Mapping)] show the spatial location of intergroup different regions, shown on lateral views.
Intergroup gyrification difference between the isolated clefts of the lip and/or palate and normal controls groups
The ICL/P group exhibited increased gyrification in the insula, middle temporal in the left hemisphere, and superiotemporal gyrus in the right hemisphere (p < 0.001, uncorrected, Figure 1 and Table 3). Gyrification of the intergroup differential regions of both groups was displayed in Supplementary Figure 2.
Table 3. Intergroup gyrification difference between the isolated clefts of the lip and/or palate (ICL/P) and normal controls (NC).
Discussion
At present, very little is known concerning brain structural changes in pregnant women with malformed fetuses. In this study, we analyzed brain structural data and compare GMV, CT, and gyrification of pregnant women carrying fetuses with ICL/P (pregnancies with ICL/P) with normal controls. We found that pregnancies with ICL/P have brain morphology changes involving language, auditory, vision, and sensory cortex. These results might reveal the structural basis for the altered functional indices shown in our previous fMRI study; besides, these results may provide clues for the early detection of fetuses with ICL/P, and may be vital for preconception and prenatal counseling with non-invasive methods.
Compared with controls, the ICL/P group had decreased CT in regions involved in the language neural circuit (the left parsopercularis and the right SMG). The language neural circuit includes the Broca area, Wernicke’s area, SMG, angular gyrus, and the main long associated fibers connecting different language centers (Hage and Nieder, 2016; Hagoort, 2017). Broca’s area classically comprises cytoarchitectonic areas 44 (parsopercularis) and 45 (pars triangularis) in the left hemispheres, complemented by some authors by area 47 (pars orbitalis) (Ardila et al., 2017). Areas 44 and 45 on the left side of the brain are instrumental for the production, or articulation, of speech and language (Friederici and Gierhan, 2013). The parsopercularis is connected with the oral cavity and tongue movement area, it participates in language production and speech processing. The triangle part (BA45) participates in semantic processing (Turken and Dronkers, 2011). In addition, the ICL/P group had decreased CT in the right SMG (p < 0.05, peak-level FWE-corrected). SMG was shown to participate in auditory memory processing and speech decision-making (Mirman et al., 2015). Our previous work showed that for adults with cleft of lip and palate (CLP) after surgery, during articulation rehabilitation (focusing on improving articulation clarity), the junction part of the right postcentral gyrus and right SMG was involved, which indicated right SMG was a key region for articulation improving (Li et al., 2020).
Besides, we found decreased CT in the right PreCG and increased gyrification in the left PreCG in the ICL/P group. Our previous fMRI study showed that the functional connection between PreCG and amygdala in both hemispheres declined in pregnant women with ICL/P fetuses compared with controls (Li et al., 2019). Decreased CT and increased gyrification of PreCG in both hemispheres may have explored its structural basis.
In summary, pregnant women with ICL/P fetuses (ICL/P group) have brain morphology changes involving language, auditory, vision (fusiform gyrus), and sensory cortex (right PreCG), which may be brain changes in the ICL/P group compared to normal pregnant women. This may lead to differences in their language or auditory or cognition, just as findings in the ICL/P cohorts (Conrad et al., 2015, 2021; Kuhlmann et al., 2021), while at present the cognitive and language status of parents of the ICL/P cohorts stay unknown. This may be a direction for future research. Meanwhile, we speculated that the altered cerebral morphology of pregnant women with cleft fetuses may be caused by carrying CL/P risk genes. This study provided an important supplement to the limited research on maternal brain structure for the ICL/P cohorts. Brain morphology changes of the ICL/P group may give support to the theory that abnormal migration of cells caused during facial development occurs concurrently with abnormal migration of neuronal cells. Large samples and genetic studies will be needed to supply more confirmation.
The present study also had several limitations. Firstly, the sample size of this study was small, besides, most of the results were at the uncorrected statistic level and they did not survive under multiple comparisons correction. As a result, they cannot be considered reliable and need to be carefully interpreted. Secondly, this study is a cross-sectional study, without the comparison of structural data before and after pregnancy. The results of the longitudinal study will be more convincing. This will be the future direction.
Data availability statement
The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.
Ethics statement
The studies involving human participants were reviewed and approved by the Medical Research Ethics Committee of Beijing Gynecology and Obstetrics Hospital (No. 2017-KY-085-02). The patients/participants provided their written informed consent to participate in this study.
Author contributions
ZL contributed to study conception, design, data analysis, and interpretation, statistical analyses, and manuscript writing. CL contributed to conception, design, interpretation, and manuscript writing. YL, KW, and LW contributed to the design, participant recruitment, and data acquisition and interpretation. XZ and QW contributed to the conception, design, data acquisition and interpretation, and manuscript review and modification. All authors approved the final manuscript.
Funding
This work was supported by the National Natural Science Foundation of China (Nos. 81771909 and 62171300).
Acknowledgments
We thank all volunteers who participated in the study, and the staff of Beijing Obstetrics and Gynecology Hospital, Capital Medical University, Beijing Maternal and Child Health Care Hospital, Department of Ultrasound, and Department of Radiology, China, for their selfless and valuable assistance.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fnhum.2022.959710/full#supplementary-material
Supplementary Figure 1 | Altered gray matter volume (GMV) between the isolated clefts of the lip and/or palate (ICL/P) group and normal controls (NC). A two-sample t-test was used; statistical significance was p < 0.001, uncorrected. Color bars indicate T-values. The brain maps show the spatial location of intergroup different regions, shown on a cross or coronal views.
Supplementary Figure 2 | Cortical thickness (CT) and gyrification of the intergroup differential regions of the isolated clefts of the lip and/or palate (ICL/P) group and normal controls (NC). The statistical level of the intergroup differential regions was p < 0.001, uncorrected. (L)/(R) means left/right hemisphere.
Footnotes
- ^ https://neuro-jena.github.io/cat//index.html#DOWNLOAD
- ^ https://www.fil.ion.ucl.ac.uk/spm/software/spm12/
References
Adamson, C. L., Anderson, V. A., Nopoulos, P., Seal, M. L., and Da Costa, A. C. (2014). Regional brain morphometric characteristics of nonsyndromic cleft lip and palate. Dev. Neurosci. 36, 490–498. doi: 10.1159/000365389
Ardila, A., Bernal, B., and Rosselli, M. (2017). Should Broca’s area include Brodmann area 47? Psicothema 29, 73–77.
Ashburner, J. (2007). A fast diffeomorphic image registration algorithm. Neuroimage 38, 95–113. doi: 10.1016/j.neuroimage.2007.07.007
Ashburner, J., and Friston, K. J. (2005). Unified segmentation. Neuroimage. 26, 839–851. doi: 10.1016/j.neuroimage.2005.02.018
Carmona, S., Martinez-Garcia, M., Paternina-Die, M., Barba-Muller, E., Wierenga, L. M., Aleman-Gomez, Y., et al. (2019). Pregnancy and adolescence entail similar neuroanatomical adaptations: A comparative analysis of cerebral morphometric changes. Hum. Brain Mapp. 40, 2143–2152. doi: 10.1002/hbm.24513
Conrad, A. L., Kuhlmann, E., van der Plas, E., and Axelson, E. (2021). Brain structure and neural activity related to reading in boys with isolated oral clefts. Child Neuropsychol. 27, 621–640. doi: 10.1080/09297049.2021.1879765
Conrad, A. L., Richman, L., and Nopoulos, P. (2015). Reading achievement in boys with non-syndromic cleft palate only: Relationship to neuropsychological skill and neurocircuitry. Dev. Neuropsychol. 40, 395–406. doi: 10.1080/87565641.2016.1142991
Dahnke, R., Yotter, R. A., and Gaser, C. (2013). Cortical thickness and central surface estimation. Neuroimage. 65, 336–348. doi: 10.1016/j.neuroimage.2012.09.050
Desikan, R. S., Ségonne, F., Fischl, B., Quinn, B. T., Dickerson, B. C., Blacker, D., et al. (2006). An automated labeling system for subdividing the human cerebral cortex on MRI scans into gyral based regions of interest. Neuroimage. 31, 968–980. doi: 10.1016/j.neuroimage.2006.01.021
Feldman, R., Braun, K., and Champagne, F. A. (2019). The neural mechanisms and consequences of paternal caregiving. Nat. Rev. Neurosci. 20, 205–224. doi: 10.1038/s41583-019-0124-6
Friederici, A. D., and Gierhan, S. M. (2013). The language network. Curr. Opin. Neurobiol. 23, 250–254. doi: 10.1016/j.conb.2012.10.002
Gonzalez, J. G., Elizondo, G., Saldivar, D., Nanez, H., Todd, L. E., and Villarreal, J. Z. (1988). Pituitary gland growth during normal pregnancy: An in vivo study using magnetic resonance imaging. Am. J. Med. 85, 217–220. doi: 10.1016/s0002-9343(88)80346-2
Hage, S. R., and Nieder, A. (2016). Dual neural network model for the evolution of speech and language. Trends Neurosci. 39, 813–829. doi: 10.1016/j.tins.2016.10.006
Hagoort, P. (2017). The core and beyond in the language-ready brain. Neurosci. Biobehav. Rev. 81, 194–204. doi: 10.1016/j.neubiorev.2017.01.048
Hoekzema, E., Barba-Muller, E., Pozzobon, C., Picado, M., Lucco, F., Garcia-Garcia, D., et al. (2017). Pregnancy leads to long-lasting changes in human brain structure. Nat. Neurosci. 20, 287–296. doi: 10.1038/nn.4458
Hoekzema, E., Tamnes, C. K., Berns, P., Barba-Muller, E., Pozzobon, C., Picado, M., et al. (2020). Becoming a mother entails anatomical changes in the ventral striatum of the human brain that facilitate its responsiveness to offspring cues. Psychoneuroendocrinology 112:104507. doi: 10.1016/j.psyneuen.2019.104507
Kuhlmann, E., van der Plas, E., Axelson, E., and Conrad, A. L. (2021). Brain developmental trajectories in children and young adults with isolated cleft lip and/or cleft palate. Dev. Neuropsychol. 46, 314–326. doi: 10.1080/87565641.2021.1946691
Li, Z., Li, C., Liang, Y., Wang, K., Zhang, W., Chen, R., et al. (2019). Altered functional connectivity and brain network property in pregnant women with cleft fetuses. Front. Psychol. 10:2235. doi: 10.3389/fpsyg.2019.02235
Li, Z., Zhang, W., Li, C., Wang, M., Wang, S., Chen, R., et al. (2020). Articulation rehabilitation induces cortical plasticity in adults with non-syndromic cleft lip and palate. Aging 12, 13147–13159. doi: 10.18632/aging.103402
Luders, E., Thompson, P. M., Narr, K. L., Toga, A. W., Jancke, L., and Gaser, C. (2006). A curvature-based approach to estimate local gyrification on the cortical surface. Neuroimage. 29, 1224–1230. doi: 10.1016/j.neuroimage.2005.08.049
Mirman, D., Chen, Q., Zhang, Y., Wang, Z., Faseyitan, O. K., Coslett, H. B., et al. (2015). Schwartz, neural organization of spoken language revealed by lesion-symptom mapping. Nat. Commun. 6:6762. doi: 10.1038/ncomms7762
Nelander, M., Wikstrom, A. K., Weis, J., Bergman, L., Larsson, A., Sundstrom-Poromaa, I., et al. (2018). Cerebral osmolytes and plasma osmolality in pregnancy and preeclampsia: A proton magnetic resonance spectroscopy study. Am. J. Hypertens 31, 847–853. doi: 10.1093/ajh/hpy019
Oatridge, A., Holdcroft, A., Saeed, N., Hajnal, J. V., Puri, B. K., Fusi, L., et al. (2002). Change in brain size during and after pregnancy: Study in healthy women and women with preeclampsia. AJNR Am. J. Neuroradiol. 23, 19–26.
Pawluski, J. L., Hoekzema, E., Leuner, B., and Lonstein, J. S. (2022). Less can be more: Fine tuning the maternal brain. Neurosci. Biobehav Rev. 133:104475. doi: 10.1016/j.neubiorev.2021.11.045
Schmaal, L., Hibar, D. P., Samann, P. G., Hall, G. B., Baune, B. T., Jahanshad, N., et al. (2017). Cortical abnormalities in adults and adolescents with major depression based on brain scans from 20 cohorts worldwide in the ENIGMA major depressive disorder working group. Mol. Psychiatry 22, 900–909. doi: 10.1038/mp.2016.60
Turken, A. U., and Dronkers, N. F. (2011). The neural architecture of the language comprehension network: Converging evidence from lesion and connectivity analyses. Front. Syst. Neurosci. 5:1. doi: 10.3389/fnsys.2011.00001
Yotter, R. A., Dahnke, R., Thompson, P. M., and Gaser, C. (2011). Topological correction of brain surface meshes using spherical harmonics. Hum. Brain Mapp. 32, 1109–1124. doi: 10.1002/hbm.21095
Zhao, L., Matloff, W., Ning, K., Kim, H., Dinov, I. D., and Toga, A. W. (2019). Age-related differences in brain morphology and the modifiers in middle-aged and older adults. Cereb. cortex 29, 4169–4193. doi: 10.1093/cercor/bhy300
Keywords: pregnant women, isolated clefts of the lip and/or palate (ICL/P), cortical thickness (CT), gray matter volume (GMV), magnetic resonance imaging (MRI), gyrification
Citation: Li Z, Li C, Liang Y, Wang K, Wang L, Zhang X and Wu Q (2022) Anomalous cerebral morphology of pregnant women with cleft fetuses. Front. Hum. Neurosci. 16:959710. doi: 10.3389/fnhum.2022.959710
Received: 02 June 2022; Accepted: 15 August 2022;
Published: 07 September 2022.
Edited by:
Wen Qin, Tianjin Medical University General Hospital, ChinaReviewed by:
Zheng Wei Min, Capital Medical University, ChinaChuanjun Zhuo, Tianjin Anding Hospital, China
Copyright © 2022 Li, Li, Liang, Wang, Wang, Zhang and Wu. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Xu Zhang, zhangxu@ccmu.edu.cn; Qingqing Wu, qingqingwu@ccmu.edu.cn