Pathway analysis of a genome-wide association study in schizophrenia

Objective

The aim of this study was to identify the candidate single nucleotide polymorphisms (SNPs) and candidate mechanisms that contribute to schizophrenia susceptibility and to generate a SNP to gene to pathway hypothesis using an analytical pathway-based approach.

Methods

We used schizophrenia GWAS data of the genotypes of 660,259 SNPs in 1378 controls and 1351 cases of European descent after quality control filtering. ICSNPathway (Identify candidate Causal SNPs and Pathways) analysis was applied to the schizophrenia GWAS dataset. The first stage involved the pre-selection of candidate SNPs by linkage disequilibrium analysis and the functional SNP annotation of the most significant SNPs found. The second stage involved the annotation of biological mechanisms for the pre-selected candidate SNPs using improved-gene set enrichment analysis.

Results

ICSNPathway analysis identified fifteen candidate SNPs, ten candidate pathways, and nine hypothetical biological mechanisms. The most strongly associated potential pathways were as follows. First, rs1644731 and rs1644730 to RDH8 to estrogen biosynthetic process (p < 0.001, FDR < 0.001). The genes involved in this pathway are RDH8 and HSD3B1 (p < 0.05). All-trans-retinol dehydrogenase (RDH8) is a visual cycle enzyme that reduces all-trans-retinal to all-trans-retinol in the presence of NADPH. The chemical reactions and pathways involved result in the formation of estrogens, which are C18 steroid hormones that can stimulate the development of female sexual characteristics. Second, rs1146031 to ACVR1 to mesoderm formation and activin binding (p < 0.001, FDR = 0.032, 0.034). Two of 15 candidate genes are known genes associated with schizophrenia: KCNQ2 and APOL2. One of the 10 candidate pathways, estrogen biosynthetic process, is known to be associated with schizophrenia (p < 0.001, FDR < 0.001). However, 13 of candidate genes (RDH8, ACVR1, PSMD9, KCNAB1, SLC17A3, ARCN1, COG7, STAB2, LRPAP1, STAB1, CXCL16, COL4A4, EXOSC3) and 9 of candidate pathways were novel.

Conclusion

By applying ICSNPathway analysis to schizophrenia GWAS data, we identified candidate SNPs, genes like KCNQ2 and APOL2 and pathways involving the estrogen biosynthetic process may contribute to schizophrenia susceptibility. Further analyses are needed to validate the results of this analysis.     

Whole-Genome Pathway Analysis on 132,497 Individuals Identifies Novel Gene-Sets Associated with Body Mass Index

Whole genome pathway analysis is a powerful tool for the exploration of the combined effects of gene-sets within biological pathways. This study applied Interval Based Enrichment Analysis (INRICH) to perform whole-genome pathway analysis of body-mass index (BMI). We used a discovery set composed of summary statistics from a meta-analysis of 123,865 subjects performed by the GIANT Consortium, and an independent sample of 8,632 subjects to assess replication of significant pathways. We examined SNPs within nominally significant pathways using linear mixed models to estimate their contribution to overall BMI heritability. Six pathways replicated as having significant enrichment for association after correcting for multiple testing, including the previously unknown relationships between BMI and the Reactome regulation of ornithine decarboxylase pathway, the KEGG lysosome pathway, and the Reactome stabilization of P53 pathway. Two non-overlapping sets of genes emerged from the six significant pathways. The clustering of shared genes based on previously identified protein-protein interactions listed in PubMed and OMIM supported the relatively independent biological effects of these two gene-sets. We estimate that the SNPs located in examined pathways explain ∼20% of the heritability for BMI that is tagged by common SNPs (3.35% of the 16.93% total).     

Harnessing expression data to identify novel candidate genes in polycystic ovary syndrome

