Genomic Convergence Analysis of Schizophrenia: mRNA Sequencing Reveals Altered Synaptic Vesicular Transport in Post-Mortem Cerebellum

Abstract

Schizophrenia (SCZ) is a common, disabling mental illness with high heritability but complex, poorly understood genetic etiology. As the first phase of a genomic convergence analysis of SCZ, we generated 16.7 billion nucleotides of short read, shotgun sequences of cDNA from post-mortem cerebellar cortices of 14 patients and six, matched controls. A rigorous analysis pipeline was developed for analysis of digital gene expression studies. Sequences aligned to approximately 33,200 transcripts in each sample, with average coverage of 450 reads per gene. Following adjustments for confounding clinical, sample and experimental sources of variation, 215 genes differed significantly in expression between cases and controls. Golgi apparatus, vesicular transport, membrane association, Zinc binding and regulation of transcription were over-represented among differentially expressed genes. Twenty three genes with altered expression and involvement in presynaptic vesicular transport, Golgi function and GABAergic neurotransmission define a unifying molecular hypothesis for dysfunction in cerebellar cortex in SCZ.

About the Authors:

Joann Mudge

Affiliation: National Center for Genome Resources, Santa Fe, New Mexico, United States of America

Neil A. Miller

Affiliation: National Center for Genome Resources, Santa Fe, New Mexico, United States of America

Irina Khrebtukova

Affiliation: Illumina Inc., Hayward, California, United States of America

Ingrid E. Lindquist

Affiliation: National Center for Genome Resources, Santa Fe, New Mexico, United States of America

Gregory D. May

Affiliation: National Center for Genome Resources, Santa Fe, New Mexico, United States of America

Jim J. Huntley

Affiliation: National Center for Genome Resources, Santa Fe, New Mexico, United States of America

Shujun Luo

Affiliation: Illumina Inc., Hayward, California, United States of America

Lu Zhang

Affiliation: Illumina Inc., Hayward, California, United States of America

Jennifer C. van Velkinburgh

Affiliation: National Center for Genome Resources, Santa Fe, New Mexico, United States of America

Andrew D. Farmer

Affiliation: National Center for Genome Resources, Santa Fe, New Mexico, United States of America

Sharon Lewis

Affiliation: National Center for Genome Resources, Santa Fe, New Mexico, United States of America

William D. Beavis

Affiliation: National Center for Genome Resources, Santa Fe, New Mexico, United States of America

Faye D. Schilkey

Affiliation: National Center for Genome Resources, Santa Fe, New Mexico, United States of America

Selene M. Virk

Affiliation: National Center for Genome Resources, Santa Fe, New Mexico, United States of America

C. Forrest Black

Affiliation: National Center for Genome Resources, Santa Fe, New Mexico, United States of America

M. Kathy Myers

Affiliation: National Center for Genome Resources, Santa Fe, New Mexico, United States of America

Lar C. Mader

Affiliation: National Center for Genome Resources, Santa Fe, New Mexico, United States of America

Ray J. Langley

Affiliation: National Center for Genome Resources, Santa Fe, New Mexico, United States of America

John P. Utsey

Affiliation: National Center for Genome Resources, Santa Fe, New Mexico, United States of America

Ryan W. Kim

Affiliation: National Center for Genome Resources, Santa Fe, New Mexico, United States of America

Rosalinda C. Roberts

Affiliation: Department of Psychiatry, University of Alabama at Birmingham, Birmingham, Alabama, United States of America

Sat Kirpal Khalsa

Affiliation: Northern New Mexico College, Española, New Mexico, United States of America

Meredith Garcia

Affiliation: Northern New Mexico College, Española, New Mexico, United States of America

Victoria Ambriz-Griffith

Affiliation: Northern New Mexico College, Española, New Mexico, United States of America

Richard Harlan

Affiliation: Northern New Mexico College, Española, New Mexico, United States of America

Wendy Czika

Affiliation: SAS Institute, Cary, North Carolina, United States of America

Stanton Martin

Affiliation: SAS Institute, Cary, North Carolina, United States of America

Russell D. Wolfinger

Affiliation: SAS Institute, Cary, North Carolina, United States of America

Nora I. Perrone-Bizzozero

Affiliation: Department of Neurosciences, University of New Mexico, Albuquerque, New Mexico, United States of America

Gary P. Schroth

Affiliation: Illumina Inc., Hayward, California, United States of America

Stephen F. Kingsmore

* E-mail: [email protected]

Affiliation: National Center for Genome Resources, Santa Fe, New Mexico, United States of America

Introduction

