First Episode Psychosis and Schizophrenia Are Systemic Neuro-Immune Disorders Triggered by a Biotic Stimulus in Individuals with Reduced Immune Regulation and Neuroprotection

There is evidence that schizophrenia is characterized by activation of the immune-inflammatory response (IRS) and compensatory immune-regulatory systems (CIRS) and lowered neuroprotection. Studies performed on antipsychotic-naïve first episode psychosis (AN-FEP) and schizophrenia (FES) patients are important as they may disclose the pathogenesis of FES. However, the protein–protein interaction (PPI) network of FEP/FES is not established. The aim of the current study was to delineate a) the characteristics of the PPI network of AN-FEP and its transition to FES; and b) the biological functions, pathways, and molecular patterns, which are over-represented in FEP/FES. Toward this end, we used PPI network, enrichment, and annotation analyses. FEP and FEP/FES are strongly associated with a response to a bacterium, alterations in Toll-Like Receptor-4 and nuclear factor-κB signaling, and the Janus kinases/signal transducer and activator of the transcription proteins pathway. Specific molecular complexes of the peripheral immune response are associated with microglial activation, neuroinflammation, and gliogenesis. FEP/FES is accompanied by lowered protection against inflammation, in part attributable to dysfunctional miRNA maturation, deficits in neurotrophin and Wnt/catenin signaling, and adherens junction organization. Multiple interactions between reduced brain derived neurotrophic factor, E-cadherin, and β-catenin and disrupted schizophrenia-1 (DISC1) expression increase the vulnerability to the neurotoxic effects of immune molecules, including cytokines and complement factors. In summary: FEP and FES are systemic neuro-immune disorders that are probably triggered by a bacterial stimulus which induces neuro-immune toxicity cascades that are overexpressed in people with reduced anti-inflammatory and miRNA protections, cell–cell junction organization, and neurotrophin and Wnt/catenin signaling.


Introduction
In 1995, Smith and Maes [1] launched the monocyte-T lymphocyte theory of schizophrenia, which incorporated neurodevelopmental factors and activation of immune pathways top five nodes with the highest degree, and top non-hub bottlenecks as nodes with the highest betweenness centrality. Together they shape the backbone of the network.
Network clustering was carried out to cluster highly interconnected genes to identify protein commonalities with similar attributes and functions. In the present study, we employed STRING to perform Markov Clustering (MCL), which robustly identifies annotated complexes, and Cytoscape plugins, namely ClusterMaker. Molecular Complex Detection (MCODE) was performed using Metascape to detect smaller components of densely connected nodes which represent molecular complexes [17].

