GWAS meta-analysis reveals novel loci and genetic correlates for general cognitive function: a report from the COGENT consortium.

This corrects the article DOI: 10.1038/mp.2016.244.

GWAS meta-analysis reveals novel loci and genetic correlates for general cognitive function: a report from the COGENT consortium.

The complex nature of human cognition has resulted in cognitive genomics lagging behind many other fields in terms of gene discovery using genome-wide association study (GWAS) methods. In an attempt to overcome these barriers, the current study utilized GWAS meta-analysis to examine the association of common genetic variation (~8M single-nucleotide polymorphisms (SNP) with minor allele frequency ⩾1%) to general cognitive function in a sample of 35 298 healthy individuals of European ancestry across 24 cohorts in the Cognitive Genomics Consortium (COGENT). In addition, we utilized individual SNP lookups and polygenic score analyses to identify genetic overlap with other relevant neurobehavioral phenotypes. Our primary GWAS meta-analysis identified two novel SNP loci (top SNPs: rs76114856 in the CENPO gene on chromosome 2 and rs6669072 near LOC105378853 on chromosome 1) associated with cognitive performance at the genome-wide significance level (P<5 × 10-8). Gene-based analysis identified an additional three Bonferroni-corrected significant loci at chromosomes 17q21.31, 17p13.1 and 1p13.3. Altogether, common variation across the genome resulted in a conservatively estimated SNP heritability of 21.5% (s.e.=0.01%) for general cognitive function. Integration with prior GWAS of cognitive performance and educational attainment yielded several additional significant loci. Finally, we found robust polygenic correlations between cognitive performance and educational attainment, several psychiatric disorders, birth length/weight and smoking behavior, as well as a novel genetic association to the personality trait of openness. These data provide new insight into the genetics of neurocognitive function with relevance to understanding the pathophysiology of neuropsychiatric illness.

GWA study data mining and independent replication identify cardiomyopathy-associated 5 (CMYA5) as a risk gene for schizophrenia.

We conducted data-mining analyses using the Clinical Antipsychotic Trials of Intervention Effectiveness (CATIE) and molecular genetics of schizophrenia genome-wide association study supported by the genetic association information network (MGS-GAIN) schizophrenia data sets and performed bioinformatic prioritization for all the markers with P-values ≤0.05 in both data sets. In this process, we found that in the CMYA5 gene, there were two non-synonymous markers, rs3828611 and rs10043986, showing nominal significance in both the CATIE and MGS-GAIN samples. In a combined analysis of both the CATIE and MGS-GAIN samples, rs4704591 was identified as the most significant marker in the gene. Linkage disequilibrium analyses indicated that these markers were in low LD (3 828 611-rs10043986, r(2)=0.008; rs10043986-rs4704591, r(2)=0.204). In addition, CMYA5 was reported to be physically interacting with the DTNBP1 gene, a promising candidate for schizophrenia, suggesting that CMYA5 may be involved in the same biological pathway and process. On the basis of this information, we performed replication studies for these three single-nucleotide polymorphisms. The rs3828611 was found to have conflicting results in our Irish samples and was dropped out without further investigation. The other two markers were verified in 23 other independent data sets. In a meta-analysis of all 23 replication samples (family samples, 912 families with 4160 subjects; case-control samples, 11 380 cases and 15 021 controls), we found that both markers are significantly associated with schizophrenia (rs10043986, odds ratio (OR)=1.11, 95% confidence interval (CI)=1.04-1.18, P=8.2 × 10(-4) and rs4704591, OR=1.07, 95% CI=1.03-1.11, P=3.0 × 10(-4)). The results were also significant for the 22 Caucasian replication samples (rs10043986, OR=1.11, 95% CI=1.03-1.17, P=0.0026 and rs4704591, OR=1.07, 95% CI=1.02-1.11, P=0.0015). Furthermore, haplotype conditioned analyses indicated that the association signals observed at these two markers are independent. On the basis of these results, we concluded that CMYA5 is associated with schizophrenia and further investigation of the gene is warranted.

Serious obstetric complications interact with hypoxia-regulated/vascular-expression genes to influence schizophrenia risk.

