Common Variants on Xq28 Conferring Risk of Schizophrenia in Han Chinese

Schizophrenia is a highly heritable, severe psychiatric disorder affecting approximately 1% of the world population. A substantial portion of heritability is still unexplained and the pathophysiology of schizophrenia remains to be elucidated. To identify more schizophrenia susceptibility loci, we performed a genome-wide association study (GWAS) on 498 patients with schizophrenia and 2025 controls from the Han Chinese population, and a follow-up study on 1027 cases and 1005 controls. In the follow-up study, we included 384 single nucleotide polymorphisms (SNPs) which were selected from the top hits in our GWAS (130 SNPs) and from previously implicated loci for schizophrenia based on the SZGene database, NHGRI GWAS Catalog, copy number variation studies, GWAS meta-analysis results from the international Psychiatric Genomics Consortium (PGC) and candidate genes from plausible biological pathways (254 SNPs).Within the chromosomal region Xq28, SNP rs2269372 in RENBP achieved genome-wide significance with a combined P value of 3.98×10⁻⁸ (OR of allele A = 1.31). SNPs with suggestive P values were identified within 2 genes that have been previously implicated in schizophrenia, MECP2 (rs2734647, P _combined = 8.78×10⁻⁷, OR = 1.28; rs2239464, P _combined = 6.71×10⁻⁶, OR = 1.26) and ARHGAP4 (rs2269368, P _combined = 4.74×10⁻⁷, OR = 1.25). In addition, the patient sample in our follow-up study showed a significantly greater burden for pre-defined risk alleles based on the SNPs selected than the controls. This indicates the existence of schizophrenia susceptibility loci among the SNPs we selected. This also further supports multigenic inheritance in schizophrenia. Our findings identified a new schizophrenia susceptibility locus on Xq28, which harbor the genes RENBP, MECP2, and ARHGAP4.Key words: schizophrenia, genome-wide association study, Han Chinese, MECP2, ARHGAP4, RENBP     

Common variants at 1p36 are associated with superior frontal gyrus volume

The superior frontal gyrus (SFG), an area of the brain frequently found to have reduced gray matter in patients with schizophrenia, is involved in self-awareness and emotion, which are impaired in schizophrenia. However, no genome-wide association studies of SFG volume have investigated in patients with schizophrenia. To identify single-nucleotide polymorphisms (SNPs) associated with SFG volumes, we demonstrated a genome-wide association study (GWAS) of gray matter volumes in the right or left SFG of 158 patients with schizophrenia and 378 healthy subjects. We attempted to bioinformatically ascertain the potential effects of the top hit polymorphism on the expression levels of genes at the genome-wide region. We found associations between five variants on 1p36.12 and the right SFG volume at a widely used benchmark for genome-wide significance (P<5.0 × 10⁻⁸). The strongest association was observed at rs4654899, an intronic SNP in the eukaryotic translation initiation factor 4 gamma, 3 (EIF4G3) gene on 1p36.12 (P=7.5 × 10⁻⁹). No SNP with genome-wide significance was found in the volume of the left SFG (P>5.0 × 10⁻⁸); however, the rs4654899 polymorphism was identified as the locus with the second strongest association with the volume of the left SFG (P=1.5 × 10⁻⁶). In silico analyses revealed a proxy SNP of rs4654899 had effect on gene expression of two genes, HP1BP3 lying 3′ to EIF4G3 (P=7.8 × 10⁻⁶) and CAPN14 at 2p (P=6.3 × 10⁻⁶), which are expressed in moderate-to-high levels throughout the adult human SFG. These results contribute to understand genetic architecture of a brain structure possibly linked to the pathophysiology of schizophrenia.     

Meta-analysis of genome-wide association studies for panic disorder in the Japanese population

