Europe PMC

This website requires cookies, and the limited processing of your personal data in order to function. By using the site you are agreeing to this as outlined in our privacy notice and cookie policy.

Characterisation of genetic variation in ST8SIA2 and its interaction region in NCAM1 in patients with bipolar disorder.

Author information

Affiliations

  • 1. Neuroscience Research Australia, Sydney, New South Wales, Australia; Schizophrenia Research Institute, Sydney, New South Wales, Australia.
      Authors
    • Shaw AD 1
    (1 author)
  • 2. Neuroscience Research Australia, Sydney, New South Wales, Australia; Schizophrenia Research Institute, Sydney, New South Wales, Australia; School of Medical Sciences, Faculty of Medicine, University of New South Wales, Sydney, New South Wales, Australia.
      Authors
    • Tiwari Y 2
    • Schofield PR 2
    • Fullerton JM 2
    (3 authors)
  • 3. Peter Wills Bioinformatic Centre, Garvan Institute, Sydney, New South Wales, Australia.
      Authors
    • Kaplan W 3
    (1 author)
  • 4. Neuroscience Research Australia, Sydney, New South Wales, Australia.
      Authors
    • Heath A 4
    (1 author)
  • 5. School of Psychiatry, University of New South Wales, Sydney, New South Wales, Australia; Black Dog Institute, Sydney, New South Wales, Australia.
      Authors
    • Mitchell PB 5
    (1 author)

ORCIDs linked to this article

Plos one, 20 Mar 2014, 9(3):e92556
https://doi.org/10.1371/journal.pone.0092556 PMID: 24651862 PMCID: PMC3961385

This article is in the Europe PMC Open access subset. Refer to the copyright information in the article for licensing details.
Free full text in Europe PMC

Abstract 

Alpha-2,8-sialyltransferase 2 (ST8SIA2) is an enzyme responsible for the transfer of polysialic acid (PSA) to glycoproteins, principally the neuronal cell adhesion molecule (NCAM1), and is involved in neuronal plasticity. Variants within ST8SIA2 have previously shown association with bipolar disorder, schizophrenia and autism. In addition, altered PSA-NCAM expression in brains of patients with schizophrenia or bipolar disorder indicates a functional dysregulation of glycosylation in mental illness. To explore the role of sequence variation affecting PSA-NCAM formation, we conducted a targeted re-sequencing study of a ∼ 100 kb region--including the entire ST8SIA2 gene and its region of interaction with NCAM1--in 48 Caucasian cases with bipolar disorder using the Roche 454 platform. We identified over 400 DNA variants, including 47 putative novel variants not described in dbSNP. Validation of a subset of variants via Sequenom showed high reliability of Roche 454 genotype calls (97% genotype concordance, with 80% of novel variants independently verified). We did not observe major loss-of-function mutations that would affect PSA-NCAM formation, either by ablating ST8SIA2 function or by affecting the ability of NCAM1 to be glycosylated. However, we identified 13 SNPs in the UTRs of ST8SIA2, a synonymous coding SNP in exon 5 (rs2305561, P207P) and many additional non-coding variants that may influence splicing or regulation of ST8SIA2 expression. We calculated nucleotide diversity within ST8SIA2 on specific haplotypes, finding that the diversity on the specific "risk" and "protective" haplotypes was lower than other non-disease-associated haplotypes, suggesting that putative functional variation may have arisen on a spectrum of haplotypes. We have identified common and novel variants (rs11074064, rs722645, 15:92961050) that exist on a spectrum of haplotypes, yet are plausible candidates for conferring the effect of risk and protective haplotypes via multiple enhancer elements. A Galaxy workflow/pipeline for sequence analysis used herein is available at: https://main.g2.bx.psu.edu/u/a-shaw-neura/p/next-generation-resources.

Free full text 

Logo of plosoneLink to Publisher's site
PLoS One. 2014; 9(3): e92556.
Published online 2014 Mar 20. https://doi.org/10.1371/journal.pone.0092556
PMCID: PMC3961385
PMID: 24651862

Characterisation of Genetic Variation in ST8SIA2 and Its Interaction Region in NCAM1 in Patients with Bipolar Disorder

Alex D. Shaw, 1 , 2 Yash Tiwari, 1 , 2 , 3 Warren Kaplan, 4 Anna Heath, 1 Philip B. Mitchell, 5 , 6 Peter R. Schofield, 1 , 2 , 3 and Janice M. Fullerton 1 , 2 , 3 , *
Jeong-Sun Seo, Editor
Author information Article notes Copyright and License information Disclaimer

Associated Data

Supplementary Materials

Abstract

Alpha-2,8-sialyltransferase 2 (ST8SIA2) is an enzyme responsible for the transfer of polysialic acid (PSA) to glycoproteins, principally the neuronal cell adhesion molecule (NCAM1), and is involved in neuronal plasticity. Variants within ST8SIA2 have previously shown association with bipolar disorder, schizophrenia and autism. In addition, altered PSA-NCAM expression in brains of patients with schizophrenia or bipolar disorder indicates a functional dysregulation of glycosylation in mental illness. To explore the role of sequence variation affecting PSA-NCAM formation, we conducted a targeted re-sequencing study of a ~100 kb region – including the entire ST8SIA2 gene and its region of interaction with NCAM1 – in 48 Caucasian cases with bipolar disorder using the Roche 454 platform. We identified over 400 DNA variants, including 47 putative novel variants not described in dbSNP. Validation of a subset of variants via Sequenom showed high reliability of Roche 454 genotype calls (97% genotype concordance, with 80% of novel variants independently verified). We did not observe major loss-of-function mutations that would affect PSA-NCAM formation, either by ablating ST8SIA2 function or by affecting the ability of NCAM1 to be glycosylated. However, we identified 13 SNPs in the UTRs of ST8SIA2, a synonymous coding SNP in exon 5 (rs2305561, P207P) and many additional non-coding variants that may influence splicing or regulation of ST8SIA2 expression. We calculated nucleotide diversity within ST8SIA2 on specific haplotypes, finding that the diversity on the specific “risk” and “protective” haplotypes was lower than other non-disease-associated haplotypes, suggesting that putative functional variation may have arisen on a spectrum of haplotypes. We have identified common and novel variants (rs11074064, rs722645, 15[ratio]92961050) that exist on a spectrum of haplotypes, yet are plausible candidates for conferring the effect of risk and protective haplotypes via multiple enhancer elements. A Galaxy workflow/pipeline for sequence analysis used herein is available at: https://main.g2.bx.psu.edu/u/a-shaw-neura/p/next-generation-resources.

Introduction

Polysialic acid (PSA) is an important carbohydrate that plays a crucial role in neural development and plasticity [1]–[3] (reviewed in [4], [5]). Homopolymers of α2,8-linked sialic acid are synthesised in the Golgi by sialyltransferase enzymes, of which there are six in the human genome. Two sialyltransferases, encoded by ST8SIA2 on human chromosome 15q26.1 and ST8SIA4 on 5q21.1, catalyse the transfer of long chains (>50 residues) of PSA onto their major substrate, neuronal cell adhesion molecule 1 (NCAM1) [2], [6] (reviewed in [5]), to form PSA-NCAM. Expression of PSA-NCAM is particularly high in early brain development [1]. PSA addition to NCAM1 increases the hydrated volume of cells and limits the adhesive properties of NCAM1, enabling neuronal migration, dendrite formation, axon targeting and synaptic plasticity (reviewed in [4]). Mice expressing only PSA-free NCAM after knockout of the relevant sialyltransferase genes (st8sia2 −/−; st8sia4 −/−), show severe malformation of major brain axon tracts, which are reminiscent of those observed in patients with schizophrenia [7], [8]. Alterations in the expression of PSA-NCAM in the adult post-mortem brain – particularly in the amygdala and hippocampal dentate gyrus & CA1 regions – have been shown in patients with schizophrenia, bipolar disorder, major depression and drug-refractory temporal lobe epilepsy [9]–[12]. Further, expression of PSA-NCAM is modulated in rats by treatment with common antipsychotics [13]–[15].

There is growing genetic evidence of the involvement of ST8SIA2 (also known as SIAT8B or STX) in conferring risk to mental illness. We and others have previously identified the 15q25–26 genomic region as harbouring a mental illness risk gene through linkage in family cohorts with bipolar spectrum disorder [16], bipolar disorder and schizophrenia [17], [18] or psychosis [19]. We subsequently found evidence of disease association with variants in ST8SIA2 in both bipolar disorder and schizophrenia [20]. Other groups have shown association with schizophrenia in Japanese and Chinese cohorts in the promoter region of ST8SIA2 [21], [22], which was selected as a candidate gene after NCAM1 was implicated in disease risk [23]. More recently, genome-wide association studies have implicated SNPs downstream of ST8SIA2 in risk of bipolar disorder [24], and variants in intron 2 have also been associated with autism spectrum disorder (verbal subtype) [25]. Furthermore, a 520 kb hemizygous deletion of three genes including ST8SIA2 has been identified in a patient with autism spectrum disorder [26]. Hence, characterisation of the nucleotide variation in ST8SIA2 in patients with a variety of mental illnesses is an important step in understanding the molecular mechanisms through which this gene increases disease risk.

Recent advances in high-throughput sequencing have allowed the rapid sequencing of whole gene loci in multiple individuals [27]. We have exploited this technique to sequence the coding, intronic and flanking sequence of ST8SIA2 and its region of interaction in NCAM1, in lymphocyte DNA from individuals with bipolar disorder. We incorporated a number of methods into our analysis pipeline for quality control in variant detection and base-calling. Altogether, we have identified 396 single nucleotide variants including 42 novel single nucleotide variants across 95 kb of genomic DNA in ST8SIA2, and 21 variants, 5 of which are novel, in a 6 kb region of NCAM1. We utilised data from the 1000 Genomes Project [28] to identify variants unique to bipolar individuals, and those which are at altered frequency in bipolar cases compared to the 1000 Genomes population panels [29] with the same racial background. We examined the haplotypic background of each variant in ST8SIA2 with respect to a previously identified risk haplotype [20] and explored evidence of functional impact by comparisons with ENCODE data [30], [31] and co-localisation with ST8SIA2-specific enhancer elements that contribute to altered target gene expression [32]. Through an integrative analysis of the above methods, we have identified and verified the most compelling candidate variants for increasing susceptibility to mental illness by modulating ST8SIA2 function.

Results

454 Sequencing and Mapping