Novel pathways in polycystic ovary syndrome (PCOS) are being identified in gene expression studies in PCOS tissues; such pathways may contain key genes in disease etiology. Previous expression studies identified both dickkopf homolog 1 (DKK1) and DnaJ (Hsp40) homolog, subfamily B, member 1 (DNAJB1) as differentially expressed in PCOS tissue, implicating them as candidates for PCOS susceptibility. To test this, we genotyped a discovery cohort of 335 PCOS cases and 198 healthy controls for three DKK1 single nucleotide polymorphisms (SNPs) and four DNAJB1 SNPs and a replication cohort of 396 PCOS cases and 306 healthy controls for 1 DKK1 SNP and 1 DNAJB1 SNP. SNPs and haplotypes were determined and tested for association with PCOS and component phenotypes. We found that no single nucleotide polymorphisms were associated with PCOS risk; however, the major allele of rs1569198 from DKK1 was associated with increased total testosterone (discovery cohort P = 0.0035) and dehydroepiandrosterone sulfate (replication cohort P = 0.05). Minor allele carriers at rs3962158 from DNAJB1 had increased fasting insulin (discovery cohort P = 0.003), increased HOMA-IR (discovery cohort P = 0.006; replication cohort P = 0.036), and increased HOMA-%B (discovery cohort P = 0.004). Carriers of haplotype 2 at DNAJB1 also had increased fasting insulin, HOMA-IR, and HOMA-%B. These findings suggest that genetic variation in DKK1 and DNAJB1 may have a role in the hyperandrogenic and metabolic dysfunction of PCOS, respectively. Our results also demonstrate the utility of gene expression data as a source of novel candidate genes in PCOS, a complex and still incompletely defined disease, for which alternative methods of gene identification are needed.     

Network-assisted investigation of combined causal signals from genome-wide association studies in schizophrenia

With the recent success of genome-wide association studies (GWAS), a wealth of association data has been accomplished for more than 200 complex diseases/traits, proposing a strong demand for data integration and interpretation. A combinatory analysis of multiple GWAS datasets, or an integrative analysis of GWAS data and other high-throughput data, has been particularly promising. In this study, we proposed an integrative analysis framework of multiple GWAS datasets by overlaying association signals onto the protein-protein interaction network, and demonstrated it using schizophrenia datasets. Building on a dense module search algorithm, we first searched for significantly enriched subnetworks for schizophrenia in each single GWAS dataset and then implemented a discovery-evaluation strategy to identify module genes with consistent association signals. We validated the module genes in an independent dataset, and also examined them through meta-analysis of the related SNPs using multiple GWAS datasets. As a result, we identified 205 module genes with a joint effect significantly associated with schizophrenia; these module genes included a number of well-studied candidate genes such as DISC1, GNA12, GNA13, GNAI1, GPR17, and GRIN2B. Further functional analysis suggested these genes are involved in neuronal related processes. Additionally, meta-analysis found that 18 SNPs in 9 module genes had P _meta<1×10⁻⁴, including the gene HLA-DQA1 located in the MHC region on chromosome 6, which was reported in previous studies using the largest cohort of schizophrenia patients to date. These results demonstrated our bi-directional network-based strategy is efficient for identifying disease-associated genes with modest signals in GWAS datasets. This approach can be applied to any other complex diseases/traits where multiple GWAS datasets are available.     

Comprehensive SNP Scan of DNA Repair and DNA Damage Response Genes Reveal Multiple Susceptibility Loci Conferring Risk to Tobacco Associated Leukoplakia and Oral Cancer

Polymorphic variants of DNA repair and damage response genes play major role in carcinogenesis. These variants are suspected as predisposition factors to Oral Squamous Cell Carcinoma (OSCC). For identification of susceptible variants affecting OSCC development in Indian population, the “maximally informative” method of SNP selection from HapMap data to non-HapMap populations was applied. Three hundred twenty-five SNPs from 11 key genes involved in double strand break repair, mismatch repair and DNA damage response pathways were genotyped on a total of 373 OSCC, 253 leukoplakia and 535 unrelated control individuals. The significantly associated SNPs were validated in an additional cohort of 144 OSCC patients and 160 controls. The rs12515548 of MSH3 showed significant association with OSCC both in the discovery and validation phases (discovery P-value: 1.43E-05, replication P-value: 4.84E-03). Two SNPs (rs12360870 of MRE11A, P-value: 2.37E-07 and rs7003908 of PRKDC, P-value: 7.99E-05) were found to be significantly associated only with leukoplakia. Stratification of subjects based on amount of tobacco consumption identified SNPs that were associated with either high or low tobacco exposed group. The study reveals a synergism between associated SNPs and lifestyle factors in predisposition to OSCC and leukoplakia.     

