A Cross-Species Systems Genetics Analysis Links APBB1IP as a Candidate for Schizophrenia and Prepulse Inhibition

Background

Schizophrenia is a severe mental disorder characterized by hallucinations, delusions, social withdrawal, and cognitive deficits such as poor executive function and working memory (Burmeister et al., 2008). It is a major cause of morbidity contributing significantly to the global burden of disease (Lopez et al., 2006) and affects about 1% of the population. While onset typically occurs after the period of brain development that follows puberty in late adolescence or early adult life (Kirkbride et al., 2012), there is now strong clinical and epidemiological evidence that schizophrenia reflects a disturbance of neurodevelopment (Owen et al., 2016). Both genetic and environmental risk factors for schizophrenia have been identified that are linked to immune function, e.g., (Murray et al., 2019), an idea that was first postulated more than a century ago (for a review see Khandaker et al., 2015). The brain is no longer considered an immunologically privileged site, isolated from the immune system via the blood-brain barrier (Carson et al., 2006). Microglia derived from the hemopoietic system beyond the central nervous system (CNS), make up about 10% of the brain cell mass, equal to that of neurons (Hanisch and Kettenmann, 2007). Animal models provide evidence that systemic inflammatory cytokines, e.g., interleukin 6 (IL6) or tumor necrosis factor-α (TNF-α), communicate with the brain, affecting synaptic plasticity, cortisol concentrations, neurotransmitters, activation of the hypothalamic-pituitary-adrenal axis, and increased oxidative stress in the brain (Dantzer et al., 2008; Khandaker et al., 2015; Murray et al., 2019). These effects could contribute to changes in mood, cognition and behavior, supported by epidemiological studies that suggest links between schizophrenia and systemic inflammation (Brown and Derkits, 2009; Miller et al., 2011; Khandaker et al., 2014).

In order to understand the etiology and biology of schizophrenia it is essential to identify the underlying causal variants. This will allow us to link the disorder to specific proteins and pathways, opening the door for therapeutic targets, and a more accurate prediction of genetic predisposition.

Twin studies estimate broad heritability of schizophrenia at ~81% (Sullivan et al., 2003; Craddock et al., 2005), and heritability of 31–44% is reported for nuclear and extended families (Light et al., 2014) indicating a substantial genetic component; however the role of genetic variation in the etiology of schizophrenia is complex. A recent, large-scale genome-wide association study (GWAS) identified 108 independent genomic loci exhibiting genome-wide significance associated with schizophrenia, suggesting that schizophrenia is highly polygenic, with many single nucleotide polymorphisms (SNPs), each of small effect (Schizophrenia Working Group of the Psychiatric Genomics Consortium, 2014). But GWAS have been unable to account for the large portion of the heritability shown by twin and family studies. To reduce type 1 error and account for association testing of ~1 million SNPs in the human genome, multiple testing corrected significance threshold of p-value ≤5 × 10⁻⁸ is required (Barsh et al., 2012). This stringent significance threshold leads to many false-negative associations, potentially preventing the discovery of additional loci that could potentially offer important insights. Despite the modest statistical power due to high multiple testing correction, loci in human GWAS are defined with high precision, potentially to the individual SNP level. In contrast, mouse model linkage studies are often statistically powered enough to detect genetic effects but have lower resolution, producing loci that include tens or hundreds of genes (e.g., Ackert-Bicknell et al., 2010; Hager et al., 2012). Mouse systems have been extensively used to as a gateway to experimental precision medicine aiming to establish the genetic basis of disease traits due to the high degree of homology between mouse and human genomes, and because of the possibility of precisely manipulating the mouse genome and experimental variables (Peters et al., 2007). For example, recombinant inbred (RI) mouse models, such as the BXD family used here, have the advantage that the same genome can be replicated an unlimited number of times, allowing both invasive tissue collection, and for tissue to be collected along a developmental timeline. To capitalize on the advantages of both GWAS and genetic mouse models we and others have developed methods that combine data from both systems, gaining power from mouse model studies and precision from human GWAS (Himes et al., 2013; Ashbrook et al., 2014, 2015; Jing et al., 2016; Mignogna et al., 2018).

Prepulse inhibition (PPI) of the startle response is a form of sensorimotor gating that occurs when a weak, subthreshold “prestimulus” presented 30–500 ms prior to an intense startling stimulus inhibits the startle response (Swerdlow et al., 2000). PPI is a highly conserved phenomenon and observed across sensory systems, but acoustic PPI is frequently used for investigations and is usually measured from the eyeblink response (Swerdlow et al., 1999). PPI deficits were first observed in schizophrenia patients over 40 years ago (for a review see Braff et al., 2001), but are also evident in their unaffected first degree relatives (Cadenhead et al., 2000), and patients with schizotypal personality disorder (Cadenhead et al., 1993), reflecting a genetically transmitted susceptibility to sensorimotor gating deficits. In fact, PPI heritability has been estimated at 32% (Greenwood et al., 2007), similar to the schizophrenia heritability estimates of 31% and 44% for nuclear and extended families. This suggests similar heritabilities for both the disease and the endophenotype (Light et al., 2014). In terms of genetic analysis, endophenotypes are heritable, measurable biomarkers that are potentially less genetically complex (Gottesman and Gould, 2003). PPI deficits are not unique to schizophrenia however, and have been observed in a number of psychiatric and neurological disorders including bipolar disorder, autism, obsessive-compulsive disorder, Tourette’s syndrome and Huntington’s disease (Kohl et al., 2013), but are most widely replicated in schizophrenia patients (Braff et al., 2001). Thus PPI meets the criteria for a viable endophenotype, outlined by Turetsky et al. (2007), which specifies that the endophenotype must be associated with a disorder, heritable, evident in unaffected family members, easily measured, with high test, re-test reliability. This, combined with the observation that inbred mouse strains demonstrate considerable phenotypic variation in PPI (for a review see Petryshen et al., 2005), renders PPI a behavioral phenotype with rich potential for discovery of genes in genetic mouse models that confer risk or resistance to schizophrenia.

Animal models have proven invaluable in genetics research and previous research has demonstrated the effectiveness of using a cross-species approach to identify genes underlying specific traits (Himes et al., 2013; Ashbrook et al., 2014, 2015; Jing et al., 2016; Mignogna et al., 2018). The benefit of this approach is that it allows the investigation of phenotypes without requiring the quantification of experimental perturbations on the biological system. Our study uses a cross-species approach to identify genes associated with PPI regulation, which thus may contribute to the PPI deficits seen in schizophrenia. We use data obtained from populations that have undergone segregation analysis for large numbers of common sequence variants and associated differences in phenotypes (Ashbrook et al., 2014, 2015). First, we used data from PPI in the most deeply phenotyped mammalian model system, the RI mouse panel BXD (Wu et al., 2014), to identify quantitative trait loci (QTL). This identifies areas of the genome containing genetic variants that influence PPI in mice and therefore may be influencing aspects of schizophrenia in humans. We then went on to investigate these genomic regions in a large human GWAS of schizophrenia, the SNP summaries of which are available in an online database (Schizophrenia Working Group of the Psychiatric Genomics Consortium, 2014). We identified one gene (amyloid-β precursor protein-binding family B member 1 interacting protein; APBB1IP) that was significantly associated with PPI in mice and the risk of schizophrenia in humans at the genome-wide level. Our systems genetics approach demonstrates that APBB1IP coexpresses with several other genes related to schizophrenia in several brain regions, providing a potential mechanism of action.

Materials and Methods

Mouse and Human Data

The BXD RI panel was derived by crossing the parental B6 and D2 strains followed by >20 generations of inbreeding, and currently consists of ~150 commercially available and densely genotyped strains (Williams et al., 2001; Peirce et al., 2004; Ashbrook et al., 2019). The BXDs are the most deeply phenotyped mammalian experimental population, and have been used as a genetic reference panel for over 40 years (Morse et al., 1979). The online resource www.genenetwork.org contains BXD genotypes and software for rapid QTL mapping (GeneNetwork). This resource also contains over 7000 BXD phenotypes uploaded by researchers that can be compared and assessed for genetic correlations and common QTL. These phenotypes include strain mRNA expression levels for many tissues that can be utilized for systems genetics analysis and prioritization of positional candidate genes in behavioral and physiological QTL (Bubier and Chesler, 2012). For our analysis, we used data on PPI startle (Bagley, 2017), available in GeneNetwork (Sloan et al., 2016). We selected any trait with at least one significant QTL (p ≤ 0.05) at the genome-wide level. Significance was calculated within GeneNetwork, using 5,000 permutations of trait values and genotypes. A p-value of 0.05 is defined as a likelihood ratio statistic (LRS) score greater than 95% of the permutated datasets (Wang et al., 2003).

