Skip to main content
Full access
New Research
Published Online: 1 July 2014

Prenatal Expression Patterns of Genes Associated With Neuropsychiatric Disorders

Abstract

Objective

Neurodevelopmental disorders presumably involve events that occur during brain development. The authors hypothesized that neuropsychiatric disorders considered to be developmental in etiology are associated with susceptibility genes that are relatively upregulated during fetal life (i.e., differentially expressed).

Method

The authors investigated the presence of prenatal expression enrichment of susceptibility genes systematically, as composite gene sets associated with six neuropsychiatric disorders in the microarray-based “BrainCloud” dorsolateral prefrontal cortex transcriptome.

Results

Using a fetal/postnatal log2-fold change threshold of 0.5, genes associated with syndromic neurodevelopmental disorders (N=31 genes, p=3.37×10–3), intellectual disability (N=88 genes, p=5.53×10–3), and autism spectrum disorder (N=242 genes, p=3.45×10–4) were relatively enriched in prenatal transcript abundance, compared with the overall transcriptome. Genes associated with schizophrenia by genome-wide association studies were not preferentially fetally expressed (N=106 genes, p=0.46), nor were genes associated with schizophrenia by exome sequencing (N=212 genes, p=0.21), but specific genes within copy-number variant regions associated with schizophrenia were relatively enriched in prenatal transcript abundance, and genes associated with schizophrenia by meta-analysis were functionally enriched for some neurodevelopmental processes. In contrast, genes associated with neurodegenerative disorders were significantly underexpressed during fetal life (N=46 genes, p=1.67×10–3).

Conclusions

The authors found evidence for relative prenatal enrichment of putative susceptibility genes for syndromic neurodevelopmental disorders, intellectual disability, and autism spectrum disorder. Future transcriptome-level association studies should evaluate regions other than the dorsolateral prefrontal cortex, at other time points, and incorporate further RNA sequencing analyses.
Adult-onset psychiatric disorders are increasingly being viewed as of neurodevelopment origin, at least in part. For example, neurodevelopmental hypotheses of schizophrenia introduced and articulated in the 1980s posited that a pathological perturbation occurs early in life and that after a latent period, as the brain matures in early adulthood and calls into action the compromised neural systems, and concurrent with environmental complexity, the perturbation becomes clinically manifest, resulting in the symptoms that constitute the syndrome of adult-onset schizophrenia (1).
The neurodevelopmental model of schizophrenia has gained considerable traction and evidentiary support since first proposed. Different versions of the neurodevelopmental model have argued that different temporal periods of brain development are most critical in pathogenesis: e.g., prenatal or postnatal events (2) or adolescent synaptic refinement (3). Epidemiological evidence for prenatal insults associated with increased schizophrenia frequency includes famine exposure (4) and maternal infection (5). Minor physical anomalies in the craniofacial region and limbs observed in patients with schizophrenia suggest disturbed prenatal ectoderm development (6). Evidence for perinatal insults includes the association of obstetric complications with schizophrenia, which leads to up to a twofold increase in the frequency of schizophrenia, with the magnitude of the effect varying by the type of obstetric complication (7, 8), and with schizophrenia risk-associated genes potentially interacting with obstetric complications (9, 10). In children who subsequently develop schizophrenia, motor abnormalities (11), social abnormalities (12), and cognitive deficits (13, 14) are present long before the onset of psychotic symptomatology. In animal models, prenatal and neonatal disruptions of cortical development lead to the emergence of abnormalities in a number of cognitive and dopamine-related behaviors and glutamatergic-gabaergic pharmacological abnormalities that do not emerge until early adulthood (1517). At the cellular level, altered postnatal maturation processes hypothesized to play a role in pathogenesis include GABA interneuron maturation, pruning of glutamate synapses, maturation of dopaminergic projections, and oligodendrocyte differentiation and myelination (18).
Yet, despite the accumulating circumstantial evidence for a neurodevelopmental pathogenesis, a definitive cellular or molecular explanation of the neurobiology of an early pathogenic process has remained elusive. Emerging molecular evidence suggests that analyzing temporal expression patterns of genes associated with schizophrenia may provide insight into neurodevelopmental pathogenesis; for several genes associated with schizophrenia, alternatively spliced transcripts that are expressed relatively abundantly in the developing fetal brain may in fact contribute differentially to the mechanism of risk that underlies the clinical association (19). Alternatively spliced transcripts that are expressed prenatally and associated with schizophrenia allelic risk variants have been demonstrated for multiple schizophrenia candidate genes to date, including DISC1, neuregulin 1, neuregulin 3, GAD1, ZNF804A, and KCNH2 (2025). Furthermore, some groups reported that genes associated with schizophrenia by exome sequencing of de novo rare variants are prenatally enriched in transcript abundance (26, 27).
More broadly, beyond schizophrenia, other psychiatric disorders have been implicated as having a developmental etiology. Obstetric complications may be an antecedent, or risk factor, for psychiatric disorders other than schizophrenia (28). ASD and intellectual disability are readily classified as neurodevelopmental disorders in the clinical and research literature (29), and their clinical onset in infancy or early childhood especially implicates an early developmental abnormality.
An intriguing hypothesis then is that a neuropsychiatric disorder considered to be developmental in etiology is associated with genes that are preferentially fetally expressed, that is, with genes that are more transcriptionally active or abundant (i.e., differentially regulated) in the prenatal period compared with the postnatal period. Furthermore, the preferentially fetally expressed genes are likely to be enriched for functions that have been associated in previous gene expression studies with the developing fetal brain: neural development and differentiation, axon guidance, transcription factor, and gene regulation activity (30, 31). This would implicate factors related to gene regulation in the fetal developing brain as critical for the genetic association with a developmental phenotype.
The present study investigated risk genes systematically, as composite gene sets associated with a specific neuropsychiatric disorder, for the presence of transcripts that are preferentially fetally expressed in the transcriptome of the dorsolateral prefrontal cortex. We report here the prenatal expression patterns of susceptibility genes associated with four neuropsychiatric disorders that are thought to have developmental components to their pathogenesis: syndromic neurodevelopmental disorders, autism spectrum disorder (ASD), intellectual disability, and schizophrenia, and one group of disorders thought not to be principally of developmental origin: neurodegenerative disorders. We also include genes associated with bipolar affective disorder in the present analysis.