Panic disorder (PD) is a moderately heritable anxiety disorder whose pathogenesis is not well understood. Due to the lack of power in previous association studies, genes that are truly associated with PD might not be detected. In this study, we conducted a genome-wide association study (GWAS) in two independent data sets using the Affymetrix Mapping 500K Array or Genome-Wide Human SNP Array 6.0. We obtained imputed genotypes for each GWAS and performed a meta-analysis of two GWAS data sets (718 cases and 1717 controls). For follow-up, 12 single-nucleotide polymorphisms (SNPs) were tested in 329 cases and 861 controls. Gene ontology enrichment and candidate gene analyses were conducted using the GWAS or meta-analysis results. We also applied the polygenic score analysis to our two GWAS samples to test the hypothesis of polygenic components contributing to PD. Although genome-wide significant SNPs were not detected in either of the GWAS nor the meta-analysis, suggestive associations were observed in several loci such as BDKRB2 (P=1.3 × 10⁻⁵, odds ratio=1.31). Among previous candidate genes, supportive evidence for association of NPY5R with PD was obtained (gene-wise corrected P=6.4 × 10⁻⁴). Polygenic scores calculated from weakly associated SNPs (P<0.3 and 0.4) in the discovery sample were significantly associated with PD status in the target sample in both directions (sample I to sample II and vice versa) (P<0.05). Our findings suggest that large sets of common variants of small effects collectively account for risk of PD.     

Genome-wide association study identifies multiple novel loci associated with disease progression in subjects with mild cognitive impairment

Alzheimer''s disease (AD) is the leading cause of dementia among the elderly population; however, knowledge about genetic risk factors involved in disease progression is limited. We conducted a genome-wide association study (GWAS) using clinical decline as measured by changes in the Clinical Dementia Rating-sum of boxes as a quantitative trait to test for single-nucleotide polymorphisms (SNPs) that were associated with the rate of progression in 822 Caucasian subjects of amnestic mild cognitive impairment (MCI). There was no significant association with disease progress for any of the recently identified disease susceptibility variants in CLU, CR1, PICALM, BIN1, EPHA1, MS4A6A, MS4A4E or CD33 following multiple testing correction. We did, however, identify multiple novel loci that reached genome-wide significance at the 0.01 level. These top variants (rs7840202 at chr8 in UBR5: P=4.27 × 10⁻¹⁴; rs11637611 with a cluster of SNPs at chr15q23 close to the Tay–Sachs disease locus: P=1.07 × 10⁻¹⁵; and rs12752888 at chr1: P=3.08 × 10⁻¹¹) were also associated with a significant decline in cognition as well as the conversion of subjects with MCI to a diagnosis of AD. Taken together, these variants define approximately 16.6% of the MCI sub-population with a faster rate of decline independent of the other known disease risk factors. In addition to providing new insights into protein pathways that may be involved with the progress to AD in MCI subjects, these variants if further validated may enable the identification of a more homogeneous population of subjects at an earlier stage of disease for testing novel hypotheses and/or therapies in the clinical setting.     

Genome-wide pharmacogenomic study of citalopram-induced side effects in STAR*D

Affecting about 1 in 12 Americans annually, depression is a leading cause of the global disease burden. While a range of effective antidepressants are now available, failure and relapse rates remain substantial, with intolerable side effect burden the most commonly cited reason for discontinuation. Thus, understanding individual differences in susceptibility to antidepressant therapy side effects will be essential to optimize depression treatment. Here we perform genome-wide association studies (GWAS) to identify genetic variation influencing susceptibility to citalopram-induced side effects. The analysis sample consisted of 1762 depression patients, successfully genotyped for 421K single-nucleotide polymorphisms (SNPs), from the Sequenced Treatment Alternatives to Relieve Depression (STAR^*D) study. Outcomes included five indicators of citalopram side effects: general side effect burden, overall tolerability, sexual side effects, dizziness and vision/hearing side effects. Two SNPs met our genome-wide significance criterion (q<0.1), ensuring that, on average, only 10% of significant findings are false discoveries. In total, 12 additional SNPs demonstrated suggestive associations (q<0.5). The top finding was rs17135437, an intronic SNP within EMID2, mediating the effects of citalopram on vision/hearing side effects (P=3.27 × 10⁻⁸, q=0.026). The second genome-wide significant finding, representing a haplotype spanning ∼30 kb and eight genotyped SNPs in a gene desert on chromosome 13, was associated with general side effect burden (P=3.22 × 10⁻⁷, q=0.096). Suggestive findings were also found for SNPs at LAMA1, AOX2P, EGFLAM, FHIT and RTP2. Although our findings require replication and functional validation, this study demonstrates the potential of GWAS to discover genes and pathways that potentially mediate adverse effects of antidepressant medications.     