Human GWAS data for schizophrenia were obtained from the Psychiatric Genomics Consortium (PGC; Schizophrenia Working Group of the Psychiatric Genomics Consortium, 2014), containing 36,989 cases and 113,075 controls. The Knowledge-Based Mining System for Genome-wide Genetic Studies (KGG¹; version 4.0) was used to convert the SNP p-values to gene p-values using the GATES method, an extension of a Simes test that integrates functional information and association evidence to rapidly combine the p-values of the SNPs within a gene, to produce valid gene-based p-values (Li et al., 2011). The GATES test is advantageous because it requires only summary GWAS and ancestral population linkage disequilibrium (LD) data; it provides effective control of the type 1 error rate regardless of gene size and LD pattern among markers, and permutation or simulation is not required to evaluate empirical significance (Li et al., 2011). Gene locations were from the Hg19 genome build. We used the complete set of data from the PGC GWAS for schizophrenia (PGC-GWAS; ckqny.scz2snpres.gz), and ancestry matched genotypes from phase 3 of 1,000 Genomes Project to adjust LD between variants (Li et al., 2014; 1 kg.phase3.v5.shapeit2.eur.hg19.tar). For each gene in Hg19 a GATES p-value was calculated for association with schizophrenia.

For the joint mouse-human analysis, the human homologs of genes within the mouse QTL were identified using Ensembl version 94 (Zerbino et al., 2017). To assess if any particular gene is associated with both PPI in mouse and schizophrenia in humans, we examined the GATES derived p-values of human homologs of genes within the significant BXD QTL. The human GWAS significance values were Bonferroni corrected for multiple comparisons using the total number of homologous genes (Ashbrook et al., 2014, 2015).

Gene Ontology

Genease (Genease) and WebGestalt (Zhang et al., 2005) tools were used for Gene Ontology (GO) annotation. For microarray-based transcriptome data used in enrichment analyses in WebGestalt, the array IDs were used, and the total probes on the array were used as background. For RNA-seq data, ENTREZ IDs were used, and all protein-coding genes were used as background.

Results

We began by identifying QTL for PPI in the BXD RI family. Next, we investigated the homologous genomic regions in human schizophrenia GWAS data and were able to identify one novel candidate gene, APBB1IP (also referred to as PREL1, RARP1, and RIAM; Entrez 54518).

QTL for Schizophrenia-Related Phenotypes in Mice

PPI deficits are associated with schizophrenia (Braff et al., 2001), and we, therefore, investigated genetic variation in this trait. Two PPI traits met our criteria in GeneNetwork: GN18755 (trait description: Acoustic startle response, PPI, change in startle response with 74 dB prepulse from response with no prepulse in females, unstressed controls), which has one significant QTL with a maximum LRS of 21.4 on chromosome 2 (Figure 1). The second trait (GN11428: Acoustic startle response, PPI at 85 dB for males; Philip et al., 2010), has one significant QTL on chromosome 19 with a maximum LRS of 19.2 (Figure 2).

FIGURE 1

Figure 1. Quantitative trait loci (QTL) map showing significant QTL on chromosome 2 between 22 and 24 Mbp, prepulse inhibition (PPI) trait (GeneNetwork ID 18755). The red horizontal line towards the top of the figures illustrates the genome-wide significance threshold, which is p-value ≤ 0.05 genome-wide corrected [significant likelihood ratio statistic (LRS) = 20.54]. The lower, blue line indicates the significance of the trait at each position along the genomic regions. The lower, green and red line shows the negative additive coefficient. Regions where DBA/2J alleles increase trait values are shown in green, and regions where the C57BL/6J alleles increase trait values are shown in red, with the scale in green on the left. The multi-colored blocks at the top of the figure represent genes, showing where they fall in the genomic region, and their length. The position of segregating single nucleotide polymorphisms (SNPs) in the BXD family are shown by the orange Seismograph track at the bottom of each map. Adapted from genenetwork.org.

FIGURE 2

Figure 2. QTL map showing a significant QTL on chromosome 19 between 31 and 33 Mbp, PPI trait (GeneNetwork ID 11428). The red horizontal line towards the top of the figures illustrates the genome-wide significance threshold, which is p-value ≤ 0.05 genome-wide corrected (significant LRS = 17.30). The lower, blue line indicates the significance of the trait at each position along the genomic regions. The lower, green and red line shows the negative additive coefficient. Regions where DBA/2J alleles increase trait values are shown in green, and regions where the C57BL/6J alleles increase trait values are shown in red, with the scale in green on the left. The multi-colored blocks at the top of the figure represent genes, showing where they fall in the genomic region, and their length. The position of segregating SNPs in the BXD family are shown by the orange seismograph track at the bottom of each map. Adapted from genenetwork.org.

To further characterize the two QTL we used Genome-wide Efficient Mixed Model Association (GEMMA) software (Zhou and Stephens, 2012) and Leave One Chromosome Out (LOCO) in GeneNetwork2². For GN18755 the most significant marker was rs13476379 at chromosome 2: 22.82 Mb. This marker has a −log(p) value of 6.16, compared to the surrounding marker values of 4.13. This marker lies within the Apbb1ip gene. For GN11428, the most significant marker was rs6197068 at chr19: 32.74 Mb. This marker has a −log(p) value of 5.69. This lies between the genes Atad1 and Pten.

The QTL on proximal chromosome 2 (22.62–22.86 Mbp) contains four genes with human homologs, and the QTL on proximal chromosome 19 (31.95–32.75 Mbp) contains eight genes with human homologs, all of which are on chromosome 10 (Table 1). We then examined the PGC Schizophrenia GWAS data to identify variants within the homologous region significantly associated with schizophrenia. Only one of the twelve genes, APBB1IP, has a nominally significant p-value (p = 0.0117).

TABLE 1

Table 1. Human homologs (genome assembly: GRCh38.p12) of the 12 mouse genes (Genome assembly: GRCm38.p6) contained in the genome-wide significant quantitative trait loci (QTL) of prepulse inhibition (PPI) phenotypic traits 18755 and 11428.

Next, we checked for associations in an independent dataset, the UKBiobank collection using PheWeb (PheWeb). This shows that two variants in APBB1IP are significantly associated with self-reported schizophrenia (UKBB ID 20002_1289). The first is rs570949810 at chr10: 26,823,566 (T/G, p = 1.8e-9) and the second is rs554788824 at chr10:26,818,850 (G/A, p = 8e-9). Both variants are also associated with usage of two antipsychotics used to treat schizophrenia, trifluoperazine (p = 1.1e-6 and 5.7e-7, respectively) and amisulpride (p = 7.2e-7 and 1.3e-5, respectively). This further supports APBB1IP as a schizophrenia associated gene.

These results thus provide strong evidence that APBB1IP influences schizophrenia-related traits in both mouse and human. We can next take advantage of massive gene expression datasets to identify plausible mechanisms by which APBB1IP may be linked to schizophrenia.

Human Gene Expression

To identify the regions in which our candidate gene is expressed we used the Genotype-Tissue Expression (GTEx) portal³. The GTEx portal indicates that APBB1IP has high expression in whole blood, in lymphocytes and in the spleen (Figure 3). This links APBB1IP to the immune system, confirming previous findings (Zhang and Wang, 2012; Su et al., 2015; Patsoukis et al., 2017). To expand on this, we used GeneNetwork to query the GTEXv5 Human Whole Blood (GN770) data and searched for genes that are significantly coexpressed with APBB1IP in whole blood. There are 3464 genes which show coexpression (r ≥ 0.5 or ≤ −0.5, p ≤ 0.05; Supplementary Table S1), and enrichment analysis of these genes shows clear links to the immune system (e.g., regulation of immune response, GO:0050776, FDR = 5.19e-10; Innate Immune System, R-HSA-168249, FDR < 1e-20; Supplementary Figure S1, Supplementary Table S2). Interestingly, the top two enriched drug targets interact with ARP2/3, which has previously been linked to schizophrenia (Gokhale et al., 2016; Yan et al., 2016; FDR = 7.96e-05 for both; Supplementary Table S3).

FIGURE 3

Figure 3. Gene Expression of APBB1IP in Human Tissues taken from www.gtexportal.org (Lonsdale et al., 2013). Expression values are shown in log-transformed Transcripts Per Million (TPM). These TPM values were calculated from a gene model with isoforms collapsed to a single gene. We did not apply any other normalization steps. Violin plots are shown as median and 25th and 75th percentiles. Outliers above or below 1.5 times the interquartile range are displayed as points. Adapted from www.gtexportal.org.

However, GTEx only uses adult tissue (20–71 years old; Lonsdale et al., 2013), and schizophrenia is known to involve neurodevelopmental deficits, which would occur before this age range. As such, we examined the Allen Brain Atlas Human Brain (Hawrylycz et al., 2012) and Developing Human Brain (Miller et al., 2014) data. The Developing Brain data suggests that APBB1IP expression increases postnatally in a number of brain areas (Supplementary Table S4) and that in the adult brain there is higher expression in the cerebellum and dentate gyrus, compared with other brain areas. Both of these structures have been previously linked to schizophrenia (Yuan et al., 2015; Mothersill et al., 2016; Nakao et al., 2017). Given these links, we then repeated the above analysis, using the cerebellar and hippocampal GTEXv5 data in GeneNetwork. In the cerebellum, 282 genes coexpress with APBB1IP (r ≥ 0.5 or ≤ −0.5, p < 0.05; Supplementary Table S5). Again, we see enrichment for immune-related genes (e.g., regulation of immune system process, GO:0002682, FDR = 1.9e-12; TYROBP Causal Network, WP3945, FDR < 1e-20; Supplementary Figure S2; Supplementary Table S6). This is of further interest as TYROBP has also been linked to schizophrenia (de Baumont et al., 2015).