To characterise genetic variation that may quantitatively impact formation of PSA-NCAM, we performed massively parallel sequencing on the 454 GS-FLX platform of two regions that we hypothesise are associated with ST8SIA2-related developmental pathology in mental illness: 1) a ~95 kb region including the entire coding and intronic regions and ~20 kb of flanking sequence of ST8SIA2 (hg19/GRCh37: chr15[ratio]92,919,255–93,013,920 bp; Figure 1); and 2) a ~6 kb region within NCAM1 (hg19/GRCh37: chr11[ratio]113,100,972–113,107,571 bp; Figure 2). This latter region contains the 5th immunoglobulin-like domain and fibronectin-1 domain, and specifically the 5th and 6th glycosylation sites and acidic patch (that are critical for enzyme recognition, docking and polysialylation by ST8SIA2 to form PSA-NCAM [33]–[35].

An external file that holds a picture, illustration, etc. Object name is pone.0092556.g001.jpg
Coverage of ST8SIA2 gene locus by 454 and Sanger sequencing.

The 96~20 kb flanking sequence (chr15[ratio]92,919,255–93,013,920 bp; UCSC Genome Browser hg19 http://genome.ucsc.edu) was divided into eight long-range PCR amplicons (454 Amplicons 1–8, in pink) for library preparation prior to 454 sequencing. The mean depth of coverage across all 48 samples for each long-range amplicon is shown in burgundy. A GC-rich region surrounding the promoter and exon 1 was sequenced via Sanger sequencing using 8 short overlapping amplicons (Sanger amplicons 1–8; inset). Pink dots on the GC percent or mean depth per sample tracks indicate loci where the level of GC richness, or read depth is above the range indicated. Regions which failed to amplify, or which contained simple repeats, long homopolymers (length > = 10), or short interspersed elements (SINEs) tended to have lower depth of coverage, and were inaccessible to next-generation sequencing in phase I of the 1000 genomes project (1000 G Accs strict).

An external file that holds a picture, illustration, etc. Object name is pone.0092556.g002.jpg
Region of NCAM1 gene covered via Roche 454 sequencing.

The sequenced region spanned 6.6(chr11[ratio]113,100,972–113,107,571 bp; hg19) and included exons 8–12 (red boxes) and less frequently used exons 8a and mini-a (grey boxes) [52]. The locations of the 5th and 6th glycosylation sites and acidic patch required for glycosylation are in exons 9, 10 and 11 and are indicated by green, blue and white lines, respectively. The locations of SNPs identified in this study are indicated with arrows. Names of known SNPs are in black text, and novel SNPs are listed as NV (novel base substitutions) or InDel (novel insertion-deletion change) followed by their base position (hg19). Two coding SNPs located in exons 8a and 11 are marked with closed boxes, with their effect on amino acid sequence shown.

We applied BWA mapping [36] to the entire genome to ensure that mapped reads were specific to regions amplified by long-range PCR, and to determine rates of target enrichment and non-specific amplification. We found that 96% of reads across all individuals mapped to a unique site within the human genome, and that 93% of reads that mapped to the human genome mapped within the target regions. The remaining mapping is likely due to some amplification of non-targeted regions by the long-range PCR primers.

Coverage

Depth and uniformity of coverage are critical for sensitive SNP detection and accurate genotype calls, and are important considerations in assessing of the completeness of variant detection in targeted resequencing. Hence, we examined the mean depth of coverage across each base-pair of the ST8SIA2 region (Figure 1). We found that depth of coverage dropped substantially and locally in regions of low complexity, which included simple repeat sequences and very long homopolymers (length ≥10). These regions are known to be difficult to base call accurately due to the deficiencies in sequencing chemistry of 454 [37]. Indeed, many of these regions were found to be less accessible to sequencing according to the “strict” stringency assessment of the human genome conducted in phase I of the 1000 Genomes Project [28]. Many of these regions also coincided with the ends of primate-specific (recently evolved) Alu repeat elements, which are rich in A and T homopolymers [38].

Aside from local sequence context, a major determinant of mean depth of coverage was inclusion within a specific long-range PCR amplicon, with the NCAM amplicon and ST8SIA2 amplicons 2, 5 and 6 having the highest mean coverage (Figure 3). We assessed the depth of coverage in each sample, and also found variation between samples (Figure 3). These differences demonstrate that amplicon pooling and sample pooling steps each contributed to variation in the depth of sequence coverage within genomic regions. Overall, although some samples contained amplicons with poor coverage (Figure 3), the mean coverage for each amplicon across all samples was >10, a level sufficient for ascertainment of diploid genotypes.

An external file that holds a picture, illustration, etc. Object name is pone.0092556.g003.jpg
Heatmap of mean depth of coverage per amplicon and per sample.

Amplicons are listed along the x-axis, and samples (x1-48) are listed along the y-axis. Depth of coverage is indicated according to the colour key, with white indicating <2 reads coverage and red indicating >35 reads coverage. The plot was created using R [63].

SNP calling

We used two different algorithms to call SNPs from 454-generated data: i) GS Reference Mapper (Refmapper), which is a vendor-supplied algorithm specific to the 454 platform, and calls SNPs based on alignment data from each individual independently; and ii) GATK [39], [40], which is not platform specific, and calls SNPs based on alignment data from all individuals combined. As raw GATK SNP calls require some level of post-hoc filtering, we initially used the GATK-recommended filtering criteria for whole genome SNP detection, and then increased the stringency to remove a subset of likely spurious variant calls (further details are supplied in Note S1; Figure S1). In order to maximise our ability to detect true SNPs for functional assessment, the union of both Refmapper and GATK algorithms (i.e., those SNPs that were called by at least one algorithm) formed our list of putative SNPs.

We compared the locus of SNPs called using each algorithm and found that 89% of called SNPs were common to both algorithms (Table 1). Taking only the subset of putative novel SNPs in each set, we found that each algorithm alone and the intersection between the algorithms generated approximately the same proportion of putative novel variants (8–9%), and ~60% of the putative novel variants were called in both algorithms (Table 1). Twenty one SNPs filtered out in GATK were called as genuine variants by Refmapper. A further 23 SNPs were identified by only one algorithm.

Table 1

Comparison of SNP detection methods from 454 data, and identification of bona fide SNPs for functional analysis.
Setn
RefmapperAll383
Novel34 (9%)
GATK (passed filter)All349
Novel28 (8%)
GATK (failed filter)All232
IntersectionAll344
Novel24 (7%)
UnionAll388
Novel38 (10%)
SangerAll24
Novel5
The number of SNPs (n) which were detected individually by Refmapper and GATK are shown. The SNPs that were called by both algorithms (and passed filtering) are given in the intersection, and the union represents SNPs that were called by at least one package. SNPs detected in the 5.3 kb region sequenced via Sanger sequencing are also listed.

We sequenced a GC-rich region that was not amplified by long-range PCR using bi-directional Sanger sequencing (Figure 1). SNPs were called using Consed and the Phred/Phrap/Polyphred software suite [41], [42], with manual inspection of chromatograms for SNP identification and genotyping. Non-reference SNP and genotype calls were only made if the non-reference base was called in both the forward and reverse reads. Using this methodology we identified 24 additional SNPs in the ST8SIA2 locus (19 known, 5 novel).

Genotype calling and accuracy

Diploid genotype calling from high throughput sequencing relies on a minimal level of read coverage in each individual at each base position to accurately distinguish true heterozgous sites from sites with random base-calling errors. Due to variation in coverage across regions and across amplicons (Figure 3), we reasoned that a proportion of the targeted region in each individual was marginal, i.e. was covered to a depth below that required for statistically confident heterozygote calling (coverage <10). In order to use data from these regions, and other regions with low local coverage, we used the probabilistic genotype calling module within GATK. This module uses read data at all sites, even those with low coverage, in the simultaneous estimation of cohort allele frequency and calculation of genotype likelihoods. The diploid genotype with the highest likelihood is assigned as the genotype call, along with a likelihood ratio calculation of the confidence of this call [denoted Genotype Quality (GQ)].