Enrichment Analysis
The list of seed and first-order genes (divided into communities by cluster analysis and up-and downregulated genes) were extended with known protein interactions from STRING, IntAct (https://www.ebi.ac.uk/intact/, as accessed 19 September 2021), GOnet (https://tools.dice-database.org/Gonet/, as accessed 19 September 2021), Metascape (http://metascape.org, as accessed 19 September 2021), inBio Discover (https://inbio-discover. com/, as accessed 19 September 2021), Enrichr (https://maayanlab.cloud/Enrichr/, as accessed 19 September 2021), or the R package ClusterProfiler 4.0 and examined for their pathway, function, or disease enrichment scores. In this study, we use the false discovery rate (FDR) corrected p-values. We searched the networks against GO biological processes (assemblies of molecular functions in pathways), GO molecular functions (protein activities at the molecular level), GO cellular component (location of the proteins) (www.geneontology. org, as accessed 19 September 2021), Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways (https://genome.jp/kegg/, as accessed 19 September 2021), REACTOME (the European Bio-informatics Institute pathway database) (https://reactome.org, as accessed 19 September 2021), Translational Regulatory Relationships (TTRUST) (www.grnpedia.org, as accessed 19 September 2021), and DOID human disease phenotypes (Disease Ontology-Institute for Genome Sciences @ University of Maryland (disease-ontology.org, as accessed 19 September 2021). Moreover, Metascape was used to delineate and visualize the GO biological pathways and molecular processes, and PANTHER (PANTHER-Gene List Analysis (pantherdb.org), as accessed 19 September 2021), REACTOME, KEGG, and Wiki (WikiPathways-WikiPathways, as accessed 19 September 2021) pathways, which are over-represented in the gene sets. Metascape automatically clusters, hierarchically, the significant terms into a tree based on Kappa-statistical similarities among their gene memberships. This is useful because GO terms heavily overlap and, therefore, output GO terms show a large degree of redundancy. The top 10 performing terms obtained with Enrichr are displayed as bar graphs produced with Appyters.

Annotation Analysis and Annotation Visualization
Annotation analysis was employed to examine GO terms annotating the genes/proteins assembled into functional sets. GOnet analysis was employed to construct interactive graphs which, at the same time, contain genes and GO terms. Furthermore, we made custom GO term annotation lists consisting of the most representative GO terms (by preference, the significant children). inBio Discover was used to show the network with the top 4 selected DOID annotations. R package ClusterProfiler was employed to make dot plots of annotated GO term leaves, with size of the dots indicating gene number and color of the dots indicating the p-values.

The PPI Network Topography of FEP and FES/FEP
All genes/proteins were used to construct an undirected network representing the protein interactions of FEP/FES. The zero-order network consisted of 23 nodes, and the number of edges (n = 88) exceeded the expected number of edges (n = 24), with p-enrichment value of 1.0 × 10 −16 , and average node degree = 7.65 and average local clustering coefficient = 0.709. There were three singletons, namely PLVAP, DROSHA, and All genes/proteins were used to construct an undirected network representing the protein interactions of FEP/FES. The zero-order network consisted of 23 nodes, and the number of edges (n = 88) exceeded the expected number of edges (n = 24), with p-enrichment value of 1.0e-16, and average node degree = 7.65 and average local clustering coefficient = 0.709. There were three singletons, namely PLVAP, DROSHA, and C1QA. Figure  1 displays the first order protein network that showed one singleton (PLVAP) and 92 nodes; the number of edges (n = 1063) exceeded the expected number of edges (n = 273) with p-enrichment value of 1.0 × 10 −16 , average node degree = 23.1, average local clustering coefficient = 0.678, average number of neighbors = 23.109, network diameter = 4 and radius = 2, characteristic path length = 1.919, network density = 0.254, and heterogeneity = 0.586.  ) was centered around BDNF, CCL11, CCL3, CSF2, IFNG, IL10, IL12A, IL13, IL1R1,  IL4, IL5, IL6, MBP, PON1, and TNF, and (2) a second adhesion-associated cluster (yellow nodes) was centered around CDH1 and CTNNB1.
The seed proteins were used to construct a PPI network representing the protein interactions in FEP only. The first-order network (first shell) shows 65 nodes, the number of edges (n = 811) exceeded the expected number of edges (n = 193) with p-enrichment value of 1.0 × 10 −16 , with average node degree = 25, average local clustering coefficient = 0.75, average number of neighbors = 24.954, network diameter = 4 and radius = 2, characteristic path length = 1.780, network density = 0.390, and heterogeneity = 0.549. The top five hubs are, in descending order of importance, IL6 (87), TNF (85), IL10 (71), IL4 (71), and IFNG (68). BDNF was the most important bottleneck (0.0563), followed by IL6 (0.0473) and TNF (0.046).  Figure 2 shows that the most important GO terms over-represented in the upregulated genes were: inflammatory response and cytokine production, and response to LPS and receptor signaling pathway via JAK-STAT. Figure 3 shows that the most important GO terms over-represented in the downregulated genes were: neurotrophin signaling pathway, production of miRNA involved in gene silencing, and neuron projection morphogenesis.   Table 1 shows the biological interpretation of the MCODE analysis performed using multiple databases (GO biological and molecular, KEGG, WikiPaths, PANTHER, and RE-ACTOME gene sets) in FEP. MCODE performed on the upregulated genes in FEP detected one cluster which represented a response to LPS, response to molecule of bacterial origin, and inflammatory response. A second MCODE analysis conducted on the upregulated genes revealed two molecular complexes, namely signaling by interleukins and a second, which is shown in Table 1, representing TNFR1-induced NFκB signaling pathway, death receptor signaling, and TNF signaling. MCODE performed on the downregulated genes in FEP detected two small complexes, which represented: (1) neurotrophin/tropomyosin related kinase B (TrkB) signaling pathway and cellular component morphogenesis; and (2) RNA interference, production of miRNAs involved in gene silencing by miRNA, and production of small RNA involved in gene silencing by RNA. The search for a regulatory    Table 1 shows the biological interpretation of the MCODE analysis performed using multiple databases (GO biological and molecular, KEGG, WikiPaths, PANTHER, and RE-ACTOME gene sets) in FEP. MCODE performed on the upregulated genes in FEP detected one cluster which represented a response to LPS, response to molecule of bacterial origin, and inflammatory response. A second MCODE analysis conducted on the upregulated genes revealed two molecular complexes, namely signaling by interleukins and a second, which is shown in Table 1, representing TNFR1-induced NFκB signaling pathway, death receptor signaling, and TNF signaling. MCODE performed on the downregulated genes in FEP detected two small complexes, which represented: (1) neurotrophin/tropomyosin related kinase B (TrkB) signaling pathway and cellular component morphogenesis; and (2) RNA interference, production of miRNAs involved in gene silencing by miRNA, and production of small RNA involved in gene silencing by RNA. The search for a regulatory  Table 1 shows the biological interpretation of the MCODE analysis performed using multiple databases (GO biological and molecular, KEGG, WikiPaths, PANTHER, and RE-ACTOME gene sets) in FEP. MCODE performed on the upregulated genes in FEP detected one cluster which represented a response to LPS, response to molecule of bacterial origin, and inflammatory response. A second MCODE analysis conducted on the upregulated genes revealed two molecular complexes, namely signaling by interleukins and a second, which is shown in Table 1, representing TNFR1-induced NFκB signaling pathway, death receptor signaling, and TNF signaling. MCODE performed on the downregulated genes in FEP detected two small complexes, which represented: (1) neurotrophin/tropomyosin related kinase B (TrkB) signaling pathway and cellular component morphogenesis; and (2) RNA interference, production of miRNAs involved in gene silencing by miRNA, and production of small RNA involved in gene silencing by RNA. The search for a regulatory relationship using TTRUST (in Metascape) shows that the top two over-represented tran-   Figure 4 shows a heat map with the top 20 biological GO enriched terms in the first FEP/FES cluster, indicating that the most important over-represented terms were: inflammatory response, positive regulation of cytokine production, I-kappa kinase/NFκB signaling, and a response to LPS. Figure 5 shows a heat map with the top 20 biological GO enriched terms in the second cluster of FEP/FES, namely Wnt signaling pathway, cell-cell junction organization, and beta-catenin-TCF complex assembly. Table 2 shows the top GP enriched terms in the complement factors, DISC1 and DROSHA.         Electronic Supplementary File (ESF) Figure S1 shows a bar graph with the top 10 performing cellular GO terms which were accumulated in FEP/FES genes, indicating that the catenin complex, and cytoplasmatic and vesicle membranes were the most significant cellular terms. ESF Figure S2 shows the top 10 performing REACTOME terms accumulated in the FEP/FES gene list, indicating that the first terms pointed towards an immune response and that the top 7-10 terms pointed towards different aspects of Toll-Like Receptor (TLR) 3/4 signaling. ESF Figure S3 displays the top 10 PANTHER terms, which were over-expressed in the FEP/FES gene list with CCKR, TLR, apoptosis, Wnt, and cadherin signaling as top pathways. ESF Figure S4 displays the top 10 WikiPathway terms which were over-expressed in the FEP/FES gene list, including TLR4 signaling and miRNA involvement.  Table 1 shows the results of a MCODE analysis performed on all FEP/FES genes and detected a molecular complex (with core genes IFNG, IL6, CCL3, IL4, IL12A, and IL13), which represented microglial cell activation, positive regulation of tyrosine phosphorylation of STAT protein, and regulation of tyrosine phosphorylation of STAT protein. Annotation analysis also shows that a neuroinflammatory response (GO:0150076, pFDR = 1.0 × 10 −7.13 ) was accumulated in the seed gene list selection. It is interesting to note that behavior (GO:0007610; log10 p = −4.941; 10 overlapping genes) and learning and memory (GO:0007611; log10 p = −4.908; 7 overlapping genes) accumulated in the FEP/FES genes.