Finally, we used GeneFriends to examine the potential roles of APBB1IP, which is agnostic to tissue type. Here, we found that 785 genes commonly coexpress with APBB1IP (r ≥ 0.5 or ≤ −0.5; Supplementary Table S7). The enrichment analysis suggests strong involvement with immune-related processes (Supplementary Figure S3).

We have thus strong evidence that, in humans, APBB1IP is linked to immune function, and is expressed in the developing brain. In the next step, we use murine gene expression data to explore the cross-species translatability of this finding, which may enable future work to develop interventions that can be tested in mouse models.

Mouse Gene Expression

We first repeated the same GeneFriends analysis as above but using mouse data. This showed 179 commonly coexpressed genes (r ≥ 0.5 or ≤ −0.5; Supplementary Table S8), showing that genes that commonly coexpress with Apbb1ip are involved in immune response (Supplementary Figure S4). Next, we sought to establish if Apbb1ip is expressed in homologous tissues to APBB1IP. The Allen Brain Atlas shows expression of Apbb1ip across a number of brain regions in young adult mice (Figures 4, 5), analogous to the human brain. In particular, we again see the expression in the dendrite gyrus, as we saw in humans. Changes in gene expression in the dendrite gyrus in schizophrenia have been reported by others (Tavitian et al., 2019).

FIGURE 4

Figure 4. Bar graph showing expression of Apbb1ip across different brain regions in young adult mice, adapted from Allen Brain Atlas (Brain-Map; Lein et al., 2007; Sunkin et al., 2013; Zaldivar and Krichmar, 2014). The data is from Apbb1ip probe RP_050111_02_C07, with data taken from a sagittal cut from a 56-day old male C57BL/6J.

FIGURE 5

Figure 5. Image of expression of Apbb1ip in the young adult mouse brain, taken from the Allen Brain Atlas (Brain-Map; Lein et al., 2007; Sunkin et al., 2013; Zaldivar and Krichmar, 2014). The data is from Apbb1ip probe RP_050111_02_C07, with data taken from a sagittal cut from a 56-day old male C57BL/6J.

In the next step of our analysis, we analyzed deep gene expression data available for the BXD family. This enables us to map expression QTL (eQTL) to identify regions of the genome that are regulated by our gene of interest. Second, we can examine gene co-expression to find groups of genes which are under similar regulatory control and are therefore likely to be involved in the same biological processes. Thirdly, we can detect low levels of gene expression, or gene expression that only occurs under specific conditions.

Because of the high expression of APBB1IP in human lymphocytes, we first used two lymphocyte cell populations from the BXD family, T-helper cells and T-regulatory cells. In both cell types there is a cis-eQTL, indicating that a variant in Apbb1ip regulates expression of this gene, although this only reaches genome-wide significance in the T-helper cells (LOD = 7.8, 1.5 LOD confidence interval chromosome 2:17.95–23.77 Mb), and reaches chromosome-wide significance in T-regulatory cells. We next looked for the probes that correlated significantly with the expression of Apbb1ip (p < 0.05) in both of these datasets, to find a co-regulated network. In T-helper cells, 3,623 probes (3,249 genes; Supplementary Table S9) were co-expressed with Apbb1ip and in T-regulatory cells 2,565 (2,367 genes; Supplementary Table S10). Of these, 507 genes were co-expressed in both cell types, and gene-set enrichment was carried out on these 507 using WebGestalt. The most highly enriched annotations were cell migration (GO:0016477, FDR = 6.73e-03) and cell motility (GO:0048870, FDR = 6.73e-03). Interestingly, cell migration has been found to be dysregulated in schizophrenia patient-derived cells (Tee et al., 2017).

We next moved on to investigate more broadly all tissues with transcriptome data in the BXD. For each dataset, we searched for genome-wide significant eQTL, and recorded their position. To reduce false positives, we concentrated only on that eQTL which occurred in more than one dataset. We found the same cis-eQTL location that we detected in T-helper cells in a number of tissues (pituitary, aged hippocampus, midbrain, spleen, gastrointestinal tract, hypothalamus), showing that a variant near Apbb1ip influences its expression across tissues. This does not mean that this variant is not expressed in tissues other than those given above, but simply that no effect could be detected in that particular dataset. For each of these issues with a confirmed Apbb1ip cis-eQTL, we carried out a coexpression analysis. If the tissue contained multiple probes for Apbb1ip then these were combined together using a PCA to produce an eigengene (Ringnér, 2008), which was then correlated against the whole dataset for that tissue.

We found that the midbrain (GN381) had a very strong cis-eQTL (chromosome 2, 17.95–23.77 Mb, peak LOD = 11.75). As this clearly showed that Apbb1ip controls its own expression, we then looked for genes that coexpress with Apbb1ip. There were 3,093 probes which correlated with Apbb1ip expression, representing 1,607 unambiguously mapped genes (Supplementary Table S11). These genes showed significant enrichment for immune annotations, including immune response (GO:0006955, FDR = 3.25e-10; Supplementary Figure S5). This shows that the immune gene coexpression seen in the human brain is also seen in the murine brain.

Further, in the gastrointestinal data (Supplementary Table S12), there was significant enrichment for a range of immune-related annotations (e.g., immune system process, GO:0002376, FDR = 0e+00; Supplementary Figure S6). When correlated with the phenotypes database on GN, the most significant correlation was with visual discrimination learning using a touchscreen assay, mean reaction time to nose-poke in 3-month-old males (GN 16210, rho = −0.8137, p = 2.06e-05, N = 17). This is of interest, as reaction time is known to be disrupted in schizophrenia patients and their relatives (Chrobak et al., 2017).

In the hypothalamus data (GN159) we saw a significant cis-eQTL (19.02–36.10 Mb). There was significant coexpression with 17,960 probes, mapping to 9,160 genes, and once again we found immune system processes in the enrichment analysis (GO:0002376, FDR < 1e-20; Supplementary Figure S7). We also identified significant enrichment for focal adhesion-related annotations (mmu04510, FDR = 1.15e-08; WP85, FDR = 3.07e-07; Supplementary Table S13), previously shown to be altered in schizophrenia (Fan et al., 2013). The highest correlating phenotype was an infectious disease, immune system phenotype (rho = 0.7326, p = 4.71e-06, n = 27; GN 15561).

Pituitary (GN335) had nine probes with strong cis-eQTLs, so these were combined into an eigengene, and remapped to produce a cis-eQTL with a max LOD of 14.55 running from 19.02 to 23.77Mb. Once again, amongst the coexpressed genes (Supplementary Table S14) we saw enrichment for the immune system process (GO:0002376, FDR < 1e-20; Supplementary Figure S8). We also saw focal adhesion as significant again (mmu04510, FDR = 7.37e-07; Supplementary Table S15).

In the spleen (GN227), 3,700 probes, mappable to 2,846 genes, were significantly correlated with Apbb1ip (Supplementary Table S16). Unlike the majority of other tissues, and more similar to the T-cell data, the most enriched annotations were associated with cell cycle (GO:0007049, FDR = 1.5e-04; Supplementary Figure S9).

Given previous research suggesting striatal expression of genes affecting mental disorders (Ashbrook et al., 2015), and strong evidence for shared genetics between mental disorders (Cross-Disorder Group of the Psychiatric Genomics Consortium, 2013), we next analyzed striatal expression of Apbb1ip using a transcriptome dataset for postnatal day 14 mice. Again, we used GeneNetwork for QTL mapping and identified three probes for Apbb1ip, which show a strong QTL on chromosome 2, between 88 and 98 Mbps. We computed the 1st principal component of these three probe sets (an eigengene, Ringnér, 2008), and mapped this. This produced a peak LOD of 5.58 at 97.32 Mb, with a confidence interval from 87.46 to 99.75 Mbp. This region contains 146 genes, 88 of which contain SNPs or indels. We then narrowed down these genes to those which showed a significant cis-eQTL (indicating they regulate their own expression), and the expression of which correlated with Apbb1ip. These genes were Morf4l1, Mthfs, 1600029I14Rik, Plscr4, 1190002N15Rik, Armc8, Zbtb38, Faim, Rbp1, C730029A08Rik, Zic1, Trpc1, Adamts7, and 2810026P18Rik. Of these, Plscr4, Armc8, Zbtb38, Rbp1, Zic1, and Trpc1 have been linked to schizophrenia or schizophrenia-related disorders in previous studies (Kim et al., 2007; Li et al., 2017). For example, Zic genes are strongly expressed in the cerebellum and have been shown to mediate cerebellar development (Aruga and Millen, 2018), and behavioral abnormalities of Zic1 mutant mice can serve as models for diseases involving sensorimotor gating abnormalities, such as schizophrenia (Ogura et al., 2001).

We also discovered a strong eQTL in erythroid cells (GN150), chromosome 2:81.02–93.32, i.e., overlapping with the eQTL which we observed in the striatum. There are 5,830 probes (Supplementary Table S17), representing 4472 unambiguously mapped genes, and these genes were enriched for annotations such as RNA processing (GO:0006396, FDR < 1e-20) and translation (GO:0006412, FDR < 1e-20).