In order to assess the overall accuracy of applying this mode of genotype calling to our 454 sequence data, we compared the GATK 454-derived genotypes with genotypes determined directly via Illumina genotyping assay for the same individuals from an earlier study [20], [43], which were assumed to be 100% accurate. The available Illumina data contained 2133 genotype calls from 50 SNPs (minor allele frequency [MAF] average 0.29; range 0.01–0.49), including 1081 (50.7%) non-homozygous-reference calls. One sample (# 10) was excluded from further analysis based on high genotype discordance between platforms (~50%). The overall genotype discordance across the remaining individuals was low (~2%; Table 2). We reasoned that some genotype calls may have been made at low confidence in regions with low coverage, and that these could be removed from all genotype calls to reduce discordance. Indeed we found that eliminating low confidence genotype calls (GQ <10), which are most subject to random error, reduced the discordance 2-fold at the expense of a decrease in the number of genotypes called (Table 2).

Table 2

Genotype concordance between 454-generated genotype calls and those directly genotyped via Illumina GoldenGate or Illumina 660 W.
SetCall accuracyCall rateHet call accuracyHet call rate
All0.980.950.950.93
Low confidence removed0.990.890.980.89
Genotype calls are divided into two categories: all calls; or confident genotype calls based on Genotype Quality (GQ) ≥ 10. Call accuracy gives the concordance between the SNP genotype derived from 454 data and Illumina-generated data across all samples excluding missing data, and the call rate incorporates missing data on an individual basis. These are calculated for all genotypes, and heterozygous (het) calls only.

Genotype calling for Sanger data was performed using Polyphred. Only two SNPs were available with known genotypes from the Illumina genotyping data for assessing the accuracy of Sanger genotype calling. We found that all known genotypes were accurately called in 44 individuals with available known genotypes (88 known genotypes, 14 non-reference).

Comparison to 1000 Genomes Project data

Having ascertained a list of putative SNPs that were observed in Caucasian individuals with bipolar disorder, we were able to assess our ability to detect SNPs present in the general population by comparing our detected variants to those of Caucasian Europeans represented in phase I of the 1000 Genomes (1 kG) project (CEU and GBR populations; n = 174).

We found that many rare SNPs (MAF <0.05) were found to be unique to either the bipolar cohort or 1 kG individuals (Figure 4), which is expected due to the inevitable sampling bias which occurs when sampling small numbers of subjects from a population. In contrast, 100% of common variants (MAF ≥0.1) found in the 1 kG cohort were also found in our bipolar cohort, confirming our ability to detect common variants in the population. A small number of common variants found in our cohort (n = 3) were apparently not detected in 1 kG: however, manual inspection showed that these SNPs were adjacent to indels, or low-complexity sequence, and therefore were probably filtered from the list of integrated variants in the 1 kG cohort. We excluded these putative SNPs from further analysis.

An external file that holds a picture, illustration, etc. Object name is pone.0092556.g004.jpg
Comparison of detected SNP variants with all SNPs detected in 1 kG data.

The minor allele frequency (MAF) of all SNPs identified in the sequenced region from 174 Caucasian Europeans from phase I of the 1 kG project (n = 519) and 48 bipolar disorder cases (n = 412) were compared. SNPs found in bipolar individuals only (yellow) and 1 kG individuals only (light blue) were predominantly rare alleles whereas those found in both sets (green) were predominantly common.

To determine which variants might be associated with illness, we next compared the allele frequencies of all SNPs detected in ST8SIA2 and NCAM1 in the 48 bipolar cases with those reported in the 174 1 kG individuals. While the 1 kG population may contain individuals with bipolar disorder, this would be expected to be at the population rate of around 1.5%, so it can be treated as an unscreened control population. We found that two SNPs in NCAM1 (rs584427 and rs646558, representing ~10% of all observed) and 53 SNPs in ST8SIA2 (~14% of all observed) showed a frequency difference of 5% or more between the two populations (Table 3, Table S1).

Table 3

Summary of genetic variation identified in NCAM1.
SNPChr11 Positionpolya Freq (BP)Freq (1 kG)FreqDIFF
rs17115204113101143T/C0.0110.0060.005
11[ratio]113101620113101620-/T0.0300.0000.030#
Exon 8a
rs117219783*113101958A/G0.0210.0120.010
11[ratio]113102363 v113102363-/C0.0100.0000.010
Exon 8
11[ratio]113102679113102679G/C0.0100.0000.010
Exon 9
rs2241450113103259T/C0.0110.0060.005
Exon 10
Exon 11
rs584427* v 113103996 T/G 0.340 0.450 −0.110#
rs148753484113104198-/CACCCCAA0.1100.0880.022
rs74800050113104456T/C0.0320.0180.014
11[ratio]113104997113104997-/T0.0200.0000.020#
11[ratio]113105039113105039-/T0.0100.0000.010
rs117660538113105158G/C0.0110.012−0.001
rs4646982 v113105405G/T0.3190.3000.019
rs113886534113105410G/A0.0210.0120.010
Exon 12
rs646558 v 113105907 A/C 0.245 0.194 0.051
rs12279261113106455G/A0.1170.159−0.042
rs17597344113106464A/G0.0430.065−0.022
rs189437769113106527T/C0.0110.012−0.001
rs11214546113106528T/G0.0430.053−0.010
Mini-exon a
rs111567425113107077A/G0.0110.035−0.025
rs652285113107204A/G0.0960.0530.043
The location of each variant (NCBI 37/hg19 build) relative to known exons is shown, with exon locations indicated, and coding SNPs indicated with asterisks. The minor allele frequency (MAF) of each variant in all samples (BP+1 kG) and the BP and 1 kG cohorts separately is given. The frequency difference (freqDIFF) was calculated relative to 1 kG allele frequency, with absolute frequency differences greater than 5% indicated in bold text, and those with p values <0.1 indicated with an asterisk.
aFor each polymorphism, the minor allele is listed first.
vSNPs verified by direct genotyping are indicated.

Nucleotide diversity on specific ST8SIA2 haplotype backgrounds

We previously identified specific haplotypes which were associated with increased or decreased disease risk, and were defined by six SNPs within the ST8SIA2 locus [20]. In order to identify potentially functional variants which lie specifically on these haplotypes, we performed haplotype phasing and imputation using the BEAGLE algorithm. The data used as input was combined from Illumina genotypes, Sanger genotypes and 454 genotype likelihoods from GATK (excluding data in regions with potential allelic bias; see Note S2). Genotypes from 1 kG phase I Caucasian European individuals (n = 174) were used as a reference panel.

After phasing and imputation, 40 bipolar chromosomes were predicted to have the risk haplotype, whereas 13 bipolar chromosomes were predicted to have the protective haplotype (Figure 5). The remaining 41 bipolar chromosomes harboured neither risk nor protective haplotypes, and were designated as “other”. A large number of common variants in the bipolar individuals appear to have arisen through co-segregation of alleles in a large haplotype block between the 1st and 5th exons of ST8SIA2 (Figure 5).

An external file that holds a picture, illustration, etc. Object name is pone.0092556.g005.jpg
Haplotype phasing of all detected nucleotide variants in ST8SIA2 gene region.

The gene structure is shown above, with the minor allele frequency of each SNP detected shown by the histogram immediately below the gene (red represents frequency of alternate allele). The red, green and yellow bars on the left indicate the haplotype on which the detected variation resides (risk = red; protective = green; other haplotypes = yellow). Detected variation with respect to the specific haplotype on which it resides is shown to the right of haplotype bars, where non-reference alleles are coloured in light blue (homozygous) and dark blue (heterozygous). The locations of the six SNP markers that were used to define the haplotypes are shown below the variation. The location of the PreSTIGE enhancer elements within the 454 sequenced region from human skeletal muscle myoblasts (HSMM) are indicated in orange, and from neuronal precursor cells (NPC) are indicated in purple. PreSTIGE enhancer locations were obtained from http://genetics.case.edu/prestige/).

The nucleotide diversity across the sequenced region, as measured under the assumptions of an infinite site neutral allele model, was similar across ST8SIA2 (θ = 5.42×10−4) to rates previously reported for other genes [44]. Surprisingly, the nucleotide diversity was lower on the designated risk and protective haplotypes (θ = 3.45×10−4 and 3.74×10−4) compared to other haplotypes (θ = 8.35×10−4). This may be indicative of allelic heterogeneity, i.e., that putative functional variants may have arisen on a spectrum of haplotypes, rather than restricted specifically to the risk haplotype previously associated with illness.

Analysis of putative function

Having ascertained the location, diploid genotype and phase of identified SNPs, we turned our attention to the predicted functional impact of SNPs in the ST8SIA2 and NCAM1 sequenced regions, looking at the coding, transcribed and non-transcribed regions.

We found no non-synonymous SNPs in coding regions of ST8SIA2, although we identified a previously known synonymous SNP in exon 5 (rs2305561, P207P). We also identified a known SNP 5 bp from the 3′ end of exon 5 (rs11637874), for which the minor T allele was predicted to improve splice donor site efficiency (BDGP score: T allele = 0.95; C allele = 0.81). Allele frequencies for these two SNPs were similar to those in the 1 kG cohort (Table S1).

In NCAM1, we found a single non-synonymous SNP in exon 8a (rs117219783, G/A, G>E; Figure 2), which resulted in a non-conservative amino acid change (non-polar glycine to polar glutamic acid), although the MAF was not different to that observed in 1 kG data (MAF(A) = 0.02 vs 0.01; χ2 = 0.37, p = 0.54; Table 3). This exon is rarely used, and was reported to exist in mRNA from testis (AK314589). We also observed a synonymous variant in exon 11 of NCAM1 (rs584427, G/T, V540V), which showed some frequency distortion in the bipolar cases compared to 1 kG individuals (MAF(T) = 0.34 vs 0.45 respectively; χ2 = 3.16, p = 0.08; Table 3). Two novel insertion/deletion SNPs (at 113,101,620 and 113,104,997 bp) showed nominal allelic association (χ2 = 5.37; p = 0.02 and χ2 = 3.57; p = 0.05, respectively). No variants were observed in the glycosylation regions or the acidic patch that would affect the ability of NCAM1 to be glycosylated.

The recent release of comprehensive functional genomics data [31] allowed us to assess the potential gene-regulatory impact of the non-coding variants we identified in ST8SIA2. To identify SNPs that may affect transcriptional regulation of ST8SIA2 we looked at the localisation of SNPs within DNAse hypersensitivity sites (DHSs) in human Embryonic Stem Cells (hESC) or neuronal cell lines, surveyed in the ENCODE project [31], [45]. DHSs indicate a region of open chromatin, and DHS peaks (DHSPs) within the DNAse-Seq signal indicate chromatin that is potentially accessible for transcription factor binding, which could be modulated by sequence variants.

We found that 20 SNPs in total (five novel) were within DHSPs, five SNPs were within DHSPs in both hESC and neuronal cell lines (one novel), six within neuronal cell lines only (two novel) and nine within hESC cell lines only (two novel; Table 4). Consistent with enrichment for functional variants, assessment of Genomic Evolutionary Rate Profiling (GERP) scores revealed that 20% of DHSP-overlapping sites were conserved (GERP >2) compared to only 3% of non-DHS-overlapping sites. The five novel SNPs which were coincident with DHSPs (at 92937845, 92938435, 92962229, 92973414, 92984673) were present in only a single sequenced individual each. Two of these SNPs were conserved with a high GERP score (4.12 at 92937845 and 5.55 at 92973414) and one SNP was predicted to introduce a splice acceptor site 69 bp downstream of exon 2 (92973414, BDGP score = 0.98 with A allele). Despite the conservation of a high proportion of DHSP-overlapping SNPs, the majority (n = 9) of conserved SNPs were neither within DHS nor transcribed. This indicates that either their transcriptional regulatory function is not captured in the cell lines surveyed, or that they impart a conserved function through another mode of regulation. A further survey of DHS clusters across 125 additional cell lines included in the ENCODE project revealed that the majority (n = 6) were not present in DHSs of any other cell lines (data not shown), suggesting that for these SNPs an alternate mode of regulation is responsible rather than a conserved transcriptional regulatory function in another cell type.

Table 4

Summary of transcribed and DHSP overlapping genetic variation identified in ST8SIA2.
SNPChr 15 position novelmA Excl. to RISKmA Excl. to PROTSetDHSPGERPpolya Freq (BP)Freq (1 kG)freqDIFF
rs117156502929269662211,2T/C0.0220.0090.014
rs28496195929269681011,2A/G0.0110.035−0.023
rs16946843 92928587 ND ND 3a 2 T/C ND 0.101 ND
rs18669914992935000NDND42,30.0100.01
rs79668601929365840043A/G0.0560.063−0.008
rs13839384692936774NDND42,34.280.0100.01
rs142883061 92937137 0 1 4 1 −2.84 C/T 0.053 0.106 −0.053
5′UTR&exon1
rs3743365*929372680041.08G/T0.10640.1523−0.0459
rs3743364*929372760041.84C/T0.10640.1523−0.0459
rs19222439592937836NDND41,32.820.0100.01
15[ratio]9293784592937845Y2241,34.12A/G0.01100.011
rs3784745 92937858 0 0 4 1,2,3 −0.45 G/C 0.0851 0.1437 −0.0586
rs3784744 v 92937884 0 1 4 1,2,3 1.02 C/G 0.053 0.106 −0.053
15[ratio]9293843592938435YNDND42A/G0.01100.011
rs73543829929388260013C/T0.01100.011
rs11637377929399980113T/C0.0220.055−0.032
rs11857396929401220113A/G0.3890.3790.01
rs722645 v 92943323 0 1 1 3 G/A 0.489 0.563 −0.074
rs28680784929433490113T/C0.4240.3910.033
rs921846 v 92943966 1 0 1 2 A/G 0.489 0.555 −0.066
rs11074064v 9294466201133.91A/G0.4570.471−0.015
rs11074065 92944835 0 1 1 3 A/G 0.371 0.48 −0.109#
rs118530839295114201130.15G/A0.4360.4110.025
rs17599821929513580112,3T/C0.410.4080.002
rs118032995929528340013T/A0.0130.0060.007
rs11637898929528500113A/G0.4380.4340.004
15[ratio]9296222992962229Y003a1,2C/TND0ND
15[ratio]9296939292969392YNDND3b3ND0ND
Exon2
15[ratio]92973414v 92973414Y22115.55A/G0.01100.011
Exon3&4
rs112804054929819032211C/A0.0110.0090.002
15[ratio]9298467392984673Y0012A/G0.01100.011
Exon5
rs2305561*92987938001−0.19G/C0.10.115−0.015
rs5958585992998684NDND3b3ND0.135ND
rs76405928 93000441 0 0 1 3 T/A 0.457 0.379 0.077#
rs79846035930009951013−2.22A/G0.0540.0490.005
rs74895969930010350113A/G0.0430.055−0.011
rs4777988930024440012G/A0.4240.431−0.007
Exon6&3′UTR
rs139149207*v 930077850015.45A/C0.0110.0030.008
rs2290492* v 93007974 0 0 1 4.64 A/G 0.278 0.221 0.057
rs12904773*93008168011G/C0.0650.0370.028
rs8035760* v 93008472 0 0 1 A/T 0.359 0.19 0.169#
rs1869774* v 93010177 0 0 1 2 T/C 0.315 0.233 0.082#
rs115781738*930101922212C/G0.01100.011
rs116928729*93010321001A/G0.0220.0170.005
rs17600420*93011135001G/A0.3040.345−0.041
rs117763930*93011237001A/G0.0220.0170.005
rs145948851*93011468221C/T0.01100.011
rs11853992 v 93012351 0 0 1 3 2.51 G/A 0.33 0.2787 0.0513
rs147920939930124260013−1.12A/G0.0110.00860.0024
The location of transcribed and DHSP overlapping variants on chromosome 15 (hg19 build), with transcribed SNPs indicated with asterisks. SNP alleles (mA) identified exclusively on the risk haplotype are shown (0 = not exclusive, 1 = present on risk and other haplotypes, 2 = present on risk haplotype only, ND = not determined). SNPs identified exclusively on the protective haplotype are shown (0 = not exclusive, 1 = present on protective and other haplotypes, 2 = present on protective haplotype only, ND = not determined). The set from which the SNP was observed is given (1 = GATK & Refmapper; 2 = GATK only; 3a = Refmapper only; 3b = Refmapper only & filtered in GATK; 4 = Sanger). Co-localisation with DNase I hypersensitivity site (DHS) peaks (neuronal = 1; hESC = 2; foetal brain = 3) are given. Genomic Evolutionary Rate Profiling (GERP) scores are provided for each variant that is within a GERP-conserved element. The nature of the polymorphism in each cohort is given (with minor allele listed first). The minor allele frequency (MAF) of each variant in the 47 bipolar cases and 174 Caucasian individuals (CEU and GBR) from the 1000 Genomes Project (1 kG) are shown separately. The frequency difference (freqDIFF) was calculated relative to 1 kG allele frequency, with absolute frequency differences greater than 5% indicated in bold text, and those with p values <0.1 indicated with an asterisk.
aFor each polymorphism, the minor allele is listed first.
vSNPs verified by direct genotyping are indicated. The annotated list of all ST8SIA2 variation identified is given in Table S1.

To examine SNPs that may affect post-transcriptional regulation, we examined SNPs that were within ST8SIA2 untranslated regions, according to Refseq annotation. These SNPs may impact RNA-binding proteins, which are known to play a crucial role in RNA localisation, stability and splicing of neuronally expressed transcripts, which is in turn critical during brain development [46]. We identified two SNPs in the 5′UTR (rs3743365, rs3743364), and 10 SNPs within the 3′UTR, including one variant rs2290492 that lies ~350 bp past the stop codon, near predicted miRNA binding sites [47]. A large proportion of the ST8SIA2 3′UTR proximal to the stop codon is conserved and we found two SNPs in the 3′UTR with GERP >2 in conserved elements (rs139149207 and rs2290492).

Haplotype-specific functional variants

Finally, we examined how SNPs segregated with haplotypes in ST8SIA2 that were associated with mental illness. We considered the 376 SNP loci genotyped by GATK or Sanger. Of these, 347 were successfully phased and thus could be assessed for whether variants segregated with the risk, the protective or other haplotypes.

To look for variants that could confer the risk or protective effect of our previously identified haplotypes in isolation, we looked for variants that resided exclusively on the risk or protective haplotypes and not other haplotypes. We found that 19 variants (seven novel) lay exclusively on the risk haplotype (MAF = 0.011–0.024), and five variants (one novel) lay exclusively on the protective haplotype (MAF = 0.011–0.025). Of those on the risk haplotype, two are transcribed (rs115781738 and rs145948851), and five were within DHSPs (rs117156502, rs112804054, rs115781738, 92937845 and 92973414). The latter two of these sites were conserved with GERP >2. The minor alleles of five SNPs were found exclusively on the protective haplotype: one of these was novel (at base position 92956314) and one localised within a DHSP (rs141457407). All of the SNPs that reside exclusively on the risk or protective haplotype are rare, and therefore cannot individually account for the increased risk or protective effect.

A recent study demonstrated that clusters of common disease-associated variants, located within multiple enhancer elements in cis, can act in combination to regulate gene expression and increase disease risk [32]. Alleles of such SNPs would contribute to a risk or protective effect when inherited together with other regulatory alleles in a specific haplotype, but would not necessarily have a discernible regulatory effect when inherited as part of other haplotypes. To attempt to capture this class of variant, we looked for alleles that segregated with either the risk or the protective haplotype in at least one chromosome, regardless of segregation with other haplotypes in additional chromosomes, and intersected these variants with ST8SIA2-specific PreSTIGE enhancer predictions [32], which existed for two cell lines: neuronal precursor cells (NPC) and human skeletal muscle myoblasts (HSMM), the locations of which are indicated in Figure 5. We expected that alleles having a discernible effect via this mechanism would be relatively common and present in the general population, so we used phased genotypes from 1 kG CEU and GBR population panels [28], which provide a larger population sample than our bipolar cohort.

Of 348 chromosomes in the 1 kG individuals, 127 carried the risk haplotype and 42 carried the protective haplotype. Seventy-eight SNPs segregated with the ST8SIA2 6-SNP risk haplotype in ≥ 1 individual and the protective haplotype in 0 individuals, including two SNPs which were used to tag the originally identified ST8SIA2 risk haplotype. Eight SNPs were tagged by two of six SNPs constituting the risk haplotype (rs2035645 and rs4777974). Of these, three are in PreSTIGE enhancer elements (rs11074066, rs921846 and rs722645), one is also in a DHSP in hESC cells (rs921846) and one is in a foetal brain tissue DHSP cluster (rs722645). Conversely, 40 SNPs segregated with the protective haplotype in ≥ 1 individual and the risk haplotype in 0 individuals. Twenty-five SNPs were tagged by four of six SNPs constituting the protective haplotype (rs4586379, rs11637898, rs11074070 and rs3784735). Of the tagged SNPs, 13 were in PreSTIGE enhancers. Although none of these 13 were in DHSPs from hESC or neuronal cell lines, two (rs11074064 and rs11074065) were in a foetal brain tissue (DHSP) cluster.

Validation and association of putative functional variants

We selected a subset of putative functional variants for validation and association analysis via direct genotyping, in a larger cohort of bipolar cases and control individuals (n = 215 cases, including the 47 sequenced individuals; and n = 165 controls). Variants were selected for Sequenom assay design from those that: 1) lie within PreSTIGE enhancer elements or DHSPs; 2) have a GERP score >2; 3) have an allele frequency difference between bipolar and 1 kG populations of >0.05. Of the 57 SNPs selected, 12 failed assay design or quality control stages, leaving 45 SNPs successfully genotyped via Sequenom (14 novel, 31 known). We also performed manual genotyping of three 3′UTR SNPs via restriction assay (rs2290492 and rs1869774) or direct sequencing (rs8035760).