Method

Gene Set Assignment

Based on literature review, susceptibility genes were classified into six neuropsychiatric disorders and 12 gene sets based on differing genetic evidence (see data supplement that accompanies the online edition of this article for methodology of gene set assignment). The gene sets were then explored in the public database of human prefrontal cortex gene expression across the lifespan, described in BrainCloud (32)(http://braincloud.jhmi.edu/ and data accessible at NCBI GEO database accession GSE30272). As described by Colantuoni et al. (32), RNA from 269 human postmortem dorsolateral prefrontal cortex samples ranging from fetal life (N=38 fetal samples) up to age 80 years, was extracted and analyzed on custom 2-color microarrays using a reference RNA comprised of a pool of all samples in each array, with hybridization at 30,176 probes reported.

Statistical Methods

Linear regression was performed at each of 30,176 probes, modeling gene expression as a continuous dependent variable, as a function of a binary “fetal” stage variable, including surrogate variables as adjustment variables to account for potential batch effects and other unwanted technical and biological variation.
More formally, the statistical model was
for probe i∈[1,…,30,176] and subject j∈[1,…,269], with indicator variable Xjfetal∈[0,1] and surrogate variable set SV(SV1,…,SV30). This regression model generated a “fetal effect” (βi), the log2-fold change of prenatal versus postnatal transcript abundance, at each probe on the microarray.
Each gene set was then evaluated by a Wilcoxon signed rank test to determine if the fetal effect distribution of the set differed significantly from the entire dorsolateral prefrontal cortex transcriptome, using the reduced statistic of one fetal effect per gene.
The significance of genes within each gene set with at least one transcript that was preferentially fetally expressed was assessed using a binomial proportion test, compared with background proportions in the genome, at increasing stringency of fetal effects on the log2 scale, 0.5, 1.0, and 1.5 (that is, a 1.4-fold, twofold, and 2.8-fold increase on a non-log2 scale,).

Results

Fetal Effect Analysis

As an example of the “fetal effect” metric, Figure 1 depicts the neurodevelopmental expression patterns of two genes across the lifespan, ANK3 and HTR2A. ANK3 is preferentially fetally expressed while HTR2A is relatively fetally underexpressed, with fetal effects of 2.42 and −2.85 respectively. Across the transcriptome, of 30,176 probes that mapped to 17,677 genes, 43% of genes expressed a transcript with a fetal effect greater than 0 (and p<0.05). Twenty-five percent of genes expressed a transcript with a fetal effect greater than 0.5, and 12% of genes expressed a transcript with a fetal effect greater than 1. All probes with a fetal postnatal log2-fold change greater than 0.5 had a p value less than 0.05 (except for two probes). See Table S1 in the online data supplement for lists of each gene in a gene set and its “fetal effect,” that is, the log2-fold change, or the calculated regression coefficient.
FIGURE 1. Fetal Effect by Individual Gene Calculated for Each Gene Expressing a Transcript in the Dorsolateral Prefrontal Cortexa
a Two examples are illustrated. The fetal effect is prenatal compared with postnatal log2-fold change or regression coefficient. N=269 control samples (including 38 fetal samples). The x axis represents time course, with fetal stage in weeks and postnatal stage in years. Plots above are obtained from BrainCloud application (http://braincloud.jhmi.edu/).
The binomial proportion tests (Table 1) reveal statistically significant preferential fetal expression for the following disease-associated gene groups: syndromic neurodevelopmental disorders (N=31 genes, 52% of genes with a fetal effect greater than 0.5, p=3.37×10−3); intellectual disability (N=88 genes, 40% of genes with a fetal effect greater than 0.5, p=5.53×10−3); and ASD (N=242 genes, 36% of genes with a fetal effect greater than 0.5, p=3.45×10−4). In contrast, as illustrated in the density plots and Wilcoxon tests in Figure 2, we found statistically significant fetal underexpression for neurodegenerative disorders (N=46 genes, one-sided p value for Wilcoxon signed-rank test in the negative direction was 9.96×10−4). The dissociation between relative fetal overexpression in genes associated with intellectual disability, syndromic neurodevelopmental disorders, and ASD and the relative fetal underexpression of genes associated with late-life neurodegenerative disorders provides conceptual support for the present approach.
TABLE 1. Fetal Effect by Gene Set and Binomial Testa
Gene SetN of GenesN of Genes Mapped to a Probe% of Genes in Gene Set with Fetal Effect >0p% of Genes in Gene Set with Fetal Effect >0.5p (Two-Sided)% of Genes in Gene Set with Fetal Effect >1p% of Genes in Gene Set with Fetal Effect >1.5p
Genome-wide 17,67843 25 12 5 
Neurodevelopmental disorders3331651.13E−02523.37E−03231.50E−01169.64E−02
Schizophrenia (SRF) Meta-analysis4036428.95E−01226.57E−0184.57E−0133.87E−01
Schizophrenia (PGC) GWAS132106501.13E−01294.58E−01153.21E−0168.02E−01
Schizophrenia CNV154113416.58E−01269.29E−01116.94E−0147.09E−01
Schizophrenia SNV/INDEL240212481.54E−01292.06E−01151.75E−0167.47E−01
ASD CNV216176415.57E−01235.60E−0192.30E−0132.11E−01
ASD database262242494.58E−02363.45E−04194.06E−03112.98E−03
Bipolar disorder GWAS146123383.00E−01248.15E−01118.95E−0167.95E−01
Intellectual disability9488561.46E−02405.53E−03222.50E−02133.70E−02
Neurodegenerative disorders5946282.90E−02111.67E−0341.37E−020 
Overlap across all categories 62457.03E−01313.61E−01163.50E−01111.26E−01
Schizophrenia/ASD CNV overlap 73456.74E−01268.87E−0182.71E−0146.55E−01
a
The percentage of genes within each gene set that are preferentially fetally expressed by increasing stringency of fetal effect: 0.5, 1.0, and 1.5. Two-sided p values for binomial proportion tests are listed. Shading indicates significant elevation compared to the genome and bolding indicates significant decrease compared to the genome [see article PDF for shading]. ASD=autism spectrum disorder; CNV=copy number variant; SNV=single nucleotide variant; INDEL=short insertion/deletion; GWAS=genome-wide association study; SRF=Schizophrenia Research Forum; PGC=Psychiatric Genomics Consortium.
FIGURE 2. Density Plot of Fetal Effect by Gene Set and Wilcoxon Signed-Rank Testa
a The figures depict the density of each gene set’s fetal effect compared to the fetal effect of the whole genome. The x axis is the fetal effect, or regression coefficient (prenatal compared with postnatal) and the y axis is the density. The black line indicates the genome and each colored line indicates the gene set. The red line indicates significant increase, the purple line indicates significant decrease, and the blue line (default) indicates no significant difference in fetal effect of gene set compared to the transcriptome. One-sided p values are shown for Wilcoxon signed-rank test for each gene set, elevation in the positive direction on the upper right side, and decrease in the negative direction on the upper left side.
Genes associated with schizophrenia by meta-analysis were not significantly different than the genome in relative overall fetal expression (N=36 genes, 22% of genes with a fetal effect greater than 0.5, p=0.66) nor were the genes reported in a series of genome-wide association study (GWAS) reports by Ripke et al. in 2011 (33) and 2013 (34) (N=106 genes, 29% of genes with a fetal effect greater than 0.5, p=0.46). Genes associated with schizophrenia by copy number variant (CNV) analysis were also not significantly different compared with the genome in relative fetal expression (N=113 genes, one-sided p value for Wilcoxon signed-rank test in the positive direction was 0.80) nor were genes associated with schizophrenia by exome sequencing and rare de novo single nucleotide variants or short insertions/deletions (N=212 genes, one-sided p value for Wilcoxon signed rank test in the positive direction was 0.29).
Genes within loci associated with ASD by CNV analysis did not show any statistically significant difference, neither elevation nor decrease in relative fetal expression, compared with the genome (N=176 genes, one-sided p value for Wilcoxon signed-rank test in the positive direction was 0.91 and 0.09 in the negative direction). The gene set associated with bipolar disorder by GWAS was also not significantly different than the genome in prenatal expression abundance (N=123 genes, 24% with a fetal effect greater than 0.5, p=0.82). Furthermore, the set of genes that overlapped diagnostic categories (listed in Table S2 in the online data supplement) was not significantly different than the genome in prenatal expression abundance (N=62 genes, 31% with a fetal effect greater than 0.5, p=0.36).

Functional Analysis of Prenatally Enriched Genes

Gene set enrichment analysis was performed for the subset of prenatally enriched genes for each gene set (defined as a gene with a fetal effect greater than 0.5), compared with the background universe of all brain-expressed genes (see Table S3 in the online data supplement for a list of significant gene ontology [GO] terms). The gene sets that demonstrated preferential fetal expression patterns were enriched for multiple functions pertinent to early cellular development and neurodevelopment. For example, genes associated with syndromic neurodevelopmental disorders were broadly enriched for “generation of neurons” (GO:0048699) and for “system development” (GO:0048731) and genes associated with intellectual disability were enriched for “central nervous system development” (GO:0007417) and “cell projection organization” (GO:0030030).
Genes from the ASD database demonstrated the most significant enrichment of all genes sets for multiple neurodevelopmental processes, including “neurogenesis” (GO:0022008), “axogenesis” (GO:0007409), and “neuron development” (GO:0048666). The group of genes associated with schizophrenia by meta-analysis were enriched for multiple developmental processes including “nervous system development” (GO:0007399) and “neuron differentiation” (GO:0048468).
Further functional investigation of the prenatally enriched genes consisted of a coexpression analysis, determining which genes are coexpressed during the second trimester with a disease-associated and prenatally enriched gene. As reported in Table S4 in the online data supplement, while some genes revealed isolated patterns of expression, others appeared to be coexpressed with numerous genes in the second trimester, indicating a potential functional concordance or coregulation.
To determine potential interactions among prenatally enriched genes, Figure S1 in the online data supplement illustrates the “top scoring” networks from the ingenuity pathway analysis for genes prenatally enriched in the ASD gene set and for the genes prenatally enriched in the combined schizophrenia data sets (GWAS, CNV association, and single nucleotide variants) along with selected functional annotation of individual genes. The networks yielded some noteworthy potential interactions and candidates, for example the 26s Proteasome complex in ASD and the NF-kB complex in schizophrenia.

Secondary Analysis of CNV Regions

We performed a secondary post hoc analysis on genes within CNV regions associated with schizophrenia and ASD, under the assumption that these broad regions involve multiple genes, some of which may demonstrate a pattern of preferential fetal expression and therefore be more pathogenically relevant (Tables 2, 3, 4 and see Table S5 in the online data supplement). Within the loci affected by CNVs associated with schizophrenia and ASD, select genes were indeed preferentially fetally expressed, with one or two genes at each locus demonstrating a comparatively greater fetal effect.
TABLE 2. Genes Within Copy Number Variant (CNV) Regions Associated With Schizophrenia and Autism Spectrum Disorder (ASD)a
CNV Loci Associated With ASD and Schizophrenia
LocusCNV TypeGeneFetal Effectp
1q21.1Dup/DelNBPF13Pn/a 
  PRKAB21.398.97E−57
  FMO5  
  CHD1Ln/a 
  BCL92.022.50E−63
  ACP6n/a 
  GJA50.252.82E−03
  GJA8−0.491.40E−09
  GPR89Bn/a 
  GPR89Cn/a 
  PDZK1P1n/a 
  NBPF8n/a 
  NBPF24n/a 
  NBPF11n/a 
  NBPF10n/a 
2p16.3Dup/DelNRXN11.316.02E−34
15q13.3Del/DupCHRFAM7A−1.545.29E−32
  ARHGAP11Bn/a 
  MTMR150.501.10E−12
  MTMR10n/a 
  TRPM1n/a 
  MIR211n/a 
  KLF13−0.532.04E−09
  OTUD7A−0.966.01E−21
 
 
CHRNA7
0.83
4.76E−08
16p11.2Del/DupSPN0.157.88E−04
  QPRT−0.863.41E−22
  C16orf54−0.019.09E−01
  ZG16n/a 
  KIF220.224.67E−04
  MAZ−0.067.63E−02
  PRRT2−0.573.02E−13
  PAGR1 (C16orf53)1.291.49E−51
  MVP−1.163.84E−28
  CDIPT−0.599.04E−32
  SEZ6L2−1.441.83E−42
  ASPHD10.741.69E−15
  KCTD130.161.68E−02
  TMEM219−0.951.04E−47
  TAOK20.507.67E−07
  HIRIP30.162.33E−02
  INO80E−0.034.98E−01
  DOC2A−0.643.67E−30
  C16orf92n/a 
  FAM57B−0.643.67E−30
  ALDOA0.321.08E−06
  PPP4C1.291.69E−49
  TBX6−0.781.59E−18
  YPEL3−0.951.35E−39
  GDPD30.242.03E−05
  MAPK3−1.222.91E−55
 
 
CORO1A
−0.97
5.34E−40
22q11.21Del/DupDGCR6−0.732.53E−31
  PRODH−1.016.36E−17
  DGCR50.654.41E−13
  DGCR9-0.603.13E−10
  DGCR100.563.37E−05
  DGCR2−0.732.09E−28
  DGCR110.009.72E−01
  DGCR140.451.60E−07
  TSSK2−0.559.14E−08
  GSC2n/a 
  SLC25A10.482.19E−15
  CLTCL10.482.19E−15
  HIRA0.589.08E−24
  MRPL40−0.501.74E−13
  C22orf39n/a 
  UFD1L0.255.60E−07
  CDC45L3.891.59E−68
  CLDN5n/a 
  SEPT5  
  GP1BBn/a 
  TBX1−0.706.57E−25
  GNB1L-0.251.14E−04
  C22orf290.071.70E−01
  TXNRD2−0.849.34E−41
  COMT−0.282.66E−07
  ARVCF0.574.53E−21
  C22orf25−1.478.40E−63
  MIR185n/a 
  DGCR80.821.92E−09
  MIR3618n/a 
  MIR1306n/a 
  TRMT2A−0.401.14E−08
  RANBP10.404.51E−16
  ZDHHC8−0.164.09E−04
  RTN4R−0.815.47E−18
  MIR1286n/a 
  DGCR6L−0.744.76E−32
  TMEM191Bn/a 
  RIMBP3n/a 
  ZNF741.572.05E−62
  SCARF2−0.172.92E−03
  KLHL220.385.11E−16
  MED15−0.491.36E−21
  PI4KA0.502.81E−08
  SERPIND10.742.23E−09
  SNAP29−0.696.90E−19
  CRKL0.053.49E−01
  AIFM3n/a 
  LZTR10.561.98E−22
  THAP7−0.137.74E−04
  P2RX60.162.00E−03
  SLC7A4−0.831.97E−26
  HIC22.321.97E−26
a
Genes with a fetal effect greater than 0.5 are in bold and genes with a fetal effect greater than 1 are shaded [see article PDF for shading]; Genes with no microarray probes are indicated as n/a.
TABLE 3. Other Schizophrenia-Associated CNV Locia
LocusCNV TypeGenesFetal Effectp
1q21.1DelHFE20.094.82E−01
  TXNIP0.862.33E−09
  POLR3GL0.052.97E−01
  ANKRD34A−0.985.24E−33
  LIX1L1.386.15E−59
  RBM8A0.902.48E−37
  GNRHR2−0.493.85E−07
  PEX11B−0.701.44E−36
  ITGA10n/a 
  ANKRD35−0.612.04E−14
  PIAS30.071.46E−01
  NUDT17−0.303.23E−08
  POLR3C0.196.81E−05
  RNF1150.082.06E−01
  CD160−0.491.18E−15
  PDZK10.233.52E−02
  GPR89An/a 
a
Genes with a fetal effect greater than 0.5 are in bold and genes with a fetal effect greater than 1 are shaded [see article PDF for shading]; Genes with no microarray probes are indicated as n/a. CNV=copy number variant.
TABLE 4. Other Autism Spectrum Disorder-Associated CNV Locia
LocusCNV TypeGeneFetal Effectp
3p14.1DelSUCLG2−0.601.49E−16
  FAM19A1n/a 
  FAM19A4n/a 
  C3orf64  
  MIR3136n/a 
  TMF10.241.03E−01
  UBA31.382.01E−56
  ARL6IP5−1.024.03E−40
  LMOD3-0.356.50E−03
  FRMD4B1.781.97E−34
 
 
MITF
0.17
3.79E−02
7q11.23DupFZD9−0.413.34E−05
  BAZ1B0.159.52E−03
  BCL7B0.211.84E−03
  TBL20.179.66E−03
  MLXIPLn/a 
  VPS37D0.818.78E−26
  DNAJC300.205.41E−04
  WBSCR22−0.372.77E−14
  STX1A−1.525.20E−49
  MIR4284n/a 
  ABHD110.046.41E−01
  CLDN3−0.273.53E−06
  CLDN40.028.86E−01
  WBSCR27−0.292.74E−02
  WBSCR280.191.21E−02
  ELN−0.037.20E−01
  LIMK1−0.502.58E−10
  MIR590n/a 
  LAT2−1.225.75E−15
  RFC20.833.10E−39
  CLIP20.187.63E−03
  GTF2IRD11.681.64E−68
 
 
GTF2I
n/a
 
7q31.1
Del/Dup
DOCK4
0.12
1.67E−01
8p23.3
Dup
DLGAP2
−2.22
1.96E−66
9p24.3Del/DupFOXD4n/a 
  CBWD10.522.04E−13
  C9orf66−0.382.01E−07
 
 
DOCK8
−1.52
4.19E−18
10q11.23-21.1DupPRKG1n/a 
  MIR605n/a 
  CSTF2T0.391.29E−06
  DKK10.019.12E−01
 
 
MBL2
0.23
2.37E−02
12q24.31DupTMEM120Bn/a 
  RHOF−1.581.67E−63
  SETD1B1.306.47E−58
  HPD−0.129.20E−02
  PSMD90.181.09E−02
  WDR66−0.585.32E−10
a
Genes with a fetal effect greater than 0.5 are in bold and genes with a fetal effect greater than 1 are shaded [see article PDF for shading]; Genes with no microarray probes are indicated as n/a. CNV=copy number variant.
At the 1q21.1 locus, PRKAB2 (fetal effect=1.39) and BCL9 (fetal effect=2.02) were most abundant in prenatal compared with postnatal transcript abundance. The single gene deletion at 2p16.3, NRXN1, also common to both schizophrenia and autism, demonstrated a fetal effect of 1.31. At the 15q13.3 locus, CHRNA7 demonstrated a fetal effect greater than other genes at 0.83, followed by MTMR15 (fetal effect=0.50). At the 16p11.2 region, PAGR1/C16orf53 (fetal effect=1.29) and PPP4C (fetal effect=1.29) were most relatively fetally abundant. Lastly at the 22q11.21 region, CDC45L (fetal effect=3.89) and HIC2 (fetal effect=2.32) were most notable in relative fetal abundance, along with ZNF74 (fetal effect=1.70). See Table S5 in the online data supplement for a list of the genes within each CNV locus for schizophrenia and ASD, respectively, to enable comparison of genes with a given disease-associated CNV locus.

Secondary Analysis of Genes in the ASD Database

Since the genes in the SFARI ASD database were of differing genetic evidence, functional biology (N=29 genes), genetic association (N=110 genes), and single rare variants (N=103 genes), we conducted a secondary analysis to determine if one subgroup was driving the fetal effect. Indeed, while the former two were not significantly different than the genome in fetal effect, the set of genes with single rare variants was highly significantly different than the overall genome (N=103 genes, one-sided p value for Wilcoxon signed rank test in the positive direction was 4.75×10–7) with 15% of the genes with a fetal effect greater than 1.5 (p=0.007)
As most of our candidate genes were sequence variants based on clinical association studies, we examined an additional ASD gene set from the study by Voineagu et al. (35) that identified differentially expressed genes between adult autism case and control subjects and also identified weighted gene coexpression networks associated with ASD. We found that the M12 module, reported to be underexpressed in autism samples, and whose eigengene was correlated with disease status, was highly significant for underexpression in fetal samples compared with the fetal expression of the overall genome (N=273 genes, p=1.59×10–34). It is important to note that only 18 of the genes in the reported M12 module overlap our combined list of 418 ASD candidate genes, and the ASD etiological candidacy of the remaining genes in the M12 module remains to be determined.

Secondary Analysis of Genes Associated With Schizophrenia by Exome Sequencing of Rare Single Nucleotide Variants

As a result of possibly discrepant findings between previously published reports of prenatal enrichment of transcripts of genes associated with schizophrenia by exome sequencing and the current analysis, we further analyzed the subset of genes recently reported by Gulsuner et al. (27). While Gulsuner et al. report that genes associated with schizophrenia by exome sequencing of rare variants are expressed during the fetal period, relative and differentially increased prenatal expression was not described in that report, which is the scope of the present inquiry. Thus, we used the publicly available BrainSpan RNA sequencing data (36) to evaluate the dorsolateral prefrontal cortex expression of 52 of the genes reported by Gulsuner et al. in fetal samples (N=17) and postnatal samples (N=18) and found that only 21 of the 52 genes had a higher fetal mean of reads per kilobase of transcript per million (RPKM) expression, that is, 40% were preferentially fetally expressed, corroborating our microarray data of a lack of a preferential fetal expression pattern for this subset of genes.

Discussion

In the present analysis, we calculated a fetal effect for each gene expressed in the dorsolateral prefrontal cortex transcriptome by comparing gene expression in the second trimester of fetal development to postnatal expression and then by comparing the disease-associated gene sets to the genome at large. The principal question we asked was whether genes associated with putatively neurodevelopmental behavioral disorders would show a relatively positive fetal effect, i.e., preferential fetal expression compared with postnatal expression.
Our finding of prenatal enrichment of expression of genes associated with syndromic neurodevelopmental disorders and intellectual disability may be interpreted as a positive control. It is probable that syndromic neurodevelopmental disorders and intellectual disability have a fetal pathogenic causality; syndromic neurodevelopmental disorders may be apparent at birth, and intellectual deficiencies are often apparent by infancy or early childhood. In contrast, the counterpoint finding of significant underexpression of genes associated with neurodegenerative disorders is consistent with an onset later in life and with increased disease-associated gene expression preceding or concurrent with gliosis, apoptosis, and other cellular processes characteristic of neurodegeneration. Thus, notwithstanding the limitations of our approach, which include the imprecision and low biological resolution of surveys of gene expression in the postmortem brain based on oligonucleotide arrays (vide infra), our positive results at the two extremes of the spectrum of neurodevelopmental disorders add credibility to the present approach.

Prenatal Expression of ASD-Associated Genes

Genes with rare single nucleotide variants that were associated with ASD demonstrated a strong pattern of preferential fetal expression. This finding, however, should be viewed in the context of the phenotypic heterogeneity of the reported ASD cases, with different studies using different criteria for ASD, some studies including syndromic forms of ASD or intellectual deficits of varying magnitude. Nevertheless, for the current putative ASD candidate gene list, there is an overall pattern of preferential fetal expression compared with the genome, as might be expected for genes contributing to the pathogenesis of a clinical phenotype present by early childhood.

Prenatal Expression of Schizophrenia-Associated Genes

The present findings do not strongly support a specific fetal pathogenic causality for the composite set of genes associated with schizophrenia, by GWAS, CNV analysis, or exome sequencing of rare variants, at least not a causality that would be reflected in differential fetal expression. We believe, however, that several caveats are important to consider in interpreting these results.
For genes implicated in genome-wide association studies, multiple genes in the linkage disequilibrium region of an associated single-nucleotide polymorphism were selected for analysis. The probable overinclusion of genes in disease-associated loci (but not true risk genes) may have decreased the probability of observing an overall fetal expression pattern. Within each CNV locus associated with schizophrenia (or ASD), one or two genes were enriched in prenatal abundance, suggesting that they may be particularly relevant to a potential fetal pathogenesis. For example, at the 22q11.2 region, CDC45L showed an especially high relative fetal abundance, the second highest of all of the genes in the present analysis. It encodes the CDC45 protein, which forms a complex required for the activation of DNA helicase and essential for DNA replication (37). At 16p11.2, PAGR1(C16orf53) was preferentially fetally expressed, and recent evidence suggests that this gene modulates glucocorticoid-receptor gene induction (38).
Overall, the gene-set enrichment analyses confirmed the critical functionality of genes that are preferentially fetally expressed and associated with syndromic neurodevelopmental disorders, intellectual disability, and ASD. The results also suggest that genes associated with schizophrenia by meta-analysis are enriched for neurodevelopmental processes, potentially supporting a fetal pathogenic mechanism for schizophrenia.

Future Transcriptome Investigations

As the present evaluation is a quantitative assessment of genes that are preferentially fetally expressed, further investigation should also consider qualitative evaluation of individual specific genes, potential splicing variation and fetal transcriptional activity, and the impact of specific genes on cellular processes and consequently neuropsychiatric disorders. Within each gene set, there may be a qualitative distinction in fetal expression, e.g., splicing or 3′ untranslated region activity, that overrides the significance of an overall quantitative measure.
It is important to consider neurodevelopmental pathogenic models beyond simply a fetal pathogenic effect. For example, a disease-associated gene may exert its pathogenic effect at multiple time points. Alternatively, the disorder may be heterogeneous with different phenotypic subtypes of differing etiology, so that one subtype may have a fetal pathogenic mechanism, while another subtype may have a different pathogenic mechanism. Furthermore, genes may be pathogenic factors in disorders like schizophrenia, but not necessarily be differentially upregulated—they may be repressed in programmed development, act at low levels in gradient patterning, or act by another mode of regulatory action in fetal development. Also, a pathogenic role for genes in fetal development may be subserved simply by their being expressed at any level in fetal developmental.
The present analysis compared prenatal to postnatal expression but did not further explore the timing of specific peaks and troughs or specific postnatal trends. Previous transcriptome-level studies of fetal postmortem samples, including those by Johnson et al. (30) and Kang et al. (31) have demonstrated that regional transcriptomes are developmentally regulated. Thus, ongoing refined examination during particular developmental time windows and at regions other than the dorsolateral prefrontal cortex may lead to further functional and ontological insight for an individual gene, gene group, or cluster of coexpressed genes.
From a methodological perspective, our postmortem analysis is potentially confounded by shifting cellular heterogeneity (postnatal samples are increasingly populated by glial cells, while prenatal samples are predominantly neuronal) and RNA integrity (higher quality in the prenatal than postnatal samples). While the issue of cellular composition warrants further investigation, the present analysis controlled for RNA integrity number (RIN) by considering “negative control genes” associated with RIN. Lastly, the microarray platform is limited by an a priori probe design, which may not include the presence of pathogenic and alternatively spliced fetal transcripts; future RNA sequencing analysis will be advantageous in permitting the investigation of novel, as yet uncharacterized transcripts. Thus, the results of the present study, while provocative, must be viewed as preliminary in exploration of the developmental aspects of psychiatric genetics. As the polygenic and developmental etiology of schizophrenia and other neuropsychiatric disorders becomes increasingly substantiated, and as the hope for gene-based diagnostics and therapeutics grows, a fetal component to pathogenesis is critical to consider, since well prior to the onset of apparent symptoms, in utero disease-associated gene expression may have prophylactic and therapeutic implications.

Footnote

Dr. Birnbaum receives support from an NIH T32, 5T32MH015330-36

Supplementary Material

Supplementary Material (758_ds001.pdf)

References

1.
Weinberger DR: Implications of normal brain development for the pathogenesis of schizophrenia. Arch Gen Psychiatry 1987; 44:660–669
2.
Murray RM, Lewis SW: Is schizophrenia a neurodevelopmental disorder? Br Med J (Clin Res Ed) 1987; 295:681–682
3.
Feinberg I: Schizophrenia: caused by a fault in programmed synaptic elimination during adolescence? J Psychiatr Res 1982-1983-1983; 17:319–334
4.
Susser E, Neugebauer R, Hoek HW, Brown AS, Lin S, Labovitz D, Gorman JM: Schizophrenia after prenatal famine: further evidence. Arch Gen Psychiatry 1996; 53:25–31
5.
Brown AS: Exposure to prenatal infection and risk of schizophrenia. Front Psychiatry 2011; 2:63
6.
Weinberg SM, Jenkins EA, Marazita ML, Maher BS: Minor physical anomalies in schizophrenia: a meta-analysis. Schizophr Res 2007; 89:72–85
7.
Lewis SW, Murray RM: Obstetric complications, neurodevelopmental deviance, and risk of schizophrenia. J Psychiatr Res 1987; 21:413–421
8.
Geddes JR, Lawrie SM: Obstetric complications and schizophrenia: a meta-analysis. Br J Psychiatry 1995; 167:786–793
9.
Cannon TD, Yolken R, Buka S, Torrey EF; Collaborative Study Group on the Perinatal Origins of Severe Psychiatric Disorders: Decreased neurotrophic response to birth hypoxia in the etiology of schizophrenia. Biol Psychiatry 2008; 64:797–802
10.
Nicodemus KK, Marenco S, Batten AJ, Vakkalanka R, Egan MF, Straub RE, Weinberger DR: Serious obstetric complications interact with hypoxia-regulated/vascular-expression genes to influence schizophrenia risk. Mol Psychiatry 2008; 13:873–877
11.
Walker EF: Developmentally moderated expressions of the neuropathology underlying schizophrenia. Schizophr Bull 1994; 20:453–480
12.
Done DJ, Crow TJ, Johnstone EC, Sacker A: Childhood antecedents of schizophrenia and affective illness: social adjustment at ages 7 and 11. BMJ 1994; 309:699–703
13.
Davidson M, Reichenberg A, Rabinowitz J, Weiser M, Kaplan Z, Mark M: Behavioral and intellectual markers for schizophrenia in apparently healthy male adolescents. Am J Psychiatry 1999; 156:1328–1335
14.
Woodberry KA, Giuliano AJ, Seidman LJ: Premorbid IQ in schizophrenia: a meta-analytic review. Am J Psychiatry 2008; 165:579–587
15.
Lipska BK, Weinberger DR: Delayed effects of neonatal hippocampal damage on haloperidol-induced catalepsy and apomorphine-induced stereotypic behaviors in the rat. Brain Res Dev Brain Res 1993; 75:213–222
16.
Lodge DJ, Grace AA: Gestational methylazoxymethanol acetate administration: a developmental disruption model of schizophrenia. Behav Brain Res 2009b; 204:306–312
17.
Lavin A, Moore HM, Grace AA: Prenatal disruption of neocortical development alters prefrontal cortical neuron responses to dopamine in adult rats. Neuropsychopharmacology 2005; 30:1426–1435
18.
Jaaro-Peled H, Hayashi-Takagi A, Seshadri S, Kamiya A, Brandon NJ, Sawa A: Neurodevelopmental mechanisms of schizophrenia: understanding disturbed postnatal brain maturation through neuregulin-1-ErbB4 and DISC1. Trends Neurosci 2009; 32:485–495
19.
Kleinman JE, Law AJ, Lipska BK, Hyde TM, Ellis JK, Harrison PJ, Weinberger DR: Genetic neuropathology of schizophrenia: new approaches to an old question and new uses for postmortem human brains. Biol Psychiatry 2011; 69:140–145
20.
Nakata K, Lipska BK, Hyde TM, Ye T, Newburn EN, Morita Y, Vakkalanka R, Barenboim M, Sei Y, Weinberger DR, Kleinman JE: DISC1 splice variants are upregulated in schizophrenia and associated with risk polymorphisms. Proc Natl Acad Sci USA 2009; 106:15873–15878
21.
Tan W, Wang Y, Gold B, Chen J, Dean M, Harrison PJ, Weinberger DR, Law AJ: Molecular cloning of a brain-specific, developmentally regulated neuregulin 1 (NRG1) isoform and identification of a functional promoter variant associated with schizophrenia. J Biol Chem 2007; 282:24343–24351
22.
Kao WT, Wang Y, Kleinman JE, Lipska BK, Hyde TM, Weinberger DR, Law AJ: Common genetic variation in Neuregulin 3 (NRG3) influences risk for schizophrenia and impacts NRG3 expression in human brain. Proc Natl Acad Sci USA 2010; 107:15619–15624
23.
Hyde TM, Lipska BK, Ali T, Mathew SV, Law AJ, Metitiri OE, Straub RE, Ye T, Colantuoni C, Herman MM, Bigelow LB, Weinberger DR, Kleinman JE: Expression of GABA signaling molecules KCC2, NKCC1, and GAD1 in cortical development and schizophrenia. J Neurosci 2011; 31:11088–11095
24.
Tao R, Cousijn H, Jaffe AE, Burnet P, Philip WJ, Edwards F, Eastwood S, Shin JH, Lane T, Walker MA, Maher BJ, Weinberger DR, Harrison P, Hyde TM, Kleinman JE: A novel transcript fetally regulated by the psychosis risk SNP rs1344706, and alterations in schizophrenia, bipolar disorder, and major depression. JAMA Psychiatry (in press)
25.
Huffaker SJ, Chen J, Nicodemus KK, Sambataro F, Yang F, Mattay V, Lipska BK, Hyde TM, Song J, Rujescu D, Giegling I, Mayilyan K, Proust MJ, Soghoyan A, Caforio G, Callicott JH, Bertolino A, Meyer-Lindenberg A, Chang J, Ji Y, Egan MF, Goldberg TE, Kleinman JE, Lu B, Weinberger DR: A primate-specific, brain isoform of KCNH2 affects cortical physiology, cognition, neuronal repolarization and risk of schizophrenia. Nat Med 2009; 15:509–518
26.
Xu B, Ionita-Laza I, Roos JL, Boone B, Woodrick S, Sun Y, Levy S, Gogos JA, Karayiorgou M: De novo gene mutations highlight patterns of genetic and neural complexity in schizophrenia. Nat Genet 2012; 44:1365–1369
27.
Gulsuner S, Walsh T, Watts AC, Lee MK, Thornton AM, Casadei S, Rippey C, Shahin H, Nimgaonkar VL, Go RC, Savage RM, Swerdlow NR, Gur RE, Braff DL, King MC, McClellan JM; Consortium on the Genetics of Schizophrenia (COGS); PAARTNERS Study Group: Spatial and temporal mapping of de novo mutations in schizophrenia to a fetal prefrontal cortical network. Cell 2013; 154:518–529
28.
Buka SL, Fan AP: Association of prenatal and perinatal complications with subsequent bipolar disorder and schizophrenia. Schizophr Res 1999; 39:113–119, discussion 160–161
29.
Diagnostic and Statistical Manual of Mental Disorders
30.
Johnson MB, Kawasawa YI, Mason CE, Krsnik Z, Coppola G, Bogdanović D, Geschwind DH, Mane SM, State MW, Sestan N: Functional and evolutionary insights into human brain development through global transcriptome analysis. Neuron 2009; 62:494–509
31.
Kang HJ, Kawasawa YI, Cheng F, Zhu Y, Xu X, Li M, Sousa AM, Pletikos M, Meyer KA, Sedmak G, Guennel T, Shin Y, Johnson MB, Krsnik Z, Mayer S, Fertuzinhos S, Umlauf S, Lisgo SN, Vortmeyer A, Weinberger DR, Mane S, Hyde TM, Huttner A, Reimers M, Kleinman JE, Sestan N: Spatio-temporal transcriptome of the human brain. Nature 2011; 478:483–489
32.
Colantuoni C, Lipska BK, Ye T, Hyde TM, Tao R, Leek JT, Colantuoni EA, Elkahloun AG, Herman MM, Weinberger DR, Kleinman JE: Temporal dynamics and genetic control of transcription in the human prefrontal cortex. Nature 2011; 478:519–523
33.
Ripke S, Schizophrenia Psychiatric Genome-Wide Association Study (GWAS) Consortium: Genome-wide association study identifies five new schizophrenia loci. Nat Genet 2011; 43:969–976
34.
Ripke S, O'Dushlaine C, Chambert K, Moran JL, Kähler AK, Akterin S, Bergen SE, Collins AL, Crowley JJ, Fromer M, Kim Y, Lee SH, Magnusson PK, Sanchez N, Stahl EA, Williams S, Wray NR, Xia K, Bettella F, Borglum AD, Bulik-Sullivan BK, Cormican P, Craddock N, de Leeuw C, Durmishi N, Gill M, Golimbet V, Hamshere ML, Holmans P, Hougaard DM, Kendler KS, Lin K, Morris DW, Mors O, Mortensen PB, Neale BM, O'Neill FA, Owen MJ, Milovancevic MP, Posthuma D, Powell J, Richards AL, Riley BP, Ruderfer D, Rujescu D, Sigurdsson E, Silagadze T, Smit AB, Stefansson H, Steinberg S, Suvisaari J, Tosato S, Verhage M, Walters JTMulticenter Genetic Studies of Schizophrenia Consortium, Levinson DF, Gejman PV, Kendler KS, Laurent C, Mowry BJ, O'Donovan MC, Owen MJ, Pulver AE, Riley BP, Schwab SG, Wildenauer DB, Dudbridge F, Holmans P, Shi J, Albus M, Alexander M, Campion D, Cohen D, Dikeos D, Duan J, Eichhammer P, Godard S, Hansen M, Lerer FB, Liang KY, Maier W, Mallet J, Nertney DA, Nestadt G, Norton N, O'Neill FA, Papadimitriou GN, Ribble R, Sanders AR, Silverman JM, Walsh D, Williams NM, Wormley BPsychosis Endophenotypes International Consortium, Arranz MJ, Bakker S, Bender S, Bramon E, Collier D, Crespo-Facorro B, Hall J, Iyegbe C, Jablensky A, Kahn RS, Kalaydjieva L, Lawrie S, Lewis CM, Lin K, Linszen DH, Mata I, McIntosh A, Murray RM, Ophoff RA, Powell J, Rujescu D, Van Os J, Walshe M, Weisbrod M, Wiersma DWellcome Trust Case Control Consortium 2, Donnelly P, Barroso I, Blackwell JM, Bramon E, Brown MA, Casas JP, Corvin AP, Deloukas P, Duncanson A, Jankowski J, Markus HS, Mathew CG, Palmer CN, Plomin R, Rautanen A, Sawcer SJ, Trembath RC, Viswanathan AC, Wood NW, Spencer CC, Band G, Bellenguez C, Freeman C, Hellenthal G, Giannoulatou E, Pirinen M, Pearson RD, Strange A, Su Z, Vukcevic D, Donnelly P, Langford C, Hunt SE, Edkins S, Gwilliam R, Blackburn H, Bumpstead SJ, Dronov S, Gillman M, Gray E, Hammond N, Jayakumar A, McCann OT, Liddle J, Potter SC, Ravindrarajah R, Ricketts M, Tashakkori-Ghanbaria A, Waller MJ, Weston P, Widaa S, Whittaker P, Barroso I, Deloukas P, Mathew CG, Blackwell JM, Brown MA, Corvin AP, McCarthy MI, Spencer CC, Bramon E, Corvin AP, O'Donovan MC, Stefansson K, Scolnick E, Purcell S, McCarroll SA, Sklar P, Hultman CM, Sullivan PF: Genome-wide association analysis identifies 13 new risk loci for schizophrenia. Nat Genet 2013; 45:1150–1159
35.
Voineagu I, Wang X, Johnston P, Lowe JK, Tian Y, Horvath S, Mill J, Cantor RM, Blencowe BJ, Geschwind DH: Transcriptomic analysis of autistic brain reveals convergent molecular pathology. Nature 2011; 474:380–384
36.
Allen Institute for Brain Science: BrainSpan: Atlas of the Developing Human Brain, 2011. http://www.brainspan.org/
37.
Kang YH, Galal WC, Farina A, Tappin I, Hurwitz J: Properties of the human Cdc45/Mcm2-7/GINS helicase complex and its action with DNA polymerase epsilon in rolling circle DNA synthesis. Proc Natl Acad Sci USA 2012; 109:6042–6047
38.
Zhang Z, Sun Y, Cho YW, Chow CC, Simons SS: PA1 protein, a new competitive decelerator acting at more than one step to impede glucocorticoid receptor-mediated transactivation. J Biol Chem 2013; 288:42–58

Information & Authors

Information

Published In

Go to American Journal of Psychiatry
Go to American Journal of Psychiatry
American Journal of Psychiatry
Pages: 758 - 767
PubMed: 24874100

History

Received: 3 November 2013
Revision received: 28 January 2014
Accepted: 21 February 2014
Published online: 1 July 2014
Published in print: July 2014

Authors

Details

Rebecca Birnbaum, M.D.
From the Lieber Institute for Brain Development, Baltimore, and the Departments of Psychiatry, Neurology, Neuroscience, and The Institute of Genetic Medicine, Johns Hopkins School of Medicine, Balitmore.
Andrew E. Jaffe, Ph.D.
From the Lieber Institute for Brain Development, Baltimore, and the Departments of Psychiatry, Neurology, Neuroscience, and The Institute of Genetic Medicine, Johns Hopkins School of Medicine, Balitmore.
Thomas M. Hyde, M.D., Ph.D.
From the Lieber Institute for Brain Development, Baltimore, and the Departments of Psychiatry, Neurology, Neuroscience, and The Institute of Genetic Medicine, Johns Hopkins School of Medicine, Balitmore.
Joel E. Kleinman, M.D., Ph.D.
From the Lieber Institute for Brain Development, Baltimore, and the Departments of Psychiatry, Neurology, Neuroscience, and The Institute of Genetic Medicine, Johns Hopkins School of Medicine, Balitmore.
Daniel R. Weinberger, M.D.
From the Lieber Institute for Brain Development, Baltimore, and the Departments of Psychiatry, Neurology, Neuroscience, and The Institute of Genetic Medicine, Johns Hopkins School of Medicine, Balitmore.

Notes

Address correspondence to Dr. Weinberger ([email protected]).

Funding Information

The authors report no financial relationships with commercial interests.
Supplementary Material

Metrics & Citations

Metrics

Citations

Export Citations

If you have the appropriate software installed, you can download article citation data to the citation manager of your choice. Simply select your manager software from the list below and click Download.

For more information or tips please see 'Downloading to a citation manager' in the Help menu.

Format
Citation style
Style
Copy to clipboard

View Options

View options

PDF/EPUB

View PDF/EPUB

Get Access

Login options

Already a subscriber? Access your subscription through your login credentials or your institution for full access to this article.

Personal login Institutional Login Open Athens login
Purchase Options

Purchase this article to access the full text.

PPV Articles - American Journal of Psychiatry

PPV Articles - American Journal of Psychiatry

Not a subscriber?

Subscribe Now / Learn More

PsychiatryOnline subscription options offer access to the DSM-5-TR® library, books, journals, CME, and patient resources. This all-in-one virtual library provides psychiatrists and mental health professionals with key resources for diagnosis, treatment, research, and professional development.

Need more help? PsychiatryOnline Customer Service may be reached by emailing [email protected] or by calling 800-368-5777 (in the U.S.) or 703-907-7322 (outside the U.S.).

Media

Figures

Other

Tables

Share

Share

Share article link

Share