Enrichment of Schizophrenia-Associated Genes

Although we have identified genes that are coexpressed with Apbb1ip and are enriched for annotations linked to schizophrenia, this does not directly link them with schizophrenia. We thus performed a permutation analysis to look specifically for enrichment for genes that have been previously associated with schizophrenia.

We used the GWAS catalog for schizophrenia (GWAS-Catalog) to identify known associations. This had information for 2008 associations, however for some of these no genes were reported (n = 708) and for many, several potential genes were reported (i.e., genes within the LD block of the associated SNP). These associations represented 1326 unique human genes, of which 959 had murine homologs. For each of the above tissues, we produced a full list of genes examined in the entire dataset (e.g., all genes on a microarray or that were recorded in an RNA-seq experiment). For each tissue, we recorded the overlap in genes coexpressed with Apbb1ip and with genes associated with schizophrenia. We then performed 100,000 permutations, taking random subsets of the same length as the coexpressed genes, and recording their overlap with the schizophrenia genes. The number of times that the random genes had a higher overlap than the co-expressed genes, divided by 100,000, produced an empirical p-value.

Table 2: enrichment for schizophrenia-associated genes in Apbb1ip coexpression genes from mouse and human tissue. The species and tissue used for coexpression are shown, the significance of the overlap between genes associated with schizophrenia and genes coexpressed with Apbb1ip, the number of genes which overlap between genes associated with schizophrenia and genes coexpressed with Apbb1ip, the total number of genes which coexpressed with Apbb1ip and the total number of schizophrenia-associated genes included in the analysis (i.e., the number of schizophrenia-associated genes which were included in the expression dataset).

TABLE 2

Table 2. Permutation analysis of the overlap between genes coexpressing with Apbb1ip in mouse tissues, and genes associated with schizophrenia in genome-wide association studies (GWAS).

Discussion

We have identified APBB1IP, on chromosome 10 in humans and chromosome 2 in mice, as a novel gene that influences schizophrenia in humans and PPI in mice, and identified a number of genes which may be involved in this association. APBB1IP is found to have high expression in the whole blood, lymphocytes and the spleen, linking it to the immune system. Coexpression and enrichment analyses in the whole blood, cerebella and hippocampus showed enrichment for immune-related genes. Similar enrichment in immune-related processes showed that genes commonly coexpressing with Apbb1ip are also involved in immune responses in the mouse.

Our results suggest that APBB1IP is involved in cell motility, the immune response and neuronal development, all of which have been linked to schizophrenia. Apbb1ip is found to be expressed in homologous tissues to ABPP1IP, particularly the dentate gyrus, previously linked to schizophrenia (Yuan et al., 2015; Nakao et al., 2017). Immune gene coexpression is seen in the human brain and the murine brain. This confirms the translatability of this gene expression and supports mouse models as useful tools in which interventions can be tested. Several brain regions of potential interest, including the hippocampus and cerebellum, did not show any eQTL for Apbb1ip in the BXD population. This could have several explanations. First, it could be that Apbb1ip is under the regulation of many loci of small effect, which cannot individually be detected with the number of strains used. Another possibility is that the strains used for the particular gene expression study did not segregate for a variant which alters Apbb1ip expression. Finally, it could be that the timing was wrong and that Apbb1ip was not being differentially expressed at the time point expression was measured.

The potential mechanism for the link between ABPP1IP and an immune response is more speculative. Evidence suggests that T cells may play an important role in the pathophysiology of schizophrenia. Acute psychosis has been associated with the increased frequency of activated lymphocytes within the cerebrospinal fluid compared with controls (Nikkilä et al., 2001). ABPP1IP has been identified as a Rap1-binding protein important for leukocyte integrin activation (Lafuente et al., 2004). Rap1 has previously been identified as interacting with candidate schizophrenia susceptibility gene, Rapgef6, disrupting amygdala function in mice (Levy et al., 2015). Inflammation triggers leukocyte recruitment into site-specific tissue. Leukocyte integrins play a critical role in leukocyte recruitment and have thus emerged as key molecules in potential therapeutic targets for immune-mediated or inflammatory diseases (Yonekawa and Harlan, 2005). Regulation of lymphocyte adhesion to endothelium and trafficking of lymphocytes to secondary lymphoid organs is critical for adaptive immunity (Mor et al., 2007). Studies in ABPP1IP mice show that they are viable, fertile and apparently healthy (Stritt et al., 2015), but peripheral lymph nodes were depleted of both B and T cells resulting in abnormal lymphocyte trafficking and defective humoral immunity to T-cell-dependent antigens (Su et al., 2015). While systematic investigations of the humoral immune system in schizophrenia are lacking, there is some research suggesting that this “unspecific” arm of the immune system may be more activated in schizophrenic patients than in controls (Wilke et al., 1996; Müller et al., 1998). Activation of this system can be found mainly in unmedicated schizophrenic patients, which suggests it may reflect the disease process rather than the antipsychotic treatment.

Conclusion

In summary, the combination of human GWAS and mouse QTL from some of the largest study systems available has enabled us to identify a novel gene, APBB1IP, which influences schizophrenia in humans and PPI in mice.

Data Availability Statement

The datasets analyzed during the current study are available in the GeneNetwork repository, http://www.genenetwork.org/webqtl/main.py; the NHGRI-EBI Catalog of published Genome-wide association studies (GWAS), https://www.ebi.ac.uk/gwas/efotraits/EFO_0000692. The GTEx Project was supported by the Common Fund of the Office of the Director of the National Institutes of Health, and by NCI, NHGRI, NHLBI, NIDA, NIMH, and NINDS. The data used for the analyses described in this manuscript were obtained from www.gtexportal.org, the GTEx Portal on 04/06/2019.

Author Contributions

RH, SC and DA conceptualized and designed the study. SC and DA carried out analyses. DA and SC wrote the manuscript, with contributions from RH. All authors read and approved the submitted version.

Funding

The financial support for this study came from Economic and Social Research Council (ESRC; ES/P000347/1) and BBSRC funded studentship for SC.

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.

Footnotes

1. ^ http://grass.cgs.hku.hk/limx/kgg/
2. ^ http://gn2.genenetwork.org
3. ^ www.gtexportal.org

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fnbeh.2019.00266/full#supplementary-material.

FIGURE S1 | Enriched Directed Acyclic Graph (DAG) taken from Webgestalt (Zhang et al., 2005; Wang et al., 2013, 2017) showing enriched Gene Ontology (GO) Biological Process (BP) categories for genes coexpressing with APBB1IP in human whole blood, taken from GTEXv5. Categories shown in red are enriched in the input dataset, whereas uncolored categories are non-enriched parent-categories. Each box lists the name of the GO category, the number of genes in that category and the significance of enrichment (p-value).

FIGURE S2 | Enriched Directed Acyclic Graph (DAG) taken from Webgestalt (Zhang et al., 2005; Wang et al., 2013, 2017) showing enriched Gene Ontology (GO) Biological Process (BP) categories for genes coexpressing with APBB1IP in human cerebellum, taken from GTEXv5. Categories shown in red are enriched in the input dataset, whereas uncolored categories are non-enriched parent-categories. Each box lists the name of the GO category, the number of genes in that category and the significance of enrichment (p-value).

FIGURE S3 | Enriched Directed Acyclic Graph (DAG) taken from Webgestalt (Zhang et al., 2005; Wang et al., 2013, 2017) showing enriched Gene Ontology (GO) Biological Process (BP) categories for genes coexpressing with APBB1IP in GeneFriends. Categories shown in red are enriched in the input dataset, whereas uncolored categories are non-enriched parent-categories. Each box lists the name of the GO category, the number of genes in that category and the significance of enrichment (p-value).

FIGURE S4 | Enriched Directed Acyclic Graph (DAG) taken from Webgestalt (Zhang et al., 2005; Wang et al., 2013, 2017) showing enriched Gene Ontology (GO) Biological Process (BP) categories for genes coexpressing with Apbb1ip in GeneFriends. Categories shown in red are enriched in the input dataset, whereas uncolored categories are non-enriched parent-categories. Each box lists the name of the GO category, the number of genes in that category and the significance of enrichment (p-value).

FIGURE S5 | Enriched Directed Acyclic Graph (DAG) taken from Webgestalt (Zhang et al., 2005; Wang et al., 2013, 2017) showing enriched Gene Ontology (GO) Biological Process (BP) categories for genes coexpressing with the Apbb1ip in midbrain, from GeneNetwork. Categories shown in red are enriched in the input dataset, whereas uncolored categories are non-enriched parent-categories. Each box lists the name of the GO category, the number of genes in that category and the significance of enrichment (p-value).

FIGURE S6 | Enriched Directed Acyclic Graph (DAG) taken from Webgestalt (Zhang et al., 2005; Wang et al., 2013, 2017) showing enriched Gene Ontology (GO) Biological Process (BP) categories for genes coexpressing with the Apbb1ip eigengene in gastrointestinal tract, from GeneNetwork. Categories shown in red are enriched in the input dataset, whereas uncolored categories are non-enriched parent-categories. Each box lists the name of the GO category, the number of genes in that category and the significance of enrichment (p-value).