Genome-wide association studies in non-anxiety individuals identified novel risk loci for depression

BackgroundDepression is a debilitating mental disorder that often coexists with anxiety. The genetic mechanisms of depression and anxiety have considerable overlap, and studying depression in non-anxiety samples could help to discover novel gene. We assess the genetic variation of depression in non-anxiety samples, using genome-wide association studies (GWAS) and linkage disequilibrium score regression (LDSC).MethodsThe GWAS of depression score and self-reported depression were conducted using the UK Biobank samples, comprising 99,178 non-anxiety participants with anxiety score <5 and 86,503 non-anxiety participants without self-reported anxiety, respectively. Replication analysis was then performed using two large-scale GWAS summary data of depression from Psychiatric Genomics Consortium (PGC). LDSC was finally used to evaluate genetic correlations with 855 health-related traits based on the primary GWAS.ResultsTwo genome-wide significant loci for non-anxiety depression were identified: rs139702470 (p = 1.54 × 10⁻⁸, OR = 0.29) locate in PIEZO2, and rs6046722 (p = 2.52 × 10⁻⁸, OR = 1.09) locate in CFAP61. These associated genes were replicated in two GWAS of depression from PGC, such as rs1040582 (p_{replication GWAS1} = 0.02, p_{replication GWAS2} = 2.71 × 10⁻³) in CFAP61, and rs11661122 (p_{replication GWAS1} = 8.16 × 10⁻³, p_{replication GWAS2} = 8.08 × 10⁻³) in PIEZO2. LDSC identified 19 traits genetically associated with non-anxiety depression (p < 0.001), such as marital separation/divorce (rg = 0.45, SE = 0.15).ConclusionsOur findings provide novel clues for understanding of the complex genetic architecture of depression.     

Systematic Integration of Brain eQTL and GWAS Identifies ZNF323 as a Novel Schizophrenia Risk Gene and Suggests Recent Positive Selection Based on Compensatory Advantage on Pulmonary Function

Genome-wide association studies have identified multiple risk variants and loci that show robust association with schizophrenia. Nevertheless, it remains unclear how these variants confer risk to schizophrenia. In addition, the driving force that maintains the schizophrenia risk variants in human gene pool is poorly understood. To investigate whether expression-associated genetic variants contribute to schizophrenia susceptibility, we systematically integrated brain expression quantitative trait loci and genome-wide association data of schizophrenia using Sherlock, a Bayesian statistical framework. Our analyses identified ZNF323 as a schizophrenia risk gene (P = 2.22×10^–6). Subsequent analyses confirmed the association of the ZNF323 and its expression-associated single nucleotide polymorphism rs1150711 in independent samples (gene-expression: P = 1.40×10^–6; single-marker meta-analysis in the combined discovery and replication sample comprising 44123 individuals: P = 6.85×10⁻¹⁰). We found that the ZNF323 was significantly downregulated in hippocampus and frontal cortex of schizophrenia patients (P = .0038 and P = .0233, respectively). Evidence for pleiotropic effects was detected (association of rs1150711 with lung function and gene expression of ZNF323 in lung: P = 6.62×10^–5 and P = 9.00×10^–5, respectively) with the risk allele (T allele) for schizophrenia acting as protective allele for lung function. Subsequent population genetics analyses suggest that the risk allele (T) of rs1150711 might have undergone recent positive selection in human population. Our findings suggest that the ZNF323 is a schizophrenia susceptibility gene whose expression may influence schizophrenia risk. Our study also illustrates a possible mechanism for maintaining schizophrenia risk variants in the human gene pool.Key words: schizophrenia, ZNF323, association, eQTL, hippocampus, positive selection     

