N-Acetyl-Cysteine Supplementation Improves Functional Connectivity Within the Cingulate Cortex in Early Psychosis: A Pilot Study

Abstract Background There is increasing evidence that redox dysregulation, which can lead to oxidative stress and eventually to impairment of oligodendrocytes and parvalbumin interneurons, may underlie brain connectivity alterations in schizophrenia. Accordingly, we previously reported that levels of brain antioxidant glutathione in the medial prefrontal cortex were positively correlated with increased functional connectivity along the cingulum bundle in healthy controls but not in early psychosis patients. In a recent randomized controlled trial, we observed that 6-month supplementation with a glutathione precursor, N-acetyl-cysteine, increased brain glutathione levels and improved symptomatic expression and processing speed. Methods We investigated the effect of N-acetyl-cysteine supplementation on the functional connectivity between regions of the cingulate cortex, which have been linked to positive symptoms and processing speed decline. In this pilot study, we compared structural connectivity and resting-state functional connectivity between early psychosis patients treated with 6-month N-acetyl-cysteine (n = 9) or placebo (n = 11) supplementation with sex- and age-matched healthy control subjects (n = 74). Results We observed that 6-month N-acetyl-cysteine supplementation increases functional connectivity along the cingulum and more precisely between the caudal anterior part and the isthmus of the cingulate cortex. These functional changes can be partially explained by an increase of centrality of these regions in the functional brain network. Conclusions N-acetyl-cysteine supplementation has a positive effect on functional connectivity within the cingulate cortex in early psychosis patients. To our knowledge, this is the first study suggesting that increased brain glutathione levels via N-acetyl-cysteine supplementation may improve brain functional connectivity.


Introduction
Mounting evidence indicates that redox dysregulation and resulting oxidative stress are key players in the pathophysiology of schizophrenia (Steullet et al., 2016). One pathway leading to redox dysregulation in schizophrenia is a deficit in glutathione (GSH), the main antioxidant and redox regulator in the brain, which has been shown to be decreased in some patients' brain, as directly measured from cerebrospinal fluid or via MR spectroscopy in the frontal lobe (Das et al., 2019;Do et al., 2000). The redox dysregulation model proposes that the combined action of a GSH synthesis deficit of genetic origin and an excess of oxidative stress caused by environmental factors during neurodevelopment, together with neuroinflammation and glutamatergic hypofunction, impairs neural connectivity and synchronization through fast-spiking parvalbumin (GABA) interneuron impairments and deficits in myelination (Do et al., 2009;Hardingham and Do, 2016;Steullet et al., 2016).
In the light of these mechanisms, there has been a great interest in N-acetyl-cysteine (NAC), a GSH precursor and antioxidant. Indeed, the first double-blinded, randomized, placebo-controlled clinical trial by Berk et al. showed that NAC is safe and effective as an add-on to antipsychotic medication to improve negative symptoms and antipsychotic-associated side effects in patients with chronic schizophrenia . Furthermore, improvements of negative symptoms were replicated in 2 other studies (Farokhnia et al., 2013;Sepehrmanesh et al., 2018). More recently, positive effects of NAC on cognition (Rapado-Castro et al., 2017;Sepehrmanesh et al., 2018) and positive symptoms (Sepehrmanesh et al., 2018) have also been reported.
In a recent randomized controlled trial by our group on the effect of NAC in early psychosis, we reported that a 6-month NAC add-on treatment increased brain GSH levels by 23% in the medial prefrontal cortex (mPFC) of early psychosis patients (EPPs) . Despite no change in negative symptoms, possibly due to low baseline levels, we observed an improvement in positive symptoms in a subgroup of patients showing high baseline peripheral oxidative status . In the same trial, we also reported improvement in cognition (particularly, processing speed)  and low-level auditory processing (Retsa et al., 2018), suggesting a more general effect of NAC intake and GSH levels on signal processing and sensory integration. Accordingly, in patients with the NAC add-on treatment, we also observed an improvement in white matter diffusion properties in the fornix. The white matter changes correlated with the augmentation of brain GSH levels, suggesting a possible restorative process along the fornix bundle (Klauser et al., 2018).
Taken together, these findings indicate that NAC supplementation can increase brain GSH levels and improve symptoms and processing speed in EPPs, possibly due to underlying changes in white matter diffusion properties. Indeed, white matter diffusion properties and synchronizations are disrupted in schizophrenia as reported in structural (Camchong et al., 2011;Whitford et al., 2014;Griffa et al., 2015;Klauser et al., 2017) and functional (Calhoun et al., 2009;Camchong et al., 2011) connectivity studies in both chronic schizophrenia and EPPs (Alonso-Solis et al., 2012;Zhang et al., 2014;Li et al., 2017). The most consistent brain network alterations evolve in a rostro-caudal fashion, from frontal regions at early stages, to a more widespread dysconnectivity involving all cerebral lobes and including the cerebellum in long-term schizophrenia (Pettersson-Yeo et al., 2011;Fornito et al., 2012;Van Den Heuvel and Fornito, 2014;Bartholomeusz et al., 2017;Skåtun et al., 2017). Some studies showed that these alterations correlate with symptom severity (Camchong et al., 2011;Whitford et al., 2014;Li et al., 2017).
Furthermore, functional and structural connectivity values in the cingulate correlate with mPFC GSH levels in healthy controls. In EPPs, only the association between GSH levels and structural connectivity values is preserved, whereas the association between GSH and functional connectivity is disrupted . These findings suggest that brain connectivity features in the cingulate may be critically related to brain GSH levels measured in the mPFC and that this association may be altered early in the time course of psychosis.
Growing evidence shows that NAC can alleviate several symptomatic dimensions of psychosis Lavoie et al., 2008;Carmeli et al., 2012;Shungu, 2012;Farokhnia et al., 2013), and it can restore changes in white matter diffusion properties. However, the effect of NAC intake on functional connectivity alterations in early psychosis has not yet been tested. In the current study, we investigated whether 6-month NAC supplementation can restore functional connectivity along the cingulum bundle. There were 3 reasons for focusing on the cingulum. First, the cingulum bundle is a key region implicated in the pathophysiology of schizophrenia and whose structural integrity has been linked with positive symptoms (Whitford et al., 2014;Kates et al., 2015) and impairment of processing