The etiology of schizophrenia is thought to include both epistasis and gene-environment interactions. We sought to test whether a set of schizophrenia candidate genes regulated by hypoxia or involved in vascular function in the brain (AKT1, BDNF, CAPON, CHRNA7, COMT, DTNBP1, GAD1, GRM3, NOTCH4, NRG1, PRODH, RGS4, TNF-alpha) interacted with serious obstetric complications to influence risk for schizophrenia. A family-based study of transmission disequilibrium was conducted in 116 trios. Twenty-nine probands had at least one serious obstetric complication (OC) using the McNeil-Sjostrom Scale, and many of the OCs reported were associated with the potential for fetal hypoxia. Analyses were conducted using conditional logistic regression and a likelihood ratio test (LRT) between nested models was performed to assess significance. Of the 13 genes examined, four (AKT1 (three SNPs), BDNF (two SNPs), DTNBP1 (one SNP) and GRM3 (one SNP)) showed significant evidence for gene-by-environment interaction (LRT P-values ranged from 0.011 to 0.037). Although our sample size was modest and the power to detect interactions was limited, we report significant evidence for genes involved in neurovascular function or regulated by hypoxia interacting with the presence of serious obstetric complications to increase risk for schizophrenia.

Allelic variation in GAD1 (GAD67) is associated with schizophrenia and influences cortical function and gene expression.

Cortical GABAergic dysfunction has been implicated as a key component of the pathophysiology of schizophrenia and decreased expression of the gamma-aminobutyric acid (GABA) synthetic enzyme glutamic acid decarboxylase 67 (GAD(67)), encoded by GAD1, is found in schizophrenic post-mortem brain. We report evidence of distorted transmission of single-nucleotide polymorphism (SNP) alleles in two independent schizophrenia family-based samples. In both samples, allelic association was dependent on the gender of the affected offspring, and in the Clinical Brain Disorders Branch/National Institute of Mental Health (CBDB/NIMH) sample it was also dependent on catechol-O-methyltransferase (COMT) Val158Met genotype. Quantitative transmission disequilibrium test analyses revealed that variation in GAD1 influenced multiple domains of cognition, including declarative memory, attention and working memory. A 5' flanking SNP affecting cognition in the families was also associated in unrelated healthy individuals with inefficient BOLD functional magnetic resonance imaging activation of dorsal prefrontal cortex (PFC) during a working memory task, a physiologic phenotype associated with schizophrenia and altered cortical inhibition. In addition, a SNP in the 5' untranslated (and predicted promoter) region that also influenced cognition was associated with decreased expression of GAD1 mRNA in the PFC of schizophrenic brain. Finally, we observed evidence of statistical epistasis between two SNPs in COMT and SNPs in GAD1, suggesting a potential biological synergism leading to increased risk. These coincident results implicate GAD1 in the etiology of schizophrenia and suggest that the mechanism involves altered cortical GABA inhibitory activity, perhaps modulated by dopaminergic function.

Significant correlation in linkage signals from genome-wide scans of schizophrenia and schizotypy.

Prior family and adoption studies have suggested a genetic relationship between schizophrenia and schizotypy. However, this has never been verified using linkage methods. We therefore attempted to test for a correlation in linkage signals from genome-wide scans of schizophrenia and schizotypy. The Irish study of high-density schizophrenia families comprises 270 families with at least two members with schizophrenia or poor-outcome schizoaffective disorder (n=637). Non-psychotic relatives were assessed using the structured interview for schizotypy (n=746). A 10-cM multipoint, non-parametric, autosomal genome-wide scan of schizophrenia was performed in Merlin. A scan of a quantitative trait comprising ratings of DSM-III-R criteria for schizotypal personality disorder in non-psychotic relatives was also performed. Schizotypy logarithm of the odds (LOD) scores were regressed onto schizophrenia LOD scores at all loci, with adjustment for spatial autocorrelation. To assess empirical significance, this was also carried out using 1000 null scans of schizotypy. The number of jointly linked loci in the real data was compared to distribution of jointly linked loci in the null scans. No markers were suggestively linked to schizotypy based on strict Lander-Kruglyak criteria. Schizotypy LODs predicted schizophrenia LODs above chance expectation genome wide (empirical P=0.04). Two and four loci yielded nonparametric LOD (NPLs) >1.0 and >0.75, respectively, for both schizophrenia and schizotypy (genome-wide empirical P=0.04 and 0.02, respectively). These results suggest that at least a subset of schizophrenia susceptibility genes also affects schizotypy in non-psychotic relatives. Power may therefore be increased in molecular genetic studies of schizophrenia if they incorporate measures of schizotypy in non-psychotic relatives.