FIGURE S7 | Enriched Directed Acyclic Graph (DAG) taken from Webgestalt (Zhang et al., 2005; Wang et al., 2013, 2017) showing enriched Gene Ontology (GO) Biological Process (BP) categories for genes coexpressing with the Apbb1ip eigengene in hypothalamus, from GeneNetwork. Categories shown in red are enriched in the input dataset, whereas uncolored categories are non-enriched parent-categories. Each box lists the name of the GO category, the number of genes in that category and the significance of enrichment (p-value).

FIGURE S8 | Enriched Directed Acyclic Graph (DAG) taken from Webgestalt (Zhang et al., 2005; Wang et al., 2013, 2017) showing enriched Gene Ontology (GO) Biological Process (BP) categories for genes coexpressing with the Apbb1ip eigengene in pituitary gland, from GeneNetwork. Categories shown in red are enriched in the input dataset, whereas uncolored categories are non-enriched parent-categories. Each box lists the name of the GO category, the number of genes in that category and the significance of enrichment (p-value).

FIGURE S9 | Enriched Directed Acyclic Graph (DAG) taken from Webgestalt (Zhang et al., 2005; Wang et al., 2013, 2017) showing enriched Gene Ontology (GO) Biological Process (BP) categories for genes coexpressing with the Apbb1ip eigengene in spleen, from GeneNetwork. Categories shown in red are enriched in the input dataset, whereas uncolored categories are non-enriched parent-categories. Each box lists the name of the GO category, the number of genes in that category and the significance of enrichment (p-value).

TABLE S1 | Genes with strong, significant (r > 0.5 or < −0.5, p < 0.05), coexpressing with APBB1IP in human whole blood, taken from GTEXv5. The table shows the record ID (ENSEMBL gene ID), the gene ID (ENTREZ gene ID), the homologene ID, the gene symbol, the gene name, the location of the gene, the mean expression of the gene, the sample r (Pearson’s correlation), number of cases, the significance of the sample r, the literature correlation (similarity of words used to describe the gene and its function, using latent semantic indexing), and the tissue correlation (correlation between the two samples across tissue types). More details available at GeneNetwork.

TABLE S2 | Reactome pathways enriched in genes that are coexpressed with APBB1IP in human whole blood. Gene expression data taken from GTEXv5, enrichment analysis carried out using Webgestalt (Zhang et al., 2005; Wang et al., 2013, 2017). Reactome pathway IDs, the name of the pathway, the number of genes in the pathway, and the false discovery rate (FDR) are shown.

TABLE S3 | DrugBank drug targets enriched in genes that are coexpressed with APBB1IP in human whole blood. Gene expression data taken from GTEXv5, enrichment analysis carried out using Webgestalt (Zhang et al., 2005; Wang et al., 2013, 2017). DrugBank IDs, the name of the drug, the number of genes that are targets for the drug, and the false discovery rate (FDR) are shown.

TABLE S4 | Expression of APBB1IP in the developing human brain. Prenatal LMD Microarray data was taken from BrainSpan for 516 brain structures at postconception days 15, 16 and 21, using three different probes (Miller et al., 2014).

TABLE S5 | Genes with strong, significant (r > 0.5 or < −0.5, p < 0.05), coexpressing with APBB1IP in human cerebellum, taken from GTEXv5. The table shows the record ID (ENSEMBL gene ID), the gene ID (ENTREZ gene ID), the homologene ID, the gene symbol, the gene name, the location of the gene, the mean expression of the gene, the sample r (Pearson’s correlation), number of cases, the significance of the sample r, the literature correlation (similarity of words used to describe the gene and its function, using latent semantic indexing), and the tissue correlation (correlation between the two samples across tissue types). More details available at (GeneNetwork).

TABLE S6 | Wikipathways enriched in genes which are coexpressed with APBB1IP in the human cerebellum. Gene expression data taken from GTEXv5, enrichment analysis carried out using Webgestalt (Zhang et al., 2005; Wang et al., 2013, 2017). Wikipathway IDs, the name of the pathway, the number of genes in the pathway, and the false discovery rate (FDR) are shown.

TABLE S7 | Genes commonly coexpressed with APBB1IP in human as determined by GeneFriends. This table shows genes that commonly coexpress with APBB1IP, as identified by GeneFriends. This coexpression is independent of tissue or treatment. Entrez gene ID, gene symbol and coexpression value with APBB1IP are given for each of the genes.

TABLE S8 | Genes commonly coexpressed with Apbb1ip in the mouse as determined by GeneFriends. This table shows genes that commonly coexpress with Apbb1ip, as identified by GeneFriends. This coexpression is independent of tissue or treatment. Entrez gene ID, gene symbol and coexpression value with Apbb1ip are given for each of the genes.

TABLE S9 | The top 10,000 probes correlated with expression of Apbb1ip (probe ID 1443145_at) in T-helper cells of the BXD family (GN319; Sloan et al., 2016). The table shows, the probe ID, Mouse gene’s Entrez gene ID, the homologous human gene’s Entrez gene ID, the mouse gene symbol and location in the mouse genome.

TABLE S10 | The top 10,000 probes correlated with expression of Apbb1ip (probe ID 1443145_at) in T-regulatory cells of the BXD family (GN122; Sloan et al., 2016). The table shows, the probe ID, Mouse gene’s Entrez gene ID, the homologous human gene’s Entrez gene ID, the mouse gene symbol and location in the mouse genome.

TABLE S11 | The top 5,000 probes correlated with expression of Apbb1ip (probe ID A_55_P1961009) in the midbrain of the BXD family (GN381; Sloan et al., 2016). The table shows, the probe ID, Mouse gene’s Entrez gene ID, the homologous human gene’s Entrez gene ID, the mouse gene symbol and location in the mouse genome.

TABLE S12 | The top 15,000 probes correlated with the expression of Apbb1ip eigengene in the gastrointestinal tract of the BXD family. The table shows, the probe ID, Mouse gene’s Entrez gene ID, the homologous human gene’s Entrez gene ID, the mouse gene symbol and location in the mouse genome.

TABLE S13 | KEGG pathways enriched in genes that are coexpressed with Apbb1ip in mouse hypothalamus. Gene expression data are taken from GeneNetwork (GN159; Sloan et al., 2016), enrichment analysis carried out using Webgestalt (Zhang et al., 2005; Wang et al., 2013, 2017). KEGG pathway IDs, the name of the pathway, the number of genes in the pathway, and the false discovery rate (FDR) are shown.

TABLE S14 | The top 16,000 probes correlated with expression of Apbb1ip eigengene in the pituitary of the BXD family (GN335; Sloan et al., 2016). The table shows, the probe ID, Mouse gene’s Entrez gene ID, the homologous human gene’s Entrez gene ID, the mouse gene symbol and location in the mouse genome.

TABLE S15 | KEGG pathways enriched in genes that are coexpressed with Apbb1ip in mouse pituitary. Gene expression data are taken from GeneNetwork (GN335; Sloan et al., 2016). enrichment analysis carried out using Webgestalt (Zhang et al., 2005; Wang et al., 2013, 2017). KEGG pathway IDs, the name of the pathway, the number of genes in the pathway, and the false discovery rate (FDR) are shown.

TABLE S16 | The top 16,000 probes correlated with expression of Apbb1ip eigengene in the spleen of the BXD family (GN227; Sloan et al., 2016). The table shows, the probe ID, Mouse gene’s Entrez gene ID, the homologous human gene’s Entrez gene ID, the mouse gene symbol and location in the mouse genome.

TABLE S17 | The top 16,000 probes correlated with expression of Apbb1ip (probe ID ILM4560441) in erythroid cells of the BXD family (GN150; Sloan et al., 2016). The table shows, the probe ID, Mouse gene’s Entrez gene ID, the homologous human gene’s Entrez gene ID, the mouse gene symbol and location in the mouse genome.

References

Ackert-Bicknell, C. L., Karasik, D., Li, Q., Smith, R. V., Hsu, Y. H., Churchill, G. A., et al. (2010). Mouse BMD quantitative trait loci show improved concordance with human genome-wide association loci when recalculated on a new, common mouse genetic map. J. Bone Miner. Res. 25, 1808–1820. doi: 10.1002/jbmr.72

PubMed Abstract | CrossRef Full Text | Google Scholar

Aruga, J., and Millen, K. J. (2018). ZIC1 function in normal cerebellar development and human developmental pathology. Adv. Exp. Med. Biol. 1046, 249–268. doi: 10.1007/978-981-10-7311-3_13

PubMed Abstract | CrossRef Full Text | Google Scholar

Ashbrook, D. G., Arends, D., Prins, P., Mulligan, M. K., Roy, S., Williams, E. G., et al. (2019). The expanded BXD family of mice: a cohort for experimental systems genetics and precision medicine. BioRxiv [Preprint]. doi: 10.1101/672097

CrossRef Full Text | Google Scholar

Ashbrook, D. G., Williams, R. W., Lu, L., and Hager, R. (2015). A cross-species genetic analysis identifies candidate genes for mouse anxiety and human bipolar disorder. Front. Behav. Neurosci. 9:171. doi: 10.3389/fnbeh.2015.00171

