Exon expression in lymphoblastoid cell lines from subjects with schizophrenia before and after glucose deprivation
© Martin et al; licensee BioMed Central Ltd. 2009
Received: 18 December 2008
Accepted: 22 September 2009
Published: 22 September 2009
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.
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.
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).
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.
Microarray studies in postmortem brains of subjects with schizophrenia have implicated altered gene expression in pathways involved in myelination , GABAergic and glutamatergic transmission [2, 3], synaptic plasticity [4, 5] and metabolic pathways [6–8] [for review see ]. 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 [10–12], transformed lymphocytes [3, 6, 13], peripheral blood mononuclear cells [14–17], and peripheral blood cells  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 , schizophrenia [3, 6, 10, 18], bipolar disorder , Huntington disease , epilepsy , Tourette syndrome , 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 . 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)  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)  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 . In addition to psychosocial stressors, metabolic stressors may also influence dopamine and stress responses. Previous work has also shown an increased dopamine response  and HPA activation  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 . 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 . 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.
Lymphoblastoid cell lines (LCLs) from ten subjects were grown under normal conditions or under glucose deprivation.
The dorsolateral prefrontal cortex (DLPFC) and anterior cingulate (AnCg) were dissected and RNA extracted from a total of 60 subjects.
DLPFC 28S/18S (Mean ± S.D.)
AnCg 28S/18S (Mean ± S.D.)
DLPFC RIN (Mean ± S.D.)
DLPFC RIN (Mean ± S.D.)
2.26 ± 0.36
2.08 ± 0.38
7.87 ± 0.38
8.09 ± 0.65
2.33 ± 0.46
2.11 ± 0.36
7.95 ± 0.31
8.04 ± 0.80
2.21 ± 0.32
2.02 ± 0.31
7.76 ± 0.53
8.01 ± 0.57
2.00 ± 0.43
1.90 ± 0.39
7.63 ± 0.80
8.36 ± 0.55
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 .
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 [28–30] 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 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.
QPCR validation of top ranked genes in cohort 1
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.
Diagnosis × Glucose (p-value)
Mean ΔCt (C-GD)
Mean ΔCt (SZ-GD)
Ingenuity pathway analysis
Resequencing to identify variants associated with altered gene expression in schizophrenia
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.
Affymetrix Probeset ID
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 . 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.
IRF5 rs10954213 allele and genotype counts are shown for Cohort 2.
IRF5 rs10954213 genotype and allele counts are shown for Cohort 3.
Brain gene expression of candidate biomarkers
The QPCR results showed a significant decrease in brain gene expression of GLS in subjects with schizophrenia compared to controls.
BPD versus C p-value
MDD versus C p-value
SZ versus C p-value
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 . 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 .
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 . Subjects with schizophrenia also had higher fasting plasma levels of glucose, insulin and cortisol and were more insulin resistant than control subjects . 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 . Insulin has been shown to increase the expression of genes that are decreased in the hippocampus in subjects with schizophrenia in neuroblastoma cells . Moreover, glucose administration improves verbal declarative memory in clozapine-treated subjects with schizophrenia , although the mechanism of action remains unknown. Previous research has implicated metabolic and mitochondrial genes [3, 7, 8, 39, 40], including ornithine aminotransferase (OAT)  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  and in DLPFC  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  and PPARGC1B participates in mitochondrial biogenesis and plays a role in fat oxidation and nonoxidative glucose metabolism , 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 [46–48]. This gene is located at 20q12, a previously identified locus of schizophrenia . This gene is interesting in light of findings of increased prevalence of diabetes in schizophrenia , 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 , 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 . Of the 54 genes identified by Maycox et al (2009) , 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  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 . Previous work has shown a high level of agreement between digital counts by deep sequencing agreement with exon array . 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  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 , 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  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 . The SNP rs10954213 variant in IRF5 is in the 3'- untranslated region (UTR) previously shown to influence IRF5 gene expression  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 , 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 . 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 . 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 [58–60] 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.
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 , and genotyping in the identification of gene variants that are associated with regulation of expression in cell lines to find disease specific associations.
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.
- 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.View ArticlePubMedPubMed CentralGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- Ben-Shachar D, Laifenfeld D: Mitochondria, synaptic plasticity, and schizophrenia. Int Rev Neurobiol. 2004, 59: 273-296. full_text.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedPubMed CentralGoogle Scholar
- 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.PubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticleGoogle Scholar
- 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.View ArticlePubMedPubMed CentralGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedPubMed CentralGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- Walker EF, Diforio D: Schizophrenia: a neural diathesis-stress model. Psychol Rev. 1997, 104 (4): 667-685. 10.1037/0033-295X.104.4.667.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.PubMedGoogle Scholar
- 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.View ArticlePubMedPubMed CentralGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedPubMed CentralGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedPubMed CentralGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticleGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedPubMed CentralGoogle Scholar
- 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.View ArticlePubMedPubMed CentralGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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
- 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.View ArticlePubMedPubMed CentralGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedPubMed CentralGoogle Scholar
- Fraser HB, Xie X: Common polymorphic transcript variation in human disease. Genome Res. 2009, 19 (4): 567-575. 10.1101/gr.083477.108.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedPubMed CentralGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- The pre-publication history for this paper can be accessed here:http://www.biomedcentral.com/1755-8794/2/62/prepub