A case control association study and cognitive function analysis of neuropilin and tolloid-like 1 gene and schizophrenia in the Japanese population

Background

Using a knock-out mouse model, it was shown that NETO1 is a critical component of the NMDAR complex, and that loss of Neto1 leads to impaired hippocampal long term potentiation and hippocampal-dependent learning and memory. Moreover, hemizygosity of NETO1 was shown to be associated with autistic-like behavior in humans.

Purpose of the Research

We examined the association between schizophrenia and the neuropilin and tolloid-like 1 gene (NETO1). First, we selected eight single nucleotide polymorphisms (SNPs) within the NETO1 locus, based on the Japanese schizophrenia genome wide association study (JGWAS) results and previously conducted association studies. These SNPs were genotyped in the replication sample comprised of 963 schizophrenic patients and 919 healthy controls. We also examined the effect of associated SNPs on scores in the Continuous Performance Test and the Wisconsin Card Sorting Test Keio version (schizophrenic patients 107, healthy controls 104).

Results

There were no significant allele-wise and haplotype-wise associations in the replication analysis after Bonferroni correction. However, in meta-analysis (JGWAS and replication dataset) three association signals were observed (rs17795324: p = 0.028, rs8098760: p = 0.017, rs17086492: p = 0.003). These SNPs were followed up but we could not detect the allele-specific effect on cognitive performance measured by the Continuous performance test (CPT) and Wisconsin Card Sorting test (WCST).

Major Conclusions

We did not detect evidence for the association of NETO1 with schizophrenia in the Japanese population. Common variants within the NETO1 locus may not increase the genetic risk for schizophrenia in the Japanese population. Additionally, common variants investigated in the current study did not affect cognitive performance, as measured by the CPT and WCST.     

Genome-wide pathway analysis of genome-wide association studies on systemic lupus erythematosus and rheumatoid arthritis

The aim of this study was to explore candidate single nucleotide polymorphisms (SNPs) and candidate mechanisms of systemic lupus erythematosus (SLE) and rheumatoid arthritis (RA). Two SLE genome-wide association studies (GWASs) datasets were included in this study. Meta-analysis was conducted using 737,984 SNPs in 1,527 SLE cases and 3,421 controls of European ancestry, and 4,429 SNPs that met a threshold of p?<?0.01 in a Korean RA GWAS dataset was used. ICSNPathway (identify candidate causal SNPs and pathways) analysis was applied to the meta-analysis results of the SLE GWAS datasets, and a RA GWAS dataset. The most significant result of SLE GWAS meta-analysis concerned rs2051549 in the human leukocyte antigen (HLA) region (p?=?3.36E?22). In the non-HLA region, meta-analysis identified 6 SNPs associated with SLE with genome-wide significance (STAT4, TNPO3, BLK, FAM167A, and IRF5). ICSNPathway identified five candidate causal SNPs and 13 candidate causal pathways. This pathway-based analysis provides three hypotheses of the biological mechanism involved. First, rs8084 and rs7192?→?HLA-DRA?→?bystander B cell activation. Second, rs1800629?→?TNF?→?cytokine network. Third, rs1150752 and rs185819?→?TNXB?→?collagen metabolic process. ICSNPathway analysis identified three candidate causal non-HLA SNPs and four candidate causal pathways involving the PADI4, MTR, PADI2, and TPH2 genes of RA. We identified five candidate SNPs and thirteen pathways, involving bystander B cell activation, cytokine network, and collagen metabolic processing, which may contribute to SLE susceptibility, and we revealed candidate causal non-HLA SNPs, genes, and pathways of RA.     

A Genome-Wide Association Study Suggests Novel Loci Associated with a Schizophrenia-Related Brain-Based Phenotype