PubMed Abstract | CrossRef Full Text | Google Scholar

Ashbrook, D. G., Williams, R. W., Lu, L., Stein, J. L., Hibar, D. P., Nichols, T. E., et al. (2014). Joint genetic analysis of hippocampal size in mouse and human identifies a novel gene linked to neurodegenerative disease. BMC Genomics 15:850. doi: 10.1186/1471-2164-15-850

PubMed Abstract | CrossRef Full Text | Google Scholar

Bagley, J. R. K. (2017). Localization of genetic variants that interact with prenatal stress to alter cocaine response and acoustic sensorimotor behaviors. BXD Phenotypes.

Barsh, G. S., Copenhaver, G. P., Gibson, G., and Williams, S. M. (2012). Guidelines for genome-wide association studies. PLoS Genet. 8:e1002812. doi: 10.1371/journal.pgen.1002812

PubMed Abstract | CrossRef Full Text | Google Scholar

Braff, D. L., Geyer, M. A., and Swerdlow, N. R. (2001). Human studies of prepulse inhibition of startle: normal subjects, patient groups and pharmacological studies. Psychopharmacology 156, 234–258. doi: 10.1007/s002130100810

PubMed Abstract | CrossRef Full Text | Google Scholar

Brain-Map. Available online at: www.brain-map.org.

BrainSpan. Available online at: http://www.brainspan.org/.

Brown, A. S., and Derkits, E. J. (2009). Prenatal infection and schizophrenia: a review of epidemiologic and translational studies. Am. J. Psychiatry 167, 261–280. doi: 10.1176/appi.ajp.2009.09030361

PubMed Abstract | CrossRef Full Text | Google Scholar

Bubier, J. A., and Chesler, E. J. (2012). Accelerating discovery for complex neurological and behavioral disorders through systems genetics and integrative genomics in the laboratory mouse. Neurotherapeutics 9, 338–348. doi: 10.1007/s13311-012-0111-3

PubMed Abstract | CrossRef Full Text | Google Scholar

Burmeister, M., McInnis, M. G., and Zöllner, S. (2008). Psychiatric genetics: progress amid controversy. Nat. Rev. Genet. 9, 527–540. doi: 10.1038/nrg2381

PubMed Abstract | CrossRef Full Text | Google Scholar

Cadenhead, K. S., Geyer, M. A., and Braff, D. L. (1993). Impaired startle prepulse inhibition and habituation in patients with schizotypal personality disorder. Am. J. Psychiatry 150, 1862–1867. doi: 10.1176/ajp.150.12.1862

PubMed Abstract | CrossRef Full Text | Google Scholar

Cadenhead, K. S., Swerdlow, N. R., Shafer, K. M., Diaz, M., and Braff, D. L. (2000). Modulation of the startle response and startle laterality in relatives of schizophrenic patients and in subjects with schizotypal personality disorder: evidence of inhibitory deficits. Am. J. Psychiatry 157, 1660–1668. doi: 10.1176/appi.ajp.157.10.1660

PubMed Abstract | CrossRef Full Text | Google Scholar

Carson, M. J., Doose, J. M., Melchior, B., Schmid, C. D., and Ploix, C. C. (2006). CNS immune privilege: hiding in plain sight. Immunol. Rev. 213, 48–65. doi: 10.1111/j.1600-065x.2006.00441.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Chrobak, A. A., Siuda-Krzywicka, K., Siwek, G. P., Tereszko, A., Janeczko, W., Starowicz-Filip, A., et al. (2017). Disrupted implicit motor sequence learning in schizophrenia and bipolar disorder revealed with ambidextrous serial reaction time task. Prog. Neuropsychopharmacol. Biol. Psychiatry 79, 169–175. doi: 10.1016/j.pnpbp.2017.06.025

PubMed Abstract | CrossRef Full Text | Google Scholar

Cross-Disorder Group of the Psychiatric Genomics Consortium. (2013). Identification of risk loci with shared effects on five major psychiatric disorders: a genome-wide analysis. Lancet 381, 1371–1379. doi: 10.1016/s0140-6736(12)62129-1

PubMed Abstract | CrossRef Full Text | Google Scholar

Craddock, N., O’Donovan, M., and Owen, M. (2005). The genetics of schizophrenia and bipolar disorder: dissecting psychosis. J. Med. Genet. 42, 193–204. doi: 10.1136/jmg.2005.030718

PubMed Abstract | CrossRef Full Text | Google Scholar

Dantzer, R., O’Connor, J. C., Freund, G. G., Johnson, R. W., and Kelley, K. W. (2008). From inflammation to sickness and depression: when the immune system subjugates the brain. Nat. Rev. Neurosci. 9, 46–56. doi: 10.1038/nrn2297

PubMed Abstract | CrossRef Full Text | Google Scholar

de Baumont, A., Maschietto, M., Lima, L., Carraro, D. M., Olivieri, E. H., Fiorini, A., et al. (2015). Innate immune response is differentially dysregulated between bipolar disease and schizophrenia. Schizophr. Res. 161, 215–221. doi: 10.1016/j.schres.2014.10.055

PubMed Abstract | CrossRef Full Text | Google Scholar

Fan, Y., Abrahamsen, G., Mills, R., Calderon, C. C., Tee, J. Y., Leyton, L., et al. (2013). Focal adhesion dynamics are altered in schizophrenia. Biol. Psychiatry 74, 418–426. doi: 10.1016/j.biopsych.2013.01.020

PubMed Abstract | CrossRef Full Text | Google Scholar

Genease. Available online at: http://research.cchmc.org/mershalab/GENEASE/analyze.html.

GeneNetwork. Available online at: http://www.genenetwork.org/reference.html.

Gokhale, A., Hartwig, C., Freeman, A. H., Das, R., Zlatic, S. A., Vistein, R., et al. (2016). The proteome of BLOC-1 genetic defects identifies the Arp2/3 actin polymerization complex to function downstream of the schizophrenia susceptibility factor dysbindin at the synapse. J. Neurosci. 36, 12393–12411. doi: 10.1523/jneurosci.1321-16.2016

PubMed Abstract | CrossRef Full Text | Google Scholar

Gottesman, I. I., and Gould, T. D. (2003). The endophenotype concept in psychiatry: etymology and strategic intentions. Am. J. Psychiatry 160, 636–645. doi: 10.1176/appi.ajp.160.4.636

PubMed Abstract | CrossRef Full Text | Google Scholar

Greenwood, T. A., Braff, D. L., Light, G. A., Cadenhead, K. S., Calkins, M. E., Dobie, D. J., et al. (2007). Initial heritability analyses of endophenotypic measures for schizophrenia: the consortium on the genetics of schizophrenia. Arch. Gen. Psychiatry 64, 1242–1250. doi: 10.1001/archpsyc.64.11.1242

PubMed Abstract | CrossRef Full Text | Google Scholar

GWAS-Catalog. Available online at: https://www.ebi.ac.uk/gwas/efotraits/EFO_0000692. Accessed December 10, 2018.

Hager, R., Lu, L., Rosen, G. D., and Williams, R. W. (2012). Genetic architecture supports mosaic brain evolution and independent brain-body size regulation. Nat. Commun. 3:1079. doi: 10.1038/ncomms2086

PubMed Abstract | CrossRef Full Text | Google Scholar

Hanisch, U.-K., and Kettenmann, H. (2007). Microglia: active sensor and versatile effector cells in the normal and pathologic brain. Nat. Neurosci. 10, 1387–1394. doi: 10.1038/nn1997

PubMed Abstract | CrossRef Full Text | Google Scholar

Hawrylycz, M. J., Lein, E. S., Guillozet-Bongaarts, A. L., Shen, E. H., Ng, L., Miller, J. A., et al. (2012). An anatomically comprehensive atlas of the adult human brain transcriptome. Nature 489, 391–399. doi: 10.1038/nature11405

PubMed Abstract | CrossRef Full Text | Google Scholar

Himes, B. E., Sheppard, K., Berndt, A., Leme, A. S., Myers, R. A., Gignoux, C. R., et al. (2013). Integration of mouse and human genome-wide association data identifies KCNIP4 as an asthma gene. PLoS One 8:e56179. doi: 10.1371/journal.pone.0056179

PubMed Abstract | CrossRef Full Text | Google Scholar

Jing, J., Pattaro, C., Hoppmann, A., Okada, Y., Fox, C. S., Köttgen, A., et al. (2016). Combination of mouse models and genomewide association studies highlights novel genes associated with human kidney function. Kidney Int. 90, 764–773. doi: 10.1016/j.kint.2016.04.004

PubMed Abstract | CrossRef Full Text | Google Scholar

Khandaker, G. M., Cousins, L., Deakin, J., Lennox, B. R., Yolken, R., and Jones, P. B. (2015). Inflammation and immunity in schizophrenia: implications for pathophysiology and treatment. Lancet Psychiatry 2, 258–270. doi: 10.1016/s2215-0366(14)00122-9

PubMed Abstract | CrossRef Full Text | Google Scholar

Khandaker, G. M., Pearson, R. M., Zammit, S., Lewis, G., and Jones, P. B. (2014). Association of serum interleukin 6 and C-reactive protein in childhood with depression and psychosis in young adult life: a population-based longitudinal study. JAMA Psychiatry 71, 1121–1128. doi: 10.1001/jamapsychiatry.2014.1332