The concordance rate for Sequenom-derived versus 454-derived genotypes was 97%, just below the concordance rate from Illumina genotyping assay validation. For manually genotyped variants, we found that rs8035760 had a low rate of concordant genotype calls (56%) compared to the other 3′UTR variants (97%). This was likely due to erroneous rs8035760 genotype calls from 454 sequencing, consistent with a low SNP call quality (QD = 6) and proximity of the site to a short homopolymer in the human reference sequence (data not shown).

Eleven of the 14 novel variants (79%) were validated in one or more samples (Table 5). Of the three SNPs that were not validated, two were called by Refmapper only whereas the third was called by both GATK and Refmapper. We also found additional cases from the larger cohort who were heterozygous for one of four of the genotyped novel SNPs, including SNP 15[ratio]92961050, which was found in an additional three cases (two bipolar type I and one bipolar type II), and no control individuals (uncorrected p = 0.079). This SNP is within a PreSTIGE putative enhancer region in intron 1. All remaining SNPs that were found in additional cases were also found in at least one unaffected control individual.

Table 5

Fourteen novel SNPs assayed by Sequenom for verification.
Base Position (hg19)Reference AlleleAlternative AlleleAC in Affected (n = 215)AC in Unaffected (n = 165)Set
11: 113102360TNP003a
15: 92926056GA201
15: 92929229GC101
15: 92956314AG101
15: 92958227CG101
15: 92958793CA101
15: 92959424ANP001
15: 92961050AG403a
15: 92961385AG111
15: 92973414GA101
15: 92992055CT201
15: 92993375CT211
15: 92994038TNP003a
15: 92994449GA111
Three SNPs that were detected in 454-generated sequence were non-polymorphic (NP) after direct genotyping. The allele count (AC) is given from a total of 380 subjects (215 cases including the 48 sequenced individuals, 165 controls). The algorithm set from which the SNP was observed is provided (1 = GATK & Refmapper; 2 = GATK only; 3a = Refmapper only; 3b = Refmapper only & filtered in GATK; 4 = Sanger).

Next, we performed a test of association with bipolar disorder (n = 215 cases, 165 controls). Thirteen subjects were excluded from association analysis due to high rates of genotype missingness (>10%; five cases and five controls). Two SNPs were nominally associated (rs11074064 & rs722645, p<0.05; Table 6) and lie in an LD block within intron 1 (chr15[ratio]92,940,731–92,944,864; Figure S2), although these SNPs did not remain significant after multiple testing correction for 11 independent tests. Interestingly, both rs11074064 and rs722645 were among those found to segregate only with the protective or risk haplotype, with the same direction of effect as the protective and risk alleles respectively. Both rs11074064 and rs722645 also lie within the same PreSTIGE enhancer element spanning the TSS and part of intron 1 and are present within each of two neighbouring DHSPs in foetal brain (data not shown). Thus, these variants are plausible candidates for functional SNPs that contribute to the effect of the protective or risk haplotypes respectively.

Table 6

SNPs nominally associated with bipolar disorder.
SNPbase Position (hg19)Alt alleleFreq (Affected)Freq (Unaffected)Ref alleleP-valueOdds ratio (95%CI)
rs1107406492944662A0.4230.509G0.0190.71 (0.53–0.95)
rs72264592943323G0.4590.381A0.0351.38 (1.02–1.85)
The single nucleotide polymorphisms (SNP) nominally associated with bipolar disorder (p<0.05) are given. The SNP position on chromosome 15 is given relative to hg19 build. The reference (Ref) and alternative (Alt) alleles are shown, with the frequency (Freq) of the alternative allele in affected individuals (n = 210 cases) and unaffected controls (n = 160) shown. Odds ratio of effect and the 95% confidence interval (CI) are given.

Discussion