Enrichment Analysis in FEP/FES
ESF Figure S5 shows a network built with OmicsNet (using InAct; 2581 nodes, 4154 edges, 86 seeds). The top three over-represented transcriptional TTRUST networks were: SP1 (132 hits), NFKB1 (95), and RELA (95).        Figure 8 shows the results of annotation analysis using R package ClusterProfiler and a custom-made GO list comprising parent and child terms exploring response, cellular, or defense response to a variety of stressors, with the aim to differentiate the type of responses, including biotic responses. We observed that the most important over-represented GO terms in the FEP/FES gene list were responses to other organisms, external biotic stimulus, response to lipid, response to bacterium, and response to LPS, whereas there was much less or no evidence for a response to a virus or other biotic (fungal, parasites) or abiotic stimuli.  Figure 8 shows the results of annotation analysis using R package ClusterProfiler and a custom-made GO list comprising parent and child terms exploring response, cellular, or defense response to a variety of stressors, with the aim to differentiate the type of responses, including biotic responses. We observed that the most important over-represented GO terms in the FEP/FES gene list were responses to other organisms, external biotic stimulus, response to lipid, response to bacterium, and response to LPS, whereas there was much less or no evidence for a response to a virus or other biotic (fungal, parasites) or abiotic stimuli.

Annotation Analysis and Visualization
Consequently, we explored the network of interrelated genes using inBio Discover to delineate which diseases (DOID) are over-represented in the query gene list. Table 3 shows the top nine DOID disease annotations. This table indicates that the proteins involved in FES/FEP are enriched in a number of immune diseases, including immunodeficiency disease and autoimmune disease, and three gastro-intestinal diseases, namely intestinal disease, inflammatory bowel disease, and colitis. ESF Figure S6 shows the extended network and those four of the nine top annotations that were over-represented. The nodes of the extended PPI network involved in immune system and autoimmune disease and inflammatory bowel and intestinal disease are colored in red, orange, blue, and green, respectively. ESF Figure S7 shows all nodes (in red color) which were accumulated in intestinal disease (DOID:5295). ESF Figure S8 shows the results of an annotation analysis for all differentially expressed proteins, using R package ClusterProfiler and the GO list consisting of only leaf children in biological processes. Consequently, we explored the network of interrelated genes using inBio Discover to delineate which diseases (DOID) are over-represented in the query gene list. Table 3 shows the top nine DOID disease annotations. This table indicates that the proteins involved in FES/FEP are enriched in a number of immune diseases, including immunodeficiency disease and autoimmune disease, and three gastro-intestinal diseases, namely intestinal disease, inflammatory bowel disease, and colitis. ESF Figure 6 shows the extended network and those four of the nine top annotations that were over-represented. The nodes of the extended PPI network involved in immune system and autoimmune disease and inflammatory bowel and intestinal disease are colored in red, orange, blue, and green, respectively. ESF Figure 7 shows all nodes (in red color) which were accumulated in intestinal disease (DOID:5295). ESF Figure 8 shows the results of an annotation analysis for all differentially expressed proteins, using R package ClusterProfiler and the GO list consisting of only leaf children in biological processes.