Significance Statement
In clinical trials, N-acetyl-cysteine as an add-on treatment to patients with schizophrenia and early psychosis patients showed significant increase of glutathione brain concentration and improvements in symptomatic expression and processing speed. We here provide a novel insight into the effect of N-acetyl-cysteine supplementation by investigating brain functional connectivity of early psychosis patients. The functional connectivity between the caudal anterior cingulate cortex and the isthmus of cingulate cortex increases after the N-acetyl-cysteine supplementation. This main finding of our study showed that this add-on might be useful for restoring functional connectivity alterations in early psychosis patients, which have been linked to symptoms in previous studies.
speed (Karbasforoushan et al., 2015;Kochunov et al., 2017), 2 symptomatic dimensions of which were improved in the study by Conus and colleagues . Second, we previously observed a correlation between brain GSH levels and functional connectivity along the cingulum in healthy individuals, indicating that higher GSH levels correspond to better synchronization of the medial prefrontal and posterior cingulate cortices. This association was disturbed in EPPs . Lastly, we recently showed that NAC has the ability to increase GSH levels in the mPFC, highlighting again that the cingulate is an important region to focus on. The current study is based on a subgroup of the cohort investigated in , which was composed of patients who agreed to participate in a very demanding imaging study. Given the increase in GSH levels in the mPFC in this same cohort, we hypothesized that NAC supplementation would lead to an increase in functional and structural connectivity along the cingulum bundle, which would correlate with changes in brain GSH levels.

Clinical Trial Protocol
NAC (2700 mg/d) and placebo were administered to EPPs for 6 months following a double-blinded randomized design. Restingstate functional magnetic resonance imaging (fMRI), diffusion spectrum imaging (DSI), and magnetic resonance spectroscopy (MRS) were performed at baseline and at the end of the study after 6 months of NAC administration (Swiss Medic [2008DR2308], ClinicalTrial.gov [NCT01354132]). The study reported here is based on a subsample of a larger clinical trial .