DEGS2 polymorphism associated with cognition in schizophrenia is associated with gene expression in brain

A genome-wide association study of cognitive deficits in patients with schizophrenia in Japan found association with a missense genetic variant (rs7157599, Asn8Ser) in the delta(4)-desaturase, sphingolipid 2 (DEGS2) gene. A replication analysis using Caucasian samples showed a directionally consistent trend for cognitive association of a proxy single-nucleotide polymorphism (SNP), rs3783332. Although the DEGS2 gene is expressed in human brain, it is unknown how DEGS2 expression varies during human life and whether it is affected by psychiatric disorders and genetic variants. To address these questions, we examined DEGS2 messenger RNA using next-generation sequencing in postmortem dorsolateral prefrontal cortical tissue from a total of 418 Caucasian samples including patients with schizophrenia, bipolar disorder and major depressive disorder. DEGS2 is expressed at very low levels prenatally and increases gradually from birth to adolescence and consistently expressed across adulthood. Rs3783332 genotype was significantly associated with the expression across all subjects (F_3,348=10.79, P=1.12 × 10⁻³), particularly in control subjects (F_1,87=13.14, P=4.86 × 10⁻⁴). Similar results were found with rs715799 genotype. The carriers of the risk-associated minor allele at both loci showed significantly lower expression compared with subjects homozygous for the non-risk major allele and this was a consistent finding across all diagnostic groups. DEGS2 expression showed no association with diagnostic status after correcting for multiple testing (P>0.05). Our findings demonstrate that a SNP showing genome-wide association study significant association with cognition in schizophrenia is also associated with regulation of DEGS2 expression, implicating a molecular mechanism for the clinical association.     

A Population-Based Study of Genetic Variation and Psychotic Experiences in Adolescents

Psychotic experiences are not uncommon in general population samples, but no studies have examined to what extent confirmed risk variants for schizophrenia are associated with such experiences. A total of 3483 children in a birth cohort study participated in semistructured interviews for psychotic experiences at ages 12 and 18. We examined whether (1) a composite measure of risk for schizophrenia conferred by common alleles (polygenic score) was associated with psychotic experiences, (2) variants with genome-wide evidence for association with schizophrenia were associated with psychotic experiences, and (3) we could identify genetic variants for psychotic experiences using a genome-wide association (GWA) approach. We found no evidence that a schizophrenia polygenic score, or variants showing genome-wide evidence of association with schizophrenia, were associated with adolescent psychotic experiences within the general population. In fact, individuals who had a higher number of risk alleles for genome-wide hits for schizophrenia showed a decreased risk of psychotic experiences. In the GWA study, no variants showed GWA for psychotic experiences, and there was no evidence that the strongest hits (P < 5 × 10⁻⁵) were enriched for variants associated with schizophrenia in large consortia. Although polygenic scores are weak tools for prediction of schizophrenia, they show strong evidence of association with this disorder. Our findings, however, lend little support to the hypothesis that psychotic experiences in population-based samples of adolescents share a comparable genetic architecture to schizophrenia, or that utilizing a broader and more common phenotype of psychotic experiences will be an efficient approach to increase understanding of the genetic etiology of schizophrenia.Key words: psychosis, schizophrenia, epidemiology, ALSPAC, GWAS, polygenic     

A common biological basis of obesity and nicotine addiction

Smoking influences body weight such that smokers weigh less than non-smokers and smoking cessation often leads to weight increase. The relationship between body weight and smoking is partly explained by the effect of nicotine on appetite and metabolism. However, the brain reward system is involved in the control of the intake of both food and tobacco. We evaluated the effect of single-nucleotide polymorphisms (SNPs) affecting body mass index (BMI) on smoking behavior, and tested the 32 SNPs identified in a meta-analysis for association with two smoking phenotypes, smoking initiation (SI) and the number of cigarettes smoked per day (CPD) in an Icelandic sample (N=34 216 smokers). Combined according to their effect on BMI, the SNPs correlate with both SI (r=0.019, P=0.00054) and CPD (r=0.032, P=8.0 × 10⁻⁷). These findings replicate in a second large data set (N=127 274, thereof 76 242 smokers) for both SI (P=1.2 × 10⁻⁵) and CPD (P=9.3 × 10⁻⁵). Notably, the variant most strongly associated with BMI (rs1558902-A in FTO) did not associate with smoking behavior. The association with smoking behavior is not due to the effect of the SNPs on BMI. Our results strongly point to a common biological basis of the regulation of our appetite for tobacco and food, and thus the vulnerability to nicotine addiction and obesity.     