There is growing evidence of the involvement of ST8SIA2 in risk of mental illness, with the association of variants within the gene with bipolar disorder [20], schizophrenia [21], [22] and autism [25] (reviewed in [48]). Post-mortem studies also implicate altered levels of PSA-NCAM in brains of people with bipolar disorder, schizophrenia and major depression [9]–[11]. Formation of PSA-NCAM is critical for fibre tract formation, neurogenesis and plasticity (reviewed in [4]), and is dependent on availability and functionality of both enzyme and substrate. Animal studies that have begun to assess combinations of mutant alleles, have shown the amount of NCAM devoid of PSA during brain development is critical, and highlights the important balance between these molecules for brain development and function [7], [49] (reviewed in [50]). Further, loss-of-function mutations have previously been reported in ST8SIA2 in schizophrenia [21], [51]. Hence, we conducted a targeted resequencing study to identify genetic variants which may impact formation of PSA-NCAM in bipolar disorder. We targeted both ST8SIA2 and the region of interaction in its principal substrate NCAM1 [33]–[35]. NCAM1 has also shown genetic association with bipolar disorder [23], schizophrenia [52] and neurocognition [53].

In any sequencing study, it is important that the coverage is adequate to be confident that all detectable nucleotide variation is identified. We performed a number of tests for quality control in our data to ensure that the nucleotide diversity we report herein is accurate and complete. We used GATK, which has so far been predominantly used for genome or exome-wide studies, as our principal method for diploid genotype calling. Using the probabilistic genotyping module within GATK resulted in high genotype accuracy, even when including genotypes called on sites with low coverage. Diploid genotyping is critically important for detailed genetic analysis, including haplotype reconstruction and accurate assessment of allele frequency within disease or general populations, and relies on genotypes that are as complete and accurate as possible given the current limitations of sequencing data. In comparison, the use of a simple threshold approach to distinguish between homozygote and heterozygote variants requires a high level of coverage for accurate heterozygote calling. With this approach many regions that have systematically low coverage with current sequencing technology would be missed, which would particularly affect the identification and genotyping of rare variants. Despite our best efforts to equalise coverage across every base pair in every individual, we acknowledge that important loss-of-function variants may exist in ST8SIA2 which were not identified in the current study, due to low coverage of some genomic regions in particular individuals.

We identified 72 variants in our study that were not observed in the 1 kG Caucasian population (41 novel and 26 known SNPs in ST8SIA2, 5 novel SNPs in NCAM1). Unfortunately, we were not able to compare the rates of novel variant detection between the bipolar and 1 kG populations due to inherent differences in coverage and platform – which is a major limitation of the study. However, as the majority of novel SNPs were observed in a small number of cases, we can deduce that they cannot individually account for the increased risk attached to the specific risk haplotype previously identified [20]. Indeed, the majority of variants may have arisen on a spectrum of haplotypes – not restricted specifically to the risk haplotype previously associated with illness – consistent with expectations of allelic heterogeneity, which confounds identification of disease-causing variants.

We observed no mutations in the protein-coding regions of either ST8SIA2 or NCAM1 that would reduce functionality of the enzyme or its substrate protein. Therefore, we focused our analysis on prediction of function in non-protein-coding gene regions. The enrichment of SNPs at conserved positions in the genome within DHSP sites in our data is consistent with rare and deleterious variants being present in ST8SIA2. Although measures of cross-species conservation are useful in highlighting variant sites that are most likely to be functional, not all functional residues are expected to have been conserved across species evolution. Despite the wealth of data available on DHSs, most conserved sites have no obvious predictable functional effect. This may be indicative of the limitations of the currently available cellular models, or current knowledge of modes of gene regulation. In the absence of any clear predicted function of these conserved sequences, forward genetics techniques could be used to examine the impact of these variants.

In addition to rare and highly deleterious variants, common variants that have a more subtle effect on ST8SIA2 function are also plausible candidates for contributing to increased risk of mental illness by regulating the timing and level of gene expression. The dose-responsive and functionally redundant nature of ST8SIA2 and ST8SIA4 is well characterised [7], and it is possible that genetic interaction between these two genes determines risk, with the dosage of each gene contributing. Sites within the 3′UTR are particularly attractive candidates for such an effect, given the important role of microRNAs in fine-tuning gene expression in the brain [54]. MicroRNA binding sites are involved in exquisite control of gene dosage during differentiation, a process that if disrupted may lead to a subtle neurodevelopmental phenotype that may precede bipolar disorder symptomatology. We found two SNPs within a conserved region at the start of the 3′UTR that are strong candidates for altering the post-transcriptional regulation of ST8SIA2. Common variants that impart an effect on sites of transcriptional regulation are also candidates for imparting subtle gene-regulatory effects. A recent finding implicated a combinatorial effect of multiple enhancers on gene expression and, in turn, genetic risk of disease [32]. Intriguingly, in that study, ST8SIA2 was identified as a bipolar GWAS locus in which variants in enhancers functionally impact expression in neuronal precursor cells [32]. Although the GWAS result [24] was not obtained in the same population as our cohort, this finding suggests that the combination of multiple DNA variants across multiple enhancers is an important mode of regulation at the ST8SIA2 locus in bipolar disorder. This pattern of risk inheritance is entirely consistent with a combinatorial effect of many risk variants that lie on alternative haplotypic backgrounds, as we indicate.

In the current study, we identified two variants with putative functional effects that are nominally associated with disease and were found to segregate with the either the protective or the risk haplotype (rs11074064 and rs722645). Under the model of a combinatorial effect from multiple enhancers, these variants, in combination with other functional variants inherited with a specific haplotype, could confer the functional effect of the risk and protective haplotypes. As these variants would be expected to alter gene expression, further work combining close genetic characterisation at the ST8SIA2 locus with gene expression data could be used to effectively pinpoint the genetic determinants of this potentially critical mode of gene regulation.

While we have not identified obvious loss-of-function mutations that would affect the formation of PSA-NCAM in this cohort of patients with bipolar disorder, we have characterised the nucleotide variation that may be associated with the disease, and importantly, in the context of the specific risk and protective haplotypes previously identified. Further analysis of the functional impact of these candidate variants on the regulation of gene expression in vitro will likely elucidate further insights into the mechanisms through which this candidate gene increases risk to mental illness.

Materials and Methods

Ethics statement

