Schizophrenia is a chronic and debilitating psychiatric condition affecting slightly more than 1% of the population worldwide and it is a multifactorial disorder with a high degree of heritability (80%) based on family and twin studies. Increasing lines of evidence suggest intermediate phenotypes/endophenotypes are more associated with causes of the disease and are less genetically complex than the broader disease spectrum. Negative symptoms in schizophrenia are attractive intermediate phenotypes based on their clinical and treatment response features. Therefore, our objective was to identify genetic variants underlying the negative symptoms of schizophrenia by analyzing two genome-wide association (GWA) data sets consisting of a total of 1,774 European-American patients and 2,726 controls. Logistic regression analysis of negative symptoms as a binary trait (adjusted for age and sex) was performed using PLINK. For meta-analysis of two datasets, the fixed-effect model in PLINK was applied. Through meta-analysis we identified 25 single nucleotide polymorphisms (SNPs) associated with negative symptoms with p<5×10(-5). Especially we detected five SNPs in the first two genes/loci strongly associated with negative symptoms of schizophrenia (P(meta-analysis)<6.22×10(-6)), which included three SNPs in the BCL9 gene: rs583583 showed the strongest association at a P(meta-analysis) of 6.00×10(-7) and two SNPs in the C9orf5 (the top SNP is rs643410 with a p = 1.29 ×10(-6)). Through meta-analysis, we identified several additional negative symptoms associated genes (ST3GAL1, RNF144, CTNNA3 and ZNF385D). This is the first report of the common variants influencing negative symptoms of schizophrenia. These results provide direct evidence of using of negative symptoms as an intermediate phenotype to dissect the complex genetics of schizophrenia. However, additional studies are warranted to examine the underlying mechanisms of these disease-associated SNPs in these genes.