Dysbindin (DTNBP1, 6p22.3) is associated with childhood-onset psychosis and endophenotypes measured by the Premorbid Adjustment Scale (PAS).

Straub et al. (2002) recently identified the 6p22.3 gene dysbindin (DTNBP1) through positional cloning as a schizophrenia susceptibility gene. We studied a rare cohort of 102 children with onset of psychosis before age 13. Standardized ratings of early development, medication response, neuropsychological and cognitive performance, premorbid dysfunction and clinical follow-up were obtained. Fourteen SNPs were genotyped in the gene DTNBP1. Family-based pairwise and haplotype transmission disequilibrium test (TDT) analysis with the clinical phenotype, and quantitative transmission disequilibrium test (QTDT) explored endophenotype relationships. One SNP was associated with diagnosis (TDT p=.01). The QTDT analyses showed several significant relationships. Four adjacent SNPs were associated (p values=.0009-.003) with poor premorbid functioning. These findings support the hypothesis that this and other schizophrenia susceptibility genes contribute to early neurodevelopmental impairment.

GAD1 (2q31.1), which encodes glutamic acid decarboxylase (GAD67), is associated with childhood-onset schizophrenia and cortical gray matter volume loss.

Postmortem brain studies have shown deficits in the cortical gamma-aminobutyric acid (GABA) system in schizophrenic individuals. Expression studies have shown a decrease in the major GABA-synthesizing enzyme (glutamic acid decarboxylase (GAD67) mRNA levels in neurons in dorsolateral prefrontal cortex in schizophrenics relative to controls. In the present study, SNPs in and around the GAD1 gene, which encodes the protein GAD67, were tested on a rare, severely ill group of children and adolescents with childhood-onset schizophrenia (COS) (n=72), in a family-based association analysis. Compared to adult-onset samples, the COS sample has evidence for more salient familial, and perhaps genetic, risk factors for schizophrenia, as well as evidence for frontal cortical hypofunction, and greater decline in cortical gray matter volume on anatomic brain MRI scans during adolescence. We performed family-based TDT and haplotype association analyses of the clinical phenotype, as well as association analyses with endophenotypes using the QTDT program. Three adjacent SNPs in the 5' upstream region of GAD1 showed a positive pairwise association with illness in these families (P=0.022-0.057). Significant transmission distortion of 4-SNP haplotypes was also observed (P=0.003-0.008). Quantitative trait TDT analyses showed an intriguing association between several SNPs and increased rate of frontal gray matter loss. These observations, when taken together with the positive results reported recently in two independent adult-onset schizophrenia pedigree samples, suggest that the gene encoding GAD67 may be a common risk factor for schizophrenia.

Clinical features of psychotic disorders and polymorphisms in HT2A, DRD2, DRD4, SLC6A3 (DAT1), and BDNF: a family based association study.

Schizophrenia is clinically heterogeneous and multidimensional, but it is not known whether this is due to etiological heterogeneity. Previous studies have not consistently reported association between any specific polymorphisms and clinical features of schizophrenia, and have primarily used case-control designs. We tested for the presence of association between clinical features and polymorphisms in the genes for the serotonin 2A receptor (HT2A), dopamine receptor types 2 and 4, dopamine transporter (SLC6A3), and brain-derived neurotrophic factor (BDNF). Two hundred seventy pedigrees were ascertained on the basis of having two or more members with schizophrenia or poor outcome schizoaffective disorder. Diagnoses were made using a structured interview based on the SCID. All patients were rated on the major symptoms of schizophrenia scale (MSSS), integrating clinical and course features throughout the course of illness. Factor analysis revealed positive, negative, and affective symptom factors. The program QTDT was used to implement a family-based test of association for quantitative traits, controlling for age and sex. We found suggestive evidence of association between the His452Tyr polymorphism in HT2A and affective symptoms (P = 0.02), the 172-bp allele of BDNF and negative symptoms (P = 0.04), and the 480-bp allele in SLC6A3 (= DAT1) and negative symptoms (P = 0.04). As total of 19 alleles were tested, we cannot rule out false positives. However, given prior evidence of involvement of the proteins encoded by these genes in psychopathology, our results suggest that more attention should be focused on the impact of these alleles on clinical features of schizophrenia.

Expression analysis of neuregulin-1 in the dorsolateral prefrontal cortex in schizophrenia.

Genetic linkage and association have implicated neuregulin-1 (NRG-1) as a schizophrenia susceptibility gene. We measured mRNA expression levels of the three major isoforms of NRG-1 (ie type I, type II, and type III) in the postmortem dorsolateral prefrontal cortex (DLPFC) from matched patients and controls using real-time quantitative RT-PCR. Expression levels of three internal controls-GAPDH, cyclophilin, and beta-actin-were unchanged in schizophrenia, and there were no changes in the absolute levels of the NRG-1 isoforms. However, type I expression normalized by GAPDH levels was significantly increased in schizophrenia DLPFC (by 23%) and positively correlated with antipsychotic medication dosage. Type II/type I and type II/type III ratios were significantly decreased (18 and 23% respectively). There was no effect on the NRG-1 mRNA levels of genotype at two SNPs previously associated with schizophrenia, suggesting that these alleles are not functionally responsible for abnormal NRG-1 expression patterns in patients. Subtle abnormalities in the expression patterns of NRG-1 mRNA isoforms in DLPFC may be associated with schizophrenia.

Genome-wide scans of three independent sets of 90 Irish multiplex schizophrenia families and follow-up of selected regions in all families provides evidence for multiple susceptibility genes.

From our linkage study of Irish families with a high density of schizophrenia, we have previously reported evidence for susceptibility genes in regions 5q21-31, 6p24-21, 8p22-21, and 10p15-p11. In this report, we describe the cumulative results from independent genome scans of three a priori random subsets of 90 families each, and from multipoint analysis of all 270 families in ten regions. Of these ten regions, three (13q32, 18p11-q11, and 18q22-23) did not generate scores above the empirical baseline pairwise scan results, and one (6q13-26) generated a weak signal. Six other regions produced more positive pairwise and multipoint results. They showed the following maximum multipoint H-LOD (heterogeneity LOD) and NPL scores: 2p14-13: 0.89 (P = 0.06) and 2.08 (P = 0.02), 4q24-32: 1.84 (P = 0.007) and 1.67 (P = 0.03), 5q21-31: 2.88 (P= 0.0007), and 2.65 (P = 0.002), 6p25-24: 2.13 (P = 0.005) and 3.59 (P = 0.0005), 6p23: 2.42 (P = 0.001) and 3.07 (P = 0.001), 8p22-21: 1.57 (P = 0.01) and 2.56 (P = 0.005), 10p15-11: 2.04 (P = 0.005) and 1.78 (P = 0.03). The degree of 'internal replication' across subsets differed, with 5q, 6p, and 8p being most consistent and 2p and 10p being least consistent. On 6p, the data suggested the presence of two susceptibility genes, in 6p25-24 and 6p23-22. Very few families were positive on more than one region, and little correlation between regions was evident, suggesting substantial locus heterogeneity. The levels of statistical significance were modest, as expected from loci contributing to complex traits. However, our internal replications, when considered along with the positive results obtained in multiple other samples, suggests that most of these six regions are likely to contain genes that influence liability to schizophrenia.

Association of the tryptophan hydroxylase gene with smoking initiation but not progression to nicotine dependence.

We investigated the association between two markers in the seventh intron of the tryptophan hydroxylase gene (TPH C218A and C779A) in a population-based case control study of 780 genotyped subjects. As in prior studies, the two markers were in strong linkage disequilibrium. The phenotypes we studied were smoking initiation and progression to nicotine dependence. Allele, genotype, and estimated haplotype frequencies for each marker were highly significantly different for smoking initiation (P < 0.0004 for each comparison) and were nonsignificant for progression to nicotine dependence. An empirical test suggested that the positive results were unlikely to have resulted from population stratification. Our results are similar to those of Lerman et al. [2001: Am J Med Genet (Neuropsychiatr Genet) 105:000-000] in associating these TPH markers with a construct related to smoking initiation but dissimilar in the variable implicated. If these results replicate in other samples, the serotonergic system may be involved in the etiology of smoking initiation given the rate-limiting role of TPH in the biosynthesis of serotonin.

Effect of COMT Val108/158 Met genotype on frontal lobe function and risk for schizophrenia.

Abnormalities of prefrontal cortical function are prominent features of schizophrenia and have been associated with genetic risk, suggesting that susceptibility genes for schizophrenia may impact on the molecular mechanisms of prefrontal function. A potential susceptibility mechanism involves regulation of prefrontal dopamine, which modulates the response of prefrontal neurons during working memory. We examined the relationship of a common functional polymorphism (Val(108/158) Met) in the catechol-O-methyltransferase (COMT) gene, which accounts for a 4-fold variation in enzyme activity and dopamine catabolism, with both prefrontally mediated cognition and prefrontal cortical physiology. In 175 patients with schizophrenia, 219 unaffected siblings, and 55 controls, COMT genotype was related in allele dosage fashion to performance on the Wisconsin Card Sorting Test of executive cognition and explained 4% of variance (P = 0.001) in frequency of perseverative errors. Consistent with other evidence that dopamine enhances prefrontal neuronal function, the load of the low-activity Met allele predicted enhanced cognitive performance. We then examined the effect of COMT genotype on prefrontal physiology during a working memory task in three separate subgroups (n = 11-16) assayed with functional MRI. Met allele load consistently predicted a more efficient physiological response in prefrontal cortex. Finally, in a family-based association analysis of 104 trios, we found a significant increase in transmission of the Val allele to the schizophrenic offspring. These data suggest that the COMT Val allele, because it increases prefrontal dopamine catabolism, impairs prefrontal cognition and physiology, and by this mechanism slightly increases risk for schizophrenia.

Analysis of epistasis in linked regions in the Irish study of high-density schizophrenia families.

Epistasis may be important in the etiology of schizophrenia. Analysis of epistasis has been important in the positional cloning of a gene involved in the etiology of type II diabetes mellitus. We investigated the importance of epistasis among six linked regions in 268 multiplex pedigrees in the Irish Study of High-Density Schizophrenia Families (ISHDSF) by computing pairwise correlations between nonparametric linkage scores for narrow, intermediate, and broad diagnostic definitions. The linked regions were on chromosomes 2, 4, 5, 6, 8, and 10. No correlation reached our a priori level of statistical significance. Using this statistical approach, we did not find evidence of important epistatic effects among these six regions in the ISHDSF.

An association study of DRD5 with smoking initiation and progression to nicotine dependence.

A large body of genetic epidemiological data strongly implicate genetic factors in the etiology of smoking behavior. Polymorphisms of genes in the dopaminergic system are plausible functional candidate genes and a linkage and an association study suggested that the type 5 dopamine receptor gene (DRD5) may be etiologically involved. We investigated the association of four DRD5 polymorphisms with smoking initiation and progression to nicotine dependence in a population-based sample of over 900 subjects. For smoking initiation, there was no significant association with the four DRD5 markers we studied; however, maximum likelihood analyses suggested the presence of a haplotype protective against smoking initiation. For progression to nicotine dependence, there were no strongly significant associations with the four DRD5 markers or for the estimated haplotypes. These data are not consistent with a strong etiological role for DRD5 in the etiology of these complex smoking behaviors.

Haplotypes of four novel single nucleotide polymorphisms in the nicotinic acetylcholine receptor beta2-subunit (CHRNB2) gene show no association with smoking initiation or nicotine dependence.

Several types of evidence, including experiments with mice that lack the nicotinic acetylcholine receptor beta2-subunit gene (CHRNB2), have suggested that a beta2-containing nicotinic receptor is necessary for at least some of the reinforcing properties of nicotine. However, sequence variations in CHRNB2 have not been reported, and its role in influencing human smoking behavior and nicotine dependence is not known. We screened most of the introns and exons and found five novel single nucleotide polymorphisms (SNPs). We tested four of these SNPs in three large, carefully selected samples: nonsmokers (n = 317) and regular smokers low levels of nicotine dependence (ND, n = 238), or smokers with high-ND (n = 317). None of the four polymorphisms we tested, nor their estimated haplotypes, were associated with smoking initiation or progression to nicotine dependence.