PubMed Abstract | CrossRef Full Text | Google Scholar

Kim, S., Choi, K. H., Baykiz, A. F., and Gershenfeld, H. K. (2007). Suicide candidate genes associated with bipolar disorder and schizophrenia: an exploratory gene expression profiling analysis of post-mortem prefrontal cortex. BMC Genomics 8:413. doi: 10.1186/1471-2164-8-413

PubMed Abstract | CrossRef Full Text | Google Scholar

Kirkbride, J. B., Errazuriz, A., Croudace, T. J., Morgan, C., Jackson, D., Boydell, J., et al. (2012). Incidence of schizophrenia and other psychoses in England, 1950–2009: a systematic review and meta-analyses. PLoS One 7:e31660. doi: 10.1371/journal.pone.0031660

PubMed Abstract | CrossRef Full Text | Google Scholar

Kohl, S., Heekeren, K., Klosterkötter, J., and Kuhn, J. (2013). Prepulse inhibition in psychiatric disorders-apart from schizophrenia. J. Psychiatr. Res. 47, 445–452. doi: 10.1016/j.jpsychires.2012.11.018

PubMed Abstract | CrossRef Full Text | Google Scholar

Lafuente, E. M., van Puijenbroek, A. A., Krause, M., Carman, C. V., Freeman, G. J., Berezovskaya, A., et al. (2004). RIAM, an Ena/VASP and Profilin ligand, interacts with Rap1-GTP and mediates Rap1-induced adhesion. Dev. Cell 7, 585–595. doi: 10.1016/j.devcel.2004.07.021

PubMed Abstract | CrossRef Full Text | Google Scholar

Lein, E. S., Hawrylycz, M. J., Ao, N., Ayres, M., Bensinger, A., Bernard, A., et al. (2007). Genome-wide atlas of gene expression in the adult mouse brain. Nature 445, 168–176. doi: 10.1038/nature05453

PubMed Abstract | CrossRef Full Text | Google Scholar

Levy, R., Kvajo, M., Li, Y., Tsvetkov, E., Dong, W., Yoshikawa, Y., et al. (2015). Deletion of Rapgef6, a candidate schizophrenia susceptibility gene, disrupts amygdala function in mice. Transl. Psychiatry 5:e577. doi: 10.1038/tp.2015.75

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, J., Zhang, W., Yang, H., Howrigan, D. P., Wilkinson, B., Souaiaia, T., et al. (2017). Spatiotemporal profile of postsynaptic interactomes integrates components of complex brain disorders. Nat. Neurosci. 20, 1150–1161. doi: 10.1038/nn.4594

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, M., Jiang, L., and Li, J. (2014). KGG: A Systematic Biological Knowledge-Based Mining System for Genome-wide Genetic Studies (Version 4.0). Available online at: http://grass.cgs.hku.hk/limx/kgg/phasedgty.html. Accessed April 24, 2018.

Li, M.-X., Gui, H.-S., Kwan, J. S., and Sham, P. C. (2011). GATES: a rapid and powerful gene-based association test using extended Simes procedure. Am. J. Hum. Genet. 88, 283–293. doi: 10.1016/j.ajhg.2011.01.019

PubMed Abstract | CrossRef Full Text | Google Scholar

Light, G., Greenwood, T. A., Swerdlow, N. R., Calkins, M. E., Freedman, R., Green, M. F., et al. (2014). Comparison of the heritability of schizophrenia and endophenotypes in the COGS-1 family study. Schizophr. Bull. 40, 1404–1411. doi: 10.1093/schbul/sbu064

PubMed Abstract | CrossRef Full Text | Google Scholar

Lonsdale, J., Thomas, J., Salvatore, M., Phillips, R., Lo, E., Shad, S., et al. (2013). The genotype-tissue expression (GTEx) project. Nat. Genet. 45, 580–585. doi: 10.1038/ng.2653

PubMed Abstract | CrossRef Full Text | Google Scholar

Lopez, A. D., Mathers, C. D., Ezzati, M., Jamison, D. T., and Murray, C. J. (2006). Global and regional burden of disease and risk factors, 2001: systematic analysis of population health data. Lancet 367, 1747–1757. doi: 10.1016/s0140-6736(06)68770-9

PubMed Abstract | CrossRef Full Text | Google Scholar

Mignogna, K. M., Bacanu, S. A., Riley, B. P., Wolen, A. R., and Miles, M. F. (2018). Cross-species alcohol dependence-associated gene networks: co-analysis of mouse brain gene expression and human genome-wide association data. PLoS One 14:e0202063. doi: 10.1371/journal.pone.0202063

PubMed Abstract | CrossRef Full Text | Google Scholar

Miller, B. J., Buckley, P., Seabolt, W., Mellor, A., and Kirkpatrick, B. (2011). Meta-analysis of cytokine alterations in schizophrenia: clinical status and antipsychotic effects. Biol. Psychiatry 70, 663–671. doi: 10.1016/j.biopsych.2011.04.013

PubMed Abstract | CrossRef Full Text | Google Scholar

Miller, J. A., Ding, S.-L., Sunkin, S. M., Smith, K. A., Ng, L., Szafer, A., et al. (2014). Transcriptional landscape of the prenatal human brain. Nature 508, 199–206. doi: 10.1038/nature13185.

PubMed Abstract | CrossRef Full Text | Google Scholar

Mor, A., Dustin, M. L., and Philips, M. R. (2007). Small GTPases and LFA-1 reciprocally modulate adhesion and signaling. Immunol. Rev. 218, 114–125. doi: 10.1111/j.1600-065x.2007.00538.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Morse, H. C., Chused, T. M., Hartley, J. W., Mathieson, B. J., Sharrow, S. O., and Taylor, B. A. (1979). Expression of xenotropic murine leukemia viruses as cell-surface gp70 in genetic crosses between strains DBA/2 and C57BL/6. J. Exp. Med. 149, 1183–1196. doi: 10.1084/jem.149.5.1183

PubMed Abstract | CrossRef Full Text | Google Scholar

Mothersill, O., Knee-Zaska, C., and Donohoe, G. (2016). Emotion and theory of mind in schizophrenia-investigating the role of the cerebellum. Cerebellum 15, 357–368. doi: 10.1007/s12311-015-0696-2

PubMed Abstract | CrossRef Full Text | Google Scholar

Müller, N., Schlesinger, B., Hadjamu, M., Riedel, M., Schwarz, M., Primbs, J., et al. (1998). Increased frequency of CD8 positive gamma/delta T-lymphocytes (CD8+ γ/δ+) in unmedicated schizophrenic patients: relation to impairment of the blood-brain barrier and HLA-DPA*02011. Schizophr. Res. 12, 69–71. doi: 10.1016/s0920-9964(98)00036-x

PubMed Abstract | CrossRef Full Text | Google Scholar

Murray, K. N., Edye, M. E., Manca, M., Vernon, A. C., Oladipo, J. M., Fasolino, V., et al. (2019). Evolution of a maternal immune activation (mIA) model in rats: early developmental effects. Brain Behav. Immun. 75, 48–59. doi: 10.1016/j.bbi.2018.09.005

PubMed Abstract | CrossRef Full Text | Google Scholar

Nakao, A., Miyazaki, N., Ohira, K., Hagihara, H., Takagi, T., Usuda, N., et al. (2017). Immature morphological properties in subcellular-scale structures in the dentate gyrus of Schnurri-2 knockout mice: a model for schizophrenia and intellectual disability. Mol. Brain 10:60. doi: 10.1186/s13041-017-0339-2

PubMed Abstract | CrossRef Full Text | Google Scholar

Nikkilä, H. V., Müller, K., Ahokas, A., Rimón, R., and Andersson, L. C. (2001). Increased frequency of activated lymphocytes in the cerebrospinal fluid of patients with acute schizophrenia. Schizophr. Res. 49, 99–105. doi: 10.1016/s0920-9964(99)00218-2

PubMed Abstract | CrossRef Full Text | Google Scholar

Ogura, H., Aruga, J., and Mikoshiba, K. (2001). Behavioral abnormalities of Zic1 and Zic2 mutant mice: implications as models for human neurological disorders. Behav. Genet. 31, 317–324. doi: 10.1023/a:1012235510600

PubMed Abstract | CrossRef Full Text | Google Scholar

Owen, M., Sawa, A., and PB, M. (2016). Schizophrenia. Lancet 388, 86–97. doi: 10.1016/S0140-6736(15)01121-6

PubMed Abstract | CrossRef Full Text | Google Scholar

Patsoukis, N., Bardhan, K., Weaver, J. D., Sari, D., Torres-Gomez, A., Li, L., et al. (2017). The adaptor molecule RIAM integrates signaling events critical for integrin-mediated control of immune function and cancer progression. Sci. Signal. 10:eaam8298. doi: 10.1126/scisignal.aam8298

PubMed Abstract | CrossRef Full Text | Google Scholar

Peirce, J. L., Lu, L., Gu, J., Silver, L. M., and Williams, R. W. (2004). A new set of BXD recombinant inbred lines from advanced intercross populations in mice. BMC Genet. 5:7. doi: 10.1186/1471-2156-5-7