Dissociation of accumulated genetic risk and disease severity in patients with schizophrenia

Genotype–phenotype correlations of common monogenic diseases revealed that the degree of deviation of mutant genes from wild-type structure and function often predicts disease onset and severity. In complex disorders such as schizophrenia, the overall genetic risk is still often >50% but genotype–phenotype relationships are unclear. Recent genome-wide association studies (GWAS) replicated a risk for several single-nucleotide polymorphisms (SNPs) regarding the endpoint diagnosis of schizophrenia. The biological relevance of these SNPs, however, for phenotypes or severity of schizophrenia has remained obscure. We hypothesized that the GWAS ‘top-10'' should as single markers, but even more so upon their accumulation, display associations with lead features of schizophrenia, namely positive and negative symptoms, cognitive deficits and neurological signs (including catatonia), and/or with age of onset of the disease prodrome as developmental readout and predictor of disease severity. For testing this hypothesis, we took an approach complementary to GWAS, and performed a phenotype-based genetic association study (PGAS). We utilized the to our knowledge worldwide largest phenotypical database of schizophrenic patients (n>1000), the GRAS (Göttingen Research Association for Schizophrenia) Data Collection. We found that the ‘top-10'' GWAS-identified risk SNPs neither as single markers nor when explored in the sense of a cumulative genetic risk, have any predictive value for disease onset or severity in the schizophrenic patients, as demonstrated across all core symptoms. We conclude that GWAS does not extract disease genes of general significance in schizophrenia, but may yield, on a hypothesis-free basis, candidate genes relevant for defining disease subgroups.     

Hostility in adolescents and adults: a genome-wide association study of the Young Finns

Hostility is a multidimensional personality trait with changing expression over the life course. We performed a genome-wide association study (GWAS) of the components of hostility in a population-based sample of Finnish men and women for whom a total of 2.5 million single-nucleotide polymorphisms (SNPs) were available through direct or in silico genotyping. Hostility dimensions (anger, cynicism and paranoia) were assessed at four time points over a 15-year interval (age range 15–30 years at phase 1 and 30–45 years at phase 4) in 982–1780 participants depending on the hostility measure. Few promising areas from chromosome 14 at 99 cM (top SNPs rs3783337, rs7158754, rs3783332, rs2181102, rs7159195, rs11160570, rs941898, P values <3.9 × 10⁻⁸ with nearest gene Enah/Vasp-like (EVL)) were found suggestively to be related to paranoia and from chromosome 7 at 86 cM (top SNPs rs802047, rs802028, rs802030, rs802026, rs802036, rs802025, rs802024, rs802032, rs802049, rs802051, P values <6.9 × 10⁻⁷ with nearest gene CROT (carnitine O-octanoyltransferase)) to cynicism, respectively. Some shared suggestive genetic influence for both paranoia and cynicism was also found from chromosome 17 at 2.8 cM (SNPs rs12936442, rs894664, rs6502671, rs7216028) and chromosome 22 at 43 cM (SNPs rs7510759, rs7510924, rs7290560), with nearest genes RAP1 GTPase activating protein 2 (RAP1GAP2) and KIAA1644, respectively. These suggestive associations did not replicate across all measurement times, which warrants further study on these SNPs in other populations.     

A genome-wide association study of antidepressant response in Koreans