Patients with schizophrenia and their siblings typically show subtle changes of brain structures, such as a reduction of hippocampal volume. Hippocampal volume is heritable, may explain a variety of cognitive symptoms of schizophrenia and is thus considered an intermediate phenotype for this mental illness. The aim of our analyses was to identify single-nucleotide polymorphisms (SNP) related to hippocampal volume without making prior assumptions about possible candidate genes. In this study, we combined genetics, imaging and neuropsychological data obtained from the Mind Clinical Imaging Consortium study of schizophrenia (n = 328). A total of 743,591 SNPs were tested for association with hippocampal volume in a genome-wide association study. Gene expression profiles of human hippocampal tissue were investigated for gene regions of significantly associated SNPs. None of the genetic markers reached genome-wide significance. However, six highly correlated SNPs (rs4808611, rs35686037, rs12982178, rs1042178, rs10406920, rs8170) on chromosome 19p13.11, located within or in close proximity to the genes NR2F6, USHBP1, and BABAM1, as well as four SNPs in three other genomic regions (chromosome 1, 2 and 10) had p-values between 6.75×10⁻⁶ and 8.3×10⁻⁷. Using existing data of a very recently published GWAS of hippocampal volume and additional data of a multicentre study in a large cohort of adolescents of European ancestry, we found supporting evidence for our results. Furthermore, allelic differences in rs4808611 and rs8170 were highly associated with differential mRNA expression in the cis-acting region. Associations with memory functioning indicate a possible functional importance of the identified risk variants. Our findings provide new insights into the genetic architecture of a brain structure closely linked to schizophrenia. In silico replication, mRNA expression and cognitive data provide additional support for the relevance of our findings. Identification of causal variants and their functional effects may unveil yet unknown players in the neurodevelopment and the pathogenesis of neuropsychiatric disorders.     

Bivariate genome-wide association study suggests that the DARC gene influences lean body mass and age at menarche

Lean body mass (LBM) and age at menarche (AAM) are two important complex traits for human health. The aim of this study was to identify pleiotropic genes for both traits using a powerful bivariate genome-wide association study (GWAS). Two studies, a discovery study and a replication study, were performed. In the discovery study, 909622 single nucleotide polymorphisms (SNPs) were genotyped in 801 unrelated female Han Chinese subjects using the Affymetrix human genome-wide SNP array 6.0 platform. Then, a bivariate GWAS was performed to identify the SNPs that may be important for LBM and AAM. In the replication study, significant findings from the discovery study were validated in 1692 unrelated Caucasian female subjects. One SNP rs3027009 that was bivariately associated with left arm lean mass and AAM in the discovery samples (P=7.26×10^?6) and in the replication samples (P=0.005) was identified. The SNP is located at the upstream of DARC (Duffy antigen receptor for chemokines) gene, suggesting that DARC may play an important role in regulating the metabolisms of both LBM and AAM.     

Genome-wide association study for atopy and allergic rhinitis in a Singapore Chinese population

Allergic rhinitis (AR) is an atopic disease which affects about 600 million people worldwide and results from a complex interplay between genetic and environmental factors. However genetic association studies on known candidate genes yielded variable results. The aim of this study is to identify the genetic variants that influence predisposition towards allergic rhinitis in an ethnic Chinese population in Singapore using a genome-wide association study (GWAS) approach. A total of 4461 ethnic Chinese volunteers were recruited in Singapore and classified according to their allergic disease status. The GWAS included a discovery stage comparing 515 atopic cases (including 456 AR cases) and 486 non-allergic non-rhinitis (NANR) controls. The top SNPs were then validated in a replication cohort consisting of a separate 2323 atopic cases (including 676 AR cases) and 511 NANR controls. Two SNPs showed consistent association in both discovery and replication phases; MRPL4 SNP rs8111930 on 19q13.2 (OR = 0.69, P_combined = 4.46×10⁻⁰⁵) and BCAP SNP rs505010 on chromosome 10q24.1 (OR = 0.64, P_combined = 1.10×10⁻⁰⁴). In addition, we also replicated multiple associations within known candidates regions such as HLA-DQ and NPSR1 locus in the discovery phase. Our study suggests that MRPL4 and BCAP, key components of the HIF-1α and PI3K/Akt signaling pathways respectively, are two novel candidate genes for atopy and allergic rhinitis. Further study on these molecules and their signaling pathways would help in understanding of the pathogenesis of allergic rhinitis and identification of targets for new therapeutic intervention.     