This study was carried out in accordance with the latest version of the Declaration of Helsinki after specific approval of the study and consent procedures by the University of New South Wales Human Research Ethics Committee (HREC # 04144 and 10078). All participants provided informed written consent prior to inclusion in the study.

Cohort selection

Individuals were selected from the Australian bipolar disorder case-control cohort, which has been described previously [20]. Diagnosis was made by structured interview using the Diagnostic Interview for Genetic Studies (DIGS), using Research Diagnostic Criteria. To enrich for individuals who may have a genetic defect in ST8SIA2, we selected one individual affected with bipolar disorder type 1 (BPI) or type 2 (BPII) (n = 39 and n = 3 respectively) from 42 families who showed evidence of linkage to chromosome 15q25–26 [16], and six additional unrelated BPI cases from our case-control cohort to ensure equivalent representation of major haplotypes to those observed in population samples [20]. A total of 45 BPI cases (94%) and three BPII cases (6%) formed the selected sample.

454 PCR enrichment and Sequencing

We used long-range PCR to amplify genomic DNA from the target regions prior to 454 sequencing, dividing the ST8SIA2 locus into eight amplicons and covering the NCAM1 region of interest in a single amplicon.

Although we initially attempted to enrich the entire ST8SIA2 locus in this manner, we could not amplify a region surrounding the ST8SIA2 transcription start site, likely owing to a high GC content at the intended primer binding site (Figure 1). This GC-rich region was also found to be less accessible to next generation sequencing in an analysis of accessible regions reported in the 1 kG project (Figure 1). After several optimisation attempts (primer sequences: F-AAAGTTACAAAAGGAATGCAAATG, R-AGGAACCTCGACTGCCAAC; F-CCCACTACCACTTCCAGCTC, R-GCAGAGCAGTTTGGGAAGTC and F-CCCCCTTCAGAGATGATGAA, R-GCAGAGCAGTTTGGGAAGTC) this gap region was divided into two subregions. The 3′ subregion was successfully amplified by long-range PCR (ST8-2), but the 5′ subregion could not be amplified (primer sequences F-CCCCCTTCAGAGATGATGAA R-AAAGAAGCAGGGAGCACTCA). In order to effectively examine nucleotide variation within this region, we divided the unamplified region into 8 overlapping short-range PCR amplicons for subsequent bi-directional Sanger sequencing (Figure 1).

Final primer sequences (5′ to 3′) for long-range PCR were as follows: ST8-1F-AGCGTCTTTTCAGGAGGTGA; ST8-1R-CTACCCTGACCCAGCAACAT; ST8-2F-CGCTTCCTGCTCTCATTTTC; ST8-2R-GTTTCCTCCTTGCCATCGT; ST8-3F-TCCCAGTGAAGAGCACAGTC; ST8-3R-TTCCCATTGCCCTGAGTATC; ST8-4F-TAAGGTGGAGCTAGGGACCA; ST8-4R- GGAGACCAAATGCCTTGAAA; ST8-5F-TCCCATCTGTGATTCCATGA; ST8-5R-GGGCAGGATCTTCTTTCTCC; ST8-6F-AATGACAAGTGCCCCATAGC; ST8-6R-ACTGCAAACTCATGCTGACAA; ST8-7F-ACCCACTTTTTCTCCCCAGT; ST8-7R-TATGAGCACGGACAGCAATC; ST8-8F-ATCAGCCTTTCCCAACAGC; ST8-8R-CCACTCCCTCTACCCAATTTC; N1F-CATCTTTCAACCAACCAGCA; and N1R-CCTATTCCCCTCCAGCTACC.

Long-range PCR was performed using the Expand Long Range dNTPack kit (Roche, Mannheim, Germany), according to manufacturer's instructions. Briefly, 50 μl reaction volumes consisted of 10 μl of 5x Expand Long Range Buffer with 12.5 mM MgCl2, 2.5 μl of 10 mM PCR Nucleotide Mix, 0.3 μM of each primer, 0.7 μl of 5 U/μl Expand Long Range Enzyme Mix, with 500 ng genomic DNA. DMSO (5%) was added to PCR mix of all amplicons except amplicons ST8-2 and ST8-3. Thermal cycling conditions were performed on a PTC-240 DNA Engine Tetrad 2 Cycler (MJ Research, Bio-Rad Laboratories, Copenhagen, Denmark) provided by the manufacturer. The annealing temperature was 60°C for all amplicons, except amplicons ST8-3 and ST8-5 which were 59°C and 58°C, respectively.

PCR products from each amplicon were quantified by comparison to DNA molecular weight marker XV (Roche, Mannheim, Germany) after gel electrophoresis, using Bio-Rad Quality One software and Molecular Imager ChemiDoc XRS System (Bio-Rad, Hercules, CA, USA). Individual PCR products from the nine amplicons for each individual were purified using MultiScreen96-PCR filter plate (Millipore, Bedford, MA), and resuspended in 30 μl Milli-Q-grade water. Quality and concentration was checked by Nanodrop 1000 Spectrophotometer (Biolab, Musgrave, VIC, Australia). After concentration normalisation, all amplicons for each individual were pooled at an equimolar ratio (approx. 5 ng/μl each) to generate 1 μg of template for library preparation, which was performed by the Ramaciotti Centre for Gene Function Analysis (University of New South Wales, Sydney, Australia). Sequencing was performed on the 454 GS-FLX (454 Life Sciences, Branford, CT, USA).

Sanger Sequencing

Overlapping primers were designed across the target region (hg19: chr15[ratio]92933439-92938793 bp), based on overlapping amplicons of at least 500 bp, using the PCR suite [55] and Primer3 [56]. Primers used were 1L-ACGTGCGGAAGAGGATATTG; 1R-ATTTTCCAAGAAGCCCGAAG; 2L-GGCCTTTAATTTGCCACCAC; 2R-AGGTGACCCACGATTACAGC; 3L-TGGGTGCAAAGAAGGAAGAC; 3R-TGCTTGAAGGGTGGTTTAGG; 4L-GCGCATCCTGTACTTTCCTC; 4R-TTATCCCCTGGACCACTCAG; 5L-GGGCACAAGCAGTACTTTGG; 5R-GGCAGATGTTTCAGCAGTTG; 6L-CGCCCACAAACTGTAGTCAG; 6R-CGGAGGGTGGAGAGTACAGA; 7L-GAAATCGGGTAAATAGCTGCTC; 7R-GCGAGGTCGGGAGAAA; 8F-GGAGACGGAGACATTTCGGG; 8R-GCTCTTGGGGAAACCGAGAA. Genomic DNA was amplified using Platinum Taq (Invitrogen, Carlsbad, CA) and purified using Multiscreen96-PCR plate. Sequencing reactions were performed in both forward and reverse directions for each amplicon, using the Big Dye Terminator v3.1 Cycle Sequencing Kit (Applied Biosystems, Foster City, CA, USA). The extension products were analysed on an ABI 3730 by the Ramaciotti Centre (University of New South Wales, Australia). The Phred/Phrap/Polyphred software suite [41] was used for base calling, sequence alignment and polymorphism identification. Consed [42] was used for viewing sequences and traces, and to confirm polymorphisms in each individual by their presence on both forward and reverse sequences.

Data Acquisition and Processing

dbSNP 135 vcf files were downloaded from the NCBI ftp server (ftp://ftp.ncbi.nih.gov/snp/organisms/human_9606/VCF/), and the subset of variants within the sequenced region was selected using SelectVariants and CombineVariants modules of the GATK. The 1 kG reference panel in BEAGLE format was downloaded from http://bochet.gcc.biostat.washington.edu/beagle/1000_Genomes.phase1_release_v3/. A subset of reference panel genotypes, including Caucasian European individuals (CEU and GBR populations) and sequence target region markers, was selected from the reference panel using BEAGLE utilities (http://faculty.washington.edu/browning/beagle_utilities/utilities.html). Genotype data for 1 kG Caucasian European Individuals was downloaded using the 1000 Genomes Data Browser.

454 Analysis: Mapping

BWA mapping: SFF files containing chromatograms were converted to fastq files using sff2fastq (https://github.com/indraniel/sff2fastq). Fastq sequences were demultiplexed using appropriate RLMID barcodes for each sample and the barcode sequences were subsequently trimmed using the FASTX toolkit (http://hannonlab.cshl.edu/fastx_toolkit/). Reads were aligned to the b37 Human Genome Reference Consortium human reference sequence, downloaded from Broad Institute (http://www.broadinstitute.org/gatk/guide/article.php?id=1213). Alignment was performed using BWA [36] in Smith-Waterman mode resulting in sam files [57]. Readgroup meta information, necessary for combined analysis of individual samples in GATK [40], was added and files were merged into a single bam file using picardtools (http://picard.sourceforge.net/).

Refmapper: Reads were aligned to fasta files of chromosomes 11 and 15 of the hg19 UCSC human reference sequence using Refmapper (version 2.6, 454 Life Sciences, Branford, CT, USA).

454 Analysis: SNP calling and Genotyping

GATK: The GATK modules CountCovariates and TableRecalibration were used for basecalling recalibration. When applying default covariates, we found that base call quality scores did not vary predictably when stratified across machine cycle, which is in contrast to findings from whole-genome data [39], and is likely due to the smaller size of our targeted sequencing dataset. To avoid potentially introducing artefact, we excluded this covariate during base call quality recalibration within our analysis.

Following recalibration, UnifiedGenotyper was used for SNP calling and genotyping. SNPs were then filtered using VariantFiltration module based on Best Practice Variant Detection V4 recommended hard filters (http://www.broadinstitute.org/gatk/guide/article?id=1186), these were QD (Quality By Depth) <2.0 (GATK recommended) or 5.0 (our analysis), MQ (Mapping Quality) <40.0, FS (Fisher Strand Bias) >60.0, HaplotypeScore (Haplotype Score) >13.0, MQRankSum (Mapping Quality Rank-Sum p value) <−12.5, ReadPosRankSum (Read Position Rank Sum Test) <−8.0. A custom Galaxy workflow incorporating these steps, which are based on the Best Practice Variant Detection v3 is available (https://main.g2.bx.psu.edu/u/a-shaw-neura/p/next-generation-resources).

Refmapper: SNP calling was performed by merging the HCDiffs.txt output from Refmapper for the 48 individually mapped individuals.

Genotyping and Allele Frequency Comparison

Illumina genotype data was prepared using PLINK-Seq (http://atgu.mgh.harvard.edu/plinkseq/) from two association studies of Australian bipolar disorder case-control cohort samples: one using an Illumina BeadArray [20] and one using an Illumina 660 W chip [43]. Overlapping sample genotypes from the two datasets were checked for exact concordance using GATK VariantEval module before merging into a single set of comparison genotypes. Concordance between Illumina genotypes and sequencing genotypes was examined using GATK VariantEval module.

1 kG Caucasian European comparison: SNPs identified in both datasets or unique to each dataset were selected using the SelectVariants GATK module. The sets of unique or overlapping SNPs were then compared using GATK VariantEval module. SNPs that were co-located with indels identified in 1 kG, and SNPs identified in Refmapper but not GATK, were not included in the AF analysis, as the AF at these sites could not be compared across the two datasets. A custom Galaxy workflow incorporating the above steps is available at https://main.g2.bx.psu.edu/u/a-shaw-neura/p/next-generation-resources.

BEAGLE Imputation and Phasing

Haplotypes were phased and genotype calls and quality scores were adjusted using GATK and BEAGLE, as per the example usage (http://gatkforums.broadinstitute.org/discussion/43/interface-with-beagle-software). Genotype likelihoods from 454 were prepared for BEAGLE input using ProduceBeagleInput GATK module. These likelihoods were then merged with genotype data from Sanger-sequenced region and Illumina genotypes, with Illumina genotypes taking priority over sequencing calls for discordant genotypes. BEAGLE was run using default parameters, with or without additional genotypes from a 1 kG reference panel as input. Phased genotypes were then converted into GATK vcf format using the BeagleOutputToVCF GATK module.

Bioinformatic analysis of predicted function

DNAse hypersensitivity site uniform peaks (DHSPs) from the ENCODE project [31], [45] were downloaded from UCSC Genome Browser [58]. Peaks from hESC (h1 and h7 lines) and neuronal (SK-N-SH RA+, SK-N-MC and BE(2)C lines) were each downloaded into the Galaxy suite. For foetal brain, DHSPs from nine library preparations from six donors were downloaded from the Roadmap Epigenomics Project [59] browser. The peaks were clustered in Galaxy to eliminate peaks that were only present in a single library preparation. Genome Evolutionary Rate Profiling (GERP) scores were downloaded from the UCSC Genome Browser for each variant site, and GERP conserved elements were obtained from http://mendel.stanford.edu/SidowLab/downloads/gerp/. Splice site predictions were performed using two web-based tools: NNSPLICE (http://www.fruitfly.org/about/index.html) [60], and Human Splicing Finder V2.4.1 (http://www.umd.be/HSF/) [61]. Variant Effect Predictor (http://asia.ensembl.org/info/docs/variation/vep/index.html) [62] was used to identify putative functional variants within the coding regions. PreSTIGE Enhancer regions for ST8SIA2 were downloaded from http://genetics.case.edu/prestige/.

Sequenom Validation

SNPs were selected for Sequenom validation based on the following criteria: Criteria 1 ([GERP >2] and [in DHS or PreSTIGE Enhancer or Frequency difference ≥ 0.05 between 454 bipolar and 1 K genomes CEU]); and Criteria 2 ([in PreSTIGE Enhancer or DHS] and [Frequency difference ≥ 0.05 or novel]). Assays for selected SNPs were designed and performed at the Australian Genome Research Facility (University of Queensland, Brisbane, Australia).

Manual Genotyping and Association Analysis

For three SNPs, manual genotyping by RFLP or Sanger sequencing was performed in Australian bipolar disorder case-control cohort individuals. Merging of manual genotypes and Sequenom genotypes and subsequent statistical analyses were performed using PLINK software. SNPs with MAF <0.05 were removed prior to association analysis. Association analysis for bipolar disorder was performed using a broad disease model where individuals diagnosed with BPI, SZMA, BPII or MDD were all classified as affected (n = 215).

Supporting Information

Figure S1

Comparison of genotype quality score and minor allele frequency to identify potentially erroneous SNPs. All points represent raw SNPs that were novel (not present in dbSNP 135). Allele frequency (AF) vs. Quality by Depth score (QD). Hollow diamonds – SNPs that failed filtering according to GATK recommended filtering, Black circles – SNPs that passed filtering. Points in dashed oval have low QD and high AF, suggesting they are not bona-fide SNPs.

(DOCX)

Figure S2

Linkage disequilibrium plot of SNPs verified by direct genotyping in larger case control cohort (n = 207 bipolar cases, 160 controls). Six of the seven SNPs in block 2 had a p>0.1. Haplotype analysis of the block 2 SNPs showed a trend for association (omnibus p = 0.077), with the major haplotype (red) more frequent in affected cases (F_A) than unaffected controls (F_U), representing a risk haplotype. The second most common haplotype (green) was more frequent in unaffected controls (F_U) than affected cases (F_A), representing a protective haplotype.

(DOCX)

Table S1

Complete list of genetic variation identified in ST8SIA2 . The location of all identified variants on chromosome 15 ST8SIA2 region (hg19 build), with transcribed SNPs indicated with asterisks. SNPs identified exclusively on the risk haplotype are shown (0 = not exclusive, 1 = present on risk and other haplotypes, 2 = present on risk haplotype only, ND = not determined). SNPs identified exclusively on the protective haplotype are shown (0 = not exclusive, 1 = present on protective and other haplotypes, 2 = present on protective haplotype only, ND = not determined). The set from which the SNP was observed is given (1 = GATK & Refmapper; 2 = GATK only; 3a = Refmapper only; 3b = Refmapper only & filtered in GATK; 4 = Sanger). Co-localisation with DNase I hypersensitivity site peaks (DHSPs; neuronal = 1; hESC = 2; foetal brain = 3) are given. Genomic Evolutionary Rate Profiling (GERP) scores are provided for each variant that is within a GERP-conserved element. The nature of the polymorphism in each cohort is given (with minor allele listed first). The minor allele frequency (mAF) of each variant in the 47 bipolar cases and 174 Caucasian individuals (CEU and GBR) from the 1000 Genomes Project (1 kG) are shown separately. The frequency difference (freqDIFF) was calculated relative to 1 kG allele frequency, and those with p values <0.1 indicated with an asterisk. aFor each polymorphism, the minor allele is listed first. v SNPs verified by direct genotyping are indicated.

(DOCX)

Note S1

GATK SNP Filtering Threshold Note S1: GATK SNP Filtering Threshold.

(DOCX)

Note S2

Removal of data in regions of PCR allelic-bias.

(DOCX)

Acknowledgments

We thank all patients who participated in this study. We also thank Mr Jason Koval for preparing the sequencing libraries and performing the 454 GS-FLX sequencing at the Ramaciotti Centre for Gene Function Analysis (University of New South Wales, Sydney, Australia). We thank Erica McAuley for providing the SNP genotypes for genotype calling accuracy studies, and Anna Heath and Kerrie Pierce for their assistance with DNA handling. DNA samples were prepared by Genetic Repositories Australia (Sydney, Australia), supported by NHMRC Enabling Grant 401184. Carol Dobson-Stone provided feedback during manuscript preparation. Finally, we thank Tamara McDonald from the Australian Genome Research Facility for her assistance in the design of Sequenom assays and generation of genotypes for the validation study.

Funding Statement

This work was supported by the Australian National Medical and Health Research Council (project grant 630574, and program grants 510135 and 1037196) and by the Schizophrenia Research Institute, utilizing infrastructure funding from NSW Health. Genetic Repositories Australia is supported by an Australian National Health and Medical Research Council (Grant # 401184). The authors used Garvan Galaxy, funded by a Cancer Institute NSW grant (11/REG/1–10). Mr Yash Tiwari was supported by an Australian Postgraduate Award. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.

References

1. Finne J (1982) Occurrence of unique polysialosyl carbohydrate units in glycoproteins of developing brain. J Biol Chem 257: 11966–11970. [Abstract] [Google Scholar]
2. Eckhardt M, Mühlenhoff M, Bethe A (1995) Molecular characterization of eukaryotic polysialyltransferase-1. Nature 373: 715–718. [Abstract] [Google Scholar]
3. Nakayama J, Fukuda M (1995) Expression cloning of a human polysialyltransferase that forms the polysialylated neural cell adhesion molecule present in embryonic brain. Proc Natl Acad Sci 92: 7031–7035. [Europe PMC free article] [Abstract] [Google Scholar]
4. Rutishauser U (2008) Polysialic acid in the plasticity of the developing and adult vertebrate nervous system. Nat Rev Neurosci 9: 26–35 10.1038/nrn2285 [Abstract] [CrossRef] [Google Scholar]
5. Mühlenhoff M, Rollenhagen M, Werneburg S, Gerardy-Schahn R, Hildebrandt H (2013) Polysialic Acid: Versatile Modification of NCAM, SynCAM 1 and Neuropilin-2. Neurochem Res 38: 1134–1143 10.1007/s11064-013-0979-2 [Abstract] [CrossRef] [Google Scholar]
6. Fukuda M (1996) A Human Polysialyltransferase Directs in Vitro Synthesis of Polysialic Acid. J Biol Chem 271: 1829–1832 10.1074/jbc.271.4.1829 [Abstract] [CrossRef] [Google Scholar]
7. Hildebrandt H, Mühlenhoff M, Oltmann-Norden I, Röckle I, Burkhardt H, et al. (2009) Imbalance of neural cell adhesion molecule and polysialyltransferase alleles causes defective brain connectivity. Brain 132: 2831–2838 10.1093/brain/awp117 [Abstract] [CrossRef] [Google Scholar]
8. Kröcher T, Malinovskaja K, Jürgenson M, Aonurm-Helm A, Zharkovskaya T, et al... (2013) Schizophrenia-like phenotype of polysialyltransferase ST8SIA2-deficient mice. Brain Struct Funct. 10.1007/s00429-013-0638-z. [Abstract]
9. Barbeau D, Liang JJ, Robitalille Y, Quirion R, Srivastava LK (1995) Decreased expression of the embryonic form of the neural cell adhesion molecule in schizophrenic brains. Proc Natl Acad Sci 92: 2785–2789. [Europe PMC free article] [Abstract] [Google Scholar]
10. Varea E, Guirado R, Gilabert-Juan J, Martí U, Castillo-Gomez E, et al. (2012) Expression of PSA-NCAM and synaptic proteins in the amygdala of psychiatric disorder patients. J Psychiatr Res 46: 189–197 10.1016/j.jpsychires.2011.10.011 [Abstract] [CrossRef] [Google Scholar]
11. Maheu ME, Davoli MA, Turecki G, Mechawar N (2013) Amygdalar expression of proteins associated with neuroplasticity in major depression and suicide. J Psychiatr Res 47: 384–390 10.1016/j.jpsychires.2012.11.013 [Abstract] [CrossRef] [Google Scholar]
12. Mikkonen M, Soininen H, Kälviänen R, Tapiola T, Ylinen A, et al. (1998) Remodeling of neuronal circuitries in human temporal lobe epilepsy: increased expression of highly polysialylated neural cell adhesion molecule in the hippocampus and the entorhinal cortex. Ann Neurol 44: 923–934 10.1002/ana.410440611 [Abstract] [CrossRef] [Google Scholar]
13. Castillo-Gómez E, Gómez-Climent MA, Varea E, Guirado R, Blasco-Ibáñez JM, et al. (2008) Dopamine acting through D2 receptors modulates the expression of PSA-NCAM, a molecule related to neuronal structural plasticity, in the medial prefrontal cortex of adult rats. Exp Neurol 214: 97–111 10.1016/j.expneurol.2008.07.018 [Abstract] [CrossRef] [Google Scholar]
14. Maćkowiak M, Dudys D, Chocyk A, Wedzony K (2009) Repeated risperidone treatment increases the expression of NCAM and PSA-NCAM protein in the rat medial prefrontal cortex. Eur Neuropsychopharmacol 19: 125–137 10.1016/j.euroneuro.2008.10.001 [Abstract] [CrossRef] [Google Scholar]
15. Frasca A, Fumagalli F, Ter Horst J, Racagni G, Murphy KJ, et al. (2008) Olanzapine, but not haloperidol, enhances PSA-NCAM immunoreactivity in rat prefrontal cortex. Int J Neuropsychopharmacol 11: 591–595 10.1017/S1461145708009061 [Abstract] [CrossRef] [Google Scholar]
16. McAuley EZ, Blair IP, Liu Z, Fullerton JM, Scimone A, et al. (2009) A genome screen of 35 bipolar affective disorder pedigrees provides significant evidence for a susceptibility locus on chromosome 15q25–26. Mol Psychiatry 14: 492–500 10.1038/sj.mp.4002146 [Abstract] [CrossRef] [Google Scholar]
17. Vazza G, Bertolin C, Scudellaro E, Vettori A, Boaretto F, et al. (2007) Genome-wide scan supports the existence of a susceptibility locus for schizophrenia and bipolar disorder on chromosome 15q26. Mol Psychiatry 12: 87–93 10.1038/sj.mp.4001895 [Abstract] [CrossRef] [Google Scholar]
18. Maziade M, Roy M, Chagnon YC, Cliche D, Fournier J-P, et al. (2005) Shared and specific susceptibility loci for schizophrenia and bipolar disorder: a dense genome scan in Eastern Quebec families. Mol Psychiatry 10: 486–499 10.1038/sj.mp.4001594 [Abstract] [CrossRef] [Google Scholar]
19. Park N, Juo SH, Cheng R, Liu J, Loth JE, et al. (2004) Linkage analysis of psychosis in bipolar pedigrees suggests novel putative loci for bipolar disorder and shared susceptibility with schizophrenia. Mol Psychiatry 9: 1091–1099 10.1038/sj.mp.4001541 [Abstract] [CrossRef] [Google Scholar]
20. McAuley EZ, Scimone A, Tiwari Y, Agahi G, Mowry BJ, et al. (2012) Identification of Sialyltransferase 8B as a Generalized Susceptibility Gene for Psychotic and Mood Disorders on Chromosome 15q25–26. PLoS One 7: e38172 10.1371/journal.pone.0038172 [Europe PMC free article] [Abstract] [CrossRef] [Google Scholar]
21. Arai M, Yamada K, Toyota T, Obata N, Haga S, et al. (2006) Association between polymorphisms in the promoter region of the sialyltransferase 8B (SIAT8B) gene and schizophrenia. Biol Psychiatry 59: 652–659 10.1016/j.biopsych.2005.08.016 [Abstract] [CrossRef] [Google Scholar]
22. Tao R, Li C, Zheng Y, Qin W, Zhang J, et al. (2007) Positive association between SIAT8B and schizophrenia in the Chinese Han population. Schizophr Res 90: 108–114 10.1016/j.schres.2006.09.029 [Abstract] [CrossRef] [Google Scholar]
23. Arai M, Itokawa M, Yamada K, Toyota T, Arai M, et al. (2004) Association of neural cell adhesion molecule 1 gene polymorphisms with bipolar affective disorder in Japanese individuals. Biol Psychiatry 55: 804–810 10.1016/j.biopsych.2004.01.009 [Abstract] [CrossRef] [Google Scholar]
24. Lee MTM, Chen CH, Lee CS, Chen CC, Chong MY, et al. (2011) Genome-wide association study of bipolar I disorder in the Han Chinese population. Mol Psychiatry 16: 548–556 10.1038/mp.2010.43 [Abstract] [CrossRef] [Google Scholar]
25. Anney R, Klei L, Pinto D, Regan R, Conroy J, et al. (2010) A genome-wide scan for common alleles affecting risk for autism. Hum Mol Genet 19: 4072–4082 10.1093/hmg/ddq307 [Europe PMC free article] [Abstract] [CrossRef] [Google Scholar]
26. Kamien B, Harraway J, Lundie B, Smallhorne L, Gibbs V, et al... (2013) Characterization of a 520 kb deletion on chromosome 15q26.1 including ST8SIA2 in a patient with behavioral disturbance, autism spectrum disorder, and epilepsy. Am J Med Genet: 1–7. 10.1002/ajmg.a.36345. [Abstract]
27. Harismendy O, Frazer K (2009) Method for improving sequence coverage uniformity of targeted genomic intervals amplified by LR-PCR using Illumina GA sequencing-by-synthesis technology. Biotechniques: 229–231. 10.2144/000113082. [Abstract]
28. The 1000 Genomes Project Consortium (2012) An integrated map of genetic variation from 1,092 human genomes. Nature 491: 56–65 10.1038/nature11632 [Europe PMC free article] [Abstract] [CrossRef] [Google Scholar]
29. America N, Accuracy G (2011) Genomes by the thousand. Nature 467: 1026–1027. [Google Scholar]
30. Dunham I, Kundaje A, Aldred SF, Collins PJ, Davis C a, et al. (2012) An integrated encyclopedia of DNA elements in the human genome. Nature 489: 57–74 10.1038/nature11247 [Europe PMC free article] [Abstract] [CrossRef] [Google Scholar]
31. Myers RM, Stamatoyannopoulos J, Snyder M, Dunham I, Hardison RC, et al. (2011) A user's guide to the encyclopedia of DNA elements (ENCODE). PLoS Biol 9: e1001046 10.1371/journal.pbio.1001046 [Europe PMC free article] [Abstract] [CrossRef] [Google Scholar]
32. Corradin O, Saiakhova A, Akhtar-Zaidi B, Myeroff L, Willis J, et al. (2014) Combinatorial effects of multiple enhancer variants in linkage disequilibrium dictate levels of gene expression to confer susceptibility to common traits. Genome Res 24: 1–13 10.1101/gr.164079.113 [Europe PMC free article] [Abstract] [CrossRef] [Google Scholar]
33. Close BE, Mendiratta SS, Geiger KM, Broom LJ, Ho L-L, et al. (2003) The minimal structural domains required for neural cell adhesion molecule polysialylation by PST/ST8Sia IV and STX/ST8Sia II. J Biol Chem 278: 30796–30805 10.1074/jbc.M305390200 [Abstract] [CrossRef] [Google Scholar]
34. Foley DA, Swartzentruber KG, Lavie A, Colley KJ (2010) Structure and mutagenesis of neural cell adhesion molecule domains: evidence for flexibility in the placement of polysialic acid attachment sites. J Biol Chem 285: 27360–27371 10.1074/jbc.M110.140038 [Europe PMC free article] [Abstract] [CrossRef] [Google Scholar]
35. Foley DA, Swartzentruber KG, Colley KJ (2009) Identification of sequences in the polysialyltransferases ST8Sia II and ST8Sia IV that are required for the protein-specific polysialylation of the neural cell adhesion molecule, NCAM. J Biol Chem 284: 15505–15516 10.1074/jbc.M809696200 [Europe PMC free article] [Abstract] [CrossRef] [Google Scholar]
36. Li H, Durbin R (2010) Fast and accurate long-read alignment with Burrows-Wheeler transform. Bioinformatics 26: 589–595 10.1093/bioinformatics/btp698 [Europe PMC free article] [Abstract] [CrossRef] [Google Scholar]
37. Harismendy O, Ng PC, Strausberg RL, Wang X, Stockwell TB, et al. (2009) Evaluation of next generation sequencing platforms for population targeted sequencing studies. Genome Biol 10: R32 10.1186/gb-2009-10-3-r32 [Europe PMC free article] [Abstract] [CrossRef] [Google Scholar]
38. Roy-Engel A, Salem A (2002) Active Alu element “A-tails”: size does matter. Genome Res 12: 1333–1344 10.1101/gr.384802 [Europe PMC free article] [Abstract] [CrossRef] [Google Scholar]
39. DePristo M a, Banks E, Poplin R, Garimella K V, Maguire JR, et al. (2011) A framework for variation discovery and genotyping using next-generation DNA sequencing data. Nat Genet 43: 491–498 10.1038/ng.806 [Europe PMC free article] [Abstract] [CrossRef] [Google Scholar]
40. McKenna A, Hanna M, Banks E, Sivachenko A, Cibulskis K, et al. (2010) The Genome Analysis Toolkit: a MapReduce framework for analyzing next-generation DNA sequencing data. Genome Res 20: 1297–1303 10.1101/gr.107524.110 [Europe PMC free article] [Abstract] [CrossRef] [Google Scholar]
41. Nickerson D a, Tobe VO, Taylor SL (1997) PolyPhred: automating the detection and genotyping of single nucleotide substitutions using fluorescence-based resequencing. Nucleic Acids Res 25: 2745–2751. [Europe PMC free article] [Abstract] [Google Scholar]
42. Gordon D, Abajian C, Green P (1998) Consed: a graphical tool for sequence finishing. Genome Res 8: 195–202. [Abstract] [Google Scholar]
43. Mühleisen T, Leber M, Schulze T, Strohmaier J, Degenhardt F, et al... (2014) Genome-wide association study reveals four new risk loci for bipolar disorder. Nat Commun. In press. [Abstract]
44. Halushka MK, Fan JB, Bentley K, Hsie L, Shen N, et al. (1999) Patterns of single-nucleotide polymorphisms in candidate genes for blood-pressure homeostasis. Nat Genet 22: 239–247 10.1038/10297 [Abstract] [CrossRef] [Google Scholar]
45. Rosenbloom KR, Sloan C a, Malladi VS, Dreszer TR, Learned K, et al. (2013) ENCODE data in the UCSC Genome Browser: year 5 update. Nucleic Acids Res 41: D56–63 10.1093/nar/gks1172 [Europe PMC free article] [Abstract] [CrossRef] [Google Scholar]
46. Darnell JC, Richter JD (2012) Cytoplasmic RNA-binding proteins and the control of complex brain function. Cold Spring Harb Perspect Biol 4: a012344 10.1101/cshperspect.a012344 [Europe PMC free article] [Abstract] [CrossRef] [Google Scholar]
47. Friedman RC, Farh KK-H, Burge CB, Bartel DP (2009) Most mammalian mRNAs are conserved targets of microRNAs. Genome Res 19: 92–105 10.1101/gr.082701.108 [Europe PMC free article] [Abstract] [CrossRef] [Google Scholar]
48. Sato C, Kitajima K (2013) Impact of structural aberrancy of polysialic acid and its synthetic enzyme ST8SIA2 in schizophrenia. Front Cell Neurosci 7: 61 10.3389/fncel.2013.00061 [Europe PMC free article] [Abstract] [CrossRef] [Google Scholar]
49. Chazal G, Durbec P, Jankovski A, Rougon G, Cremer H (2000) Consequences of neural cell adhesion molecule deficiency on cell migration in the rostral migratory stream of the mouse. J Neurosci 20: 1446–1457. [Europe PMC free article] [Abstract] [Google Scholar]
50. Hildebrandt H, Mühlenhoff M, Weinhold B, Gerardy-Schahn R (2007) Dissecting polysialic acid and NCAM functions in brain development. J Neurochem 103 Suppl: 56–6410.1111/j.1471-4159.2007.04716.x [Abstract] [CrossRef] [Google Scholar]
51. Isomura R, Kitajima K, Sato C (2011) Structural and functional impairments of polysialic acid by a mutated polysialyltransferase found in schizophrenia. J Biol Chem 286: 21535–21545 10.1074/jbc.M111.221143 [Europe PMC free article] [Abstract] [CrossRef] [Google Scholar]
52. Atz M, Rollins B, Vawter M (2007) NCAM1 association study of bipolar disorder and schizophrenia: polymorphisms and alternatively spliced isoforms lead to similarities and differences. Psychiatr Genet 17: 55–67. [Europe PMC free article] [Abstract] [Google Scholar]
53. Sullivan PF, Keefe RSE, Lange L a, Lange EM, Stroup TS, et al. (2007) NCAM1 and neurocognition in schizophrenia. Biol Psychiatry 61: 902–910 10.1016/j.biopsych.2006.07.036 [Abstract] [CrossRef] [Google Scholar]
54. Qureshi I a, Mehler MF (2012) Emerging roles of non-coding RNAs in brain evolution, development, plasticity and disease. Nat Rev Neurosci 13: 528–541 10.1038/nrn3234 [Europe PMC free article] [Abstract] [CrossRef] [Google Scholar]
55. Van Baren MJ, Heutink P (2004) The PCR suite. Bioinformatics 20: 591–593 10.1093/bioinformatics/btg473 [Abstract] [CrossRef] [Google Scholar]
56. Rozen S, Skaletsky HJ (2000) Primer3 on the WWW for general users and for biologist programmers. Bioinformatics Methods and Protocols. pp. 365–386. [Abstract]
57. Li H, Handsaker B, Wysoker A, Fennell T, Ruan J, et al. (2009) The Sequence Alignment/Map format and SAMtools. Bioinformatics 25: 2078–2079 10.1093/bioinformatics/btp352 [Europe PMC free article] [Abstract] [CrossRef] [Google Scholar]
58. Karolchik D, Barber GP, Casper J, Clawson H, Cline MS, et al. (2014) The UCSC Genome Browser database: 2014 update. Nucleic Acids Res 42: D764–70 10.1093/nar/gkt1168 [Europe PMC free article] [Abstract] [CrossRef] [Google Scholar]
59. Bernstein BE, Stamatoyannopoulos J a, Costello JF, Ren B, Milosavljevic A, et al. (2010) The NIH Roadmap Epigenomics Mapping Consortium. Nat Biotechnol 28: 1045–1048 10.1038/nbt1010-1045 [Europe PMC free article] [Abstract] [CrossRef] [Google Scholar]
60. Reese MG, Eeckman FH, Kulp D, Haussler D (1997) Improved splice site detection in Genie. J Comput Biol 4: 311–323 10.1089/cmb.1997.4.311 [Abstract] [CrossRef] [Google Scholar]
61. Desmet F-O, Hamroun D, Lalande M, Collod-Béroud G, Claustres M, et al. (2009) Human Splicing Finder: an online bioinformatics tool to predict splicing signals. Nucleic Acids Res 37: e67 10.1093/nar/gkp215 [Europe PMC free article] [Abstract] [CrossRef] [Google Scholar]
62. McLaren W, Pritchard B, Rios D, Chen Y, Flicek P, et al. (2010) Deriving the consequences of genomic variants with the Ensembl API and SNP Effect Predictor. Bioinformatics 26: 2069–2070 10.1093/bioinformatics/btq330 [Europe PMC free article] [Abstract] [CrossRef] [Google Scholar]
63. R Core Team (2012) R: A Language and Environment for Statistical Computing. Vienna, Austria: R Foundation for Statistical Computing.
Articles from PLOS ONE are provided here courtesy of PLOS

Full text links 

Read article at publisher's site: https://doi.org/10.1371/journal.pone.0092556

Read article for free, from open access legal sources, via Unpaywall: https://journals.plos.org/plosone/article/file?id=10.1371/journal.pone.0092556&type=printableUnpaywall

Citations & impact 

Impact metrics

25
Citations
Jump to Citations
7
Data citations
Jump to Data

Citations of article over time

Alternative metrics

Altmetric item for https://www.altmetric.com/details/2661698
Altmetric
Discover the attention surrounding your research
https://www.altmetric.com/details/2661698

Smart citations by scite.ai
Smart citations by scite.ai include citation statements extracted from the full text of the citing article. The number of the statements may be higher than the number of citations provided by EuropePMC if one paper cites another multiple times or lower if scite has not yet processed some of the citing articles.
Explore citation contexts and check if this article has been supported or disputed.
https://scite.ai/reports/10.1371/journal.pone.0092556

Supporting
Mentioning
Contrasting
0
17
0

Article citations

Go to all (25) article citations

Data 

Data behind the article

This data has been text mined from the article, or deposited into data resources.

BioStudies: supplemental material and supporting data

SNPs (Showing 54 of 54)

Similar Articles 

To arrive at the top five similar articles we use a word-weighted algorithm to compare words from the Title and Abstract of each citation.

Funding 

Funders who supported this work.

Medical Research Council