PubMed Abstract | CrossRef Full Text | Google Scholar

Peters, L. L., Robledo, R. F., Bult, C. J., Churchill, G. A., Paigen, B. J., and Svenson, K. L. (2007). The mouse as a model for human biology: a resource guide for complex trait analysis. Nat. Rev. Genet. 8:58. doi: 10.1038/nrg2025

PubMed Abstract | CrossRef Full Text | Google Scholar

Petryshen, T. L., Kirby, A., Hammer, R. P.Jr., Purcell, S., Singer, J. B., Hill, A. E., et al. (2005). Two quantitative trait loci for prepulse inhibition of startle identified on mouse chromosome 16 using chromosome substitution strains. Genetics 171, 1895–1904. doi: 10.1534/genetics.105.045658

PubMed Abstract | CrossRef Full Text | Google Scholar

PGC-GWAS. Available online at: http://www.med.unc.edu/pgc/results-and-downloads/downloads.

PheWeb. Available online at: http://pheweb.sph.umich.edu/.

Philip, V. M., Duvvuru, S., Gomero, B., Ansah, T. A., Blaha, C. D., Cook, M. N., et al. (2010). High-throughput behavioral phenotyping in the expanded panel of BXD recombinant inbred strains. Genes Brain Behav. 9, 129–159. doi: 10.1111/j.1601-183x.2009.00540.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Ringnér, M. (2008). What is principal component analysis? Nat. Biotechnol. 26:303. doi: 10.1038/nbt0308-303

PubMed Abstract | CrossRef Full Text | Google Scholar

Schizophrenia Working Group of the Psychiatric Genomics Consortium (2014). Biological insights from 108 schizophrenia-associated genetic loci. Nature 511, 421–427. doi: 10.1038/nature13595

PubMed Abstract | CrossRef Full Text | Google Scholar

Sloan, Z., Arends, D., Broman, K. W., Centeno, A., Furlotte, N., Nijveen, H., et al. (2016). GeneNetwork: framework for web-based genetics. J. Open Source Software 1:2. doi: 10.21105/joss.00025

CrossRef Full Text | Google Scholar

Stritt, S., Wolf, K., Lorenz, V., Vögtle, T., Gupta, S., Bösl, M. R., et al. (2015). Rap1-GTP-interacting adaptor molecule (RIAM) is dispensable for platelet integrin activation and function in mice. Blood 125, 219–222. doi: 10.1182/blood-2014-08-597542

PubMed Abstract | CrossRef Full Text | Google Scholar

Su, W., Wynne, J., Pinheiro, E. M., Strazza, M., Mor, A., Montenont, E., et al. (2015). Rap1 and its effector RIAM are required for lymphocyte trafficking. Blood 126, 2695–2703. doi: 10.1182/blood-2015-05-644104

PubMed Abstract | CrossRef Full Text | Google Scholar

Sullivan, P. F., Kendler, K. S., and Neale, M. C. (2003). Schizophrenia as a complex trait: evidence from a meta-analysis of twin studies. Arch. Gen. Psychiatry 60, 1187–1192. doi: 10.1001/archpsyc.60.12.1187

PubMed Abstract | CrossRef Full Text | Google Scholar

Sunkin, S. M., Ng, L., Lau, C., Dolbeare, T., Gilbert, T. L., Thompson, C. L., et al. (2013). Allen brain atlas: an integrated spatio-temporal portal for exploring the central nervous system. Nucleic Acids Res. 41, D996–D1008. doi: 10.1093/nar/gks1042

PubMed Abstract | CrossRef Full Text | Google Scholar

Swerdlow, N. R., Braff, D., and Geyer, M. (2000). Animal models of deficient sensorimotor gating: what we know, what we think we know and what we hope to know soon. Behav. Pharmacol. 11, 185–204. doi: 10.1097/00008877-200006000-00002

PubMed Abstract | CrossRef Full Text | Google Scholar

Swerdlow, N. R., Braff, D. L., and Geyer, M. A. (1999). Cross-species studies of sensorimotor gating of the startle reflex. Ann. N Y Acad. Sci. 877, 202–216. doi: 10.1111/j.1749-6632.1999.tb09269.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Tavitian, A., Song, W., and Schipper, H. M. (2019). Dentate gyrus immaturity in schizophrenia. Neuroscientist. doi: 10.1177/1073858418824072 [Epub ahead of print].

PubMed Abstract | CrossRef Full Text | Google Scholar

Tee, J. Y., Sutharsan, R., Fan, Y., and Mackay-Sim, A. (2017). Cell migration in schizophrenia: patient-derived cells do not regulate motility in response to extracellular matrix. Mol. Cell. Neurosci. 80, 111–122. doi: 10.1016/j.mcn.2017.03.005

PubMed Abstract | CrossRef Full Text | Google Scholar

Turetsky, B. I., Calkins, M. E., Light, G. A., Olincy, A., Radant, A. D., and Swerdlow, N. R. (2007). Neurophysiological endophenotypes of schizophrenia: the viability of selected candidate measures. Schizophr. Bull. 33, 69–94. doi: 10.1093/schbul/sbl060

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, J., Duncan, D., Shi, Z., and Zhang, B. (2013). WEB-based GEne SeT AnaLysis Toolkit (WebGestalt): update 2013. Nucleic Acids Res. 41, W77–W83. doi: 10.1093/nar/gkt439

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, J., Williams, R. W., and Manly, K. F. (2003). WebQTL: web-based complex trait analysis. Neuroinformatics 1, 299–308. doi: 10.1385/ni:1:4:299

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, J., Vasaikar, S., Shi, Z., Greer, M., and Zhang, B. (2017). WebGestalt 2017: a more comprehensive, powerful, flexible and interactive gene set enrichment analysis toolkit. Nucleic Acids Res. 45, W130–w137. doi: 10.1093/nar/gkx356

PubMed Abstract | CrossRef Full Text | Google Scholar

Wilke, I., Arolt, V., Rothermundt, M., Weitzsch, C., Hornberg, M., and Kirchner, H. (1996). Investigations of cytokine production in whole blood cultures of paranoid and residual schizophrenic patients. Eur. Arch. Psychiatry Clin. Neurosci. 246, 279–284. doi: 10.1007/bf02190280

PubMed Abstract | CrossRef Full Text | Google Scholar

Williams, R. W., Gu, J., Qi, S., and Lu, L. (2001). The genetic structure of recombinant inbred mice: high-resolution consensus maps for complex trait analysis. Genome Biol. 2:RESEARCH0046.0041. doi: 10.1186/gb-2001-2-11-research0046

PubMed Abstract | CrossRef Full Text | Google Scholar

Wu, Y., Williams, E. G., Dubuis, S., Mottis, A., Jovaisaite, V., Houten, S. M., et al. (2014). Multilayered genetic and omics dissection of mitochondrial activity in a mouse reference population. Cell 158, 1415–1430. doi: 10.1016/j.cell.2014.07.039

PubMed Abstract | CrossRef Full Text | Google Scholar

Yan, Z., Kim, E., Datta, D., Lewis, D. A., and Soderling, S. H. (2016). Synaptic actin dysregulation, a convergent mechanism of mental disorders? J. Neurosci. 36, 11411–11417. doi: 10.1523/jneurosci.2360-16.2016

PubMed Abstract | CrossRef Full Text | Google Scholar

Yonekawa, K., and Harlan, J. M. (2005). Targeting leukocyte integrins in human diseases. J. Leukoc. Biol. 77, 129–140. doi: 10.1189/jlb.0804460

PubMed Abstract | CrossRef Full Text | Google Scholar

Yuan, Y., Wang, H., Wei, Z., and Li, W. (2015). Impaired autophagy in hilar mossy cells of the dentate gyrus and its implication in schizophrenia. J. Genet. Genomics 42, 1–8. doi: 10.1016/j.jgg.2014.12.001

PubMed Abstract | CrossRef Full Text | Google Scholar

Zaldivar, A., and Krichmar, J. L. (2014). Allen brain atlas-driven visualizations: a web-based gene expression energy visualization tool. Front. Neuroinform. 8:51. doi: 10.3389/fninf.2014.00051

PubMed Abstract | CrossRef Full Text | Google Scholar

Zerbino, D. R., Achuthan, P., Akanni, W., Amode, M. R., Barrell, D., Bhai, J., et al. (2017). Ensembl 2018. Nucleic Acids Res. 46, D754–D761. doi: 10.1093/nar/gkx1098

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhang, B., Kirov, S., and Snoddy, J. (2005). WebGestalt: an integrated system for exploring gene sets in various biological contexts. Nucleic Acids Res. 33, W741–W748. doi: 10.1093/nar/gki475

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhang, Y., and Wang, H. (2012). Integrin signalling and function in immune cells. Immunology 135, 268–275. doi: 10.1111/j.1365-2567.2011.03549.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhou, X., and Stephens, M. (2012). Genome-wide efficient mixed-model analysis for association studies. Nat. Genet. 44, 821–824. doi: 10.1038/ng.2310

PubMed Abstract | CrossRef Full Text | Google Scholar