A Major Histocompatibility Class I Locus Contributes to Multiple Sclerosis Susceptibility Independently from HLA-DRB1*15:01

Background

In Northern European descended populations, genetic susceptibility for multiple sclerosis (MS) is associated with alleles of the human leukocyte antigen (HLA) Class II gene DRB1. Whether other major histocompatibility complex (MHC) genes contribute to MS susceptibility is controversial.

Methodology/Principal Findings

A case control analysis was performed using 958 single nucleotide polymorphisms (SNPs) spanning the MHC assayed in two independent datasets. The discovery dataset consisted of 1,018 cases and 1,795 controls and the replication dataset was composed of 1,343 cases and 1,379 controls. The most significantly MS-associated SNP in the discovery dataset was rs3135391, a Class II SNP known to tag the HLA-DRB1*15:01 allele, the primary MS susceptibility allele in the MHC (O.R. = 3.04, p<1×10⁻⁷⁸). To control for the effects of the HLA-DRB1*15:01 haplotype, case control analysis was performed adjusting for this HLA-DRB1*15:01 tagging SNP. After correction for multiple comparisons (false discovery rate = .05) 52 SNPs in the Class I, II and III regions were significantly associated with MS susceptibility in both datasets using the Cochran Armitage trend test. The discovery and replication datasets were merged and subjects carrying the HLA-DRB1*15:01 tagging SNP were excluded. Association tests showed that 48 of the 52 replicated SNPs retained significant associations with MS susceptibility independently of the HLA-DRB1*15:01 as defined by the tagging SNP. 20 Class I SNPs were associated with MS susceptibility with p-values ≤1×10⁻⁸. The most significantly associated SNP was rs4959039, a SNP in the downstream un-translated region of the non-classical HLA-G gene (Odds ratio 1.59, 95% CI 1.40, 1.81, p = 8.45×10⁻¹³) and is in linkage disequilibrium with several nearby SNPs. Logistic regression modeling showed that this SNP''s contribution to MS susceptibility was independent of the Class II and Class III SNPs identified in this screen.

Conclusions

A MHC Class I locus contributes to MS susceptibility independently of the HLA-DRB1*15:01 haplotype.     

Possible Associations of NTRK2 Polymorphisms with Antidepressant Treatment Outcome: Findings from an Extended Tag SNP Approach

Background

Data from clinical studies and results from animal models suggest an involvement of the neurotrophin system in the pathology of depression and antidepressant treatment response. Genetic variations within the genes coding for the brain-derived neurotrophic factor (BDNF) and its key receptor Trkb (NTRK2) may therefore influence the response to antidepressant treatment.

Methods

We performed a single and multi-marker association study with antidepressant treatment outcome in 398 depressed Caucasian inpatients participating in the Munich Antidepressant Response Signature (MARS) project. Two Caucasian replication samples (N = 249 and N = 247) were investigated, resulting in a total number of 894 patients. 18 tagging SNPs in the BDNF gene region and 64 tagging SNPs in the NTRK2 gene region were genotyped in the discovery sample; 16 nominally associated SNPs were tested in two replication samples.

Results

In the discovery analysis, 7 BDNF SNPs and 9 NTRK2 SNPs were nominally associated with treatment response. Three NTRK2 SNPs (rs10868223, rs1659412 and rs11140778) also showed associations in at least one replication sample and in the combined sample with the same direction of effects (P_corr = .018, P_corr = .015 and P_corr = .004, respectively). We observed an across-gene BDNF-NTRK2 SNP interaction for rs4923468 and rs1387926. No robust interaction of associated SNPs was found in an analysis of BDNF serum protein levels as a predictor for treatment outcome in a subset of 93 patients.