The PPI Network of FEP/FES
The first major finding of this study is that we were able to construct zero and firstorder PPI networks of FES/FEP which show high connectivity and some unexpected interactions, including the central role of BDNF. We found that the backbone of the FEP network is shaped by six genes (IL6, TNF, IL10, IL4, IFNG, and BDNF) and the backbone of FEP/FES by IL6, TNF, IL10, CCTNB1, CDH1, and BDNF. Moreover, our network analysis revealed two major protein communalities, namely one centered around immuneinflammatory genes and BDNF, and a second cluster around cell-cell junction genes (CTNNB1 and CDH1). Most importantly, BDNF appears to function as a switch between both communalities, as it shows interconnections with immune genes and CTNNB1 and the high-affinity receptor of BDNF, namely neurotrophic tyrosine kinase receptor 2 (NTRK2), which was allocated to the cell-cell junction cluster. Another, smaller, bridge between both communalities is CDH1, which is interconnected with CTNNB1 and with CSF2, IL4, TNF, IL10, and IL6.
Therefore, it appears that an immune response, as observed in FEP/FES, comprises not only IRS and CIRS components [4], but also different neurotrophic (BDNF, NTF4, and NGFR) and cell adhesion (CTNNB1 and CDH1) factors. By inference, the neurotoxic effects of the immune response may be counterbalanced by CIRS (Th-2 and Treg) [7], neuroprotective, and cell-cell adhesion genes. It is tempting to speculate that the participation of neuroprotective and cell-cell adhesion genes in this immune network is another conserved regulatory process protecting against detrimental neuro-immune effects.

Trigger Factors in FEP and FEP/FES
GO annotation and MCODE analysis revealed that the upregulated genes in FEP and the cluster-1 genes in FEP/FES were highly significantly associated with a response to a molecule of bacterial origin and a response to LPS. In fact, a response to LPS was the second most important GO term overexpressed in the upregulated proteins in FEP. Moreover, annotation analysis using a custom-made GO list with possible "responses to . . . " or "cellular responses to . . . " or a "defense response to . . . " showed that the most significant paths enriched in the network were responses to an organic substance (determined as cytokine or lipid), an external stimulus, response to other organisms, or an external biotic stimulus. Further analysis showed that a response to a bacterium and LPS were the most significant paths enriched in the FEP/FES network, whereas a response to a virus was less significant and did not appear in the GO and MCODE enrichment analysis. Moreover, an abiotic stimulus and fungal or parasite stimuli were not significantly enriched.
All in all, it appears that FEP and FEP/FES may be triggered or maintained by a response to LPS of Gram-negative bacteria. These findings corroborate a recent study which reported increased IgA/IgM responses to LPS of Gram-negative bacteria in FES and deficit schizophrenia, and a significant association between these indicants of bacterial load and the G-CoDe (general cognitive impairments) and symptom profiles, especially psychosis [7,18]. Such findings may indicate increased bacterial translocation in FES/FEP through leaky gut (increased gut permeability) with deficits in tight and adherens junctions and the vascular barrier, as well [7]. Interestingly, three of the top nine DOID annotations which were overrepresented in the PPI network of FEP/FES comprised intestinal annotations, including intestinal disease, inflammatory bowel disease, and colitis.