MRI Study Participants
Participants were recruited from the Treatment and Early Intervention in Psychosis Program (Lausanne University Hospital, Switzerland) (Baumann et al., 2013). Patients meeting criteria for early psychosis were included in the trial. The Lausanne Psychosis cohort was created under the supervision of ethics boards with representatives at different levels (CHUV/Canton/Fédéral). The study was approved by the local research ethic committee (Commission cantonale d'éthique de la recherche sur l'être humain), and written informed consent was obtained from all participants. A complete description of the study and of the cohort can be found elsewhere . Among the 63 randomized participants of this study (32 for the NAC arm and 31 for the placebo arm), 20 patients (9 NAC, 11 placebo) agreed to participate in a complementary MRI study. These 20 patients (7 women, 13 men; 15 right-handed, 4 left-handed, 1 ambidextrous; aged 25 ± 6 years) were scanned at baseline and at the 6-month follow-up after NAC or placebo add-on. A total of 74 age and gender-matched healthy controls (27 women, 47 men; 62 right-handed, 9 left-handed, 2 ambidextrous; aged 26 ± 6 years) were recruited and assessed by the Diagnostic Interview for Genetic Studies (Preisig et al., 1999). Major mood, psychotic or substance-use disorders, and having a first-degree relative with a psychotic disorder were exclusion criteria for the healthy subjects. Neurological disorders and severe head trauma were exclusion criteria for all subjects.

Structural, Diffusion, and Functional MRI Acquisitions
T1-weighted volumes and fMRI recordings were acquired on a 3-Tesla scanner (Trio, Siemens Medical, Germany) equipped with a 32-channel head coil. For the T1-weighted volumes, a magnetization-prepared rapid acquisition gradient echo sequence was acquired with 1-mm in-plane resolution and 1.2mm slice thickness, covering 240 ± 257 ± 160 voxels. The TR, TE, and TI were 2300, 2.98, and 900 ms, respectively. For the fMRI recordings, each subject was scanned in resting-state conditions using a standard gradient echo planar imaging sequence sensitive to blood oxygen-level-dependent contrast. An axial plane was used with a 64-± 58-voxel matrix (resolution 3.3 ± 3.3 mm 3 ). Thirty-two slices of 3.3-mm thickness with a 0.3-mm gap were acquired. Acquisition and repetition times were 9 minutes and 1920 milliseconds, respectively. A DSI sequence was also acquired during the same MRI session with a 2.2-± 2.2-± 3-mm 3 resolution, covering 96 ± 96 ± 34 voxels. The TR and TE were 6100 and 144 ms, respectively. A q4-half acquisition scheme was used with a maximum b-value of 8000 s/mm 2 and 1 b o volume (Wedeen et al., 2008). The same MRI protocols were applied to both the EPP and control groups.

MRS Acquisitions
The MRS experiment was performed with a transverse electromagnetic volume head coil (MR Instruments, Minneapolis, MN). The magnetic field homogeneity was optimized by adjusting first-and second-order shims using FAST(EST)MAP (Gruetter, 1993). In vivo proton nuclear magnetic resonance spectra were acquired from a volume of interest positioned in the medial bilateral prefrontal lobe using the short-TE spin-echo fullintensity acquired localized single voxel spectroscopy technique (Mekle et al., 2009). The scan covered a volume of interest of 20 ± 25 ± 25 mm 3 with TR/TE = 4000/6 ms, an acquisition bandwidth of 2 kHz, an averaged number of 148, and a vector size of 2048. The GSH concentration was quantified via the watersuppressed in vivo proton nuclear magnetic resonance imaging resonance spectra using LCModel (Provencher, 1993) with the unsuppressed water proton nuclear magnetic resonance spectra as an internal reference. The spectral range for analysis was set to 0.2 to 4.2 ppm. GSH was quantified with a Cramer-Rao lower bounds of 10 ± 3%. More details can be found in Xin and al. (Xin et al., 2016).