Conclusions/Limitations

Although not all associations in the discovery analysis could be unambiguously replicated, the findings of the present study identified single nucleotide variations in the BDNF and NTRK2 genes that might be involved in antidepressant treatment outcome and that have not been previously reported in this context. These new variants need further validation in future association studies.     

A high-density genome-wide association screen of sporadic ALS in US veterans

Following reports of an increased incidence of amyotrophic lateral sclerosis (ALS) in U.S. veterans, we have conducted a high-density genome-wide association study (GWAS) of ALS outcome and survival time in a sample of U.S. veterans. We tested ∼1.3 million single nucleotide polymorphisms (SNPs) for association with ALS outcome in 442 incident Caucasian veteran cases diagnosed with definite or probable ALS and 348 Caucasian veteran controls. To increase power, we also included genotypes from 5909 publicly-available non-veteran controls in the analysis. In the survival analysis, we tested for association between SNPs and post-diagnosis survival time in 639 Caucasian veteran cases with definite or probable ALS. After this discovery phase, we performed follow-up genotyping of 299 SNPs in an independent replication sample of Caucasian veterans and non-veterans (ALS outcome: 183 cases and 961 controls; survival: 118 cases). Although no SNPs reached genome-wide significance in the discovery phase for either phenotype, three SNPs were statistically significant in the replication analysis of ALS outcome: rs6080539 (177 kb from PCSK2), rs7000234 (4 kb from ZNF704), and rs3113494 (13 kb from LOC100506746). Two SNPs located in genes that were implicated by previous GWA studies of ALS were marginally significant in the pooled analysis of discovery and replication samples: rs17174381 in DPP6 (p = 4.4×10⁻⁴) and rs6985069 near ELP3 (p = 4.8×10⁻⁴). Our results underscore the difficulty of identifying and convincingly replicating genetic associations with a rare and genetically heterogeneous disorder such as ALS, and suggest that common SNPs are unlikely to account for a substantial proportion of patients affected by this devastating disorder.     

Identification of candidate genes for dyslexia susceptibility on chromosome 18

Background

Six independent studies have identified linkage to chromosome 18 for developmental dyslexia or general reading ability. Until now, no candidate genes have been identified to explain this linkage. Here, we set out to identify the gene(s) conferring susceptibility by a two stage strategy of linkage and association analysis.

Methodology/Principal Findings

Linkage analysis: 264 UK families and 155 US families each containing at least one child diagnosed with dyslexia were genotyped with a dense set of microsatellite markers on chromosome 18. Association analysis: Using a discovery sample of 187 UK families, nearly 3000 SNPs were genotyped across the chromosome 18 dyslexia susceptibility candidate region. Following association analysis, the top ranking SNPs were then genotyped in the remaining samples. The linkage analysis revealed a broad signal that spans approximately 40 Mb from 18p11.2 to 18q12.2. Following the association analysis and subsequent replication attempts, we observed consistent association with the same SNPs in three genes; melanocortin 5 receptor (MC5R), dymeclin (DYM) and neural precursor cell expressed, developmentally down-regulated 4-like (NEDD4L).

Conclusions

Along with already published biological evidence, MC5R, DYM and NEDD4L make attractive candidates for dyslexia susceptibility genes. However, further replication and functional studies are still required.     

Genetic variants in toll-like receptors are not associated with rheumatoid arthritis susceptibility or anti-tumour necrosis factor treatment outcome

Background

Several studies point to a role of Toll-like receptors (TLRs) in the development of rheumatoid arthritis (RA). We investigated if genetic variants in TLR genes are associated with RA and response to tumour necrosis factor blocking (anti-TNF) medication.

Methodology and Principal Findings