We conducted a three-stage genome-wide association study (GWAS) of response to antidepressant drugs in an ethnically homogeneous sample of Korean patients in untreated episodes of nonpsychotic unipolar depression, mostly of mature onset. Strict quality control was maintained in case selection, diagnosis, verification of adherence and outcome assessments. Analyzed cases completed 6 weeks of treatment with adequate plasma drug concentrations. The overall successful completion rate was 85.5%. Four candidate single-nucleotide polymorphisms (SNPs) on three chromosomes were identified by genome-wide search in the discovery sample of 481 patients who received one of four allowed selective serotonin reuptake inhibitor (SSRI) antidepressant drugs (Stage 1). In a focused replication study of 230 SSRI-treated patients, two of these four SNP candidates were confirmed (Stage 2). Analysis of the Stage 1 and Stage 2 samples combined (n=711) revealed GWAS significance (P=1.60 × 10^-8) for these two SNP candidates, which were in perfect linkage disequilibrium. These two significant SNPs were confirmed also in a focused cross-replication study of 159 patients treated with the non-SSRI antidepressant drug mirtazapine (Stage 3). Analysis of the Stage 1, Stage 2 and Stage 3 samples combined (n=870) also revealed GWAS significance for these two SNPs, which was sustained after controlling for gender, age, number of previous episodes, age at onset and baseline severity (P=3.57 × 10^-8). For each SNP, the response rate decreased (odds ratio=0.31, 95% confidence interval: 0.20–0.47) as a function of the number of minor alleles (non-response alleles). The two SNPs significantly associated with antidepressant response are rs7785360 and rs12698828 of the AUTS2 gene, located on chromosome 7 in 7q11.22. This gene has multiple known linkages to human psychological functions and neurobehavioral disorders. Rigorous replication efforts in other ethnic populations are recommended.     

Follow-up of loci from the International Genomics of Alzheimer's Disease Project identifies TRIP4 as a novel susceptibility gene

To follow-up loci discovered by the International Genomics of Alzheimer''s Disease Project, we attempted independent replication of 19 single nucleotide polymorphisms (SNPs) in a large Spanish sample (Fundació ACE data set; 1808 patients and 2564 controls). Our results corroborate association with four SNPs located in the genes INPP5D, MEF2C, ZCWPW1 and FERMT2, respectively. Of these, ZCWPW1 was the only SNP to withstand correction for multiple testing (P=0.000655). Furthermore, we identify TRIP4 (rs74615166) as a novel genome-wide significant locus for Alzheimer''s disease risk (odds ratio=1.31; confidence interval 95% (1.19–1.44); P=9.74 × 10⁻⁹).     

The morphometry of left cuneus mediating the genetic regulation on working memory

Working memory is a basic human cognitive function. However, the genetic signatures and their biological pathway remain poorly understood. In the present study, we tried to clarify this issue by exploring the potential associations and pathways among genetic variants, brain morphometry and working memory performance. We first carried out association analyses between 2‐back accuracy and 212 image‐derived phenotypes from 1141 Human Connectome Project (HCP) subjects using a linear mixed model (LMM). We found a significantly positive correlation between the left cuneus volume and 2‐back accuracy (T = 3.615, p = 3.150e−4, Cohen''s d = 0.226, corrected using family‐wise error [FWE] method). Based on the LMM‐based genome‐wide association study (GWAS) on the HCP dataset and UK Biobank 33 k GWAS summary statistics, we identified eight independent single nucleotide polymorphisms (SNPs) that were reliably associated with left cuneus volume in both UKB and HCP dataset. Within the eight SNPs, we found a negative correlation between the rs76119478 polymorphism and 2‐back accuracy accuracy (T = −2.045, p = .041, Cohen''s d = −0.129). Finally, an LMM‐based mediation analysis elucidated a significant effect of left cuneus volume in mediating rs76119478 polymorphism on the 2‐back accuracy (indirect effect = −0.007, 95% BCa CI = [−0.045, −0.003]). These results were also replicated in a subgroup of Caucasians in the HCP population. Further fine mapping demonstrated that rs76119478 maps on intergene CTD‐2315A10.2 adjacent to protein‐encoding gene DAAM1, and is significantly associated with L3HYPDH mRNA expression. Our study suggested this new variant rs76119478 may regulate the working memory through exerting influence on the left cuneus volume.