Tiannan Shao1†
Weiyan Wang1†
Gangrui Hei2
Ye Yang1
Yujun Long1
Xiaoyi Wang1
Jingmei Xiao1
Yuyan Huang1
Xueqin Song2
Xijia Xu3
Shuzhan Gao3
Jing Huang1
Ying Wang1
Jingping Zhao1
Renrong Wu1*- 1Department of Psychiatry, and National Clinical Research Center for Mental Disorders, The Second Xiangya Hospital of Central South University, Changsha, China
- 2Department of Psychiatry, The First Affiliated Hospital of Zhengzhou University, Zhengzhou, China
- 3Department of Psychiatry, Affiliated Nanjing Brain Hospital, Nanjing Medical University, Nanjing, China
Background: Cognitive subtypes of schizophrenia may exhibit different neurobiological characteristics. This study aimed to reveal the underlying neurobiological features between cognitive subtypes in the early course of schizophrenia (ECS). According to prior studies, we hypothesized to identify 2–4 distinct cognitive subtypes. We further hypothesized that the subtype with relatively poorer cognitive function might have lower brain spontaneous neural activity than the subtype with relatively better cognitive function.
Method: Cognitive function was assessed by the MATRICS Consensus Cognitive Battery (MCCB). Resting-state functional magnetic resonance imaging scanning was conducted for each individual. There were 155 ECS individuals and 97 healthy controls (HCs) included in the subsequent analysis. Latent profile analysis (LPA) was used to identify the cognitive subtypes in ECS individuals, and amplitude of low-frequency fluctuations (ALFFs) was used to measure brain spontaneous neural activity in ECS individuals and HCs.
Results: LPA identified two cognitive subtypes in ECS individuals, containing a severely impaired subtype (SI, n = 63) and a moderately impaired subtype (MI, n = 92). Compared to HCs, ECS individuals exhibited significantly increased ALFF in the left caudate and bilateral thalamus and decreased ALFF in the bilateral medial prefrontal cortex and bilateral posterior cingulate cortex/precuneus (PCC/PCu). In ECS cognitive subtypes, SI showed significantly higher ALFF in the left precentral gyrus (PreCG) and lower ALFF in the left PCC/PCu than MI. Furthermore, ALFFs of left PreCG were negatively correlated with several MCCB cognitive domains in ECS individuals, while ALFF of left PCC/PCu presented opposite correlations.
Conclusion: Our findings suggest that differences in the brain spontaneous neural activity of PreCG and PCC/PCu might be the potential neurobiological features of the cognitive subtypes in ECS, which may deepen our understanding of the role of PreCG and PCC/PCu in the pathogenesis of cognitive impairment in schizophrenia.
Introduction
Cognitive impairment is a core component of schizophrenia, which exists not only in the acute phase but also in the ultra-high risk stage of psychosis and persists during the clinical remission of psychiatric symptoms (Keefe et al., 2006; Bora and Murray, 2014). In schizophrenia, cognitive impairment is manifested in multiple domains, including verbal and visual learning, working memory, processing speed, problem solving, attention, and executive function (Nuechterlein et al., 2004; Sheffield et al., 2018), which could contribute to poor functional outcomes (Green, 2016). Although there is a generalized cognitive deficit in schizophrenia, not all individuals exhibit the same pattern of cognitive deficits, i.e., degrees of impairment in different cognitive dimensions from different sample groups are manifested differently (Carruthers et al., 2019). These heterogeneous patterns of cognitive impairment make it difficult to understand the pathophysiology of schizophrenia.
Accumulating evidence has successfully classified individuals with psychiatric disorders with similar cognitive characteristics and compared differences between distinct cognitive features. Participants could be classified according to their current cognitive function scores based on artificial classification (Ammari et al., 2010; Ortiz-Gil et al., 2011), or could be classified by data-driven approaches (Lim et al., 2020; Smucny et al., 2020). Instead of appearing to be limited to dividing participants into relatively intact or cognitively impaired subtypes by prior defined classification criteria, the data-driven approach could provide valuable insight into multiple cognitive subtypes that exist in psychiatric disorders (Carruthers et al., 2019).
Abnormal cognitive function also has relevant neuroimaging features. Growing attention has been paid to identifying the neurobiological changes of cognitive subtypes in schizophrenia through neuroimaging (Van Rheenen et al., 2018; Lewandowski et al., 2019). Several studies using data-driven methods have illustrated the different patterns of brain structure or resting-state functional connectivity among cognitive subtypes in individuals with schizophrenia. For brain structure, near-normal cognitive groups and impaired cognitive groups of individuals with schizophrenia or schizoaffective disorders were found at varying levels of cortex thickness (Cobia et al., 2011), or different gray and white matter volumes (Wexler et al., 2009; Van Rheenen et al., 2018). For brain function connectivity networks, three cognitive subtypes classified by k-means clustering method from 67 individuals with schizophrenia spectrum disorders showed unique hyper- or hypo-connectivity in specific functional networks (Rodriguez et al., 2019). However, previous neuroimaging studies have used structural measurements and functional connectivity approaches to investigate the neurobiological features between cognitive subtypes in schizophrenia; to the best of our knowledge, there was no study investigating differences in brain spontaneous neural activity between cognitive subtypes of ECS yet.
The amplitude of low-frequency fluctuations (ALFFs) is often used to characterize brain spontaneous neural activity (Zang et al., 2007). Several studies have investigated the association between ALFF and cognitive function in schizophrenia (Zhou et al., 2014; Wang et al., 2019). Therefore, the purpose of this study was to determine the differences in brain spontaneous neuronal activity measured by ALFF between cognitive subtypes of individuals with ECS. We compared the cognitive subtypes derived from latent profile analysis (LPA) with healthy controls (HCs) to assess the condition of cognitive impairment (Miettunen et al., 2016). Based on two recent systematic reviews on cognitive subtypes in schizophrenia, we expected to find 2–4 distinct cognitive subtypes: a severely impaired subgroup, 1–2 moderately impaired subgroups, and a relatively intact cognitive subgroup (Carruthers et al., 2019, 2021). We further hypothesized that cognitive subtypes might exhibit different brain spontaneous neuronal activities, and the brain spontaneous neuronal activity may be lower in the subtype with poorer cognitive function than in the subtype with better cognitive function.
Methods
Participants
We recruited ECS individuals and HCs from three clinical medical centers, the Second Xiangya Hospital of Central South University (Center 1), the Affiliated Nanjing Brain Hospital of Nanjing Medical University (Center 2), and the First Affiliated Hospital of Zhengzhou University (Center 3). All recruited ECS individuals met the criteria of the Diagnostic and Statistical Manual of Mental Disorders, Fourth Edition (DSM-IV) based on a Structured Clinical Interview for DSM-IV Axis-I Disorder (First et al., 1997), which were assessed by two well-trained psychiatrists in each center.
For ECS individuals, the further inclusion criteria were as follows: (1) aged 16–60; (2) experienced the first episode of psychiatric symptoms or illness duration within 3 years but in the acute phase currently; (3) Positive and Negative Syndrome Scale (PANSS) total score >60. The exclusion criteria were: (1) comorbidity of physical diseases, or other psychiatric disorders that met the DSM-IV criteria; (2) comorbidity of substance abuse or addiction; (3) unable to complete cognitive function tests and magnetic resonance imaging (MRI) examinations; (4) pregnant or lactating women. In addition, HCs with first-degree relatives that had any psychiatric disorders were also excluded.
Written informed consents were obtained from all participants or first-degree relatives of ECS individuals before participating in the study. This study was approved by the Ethics Committee of the Second Xiangya Hospital of Central South University, the Ethics Committee of the Affiliated Nanjing Brain Hospital of Nanjing Medical University, and the Ethics Committee of the First Affiliated Hospital of Zhengzhou University.
Our study totally enrolled 181 ECS individuals and 107 HCs. After neuroimaging preprocessing, 155 ECS individuals and 97 HCs were included in the subsequent analysis (subjects with head motion exceeding 2 mm or head rotation exceeding 2° were excluded). Among 155 ECS individuals, 119 (76.8%) individuals received risperidone, 16 (10.3%) olanzapine, 11 (7.1%) amisulpride, seven (4.5%) aripiprazole, and two (1.3%) paliperidone. The doses of antipsychotic drugs were equivalent to chlorpromazine dose by the defined daily doses (DDDs) method (Leucht et al., 2016). The remaining subjects of each center were as follows: 34 ECS individuals and 18 HCs (Center 1); 38 ECS individuals and 36 HCs (Center 2); 83 ECS individuals and 43 HCs (Center 3). Demographic and clinical data of ECS individuals and HCs from three centers are shown in Supplementary Table 1.
Cognitive assessment
We used Measurement and Treatment Research to Improve Cognition in Schizophrenia Consensus Cognitive Battery (MCCB) to evaluate cognitive function (Kern et al., 2008; Nuechterlein et al., 2008). The MCCB was widely used and translated into Chinese version (Shi et al., 2015). It includes nine tasks across seven cognitive domains: (1) speed of processing (Trail Making Test, part A; Brief Assessment of Cognition in Schizophrenia, Symbol Coding; Category Fluency Test, Animal); (2) attention and vigilance (Continuous Performance Test-Identical Pairs); (3) working memory (Wechsler Memory Scale, spatial span); (4) verbal learning (Hopkins Verbal Learning Test-Revised); (5) visual learning (Brief Visuospatial Memory Test); (6) reasoning and problem solving (Neuropsychological Assessment Battery, mazes test); (7) social cognition (Mayer–Salovey–Caruso Emotional Intelligence Test, managing emotions test). It takes 1–1.5 h for each subject to finish MCCB. Raw scores were converted to scale scores, then to Chinese-normalized T scores. T scores of seven cognitive domains and composite scores (the average T score of nine tasks) were calculated. In this study, all ECS individuals and HCs completed the MCCB.
Neuroimaging data acquisition
The thorough description of this section is described in the Supplementary material.
Neuroimaging data preprocessing
Image preprocessing was performed by the Data Processing & Analysis of Brain Imaging toolbox (DPABI, V4.2, http://rfmri.org/dpabi) running on MATLAB software (The MathWorks, Inc., Natick, MA, USA; Yan et al., 2016). We removed the first 10 scanning volumes in order to stabilize the magnetic resonance signal and reduce the impact of subjects not adapting to the scanning environment. Slice timing and head motion were corrected for each subject, and those whose head motion exceeded 2 mm or head rotation exceeded 2° were excluded. Afterward, images were spatially normalized to the standard Montreal Neurological Institute template by using warping parameters estimated from T1 images with a resampling standard voxel size of 3 mm × 3 mm × 3 mm. We used a 6 mm full-width at half-maximum Gaussian kernel to spatially smooth images and performed linear detrending. Then the nuisance signals were regressed out, including head motion effects (Friston 24-parameter model; Friston et al., 1996), white matter, and cerebrospinal fluid. Finally, band-pass filtering (0.01–0.08 Hz) was applied for the time series of each voxel to remove the effects of very-low-frequency drifts and high-frequency noise (Zang et al., 2007). After that, the time series of each voxel for each subject was transformed into the frequency domain under a fast Fourier transformation way in order to get the power spectrum. To standardize the ALFF values, the ALFF of each voxel would be converted into Z scores by subtracting the global mean and then dividing the global standard deviation.
Multi-site effect harmonization
Before statistical analysis, we used the ComBat Harmonization method (http://github.com/Jfortin1/ComBatHarmonization) for preprocessed data to eliminate the inter-site effects (Fortin et al., 2017). This widely used method could effectively remove unwanted variation introduced by the site, and increase statistical power (Fortin et al., 2018; Radua et al., 2020). In addition, group, age, sex, and education level were protected during the removal of inter-site effects.
Latent profile analysis
Latent profile analysis is used to classify individuals into heterogeneous subtypes based on latent variable models. It could explain the associations between the observed continuous indicator variables by regressing the continuous indicator variables onto a set of one or more latent class variables (Miettunen et al., 2016). LPA is a model-based approach, and thus has fewer prerequisites for application, more reasonable clustering criteria and result testing, and less arbitrariness than traditional clustering methods (e.g., k-means; Brusco et al., 2017; Schreiber, 2017). The flexibility of LPA makes it adaptable to the heterogeneous study of complex psychiatric and psychological phenomena with effective classification of cognitive subtypes (Lim et al., 2020; Smucny et al., 2020; De Meo et al., 2021).
Latent profile analysis was conducted by Mplus version 7.11 to identify potential homogenous subtypes of ECS individuals based on cognitive performance in seven MCCB cognitive domains (Muthén and Muthén, 2015). The number of classes was determined from an examination of models fit statistics rather than hypothesized. These model fit indices included log-likelihood ratio (LLR; Woolf, 1957), Akaike's information criteria (AIC; Akaike, 1987), Bayesian information criteria (BIC; Schwarz, 1978), sample-size adjusted BIC (ABIC; Sclove, 1987), and entropy (Celeux and Soromenho, 1996). Lo-Mendell-Rubin (LMR) tests and bootstrapped likelihood ratio tests (BLRTs) were also conducted to evaluate the significance of model improvement between n and n – 1 number of classes (McLachlan, 1978; Lo et al., 2001). A total of five models were estimated specifying from 1 to 5 latent classes.
Statistical analysis
Once the potential subtypes were identified, group differences in demographic data, clinical data, PANSS scores, and MCCB scores were analyzed by SPSS version 22.0 (IBM, Armonk, NY, USA) by using one-way analysis of variance (ANOVA), two-sample t-test, or chi-squared test. Furthermore, eta-squared (η2) was used to calculate the effect size of comparisons of each MCCB cognitive domain (Cohen, 1988). Post-hoc comparisons were carried out by Bonferroni correction if ANOVA showed significant differences between subtypes. Generally, p-values of < 0.05 were accepted as statistically significant.
ALFF analysis was conducted in the Statistical Parametric Mapping 12 toolbox (SPM12, https://www.fil.ion.ucl.ac.uk/spm/software/spm12). Two-sample t-test was designed for the comparison of ALFF maps between ECS individuals and HCs, and in ECS cognitive subtypes, with age, sex, education, and mean frame-wise displacement Jenkinson as covariates. In addition, PANSS total score was also controlled for subtypes comparison to investigate whether the result was consistent after subtracting out the effect of symptom severity. Multiple comparisons were corrected using the cluster-wise family-wise error (FWE) rates correction (cluster-wise FWE p < 0.05) with a combined individual voxel threshold of p < 0.001. Significant brain regions with discrepant ALFF between ECS cognitive subtypes were regarded as regions of interest (ROI) to extract ALFF values for subsequent correlation analysis in SPSS. To further explore the specific associations between ALFF in ROIs and MCCB scores, their correlations coefficient between ALFF values in ROIs and MCCB scores in ECS and HCs were calculated, respectively. Pearson's r was used to calculate the effect size of correlations between ALFF in ROIs and MCCB cognitive scores (Cohen, 1988). The significant level of correlations was corrected by the false discovery rate (FDR) at q < 0.05.
Results
LPA results based on MCCB cognitive domains
According to the results of model estimation, the 2-class solution presented a better fit than the 1-class solution. Though the other solutions suggested better fits than the 2-class solution based on AIC, BIC, and ABIC, there was no significant improvement over the 2-class solution according to the p-value of LMR and the p-value of BLRT (Table 1, Supplementary Figure 1). The 2-class solution classified 40.6% of the ECS individuals into the severely impaired subtype (SI, n = 63) and 59.4% into the moderately impaired subtype (MI, n = 92), respectively.
Table 1. Model estimations of latent profile analysis based on MCCB cognitive domains.
Demographic and clinical characteristics
Demographic and clinical characteristics of two ECS cognitive subtypes and HCs are shown in Table 2. The sex composition of SI significantly differed from MI. Besides, education levels were similar in both subtypes but significantly lower than the HCs. For symptom severity, PANSS scores of SI were significantly higher than MI except for positive score. Significant pairwise differences in MCCB domains were exhibited among the three groups, while differences in the working memory domain and reasoning/problem-solving domain between MI and HCs were not significant (Table 2). In addition, the MCCB performance of HC were taken as the norm to standardize the cognitive score of ECS into Z scores (see Figure 1 and Supplementary Table 2).
Table 2. Demographic and clinical data of ECS cognitive subtypes and HCs.
Figure 1. MCCB cognitive scores in ECS cognitive subtypes (SI and MI) and HCs. Z scores of ECS cognitive subtypes are standardized against means and SDs of HC (mean = 0, SD = 1). Abbreviations: ECS, early course schizophrenia; SI, severely impaired subtype; MI, moderately impaired subtype; HC, healthy control; MCCB, MATRICS Consensus Cognitive Battery; SD, standard deviation.
ALFF differences between groups
Results of the independent two-sample t-test (cluster-wise FWE p < 0.05) and cluster information are shown in Table 3 and Figure 2. Compared with HCs, ECS exhibited higher ALFF in the left caudate and bilateral thalamus, while lower ALFF was observed in the bilateral medial prefrontal cortex (MPFC) and bilateral posterior cingulate cortex/precuneus (PCC/PCu). For the subtypes comparison, ALFF of the left precentral gyrus (PreCG) was significantly higher in SI. On the contrary, the ALFF of left PCC/PCu in SI was significantly lower than MI. To control for the effect of PANSS and medication effects (chlorpromazine dose equivalence), we conducted the additional analysis with PANSS total score and medication effects as covariates. This result was similar to the previous result without controlling for the two covariates (see Supplementary material).
Table 3. Brain regions with ALFF differences in ECS and HCs and in ECS cognitive subtypes.
Figure 2. ALFF differences in ECS, HCs, and ECS cognitive subtypes (SI and MI). (A) Abnormal ALFF in ECS relative to HCs. (B) Increased ALFF in left PreCG and decreased ALFF in left PCC/PCu in SI relative to MI. Results were cluster-wise FWE corrected. Abbreviations: ALFF, amplitude of low-frequency fluctuations; ECS, early course schizophrenia; SI, severely impaired subtype; MI, moderately impaired subtype; HC, healthy control; MPFC, medial prefrontal cortex; PCC/PCu, posterior cingulate cortex/precuneus; PreCG, precentral gyrus; L, left; B, bilateral; FWE, family-wise error.
Correlations between ALFF of ROIs and cognition
For ECS individuals, there were significantly positive correlations between ALFF values of left PCC/PCu with MCCB speed of processing, attention/vigilance, verbal learning, reasoning/problem solving, and composite score (Figure 3). The effect sizes of the correlations were small to moderate. Moreover, ALFF in left PreCG was significantly negatively correlated with MCCB seven cognitive domains and composite score (Figure 3). Importantly, some correlations of ALFF in left PreCG had moderate to large effect sizes. No significant correlations were found in HCs.
Figure 3. Correlations between ALFF values of two ROIs and MCCB cognitive scores in ECS. q, p-values corrected by FDR correction. Abbreviations: ALFF, amplitude of low-frequency fluctuations; ECS, early course schizophrenia; MCCB, MATRICS Consensus Cognitive Battery; PCC/PCu, posterior cingulate cortex/precuneus; PreCG, precentral gyrus; ROI, region of interest; FDR, false discovery rate; L, left.
Discussion
The present study revealed the relationship between cognitive impairment and neurobiology in ECS individuals grouped according to their different patterns of cognitive deficit using a data-driven approach. We reported two patterns of cognitive impairment in ECS individuals. Compared to HCs, ECS individuals showed significantly increased ALFF in the left caudate and bilateral thalamus and decreased ALFF in bilateral MPFC and bilateral PCC/PCu. Interestingly, we first found that ALFF of left PreCG and left PCC/PCu were different in the two identified impaired cognitive subtypes. In addition, ALFF of left PreCG and left PCC/PCu exhibited significant correlations with MCCB cognitive domains in ECS individuals. These findings suggest that the different spontaneous neural activities of PreCG and PCC/PCu at resting-state may be the potential neurobiological features of cognitive impairment subtypes in ECS.
Consistent with the grouping results of previous studies (Morar et al., 2011; Green et al., 2013), our LPA recognized two different patterns of cognitive impairment in ECS individuals. SI exhibited the worst performance in all MCCB cognitive scores, which manifested wide deficits in whole cognitive domains, ranging from 0.7 to 3.5 standard deviations below HCs (Figure 1 and Supplementary Table 2). MI showed an intermediate deficient speed of processing, attention/vigilance, verbal learning, visual learning, and social cognition, but maintained intact working memory and reasoning/problem solving, within ~1.2 standard deviations of HCs (Figure 1 and Supplementary Table 2).
In line with previous studies, a common phenomenon in both subtypes is that the speed of processing was the most severely impaired cognitive domain (η2 = 0.655, Table 2) except for the composite score. Speed of processing appears to be a core feature of cognition, it underlies other cognitive impairments such as executive functioning and working memory (Dickinson et al., 2007). Our study further supported that impaired speed of processing plays a key role in the cognitive impairment of ECS (Sheffield et al., 2018; Lim et al., 2020). Despite ECS individuals having a lower level of education, there was no difference in the educational years between the two impaired cognitive subtypes. Moreover, severely impaired individuals showed higher negative symptoms than moderately impaired individuals but with no differences in positive symptoms, which is consistent with the clinical observation that negative symptoms and cognitive impairment might share common pathophysiological substrates (Bowins, 2011; Lincoln et al., 2017).
Working memory is described as the ability to temporarily reserve and manipulate information for further cognitive processing (Baddeley, 1992). High-level cognitive processes require the support of working memory, such as reasoning, learning, and comprehension (Baddeley, 2003). The deficit in reasoning/problem-solving domain indicates executive dysfunction (Lis et al., 2005). Impaired executive function is found to predict poor functional outcomes, failure of interventions, and restricted recovery (Green et al., 2000). By using the LPA method, we identified a subtype with intact working memory and reasoning/problem solving from 155 ECS individuals, suggesting that executive function in the subtype with mild to moderate cognitive deficits is, indeed, comparable to the healthy population. In contrast to our study hypothesis, our results found only two cognitive subtypes. A possible reason is that due to our sample size limitation, it was not possible to cluster a subtype with a small proportion of ECS individuals. Another possible reason may be that our ECS individuals had a shorter duration and more severe psychiatric symptoms compared to other studies (Carruthers et al., 2019, 2021; Lim et al., 2020).
Before analyzing brain regions with ALFF differences in the two cognitive subtypes, we first examined regions with abnormal ALFF in ECS individuals. Resting-state functional MRI analyses showed four brain regions with abnormal neural activity in ECS, including two brain regions with increased ALFF (left caudate and bilateral thalamus) and two with decreased ALFF (bilateral MPFC and bilateral PCC/PCu).
Caudate is a part of the subcortical structure, which is responsible for several neurobiological processes such as planning the behavioral execution (Grahn et al., 2008). It has been found that there was hyperactivity in caudate during working memory tasks in individuals at clinical high risk for psychosis (Thermenos et al., 2016). Individuals with schizophrenia also exhibited higher ALFF in the left caudate than HCs (Zhang et al., 2021). Thalamus is also a subcortical region, which is involved in transmitting sensory information to the cerebral cortex and regulating emotion and cognitive attention control (Sherman, 2016; Wolff and Vann, 2019). Individuals with schizophrenia exhibited reduced thalamic gray matter volume (Alemán-Gómez et al., 2020), and abnormal activation during task-related functional MRI (Byne et al., 2009). As two parts of the cortico-striatal-thalamic-cortical (CSTC) sub-circuit in the salience network (SN), abnormal functional connectivity of caudate and thalamus has been reported in previous studies (Peters et al., 2016; Huang et al., 2020). Our findings showed that increased spontaneous neural activity appeared in the CSTC sub-circuit of the SN in the early course of schizophrenia, which may suggest a compensatory mechanism to maintain normal functioning performance (Gong et al., 2020). In addition, the long-term examination should be performed to investigate whether the changes in these two brain regions still exist.
Another important finding in our study is that ECS individuals showed decreased ALFF in bilateral MPFC and bilateral PCC/PCu compared with HCs. Both MPFC and PCC/PCu are core regions of the default mode network (DMN), a crucial brain network that associates with many neurophysiological functions (Raichle, 2015). The function of MPFC and PCC/PCu was involved in introspective processes that were attenuated when attention was turned to external events (Gusnard et al., 2001). There is increasing evidence that ECS individuals exhibited decreased ALFF in MPFC and PCC/PCu compared with HCs (Ren et al., 2013; Gong et al., 2020), which is in line with our results. Taken together, our findings further support that SN and DMN play a critical role in the pathogenesis of ECS.
Although we found ALFF abnormalities in several brain regions in ECS, not all the abnormalities could reflect the neuropathological changes in cognitive impairment. We found two regions with different ALFF values between SI and MI. Compared with MI, SI presented higher ALFF values in the left precentral gyrus (PreCG) and lower ALFF values in the left PCC/PCu. After controlling the PANSS total score and medication effects, these results remained significant. Moreover, ALFF of left PCC/PCu was positively correlated with MCCB cognitive scores in ECS, while ALFF of left PreCG showed negative correlations.
Several studies focused on PCC/PCu and cognitive function have been reported. A cortical morphometric study revealed that the structural volume of PCC/PCu was associated with cognitive impairment in first-episode schizophrenia (Wang et al., 2021). Furthermore, PCC/PCu showed significant activation during the episodic memory search task (Sestieri et al., 2011). In schizophrenia individuals, the activation in PCC/PCu was not significantly enhanced during the virtual maze task, while HCs exhibited significantly enhanced activation in the same region (Siemerkus et al., 2012). PreCG, also known as the primary motor cortex, is participated in motor information processing and emotional perception (Mesulam, 1998; Watanuki et al., 2016). Abnormal regional homogeneity and voxel-mirrored homotopic connectivity of PreCG have been found in first-episode schizophrenia (Liu et al., 2018a,b). Recent evidence confirmed that the excessive activity of PreCG could result in the impairment of emotional processing in schizophrenia (Watanuki et al., 2016). Additionally, another study indicated that PreCG might be involved in some cognitive processes, such as word recognition and phonological processing (Xu et al., 2019). Based on the available evidence, our results suggest that higher spontaneous neural activity in PCC/PCu is beneficial for the preservation of cognitive function in schizophrenia individuals, whereas higher activity in PreCG plays an opposite role. In conclusion, abnormal ALFF in caudate, thalamus, MPFC, and PCC/PCu reflected the pathophysiology of ECS, with abnormal ALFF in the PCC/PCu also indicating the pathophysiology of cognitive impairment in ECS. Furthermore, the abnormal spontaneous neural activity in the PCC/PCu in schizophrenia and between cognitive subtypes provides compelling evidence for the vital effect of DMN on schizophrenia and cognitive function.
Several limitations should be considered in our study. MRI data of subjects in our study were restricted by the fact that images were collected from different acquisition centers and scanners (Jovicich et al., 2006). We used the ComBat Harmonization method to minimize the site-specific confounds and enhanced statistical power (Fortin et al., 2017). However, the ComBAt Harmonization method cannot exclude the non-linear site and scanner effects, which is one of its drawbacks. For now, this method is still the best way to minimize linear site and scanner effects (Sun et al., 2022). Second, our sample size is relatively small, we propose to expand the sample size and make validation in the independent sample population in future. Third, our study is a cross-sectional study, and further follow-up study should be conducted to validate subtype stability and outcomes. Moreover, our study did not collect full-scale IQ from schizophrenia individuals and HCs, we will collect full-scale IQ in future studies, use the subtype classification based on MCCB to compare to the subtype classification based on full-scale IQ, and investigate the differences in full-scale IQ across cognitive subtypes.
In general, our results showed that cognitive impairment in ECS might be described as two subtypes: a severely impaired group with compromised cognition across all cognitive domains and a moderately impaired group with preserved cognition in working memory and reasoning/problem solving. Furthermore, our study identified four neurobiological features of ECS and two neurobiological features of cognitive subtypes in ECS. These brain regions associated with schizophrenia and cognitive function could be potential targets for the treatment of schizophrenia and its cognitive impairment. Meanwhile, differentiating individuals into subtypes based on cognitive function could help clinicians better understand the prognosis and recovery of social function, as well as carry out individualized interventions by combining the neuroimaging features derived from the subtypes. Our findings contribute to understanding the pathophysiology of cognitive impairment in ECS from the perspective of brain spontaneous neural activity.
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 Ethics Committee of the Second Xiangya Hospital of Central South University, the Ethics Committee of the Affiliated Nanjing Brain Hospital of Nanjing Medical University, the Ethics Committee of the First Affiliated Hospital of Zhengzhou University. Written informed consent to participate in this study was provided by the participants' legal guardian/next of kin.
Author contributions
RW and JZ designed this study and wrote the protocol. RW, XS, XX, GH, YY, YL, XW, JX, YH, SG, JH, and YW performed the research. TS and WW analyzed the data and wrote the manuscript. RW, XS, and XX were involved in the revision and completion of the work. All authors contributed to the final manuscript.
Funding
This work was supported by the National Natural Science Foundation of China (Grant Nos. 81622018 and 82072096).
Acknowledgments
We would like to thank all the participants that participated in this study.
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/fnmol.2022.983995/full#supplementary-material
References
Alemán-Gómez, Y., Najdenovska, E., Roine, T., Fartaria, M. J., Canales-Rodríguez, E. J., Rovó, Z., et al. (2020). Partial-volume modeling reveals reduced gray matter in specific thalamic nuclei early in the time course of psychosis and chronic schizophrenia. Hum. Brain Mapp. 41, 4041–4061. doi: 10.1002/hbm.25108
Ammari, N., Heinrichs, R. W., and Miles, A. A. (2010). An investigation of 3 neurocognitive subtypes in schizophrenia. Schizophr Res 121, 32–38. doi: 10.1016/j.schres.2010.04.014
Baddeley, A. (2003). Working memory: looking back and looking forward. Nat. Rev. Neurosci. 4, 829–839. doi: 10.1038/nrn1201
Bora, E., and Murray, R. M. (2014). Meta-analysis of cognitive deficits in ultra-high risk to psychosis and first-episode psychosis: do the cognitive deficits progress over, or after, the onset of psychosis? Schizophr. Bull. 40, 744–755. doi: 10.1093/schbul/sbt085
Bowins, B. (2011). A cognitive regulatory control model of schizophrenia. Brain Res. Bull. 85, 36–41. doi: 10.1016/j.brainresbull.2011.02.004
Brusco, M. J., Shireman, E., and Steinley, D. (2017). A comparison of latent class, K-means, and K-median methods for clustering dichotomous data. Psychol. Methods 22, 563–580. doi: 10.1037/met0000095
Byne, W., Hazlett, E. A., Buchsbaum, M. S., and Kemether, E. (2009). The thalamus and schizophrenia: current status of research. Acta Neuropathol. 117, 347–368. doi: 10.1007/s00401-008-0404-0
Carruthers, S. P., Van Rheenen, T. E., Gurvich, C., Sumner, P. J., and Rossell, S. L. (2019). Characterising the structure of cognitive heterogeneity in schizophrenia spectrum disorders. A systematic review and narrative synthesis. Neurosci. Biobehav. Rev. 107, 252–278. doi: 10.1016/j.neubiorev.2019.09.006
Carruthers, S. P., Van Rheenen, T. E., Karantonis, J. A., and Rossell, S. L. (2021). Characterising demographic, clinical and functional features of cognitive subgroups in schizophrenia spectrum disorders: a systematic review. Neuropsychol. Rev. 25, 1–21. doi: 10.1007/s11065-021-09525-0
Celeux, G., and Soromenho, G. (1996). An entropy criterion for assessing the number of clusters in a mixture model. J. Classif. 13, 195–212.
Cobia, D. J., Csernansky, J. G., and Wang, L. (2011). Cortical thickness in neuropsychologically near-normal schizophrenia. Schizophr. Res. 133, 68–76. doi: 10.1016/j.schres.2011.08.017
Cohen, J. (1988). Statistical Power Analysis for the Behavioural Sciences. Hillside, NJ: Lawrence Erlbaum Associates.
De Meo, E., Portaccio, E., Giorgio, A., Ruano, L., Goretti, B., Niccolai, C., et al. (2021). Identifying the distinct cognitive phenotypes in multiple sclerosis. JAMA Neurol. 78, 414–425. doi: 10.1001/jamaneurol.2020.4920
Dickinson, D., Ramsey, M. E., and Gold, J. M. (2007). Overlooking the obvious: a meta-analytic comparison of digit symbol coding tasks and other cognitive measures in schizophrenia. Arch. Gen. Psychiatry 64, 532–542. doi: 10.1001/archpsyc.64.5.532
First, M. B., Spitzer, R., Gibbon, M., and Williams, J. B. (1997). Structured Clinical Interview for DSM-IV Axis I Disorders, Clinician Version (SCID-CV). Washington, DC: American Psychiatric Press.
Fortin, J. P., Cullen, N., Sheline, Y. I., Taylor, W. D., Aselcioglu, I., Cook, P. A., et al. (2018). Harmonization of cortical thickness measurements across scanners and sites. Neuroimage 167, 104–120. doi: 10.1016/j.neuroimage.2017.11.024
Fortin, J. P., Parker, D., Tunç, B., Watanabe, T., Elliott, M. A., Ruparel, K., et al. (2017). Harmonization of multi-site diffusion tensor imaging data. Neuroimage 161, 149–170. doi: 10.1016/j.neuroimage.2017.08.047
Friston, K. J., Williams, S., Howard, R., Frackowiak, R. S., and Turner, R. (1996). Movement-related effects in fMRI time-series. Magn. Reson. Med. 35, 346–355. doi: 10.1002/mrm.1910350312
Gong, J., Wang, J., Luo, X., Chen, G., Huang, H., Huang, R., et al. (2020). Abnormalities of intrinsic regional brain activity in first-episode and chronic schizophrenia: a meta-analysis of resting-state functional MRI. J. Psychiatry Neurosci. 45, 55–68. doi: 10.1503/jpn.180245
Grahn, J. A., Parkinson, J. A., and Owen, A. M. (2008). The cognitive functions of the caudate nucleus. Prog. Neurobiol. 86, 141–155. doi: 10.1016/j.pneurobio.2008.09.004
Green, M. F. (2016). Impact of cognitive and social cognitive impairment on functional outcomes in patients with schizophrenia. J. Clin. Psychiatry 77(Suppl 2), 8–11. doi: 10.4088/JCP.14074su1c.02
Green, M. F., Kern, R. S., Braff, D. L., and Mintz, J. (2000). Neurocognitive deficits and functional outcome in schizophrenia: are we measuring the “right stuff”? Schizophr. Bull. 26, 119–136. doi: 10.1093/oxfordjournals.schbul.a033430
Green, M. J., Cairns, M. J., Wu, J., Dragovic, M., Jablensky, A., Tooney, P. A., et al. (2013). Genome-wide supported variant MIR137 and severe negative symptoms predict membership of an impaired cognitive subtype of schizophrenia. Mol. Psychiatry 18, 774–780. doi: 10.1038/mp.2012.84
Gusnard, D. A., Raichle, M. E., and Raichle, M. E. (2001). Searching for a baseline: functional imaging and the resting human brain. Nat. Rev. Neurosci. 2, 685–694. doi: 10.1038/35094500
Huang, H., Botao, Z., Jiang, Y., Tang, Y., Zhang, T., Tang, X., et al. (2020). Aberrant resting-state functional connectivity of salience network in first-episode schizophrenia. Brain Imag. Behav. 14, 1350–1360. doi: 10.1007/s11682-019-00040-8
Jovicich, J., Czanner, S., Greve, D., Haley, E., van der Kouwe, A., Gollub, R., et al. (2006). Reliability in multi-site structural MRI studies: effects of gradient non-linearity correction on phantom and human data. Neuroimage 30, 436–443. doi: 10.1016/j.neuroimage.2005.09.046
Keefe, R. S., Bilder, R. M., Harvey, P. D., Davis, S. M., Palmer, B. W., Gold, J. M., et al. (2006). Baseline neurocognitive deficits in the CATIE schizophrenia trial. Neuropsychopharmacology 31, 2033–2046. doi: 10.1038/sj.npp.1301072
Kern, R. S., Nuechterlein, K. H., Green, M. F., Baade, L. E., Fenton, W. S., Gold, J. M., et al. (2008). The MATRICS consensus cognitive battery, part 2: co-norming and standardization. Am. J. Psychiatry 165, 214–220. doi: 10.1176/appi.ajp.2007.07010043
Leucht, S., Samara, M., Heres, S., and Davis, J. M. (2016). Dose equivalents for antipsychotic drugs: the DDD method. Schizophr. Bull. 42(Suppl 1), S90–94. doi: 10.1093/schbul/sbv167
Lewandowski, K. E., McCarthy, J. M., Öngür, D., Norris, L. A., Liu, G. Z., Juelich, R.J., et al. (2019). Functional connectivity in distinct cognitive subtypes in psychosis. Schizophr. Res. 204, 120–126. doi: 10.1016/j.schres.2018.08.013
Lim, K., Smucny, J., Barch, D. M., Lam, M., Keefe, R. S. E., and Lee, J. (2020). Cognitive subtyping in schizophrenia: a latent profile analysis. Schizophr. Bull. 47, 712–721. doi: 10.1093/schbul/sbaa157
Lincoln, T. M., Dollfus, S., and Lyne, J. (2017). Current developments and challenges in the assessment of negative symptoms. Schizophr. Res. 186, 8–18. doi: 10.1016/j.schres.2016.02.035
Lis, S., Krieger, S., Wilhelm, J., and Gallhofer, B. (2005). Feedback about previous action improves executive functioning in schizophrenia: an analysis of maze solving behaviour. Schizophr. Res. 78, 243–250. doi: 10.1016/j.schres.2005.02.008
Liu, Y., Guo, W., Zhang, Y., Lv, L., Hu, F., Wu, R., et al. (2018a). Decreased resting-state interhemispheric functional connectivity correlated with neurocognitive deficits in drug-naive first-episode adolescent-onset schizophrenia. Int. J. Neuropsychopharmacol. 21, 33–41. doi: 10.1093/ijnp/pyx095
Liu, Y., Zhang, Y., Lv, L., Wu, R., Zhao, J., and Guo, W. (2018b). Abnormal neural activity as a potential biomarker for drug-naive first-episode adolescent-onset schizophrenia with coherence regional homogeneity and support vector machine analyses. Schizophr. Res. 192, 408–415. doi: 10.1016/j.schres.2017.04.028
Lo, Y. T., Mendell, N. R., and Rubin, D. B. (2001). Testing the number of components in a normal mixture. Biometrika 88, 767–778. doi: 10.1093/biomet/88.3.767
McLachlan, G. J. (1978). On bootstrapping the likelihood ratio test statistic for the number of components in a normal mixture. Appl. Stat. 23, 318–324.
Mesulam, M. M. (1998). From sensation to cognition. Brain 121(Pt 6), 1013–1052. doi: 10.1093/brain/121.6.1013
Miettunen, J., Nordström, T., Kaakinen, M., and Ahmed, A. O. (2016). Latent variable mixture modeling in psychiatric research: a review and application. Psychol. Med. 46, 457–467. doi: 10.1017/S0033291715002305
Morar, B., Dragović, M., Waters, F. A., Chandler, D., Kalaydjieva, L., and Jablensky, A. (2011). Neuregulin 3 (NRG3) as a susceptibility gene in a schizophrenia subtype with florid delusions and relatively spared cognition. Mol. Psychiatry 16, 860–866. doi: 10.1038/mp.2010.70
Muthén, L. K., and Muthén, B. O. (2015). Mplus User's Guide, 7th Edn. Los Angeles, CA: Muthén and Muthén.
Nuechterlein, K. H., Barch, D. M., Gold, J. M., Goldberg, T. E., Green, M. F., and Heaton, R. K. (2004). Identification of separable cognitive factors in schizophrenia. Schizophr. Res. 72, 29–39. doi: 10.1016/j.schres.2004.09.007
Nuechterlein, K. H., Green, M. F., Kern, R. S., Baade, L. E., Barch, D. M., Cohen, J. D., et al. (2008). The MATRICS consensus cognitive battery, part 1: test selection, reliability, and validity. Am. J. Psychiatry 165, 203–213. doi: 10.1176/appi.ajp.2007.07010042
Ortiz-Gil, J., Pomarol-Clotet, E., Salvador, R., Canales-Rodríguez, E. J., Sarró, S., Gomar, J. J., et al. (2011). Neural correlates of cognitive impairment in schizophrenia. Br. J. Psychiatry 199, 202–210. doi: 10.1192/bjp.bp.110.083600
Peters, S. K., Dunlop, K., and Downar, J. (2016). Cortico-striatal-thalamic loop circuits of the salience network: a central pathway in psychiatric disease and treatment. Front. Syst. Neurosci. 10, 104. doi: 10.3389/fnsys.2016.00104
Radua, J., Vieta, E., Shinohara, R., Kochunov, P., Quidé, Y., Green, M. J., et al. (2020). Increased power by harmonizing structural MRI site differences with the ComBat batch adjustment method in ENIGMA. Neuroimage 218, 116956. doi: 10.1016/j.neuroimage.2020.116956
Raichle, M. E. (2015). The brain's default mode network. Annu. Rev. Neurosci. 38, 433–447. doi: 10.1146/annurev-neuro-071013-014030
Ren, W., Lui, S., Deng, W., Li, F., Li, M., Huang, X., et al. (2013). Anatomical and functional brain abnormalities in drug-naive first-episode schizophrenia. Am. J. Psychiatry 170, 1308–1316. doi: 10.1176/appi.ajp.2013.12091148
Rodriguez, M., Zaytseva, Y., Cvrčková, A., Dvoraček, B., Dorazilová, A., Jonáš, J., et al. (2019). Cognitive profiles and functional connectivity in first-episode schizophrenia spectrum disorders: linking behavioral and neuronal data. Front. Psychol. 10, 689. doi: 10.3389/fpsyg.2019.00689
Schreiber, J. B. (2017). Latent class analysis: an example for reporting results. Res. Soc. Adm. Pharm. 13, 1196–1201. doi: 10.1016/j.sapharm.2016.11.011
Schwarz, G. (1978). Estimating the dimension of a model. Ann. Stat. 6, 461–464. doi: 10.1214/aos/1176344136
Sclove, S. L. (1987). Application of model-selection criteria to some problems in multivariate analysis. Psychometrika 52, 333–343.
Sestieri, C., Corbetta, M., Romani, G. L., and Shulman, G. L. (2011). Episodic memory retrieval, parietal cortex, and the default mode network: functional and topographic analyses. J. Neurosci. 31, 4407–4420. doi: 10.1523/JNEUROSCI.3335-10.2011
Sheffield, J. M., Karcher, N. R., and Barch, D. M. (2018). Cognitive deficits in psychotic disorders: a lifespan perspective. Neuropsychol. Rev. 28, 509–533. doi: 10.1007/s11065-018-9388-2
Sherman, S. M. (2016). Thalamus plays a central role in ongoing cortical functioning. Nat. Neurosci. 19, 533–541. doi: 10.1038/nn.4269
Shi, C., Kang, L., Yao, S., Ma, Y., Li, T., Liang, Y., et al. (2015). The MATRICS consensus cognitive battery (MCCB): co-norming and standardization in China. Schizophr. Res. 169, 109–115. doi: 10.1016/j.schres.2015.09.003
Siemerkus, J., Irle, E., Schmidt-Samoa, C., Dechent, P., and Weniger, G. (2012). Egocentric spatial learning in schizophrenia investigated with functional magnetic resonance imaging. Neuroimage Clin. 1, 153–163. doi: 10.1016/j.nicl.2012.10.004
Smucny, J., Iosif, A. M., Eaton, N. R., Lesh, T. A., Ragland, J. D., Barch, D. M., et al. (2020). Latent profiles of cognitive control, episodic memory, and visual perception across psychiatric disorders reveal a dimensional structure. Schizophr. Bull. 46, 154–162. doi: 10.1093/schbul/sbz025
Sun, D., Rakesh, G., Haswell, C. C., Logue, M., Baird, C. L., O'Leary, E. N., et al. (2022). A comparison of methods to harmonize cortical thickness measurements across scanners and sites. Neuroimage 261, 119509. doi: 10.1016/j.neuroimage.2022.119509
Thermenos, H. W., Juelich, R. J., DiChiara, S. R., Mesholam-Gately, R. I., Woodberry, K. A., Wojcik, J., et al. (2016). Hyperactivity of caudate, parahippocampal, and prefrontal regions during working memory in never-medicated persons at clinical high-risk for psychosis. Schizophr. Res. 173, 1–12. doi: 10.1016/j.schres.2016.02.023
Van Rheenen, T. E., Cropley, V., Zalesky, A., Bousman, C., Wells, R., Bruggemann, J., et al. (2018). Widespread volumetric reductions in schizophrenia and schizoaffective patients displaying compromised cognitive abilities. Schizophr. Bull. 44, 560–574. doi: 10.1093/schbul/sbx109
Wang, C., Oughourlian, T., Tishler, T. A., Anwar, F., Raymond, C., Pham, A. D., et al. (2021). Cortical morphometric correlational networks associated with cognitive deficits in first episode schizophrenia. Schizophr. Res. 231, 179–188. doi: 10.1016/j.schres.2021.04.001
Wang, P., Yang, J., Yin, Z., Duan, J., Zhang, R., Sun, J., et al. (2019). Amplitude of low-frequency fluctuation (ALFF) may be associated with cognitive impairment in schizophrenia: a correlation study. BMC Psychiatry 19, 30. doi: 10.1186/s12888-018-1992-4
Watanuki, T., Matsuo, K., Egashira, K., Nakashima, M., Harada, K., Nakano, M., et al. (2016). Precentral and inferior prefrontal hypoactivation during facial emotion recognition in patients with schizophrenia: a functional near-infrared spectroscopy study. Schizophr. Res. 170, 109–114. doi: 10.1016/j.schres.2015.11.012
Wexler, B. E., Zhu, H., Bell, M. D., Nicholls, S. S., Fulbright, R. K., Gore, J. C., et al. (2009). Neuropsychological near normality and brain structure abnormality in schizophrenia. Am. J. Psychiatry 166, 189–195. doi: 10.1176/appi.ajp.2008.08020258
Wolff, M., and Vann, S. D. (2019). The cognitive thalamus as a gateway to mental representations. J. Neurosci. 39, 3–14. doi: 10.1523/JNEUROSCI.0479-18.2018
Woolf, B. (1957). The log likelihood ratio test (the G-test); methods and tables for tests of heterogeneity in contingency tables. Ann. Hum. Genet. 21, 397–409. doi: 10.1111/j.1469-1809.1972.tb00293.x
Xu, X. M., Jiao, Y., Tang, T. Y., Lu, C. Q., Zhang, J., Salvi, R., et al. (2019). Altered spatial and temporal brain connectivity in the salience network of sensorineural hearing loss and tinnitus. Front. Neurosci. 13, 246. doi: 10.3389/fnins.2019.00246
Yan, C. G., Wang, X. D., Zuo, X. N., and Zang, Y. F. (2016). DPABI: data processing and analysis for (resting-state) brain imaging. Neuroinformatics 14, 339–351. doi: 10.1007/s12021-016-9299-4
Zang, Y. F., He, Y., Zhu, C. Z., Cao, Q. J., Sui, M. Q., Liang, M., et al. (2007). Altered baseline brain activity in children with ADHD revealed by resting-state functional MRI. Brain Dev. 29, 83–91. doi: 10.1016/j.braindev.2006.07.002
Zhang, Y., Liao, J., Li, Q., Zhang, X., Liu, L., Yan, J., et al. (2021). Altered resting-state brain activity in schizophrenia and obsessive-compulsive disorder compared with non-psychiatric controls: commonalities and distinctions across disorders. Front. Psychiatry 12, 681701. doi: 10.3389/fpsyt.2021.681701
Keywords: schizophrenia, cognitive subtypes, latent profile analysis, MCCB, MRI, ALFF
Citation: Shao T, Wang W, Hei G, Yang Y, Long Y, Wang X, Xiao J, Huang Y, Song X, Xu X, Gao S, Huang J, Wang Y, Zhao J and Wu R (2022) Identifying and revealing different brain neural activities of cognitive subtypes in early course schizophrenia. Front. Mol. Neurosci. 15:983995. doi: 10.3389/fnmol.2022.983995
Received: 01 July 2022; Accepted: 07 September 2022;
Published: 03 October 2022.
Edited by:
Arianna Bellucci, University of Brescia, ItalyReviewed by:
Johanna Seitz-Holland, Brigham and Women's Hospital and Harvard Medical School, United StatesTianhong Zhang, Shanghai Jiao Tong University, China
Li Hui, Suzhou Guangji Hospital, China
Bing Liu, Beijing Normal University, China
Copyright © 2022 Shao, Wang, Hei, Yang, Long, Wang, Xiao, Huang, Song, Xu, Gao, Huang, Wang, Zhao 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: Renrong Wu, d3VyZW5yb25nJiN4MDAwNDA7Y3N1LmVkdS5jbg==
†These authors have contributed equally to this work