22 single nucleotide polymorphisms (SNPs) in seven TLR genes were genotyped in a Dutch cohort consisting of 378 RA patients and 294 controls. Significantly associated variants were investigated in replication cohorts from The Netherlands, United Kingdom and Sweden (2877 RA patients and 2025 controls). 182 of the Dutch patients were treated with anti-TNF medication. Using these patients and a replication cohort (269 Swedish patients) we analysed if genetic variants in TLR genes were associated with anti-TNF outcome. In the discovery phase of the study we found a significant association of SNPs rs2072493 in TLR5 and rs3853839 in TLR7 with RA disease susceptibility. Meta-analysis of discovery and replication cohorts did not confirm these findings. SNP rs2072493 in TLR5 was associated with anti-TNF outcome in the Dutch but not in the Swedish population.

Conclusion

We conclude that genetic variants in TLRs do not play a major role in susceptibility for developing RA nor in anti-TNF treatment outcome in a Caucasian population.     

Significant associations of CHRNA2 and CHRNA6 with nicotine dependence in European American and African American populations

The direct physiological effects that promote nicotine dependence (ND) are mediated by nicotinic acetylcholine receptors (nAChRs). In line with the genetic and pharmacological basis of addiction, many previous studies have revealed significant associations between variants in the nAChR subunit genes and various measures of ND in different ethnic samples. In this study, we first examined the association of variants in nAChR subunits α2 (CHRNA2) and α6 (CHRNA6) genes on chromosome 8 with ND using a family sample consisting of 1,730 European Americans (EAs) from 495 families and 1,892 African Americans (AAs) from 424 families (defined as the discovery family sample). ND was assessed by two standard quantitative measures: smoking quantity (SQ) and the Fagerström Test for ND (FTND). We found nominal associations for all seven tested SNPs of the genes with at least one ND measure in the EA sample and for two SNPs in CHRNA2 in the AA sample. Of these, associations of SNPs rs3735757 with FTND (P = 0.0068) and rs2472553 with both ND measures (with a P value of 0.0043 and 0.00086 for SQ and FTND, respectively) continued to be significant in the EA sample even after correction for multiple tests. Further, we found several haplotypes that were significantly associated with ND in the EA sample in CHRNA6 and in the both EA and AA samples in CHRNA2. To confirm the associations of the two genes with ND, we conducted a replication study with an independent case–control sample from the SAGE study, which showed a significant association of the two genes with ND, although the significantly associated SNPs were not always the same in the two samples. Together, these findings indicate that both CHRNA2 and CHRNA6 play a significant role in the etiology of ND in AA and EA smokers. Further replication in additional independent samples is warranted.     

The Impact of Genome-Wide Supported Schizophrenia Risk Variants in the Neurogranin Gene on Brain Structure and Function

The neural mechanisms underlying genetic risk for schizophrenia, a highly heritable psychiatric condition, are still under investigation. New schizophrenia risk genes discovered through genome-wide association studies (GWAS), such as neurogranin (NRGN), can be used to identify these mechanisms. In this study we examined the association of two common NRGN risk single nucleotide polymorphisms (SNPs) with functional and structural brain-based intermediate phenotypes for schizophrenia. We obtained structural, functional MRI and genotype data of 92 schizophrenia patients and 114 healthy volunteers from the multisite Mind Clinical Imaging Consortium study. Two schizophrenia-associated NRGN SNPs (rs12807809 and rs12541) were tested for association with working memory-elicited dorsolateral prefrontal cortex (DLPFC) activity and surface-wide cortical thickness. NRGN rs12541 risk allele homozygotes (TT) displayed increased working memory-related activity in several brain regions, including the left DLPFC, left insula, left somatosensory cortex and the cingulate cortex, when compared to non-risk allele carriers. NRGN rs12807809 non-risk allele (C) carriers showed reduced cortical gray matter thickness compared to risk allele homozygotes (TT) in an area comprising the right pericalcarine gyrus, the right cuneus, and the right lingual gyrus. Our study highlights the effects of schizophrenia risk variants in the NRGN gene on functional and structural brain-based intermediate phenotypes for schizophrenia. These results support recent GWAS findings and further implicate NRGN in the pathophysiology of schizophrenia by suggesting that genetic NRGN risk variants contribute to subtle changes in neural functioning and anatomy that can be quantified with neuroimaging methods.