Upregulated Pathways, Molecular and Cellular Processes in FEP
The third major finding of this study is that the upregulated genes in FEP were enriched in key pathways and cellular functions that play a key role in immune-inflammatory signaling, namely the receptor signaling pathway via Janus kinases/signal transducer and activator of transcription proteins (JAK-STAT) pathway, the TNFR1-induced NFκB signaling pathway, and TNF-and death receptor signaling. Moreover, STAT3 and STAT6 were prominent nodes in the first order network, with STAT3 occupying a central position. Many cytokines (e.g., IL-2, IL-4, IL-6, IL-10, IL-12, IFN-γ) signal via the JAK-STAT pathway, thereby transactivating Janus kinases leading to translocation of STATs to the nucleus and upregulation of cytokine-modifiable genes [19]. The JAK-STAT pathway plays a key role in inflammation, cell death, cell division, and polarization of T cells; STAT3 is associated with autoimmune responses, and STAT6 plays a key role in M2 macrophage activation and Th-2 differentiation with production of IL-4, IL-5, IL-9, and IL-13 [19][20][21]. Recently, Sharma et al. [22] reported that a subgroup of schizophrenia patients showed increased STAT1 levels.
Our enrichment analyses show that TNF plays an important role in FEP and suggests that TNF-induced IκB kinase (IKK) activation with consequent NFκB translocation to the nucleus is one of the key processes in FEP. Moreover, we found that both NFκB1 and RELA (NFκB p65 unit or transcription factor p65) are the most prominent transcription factors in the FEP network. RELA plays a role in NFκB activation and translocation of the released NFκB complex to the nucleoplasm and contributes to DNA binding in the NFκB complex [23]. This NFκB heterodimeric RELA-NFKB1 complex functions as a transcriptional activator and plays a key role in gene expression of multiple cytokines [23]. In schizophrenia, it was shown that increased NFκB activity may contribute to cortical immune activation [24]. In a Japanese population, schizophrenia is associated with variants of the RELA gene and has a significant effect on pre-pulse inhibition [25]. Importantly, this pathway is not only induced by TNF, but also by other cytokines and LPS [26]. Previously, we have reviewed the many neurotoxic effects of increased TNF-α levels in schizophrenia, and especially in deficit schizophrenia [27]. Moreover, the frequency of the TNF2(A) allele, which affects plasma TNF levels, is significantly increased in schizophrenia and TNF2 homozygotes are detected in schizophrenia only [6].

Downregulated Pathways, Molecular and Cellular Processes in FEP
The fourth major finding is that the downregulated genes were enriched in the receptor protein tyrosine kinase (RTK) and neurotrophin/Trk receptor signaling pathways, cellular component morphogenesis, and production of miRNAs involved in gene silencing. RTK is a family of high-affinity cell surface receptors, including Trk, which is activated by neurotrophins such as BDNF [28]. Interestingly, RTKs regulate the threshold for macrophage activation, thereby promoting homeostasis and protecting tissues from inflammatory damage [29]. BDNF/Trk signaling contributes to axonal growth, axonal guidance, plasticity, dendritic arborization, synapse structure and formation and connections, neurogenesis, differentiation of new neurons and synapses, and axonal and dendritic sprouting [30]. Moreover, BDNF may have anti-inflammatory and anti-apoptopic effects via modulation of MyD88/NFκB and PI3K/AKT-signaling pathways [31]. Following bacterial infection, BDNF pretreatment reduces the expression of TNF-α, IL-16, IL-1β, and the NFκB pathways, and increases IL-10 and Trk expression [31]. LPS-associated inflammation alters BDNF/TrkB signaling in the hippocampus, nucleus accumbens, and prefrontal cortex in association with the onset of depressive behaviors [32]. A recent meta-analysis showed that schizophrenia is associated with reduced BDNF with a medium effect size (Hedges g = −0.458, p < 0.004), and that these effects are not influenced by the drug state of the patients [33]. Moreover, the interaction between BDNF and NTRK2 gene polymorphisms may increase susceptibility to paranoid schizophrenia [34].
We found that the downregulated DISC1 gene is enriched in various GO biological terms indicating neurogenesis, axonogenesis, and axon extension. DISC1 is now established as a risk factor for schizophrenia and other major psychiatric illness and is involved in aspects of adult progenitor proliferation, neurogenesis, neurite outgrowth, cytoskeletal modulation, signal transduction, and CTNNB1 abundance [35,36]. It is interesting to note that inflammatory signaling via TLR3 is accompanied by impairments in dendritic spine and growth via effects on MYD88 and, consequently, DISC1 expression [37].
Another highly significant GO path enriched in the downregulated network was the "production of miRNAs involved in gene silencing by miRNA", which may, at least in part, be ascribed to downregulation of DROSHA, which is a ribonuclease (RNase) III family enzyme and plays a key role in miRNA maturation [38]. In mammals, the miRNA network comprises 5000-10,000 miRNAs which regulate the expression of 60% of the protein-coding genes through translational silencing and mRNA destabilization [39,40]. Importantly, miRNA regulate the adaptive and innate immune response and act as fine-tuning regulators, preventing an overzealous inflammatory response and thereby maintaining homeostasis [40]. Many of the miRNAs which are associated with schizophrenia phenotypes [41,42] display immune regulatory effects. For example, miR-9 exerts a negative feedback on NFκB and is dysregulated in neural progenitor cells of schizophrenia patients [43]; miR-132 inhibits inflammation signaling (via acetylcholine, STAT3, and NFKB) and is dysregulated in schizophrenia [44]; miR-146 inhibits inflammatory responses and is downregulated in monocytes of postpartum psychosis patients [45]; and miR-149 inhibits LPS-induced inflammation (via STAT3, NFκB, TNF, IL-6) and is a candidate biomarker of psychiatric disease including bipolar disorders [46].