Data Processing
Magnetization-prepared rapid acquisition gradient echo volumes were segmented into white matter, gray matter, and cerebrospinal fluid compartments. The gray matter volume was parcellated into 68 cortical and 14 subcortical anatomical regions and the brainstem according to the Desikan-Killiany atlas (Desikan et al., 2006) using FreeSurfer software (Fischl et al., 2002) (https://surfer.nmr.mgh.harvard.edu/).
Whole-brain tractography was performed on reconstructed DSI data using CMTK software (Daducci et al., 2012) (http:// connectomics.org/), which allowed individual structural connectivity matrices to be estimated. The number of streamlines connecting each pair of gray matter regions was used as a measure of structural connectivity strength. Furthermore, we built a group-representative structural connectome by averaging the individual structural connectivity matrices over the 74 healthy subjects. Structural connections that were not present in all the healthy subjects were discarded from further analyses.
fMRI data were processed according to a state-of-the-art pipeline including discarding the first 5 time points for signal stabilization, slice-timing correction, motion correction, regressions of 6 motion parameters, averaging white matter and cerebrospinal fluid signals, linear detrending, and bandpass filtering (0.01-0.1 Hz) using CMTK software. The structural volumes were linearly registered to the individual fMRI spaces. Voxelwise fMRI time series were averaged over the 83 gray matter regions. Functional connectivity was assessed by computing the Pearson's correlation coefficient between the temporal signals of brain region pairs. Negative values were not discarded, but absolute values were instead taken. Indeed, recent studies confirm that negative correlation values may have biological origins (Chai et al., 2012) and present specific features in schizophrenia (Parente et al., 2018). We filtered the functional connections using the group-representative structural connectome. We performed the filtering to discard the weakest functional connections (Goñi et al., 2014) and reduce the number of false positives without introducing a threshold bias (Drakesmith et al., 2015;van den Heuvel et al., 2017).

Connectivity Analysis
The resting-state functional network of each subject was defined using the functional connectivity matrix, with the 83 brain regions corresponding to the 83 nodes of the network and the pairwise Pearson's correlation values representing the edges' strength. The structural network was defined using the number of streamlines as the edge weight. As a first step, we investigated the effect of NAC or placebo intake on the functional and structural connectivity of the cingulate cortex, which have been shown to be impaired in schizophrenia (Van Den Heuvel and Fornito, 2014;Whitford et al., 2014;Monin et al., 2015). More specifically, we considered the functional and structural connectivity between 3 anterior (medial orbitofrontal, rostral anterior cingulate, and caudal anterior cingulate cortices) and 2 posterior (posterior cingulate and isthmus of cingulate cortices) regions of the cingulate cortex (Figure 1), resulting in 6 connections to be tested. The connectivity values of the left and right hemispheres were averaged. We compared the baseline connectivity values with the 6-month connectivity values for both the NAC and placebo EPP groups. Furthermore, we compared the brain connectivity of each EPP group with that of the control subjects. Next, we characterized integration and segregation properties of the functional and structural brain networks using graph theoretical measures. At the global level, we computed the global efficiency and average clustering coefficient of each individual (structural and functional) brain network. The global efficiency represents the efficiency of communication through shortest paths in the brain network, while the average clustering coefficient quantifies the level of local connectedness in the network. At the local level, we computed the node betweenness centrality and the edge betweenness centrality for the regions and connections of the cingulate cortex (Figure 1). The betweenness centrality quantifies the number of shortest paths in the network passing through a given node or connection. Here, the betweenness centrality quantifies the centrality of the cingulate cortex regions and cingulum connections within the overall brain network. The graph measures were computed with the MATLAB Brain Connectivity Toolbox (Rubinov and Sporns, 2010). The formal definitions of these measures can be found in (Rubinov and Sporns, 2010).

Statistical Analysis
The relationship between the brain connectivity values and the GSH levels was assessed with the Spearman correlation coefficient. Statistical analyses were performed using MATLAB functions (MATLAB R2017a) and in-house code.

Results
First, we investigated the effect of NAC on functional connectivity in the cingulate cortex. To do this, we analyzed the relative variation of functional connectivity values between the anterior and posterior cingulate regions before and after the supplementation with NAC or placebo. The median of functional connectivity differences before and after supplementation (ΔFC) was significantly higher in NAC compared with placebo patients for the caudal anterior cingulate-isthmus of cingulate connection (1-sided Wilcoxon: P = .006, P = .03 adjusted with Bonferroni correction) (Figure 2a). There were no significant changes for the other connections between the anterior and posterior parts of the cingulate cortex. For this reason, we focused our further analyses on characterizing the functional connectivity changes between the caudal anterior cingulate and isthmus of cingulate cortices. In a second step, we extended our investigation to possible differences in functional connectivity between the caudal anterior cingulate and isthmus of cingulate cortices in patients compared with the healthy controls. The distributions of functional connectivity values for both patient's baseline and follow-up data were calculated. Post-hoc tests were performed to compare the NAC and placebo groups with a control group. No significant changes were found between EPPs at baseline and healthy controls (one-sided Wilcoxon P = .27). We observed a heterogeneity in the 2 EPPs baseline groups, NAC and placebo, which was not significant (Wilcoxon P = .068). However, the results showed a significant increase in functional connectivity at the follow-up compared with baseline for the NAC patients (1-sided paired Wilcoxon, P = .0039), whereas no change was observed for the placebo patients (1-sided paired Wilcoxon P = .41), as expected from the results of Figure 2A. Moreover, the functional connectivity between the caudal and isthmus of the cingulate cortices in the NAC group after 6-month supplementation was larger than that of the control group (1-sided Wilcoxon, P = .0088) (Figure 2b).
Similarly, we explored possible structural connectivity alterations between the anterior and the posterior parts of the cingulate cortex (5 connections of interest; Figure 1); as such we compared the NAC, placebo, and healthy control groups. We used the number of streamlines as an edge weight to characterize the structural connectivity. No significant differences were found after correction for multiple comparisons. Nevertheless, in view of the findings regarding functional connectivity impairments, we focused our analysis on the caudalisthmus connection. After 6-month supplementation of NAC, this connection showed a trend of an increased number of streamlines connecting the caudal and isthmus of the cingulate cortices (1-sided paired Wilcoxon, P = .082, P = .41 adjusted with Bonferroni correction) (Figure 3).
We hypothesized that the increase under NAC supplementation of the functional connectivity between the caudal cingulate and isthmus of the cingulate cortices could imply a reorganization of the brain functional network around this specific connection. First, we characterized the global topological changes of the brain functional networks by computing the global efficiency and average clustering coefficient of the network. No significant changes were found between the baseline and follow-up data for these global measures. Second, we assessed the nodal betweenness centrality of the caudal cingulate and isthmus of cingulate regions, and the edge betweenness centrality of their functional connection. In regard to these local network measures, the change (between baseline and 6-month follow-up) in edge between centrality was significantly larger in NAC patients than in placebo patients (1-sided Wilcoxon, P = .0073) (Figure 4a). Moreover, the change in nodal betweenness centrality for the isthmus of cingulate region was significantly higher for NAC patients than for placebo patients (1-sided Wilcoxon, P = .0049) (Figure 4b). The same trend was found for the caudal anterior cingulate region (1-sided Wilcoxon, P = .11 for average left-right hemispheres, significant in the right hemisphere, P = .01).
On one hand, our results showed an increase of functional connectivity in the cingulate regions. On the other hand, Conus and al. (Conus et al., 2018) already showed in the larger cohort of this clinical trial a significant increase of GSH brain concentration, suggesting that the supplementation of NAC, precursor of GSH, was able to help restore the redox dysregulation mechanism responsible for connectivity alteration in schizophrenia (Steullet et al., 2016). To verify this hypothesis, we also investigated the relationship between ΔFC and GSH changes (ΔGSH) in NAC patients. The correlation between ΔFC and ΔGSH (Spearman's correlation coefficient, r = 0.33, P = .17) was not significant. Then, a linear discriminant analysis was performed on the (ΔFC, ΔGSH) values with respect to the 2 groups of NAC patients and placebo patients ( Figure 5). The classification error on our dataset was 11%. We estimated the expected classification error on a new dataset at 11% using a k-fold cross-validation procedure based on the resampling of our data.

Discussion
The main finding of this study is that functional connectivity along the cingulum bundle was increased after 6-month NAC supplementation in EPPs. We observed that functional connectivity between the caudal anterior cingulate and isthmus of cingulate cortices increased with NAC supplementation. We further hypothesized that baseline functional connectivity was decreased between EPPs compared with healthy controls along the cingulum bundle. However, no significant differences in functional connectivity were found between the 2 groups, which may be due to the effect of the disease in early psychosis not being strong enough or the sample being too small. The functional connectivity increase after NAC supplementation might be related to compensatory mechanisms, considering that functional connectivity loss in the cingulate regions have been previously observed in patients with psychosis (Van Den Heuvel and Fornito, 2014;Skåtun et al., 2017b). Changes in functional connectivity after NAC intake has been studied in the context of smoking cessation (Froeliger et al., 2015), where NAC restores fronto-striatal resting-state functional connectivity in preventing relapse during smoking cessation. In summary, this is the first study examining the impact of NAC on functional connectivity in EPPs. In spite of a heterogeneity between our 2 EPPs groups at baseline, the relative difference of functional connectivity between the anterior and posterior regions of the cingulate cortex was strengthened in the NAC group after 6 months of treatment compared with the placebo group and baseline values in healthy controls.
Changes in functional connectivity are difficult to interpret quantitatively because of the complex physiology of the BOLD signal (Gauthier and Fan, 2019). However, evidence relates the BOLD signal to the electrophysiological brain activity (Logothetis and Wandell, 2004) and suggests that an increase in FC could reflect higher synchronization between the corresponding brain regions. This higher synchronization could have different origins and could be related to, for instance, the reorganization of functional brain networks due to white matter alterations.
We thus provide evidence that the functional connections along the cingulum may be rearranged by NAC supplementation.
Indeed, edge betweenness centrality was significantly increased for the same caudal-isthmus of cingulate connection. A higher number of shortest paths passed through this connection after NAC supplementation, suggesting that this connection becomes more efficient in linking rostral to caudal regions of the cingulate cortex. This functional network reorganization could be a local compensatory process to alleviate the overall connectivity disruptions reported in the literature by exploiting new communication pathways and could reflect underlying white matter changes.
Functional connectivity abnormalities in schizophrenia may be underpinned by white matter alterations (Zalesky et al., 2011;Klauser et al., 2017). Thus, we next examined whether changes in structural connectivity may underlie changes in functional connectivity. Indeed, a restoration of structural alterations could explain improvements in functional connectivity by rehabilitating the structural pathways that convey the transmission of information. Even if the relationship between structural and functional connectivity is not yet fully understood, they are strongly related and functional connectivity to some extent reflects the underlying structural scaffold (Honey et al., 2010;Cabral et al., 2013;Goñi et al., 2014). In our sample, the structural connectivity (as quantified by the number of streamlines) between the caudal and isthmus of cingulate regions was increased to a trend level after 6-month NAC supplementation. A longer period of molecule supplementation or a larger cohort might be necessary to detect statistically significant changes in the white matter organization. In this regard, it is interesting to note that in the same cohort, 6-month NAC supplementation improved the white matter properties of the body of the fornix (as assessed with diffusion MRI and generalized fractional anisotropy) (Klauser et al., 2018). This white matter region is much thinner than the cingulum and has been shown to be very vulnerable to oxidative stress (Corcoba et al., 2015).
One of the possible mechanisms for explaining the observed functional connectivity alterations in schizophrenia is the hypothesis that a redox dysregulation along with neuroinflammation and glutamatergetic system NMDAR hypofunction would cause macro-and microcircuit impairments resulting in altered functional connectivity (Steullet Figure 5. Linear discriminant analysis. Δ Glutathione (GSH) and median of the functional connectivity differences before and after supplementation (ΔFC) are used as features of a linear discriminant analysis to classify the subjects into the 2 arms. The misclassification error was estimated at 5% via a resubstitution validation and 11% using a k-fold cross-validation procedure. et al., 2016). To explore this hypothesis, we investigated the relationship between the brain GSH concentration and the functional connectivity changes observed in our sample. However, the direct relationship between GSH levels and functional connectivity values between the caudal and isthmus cortices was not significant. Nevertheless, we demonstrated a dependence between these parameters by using linear discriminant analysis to predict a subject's alliance to the NAC or placebo groups. Indeed, ΔFC and ΔGSH values allowed for accurate classification of the patients into the 2 groups. Replication of these results in a larger cohort could reinforce the evidence for the redox dysregulation model. Moreover, cingulate functional connectivity values and GSH levels may be useful markers in monitoring the NAC treatment response as these 2 markers have already been related to symptoms and cognitive scores (Yan et al., 2012;Conus et al., 2018).
Our study focused only on the connectivity within the cingulate cortex. Indeed, the human cingulate cortex is a functionally heterogeneous structure (Margulies et al., 2007), which can be subdivided into regions with specific functional and neurobiological profiles (Yu et al., 2011). Moreover, diffusion properties of the cingulate subdivisions allow for differentiation of the uniform cingulum white matter tract into more specific tracts (Jones et al., 2013). In schizophrenia, functional connectivity between these cingulate subdivisions and the other brain regions is significantly altered (Wang et al., 2015). However, the nature of the alteration is highly dependent on the fMRI signal processing steps, particularly the regression of the global signal (Wang et al., 2015). Similarly, aberrant functional connectivity has been found in schizophrenia patients between the default mode network, notably involving the medial orbitofrontal cortex, the posterior cingulate cortex, and/or some regions of the cingulate cortex and other functional systems. These aberrations describe either a hyperconnectivity of the default mode network (Whitfield-Gabrieli et al., 2009) or a mix of hyper-and hypoconnectivity dependent on these regions (Yan et al., 2012). However, the connectivity within the cingulate cortex, as studied in our analysis, is sparsely studied. The functional connectivity within the default mode network tends to be lower in patients with schizophrenia compared with healthy populations (Bastos-Leite et al., 2015;Du et al., 2016). The changes of functional connectivity within the cingulate cortex observed after NAC supplementation showed local modification within the cingulate subnetwork, but the global effect on the whole brain remains to be studied in a larger cohort. Further analyses could explore how the alterations of functional connectivity within the cingulate subnetwork may impact the connectivity between other brain subnetworks and connectivity of other functional systems.
Betweenness centrality values are highly dependent on the reliability of the functional network estimation. One of the main limitations in this estimation is to find a balance between sensitivity and specificity by choosing a threshold for the connectivity values ). An excess of false positive or false negative connections could significantly affect a range of network measures, such as the global efficiency and clustering coefficient. However, the application of a threshold can result in artifactual differences between the average degree and the number of nodes of the different networks, thus changing the essential conditions for network comparison (van Wijk et al., 2010). The application of a proportional threshold could reduce this problem by maintaining the same density over the different networks. However, fixing the density of the network could increase the number of spurious edges in pathological networks suffering from connectivity alterations and could bias the comparison between the different arms (van den Heuvel et al., 2017;Hallquist and Hillary, 2019). In our analysis, we filtered our functional connectivity matrix by the structural scaffold without applying any threshold. This method has been developed based on the strong structure-function relationship and is designed to reduce the number of false positive connections in the functional network.
The limitations of the current study deserve some further considerations. First, the small sample size of our study limits the generalizability of the findings, which should be replicated in a larger sample. Second, the nature of the connectivity alterations in fMRI studies depend on the functional signal processing pipeline and notably the regression of the global signal. A larger cohort and a comparison of different signal processing steps would allow for a deeper exploration of the effect of NAC on connectivity of functional brain networks.
In conclusion, we investigated the connectivity changes within the cingulate cortex of EPPs after NAC supplementation. We found a significant change in functional connectivity between the caudal anterior cingulate cortex and the isthmus of cingulate cortex after NAC supplementation. Further studies are needed to assess to what extent NAC may be useful for restoring functional connectivity in the cingulate cortex in early psychosis.