Schizophrenia (SCZ [MIM #181500]) is a common, severe psychiatric disorder with a strong genetic component and complex inheritance [1]–[5]. Cytogenetic, linkage, positional cloning, candidate gene association and genome-wide studies have identified several credible SCZ risk genes. However, few have yet been replicated or translated into causal alleles, in vitro diagnostics or therapeutics [6]–[8]. In many studies of SCZ, genetic analysis has been impeded by phenotypic definition based upon multiple, subjectively ascertained, behavioral parameters that lack reference to biological mechanism [9]–[12]. In addition, more than 20 whole-genome linkage scans have demonstrated heterogeneity of linkage [13], suggesting the existence of genocopies (similar phenotypes that are determined by distinct risk loci) [12]. Evidence exists that the genetic architecture of SCZ may be further obscured by allelic heterogeneity (additive genetic variance segregating in the population at causative loci), epistasis (different combinations of loci producing a phenotype in different pedigrees), pleiotropy (loci that affect more than one phenotype) and phenocopies (heterogeneous environmental factors that mimic allelic effects) [14]. In addition, contributions of risk alleles to complex traits may not fit basal multiplicative and/or threshold models, and studies performed to date may not have had sufficient power, or appropriate theory, to assess non-linear (i.e., epistatic and genotype-by-environment) models. As a result, case-control association studies have identified numerous significantly associated susceptibility loci, but lack of replication among studies is widespread [3], [4], [8], [15]–[18]. Furthermore, the most validated loci were largely selected based on involvement in networks implicated in SCZ (such as dopaminergic and glutamatergic neurotransmission), introducing bias and limiting identification of novel risk factors.

Absence of a clear understanding of the molecular basis of SCZ imposes significant challenges to timely diagnosis and prognostic or therapeutic categorization [12], [19]. Supplementation of diagnostic criteria with biomarkers that are causally related to SCZ or endophenotypes may allow definition of homogeneous subgroups that are predictive of progression and therapeutic response in individual patients [20] and would serve as a starting point for development of therapeutics directed at causal variants. An alternative approach for molecular dissection of SCZ is identification of altered gene expression in affected tissues. Because gene expression reflects both genetic and environmental influences, it may be particularly useful for identifying risk factors for a complex disorder such as SCZ, which is believed to have a multifactorial etiology [21]. Two factors have hitherto limited the effectiveness of gene expression analysis in SCZ: Firstly, mRNA analyses in post-mortem brains in SCZ is challenging due to type I and type II errors resulting from variation in cause of death (affecting agonal gene expression), postmortem interval (affecting RNA quality), concurrent medication, substance abuse, age, sex, race and duration of illness [22]. Secondly, in common with genome-wide association studies, gene expression comparisons employing available cohort sizes sail between the Scylla of many false-positives due to multiple comparisons and the Charybdis of insufficient power to detect true-positives following statistical correction [23]. Recently, however, studies of mRNA expression in post mortem brains in SCZ that account for these variables have started to be reported [7], [21], [22].

An elegant, new approach to navigate Scylla and Charybdis and, thereby, accomplish molecular definition of SCZ may be genomic convergence analysis [24]. Predicated on an implication of the central dogma of molecular biology, genomic convergence analysis posits that clinically relevant nucleotide variation should result in detectable cis- and trans-effects in messenger RNA (mRNA) that amalgamate into functional changes in networks and pathways. Importantly, genomic convergence analysis provides a strategy to collectively interpret and employ the massive, disease-related data sets produced by unbiased (i.e. non-hypothesis driven) linkage and expression studies. Indeed, integration of gene expression and genetic linkage data has shown promise in several neurologic disorders [24]–[27] and has started to be applied to SCZ [7], [28]–[30].

mRNA sequencing with shotgun, massively parallel sequencing platforms has recently shown utility for measurement of transcript abundance, splice isoforms and allelic influence on gene expression [31]–[50]. mRNA abundance is determined by sequencing either 3′ end tags or random cDNA fragments (digital transcript expression, DTE), followed by read alignment to reference databases and calculation of aligned read frequencies. Potential advantages of DTE in comparison to array hybridization include: single molecule sensitivity (corresponding to approximately 1 mRNA molecule per 30 cells; Hayashizake, personal communication); absolute, rather than relative, measurement of transcript abundance; sequence verification for each measurement; comprehensive detection of both known and novel, unannotated transcripts and isoforms; applicability to any eukaryotic species; very little technical imprecision; absence of interference from abundant transcripts (e.g. globin); and extensibility to concomitant measurement of non-coding RNA and to detection of nucleotide and structural variation.

As the first stage of a genomic convergence analysis of SCZ, we describe shotgun mRNA sequencing of an affected tissue (post-mortem cerebellar cortex) of patients and controls, together with a rigorous and systematic approach to DTE. The approach presented is conservative due to application of statistical and bioinformatic methods that substantially reduce type I error rates by using statistical significance criteria rather than fold-change values and by incorporation of potential confounding variables, such as psychotropic medication and cause of death. It is also more comprehensive than prior microarray studies of SCZ. The stringency of this approach, in combination with genome-wide genotypes from the same individuals, is anticipated to enable a genomic convergence analysis that advances understanding of the molecular basis of SCZ.

Results

Sequencing-by-synthesis of cerebellar mRNA

16.7 billion nucleotides of shotgun, full length cDNA sequence data were generated using Illumina Genome Analyzer platforms with sequencing-by-synthesis (SBS) chemistry from 20 mRNA samples [39], [42], [43]. mRNA samples were isolated post-mortem from the lateral hemispheres of the cerebellar cortices of 14 patients with SCZ and 6 control individuals [51], [52] (Table 1). Unrelated subjects were chosen to facilitate sampling of genetic heterogeneity [53], [54]. Cases and controls were approximately matched for age (cases, 45.2±11.8 years; controls 41.3±9.2 years), sex (all male), race, post-mortem interval (cases, 12.2±5.0 hours; controls 17.7±3.3 hours), cause of death, autopsy brain pH (cases, 6.54±0.19; controls 6.46±0.10) and RNA integrity number (cases, 8.06±0.53; controls, 7.87±0.41) (Table 1). 12.5–38.7 million, high quality sequences of length 32–36 bp were generated per sample (Table 2). Sequences were aligned to the human genome and RefSeq transcript databases using the algorithm GMAP, which allowed <2 (≤6%) mismatches [55]. There was little intra-sample variability in the number of sequences aligned to each locus from run-to-run or instrument-to-instrument (Fig. 1A; all source coefficient of variation 3.4%). 43.5±6.7% of sequences aligned to a transcript and 69.4±9.6% to the genome, evidence that annotation of mRNA isoforms in Homo sapiens is incomplete [56] (Table 2). 91% of alignments were unique (Table 2). Reads aligning to more than one location contained repetitive, paralogous, polymorphic or low complexity sequences and primarily mapped to untranslated regions or highly polymorphic gene families, such as major histocompatibility genes [48]. Unmapped sequences did not align to mitochondrial or 1879 viral genomes, offering negative evidence of chronic viral etiology for SCZ in these patients.

[Figure omitted. See PDF.]

Figure 1. Characteristics of cDNA sequencing-by-synthesis with Illumina/Solexa instrument.

(A) Run-to-run comparisons of the number of reads aligned per reference transcript for 2 SBS runs of cerebellar cortex #41 (21.6 and 17.0 million reads, respectively). Runs used the same library but were performed on different sequencing instruments on different days. Coefficient of variation was 3.4%. (B) Histogram showing coverage along an “averaged” reference transcript for 1.2 Gb of cerebellar cortex #41 cDNA sequences. Coverage was calculated at 1% intervals along each transcript to which reads aligned. At each interval coverage was averaged across all transcripts and plotted. “Short transcripts” are all transcripts of ≤500 bp to which reads were aligned. “Long transcripts” are all transcripts ≥10 kb to which reads were aligned. Numbers in parentheses are the number of transcripts represented by each category. (C) Histogram showing greater than exponential decline in frequency with increasing SBS coverage on reference transcripts for cerebellar cortex #41 (1.2 Gb). Of 33,938 transcripts, 111 had more than 300× coverage. Maximum coverage was 2323×. Best fit power trendline is shown. (D) Histogram showing the number of reference transcripts covered as a function of the amount of sequence generated. Different levels of minimum average coverage were examined: at least one read aligned (black circles), 1× average coverage (red triangles), 2× average coverage (green “+”), 4× average coverage (dark blue “×”), and 8× average coverage (light blue squares). Data are from SBS of cerebellar cortex #41 (1.2 Gb).

https://doi.org/10.1371/journal.pone.0003625.g001

[Figure omitted. See PDF.]

Table 1. Clinical Features of Patients and mRNA Samples.

https://doi.org/10.1371/journal.pone.0003625.t001

[Figure omitted. See PDF.]

Table 2. cDNA Sequencing and Alignment Statistics.

https://doi.org/10.1371/journal.pone.0003625.t002

In order to further investigate differences in proportions of sequences aligning to different reference databases, an additional 2.2 billion nucleotides of shotgun, full length, paired cDNA sequence data were generated from mRNA sample 3S using an Illumina Genome Analyzer II platform and aligned to three reference databases (the human genome, RefSeq transcript and UniGene transcript databases): 46.1% of reads aligned to all three databases, representing known exons. 31.4% of reads aligned to the genome alone, representing novel exons, splice isoforms [57] or transcripts [56]. 5.5% mapped to RefSeq transcript alone, and 6.5% to UniGene transcript alone, representing sequences that span known exon junctions [57]. 10.6% of reads did not align to any of these reference sets, representing reads that span boundaries of novel exons or splice isoforms, novel genomic sequence, or poor quality sequence.

Coverage and Composition of cDNA Sequences

The number of transcripts detected in cerebellar cortex mRNA differed little between samples (33,200±1,000; Table 2), corresponding to 85±3% of RefSeq transcript entries. While 12.5 million sequences per sample was sufficient to reach a plateau in the number of transcripts detected, deeper sequence generation resulted in linear increase in average depth of coverage (Fig. 1D). As anticipated with hexamer-primed cDNA synthesis, the distribution of sequence alignments along transcripts appeared random. Transcript coverage showed a moderate 3′ bias, particularly among long transcripts, attributable to only very slight levels of degradation in these samples together with two rounds of poly-A+ RNA selection (Fig. 1B). Decreased coverage was observed at 5′ and 3′ termini due to edge effects of random priming very close (<30 nts) to the ends of the mRNA. No compositional bias was detected. The average depth of coverage achieved was 8.0-fold (Table 2). The distribution of transcript abundance in cerebellar cortex was interesting: 76–81% of expressed transcripts had an abundance of at least 1 read per million. However, as abundance increased, the number of transcripts declined greater than exponentially (Fig. 1C, r2 = 0.97). Thus, only 0.3–0.5% of expressed transcripts had an abundance of ≥1000 reads per million.

Comparison of Digital Transcript Expression and Array Hybridization

Gene expression was evaluated in mRNA samples from the 20 cerebellar cortices using both oligonucleotide array hybridization, the current standard, and tag frequencies from the Illumina cDNA sequencing assay (aligned reads per million). The dynamic range of read frequencies was two orders of magnitude greater than array hybridization (log10 dynamic range DTE–4.37; log10 dynamic range array hybridization–2.33; Fig. 2). While 41% of arrayed oligonucleotide probes had hybridization intensities greater than the conventional signal: noise threshold (Fig. 2), 85±3% transcripts had aligned reads. Array hybridization signals are transformed and often normalized prior to evaluation of inter-sample differences. For read frequencies, log transformation, but not normalization, improved overlaid kernel density estimates, univariate distribution results, and Mahalanobis distances (Fig. 3 and 4). Therefore, log transformed values were used in all subsequent expression analyses. Aligned read frequencies, adjusted for transcript length, correlated weakly with oligonucleotide array hybridization signals (r2 = 0.35; Fig. 2). Read frequencies exhibited much higher correlation coefficients in pair wise sample comparisons (r2 = 0.93–0.99) than array hybridization (r2 = 0.83–0.88; Fig. 5 and 6). Furthermore, pair wise sample correlations of genome- and transcript-aligned read frequencies had very similar correlation coefficients (data not shown). SCZ patients could readily be distinguished from controls by unsupervised principal component analysis (by Pearson product-moment correlation, Fig. 7B) or Ward hierarchical clustering of Pearson product-moment correlations of read frequencies (data not shown). In contrast, patients and controls could not be separated with these methods on the basis of array hybridization signals (Fig. 7A and data not shown). Log transformed genome- and transcript-aligned read frequencies showed identical hierarchical clustering and Pearson product-moment correlations of samples (data not shown).

[Figure omitted. See PDF.]

Figure 2. Comparison of gene expression levels measured by two methods.

Comparison of gene expression levels detected by average shotgun, mRNA SBS coverage per transcript (normalized for transcript length) versus Affymetrix oligonucleotide microarray hybridization for cerebellar cortex sample #41 (38.7 million reads). Microarray images were scaled to an average hybridization intensity of 200, and the threshold for expression was an average hybridization intensity of ≥50.

https://doi.org/10.1371/journal.pone.0003625.g002

[Figure omitted. See PDF.]

Figure 3. Comparison of Mahalanobis Distances of Gene Expression by Array Hybridization and Sequence Read Frequencies.

14 SCZ samples are indicated by blue circles, and 6 control samples by red circles. The Y-axis shows Mahalanobis distances of log transformed gene expression values. The dotted blue line indicates the cutoff value for outliers. Panel A: Log10 transformed Affymetrix array hybridization signals. Panel B: Log10 transformed genome-aligned read frequencies. Panel C: Log10 transformed transcript-aligned read frequencies. Log10 transformed array hybridization values (A) had a wider distribution of distances than Log10 transformed sequence read frequencies (B,C). Without log transformation, distances were greater and several samples represented outliers (data not shown).

https://doi.org/10.1371/journal.pone.0003625.g003

[Figure omitted. See PDF.]

Figure 4. Overlayed kernel density estimates of Gene Expression by Array Hybridization and Sequence Read Frequencies.

14 SCZ samples are indicated by blue lines, and 6 control samples by red lines. The X-axis shows log transformed gene expression values while the Y-axis shows kernel densities. Panel A: Log10 transformed Affymetrix array hybridization signals. Panel B: Log10 transformed genome-aligned read frequencies. Panel C: Log10 transformed transcript-aligned read frequencies. Log10 transformed array hybridization values less than 1.69 (equivalent to a calibrated hybridization signal of 50) are considered noise. Without log transformation, samples showed greater variability in kernel densities and sequence read frequencies showed a near exponential decay (Fig. 1C, data not shown).

https://doi.org/10.1371/journal.pone.0003625.g004

[Figure omitted. See PDF.]

Figure 5. Pairwise sample correlations of log10-transformed, genome-aligned read frequencies, showing pairwise correlation coefficients.

https://doi.org/10.1371/journal.pone.0003625.g005

[Figure omitted. See PDF.]

Figure 6. Pairwise sample correlations of log10-transformed, array hybridization signals, showing pairwise correlation coefficients.

https://doi.org/10.1371/journal.pone.0003625.g006

[Figure omitted. See PDF.]

Figure 7. Unsupervised principal component analysis of log10 transformed array hybridization signals (A) and log10 transformed read frequencies (B).

Three dimensional plots of principal component analysis by Pearson product-moment correlation. 14 SCZ samples are indicated by blue circles, and 6 control samples by red circles.

https://doi.org/10.1371/journal.pone.0003625.g007

Patient, sample, and experimental parameters were examined to quantify sources of variability in read-frequency-based DTE. Decomposition of principal components of variance showed that DTE attributed a greater proportion of the total variance to diagnosis (SCZ versus control, 45.3%) than array hybridization (14.1%) (Fig 8). The largest sources of variability of array hybridization results were mRNA quality metrics (Fig. 8A), whereas DTE results were only marginally influenced by mRNA quality. In contrast, large sources of variability in DTE results were year of sequence generation, sequencing instrument-to-instrument variation and cause of death (Fig. 8B). Variance component categories were identical in genome- and transcript-aligned read frequencies (data not shown).

[Figure omitted. See PDF.]

Figure 8. Principal components of variance of log10 transformed array hybridization signals (A) and log10 transformed read frequencies (B).

Variance components decomposition of principal components (with Pearson correlation), with partitioning of variability in terms of known effects. Patient (Diagnosis, Cause of death [Death_cause], Age, Race, Medication), sample (Post-mortem interval [PMI], brain pH [pH_2], RNA integrity number [RIN], RNA isolation date [RNA_isol_date]), and experimental (Year Sequenced, Average Read Quality, Average Read length, Cluster Station, Illumina 1G instrument [SNPster], Library Creator, % reads aligning, number of reads) parameters were examined to quantify sources of variability in read-frequency- and array hybridization-based gene expression.

https://doi.org/10.1371/journal.pone.0003625.g008

Digital Transcript Expression Results

Differences between SCZ patients and controls in gene expression in cerebellar cortex were identified with analysis of variance with diagnosis as the discriminatory effect and the major non-diagnosis components of variance as fixed effects (brain pH, post-mortem interval, RNA integrity number and RNA isolation date for array hybridization, and cause of death, post-mortem interval, sequencing instrument and year sequenced for read frequencies). Following FDR correction, no differences were identified between SCZ patients and controls by array hybridization, due to the magnitude of non-diagnosis components of variance. In contrast, 88 genes exhibited FDR-corrected differences in cerebellar expression between SCZ patients and controls in genome-aligned read frequencies and 152 genes differed significantly in transcript-aligned reads (Fig. 9; Table 3). Between these two sets of genes, 25 were identified in common to both genome- and transcript-aligned reads. 24 of these 25 genes showed congruent direction of change in genome- and transcript-aligned reads. 95% of genes exhibiting significant change in expression with only one alignment had congruent but non-significant change in the second alignment or were absent in one of the reference sets. The correlation coefficient of the case-control ratio of the 215 differentially expressed genes between genome- and transcript-aligned reads was 0.73 and between transcript-aligned reads and array hybridization was 0.46 (data not shown). A majority of genes with disparity in the magnitude of gene expression change between alignments had greater read counts for genome- than transcript-alignments, which, in some cases, appeared to be due to unannotated exons and novel splice isoforms (see below).

[Figure omitted. See PDF.]

Figure 9. Volcano plot of analysis of variance of log10-transformed, transcript- (A) and genome-aligned read frequencies (B), showing genes with FDR-corrected differences in LSMeans.

The x-axis shows the magnitude of the difference between 14 cases (S) and 6 controls (C), while the Y-axis shows the −log10(p value) for those differences. The red dashed line indicates the significance cutoff for differences with a control (−log10(p value) 3.17 for genome alignments and 2.94 for transcript alignments).

https://doi.org/10.1371/journal.pone.0003625.g009

[Figure omitted. See PDF.]

Table 3. Genes Exhibiting Significantly Altered Aligned Read Frequencies between SCZD Cases and Controls.

https://doi.org/10.1371/journal.pone.0003625.t003

Genes with the largest increase in expression in SCZ samples were the ISL2 transcription factor, LIM/homeodomain, >800% in genome and transcript alignments; AMME complex, gene 1, 580% in transcript alignments; 3-hydroxy-3-methylglutaryl-Coenzyme A synthase, 372% in transcript alignments; kelch domain containing 1, 360% in transcript alignments; interphotoreceptor matrix proteoglycan 2, 321% in transcript alignments; and NK2 transcription factor related, locus 3, 320% in genome alignments. Genes with the largest decrease in expression in SCZ samples were the solute carrier family 25 (ornithine transporter), member 2, 59% in transcript alignments and LOC126170, similar to peptidylprolyl isomerase A isoform 1, 73% in transcript alignments.

An SCZ candidate gene showing differential expression by both measures was the gamma-aminobutyric acid (GABA)-mediated neuroinhibitory receptor, GABRA1 (increased expression in SCZ by 52% [transcript alignment] and 46% [genome alignment], Table 3). Eleven other genes involved in GABA neurotransmission showed non-significant increases in expression: GABRA2 (59%), GABRA4 (72%), GABRR1 (69%), GABRB1 (168%), GABRG2 (42%), GABRE (35%), GABRB2 (33%), GABBR2 (27%), GABARAPL2 (35%), GABRA6 (14%) and GABRB3 (8%). Four GABAergic genes showed non-significant decreases in expression: GABBR1 (22%), GABRG3 (54%), GABRR2 (26%) and SLC6A1 (GAT-1; 16%).

GO annotation of genes with significantly altered expression revealed over-representation of membrane associated genes (Χ2 test, p<0.01), genes involved in Zinc binding or transport (p<0.02), regulation of transcription (p<e−6), Golgi apparatus (p<e−6) and vesicle-mediated transport (p<0.01). 13 of 20 genes involved in vesicular transport or Golgi apparatus were up-regulated.

Further annotation of these 20 genes revealed up-regulation of:

1. ▪ Nine genes encoding proteins involved in transport from the trans-Golgi network to the synaptic vesicle (golgi autoantigen, subfamily a, member 1, GOLGA1, 28% in transcript alignments, 14% in genome alignments; solute carrier family 35 (UDP-N-acetylglucosamine transporter), member A3, SLC35A3, 69% in transcript alignments, 18% in genome alignments; component of oligomeric golgi complex 6, COG6, 30% in transcript alignments, 37% in genome alignments; thyroid hormone receptor interactor 11, TRIP11, 48% in transcript alignments, 23% in genome alignments; adaptor-related protein complex 1, gamma 1 subunit, AP1G1, 6% in transcript alignments,−13% in genome alignments; ADP-ribosylation factor guanine nucleotide-exchange factor 2, ARFGEF2, 2% in transcript alignments,7% in genome alignments; vesicle docking protein p115, USO1, 36% in transcript alignments, 25% in genome alignments; Rho-associated, coiled-coil containing protein kinase 1, ROCK1, 85% in transcript alignments, 60% in genome alignments; RAB9B, 35% in transcript alignments, 23% in genome alignments and vacuolar protein sorting 35 homolog,

2. ▪ VPS35 (52% in transcript alignments, 31% in genome alignments), which is involved in retrograde transport to the Golgi apparatus,

3. ▪ Early endosome antigen 1, 162 kD, EEA1 (88% in transcript alignments, −1% in genome alignments), which is involved in homotypic fusion of early endosomes,

4. ▪ Synaptotagmin I, SYT1 (8% in transcript alignments, 15% in genome alignments), which is involved with synaptic vesicle exocytosis, and

5. ▪ AP2 associated kinase 1, AAK1 (64% in transcript alignments, 34% in genome alignments), which is involved in receptor-mediated endocytosis (Fig. 10, Table 3).

[Figure omitted. See PDF.]

Figure 10. Cartoon illustrating functions and/or synaptic locations of 23 proteins corresponding to genes with altered expression in SCZ.

15 genes were upregulated (green), whereas 8 were downregulated (red). Underlined genes had >30% change in expression. Two genes involved in transport from the endoplasmic reticulum to the Golgi (GOLM1 and GPSN2) were downregulated, ten involved in transport from the trans-Golgi network to the synaptic vesicle were upregulated (GOLGA1, SLC35A3, COG6, TRIP11, AP1G1, ARFGEF2, USO1, ROCK1, RAB9B and VPS35) and two were downregulated (STX10 and ARFRP1), two genes involved with synaptic vesicle exocytosis (EEA1 and SYT1) were upregulated and two were downregulated (SV2A and NCDN), one gene involved in receptor-mediated endocytosis was upregulated (AAK1) and one involved in retrograde transport back to the Golgi apparatus was downregulated (SNX17). In addition, three post-synaptic membrane genes showed altered expression: GABRA1 (upregulated), ZACN (downregulated) and CACNG2 (upregulated).

https://doi.org/10.1371/journal.pone.0003625.g010

Down-regulated transcripts corresponded to:

1. ▪ Two genes involved in transport from the trans-Golgi network to the synaptic vesicle (syntaxin 10, STX10, −17% in transcript alignments, −19% in genome alignments and ADP-ribosylation factor related protein 1, ARFRP1, −10% in transcript alignments, −12% in genome alignments),

2. ▪ Two genes involved in transport from the endoplasmic reticulum to the Golgi (golgi phosphoprotein 2, GOLM1, −14% in transcript alignments, −17% in genome alignments and synaptic glycoprotein 2, GPSN2, −8% in transcript alignments, −5% in genome alignments),

3. ▪ Two genes involved with synaptic vesicle exocytosis (synaptic vesicle glycoprotein 2A, SV2A, −23% in transcript alignments, −26% in genome alignments and neurochondrin, NCDN, −12% in transcript alignments, −17% in genome alignments), and

4. ▪ Sorting nexin 17, SNX17 (−13% in transcript alignments, −21% in genome alignments) which is involved in retrograde transport back to the Golgi apparatus (Fig. 10, Table 3).

In addition, three post-synaptic membrane genes showed altered expression: GABRA1 (see above), ZACN (ligand-gated ion channel, zinc activated 1; down-regulated 21% in transcript alignments and 5% in genome alignments) and CACNG2 (calcium channel, voltage-dependent, gamma subunit 2; up-regulated 169% in transcript alignments and 10% in genome alignments) (Fig 10, Table 3).

Identification of Novel Splice Isoforms

10,022 genes exhibited ≥20% more reads aligned to genomic loci than to the corresponding, annotated transcripts, evidence for novel, unannotated exons. Putative use of alternative 5′ and 3′ terminal exons, novel internal exons and read-through of introns were all detected [57]. A SCZ susceptibility locus that exhibited intron-read-through in cerebellar mRNA samples was proline oxidase (PRODH, EC 1.5.99.8, SCZ4): 421 cerebellar cortex 41 reads aligned to the sole annotated PRODH transcript, while 838 aligned to the genomic locus. 33 of the latter mapped to intron 14 (nucleotides 17280876–17280950 on Chromosome 22), resulting in insertion of 25 new, in-frame amino acids adjacent to the N-terminus (Fig. 11). In addition to reads mapping within this intron, reads showing splicing out of the intron were also observed (Fig. 11). A further 112 reads aligned to 660 nucleotides of intron 13, generating an alternative N-terminus and premature stop codon (Fig. 11). Presence of these intronic PRODH sequences was observed in all cerebellum mRNA samples and is supported by independently isolated Sanger EST sequences (Genbank accession numbers CN362766, AA300535, AA322439, CD671570, CN362770).

[Figure omitted. See PDF.]

Figure 11. Novel alternative splicing of the PRODH locus.

A. Exons (green boxes) and introns (yellow lines) of the PRODH locus are shown. A 100 bp shaded box covers part of intron 13, intron14 and exon 14 and part of exon 15. Two sSNPs are illustrated by vertical blue hash marks. B. Alignments of SBS reads to the introns and exons within the 1000 bp shaded box are shown. Sequence reads are shown as arrows pointing in the direction of their orientation relative to the genomic reference. Yellow arrows represent sequence reads that map to more than one region on the genomic reference, and green arrows represent uniquely mapping reads. Three cDNA reads that omit intron 14 and contain exon 14 and 15 sequences are highlighted. Also highlighted are cDNA reads that map within intron 14. C. The sequence of the PRODH genomic region shown in B. Green highlights represent exonic nucleotides with aligned sequence reads and yellow represent intronic nucleotides with aligned sequence reads. Reads aligning uniquely to intron 14 (75 nts) indicate the existence of an alternative splice isoform that reads through this intron. Similarly, the 3′ region of intron 13 appears to be included in a novel splice isoform(s).

https://doi.org/10.1371/journal.pone.0003625.g011

As noted above, the presence of unannotated exons and novel splice isoforms appeared to explain disparity in the magnitude of gene expression change between cases and controls when using genome- and transcript-alignments. Dedicator of cytokinesis 3 (DOCK3), for example, was significantly up-regulated in SCZ using genome alignments (225 reads/million in cases versus 200 reads per million in controls, −log10(p-value) = 3.82) but non-significantly decreased using transcript alignments (138 reads/million in cases versus 147 reads per million in controls). Inspection of genomic read alignments revealed the existence of several, putative, unannotated exons. For example 26 reads in SCZ case 5S and 17 reads in SCZ case 5 aligned uniquely to DOCK3 intronic sequences corresponding to nucleotides 5,130,7158–51307328 on Chromosome 3 (data not shown).

Discussion

As the first stage of a genomic convergence analysis of SCZ, we sought gene expression changes between post-mortem cerebellar cortices of 14 patients and 6 controls using shotgun, clonal mRNA sequencing, together with a rigorous and systematic analysis approach. Cerebellum was chosen for study since it has been shown to be affected in SCZ: Patients exhibit impaired sensory-motor integration, reduced cerebellar volume and grey matter, increased cerebellar blood flow and glucose uptake, both at rest and following cognitive challenge, greatest approximate conditional likelihood score in cerebellar vermis, and gene expression changes in the lateral hemispheres of the cerebellar cortex consistent with decreased GABAergic and increased glutamatergic neurotransmission [32], [51], [52], [58]–[66]. In addition, the simple organization and lack of perturbation by antipsychotic medication of cerebellar cortical circuits relative to other brain regions affected by SCZ was appealing for the current study. Specifically, granule cells, which are important in processing sensory information, are by far the most abundant cerebellar glutamatergic neuron, comprising the majority of neurons in the brain [67].

Cases and controls were all male and were matched reasonably well for age, race, cause of death and RNA quality metrics. Since only cases had received psychotropic medication, this variable could not be matched. There is currently considerable excitement about stratification of patients with SCZ on the basis of clinical or neuroimaging endophenotypes [20]. Unfortunately, that was not possible in this series.

Approximately 22 million random, cDNA sequence tags were generated from each of 20 samples, providing the deepest mRNA coverage reported to date and providing an unparalleled assessment of the transcriptional complexity of the human cerebellar cortex. Sequences did not align to 1879 viral genomes, offering negative evidence of chronic viral etiology for SCZ in these patients. 85% of annotated transcripts were expressed (based on molecular counting, Table 2), of which 63% were detected at a level of at least 1-fold coverage (equivalent to 56 reads per million or 1.1 copies per million). This pervasive transcription accords with other recent studies [32], [44], [56], [68]. In contrast to 3′ end tag sequencing, random-primed mRNA sequencing provided good coverage of coding domains (Fig. 1B), allowing assessment of exon utilization in cerebellar cortex. Alternative splicing is thought to generate more than 5 transcripts per human locus, most of which have not yet been annotated [56]. 85% of loci expressed at a level of at least 1-fold coverage showed evidence for utilization of unannotated exons, based on excess genomic read alignments. Proline oxidase, for example, is a very well studied gene that appeared to utilize hitherto unannotated exons in human cerebellum (Fig. 11). For this reason, sequences were aligned both to a reference transcript database (for estimation of expression of well-annotated transcripts) and to the human genome (for estimation of total expression of all exons, whether annotated or not). With additional improvements in library preparation and analysis software it should be possible to measure digital gene expression on an exon-by-exon basis in a strand-specific manner [31]–[50], [69]. The number of transcripts with aligned reads declined greater than exponentially as a function of level of expression (Fig 1B). As a result, 29 million aligned reads provided an impressive linear dynamic range (4.37×log10, Fig. 2). Despite lower read lengths and quality scores in this study (Table 2) than currently obtained with the Illumina GA II instrument, only 9% of reads mapped non-uniquely. Finally, technical imprecision of DTE measurements was paltry (CV 3.4%). In summary, our technical assessment indicated deep mRNA sequencing to be a powerful tool for digital gene expression and empiric annotation of transcribed elements in genomes, in accord with a rapidly growing literature [31]–[50], [69].

Gratifyingly, digital transcript expression values, expressed as aligned read frequencies, were amenable to very similar data transformation, quality control, pattern discovery, row-by-row modeling and annotation analyses as array hybridization datasets. This enabled direct comparison of digital transcript expression and array hybridization data metrics (Figs. 4–9). Like array hybridization signals, aligned read frequencies benefited from log transformation as evidenced by improved Mahalanobis distances, overlaid kernel density estimates, univariate distribution results, and unsupervised principal component analysis. Aligned read frequencies had greater correlation coefficients in pair wise sample comparisons than array hybridization, and much greater ability to distinguish SCZ cases from controls, as evaluated by unsupervised principal component analysis (PCA), Ward hierarchical clustering of Pearson product-moment correlations of read frequencies, and the magnitude of the component of variance attributable to diagnosis. DTE was clearly better than array hybridization for PCA-based segregation of cases and controls. Read frequencies based upon alignment to a reference transcript database and to the human genome were almost indistinguishable by these metrics. Partitioning of variance components of clinical, sample and experimental metadata revealed the largest components of array hybridization variability to be metrics of mRNA quality, in accord with recent studies of post-mortem brain tissue [70]. DTE was much less influenced by RNA quality. Non-diagnosis DTE variability included year of sequence generation, sequencing instrument variation and cause of death (Fig. 8). This is not surprising since substantial improvements occurred in sequencing instrument specifications, reagents and base-calling software during the project. Given continued, rapid evolution of generation II sequencing technologies, it is important to keep instrument specification, reagents and base-calling software relatively stable during mRNA sequencing projects to minimize non-hypothesis-related variability. Importantly, variance decomposition allowed incorporation of large clinical, sample and experimental effects as fixed effects in analyses of variance, minimizing detection of gene expression changes associated with confounding variables. Incorporation of fixed effects in array hybridization, but not DTE, analyses of variance resulted in no expression changes meeting significance cutoffs. Coefficients of correlation between aligned read frequencies, adjusted for transcript length, and array hybridization signals were 0.35 for all values and 0.46 for fold-change in genes with significant differences in expression. These values are somewhat lower than similar comparisons between array hybridization platforms and are probably reflective of the exacting RNA source [71]. It should be noted that aligned read frequencies correlate very well with quantitative PCR results [57], [72]. In conclusion, mRNA sequencing appears to be superior to array hybridization for gene expression analysis and is amenable to standard statistical procedures for quality assessment and gene expression analysis.

Expression of 215 genes differed significantly between cases and controls in post-mortem cerebellar cortex by analysis of variance (Table 3). Of these, 88 were identified in genome-aligned reads and 152 in transcript-aligned reads. Only 25 genes were common to both alignments, but 96% of the 215 had congruent direction of change and correlation between genome- and transcript-aligned fold-change was 0.73. Major contributors to disparity between alignments appeared to be detection of unannotated exons in genome alignments and forced mis-alignments in incomplete transcript datasets. The latter may largely be avoided by removal of non-unique alignments from calculations. As RefSeq transcript becomes more complete, this disparity should decrease. In the interim, it is wise to align both to a well-annotated transcript dataset and to the genome.

17 of the 215 genes with altered cerebellar expression have previously shown changes in SCZ in the dorsolateral prefrontal cortex or superior temporal gyrus (APBA2, BTG1, CACNG2, CAP2, DKFZP434A0131, GABRA1, GOLGA1, HSPBP1, KIAA0256, KPNA1, RPL3, SLC35A3, SLC39A7, SRRM2, TCF4, ZNF148 and ZNF195). Despite substantial differences in gene expression in cerebellar cortex and cerebral cortex [73], 13 genes were congruent with the direction of change reported previously. Three genes involved in GABAergic neurotransmission (GABRA6, GABRB3 and SLC6A1) showed changes in cerebellar cortex in SCZ that agreed with a previous report [52]. In the context of cerebellar cortical function in SCZ, alteration in expression of genes may be causal, consequential, compensatory, or the result of confounding factors. Because the substantial heritability of schizophrenia appears to be polygenic [1]–[5], expression differences might reflect SCZ-associated nucleotide variants that alter expression in cis (e.g. eSNPs). Indeed, seven of these genes have previously shown association with SCZ in genetic association studies (CACNG2, GABRA1, GPSN2, HIRA, PSAP, RANBP5 and TCF4) [74]. Alternatively, expression differences may represent compensatory changes in pathways and networks. GO annotation of genes with altered expression revealed over-representation of membrane association, Zinc binding or transport, regulation of transcription, Golgi apparatus and vesicle-mediated transport.

Most striking were 23 genes involved in presynaptic vesicular transport / Golgi apparatus or post-synaptic neurotransmission, 15 of which were up-regulated and 8 down-regulated. Up-regulated genes included nine involved in transport from the trans-Golgi network to the synaptic vesicle (golgi autoantigen A1, UDP-N-acetylglucosamine transporter (SLC35A3), component of oligomeric golgi complex 6, thyroid hormone receptor interactor 11, adaptor-related protein complex 1G1, ADP-ribosylation factor guanine nucleotide-exchange factor 2, vesicle docking protein p115, Rho-associated, coiled-coil containing protein kinase 1 (ROCK1) and RAB9B), one involved in synaptic vesicle exocytosis (synaptotagmin I), one involved in clathrin-mediated endocytosis (AP2 associated kinase 1 (AAK1)), one involved in homotypic fusion of early endosomes (early endosome antigen 1 (EEA1), and one involved in retrograde transport to the Golgi apparatus (vacuolar protein sorting 35 (VPS35)). Down-regulated genes included two involved in transport from the trans-Golgi network to the synaptic vesicle (syntaxin 10 and ADP-ribosylation factor related protein 1), two involved with synaptic vesicle exocytosis (synaptic vesicle glycoprotein 2A and neurochondrin), two involved in transport from the endoplasmic reticulum to the Golgi (golgi phosphoprotein 2 and synaptic glycoprotein 2) and one involved in retrograde transport to the Golgi apparatus (sorting nexin 17) (Fig. 10, Table 3). Up-regulated post-synaptic membrane genes included GABRA1 (see above) and CACNG2 (calcium channel, voltage-dependent, gamma subunit 2). Several other genes involved in GABAergic neurotransmission showed non-significant changes in expression in cerebellar cortex in SCZ, in agreement with previous reports [52]. One down-regulated post-synaptic membrane gene was ligand-gated ion channel, zinc activated 1. Most of the up-regulated genes, but none of the down-regulated genes in this set, had change in expression of greater than 30%. The largest alterations were CACNG2 (+169%), EEA1 (+88%), ROCK1 (+85%), SLC35A3 (+69%), AAK1 (64%) and GABRA1 and VPS35 (both +52%). Several previous studies have shown decreased expression of genes involved in presynaptic vesicular transport in the dorsolateral prefrontal cortex and hippocampus in SCZ [75]–[78]. While the molecular determinants of synaptic vesicular transport have been intensively studied, it is not yet possible to integrate these transcriptional changes into a single biological outcome. Isolated elevation of AAK1, for example, should decrease neurotransmitter endocytosis [79], while elevated EEA1 should increase early endosomal fusion [80]. However, activity of most of these proteins is dependent either upon post-translational modification or presence of interacting proteins. Furthermore, there is likely to be heterogeneous expression of these proteins among cerebellar cortical neurons. Nevertheless, it appears clear that presynaptic vesicular transport, Golgi function and GABAergic neurotransmission are perturbed in cerebellar cortex in SCZ, in agreement with previous analyses of gene expression in post-mortem brain in SCZ [28], [52], [75] and genetic associations of nucleotide variants in synaptic vesicular transport genes with SCZ (namely CACNG2, COG2, STX1A, SNAP29, MUTED, NRG1, PLDN, NRG3, HOMER3, RTN4R, NRG2, RTN4, CPLX2, BLOC1S1, SYNJ1, HOMER1, CLINT1, SYP, BLOC1S2, SYN3, SYT11, SYNGR1, SNAPAP, STX7, BLOC1S3, DISC1, SNAP25, DLG4, CPLX1, CNO, SYN2, DTNBP1 and DAOA) [74]. Replication of alterations in presynaptic vesicular transport and Golgi function, particularly as affecting GABAergic neurotransmission, in additional cases or other brain regions–such as prefrontal cortex–would substantiate novel molecular mechanisms underpinning SCZ and establish novel targets for therapeutic intervention.

The current study reports the first phase of a genomic convergence analysis of SCZ, which is intended to integrate linkage and expression studies. In a subsequent manuscript, we will present results of expressed nucleotide variants that differ in allele frequency between cases and controls, and integrate nucleotide variant and expression analyses. In addition, a recent study has shown substantial variation in alternative splice isoform expression and alternative polyadenylation in cerebellar cortex between normal individuals [57], and it will be interesting to ascertain whether specific examples of alternative isoform expression show association with SCZ. Confirmation of causal components of the molecular basis of SCZ is anticipated to have significant impact on clinical practice, particularly with respect to timely diagnosis and prognostic or therapeutic categorization [12], [19].

Materials and Methods

Sample preparation

was as previously described [51]. Anonymized cerebellar samples were provided from the Maryland Brain Collection with permission from the Maryland Brain Collection Steering Committee. Permission for the study of brain tissues was provided by families post-mortem in accordance with the guidelines of the Uniform Anatomical Gift Act. Fourteen samples were from patients with a diagnosis of SCZ according to DSM-IV criteria and six were controls. Cortical areas corresponding to crus I/VIIa of cerebellar hemispheres were dissected at −20°C and frozen at −80°C, as described [52]. The average pH of the samples was 6.6±0.17. Cerebellum specimens were obtained according to NIH guidelines for confidentially and privacy. Genomic DNA and total RNA were isolated from samples using standard techniques (Qiagen, Valencia, CA).

Sequencing-by-synthesis

Cerebellum samples were sequenced using Illumina Genome Analyzers with modifications for mRNA samples [39], [42], [43], [52], [57]: Following quality assessment using a Bioanalyzer 2100 (Agilent Inc., Santa Clara, CA; Table 1), poly A+ RNA was isolated from 5–10 µg total RNA by two rounds of oligo-dT selection (Invitrogen Inc., Santa Clara, CA). mRNA was annealed to high concentrations of random hexamers and reverse transcribed. Following second strand synthesis, end repair, and A-tailing, adapters complementary to sequencing primers were ligated to cDNA fragment ends. Resultant cDNA libraries were size fractionated on agarose gels, 200 bp fragments excised, and amplified by 15 cycles of polymerase chain reaction. Following quality assessment using a Bioanalyzer 2100, single-stranded cDNA-adapter fragments were randomly annealed to the surface of a flow cell in a cluster station (Illumina Inc., San Diego, CA) via primers complementary to the adapters and incubated under conditions fostering annealing of the ends of cDNA-adapter fragments to adjacent complementary primers. Primers with annealed cDNA-adapter fragments were extended with DNA polymerase and unlabelled dNTPs in a solid-phase “bridge amplification”. Resultant double-stranded products were denatured, yielding 2 single stranded fragments. The latter steps were repeated for 35 cycles, generating ∼40 million clusters of clonal amplicons. Subsequently, 32–36 cycles of sequencing-by-synthesis chemistry were performed in Illumina Genome Analyzer instruments with 4 dNTPs featuring cleavable dyes and reversible terminators. Following each base extension, 4 images (one for each nucleotide) are taken upon laser excitation. Incorporation of the next base occurred after removal of the blocked 3′ terminus and fluorescent tag of the previously incorporated nucleotide. Sequences were retained if of high quality (defined as those passing the default quality filtering parameters used in the Illumina GA Pipeline GERALD stage, i.e. clusters with intensities greater than 0.6-times the average of the highest and the sum of the two highest intensities for the first 12 cycles). Furthermore, sequencing runs with poor average nucleotide quality score graphs (associated with a minority of reads passing the default quality filtering parameters) were discarded.

Read Alignment-based gene expression profiling

High-quality reads were aligned to the human genome, Build 36.2, RefSeq transcript database, Release 22 [81], Unigene Hs, build 215 [82] using the algorithm GMAP [55] and the software system Alpheus® [48], with adjustments for short SBS reads (oligomer overlap interval = 3 nt, identities ≥34/36 or 94%). A read was denoted ‘aligned’ to a locus if its sequence alignment to the genomic reference sequence (NCBI build 36.2) fell within the boundaries of the locus coordinates on the chromosome. Locus boundaries on the genome were defined by NCBI annotations, as reported through the Nucleotide database. Annotation data was downloaded on 2/2/2007 using the Batch Entrez query tools (http://www.ncbi.nlm.nih.gov/sites/batchentrez?dbNucleotide). Only the highest scoring alignment(s) was retained. Reads with a single best alignment or with equally good alignments to alternative transcripts of the same locus were considered uniquely aligned. Aligned read frequencies (per million aligned reads) were calculated for each sample and locus using Alpheus® [48].

Array-based gene expression profiling was performed using Affymetrix Human Genome U133 Plus 2.0 oligonucleotide arrays and standard procedures (Affymetrix, Santa Clara, CA). Array-based gene expression profiling of cerebellar cortex of unaffected individuals was performed using Affymetrix Human U95A oligonucleotide arrays as previously reported [83]. Primary image analysis was performed with Affymetrix GENECHIP 3.2. Images were scaled to an average hybridization intensity of 200, and the threshold for expression was an average hybridization intensity of ≥50.

Statistical Analysis

Array hybridization signals and read frequencies were log10 transformed prior to evaluation of inter-sample differences. Overlaid kernel density estimates, univariate distribution results, Mahalanobis distances, correlation coefficients of pair wise sample comparisons, unsupervised principal component analysis (by Pearson product-moment correlation) and Ward hierarchical clustering of Pearson product-moment correlations of read frequencies were performed with JMP Genomics, Version 3.2 (SAS Institute, Cary, NC). Decomposition of principal components of variance of array signals and read frequencies, FDR corrected analysis of variance (with inclusion of the non-diagnosis components of variance discussed above as fixed effects) and chi squared comparisons of GO annotations were also performed with JMP Genomics, Version 3.2. Patient, sample, and experimental parameters examined to quantify sources of variability in measurements were age, alignment percentage, brain pH, cause of death, cluster station, diagnosis, library creator, number of reads, post-mortem interval, psychotropic medication, race, read length, read quality score, RNA isolation date, RNA integrity number, sample storage duration, sequencing instrument and year sequenced.

Acknowledgments

A Deo lumen, ab amicis auxilium. We thank Steven Day, Dan Weems, Melodie J. Rice, Terri L. Gomez, Linda Julien, Kamal Gajendran for assistance with Alpheus software development. We thank Becky Archuleta, Dawn Doyle, Pauline Brito, and Leisha Pak for assistance with variant annotation. We thank Eric Vermaas, Victor Quijano, Juying Yan, Leonda Clendenen and Jimmy Woodward for assistance in generating some of the sequencing data. We thank the NIH Neuroscience DNA Microarray Consortium and Translational Genomics Research Institute (Phoenix, AZ) for help with gene expression profiling of human cerebellar tissue. Irina Khrebtukova, Shujun Luo, Lu Zhang and Gary P. Schroth are employees of Illumina Inc. Study data is available at NCBI Gene Expression Omnibus with accession number GSE12297 and at the Short Read Archive with accession number SRA002355.1 and is freely available for query at http://schizophrenia.ncgr.org.

Author Contributions

Conceived and designed the experiments: GPS SFK. Performed the experiments: JM NAM IK IEL GM JJH LS LZ JCvV ADF RJL RWK. Analyzed the data: JM NAM IK IEL GM JJH JCvV ADF SL WB FDS SMV CFB MKM LCM RJL SKK MG VAG REH WC SM RDW NPB GPS SFK. Contributed reagents/materials/analysis tools: NAM GM ADF SL WB FDS SMV CFB MKM LCM JPU RWK RCR WC SM NPB GPS. Wrote the paper: JM SKK MG REH RDW NPB GPS SFK.

Citation: Mudge J, Miller NA, Khrebtukova I, Lindquist IE, May GD, Huntley JJ, et al. (2008) Genomic Convergence Analysis of Schizophrenia: mRNA Sequencing Reveals Altered Synaptic Vesicular Transport in Post-Mortem Cerebellum. PLoS ONE3(11): e3625. https://doi.org/10.1371/journal.pone.0003625

References

1. Faraone SV, Seidman LJ, Kremen WS, Toomey R, Pepple JR, et al. (2000) Neuropsychologic functioning among the nonpsychotic relatives of schizophrenic patients: the effect of genetic loading. Biol Psychiatry 48: 120–126.SV FaraoneLJ SeidmanWS KremenR. ToomeyJR Pepple2000Neuropsychologic functioning among the nonpsychotic relatives of schizophrenic patients: the effect of genetic loading.Biol Psychiatry48120126

2. McGuffin P, Owen MJ, O'Donovan MC, Thapar A, Gottesman (1994) London: Geskell Press. P. McGuffinMJ OwenMC O'DonovanA. ThaparGottesman II1994LondonGeskell Press

3. Lichtermann D, Karbe E, Maier W (2000) The genetic epidemiology of schizophrenia and of schizophrenia spectrum disorders. Eur Arch Psychiatry Clin Neurosci 250: 304–310.D. LichtermannE. KarbeW. Maier2000The genetic epidemiology of schizophrenia and of schizophrenia spectrum disorders.Eur Arch Psychiatry Clin Neurosci250304310

4. Ross CA, Margolis RL, Reading SA, Pletnikov M, Coyle JT (2006) Neurobiology of schizophrenia. Neuron 52: 139–153.CA RossRL MargolisSA ReadingM. PletnikovJT Coyle2006Neurobiology of schizophrenia.Neuron52139153

5. Saha S, Chant D, Welham J, McGrath J (2005) A systematic review of the prevalence of schizophrenia. PLoS Med 2: e141.S. SahaD. ChantJ. WelhamJ. McGrath2005A systematic review of the prevalence of schizophrenia.PLoS Med2e141

6. Abusaad I, MacKay D, Zhao J, Stanford P, Collier DA, et al. (1999) Stereological estimation of the total number of neurons in the murine hippocampus using the optical disector. J Comp Neurol 408: 560–566.I. AbusaadD. MacKayJ. ZhaoP. StanfordDA Collier1999Stereological estimation of the total number of neurons in the murine hippocampus using the optical disector.J Comp Neurol408560566

7. Harrison PJ, Weinberger DR (2005) Schizophrenia genes, gene expression, and neuropathology: on the matter of their convergence. Mol Psychiatry 10: 40–68; image 45.PJ HarrisonDR Weinberger2005Schizophrenia genes, gene expression, and neuropathology: on the matter of their convergence.Mol Psychiatry104068; image 45

8. Riley B, Kendler KS (2006) Molecular genetic studies of schizophrenia. Eur J Hum Genet 14: 669–680.B. RileyKS Kendler2006Molecular genetic studies of schizophrenia.Eur J Hum Genet14669680

9. (1994) Diagnostic and statistical manual of mental disorders. Washington, DC: American Psychiatric Association. 1994Diagnostic and statistical manual of mental disordersWashington, DCAmerican Psychiatric Association

10. Cardon LR, Palmer LJ (2003) Population stratification and spurious allelic association. Lancet 361: 598–604.LR CardonLJ Palmer2003Population stratification and spurious allelic association.Lancet361598604

11. Sillanpaa MJ, Auranen K (2004) Replication in genetic studies of complex traits. Ann Hum Genet 68: 646–657.MJ SillanpaaK. Auranen2004Replication in genetic studies of complex traits.Ann Hum Genet68646657

12. Tsuang MT, Stone WS, Faraone SV (2000) Toward reformulating the diagnosis of schizophrenia. Am J Psychiatry 157: 1041–1050.MT TsuangWS StoneSV Faraone2000Toward reformulating the diagnosis of schizophrenia.Am J Psychiatry15710411050

13. Lewis CM, Levinson DF, Wise LH, DeLisi LE, Straub RE, et al. (2003) Genome scan meta-analysis of schizophrenia and bipolar disorder, part II: Schizophrenia. Am J Hum Genet 73: 34–48.CM LewisDF LevinsonLH WiseLE DeLisiRE Straub2003Genome scan meta-analysis of schizophrenia and bipolar disorder, part II: Schizophrenia.Am J Hum Genet733448

14. Rutter M, Moffitt TE, Caspi A (2006) Gene-environment interplay and psychopathology: multiple varieties but real effects. J Child Psychol Psychiatry 47: 226–261.M. RutterTE MoffittA. Caspi2006Gene-environment interplay and psychopathology: multiple varieties but real effects.J Child Psychol Psychiatry47226261

15. Carter CJ (2006) Schizophrenia susceptibility genes converge on interlinked pathways related to glutamatergic transmission and long-term potentiation, oxidative stress and oligodendrocyte viability. Schizophr Res 86: 1–14.CJ Carter2006Schizophrenia susceptibility genes converge on interlinked pathways related to glutamatergic transmission and long-term potentiation, oxidative stress and oligodendrocyte viability.Schizophr Res86114

16. Crespi B, Summers K, Dorus S (2007) Adaptive evolution of genes underlying schizophrenia. Proc Biol Sci 274: 2801–2810.B. CrespiK. SummersS. Dorus2007Adaptive evolution of genes underlying schizophrenia.Proc Biol Sci27428012810

17. Kirov G, O'Donovan MC, Owen MJ (2005) Finding schizophrenia genes. J Clin Invest 115: 1440–1448.G. KirovMC O'DonovanMJ Owen2005Finding schizophrenia genes.J Clin Invest11514401448

18. Owen MJ, Craddock N, O'Donovan MC (2005) Schizophrenia: genes at last? Trends Genet 21: 518–525.MJ OwenN. CraddockMC O'Donovan2005Schizophrenia: genes at last?Trends Genet21518525

19. Kendler KS, Neale MC, Walsh D (1995) Evaluating the spectrum concept of schizophrenia in the Roscommon Family Study. Am J Psychiatry 152: 749–754.KS KendlerMC NealeD. Walsh1995Evaluating the spectrum concept of schizophrenia in the Roscommon Family Study.Am J Psychiatry152749754

20. Gottesman , Gould TD (2003) The endophenotype concept in psychiatry: etymology and strategic intentions. Am J Psychiatry 160: 636–645.Gottesman IITD Gould2003The endophenotype concept in psychiatry: etymology and strategic intentions.Am J Psychiatry160636645

21. Glatt SJ, Everall IP, Kremen WS, Corbeil J, Sasik R, et al. (2005) Comparative gene expression analysis of blood and brain provides concurrent validation of SELENBP1 up-regulation in schizophrenia. Proc Natl Acad Sci U S A 102: 15533–15538.SJ GlattIP EverallWS KremenJ. CorbeilR. Sasik2005Comparative gene expression analysis of blood and brain provides concurrent validation of SELENBP1 up-regulation in schizophrenia.Proc Natl Acad Sci U S A1021553315538

22. Olsen L, Hansen T, Jakobsen KD, Djurovic S, Melle I, et al. (2008) The estrogen hypothesis of schizophrenia implicates glucose metabolism: association study in three independent samples. BMC Med Genet 9: 39.L. OlsenT. HansenKD JakobsenS. DjurovicI. Melle2008The estrogen hypothesis of schizophrenia implicates glucose metabolism: association study in three independent samples.BMC Med Genet939

23. Cheng C, Pounds S (2007) False discovery rate paradigms for statistical analyses of microarray gene expression data. Bioinformation 1: 436–446.C. ChengS. Pounds2007False discovery rate paradigms for statistical analyses of microarray gene expression data.Bioinformation1436446

24. Hauser MA, Li YJ, Takeuchi S, Walters R, Noureddine M, et al. (2003) Genomic convergence: identifying candidate genes for Parkinson's disease by combining serial analysis of gene expression and genetic linkage. Hum Mol Genet 12: 671–677.MA HauserYJ LiS. TakeuchiR. WaltersM. Noureddine2003Genomic convergence: identifying candidate genes for Parkinson's disease by combining serial analysis of gene expression and genetic linkage.Hum Mol Genet12671677

25. Li YJ, Oliveira SA, Xu P, Martin ER, Stenger JE, et al. (2003) Glutathione S-transferase omega-1 modifies age-at-onset of Alzheimer disease and Parkinson disease. Hum Mol Genet 12: 3259–3267.YJ LiSA OliveiraP. XuER MartinJE Stenger2003Glutathione S-transferase omega-1 modifies age-at-onset of Alzheimer disease and Parkinson disease.Hum Mol Genet1232593267

26. Noureddine MA, Li YJ, van der Walt JM, Walters R, Jewett RM, et al. (2005) Genomic convergence to identify candidate genes for Parkinson disease: SAGE analysis of the substantia nigra. Mov Disord 20: 1299–1309.MA NoureddineYJ LiJM van der WaltR. WaltersRM Jewett2005Genomic convergence to identify candidate genes for Parkinson disease: SAGE analysis of the substantia nigra.Mov Disord2012991309

27. Oliveira SA, Li YJ, Noureddine MA, Zuchner S, Qin X, et al. (2005) Identification of risk and age-at-onset genes on chromosome 1p in Parkinson disease. Am J Hum Genet 77: 252–264.SA OliveiraYJ LiMA NoureddineS. ZuchnerX. Qin2005Identification of risk and age-at-onset genes on chromosome 1p in Parkinson disease.Am J Hum Genet77252264

28. Le-Niculescu H, Balaraman Y, Patel S, Tan J, Sidhu K, et al. (2007) Towards understanding the schizophrenia code: an expanded convergent functional genomics approach. Am J Med Genet B Neuropsychiatr Genet 144B: 129–158.H. Le-NiculescuY. BalaramanS. PatelJ. TanK. Sidhu2007Towards understanding the schizophrenia code: an expanded convergent functional genomics approach.Am J Med Genet B Neuropsychiatr Genet144B129158

29. Lipska BK, Mitkus SN, Mathew SV, Fatula R, Hyde TM, et al. (2006) Functional genomics in postmortem human brain: abnormalities in a DISC1 molecular pathway in schizophrenia. Dialogues Clin Neurosci 8: 353–357.BK LipskaSN MitkusSV MathewR. FatulaTM Hyde2006Functional genomics in postmortem human brain: abnormalities in a DISC1 molecular pathway in schizophrenia.Dialogues Clin Neurosci8353357

30. Middleton FA, Pato CN, Gentile KL, McGann L, Brown AM, et al. (2005) Gene expression analysis of peripheral blood leukocytes from discordant sib-pairs with schizophrenia and bipolar disorder reveals points of convergence between genetic and functional genomic approaches. Am J Med Genet B Neuropsychiatr Genet 136B: 12–25.FA MiddletonCN PatoKL GentileL. McGannAM Brown2005Gene expression analysis of peripheral blood leukocytes from discordant sib-pairs with schizophrenia and bipolar disorder reveals points of convergence between genetic and functional genomic approaches.Am J Med Genet B Neuropsychiatr Genet136B1225

31. Bainbridge MN, Warren RL, Hirst M, Romanuik T, Zeng T, et al. (2006) Analysis of the prostate cancer cell line LNCaP transcriptome using a sequencing-by-synthesis approach. BMC Genomics 7: 246.MN BainbridgeRL WarrenM. HirstT. RomanuikT. Zeng2006Analysis of the prostate cancer cell line LNCaP transcriptome using a sequencing-by-synthesis approach.BMC Genomics7246

32. Kim JB, Porreca GJ, Song L, Greenway SC, Gorham JM, et al. (2007) Polony multiplex analysis of gene expression (PMAGE) in mouse hypertrophic cardiomyopathy. Science 316: 1481–1484.JB KimGJ PorrecaL. SongSC GreenwayJM Gorham2007Polony multiplex analysis of gene expression (PMAGE) in mouse hypertrophic cardiomyopathy.Science31614811484

33. Butz JA, Roberts KG, Edwards JS (2004) Detecting changes in the relative expression of KRAS2 splice variants using polymerase colonies. Biotechnol Prog 20: 1836–1839.JA ButzKG RobertsJS Edwards2004Detecting changes in the relative expression of KRAS2 splice variants using polymerase colonies.Biotechnol Prog2018361839

34. Butz JA, Yan H, Mikkilineni V, Edwards JS (2004) Detection of allelic variations of human gene expression by polymerase colonies. BMC Genet 5: 3.JA ButzH. YanV. MikkilineniJS Edwards2004Detection of allelic variations of human gene expression by polymerase colonies.BMC Genet53

35. Cheung F, Haas BJ, Goldberg SM, May GD, Xiao Y, et al. (2006) Sequencing Medicago truncatula expressed sequenced tags using 454 Life Sciences technology. BMC Genomics 7: 272.F. CheungBJ HaasSM GoldbergGD MayY. Xiao2006Sequencing Medicago truncatula expressed sequenced tags using 454 Life Sciences technology.BMC Genomics7272

36. Chiu KP, Ariyaratne P, Xu H, Tan A, Ng P, et al. (2007) Pathway aberrations of murine melanoma cells observed in Paired-End diTag transcriptomes. BMC Cancer 7: 109.KP ChiuP. AriyaratneH. XuA. TanP. Ng2007Pathway aberrations of murine melanoma cells observed in Paired-End diTag transcriptomes.BMC Cancer7109

37. Emrich SJ, Barbazuk WB, Li L, Schnable PS (2007) Gene discovery and annotation using LCM-454 transcriptome sequencing. Genome Res 17: 69–73.SJ EmrichWB BarbazukL. LiPS Schnable2007Gene discovery and annotation using LCM-454 transcriptome sequencing.Genome Res176973

38. Jones-Rhoades MW, Borevitz JO, Preuss D (2007) Genome-wide expression profiling of the Arabidopsis female gametophyte identifies families of small, secreted proteins. PLoS Genet 3: 1848–1861.MW Jones-RhoadesJO BorevitzD. Preuss2007Genome-wide expression profiling of the Arabidopsis female gametophyte identifies families of small, secreted proteins.PLoS Genet318481861

39. Lister R, O'Malley RC, Tonti-Filippini J, Gregory BD, Berry CC, et al. (2008) Highly integrated single-base resolution maps of the epigenome in Arabidopsis. Cell 133: 523–536.R. ListerRC O'MalleyJ. Tonti-FilippiniBD GregoryCC Berry2008Highly integrated single-base resolution maps of the epigenome in Arabidopsis.Cell133523536

40. Marioni J, Mason C, Mane S, Stephens M, Gilad Y (2008) RNA-seq: An assessment of technical reproducibility and comparison with gene expression arrays. Genome Res. J. MarioniC. MasonS. ManeM. StephensY. Gilad2008RNA-seq: An assessment of technical reproducibility and comparison with gene expression arrays.Genome Res

41. Mikkilineni V, Mitra RD, Merritt J, DiTonno JR, Church GM, et al. (2004) Digital quantitative measurements of gene expression. Biotechnol Bioeng 86: 117–124.V. MikkilineniRD MitraJ. MerrittJR DiTonnoGM Church2004Digital quantitative measurements of gene expression.Biotechnol Bioeng86117124

42. Morin RD, O'Connor MD, Griffith M, Kuchenbauer F, Delaney A, et al. (2008) Application of massively parallel sequencing to microRNA profiling and discovery in human embryonic stem cells. Genome Res 18: 610–621.RD MorinMD O'ConnorM. GriffithF. KuchenbauerA. Delaney2008Application of massively parallel sequencing to microRNA profiling and discovery in human embryonic stem cells.Genome Res18610621

43. Mortazavi A, Williams BA, McCue K, Schaeffer L, Wold B (2008) Mapping and quantifying mammalian transcriptomes by RNA-Seq. Nat Methods 5: 621–628.A. MortazaviBA WilliamsK. McCueL. SchaefferB. Wold2008Mapping and quantifying mammalian transcriptomes by RNA-Seq.Nat Methods5621628

44. Nagalakshmi U, Wang Z, Waern K, Shou C, Raha D, et al. (2008) The transcriptional landscape of the yeast genome defined by RNA sequencing. Science 320: 1344–1349.U. NagalakshmiZ. WangK. WaernC. ShouD. Raha2008The transcriptional landscape of the yeast genome defined by RNA sequencing.Science32013441349

45. Ng P, Tan JJ, Ooi HS, Lee YL, Chiu KP, et al. (2006) Multiplex sequencing of paired-end ditags (MS-PET): a strategy for the ultra-high-throughput analysis of transcriptomes and genomes. Nucleic Acids Res 34: e84.P. NgJJ TanHS OoiYL LeeKP Chiu2006Multiplex sequencing of paired-end ditags (MS-PET): a strategy for the ultra-high-throughput analysis of transcriptomes and genomes.Nucleic Acids Res34e84

46. Nielsen KL, Hogh AL, Emmersen J (2006) DeepSAGE–digital transcriptomics with high sensitivity, simple experimental protocol and multiplexing of samples. Nucleic Acids Res 34: e133.KL NielsenAL HoghJ. Emmersen2006DeepSAGE–digital transcriptomics with high sensitivity, simple experimental protocol and multiplexing of samples.Nucleic Acids Res34e133

47. Ohtsu K, Smith MB, Emrich SJ, Borsuk LA, Zhou R, et al. (2007) Global gene expression analysis of the shoot apical meristem of maize (Zea mays L.). Plant J 52: 391–404.K. OhtsuMB SmithSJ EmrichLA BorsukR. Zhou2007Global gene expression analysis of the shoot apical meristem of maize (Zea mays L.).Plant J52391404

48. Sugarbaker DJ, Richards WG, Gordon GJ, Dong L, De Rienzo A, et al. (2008) Transcriptome sequencing of malignant pleural mesothelioma tumors. Proc Natl Acad Sci U S A 105: 3521–3526.DJ SugarbakerWG RichardsGJ GordonL. DongA. De Rienzo2008Transcriptome sequencing of malignant pleural mesothelioma tumors.Proc Natl Acad Sci U S A10535213526

49. Toth AL, Varala K, Newman TC, Miguez FE, Hutchison SK, et al. (2007) Wasp gene expression supports an evolutionary link between maternal behavior and eusociality. Science 318: 441–444.AL TothK. VaralaTC NewmanFE MiguezSK Hutchison2007Wasp gene expression supports an evolutionary link between maternal behavior and eusociality.Science318441444

50. Weber AP, Weber KL, Carr K, Wilkerson C, Ohlrogge JB (2007) Sampling the Arabidopsis transcriptome with massively parallel pyrosequencing. Plant Physiol 144: 32–42.AP WeberKL WeberK. CarrC. WilkersonJB Ohlrogge2007Sampling the Arabidopsis transcriptome with massively parallel pyrosequencing.Plant Physiol1443242

51. Paz RD, Andreasen NC, Daoud SZ, Conley R, Roberts R, et al. (2006) Increased expression of activity-dependent genes in cerebellar glutamatergic neurons of patients with schizophrenia. Am J Psychiatry 163: 1829–1831.RD PazNC AndreasenSZ DaoudR. ConleyR. Roberts2006Increased expression of activity-dependent genes in cerebellar glutamatergic neurons of patients with schizophrenia.Am J Psychiatry16318291831

52. Bullock WM, Cardon K, Bustillo J, Roberts R, Perrone-Bizzozero NI (2008) Altered expression of genes involved in GABAergic transmission and neuromodulation of granule cell activity in the cerebellum of patients with schizophrenia. Am J Psychiatry In press. WM BullockK. CardonJ. BustilloR. RobertsNI Perrone-Bizzozero2008Altered expression of genes involved in GABAergic transmission and neuromodulation of granule cell activity in the cerebellum of patients with schizophrenia.Am J Psychiatry In press

53. Freimer N, Sabatti C (2004) The use of pedigree, sib-pair and association studies of common diseases for genetic mapping and epidemiology. Nat Genet 36: 1045–1051.N. FreimerC. Sabatti2004The use of pedigree, sib-pair and association studies of common diseases for genetic mapping and epidemiology.Nat Genet3610451051

54. McClellan JM, Susser E, King MC (2007) Schizophrenia: a common disease caused by multiple rare alleles. Br J Psychiatry 190: 194–199.JM McClellanE. SusserMC King2007Schizophrenia: a common disease caused by multiple rare alleles.Br J Psychiatry190194199

55. Wu TD, Watanabe CK (2005) GMAP: a genomic mapping and alignment program for mRNA and EST sequences. Bioinformatics 21: 1859–1875.TD WuCK Watanabe2005GMAP: a genomic mapping and alignment program for mRNA and EST sequences.Bioinformatics2118591875

56. Birney E, Stamatoyannopoulos JA, Dutta A, Guigo R, Gingeras TR, et al. (2007) Identification and analysis of functional elements in 1% of the human genome by the ENCODE pilot project. Nature 447: 799–816.E. BirneyJA StamatoyannopoulosA. DuttaR. GuigoTR Gingeras2007Identification and analysis of functional elements in 1% of the human genome by the ENCODE pilot project.Nature447799816

57. Wang ET, Sandberg R, Luo S, Khrebtukova I, Zhang L, et al. (2008) Alternative Isoform Regulation in Human Tissue Transcriptomes. ET WangR. SandbergS. LuoI. KhrebtukovaL. Zhang2008Alternative Isoform Regulation in Human Tissue Transcriptomes.Nature Oct 22 Advance Online Publication. Nature Oct 22 Advance Online Publication.

58. Andreasen NC, Calage CA, O'Leary DS (2008) Theory of mind and schizophrenia: a positron emission tomography study of medication-free patients. Schizophr Bull 34: 708–719.NC AndreasenCA CalageDS O'Leary2008Theory of mind and schizophrenia: a positron emission tomography study of medication-free patients.Schizophr Bull34708719

59. Andreasen NC, Pierson R (2008) The role of the cerebellum in schizophrenia. Biol Psychiatry 64: 81–88.NC AndreasenR. Pierson2008The role of the cerebellum in schizophrenia.Biol Psychiatry648188

60. Bottmer C, Bachmann S, Pantel J, Essig M, Amann M, et al. (2005) Reduced cerebellar volume and neurological soft signs in first-episode schizophrenia. Psychiatry Res 140: 239–250.C. BottmerS. BachmannJ. PantelM. EssigM. Amann2005Reduced cerebellar volume and neurological soft signs in first-episode schizophrenia.Psychiatry Res140239250

61. Edwards CR, Newman S, Bismark A, Skosnik PD, O'Donnell BF, et al. (2008) Cerebellum volume and eyeblink conditioning in schizophrenia. Psychiatry Res 162: 185–194.CR EdwardsS. NewmanA. BismarkPD SkosnikBF O'Donnell2008Cerebellum volume and eyeblink conditioning in schizophrenia.Psychiatry Res162185194

62. Ho BC, Mola C, Andreasen NC (2004) Cerebellar dysfunction in neuroleptic naive schizophrenia patients: clinical, cognitive, and neuroanatomic correlates of cerebellar neurologic signs. Biol Psychiatry 55: 1146–1153.BC HoC. MolaNC Andreasen2004Cerebellar dysfunction in neuroleptic naive schizophrenia patients: clinical, cognitive, and neuroanatomic correlates of cerebellar neurologic signs.Biol Psychiatry5511461153

63. Malaspina D, Harkavy-Friedman J, Corcoran C, Mujica-Parodi L, Printz D, et al. (2004) Resting neural activity distinguishes subgroups of schizophrenia patients. Biol Psychiatry 56: 931–937.D. MalaspinaJ. Harkavy-FriedmanC. CorcoranL. Mujica-ParodiD. Printz2004Resting neural activity distinguishes subgroups of schizophrenia patients.Biol Psychiatry56931937

64. Okugawa G, Nobuhara K, Takase K, Kinoshita T (2007) Cerebellar posterior superior vermis and cognitive cluster scores in drug-naive patients with first-episode schizophrenia. Neuropsychobiology 56: 216–219.G. OkugawaK. NobuharaK. TakaseT. Kinoshita2007Cerebellar posterior superior vermis and cognitive cluster scores in drug-naive patients with first-episode schizophrenia.Neuropsychobiology56216219

65. Potkin SG, Alva G, Fleming K, Anand R, Keator D, et al. (2002) A PET study of the pathophysiology of negative symptoms in schizophrenia. Positron emission tomography. Am J Psychiatry 159: 227–237.SG PotkinG. AlvaK. FlemingR. AnandD. Keator2002A PET study of the pathophysiology of negative symptoms in schizophrenia. Positron emission tomography.Am J Psychiatry159227237

66. Puri BK, Counsell SJ, Saeed N, Bustos MG, Treasaden IH, et al. (2008) Regional grey matter volumetric changes in forensic schizophrenia patients: an MRI study comparing the brain structure of patients who have seriously and violently offended with that of patients who have not. BMC Psychiatry 8: Suppl 1S6.BK PuriSJ CounsellN. SaeedMG BustosIH Treasaden2008Regional grey matter volumetric changes in forensic schizophrenia patients: an MRI study comparing the brain structure of patients who have seriously and violently offended with that of patients who have not.BMC Psychiatry8Suppl 1S6

67. Chadderton P, Margrie TW, Hausser M (2004) Integration of quanta in cerebellar granule cells during sensory processing. Nature 428: 856–860.P. ChaddertonTW MargrieM. Hausser2004Integration of quanta in cerebellar granule cells during sensory processing.Nature428856860

68. Kapranov P, Cheng J, Dike S, Nix DA, Duttagupta R, et al. (2007) RNA maps reveal new RNA classes and a possible function for pervasive transcription. Science 316: 1484–1488.P. KapranovJ. ChengS. DikeDA NixR. Duttagupta2007RNA maps reveal new RNA classes and a possible function for pervasive transcription.Science31614841488

69. de Hoon M, Hayashizaki Y (2008) Deep cap analysis gene expression (CAGE): genome-wide identification of promoters, quantification of their expression, and network inference. Biotechniques 44: 627–628.630, 632M. de HoonY. Hayashizaki2008Deep cap analysis gene expression (CAGE): genome-wide identification of promoters, quantification of their expression, and network inference.Biotechniques44627628630, 632

70. Popova T, Mennerich D, Weith A, Quast K (2008) Effect of RNA quality on transcript intensity levels in microarray analysis of human post-mortem brain tissues. BMC Genomics 9: 91.T. PopovaD. MennerichA. WeithK. Quast2008Effect of RNA quality on transcript intensity levels in microarray analysis of human post-mortem brain tissues.BMC Genomics991

71. Chen JJ, Hsueh HM, Delongchamp RR, Lin CJ, Tsai CA (2007) Reproducibility of microarray data: a further analysis of microarray quality control (MAQC) data. BMC Bioinformatics 8: 412.JJ ChenHM HsuehRR DelongchampCJ LinCA Tsai2007Reproducibility of microarray data: a further analysis of microarray quality control (MAQC) data.BMC Bioinformatics8412

72. Canales RD, Luo Y, Willey JC, Austermiller B, Barbacioru CC, et al. (2006) Evaluation of DNA microarray results with quantitative gene expression platforms. Nat Biotechnol 24: 1115–1122.RD CanalesY. LuoJC WilleyB. AustermillerCC Barbacioru2006Evaluation of DNA microarray results with quantitative gene expression platforms.Nat Biotechnol2411151122

73. Evans SJ, Choudary PV, Vawter MP, Li J, Meador-Woodruff JH, et al. (2003) DNA microarray analysis of functionally discrete human brain regions reveals divergent transcriptional profiles. Neurobiol Dis 14: 240–250.SJ EvansPV ChoudaryMP VawterJ. LiJH Meador-Woodruff2003DNA microarray analysis of functionally discrete human brain regions reveals divergent transcriptional profiles.Neurobiol Dis14240250

74. Allen NC, Bagade S, McQueen MB, Ioannidis JP, Kavvoura FK, et al. (2008) Systematic meta-analyses and field synopsis of genetic association studies in schizophrenia: the SzGene database. Nat Genet 40: 827–834.NC AllenS. BagadeMB McQueenJP IoannidisFK Kavvoura2008Systematic meta-analyses and field synopsis of genetic association studies in schizophrenia: the SzGene database.Nat Genet40827834

75. Mirnics K, Middleton FA, Marquez A, Lewis DA, Levitt P (2000) Molecular characterization of schizophrenia viewed by microarray analysis of gene expression in prefrontal cortex. Neuron 28: 53–67.K. MirnicsFA MiddletonA. MarquezDA LewisP. Levitt2000Molecular characterization of schizophrenia viewed by microarray analysis of gene expression in prefrontal cortex.Neuron285367

76. Thompson PM, Egbufoama S, Vawter MP (2003) SNAP-25 reduction in the hippocampus of patients with schizophrenia. Prog Neuropsychopharmacol Biol Psychiatry 27: 411–417.PM ThompsonS. EgbufoamaMP Vawter2003SNAP-25 reduction in the hippocampus of patients with schizophrenia.Prog Neuropsychopharmacol Biol Psychiatry27411417

77. Vawter MP, Crook JM, Hyde TM, Kleinman JE, Weinberger DR, et al. (2002) Microarray analysis of gene expression in the prefrontal cortex in schizophrenia: a preliminary study. Schizophr Res 58: 11–20.MP VawterJM CrookTM HydeJE KleinmanDR Weinberger2002Microarray analysis of gene expression in the prefrontal cortex in schizophrenia: a preliminary study.Schizophr Res581120

78. Vawter MP, Thatcher L, Usen N, Hyde TM, Kleinman JE, et al. (2002) Reduction of synapsin in the hippocampus of patients with bipolar disorder and schizophrenia. Mol Psychiatry 7: 571–578.MP VawterL. ThatcherN. UsenTM HydeJE Kleinman2002Reduction of synapsin in the hippocampus of patients with bipolar disorder and schizophrenia.Mol Psychiatry7571578

79. Henderson DM, Conner SD (2007) A novel AAK1 splice variant functions at multiple steps of the endocytic pathway. Mol Biol Cell 18: 2698–2706.DM HendersonSD Conner2007A novel AAK1 splice variant functions at multiple steps of the endocytic pathway.Mol Biol Cell1826982706

80. Haas AK, Fuchs E, Kopajtich R, Barr FA (2005) A GTPase-activating protein controls Rab5 function in endocytic trafficking. Nat Cell Biol 7: 887–893.AK HaasE. FuchsR. KopajtichFA Barr2005A GTPase-activating protein controls Rab5 function in endocytic trafficking.Nat Cell Biol7887893

81. Pruitt KD, Tatusova T, Maglott DR (2005) NCBI Reference Sequence (RefSeq): a curated non-redundant sequence database of genomes, transcripts and proteins. Nucleic Acids Res 33: D501–504.KD PruittT. TatusovaDR Maglott2005NCBI Reference Sequence (RefSeq): a curated non-redundant sequence database of genomes, transcripts and proteins.Nucleic Acids Res33D501504

82. Pontius JU, Wagner L, Schuler GD (2003) UniGene: a unified view of the transcriptome. The NCBI Handbook. Bethesda, MD: National Center for Biotechnology Information. JU PontiusL. WagnerGD Schuler2003UniGene: a unified view of the transcriptome. The NCBI HandbookBethesda, MDNational Center for Biotechnology Information

83. Su AI, Cooke MP, Ching KA, Hakak Y, Walker JR, et al. (2002) Large-scale analysis of the human and mouse transcriptomes. Proc Natl Acad Sci U S A 99: 4465–4470.AI SuMP CookeKA ChingY. HakakJR Walker2002Large-scale analysis of the human and mouse transcriptomes.Proc Natl Acad Sci U S A9944654470