Pathways, Molecular and Cellular Processes in FEP/FES
The enrichment and annotation analysis revealed other important drug targets in FEP/FES. Firstly, we found that a neuroinflammatory response was enriched in the seed gene FEP/FES list, whilst MCODE showed that a cytokine/chemokine complex of IFN-γ, IL-6, IL-12A, CCL3, IL-4, and IL-13 was strongly associated with microglial cell activation and tyrosine phosphorylation of STAT proteins. These results extend the findings that schizophrenia is accompanied by microglial activation [47]. Moreover, the upregulated genes in FEP/FES were enriched in "the positive regulation of gliogenesis". In adulthood, gliogenesis is maintained to renew oligodendrocytes; however, following inflammatory disease and injuries, gliogenesis becomes more active (reactive astroytosis or astrogliosis) and may have negative consequences, thereby contributing to immune-inflammatory responses and altering the balance between neurogenesis and gliogenesis [48].
Secondly, WikiPathway and PANTHER enrichment analysis revealed that the upregulated genes were strongly associated with the TLR signaling (especially TLR4) and tolerance pathways. These findings extend those of previous publications indicating activation of the TLR4 proinflammatory pathway in schizophrenia [49].
CTNNB1 is a component of the Wnt/β-catenin signaling pathway and the E-cadherincatenin adhesion complex, which play a key role in epithelial integrity and tissue architecture maintenance [50,51]. The Wnt/catenin pathway is strongly involved in neurogenesis, axonal spreading and branching, connectivity between pre-and post-synaptic neuronal regions, regulation of synaptic functions and modeling of synaptic structures, modulation of excitatory synaptic transmission, LTP, and post-synaptic protein assembly [52]. The Wnt/β-catenin signaling pathway also regulates immune-inflammatory responses and T-cell-inflammation [53][54][55]. Inflammatory responses, due to infections with pathogenic bacteria, may affect the Wnt/β-catenin signaling pathway [55] and the E-cadherin-catenin adhesion complex [56]. For example, in inflammatory bowel disease, impairments in the latter complex are affected by the inflammatory milieu and may cause dysregulations of the actin cytoskeleton leading to aberrations in intracellular signaling and transcriptional regulation [57]. Moreover, the Wnt/catenin pathway may regulate BDNF expression while these two pathways may have common effector actions [58], and BDNF polymorphisms are associated with changes in the Wnt/β-catenin pathway [59]. In hippocampal neurons, BDNF-disruption of cadherin-β-catenin complexes is associated with increased synapse density [60].
Disorders in the Wnt/catenin pathway were previously described in FES and schizophrenia [7,61,62] and alterations in E-cadherin and beta-catenin levels in FES are strongly associated with increased bacterial translocation [7]. This more generalized disorder in paracellular and cell-cell junctions in FEP/FES may, at least in part, be related with the increased frequency of the Hp2 allele and the Hp2.2 genotype (prehaptoglobin-2 or zonulin) [63] and increased zonulin levels [64]. Moreover, mutations in the CTNNB1 gene (c.1943 A>G) are associated with schizophrenia [65], whilst CTNNB1 KO mice display anxiety behaviors and CTNNB1 KO in paraventricular interneurons accompanied by impairments in social interactions, repetitive behaviors, and object recognition [66]. Association, candidate gene, and genomewide association studies show that cadherins may be involved in the pathophysiology of schizophrenia [67,68].
Fourth, exploration of transcriptional regulation showed that SP1, NFκB1, and RELA were the most prominent transcription factors in the FEP/FES network. SP1 or specificity protein 1 (or transcription factor Sp1) is a ubiquitously expressed transcription factor, which regulates the expression of a variety of house-keeping and tissue-restricted genes frequently involved in immune responses, response to DNA damage, and apoptosis [69]. This explains that SP1 is associated with the pathophysiology of some neurodegenerative and neuroinflammatory disorders, including Alzheimer's and Huntington's disease and multiple sclerosis [69]. Interestingly, neurons have a decreased capacity of activating NFκB, but κB cis elements may bind to SP1 [70] and SP1 interacts with RELA to form a complex [71]. In cortical neurons, SP1 is an oxidatively-induced transcription factor which regulates neuronal survival [72]. In Huntington's disease, pathogenic SP1 cascades cause repression of neuronal genes [73]. Interestingly, both NFκB and SP1 modulate antimicrobial activity against Gram-negative bacteria [74].
Fifth, biological GO term classifications showed that complement factors were enriched not only in microglial activation and humoral immune responses, but also in synapse pruning and organization. The involvement of complement in FEP/FES agrees with previous findings showing increased plasma C3C and C4 and CSF and brain C4 complement factors in schizophrenia [2,75,76]. Moreover, increased C1qA, C3, and C4 transcripts were reported to be associated with microglial activation in the midbrain of schizophrenia patients [77]. It should be added that FES is accompanied by increased formation of the IgA C1qA CIC, which may have detrimental effects in its own right [7]. Nevertheless, the complement pathways established in FEP/FES are not a key component of the interactome and, therefore, should be considered as secondary phenomena. In addition, also PON1 was not a key component of the interactome and showed only a few interactions. Other CIRS components could not be included, such as natural IgM responses to oxidatively specific epitopes [7]. Both PON1 and natural IgM are first line innate immune defenses against bacterial infections and display strong anti-inflammatory and antioxidant properties, explaining that impairments in both systems contribute to FEP/FES [7]. Figure 9 shows a summary of the enrichment and annotation analysis reported in our study. All in all, the analyses suggest that increased bacterial translocation, in particular that of Gram-negative bacteria, is the most probable trigger and maintaining factor of FEP and FEP/FES, causing an immune-inflammatory milieu with involvement of TNF, NFκB/RELA, SP1, and JAK-STAT signaling, including tyrosine phosphorylation of STAT proteins, and death receptor and TLR4 signaling. A specific cytokine/chemokine complex consisting of IFN-γ, IL-6, IL-12A, CCL3, IL-4, and IL-13 is strongly associated with microglial activation and neuroinflammation, and this upregulated immune network additionally positively regulates gliogenesis. FEP/FES is also accompanied by lowered immune regulation and thus protection against an overzealous inflammatory process due to dysfunctional miRNA maturation, and deficits in neurotrophin/Trk, RTK and Wnt/catenin signaling. regulation and thus protection against an overzealous inflammatory process due to dysfunctional miRNA maturation, and deficits in neurotrophin/Trk, RTK and Wnt/catenin signaling. Figure 9. Summary of the findings of the enrichment/annotation analysis using upregulated and downregulated genes in first episode psychosis and schizophrenia.

