Skip to main content
  • Research article
  • Open access
  • Published:

Exon expression in lymphoblastoid cell lines from subjects with schizophrenia before and after glucose deprivation

Abstract

Background

The purpose of this study was to examine the effects of glucose reduction stress on lymphoblastic cell line (LCL) gene expression in subjects with schizophrenia compared to non-psychotic relatives.

Methods

LCLs were grown under two glucose conditions to measure the effects of glucose reduction stress on exon expression in subjects with schizophrenia compared to unaffected family member controls. A second aim of this project was to identify cis-regulated transcripts associated with diagnosis.

Results

There were a total of 122 transcripts with significant diagnosis by probeset interaction effects and 328 transcripts with glucose deprivation by probeset interaction probeset effects after corrections for multiple comparisons. There were 8 transcripts with expression significantly affected by the interaction between diagnosis and glucose deprivation and probeset after correction for multiple comparisons. The overall validation rate by qPCR of 13 diagnosis effect genes identified through microarray was 62%, and all genes tested by qPCR showed concordant up- or down-regulation by qPCR and microarray. We assessed brain gene expression of five genes found to be altered by diagnosis and glucose deprivation in LCLs and found a significant decrease in expression of one gene, glutaminase, in the dorsolateral prefrontal cortex (DLPFC). One SNP with previously identified regulation by a 3' UTR SNP was found to influence IRF5 expression in both brain and lymphocytes. The relationship between the 3' UTR rs10954213 genotype and IRF5 expression was significant in LCLs (p = 0.0001), DLPFC (p = 0.007), and anterior cingulate cortex (p = 0.002).

Conclusion

Experimental manipulation of cells lines from subjects with schizophrenia may be a useful approach to explore stress related gene expression alterations in schizophrenia and to identify SNP variants associated with gene expression.

Peer Review reports

Background

Microarray studies in postmortem brains of subjects with schizophrenia have implicated altered gene expression in pathways involved in myelination [1], GABAergic and glutamatergic transmission [2, 3], synaptic plasticity [4, 5] and metabolic pathways [68] [for review see [9]]. While examining gene expression in brain tissue from subjects with psychiatric or neurological disorders might be more relevant to disease pathology, assessment of whole blood or lymphocyte gene expression has several advantages, including practical availability of blood compared to brain tissue in human subjects and reduced confounds such as post-mortem interval (PMI), pH and agonal factors.

A number of studies have demonstrated the utility of using whole blood [1012], transformed lymphocytes [3, 6, 13], peripheral blood mononuclear cells [1417], and peripheral blood cells [18] as biomarkers in psychiatric and neurological disorders. These studies, while employing a variety of cell types, have been used to study panic disorder, nicotine dependence, post-traumatic stress disorder [15], schizophrenia [3, 6, 10, 18], bipolar disorder [10], Huntington disease [11], epilepsy [12], Tourette syndrome [12], Down syndrome and Alzheimer's disease [13, 16]. In addition to increased availability of lymphocytes and reduced confounds such as PMI, pH and agonal factors, transformed lymphocytes have the advantage of being more readily amenable to experimental manipulations than postmortem brain tissue. While it may be argued that state-dependent and brain tissue-dependent features are avoided by use of the cell lines, more stable genetic features may be prominent. Furthermore, previous work found a moderate correlation between blood and gene expression from 17 brain regions in humans, and reported that about half of a predetermined set of candidate genes relevant to schizophrenia were expressed in both whole blood and prefrontal cortex [19]. Experimental manipulation of cell lines from subjects with schizophrenia may be a useful approach to explore stress related gene expression alterations in schizophrenia and to identify SNP variants associated with gene expression. The purpose of this study was to examine the effects of glucose reduction stress on LCL gene expression in subjects with schizophrenia compared to non-psychotic relatives.

Recently, Naydenov et al (2007) [20] examined the effects of glucose deprivation in lymphocytes to study stress-induced changes in gene expression in bipolar disorder compared to controls. Naydenov et al (2007) [20] reported significant differences between subjects with bipolar disorder and controls in expression of genes related to the mitochondrial respiratory chain. This study demonstrated the advantages of combining a cellular model with experimental manipulation of glucose deprivation, which identified gene expression that was blunted in bipolar disorder but not controls.

The neural diathesis-stress model of schizophrenia proposes that in vulnerable individuals, psychosis can be triggered or worsened by stress and the interaction of stress with the dopamine system [21]. In addition to psychosocial stressors, metabolic stressors may also influence dopamine and stress responses. Previous work has also shown an increased dopamine response [22] and HPA activation [23] in subjects with schizophrenia compared to controls following metabolic stress due to 2-deoxyglucose treatment, although the molecular basis for this differential response is not understood. In the present study, it was hypothesized that a subset of genes would respond differentially to environmental manipulation (glucose reduction) in individuals with schizophrenia compared to controls.

