Background: Personal variations in genetic risk for schizophrenia relate to its phenotypic heterogeneity—both in disorder development and clinical manifestations. Abnormal glutamatergic neurotransmitter system functioning is integrated in the pathogenesis of schizophrenia. Methods: A sample of 805 Russian schizophrenia patients from the Siberian Federal region was investigated. We examined the association of 39 single nucleotide polymorphisms in eight genes (GRIN2A, GRIN2B, SLC1A2, SLC1A3, SLC17A7, GRM3, GRM7, and GRM8) involved in the glutamatergic system with the development of clinical heterogeneity of schizophrenia. The MassARRAY Analyzer 4 was used for genotyping. Results: GRIN2A rs11644461, rs8057394 and GRIN2B rs7313149 are associated with the continuous type of schizophrenia. The GRIN2A rs8057394*G allele is a relative risk factor (р = 0.019) for developing the continuous type of schizophrenia. We found a nominally significant association between negative symptoms of schizophrenia and SLC17A7 rs62126236. The SLC17A7 rs62126236*T allele has a protective effect (р = 0.039) against predominant negative symptoms in schizophrenia. The total Positive and Negative Syndrome Scale (PANSS) scores were significantly associated with GRIN2A rs9788936 after adjusting for multiple testing (р = 0.001). Conclusions: In this study the contribution of the glutamatergic gene polymorphisms to the clinical heterogeneity of schizophrenia has been demonstrated.