Conclusions
By inference, individuals with reduced CIRS [4], including miRNA maturation and deficits in neurotrophin/Trk, RTK, and Wnt/catenin signaling, neuroprotection (including neurotrophin/Trk and Wnt/catenin signaling, lowered DISC1 expression, and interactions between reduced BDNF, CDH1, CTNNB, and DISC1 may be at an increased risk to develop FEP and FES as a consequence of an immune response following bacterial infections and the ensuing neuro-immune toxicity. The latter may cause a multitude of neuronal dysfunctions and cognitive impairments and schizophrenia behaviors [7], as shown in Figure 9. Interestingly, almost all pathways or molecular patterns enriched in the interactome of FEP/FES are directly (NFκB/RELA, TLRs, TNF, adherens junctions and cell-cell junctions, complement activation with C1QA CIC formation, neurotrophin/Trk pathway, Wnt signaling, microglial activation, gliogenesis, neurotoxicity) or indirectly (JAK-STAT, DISC1, SP1) affected by LPS. As such, the increased LPS levels in FES may maintain peripheral and central immune activation resulting in neurotoxic effects on central neuronal circuitry, neurogenesis, and synapse functions [7]. Disorders in cell-cell junction organization may contribute to the maintenance of increased bacterial translocation and activation of complement cascades with increased IgA C1QA CIC formation, which may further fuel the immune response [7]. It should be added that LPS-induced maternal immune activation may have increased the vulnerability to such immune hits by inducing neurodevelopmental disorders with sensitized immune-inflammatory pathways (including in TLR4 and inflammasome) and lowered neuroprotection [78,79].
Future research in FEP should examine the microbiome, gut dysbiosis, and consequent bacterial translocation due to leaky gut as well as alterations in the microbiota-gut- Figure 9. Summary of the findings of the enrichment/annotation analysis using upregulated and downregulated genes in first episode psychosis and schizophrenia.
By inference, individuals with reduced CIRS [4], including miRNA maturation and deficits in neurotrophin/Trk, RTK, and Wnt/catenin signaling, neuroprotection (including neurotrophin/Trk and Wnt/catenin signaling, lowered DISC1 expression, and interactions between reduced BDNF, CDH1, CTNNB, and DISC1 may be at an increased risk to develop FEP and FES as a consequence of an immune response following bacterial infections and the ensuing neuro-immune toxicity. The latter may cause a multitude of neuronal dysfunctions and cognitive impairments and schizophrenia behaviors [7], as shown in Figure 9. Interestingly, almost all pathways or molecular patterns enriched in the interactome of FEP/FES are directly (NFκB/RELA, TLRs, TNF, adherens junctions and cell-cell junctions, complement activation with C1QA CIC formation, neurotrophin/Trk pathway, Wnt signaling, microglial activation, gliogenesis, neurotoxicity) or indirectly (JAK-STAT, DISC1, SP1) affected by LPS. As such, the increased LPS levels in FES may maintain peripheral and central immune activation resulting in neurotoxic effects on central neuronal circuitry, neurogenesis, and synapse functions [7]. Disorders in cell-cell junction organization may contribute to the maintenance of increased bacterial translocation and activation of complement cascades with increased IgA C1QA CIC formation, which may further fuel the immune response [7]. It should be added that LPS-induced maternal immune activation may have increased the vulnerability to such immune hits by inducing neurodevelopmental disorders with sensitized immune-inflammatory pathways (including in TLR4 and inflammasome) and lowered neuroprotection [78,79].
Future research in FEP should examine the microbiome, gut dysbiosis, and consequent bacterial translocation due to leaky gut as well as alterations in the microbiota-gut-brain axis [18], which may ultimately activate or dysregulate the many pathways discussed in the current review. In addition, it would be most informative to examine differences in the gut-brain axis between FEP, deficit schizophrenia, and other major psychiatric and psychosomatic disorders. There is some evidence that increased leaky gut and bacterial translocation and their impact on neuro-immune pathways may be a transdiagnostic phenomenon [80] which is detected in affective disorders [81] and Myalgic Encephalomyeli-tis [82]. Future research should examine the effects of new or repurposing drugs and natural anti-inflammatory and antioxidant compounds on the gut-brain axis and intracellular pathways in FEP (see Figure 9). For example, one could aim to attenuate the primary bacterial translocation in FEP by using minocycline, which has, additionally, anti-inflammatory and antioxidant capacities [83]. One could also aim to target NFκB using specific kinase inhibitors with or without natural NFκB inhibitors, such as curcuma [84]. However, a novel approach to treating FEP and, thus, preventing schizophrenia, its worsening course, and the onset of deficit syndrome appears to require a highly complex approach, as deficiencies in the CIRS and neurotrophic factors appear to be critical. As such, combination therapies with immunotherapies (including intravenous immunoglobulins, which have some efficacy in schizophrenia in case reports [85]) and neurotrophic factors [84] appear to be warranted.

Supplementary Materials:
The following are available online at https://www.mdpi.com/article/ 10.3390/cells10112929/s1, Figure S1: Heatmap (top-10) of enriched GO cellular component terms accumulated in the differently expressed proteins in first episode psychosis and schizophrenia., Figure S2: Heatmap (top-10) of enriched Reactome terms over-represented in the differently expressed proteins in first episode psychosis and schizophrenia, Figure S3: Heatmap (top-10) of enriched Panther Pathways that were over-represented in the differently expressed proteins in first episode psychosis and schizophrenia. Figure S4: Heatmap (top-10) of enriched WikiPathways that were over-represented in the differently expressed proteins in first episode psychosis and schizophrenia. Figure S5: Protein-protein interaction network (built with OmicsNet and IntAct) of first episode psychosis and schizophrenia. The top-3 transcriptional factors are shown as blue spots (SP1, NFκB1, and RELA). Proteins in green and bundles in pink. Figure S6. Protein-protein interaction network (built with inBio Discover) of differently expressed proteins in first episode psychosis and schizophrenia. Shown are 4 top annotations that are over-represented in the gene list, namely: red color: immune system disease (DOID:2914); blue color: inflammatory bowel disease (DOID:0050589); green color: intestinal disease (DOID:5295); yellow color: autoimmune disease (DOID:417) Figure S7: Differently expressed proteins in the extended protein-protein interaction network that accumulated in intestinal disease (DOID:5295). Figure S8. GO functional enrichment analysis of all differently expressed proteins in first episode psychosis and schizophrenia using R package ClusterProfiler. Only annotated GO term leaves are shown in this figure. The x-axis shows the gene ratio and the y-axis the annotated GO biological terms. The size of the dots is proportional to the gene number; p-values of all GO terms are colored as indicated in the figure. Institutional Review Board Statement: Not applicable. This study is a secondary data analysis on existing data (which were openly available to the public before the initiation of the study) using open, deidentified and non-coded datasets and, therefore, this is non-human subjects research, which is not subject to IRB approval.