Affymetrix 1.0 ST human exon arrays were utilized for this study to interrogate exon-level expression across the entire transcriptome. Exon arrays differ significantly from 3' expression arrays in the number and placement of the oligonucleotide probes and in the design of control probes for background correction. The Affymetrix 1.0 ST human exon array contains more than 5 million 25 bp oligomer probes, forming 1.4 million probe sets that together interrogate 1 million known and predicted exons up to four probes are selected from each putative exonic region. The coverage of the human 1.0 exon array represents over a 4-fold increase in probe density and an 8-fold increase in the number of perfect-match targets compared to the human U133 Plus 2.0 array. One advantage of the human 1.0 exon array platform is more robust measurements at the transcript level due to increased number of probes per transcript (30-40 probes across each entire RefSeq transcript compared to 11 probes mostly at the 3' transcript end in the U133 array). The probe sets are distributed along the entire length of the gene (as opposed to a concentration at the 3' end) allowing the potential to distinguish between different isoforms of a gene at the level of individual exons [24]. Therefore, this platform may identify alternative splicing events, which might play an important role in the pathophysiology of schizophrenia. Exon arrays have been demonstrated to attain improved sensitivity and specificity of absence/presence calls and may allow more accurate quantitative measurement of the level of gene expression compared to traditional 3' based arrays [25]. Quantitative real time PCR was used to validate exon array findings for 13 genes in lymphocyte samples. A subset of the validated genes was then assessed in brain samples from a separate cohort of subjects. Resequencing and TaqMan genotyping assays were conducted on one gene selected from exon array results in a larger cohort to identify gene variants associated with altered expression.

Methods

Subjects

Three cohorts were selected for this case-control study. Cohort 1 was used in the exon array analysis and consisted of LCLs from ten male subjects (five subjects with schizophrenia and five related non-psychotic subjects analyzed under normal and glucose deprivation conditions) for a total of 20 exon arrays (Table 1). The average age of controls and subjects with schizophrenia was 45 and 42, respectively. Cohort 2 was used for TaqMan genotyping assays on a total of 541 subjects (111 subjects with bipolar disorder, 19 subjects with major depressive disorder, 172 subjects with schizophrenia and 239 controls). Cohort 3 consisted of mRNA samples from two brain regions, dorsolateral prefrontal cortex (DLPFC) and anterior cingulate (AnCg), from 112 subjects, schizophrenia (n = 15), bipolar disorder (n = 17), major depressive disorder (n = 33) and controls (n = 47). Subjects with bipolar disorder and major depressive disorder were included in the analysis to test for the specificity of findings to schizophrenia. All three cohorts were obtained with IRB approval at the University of California, Irvine and cohorts 1 and 2 were obtained with IRB approval in Costa Rica. In cohort 3, 60 subjects were also genotyped (Table 2). Written, informed consent was obtained from all the individuals or from the surviving next-of-kin for all blood and brain samples reported in this study. Both cohorts 1 and 2 samples were collected from the Central Valley of Costa Rica in a relatively isolated population [26], and cohort 3 was collected in Southern California.

Table 1 Lymphoblastoid cell lines (LCLs) from ten subjects were grown under normal conditions or under glucose deprivation.
Table 2 The dorsolateral prefrontal cortex (DLPFC) and anterior cingulate (AnCg) were dissected and RNA extracted from a total of 60 subjects.

Ficoll extraction and Epstein-Barr virus transformation

Approximately 10 mL of blood was centrifuged at 1500 rpm for 10 min at RT (25°C). The upper layer was transferred to a 15 ml Falcon tube containing 5 ml Ficoll-Paque Plus (GE Healthcare, Piscataway, NJ) and centrifuged for 20 min at 2500 rpm at RT. The buffy coat layer was transferred to a new tube containing 10 ml of 1.0 M PBS (pH = 7.4) and centrifuged at 1000 rpm for 10 min at RT. The pellet was resuspended in 400 μl of supernatant and an equal volume was added to two wells of a culture plate. One ml of concentrated Epstein-Barr virus stock was added to each plate well and the plate was incubated for 45 min at 5% CO2 and 37°C. Next, 10 μl of phytohaemagglutinin (PHA) was added to each well and 400 μl of transforming media RPMI 1640 containing 15% fetal bovine serum, 2 mM glutamine and 25 mg of gentamicin. The cell lines were grown until confluent and then stored frozen in RPMI supplemented media and 10% dimethylsulphoxide (DMSO) until being grown for the experiment.

Glucose deprivation of lymphoblastoid cell lines

LCLs from cohort 1 subjects were grown under two environmental conditions (glucose deprived and normal glucose) to measure the effects of stress on exon expression and to examine state and trait-dependent alterations in exon expression in schizophrenia. All cell lines were thawed on the same date and grown with the same batch of RPMI-1640 supplemented media at 5% CO2 and 37°C to approximately the same cell density of 107 cells in a 25-cm2 flask. The cells from the 25-cm2 flask were split into two separate 75-cm2 flasks and grown to a confluency of approximately 5×107 cells/flask. One flask from each subject was left to grow in normal glucose (2%) while the other flask of cells was transferred into complete media containing 5.55 mM d-glucose for 24 hours. These cells were then placed into low glucose media containing 2.77 mM d-glucose (0.5%) for 48 hours. Microscopic cell death was not apparent when cell counts were taken using trypan blue stain.

Brain and lymphoblastoid cell line RNA extraction

Coronal slices of the brain were rapidly frozen on pre-cooled (to -120°C) aluminum plates, and stored at -80°C. The DLPFC and AnCg were later dissected and total RNA was extracted separately from each sample with TRIzol Reagent (Invitrogen, Carlsbad, CA) according to the supplier's protocol. Briefly, 1 ml of TRIzol was added to 5-10×106 cell pellets [10 ml of TRIzol per 50 ml confluent T75 flask of cells (pelleted)] or 50 mg DLPFC tissue. All final RNA pellets were resuspended in 100 μl of diethyl pyrocarbonate (DEPC) treated water. RNA degradation was assessed by quantifying the resulting 28S and 18S ribosomal band peak height ratios using the Agilent Bioanalyzer 2100 (Agilent, Palo Alto, CA) (see Tables 1 and 2).

Affymetrix Human Exon 1.0 ST Array

The Affymetrix Human Gene Chip Exon 1.0 ST Array interrogates over one million exons representing over 17,868 NCBI Reference Sequence (RefSeq) transcripts. Arrays were run using the manufacturer's technical protocol (Affymetrix, Santa Clara, CA). Briefly, 2 μg total RNA was subjected to a ribosomal RNA removal procedure (RiboMinus Human/Mouse Transcriptome Isolation Kit, Invitrogen) to reduce the 28S and 18S rRNA population in order to minimize background and increase sensitivity of the assay. Reduced RNA was then reverse-transcribed to cDNA using random hexamers tagged with a T7 promoter sequence followed by a second strand cDNA synthesis using DNA polymerase (GeneChip WT cDNA Synthesis Kit, Affymetrix). The resulting double stranded cDNA was then used for amplification of antisense cRNA and cleaned using the Gene Chip Sample Cleanup Module (Affymetrix). A second cycle cDNA synthesis was then performed using random primers to reverse transcribe the cRNA into sense single stranded DNA which was then fragmented, labeled, and hybridized to Affymetrix Human Gene Chip Exon 1.0 ST Arrays. Arrays were washed, stained, and scanned on the Affymetrix Fluidics Station and G7 Affymetrix high-resolution scanner using GCOS 1.3. Gene expression traits were derived from the cel files and analyzed with a Robust multi-array average (RMA) procedure [27].

Quantitative polymerase chain reaction (qPCR) of candidate diagnosis genes

The 100 μl reverse transcription reaction consisted of RNase-free water, TaqMan reverse transcriptase buffer (10 μl), 25 mM MgCl2 (22 μl), dNTPs (20 μl), Oligo d(T)16 (5 μl), RNase inhibitor (2 μl), MultiScribe Reverse Transcriptase, and 2 μg RNA from brain or lymphocytes. Real-time PCR was performed using the ABI Prism 7000 Sequence Detection System® (Applied Biosystems™), using the default thermal cycler program for all genes: 10 min of pre-incubation at 95°C followed by 40 cycles for 15 sec at 95°C and 1 min at 60°C. Individual real-time PCR reactions were carried out in a volume of 25 μl in 96-well plates (Applied Biosystems™) containing 12.5 μl SYBR Green, 5 μl cDNA (4 ng/μl), and 7.5 μl H2O. At the end of each reaction, the cycle threshold (Ct) was manually set at the level that reflected the best kinetic PCR parameters. Although each RNA sample was either DNase treated or Qiagen column purified, to further increase gene specificity one primer in each pair was designed (Primer Express, ABI) to span two exons and primers were not 3' UTR biased. The primers were tested for amplification of any residual genomic DNA contamination. The dissociation curves of real time PCR were monitored for primer-dimer pairings, which interfere with SybrGreen fluorescence measurements. A melting curve analysis was used for all primers and any primers that displayed amplification of genomic DNA less than 35 Ct were discarded and redesigned. (See Additional File 1 for primer sequences.)

TaqMan SNP genotyping

DNA samples were genotyped for an IRF5 variant (rs10954213) using TaqMan 5'-allele discrimination assays (Applied Biosystems, Foster City, CA). Allele-specific probes were labeled with the fluorescent dyes VIC and FAM, respectively. The PCR reaction was carried out in a total reaction of 5 μl (1 μl DNA, 2.5 μl TaqMan master mix, 0.125 μl 40× a IRF5 assay, 1.375 μl H2O) using the following amplification protocol: denaturation at 95° for 10 minutes, followed by 50 cycles at 92° for 15 seconds and at 58°C for 1.5 minutes. The genotype of each sample was determined by measuring allelic-specific fluorescence using the automated ABI Prism 7900 Sequence Detection system, with SDS 2.3 software for allelic discrimination (Applied Biosystems). Duplicate samples and negative controls were included to check the accuracy of genotyping.

Identification of SNPs associated with candidate gene expression

A second aim of this project was to identify cis-regulated transcripts associated with diagnosis. The mRNA by SNP Browser v 1.0.1 software was used to identify SNPs within the interferon regulatory factor 5 (IRF5) that were strongly related to expression and the regions around these SNPs were sequenced in subjects with exon array data. The mRNA by SNP Browser software [2830] contains association results of 54,675 transcripts with 406,912 SNPs (p < 0.05) and allows SNPs to be visualized in their genomic context along with linkage disequilibrium maps and putative haplotype blocks derived from the analysis of over 3 million SNPs genotyped in several populations by the International HapMap project. The IRF5 SNP with the most significant p-value related to expression in cohort 1 was followed up to examine whether there was an association between genotype and diagnosis in cohort 2, and genotype and expression in cohort 3.

Data analysis

Data were analyzed in Partek Genomic Solutions (St. Louis, MO) by a repeated measure ANOVA (mixed model) with Diagnosis × Probeset, Glucose deprivation × Probeset, Diagnosis × Probeset as main effects and corresponding interactions. Transcript expression was estimated by an average of the expression levels across all probesets in each RefSeq transcript. Probeset effects were examined to identify exon specific effects of diagnosis or glucose deprivation. The threshold for significance was obtained using a Benjamini-Hochberg step-down FDR cutoff of 25%.

Ingenuity pathway analysis

Genes with a step-down FDR less than or equal to 25% were used as the input variables for the data set to query the Ingenuity Pathways Analysis (IPA) software v7.1 canonical pathway analysis. Gene symbols were mapped to corresponding gene objects in the Ingenuity Pathways Knowledge Base. Each network or pathway was set to have a maximum of 35 focus genes and IPA identified those pathways that were most significant to the input data set. The significance of the association between the data set and the canonical pathway was determined based on the Benjamini-Hochberg step-down FDR calculated with the Fisher's exact test by calculating the probably that the association between the genes in the data set and the canonical pathway is due to chance alone.

Results

We tested the effects of glucose deprivation on LCL transcript profiles in subjects with schizophrenia compared to non-affected relatives. Diagnosis and glucose deprivation by probeset interaction effects were examined to identify exon-specific changes in expression. Transcripts were identified with significantly altered expression using an FDR cut-off of 25%. This cut-off was chosen based upon previous work by Benes et al, 2006 [31], which suggests an FDR cutoff of 5% may be too stringent for pathway analyses. There were no overall effects of diagnosis using a step-down FDR (Benjamini-Hochberg) cutoff of 25%. Step-down FDR is more conservative than step-up FDR and gives similar numbers of findings compared to Bonferroni correction [32]. Using a Benjamini-Hochberg step-down FDR cutoff of 25% there were 122 significant Diagnosis × Probeset effects, 328 Glucose deprivation × Probeset effects and 8 Diagnosis × Glucose deprivation × Probeset effects (see Additional file 2 for further detail; see Figures 1, 2 and 3 and Additional file 3 for exon array expression levels of representative genes GLS, DSC2 and ERO1L). Bonferroni correction of the mixed model ANOVA results was initially considered, but this approach assumed independent relationship among transcripts and these assumptions would be violated. Nevertheless, with a Bonferroni correction of 5%, no overall diagnosis effects, 72 Diagnosis × Probeset effects, 222 Glucose deprivation × Probesets effects and 5 Diagnosis × Glucose deprivation × Probeset effects were significant.

Figure 1
figure 1

(a) A significant Diagnosis × Glucose deprivation × Probeset interaction effect on GLS expression was detected by exon array. Exon array probeset expression levels are shown for the entire transcript. The probe set number is shown on the x-axis and the average expression for each diagnosis group is depicted on the y-axis. There was a significant decrease in probeset expression in subjects with schizophrenia (shown in blue) compared to controls (shown in red), under normal glucose conditions (b) The probeset with the most significant Diagnosis × Glucose interaction was probeset 2520346 (p = 0.001). Posthoc analysis revealed a significant decrease in GLS probeset expression in subjects with schizophrenia before and after glucose deprivation compared to normal controls (p = 0.001 and 0.0002, respectively).

Figure 2
figure 2

A significant Diagnosis × Glucose deprivation × Probeset interaction effect on DSC2 expression was detected by exon array. (a) Exon array probeset expression levels are shown for the entire transcript. The probe set number is shown on the x-axis and the average expression for each diagnosis group is depicted on the y-axis. The probe set number is shown on the x-axis, and the average expression for each diagnosis group is depicted on the y-axis. There was a significant Diagnosis × Probeset interaction effect on DSC2 expression, with a significant decrease in DSC2 expression in subjects with schizophrenia (shown in blue) compared to controls (shown in red). (b) The DSC2 probeset with the most significant diagnosis effect was probeset 3803005 (p = 2.06 × 10-5) and is depicted by dot plot. These findings were validated by qPCR (p = 0.0038).

Figure 3
figure 3

A significant Diagnosis × Glucose deprivation × Probeset interaction effect on ERO1L expression was detected by exon array. (a) Exon array probeset expression levels are shown for the entire transcript. The probe set number is shown on the x-axis and the average expression for each diagnosis group is depicted on the y-axis. There was a significant Diagnosis × Probeset interaction effect on ERO1L expression, with a significant decrease in ERO1L expression in subjects with schizophrenia (shown in blue) compared to controls (shown in red). (b) The probeset with the most significant diagnosis effect was probeset 3564812 (p = 0.0025) and is depicted by dot plot. These findings were validated by qPCR (p = 0.0037).

QPCR validation of top ranked genes in cohort 1

Thirteen candidate genes were chosen for qPCR validation based on p-value and fold-change of probesets with altered expression between schizophrenic and non-psychotic relatives' cell lines. Eight of 13 schizophrenia candidate genes were validated by qPCR (Table 3) for an overall validation rate by qPCR of 62% (p < 0.05, two tailed t-test) while the remaining 5 genes showed a trend for validation by qPCR (p < 0.1, two tailed t-test). All qPCR results were concordant compared to microarray in terms of the direction of fold changes. The exon array expression of three genes validated by qPCR (GLS, ERO1L and DSC2) is shown in Figures 1, 2 and 3).

Table 3 QPCR of LCLs was performed to confirm altered gene expression in twelve top schizophrenia candidate genes identified by exon array analysis from Diagnosis × Probeset interactions.

Ingenuity pathway analysis

Ingenuity pathway analysis was conducted using all transcripts with p-values below a Benjamini-Hochberg step-down FDR cutoff of 25%. For the 8 Diagnosis × Glucose deprivation × Probeset effect genes, IPA identified a network involving the transcription factor HFN4A, which was associated with disease specific responses to glucose deprivation. To increase the number of genes for an exploratory analysis, we examined a larger set of 67 Diagnosis × Glucose deprivation × Probeset effect genes identified by step-up FDR cutoff of 25% and a larger network involving transcription factor HFN4A was identified (see Figure 4). The three most significant functional categories were molecular transport, protein trafficking, and cellular function and maintenance (p < 0.05). The 122 genes with significant Diagnosis × Probeset effects based on a Benjamini-Hochberg step-down FDR cutoff of 25% were associated with no significant over-represented functional categories. A gene network of significant diagnosis effect genes is shown in Figure 5. The network includes 22 genes that showed diagnosis effects on gene expression. Calmodulin was a central complex located within this network. There were no over-represented functional categories for the 328 genes with significant Glucose deprivation × Probeset effects based on a Benjamini-Hochberg step-down FDR cutoff of 25%. There was not an over-representation of diagnosis-related genes implicated in schizophrenia.

Figure 4
figure 4

IPA analysis of genes with Diagnosis × Glucose deprivation × Probeset effects identified a network of 18 genes. Red indicates genes upregulated in schizophrenia in the glucose-deprived condition, green indicates genes down-regulated in schizophrenia in the glucose deprived condition, compared to normal controls in the glucose normal condition. White indicates the genes were not part of the dataset file genes. Solid lines indicate direct interactions, while dashed lines indicate indirect interactions. The network shapes are indicative of the molecular class (oval -- transcription regulator; square -- cytokine, trapezoid -- transporter; diamond -- enzyme; double circle -- complex/group; circle -- other). The network included 18 genes that showed diagnosis-dependent expression alterations following glucose deprivation, with Hepatocyte nuclear factor 4, alpha (HNF4A) as one node in this gene network.

Figure 5
figure 5

IPA analysis of genes with Diagnosis × Probeset effects identified a network of 22 genes. Red indicates genes upregulated in schizophrenia, green indicates genes down-regulated in schizophrenia. White indicates the genes were not part of the dataset file genes. Solid lines indicate direct interactions, while dashed lines indicate indirect interactions. See Figure 4 legend for description of molecular classes indicated by each network shape. The network included 22 genes that showed diagnosis effects on gene expression. Interestingly, calcium-activated calmodulin, which has previously been shown to be altered in expression in schizophrenia [49], was a centrally located complex within this network.

Resequencing to identify variants associated with altered gene expression in schizophrenia

A second aim of this study was to identify cis-regulated transcripts that were altered in schizophrenia to investigate gene variants that might be associated with genetic risk for schizophrenia. For all genes validated by qPCR, SNPs associated with expression were identified using the mRNA by SNP Browser [2830]. SNPs with highly significant LOD scores associated with expression of IRF5 were identified within the IRF5 gene using the mRNA by SNP Browser, but not for other Diagnosis × Probeset effect genes that were validated by qPCR (Table 4). Regions (~300 bps) around SNPs of interest [33] were resequenced in 10 subjects from cohort 1. Two SNPs of interest (rs4728142 and rs10488630) adjacent to exon array probeset 3023264 showed a significant association between genotypes (G and A alleles, respectively) and increased microarray expression of IRF5 (Table 5). These SNPs showed high linkage disequilibrium. Three additional SNPs not listed in the mRNA by SNP Browser software (rs11770589, rs10954214 and rs10954213) were identified within IRF5 probe set 3023264 by re-sequencing. These SNPs were included in our analysis because of their proximity to SNPs shown to be associated IRF5 expression. There was a significant effect of genotype of SNPs rs4728142, rs10954214 and rs10954213 on IRF5 expression in cohort 1 LCLs (Table 5).

Table 4 There was a significant effect of genotype at two SNPs on IRF5 expression according to the mRNA by SNP Browser data which supported the investigation of the effects of these SNPs on probeset expression in the present study.
Table 5 Candidate SNPs within probesets were genotyped by resequencing and tested for the effect of genotype on exon array expression obtained in LCLs.

SNP rs10954213 association with IRF5 expression in cohorts 2 and 3

TaqMan SNP genotyping was performed for IRF5 SNP rs10954213 on DNA extracted from lymphocyte from 541 individuals (cohort 2) and brain samples from 60 individuals (cohort 3). This SNP was selected because it was the most significantly associated with expression in our study and was also previously shown to be associated with expression [33]. Cohorts 2 and 3 included subjects with bipolar disorder and major depressive disorder in addition to subjects with schizophrenia and controls to investigate the specificity of an association of the genotype with schizophrenia.

The association of the rs10954213 SNP was also investigated in a larger sample of subjects (cohort 2, N = 541) with schizophrenia, bipolar disorder and major depressive disorder. There was no significant association between diagnosis and this IRF5 SNP (chi-square = 0.26, p = 0.99). (See tables 6 and 7 for genotype frequencies within each cohorts 2 and 3). The genotype frequencies for each diagnosis group were in Hardy-Weinberg equilibrium.

Table 6 IRF5 rs10954213 allele and genotype counts are shown for Cohort 2.
Table 7 IRF5 rs10954213 genotype and allele counts are shown for Cohort 3.

The expression of IRF5 was studied in cohort 3 in the DLPFC. The relationship between rs10954213 genotype and IRF5 expression was significant in these subjects in the DLPFC (F = 11.76; df = 2, 58; p = 8.11 × 10-5) and anterior cingulate cortex (F = 5.55; df = 2, 58; p = 0.007) from cohort 3 (see Figure 6). However, there were no significant associations between diagnosis and IRF5 expression in the DLPFC (F = 0.04; df = 1,59; p = 0.996) or the anterior cingulate (F = 0.28; df = 1,59; p = 0.84), nor Genotype × Diagnosis interaction effect on expression in the DLPFC (F = 0.90; df = 5, 55; p = 0.52) or the anterior cingulate (F = 0.90; df = 5, 55; p = 0.84).

Figure 6
figure 6

The relationship between IRF5 expression (measured by qPCR of the 3'-UTR of IRF5 in two regions of brain in cohort 3) and the IRF5 SNP rs10954213 genotype is shown. The ΔCt was calculated by subtracting from the Ct IRF5 the average Ct of two housekeeping genes (SNX16 and TSR1 in the anterior cingulate and SDHA and CAMK2A in the dorsolateral prefrontal cortex). A larger ΔCt indicates less expression of the target gene. The relationship between rs10954213 genotype and IRF5 probeset 3026432 expression was significant in these subjects in the DLPFC (F = 11.76; df = 1, 58;p = 8.11 × 10-5) and anterior cingulate cortex (F = 5.55; df = 1,58; p = 0.007) from cohort 2. The fold changes of expression in the anterior cingulate of subjects with the AG and GG genotypes compared to those with the AA genotypes were 0.52 and 0.50, respectively. The fold changes of expression in the anterior cingulate of subjects with the AG and GG genotypes compared to those with the AA genotypes were 0.33 and 0.38, respectively. Both regions showed a significant down-regulation in subjects with the AA genotype compared to the AG or GG genotype, consistent earlier reports of lower expression of the long 3'UTR in subjects carrying G alleles [33].

Brain gene expression of candidate biomarkers

Five additional genes with significantly altered probeset expression in LCLs between schizophrenia and unaffected relatives (ADCY1, DSC2, DSC3, GLS and PDE4D) were examined by qPCR expression in DLPFC tissue from subjects with schizophrenia and non-psychotic controls in cohort 3. The direction of change of GLS in both LCLs and the prefrontal cortex was consistent and decreased in schizophrenia in LCLs under glucose deprivation and the DLPFC. For comparison, subjects with bipolar disorder (BPD) and major depressive disorder (MDD) were included in the analysis. There was a significant decrease in expression of GLS in subjects with schizophrenia compared to controls (p = 0.04, Fold Change = 0.73), but no significant alterations of brain gene expression in schizophrenia of any other transcript examined. There were no significant alterations in brain expression of any of the other transcripts examined p-values for the in schizophrenia (see Table 8), although ADCY1 was down-regulated in major depressive disorder and bipolar disorder.

Table 8 The QPCR results showed a significant decrease in brain gene expression of GLS in subjects with schizophrenia compared to controls.

Discussion

The results of this study suggest exon arrays in conjunction with cell lines and brain samples may be a valuable tool to identify gene expression associated with genotype in the context of psychiatric disorders. Despite a relatively small sample of cell lines grown under two different conditions, this study identified eight differentially expressed transcripts between cases and controls when compared under two glucose conditions, 328 transcripts with significant Glucose deprivation × Probeset effects on expression, and a total of 122 transcripts with significant Diagnosis × Probeset effects on expression demonstrating probeset-specific effects of diagnosis and glucose deprivation on expression. These transcripts passed a 25% step-down Benjamini Hochberg false discovery and the complete dataset is available by request.

There were no significant main effects of diagnosis. This might be due to the fact analyses were performed using summarized transcript data and gene transcript expression was estimated using the mean expression of all probesets, which for some transcripts could be 50 - 100 probesets. On the other hand, the Glucose deprivation × Probeset and Diagnosis × Probeset interaction effects were based on probeset data summarized at the exon level which was 4 - 8 probesets. The finding that expression differences within closely matched schizophrenia and unaffected relatives are restricted within certain exons of a transcript is consistent with the alternatively spliced transcripts being an over-represented significant category between schizophrenia and controls. As the exon array does not contain junction probes, this conclusion is speculative, as the exon array would measure both spliced and unspliced transcripts, and does not actually differentiate spliced and unspliced transcripts except in situations where splicing is an all or none occurrence. Over 70% of human multi-exon genes are alternatively spliced [34]. Also, limited power may have played a role in our lack of diagnosis or glucose deprivation effects. The effect sizes detected for 100 most significant diagnosis effects ranged from 0.05 to 1.19. A power analysis done in G*Power showed 80% power to detect an effect size of 0.55 [35].

There is a growing body of evidence of glucose metabolism abnormalities in schizophrenia. In a study by Ryan et al (2003), more than 15% of subjects with schizophrenia had impaired fasting glucose tolerance [36]. Subjects with schizophrenia also had higher fasting plasma levels of glucose, insulin and cortisol and were more insulin resistant than control subjects [36]. There is also evidence for linkage between enzymes that control glycolysis and schizophrenia, including PFKFB2, a rate limiting enzyme in glycolysis, and hexokinase, a regulatory enzyme in the metabolic pathway [37]. Insulin has been shown to increase the expression of genes that are decreased in the hippocampus in subjects with schizophrenia in neuroblastoma cells [9]. Moreover, glucose administration improves verbal declarative memory in clozapine-treated subjects with schizophrenia [38], although the mechanism of action remains unknown. Previous research has implicated metabolic and mitochondrial genes [3, 7, 8, 39, 40], including ornithine aminotransferase (OAT) [7] and metabolic enzyme cytosolic malate dehydrogenase (MDH1) [3, 41] in schizophrenia, which have not been shown to be influenced by cigarette smoking or antipsychotic treatment. For these reasons, we examined the effects of glucose deprivation on gene expression in subjects with schizophrenia compared to controls.

Of the significantly altered transcripts identified in this study, only GLS showed significant Diagnosis × Probeset, significant Glucose deprivation × Probeset, and significant Diagnosis × Glucose effects. This gene and four other Diagnosis × Probeset effect genes were selected for further study in the DLPFC and only GLS was significant in DLPFC. Glutaminase, which converts glutamine to glutamate, was decreased in expression in glucose-deprived LCLs and decreased in expression in the DLPFC in subjects with schizophrenia compared to controls. This is in contrast to previous findings showing increases in glutaminase in the thalamus [42] and in DLPFC [43] of subjects with schizophrenia. In addition to GLS expression being regulated by cellular environment, the effects of diagnosis and glucose on expression varied by exon and may be isoform-specific. This gene may be worth further association study by SNP analysis or resequencing. In addition to glutaminase, there was a significant Diagnosis × Glucose deprivation × Probeset effect on PPARGC1B expression. PPARGC1B is located at 5q33, which has been previously linked to schizophrenia [44] and PPARGC1B participates in mitochondrial biogenesis and plays a role in fat oxidation and nonoxidative glucose metabolism [45], which is interesting in light of evidence of glucose metabolism abnormalities and mitochondrial dysfunction in schizophrenia. None of the other Diagnosis × Glucose deprivation × Probeset effect genes that passed a step-down FDR cutoff of 25% is known to play a role in glucose metabolism or to be sensitive to metabolic stressors. Thus, induction of gene expression changes in vitro by glucose reduction could be compared to increased glucose concentrations which may alter mitochondrial and metabolic function.

IRF5 expression was down-regulated in cohort 1. We examined the relationship between Interferon Regulatory Factor 5 (IRF5) expression and genotype and found expression of IRF5 to variation at SNP rs10954213 in all cohorts. However, IRF5 gene expression alteration in LCLs or brain was not observed in schizophrenia in a larger sample of subjects by RT-PCR and no SNP association with schizophrenia was observed.

IPA network analysis of transcripts with significant Diagnosis × Glucose deprivation × Probeset effects revealed a network of 18 genes that showed diagnosis-dependent expression alterations following glucose deprivation. D-glucose and hepatocyte nuclear factor 4, alpha (HNF4A) were two nodes in this gene network. HNF4A mutations have been associated with type II diabetes [4648]. This gene is located at 20q12, a previously identified locus of schizophrenia [44]. This gene is interesting in light of findings of increased prevalence of diabetes in schizophrenia [28], and further investigation of HNF4A variants in schizophrenia may be warranted. IPA network analysis of transcripts with significant Diagnosis × Probeset effects revealed a network of 18 down-regulated and four up-regulated genes. Interestingly, calmodulin, which has previously been shown to be altered in expression in schizophrenia [49], was a centrally located complex within this network. If replicated, genes in these networks may represent novel targets for further study in schizophrenia.

The results of this study were compared with a recent study of microarray gene expression in the DLPFC of subjects with schizophrenia and control subjects [50]. Of the 54 genes identified by Maycox et al (2009) [50], there were 12 genes with diagnosis or Diagnosis × Probeset effects displaying nominal p-values < 0.05 in the current study (see Additional file 1). These results indicate that case and control derived cell lines might be a reliable and easily maintained resource for the study of the expression phenotype that is disease associated.

The human transcriptome was recently shown to have over 90% alternatively spliced variants [51] and it will be increasingly important to identify those variants at the exon level that differentiate cases and controls. We have shown that differential exon usage is associated with case - control differences, while few transcripts showed a complete differentiation across the entire gene; instead, differential expression was restricted to specific exons within a transcript. To fully pursue individual splice variants requires more validation experiments, but having a profile across the transcript represents a first approximation.

The use of exon arrays allowed for the possibility of detecting exon-specific changes in expression [24]. Previous work has shown a high level of agreement between digital counts by deep sequencing agreement with exon array [52]. In identifying gene expression alterations in LCLs, it may be important to combine theoretically valid experimental manipulations of cell lines from subjects with schizophrenia to induce gene expression alterations relevant to the disorder. An example would be to silence gene expression through siRNA transfection, and measure effects on known interactors.

Expression of schizophrenia candidate genes in lymphoblastoid cell lines

One aspect of the present work was that many of the top candidate genes did not exhibit diagnosis or glucose deprivation effects. Although many schizophrenia candidate genes were expressed above background levels in LCLs, there was no overlap between the candidate genes identified in a recent meta-analysis of gene association studies [53] and genes with significant Diagnosis × Probeset × Glucose deprivation effects on expression found in this study. This may be due to tissue specificity (brain versus lymphocyte) of gene expression alterations in schizophrenia, or due to the selection of the control group, which consisted of relatives of subjects with schizophrenia and who likely carry many of the same genetic risk factors as their affected family members. Alternatively, reactive oxygen species generation might be increased more by increased glucose utilization [54], as we did not detect an over representation of mitochondrial or reactive oxygen genes by lowering glucose concentrations in cells. It is to be expected that increased glucose concentration would lead to more alterations in mitochondrial related genes.

A second aim of this study was to identify cis-regulated alterations in gene expression from the mRNA by SNP Browser [30] as one criterion for further study in the context of an association with the disease. There was a significant association between an IRF5 genotype and expression. A genotype effect of SNP rs10954213 on IRF5 expression was significant in LCLs, and was validated in the prefrontal cortex and anterior cingulate. However, there was no significant Diagnosis × Genotype effect on expression or a significant association of genotype with diagnosis in a larger sample for rs10954213. IRF5 exhibits genetically controlled polymorphic transcript variation [55]. The SNP rs10954213 variant in IRF5 is in the 3'- untranslated region (UTR) previously shown to influence IRF5 gene expression [33] in LCLs. This SNP creates a functional polyadenylation site and the A genotype correlates with increased expression of a transcript variant containing a shorter 3'-UTR [33], which was the same association found in LCLs, and brain tissue in the present study.

The 3'-UTR of mRNA has important functions in the stability, localization, and translation of mRNA [56]. MicroRNAs (miRNAs), are small noncoding RNAs that recruit an argonaute protein complex to a complementary target mRNA. Each miRNA-Ago complex interacts with a specific mRNA, typically through pairing of nucleotide bases between the miRNA sequence and complementary sequences in the mRNA's 3-untranslated region (3'UTR), which results in translation repression or degradation of the mRNA. A miRNA binding site for hsa-miR-185 was identified 1 bp downstream of SNP rs10954213 using the web based tool, MicroInspector [57]. It is possible that the effect of SNP rs10954213 variant on IRF5 protein levels may be mediated via this shortened 3'-UTR lacking this miRNA binding site. A further experiment to correlate reduction in protein expression with genotype would address this speculation.

We cannot rule out the possible effects of antipsychotics or cigarette smoking on gene expression. Recent research suggests that cigarette smoke and antipsychotic treatment may lead to gene-silencing effects through DNA methylation changes [5860] which could lead to long-term changes in cell lines derived from smokers. The overlap between genes included in the IPA analysis and those known to be altered by antipsychotics or nicotine were assessed with Ingenuity and no overlap was observed.

Conclusion

This study supports the use of experimental manipulation on LCLs to examine alterations in gene expression at the exon level in psychiatric disorders [3, 61]. These results support the combined use of exon array, mRNA by SNP Browser software [30], and genotyping in the identification of gene variants that are associated with regulation of expression in cell lines to find disease specific associations.

References

  1. Hakak Y, Walker JR, Li C, Wong WH, Davis KL, Buxbaum JD, Haroutunian V, Fienberg AA: Genome-wide expression analysis reveals dysregulation of myelination-related genes in chronic schizophrenia. Proc Natl Acad Sci USA. 2001, 98 (8): 4746-4751. 10.1073/pnas.081071198.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  2. Pongrac J, Middleton FA, Lewis DA, Levitt P, Mirnics K: Gene expression profiling with DNA microarrays: advancing our understanding of psychiatric disorders. Neurochem Res. 2002, 27 (10): 1049-1063. 10.1023/A:1020904821237.

    Article  CAS  PubMed  Google Scholar 

  3. Vawter MP, Ferran E, Galke B, Cooper K, Bunney WE, Byerley W: Microarray screening of lymphocyte gene expression differences in a multiplex schizophrenia pedigree. Schizophr Res. 2004, 67 (1): 41-52. 10.1016/S0920-9964(03)00151-8.

    Article  PubMed  Google Scholar 

  4. Ben-Shachar D, Laifenfeld D: Mitochondria, synaptic plasticity, and schizophrenia. Int Rev Neurobiol. 2004, 59: 273-296. full_text.

    Article  CAS  PubMed  Google Scholar 

  5. Mirnics K, Middleton FA, Lewis DA, Levitt P: Analysis of complex brain disorders with gene expression microarrays: schizophrenia as a disease of the synapse. Trends Neurosci. 2001, 24 (8): 479-486. 10.1016/S0166-2236(00)01862-2.

    Article  CAS  PubMed  Google Scholar 

  6. Vawter MP, Atz ME, Rollins BL, Cooper-Casey KM, Shao L, Byerley WF: Genome scans and gene expression microarrays converge to identify gene regulatory loci relevant in schizophrenia. Hum Genet. 2006, 119 (5): 558-570. 10.1007/s00439-006-0172-7.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  7. Middleton FA, Mirnics K, Pierri JN, Lewis DA, Levitt P: Gene expression profiling reveals alterations of specific metabolic pathways in schizophrenia. J Neurosci. 2002, 22 (7): 2718-2729.

    CAS  PubMed  Google Scholar 

  8. Prabakaran S, Swatton JE, Ryan MM, Huffaker SJ, Huang JT, Griffin JL, Wayland M, Freeman T, Dudbridge F, Lilley KS, et al: Mitochondrial dysfunction in schizophrenia: evidence for compromised brain metabolism and oxidative stress. Mol Psychiatry. 2004, 9 (7): 684-697. 10.1038/sj.mp.4001532.

    Article  CAS  PubMed  Google Scholar 

  9. Altar CA, Hunt RA, Jurata LW, Webster MJ, Derby E, Gallagher P, Lemire A, Brockman J, Laeng P: Insulin, IGF-1, and muscarinic agonists modulate schizophrenia-associated genes in human neuroblastoma cells. Biol Psychiatry. 2008, 64 (12): 1077-1087. 10.1016/j.biopsych.2008.08.031.

    Article  CAS  PubMed  Google Scholar 

  10. Tsuang MT, Nossova N, Yager T, Tsuang MM, Guo SC, Shyu KG, Glatt SJ, Liew CC: Assessing the validity of blood-based gene expression profiles for the classification of schizophrenia and bipolar disorder: a preliminary report. Am J Med Genet B Neuropsychiatr Genet. 2005, 133 (1): 1-5.

    Article  Google Scholar 

  11. Borovecki F, Lovrecic L, Zhou J, Jeong H, Then F, Rosas HD, Hersch SM, Hogarth P, Bouzou B, Jensen RV, et al: Genome-wide expression profiling of human blood reveals biomarkers for Huntington's disease. Proc Natl Acad Sci USA. 2005, 102 (31): 11023-11028. 10.1073/pnas.0504921102.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  12. Tang Y, Gilbert DL, Glauser TA, Hershey AD, Sharp FR: Blood gene expression profiling of neurologic diseases: a pilot microarray study. Arch Neurol. 2005, 62 (2): 210-215. 10.1001/archneur.62.2.210.

    Article  PubMed  Google Scholar 

  13. Kálmán J, Kitajka K, Pákáski M, Zvara Á, Juhász A, Vincze G, Janka Z, Puskás L: Gene expression profile analysis of lymphocytes from Alzheimer's patients. Psychiatric Genetics. 2005, 15: 1-6. 10.1097/00041444-200503000-00001.

    Article  PubMed  Google Scholar 

  14. Achiron A, Gurevich M: Peripheral blood gene expression signature mirrors central nervous system disease: The model of multiple sclerosis. Autoimmunity Reviews. 2006, 5: 517-522. 10.1016/j.autrev.2006.02.009.

    Article  CAS  PubMed  Google Scholar 

  15. Segman RH, Shefi N, Goltser-Dubner T, Friedman N, Kaminski N, Shalev AY: Peripheral blood mononuclear cell gene expression profiles identify emergent post-traumatic stress disorder among trauma survivors. Mol Psychiatry. 2005, 10 (5): 500-513. 10.1038/sj.mp.4001636. 425.

    Article  CAS  PubMed  Google Scholar 

  16. Maes OC, Xu S, Yu B, Chertkow HM, Wang E, Schipper HM: Transcriptional profiling of Alzheimer blood mononuclear cells by microarray. Neurobiol Aging. 2007, 28 (12): 1795-1809. 10.1016/j.neurobiolaging.2006.08.004.

    Article  CAS  PubMed  Google Scholar 

  17. Konradi C: Gene expression microarray studies in polygenic psychiatric disorders: applications and data analysis. Brain Res Brain Res Rev. 2005, 50 (1): 142-155. 10.1016/j.brainresrev.2005.05.004.

    Article  CAS  PubMed  Google Scholar 

  18. Glatt SJ, Everall IP, Kremen WS, Corbeil J, Sasik R, Khanlou N, Han M, Liew CC, Tsuang MT: Comparative gene expression analysis of blood and brain provides concurrent validation of SELENBP1 up-regulation in schizophrenia. Proc Natl Acad Sci USA. 2005, 102 (43): 15533-15538. 10.1073/pnas.0507666102.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  19. Sullivan PF, Fan C, Perou CM: Evaluating the comparability of gene expression in blood and brain. Am J Med Genet B Neuropsychiatr Genet. 2006, 141B (3): 261-268. 10.1002/ajmg.b.30272.

    Article  PubMed  Google Scholar 

  20. Naydenov A, MacDonald B, Ongur D, Konradi C: Differences in lymphocytes electron transport gene expression levels between subjects with bipolar disorder and normal controls in response to glucose deprivation stress. Archives of General Psychiatry. 2007, 64: 555-564. 10.1001/archpsyc.64.5.555.

    Article  PubMed  Google Scholar 

  21. Walker EF, Diforio D: Schizophrenia: a neural diathesis-stress model. Psychol Rev. 1997, 104 (4): 667-685. 10.1037/0033-295X.104.4.667.

    Article  CAS  PubMed  Google Scholar 

  22. Breier A, Davis OR, Buchanan RW, Moricle LA, Munson RC: Effects of metabolic perturbation on plasma homovanillic acid in schizophrenia. Relationship to prefrontal cortex volume. Arch Gen Psychiatry. 1993, 50 (7): 541-550.

    Article  CAS  PubMed  Google Scholar 

  23. Elman I, Adler CM, Malhotra AK, Bir C, Pickar D, Breier A: Effect of acute metabolic stress on pituitary-adrenal axis activation in patients with schizophrenia. Am J Psychiatry. 1998, 155 (7): 979-981.

    Article  CAS  PubMed  Google Scholar 

  24. Zhang X, Liu G, Lenburg ME, Spira A: Comparison of smoking-induced gene expression on Affymetrix Exon and 3'-based expression arrays. Genome Inform. 2007, 18: 247-257. full_text.

    PubMed  Google Scholar 

  25. Kapur K, Xing Y, Ouyang Z, Wong WH: Exon arrays provide accurate assessments of gene expression. Genome Biol. 2007, 8 (5): R82-10.1186/gb-2007-8-5-r82.

    Article  PubMed  PubMed Central  Google Scholar 

  26. DeLisi LE, Mesen A, Rodriguez C, Bertheau A, LaPrade B, Llach M, Riondet S, Razi K: Clinical characteristics of schizophrenia in multiply affected Spanish origin families from Costa Rica. Psychiatr Genet. 2001, 11 (3): 145-152. 10.1097/00041444-200109000-00006.

    Article  CAS  PubMed  Google Scholar 

  27. Irizarry RA, Hobbs B, Collin F, Beazer-Barclay YD, Antonellis KJ, Scherf U, Speed TP: Exploration, normalization, and summaries of high density oligonucleotide array probe level data. Biostatistics. 2003, 4 (2): 249-264. 10.1093/biostatistics/4.2.249.

    Article  PubMed  Google Scholar 

  28. Dixon AL, Liang L, Moffatt MF, Chen W, Heath S, Wong KC, Taylor J, Burnett E, Gut I, Farrall M, et al: A genome-wide association study of global gene expression. Nat Genet. 2007, 39 (10): 1202-1207. 10.1038/ng2109.

    Article  CAS  PubMed  Google Scholar 

  29. Moffatt MF, Kabesch M, Liang L, Dixon AL, Strachan D, Heath S, Depner M, von Berg A, Bufe A, Rietschel E, et al: Genetic variants regulating ORMDL3 expression contribute to the risk of childhood asthma. Nature. 2007, 448 (7152): 470-473. 10.1038/nature06014.

    Article  CAS  PubMed  Google Scholar 

  30. Cookson W, Liang L, Abecasis G, Moffatt M, Lathrop M: Mapping complex disease traits with global gene expression. Nature Reviews Genetics. 2009, 10: 184-194. 10.1038/nrg2537.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  31. Benes FM, Matzilevich D, Burke RE, Walsh J: The expression of proapoptosis genes is increased in bipolar disorder, but not in schizophrenia. Mol Psychiatry. 2006, 11 (3): 241-251. 10.1038/sj.mp.4001758.

    Article  CAS  PubMed  Google Scholar 

  32. Qian HR, Huang S: Comparison of false discovery rate methods in identifying genes with differential expression. Genomics. 2005, 86 (4): 495-503. 10.1016/j.ygeno.2005.06.007.

    Article  CAS  PubMed  Google Scholar 

  33. Graham RR, Kyogoku C, Sigurdsson S, Vlasova IA, Davies LR, Baechler EC, Plenge RM, Koeuth T, Ortmann WA, Hom G, et al: Three functional variants of IFN regulatory factor 5 (IRF5) define risk and protective haplotypes for human lupus. Proc Natl Acad Sci USA. 2007, 104 (16): 6758-6763. 10.1073/pnas.0701266104.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  34. Johnson JM, Castle J, Garrett-Engele P, Kan Z, Loerch PM, Armour CD, Santos R, Schadt EE, Stoughton R, Shoemaker DD: Genome-wide survey of human alternative pre-mRNA splicing with exon junction microarrays. Science. 2003, 302 (5653): 2141-2144. 10.1126/science.1090100.

    Article  CAS  PubMed  Google Scholar 

  35. Faul F, Erdfelder E, Lang AG, Buchner A: G*Power 3: a flexible statistical power analysis program for the social, behavioral, and biomedical sciences. Behav Res Methods. 2007, 39 (2): 175-191.

    Article  PubMed  Google Scholar 

  36. Ryan MC, Collins P, Thakore JH: Impaired fasting glucose tolerance in first-episode, drug-naive patients with schizophrenia. Am J Psychiatry. 2003, 160 (2): 284-289. 10.1176/appi.ajp.160.2.284.

    Article  PubMed  Google Scholar 

  37. Stone WS, Faraone SV, Su J, Tarbox SI, Van Eerdewegh P, Tsuang MT: Evidence for linkage between regulatory enzymes in glycolysis and schizophrenia in a multiplex sample. Am J Med Genet B Neuropsychiatr Genet. 2004, 127B (1): 5-10. 10.1002/ajmg.b.20132.

    Article  PubMed  Google Scholar 

  38. Stone WS, Seidman LJ, Wojcik JD, Green AI: Glucose effects on cognition in schizophrenia. Schizophr Res. 2003, 62 (1-2): 93-103. 10.1016/S0920-9964(02)00406-1.

    Article  PubMed  Google Scholar 

  39. Vawter MP, Tomita H, Meng F, Bolstad B, Li J, Evans S, Choudary P, Atz M, Shao L, Neal C, et al: Mitochondrial-related gene expression changes are sensitive to agonal-pH state: implications for brain disorders. Mol Psychiatry. 2006, 11 (7): 663-679. 10.1038/sj.mp.4001830.

    Article  CAS  Google Scholar 

  40. Iwamoto K, Bundo M, Kato T: Altered expression of mitochondria-related genes in postmortem brains of patients with bipolar disorder or schizophrenia, as revealed by large-scale DNA microarray analysis. Hum Mol Genet. 2005, 14 (2): 241-253. 10.1093/hmg/ddi022.

    Article  CAS  PubMed  Google Scholar 

  41. Vawter MP, Barrett T, Cheadle C, Sokolov BP, Wood WH, Donovan DM, Webster M, Freed WJ, Becker KG: Application of cDNA microarrays to examine gene expression differences in schizophrenia. Brain Res Bull. 2001, 55 (5): 641-650. 10.1016/S0361-9230(01)00522-6.

    Article  CAS  PubMed  Google Scholar 

  42. Bruneau EG, McCullumsmith RE, Haroutunian V, Davis KL, Meador-Woodruff JH: Increased expression of glutaminase and glutamine synthetase mRNA in the thalamus in schizophrenia. Schizophr Res. 2005, 75 (1): 27-34. 10.1016/j.schres.2004.12.012.

    Article  PubMed  Google Scholar 

  43. Gluck MR, Thomas RG, Davis KL, Haroutunian V: Implications for altered glutamate and GABA metabolism in the dorsolateral prefrontal cortex of aged schizophrenic patients. Am J Psychiatry. 2002, 159 (7): 1165-1173. 10.1176/appi.ajp.159.7.1165.

    Article  PubMed  Google Scholar 

  44. Gurling HM, Kalsi G, Brynjolfson J, Sigmundsson T, Sherrington R, Mankoo BS, Read T, Murphy P, Blaveri E, McQuillin A, et al: Genomewide genetic linkage analysis confirms the presence of susceptibility loci for schizophrenia, on chromosomes 1q32.2, 5q33.2, and 8p21-22 and provides support for linkage to schizophrenia, on chromosomes 11q23.3-24 and 20q12.1-11.23. Am J Hum Genet. 2001, 68 (3): 661-673. 10.1086/318788.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  45. Ling C, Poulsen P, Carlsson E, Ridderstrale M, Almgren P, Wojtaszewski J, Beck-Nielsen H, Groop L, Vaag A: Multiple environmental and genetic factors influence skeletal muscle PGC-1alpha and PGC-1beta gene expression in twins. J Clin Invest. 2004, 114 (10): 1518-1526.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  46. Hara K, Noda M, Waki H, Tobe K, Yamauchi T, Kadowaki H, Satou H, Tsukamoto K, Nagamatsu S, Yamagata K, et al: Maturity-onset diabetes of the young resulting from a novel mutation in the HNF-4alpha gene. Intern Med. 2002, 41 (10): 848-852. 10.2169/internalmedicine.41.848.

    Article  CAS  PubMed  Google Scholar 

  47. Weedon MN, Owen KR, Shields B, Hitman G, Walker M, McCarthy MI, Love-Gregory LD, Permutt MA, Hattersley AT, Frayling TM: Common variants of the hepatocyte nuclear factor-4alpha P2 promoter are associated with type 2 diabetes in the U.K. population. Diabetes. 2004, 53 (11): 3002-3006. 10.2337/diabetes.53.11.3002.

    Article  CAS  PubMed  Google Scholar 

  48. Gupta RK, Kaestner KH: HNF-4alpha: from MODY to late-onset type 2 diabetes. Trends Mol Med. 2004, 10 (11): 521-524. 10.1016/j.molmed.2004.09.004.

    Article  CAS  PubMed  Google Scholar 

  49. Broadbelt K, Jones LB: Evidence of altered calmodulin immunoreactivity in areas 9 and 32 of schizophrenic prefrontal cortex. J Psychiatr Res. 2008, 42 (8): 612-621. 10.1016/j.jpsychires.2007.07.006.

    Article  PubMed  Google Scholar 

  50. Maycox PR, Kelly F, Taylor A, Bates S, Reid J, Logendra R, Barnes MR, Larminie C, Jones N, Lennon M, et al: Analysis of gene expression in two large schizophrenia cohorts identifies multiple changes associated with nerve terminal function. Mol Psychiatry. 2009.

    Google Scholar 

  51. Wang W, Carvalho B, Miller ND, Pevsner J, Chakravarti A, Irizarry RA: Estimating genome-wide copy number using allele-specific mixture models. J Comput Biol. 2008, 15 (7): 857-866. 10.1089/cmb.2007.0148.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  52. Sultan M, Schulz MH, Richard H, Magen A, Klingenhoff A, Scherf M, Seifert M, Borodina T, Soldatov A, Parkhomchuk D, et al: A global view of gene activity and alternative splicing by deep sequencing of the human transcriptome. Science. 2008, 321 (5891): 956-960. 10.1126/science.1160342.

    Article  CAS  PubMed  Google Scholar 

  53. Allen NC, Bagade S, McQueen MB, Ioannidis JP, Kavvoura FK, Khoury MJ, Tanzi RE, Bertram L: Systematic meta-analyses and field synopsis of genetic association studies in schizophrenia: the SzGene database. Nat Genet. 2008, 40 (7): 827-834. 10.1038/ng.171.

    Article  CAS  PubMed  Google Scholar 

  54. Leloup C, Tourrel-Cuzin C, Magnan C, Karaca M, Castel J, Carneiro L, Colombani AL, Ktorza A, Casteilla L, Penicaud L: Mitochondrial reactive oxygen species are obligatory signals for glucose-induced insulin secretion. Diabetes. 2009, 58 (3): 673-681. 10.2337/db07-1056.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  55. Fraser HB, Xie X: Common polymorphic transcript variation in human disease. Genome Res. 2009, 19 (4): 567-575. 10.1101/gr.083477.108.

    Article  CAS  PubMed  Google Scholar 

  56. Vasudevan S, Tong Y, Steitz JA: Switching from repression to activation: microRNAs can up-regulate translation. Science. 2007, 318 (5858): 1931-1934. 10.1126/science.1149460.

    Article  CAS  PubMed  Google Scholar 

  57. Rusinov V, Baev V, Minkov IN, Tabler M: MicroInspector: a web tool for detection of miRNA binding sites in an RNA sequence. Nucleic Acids Res. 2005, W696-700. 10.1093/nar/gki364. 33 Web Server

  58. Toyooka S, Suzuki M, Tsuda T, Toyooka KO, Maruyama R, Tsukuda K, Fukuyama Y, Iizasa T, Fujisawa T, Shimizu N, et al: Dose effect of smoking on aberrant methylation in non-small cell lung cancers. Int J Cancer. 2004, 110 (3): 462-464. 10.1002/ijc.20125.

    Article  CAS  PubMed  Google Scholar 

  59. Dong E, Nelson M, Grayson DR, Costa E, Guidotti A: Clozapine and sulpiride but not haloperidol or olanzapine activate brain DNA demethylation. Proc Natl Acad Sci USA. 2008, 105 (36): 13614-13619. 10.1073/pnas.0805493105.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  60. Liu H, Zhou Y, Boggs SE, Belinsky SA, Liu J: Cigarette smoke induces demethylation of prometastatic oncogene synuclein-gamma in lung cancer cells by downregulation of DNMT3B. Oncogene. 2007, 26 (40): 5900-5910. 10.1038/sj.onc.1210400.

    Article  CAS  PubMed  Google Scholar 

  61. Sawa A, Wiegand GW, Cooper J, Margolis RL, Sharp AH, Lawler JF, Greenamyre JT, Snyder SH, Ross CA: Increased apoptosis of Huntington disease lymphoblasts associated with repeat length-dependent mitochondrial depolarization. Nat Med. 1999, 5 (10): 1194-1198. 10.1038/13518.

    Article  CAS  PubMed  Google Scholar 

Pre-publication history

Download references

Acknowledgements

We acknowledge the Pritzker Neuropsychiatric Disorders Research Consortium, NIMH Biomarker grant MH40730, NIMH Conte Center Grant P50 MH60398, and the William Lion Penzner Foundation for support of this work. The authors are members of the Pritzker Neuropsychiatric Disorders Research Consortium, which is supported by the Pritzker Neuropsychiatric Disorders Research Fund L.L.C. A shared intellectual property agreement exists between this philanthropic fund and the University of Michigan, Stanford University, the Weill Medical College of Cornell University, Hudson Alpha Institute of Biotechnology, the Universities of California at Davis, and at Irvine, to encourage the development of appropriate findings for research and clinical applications.

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Marquis P Vawter.

Additional information

Competing interests

These authors (AM, WB, RS, HA, EGJ, SJW, JB, LED, RMM, AS, WEB, MPV) have received funding from Pritzker Neuropsychiatric Disorders Research Fund LLC for this work. RMM is a consultant for Bay City Capital and received reimbursement for his consulting work which was not related with the topic of this paper. All authors indicated (AM, WB, RS, HA, EGJ, SJW, JB, LED, RMM, AS, WEB, MPB) are applying for a patent related to the content of the manuscript (Patent Application 20080199866). The authors have no other competing interests.

Authors' contributions

MVM performed computational analysis and generated qPCR, sequencing and TaqMan data. BR generated exon array data. PAS provided bioinformatics and statistical support. AM provided blood samples and diagnosis of living subjects. WB provided blood samples. RS provided diagnosis of postmortem subjects. EAM transformed and organized cell lines and assisted with genotyping. HA, EGJ, SJW and JB provided postmortem samples and critique of results. LED provided blood samples and diagnosis of living subjects. AS, RMM, WEB provided postmortem samples and critique of results. MPV contributed to the overall design of the project, provided critique of the results, participated in manuscript conceptual development and editing, and provided bioinformatics and statistical support. All authors read and approved the final manuscript.

Electronic supplementary material

Additional File 1: Primer sequences. The tables lists primer sequences for all primers used in this study. (DOC 24 KB)

12920_2008_126_MOESM2_ESM.DOC

Additional File 2: significant diagnosis, glucose deprivation and probeset interaction effects on transcript expression overlap with previous findings in the prefrontal cortex of subjects with schizophrenia compared to controls. Tables list transcripts with significant diagnosis × probeset, glucose deprivation × probeset and diagnosis × glucose deprivation × probeset effects, and overlap with previous findings in the prefrontal cortex of subjects with schizophrenia compared to controls. (DOC 686 KB)

12920_2008_126_MOESM3_ESM.DOC

Additional File 3: PPARGC1B and FBXW5 expression. Figures displaying significant Diagnosis × Glucose deprivation × Probeset interaction effects on transcripts expression for genes PPARGC1B and FBXW5. (DOC 279 KB)

Authors’ original submitted files for images

Rights and permissions

Open Access This article is published under license to BioMed Central Ltd. This is an Open Access article is distributed under the terms of the Creative Commons Attribution License ( https://creativecommons.org/licenses/by/2.0 ), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Reprints and permissions

About this article

Cite this article

Martin, M.V., Rollins, B., Sequeira, P.A. et al. Exon expression in lymphoblastoid cell lines from subjects with schizophrenia before and after glucose deprivation. BMC Med Genomics 2, 62 (2009). https://doi.org/10.1186/1755-8794-2-62

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/1755-8794-2-62

Keywords