Are There Schizophrenia Genetic Markers and Mutations ? A Systematic Review and Meta-Analyses

Background: Schizophrenia is a severe psychiatric disorder with a complex genetic factor determining its disease onset. Nevertheless, it is not clear in this mental disorder. Objective: To conduct a systematic review of articles regarding the genetic markers and mutations in schizophrenia. Methods: A systematic review of articles on genetic markers and mutations in schizophrenia, published from January 1, 2011, to September 7, 2015, on SCOPUS database was carried out. Search terms were “Genetic markers”, “Mutation”, and “Schizophrenia”. Results: Of the 527 retrieved studies, 31 met the eligibility criteria. Genetic polymorphism, Immune-associated genes, TCF4 and ZNF804A association with microRNA, Neuregulin gene, Chromosome 13q32 and 11p15.4, genes involved in glutamatergic via schizophrenia and brain structure, appeared to be associated with the origin of schizophrenia. Conclusion: Some studies show genes involved in several pathways leading to the disease pathogenesis such as that one related with the dopaminergic and immune system, or rare alleles. Some genes present no involvement in the etiology of this mental disorder. These findings clarify the genetic complexity of schizophrenia and affirm that together, the genes have an overall effect greater than the sum of the individual effect of each gene.


Introduction
Schizophrenia (MIM 181500) is a complex disease that has a lifetime risk of approximately 1% and is characterized by delusions, hallucinations, altered cognition, emotional reactivity and disorganized behavior.Genetic factors account for more than 80% of the variance in susceptibility, and risk likely results from multiple loci of small effect [1] [2].Supporting evidence has linked the high heritability of schizophrenia (SCZ) to a combination of relatively common alleles of small effect and to a few rare alleles with relatively large effects [3] [4].Genes like GRIK2, GRIA2, [5] ADAMTSL3 [6] and polymorphism present in the MAOA promoter region was associated with schizophrenia [7].Due to its complex genetic architecture and joint effects among these genes, the overall effect of a gene network is expected to have a greater effect than the sum of individual effect of each gene [8].Some of the most investigated genes in studies of susceptibility to schizophrenia are those that encode proteins of the dopaminergic system [4].The dysfunction of the dopamine system is believed to be a fundamental component in SCZ development [9].Abnormal transmission of dopamine may be involved in the pathogenesis of SCZ [10] [11] [12].Beyond that, there is compelling evidence that abnormal brain development and disturbed neuroplasticity are major factors in schizophrenia (SCZ) pathology [1] [2] [3] [13].
Also, was found that the genes involved are mainly neural-and immune-related and more likely to interact and take part in the same or related pathways [8].The abnormal functioning of the immune system including the irregulation and malfunction of different parameters of immune system are strongly associated with pathogenesis of schizophrenia.There is a strong relationship between the clinical characteristics of schizophrenia and disturbances of the immune system [14] [15].The major histocompatibility complex (MHC), an immune response gene locus chromosome 6, is the most extensively associated locus for schizophrenia in genome-wide association study (GWAS) [16] [17] [18].However, clinical and genetic heterogeneity and overlapping with other neurodevelopmental disorders complicate our understanding of the etiology of schizophrenia [18].
This research was conducted through a systematic review aiming to identify genetic markers and mutations for schizophrenia.Therefore, this work produces a synthesis of the latest findings about the genetic influence in this mental disorder.

Methods
A meta-analyse and systematic review was performed of articles about genetic markers, mutations, and their relations with schizophrenia.
It was conducted a search in the literature through the online database Scopus, in September 2015, by limiting itself to articles published between January 1, 2011 to September 7, 2015.The reason to limit the search between 2011 and results of research involving genetic mutations and markers for schizophrenia.
Therefore, genetic factors involving schizophrenia had great relevance in the scientific community.
Initially, the search terms browsed in Scopus database were: #1: "genetic markers" (MeSH term); #2: "mutation" (MeSH term); and #3: "schizophrenia" (MeSH term).Analysis of the article followed predetermined eligibility criteria.The survey was carried out in two phases: 1 AND 2 AND 3, 1 AND 3. It was made for those combinations by using the filter "subject descriptor".The index utilized was: Title, Abstract and Keywords.We adopted the following inclusion criteria: 1) written publications in English; 2) studies pertaining schizophrenia and mutations or genetic markers; 3) original articles with available full text online; 4) articles that included in the title at least one combination of terms described in the search strategy; and 5) observational (analytical or descriptive, except case reports), experimental or quasi-experimental studies (except animal studies), both prospective and retrospective.
The exclusion criteria were as follows: 1) other study designs, e.g.case reports, case series, literature reviews and comments; 2) non-original studies, including editorials, reviews, forewords, short communications and letters to the editor.
Each article was read in its entirety, and the information was entered in a spreadsheet that included authors, year of publication, the study sample description and key data.Some works found about schizophrenia differed from the proposed theme and were not included because they treat schizophrenia associated with other psychological disorders such as bipolar disorder, by addressing only symptomatic changes or because they are related to disease treatment.BioEstat 5.0 program was used to calculate the Mantel-Haenszel test and the Odds Ratio (Figure 1) in order to verify the association between the presence of the gene and the chance of developing schizophrenia; the Chi-Square Distribution (Figure 2), which is a probability distribution of the event actually occurs; and the Pearson's correlation coefficient (Figure 3), which represents the degree of association between certain variables.
The following Figure 1 shows the chances ratio of schizophrenia development associated to the genetic factor.
The following Figure 2 represents the probability distribution according to the associated event.
The following Figure 3 represents the Pearson's correlation coefficient, which highlights the genetic markers association with the schizophrenia development according its correlation strength.
To better analyze the data, the following stage involved the comparison between the articles and the division of the results obtained from the Reading of each one of them in eight categories: GENETIC POLYMORPHISM (ANNEXIN A5, THBS1, PDLIM5, HDAC, DISC 1 gene, KCNJ3, RGS2, MAOA gene, Figure 1.Mantel-Haenszel test and odds ratio.Prepared by the authors.The Mantel-Haenszel test and the odds ratio were used in order to calculate the genetic association and the disease establishment.Each set of studies is represented by a line and a square.The former represents the confidence interval and the latter the studies' effect.The square size represents the importance of each set of studies to the meta-analyses.The vertical line shows the effect absence while the diamond symbolizes the final meta-analysis result.The joint analysis involving all study patients had odds ratio 9.07 CI [8.540 -9.64].This means that if someone has a genetic load compatible with schizophrenia, the chance of developing this disease is 9 times higher. Figure 2. Chi-square distribution, prepared by the authors.The chi-square distribution was used in order to calculate the probability distribution of the event actually occurs.This involves the freedom degree and it is based on the odds ratio and on the chi-square calculated together and also according to the odds ratio set.With the p = 0.0430, the hypothesis of genetic association on the schizophrenia onset is not rejected.ANKK1, DRD3, GRIK1); IMMUNE-ASSOCIATED GENES (NKAPL, PGBD1, RELA), TCF4 and ZNF804A; ASSOCIATION WITH MICRORNA; NEUREGULIN GENE (NRG1, ERBB4); CHROMOSSOME 13q32 and 11p15.4;GENES INVOLVED IN GLUTAMATERGIC VIA (PTPN5 gene, GRIN2B, GRIK1, GRIK2, GRIA2, DTNBP1 gene, FXYD6); SCHIZOPHRENIA AND BRAIN STRUCTURE (B3GAT2, ADAMTSL3, DTNBP1 gene); and finally, NEGATIVE RESULTS (PEA15, ENTPD4, GAS2L1, ZIC2, SLC15A1, and FGF14).

Results
At first, the aforementioned search strategies resulted in 527 references.After searching the title and abstract of the considered citations for eligibility based on study inclusion criteria, 496 articles were excluded and 31 articles were further retrieved and included in the final sample (Figure 4).three studies).Among the 31 studies, some works were referenced in more than one category.The categorization of studies aims to a better organizational quality systematic review and it is not compulsory that each article must be referenced only in their respective category.Table 1 provides an overview of all studies included in the final sample and of all data elements used during the data analysis process.

Discussion
Schizophrenia is a heterogeneous and complex disorder, the etiology of which remains unknown [19].However, it is considered a complex trait resulting from both genetic and shared environmental etiological influences as reviewed by Sullivan [20].The genetic contribution to risk is high and heritability estimates based on clinical ascertainment are usually given as over 80% [21] [22] [23] (Figure 5).

Annexin A5 and Apoptotic Processes
It is proposed that both pre and postnatal as well as genetically determined abnormalities of the apoptotic processes might be among factors responsible for the development of schizophrenia [24].
It was shown that annexin A5 coupled to membrane of apoptotic cells is a ligand for C1q protein [25].
The increased levels of annexin A5 and Hficolin in the blood of schizophrenia patients might result from the intensification of the processes of programmed death of cells in the central blood circulation or/and might reflect alterations at the level of neuronal apoptosis due to increase of the bloodbrain barrier permeability in this disorder [26].
The results on genotyping the rs1157945 polymorphism of the annexin A5 gene definitely suggest the association between this mutation and schizophrenia, which enable one to consider the rs1157945*T minor allele of the annexin A5 gene to be a risk factor for disease development [27].Results confirmed previous findings from different ethnic populations, in that the rs1800497 and rs909253 polymorphisms were both associated with risk of schizophrenia.Differences between the genotypes of cases and controls were strongly significant (P = 0.0005 for rs1800497 and P = 0.001 for rs909253).The relative risk to schizophrenia was 1.2 (P = 0.01) for the C allele and 0.8 (P = 0.04) for the G allele.The CC, GG, and combined CC/AA genotypes were all more frequent in cases than in controls.These results support an association between ANKK1 and LTA genetic markers and vulnerability to schizophrenia and show the potential influence of just one copy of the mutant C or G allele in the Egyptian population.
This study reports the upregulation of JUN, HIST2H2BE, FOSB, FOS, EGR1 and TCF4 in the fibroblasts of SCZ patients.A significant alteration in EGR1 expression is also present in SCZ PBCs compared to controls and to MDD and BD patients, suggesting that this gene could be a specific biomarker helpful in the differential diagnosis of major psychoses.
It was found a secondary association between PDLIM5 variants and the paranoid subtype of schizophrenia in Emirati Arabs suggesting that PDLIM5 may represent a determinate/marker for schizophrenia subtype specification.However, no associations were found with variants in PICK1, NRG3 or DISC1 genes.
The study supports the notion that the DISC1 gene is a marker of schizophrenia susceptibility and that rs2509382 in the mutual DISC1 translocation region is a susceptibility marker for schizophrenia among males in Malaysia.

Dai et al. (2014) [11]
A total of 30 paranoid (15 males and 15 females) and 29 undifferentiated (15 males and 14 females) SCZ patients were collected as cases from Ningbo Kangning Hospital.In addition, 26 age-and gender-matched volunteers in Ningbo Kangning Hospital were enrolled as healthy controls (12 males and 15 females).
Their findings supported that DRD3 gene body hypermethylation was significantly associated with the risk of SCZ.

Gareeva et al. (2014) [5]
A test sample included 257 patients diagnosed with PSZ F.20.0xx according to the 10th revision of the International Classification of Diseases (ICD10), of which 108 ( The findings support the hypothesis that glutamate receptor genes are involved in the etiology and pathogenesis of schizophrenia.
Joshi et al.
(2014) [40] Postmortem frozen tissues from 30 schizophrenia cases, 7 schizoaffective cases, and 37 well-matched control individuals were provided by the New South Wales Tissue Resource Centre (TRC cohort).
Their findings demonstrate that ErbB4-pan laminar mRNA expression is elevated (layers I, II, V) in schizophrenia.At the cellular level, ErbB4-pan mRNA+ signal was detected predominantly in interneuron-like neurons.They provide evidence from this independent Australian postmortem cohort that ErbB4-JMa expression is elevated in schizophrenia and is linked to deficits in dendrite-targeting somatostatin, neuropeptide Y and vasoactive intestinal peptide interneurons.
Kebir et al.
(2014) [41] Families were included in this study based on the presence of one proband diagnosed with schizophrenia or schizoaffective disorder.The global sample included 951 Caucasian subjects in 313 nuclear families (325 cases and their 626 parents forming This first exploratory systematic study of the HDAC genes provides consistent support for the involvement of the HDAC3 gene in the etiology of schizophrenia.A statistical epistatic interaction between Continued 325 case-parent trios).It was composed of four samples: 1) The REFAPSY cohort: 274 subjects in 88 nuclear families (98 cases and 176 parents forming 98 case-parent trios) (multi-center enrollment, France), 2) The "Colombes" cohort: 182 subjects in 60 nuclear families (62 cases and 120 parents forming 62 case-parent trios) (Paris, France), 3) The "Barcelona" cohort: 411 subjects in 137 nuclear families (137 cases and 274 parents forming 137 case-parent trios) (multi-center enrollment, Spain), and 4) The "Montreal" cohort: 84 subjects in 28 nuclear families (28 cases and 56 founders.forming 28 case-parent trios) (Montreal, Canada).HDAC9, HDAC10, and HDAC11 was detected and seems biologically plausible.
As a result, it was identified one gene set with a joint effect significantly associated with schizophrenia and gene expression profiling analysis suggested that they were mainly neuro-and immune-related genes, such as glutamatergic gene (GRM5), GABAergic genes (GABRB1, GABARAP) and genes located in the MHC region (HLA-C, TAP2, HIST1H1B).Further pathway enrichment analysis suggested that these genes are involved in processes related to neuronal and immune systems, such as the Adherens junction pathway, the Neurotrophin signaling pathway and the Toll-like receptor signaling pathway.

Gareeva et al. (2013) [52]
A group of 258 patients diagnosed with paranoid schizophrenia F20.0xx (according to the international classification of diseases ICD10, paranoid schizophrenia in subjects of different ethnic backgrounds) including 110 Russians (50 women and 60 men) and 148 Tatars (68 women and 80 men) undergoing treatment in the Republican Psychiatric Hospital No. 1 of the Ministry of Health of the Republic Bashkortostan participated in the study.
The results are consistent with those obtained previously and support the hypothesis concerning the association of RGS2 gene polymorphisms with the risk of extrapyramidal syndrome development during haloperidol therapy and their involvement in the etiology and pathogenesis of schizophrenia.factor", 2p24.3 and 7q31.1 for "prodromal impairment factor", 1q32.1, 9p21.3, and 9q31.2 for "negative symptom factor", and 10p13 for "disorganization factor".Among these loci, chromosome 2p24.3and 1q32.1 overlap with susceptibility loci of schizophrenia identified in our previous linkage studies.This study suggests the existence of genetic loci related to various clinical features of schizophrenia.Further genetic analyses for these dimensional phenotypes are warranted.
Cao et al.
(2013) [45] The case samples were 310 random unrelated schizophrenia patients (144 male and 166 female) Their results provide further evidence for an effect of the DISC1 gene on the etiology of schizophrenia and suggest that STRs in the DISC1 gene may be genetic risk factors for schizophrenia.genotype (p = 0.038; OR = 2.12) and GRIN2B*T allele (p = 0.028; OR = 2.03), rs1805247, genotype GRIN2B*A/*A (p = 0.042; OR = 2.12), rs1805476; 2) genetic markers of the reduced risk of developing paranoid schizophrenia; 3) genetic markers of therapy response and the risk of side effects development during neuroleptics (haloperidol) treatment in Bashkortostan.

Boyajyan et al. (2013)
[27] Study subjects were patients with the paranoid form of schizophrenia diagnosed by psychiatrists at the Psychiatric Medical Center of the Ministry of Health of the Republic of Armenia based on the criteria of the International Classification of Diseases (ICD10 code: F20.0).In total, 225 chronic patients treated with typical neuroleptic haloperidol, 25 firstepisode neurolepticfree patients, and 225 healthy subjects (control group) were investigated.In the group of chronic patients 154 males and 71 females were enrolled.
It was shown that the pathogenesis of schizophrenia is characterized by an increased rate of apoptosis, which is more pronounced in the case of the first-episode neuroleptic-free patients than in the case of chronic patients that receive typical neuroleptic haloperidol.It was also shown that the rs11575945 polymorphism of the annexin A5 gene is associated with schizophrenia, and its minor allele is responsible for higher levels of the annexin A5 protein in the blood and represents one of the risk factors for the development of this disease.
13q32 has been reported to be linked to schizophrenia by multiple different studies.Thus, this study provides additional supporting evidence for an aetiological role of variants at 13q32 in schizophrenia.
Zhang et al.
(2012) [12] 768 patients and 1348 controls.A second sample used 963 patients and 992 controls in the Han Chinese population.
Association analysis of nine SNPs failed to detect any positive markers or haplotypes.Then, was tested rs4680 in a validation and no significant association was observed, but the cases significantly deviated from Hardy-Weinberg equilibrium (p = 5.7e-4).There was no association of rs4680 with SZ in the combined sample (n = 4071, p = 0.110, odds ratio = 1.08).
Hirata et al.
(2012) [104] 202 unrelated patients with a DSM-III-R or DSM-IV diagnosis of schizophrenia, as well as for the schizophrenia probands of 108 small nuclear families.Case-control sample and family based sample are completely independent.Each schizophrenia case was matched with a healthy control (N = 211) on the basis of age, ethnicity, and sex.There were 128 male and 74 female cases.In the case-control study, the Caucasian-only sample was 169 of 202 cases (83.7%).In the family-based study, the Caucasian subset was 73 of 108 probands (67.5%).
Was analyzed eight SNPs across the GRIK1 gene.The marker rs469472 was associated with schizophrenia when we combined the case-control and family samples (p = 0.027).
Zhang et al.
Was found association between ZNF804A and schizophrenia.ZNF804A single nucleotide polymorphisms (SNPs) encompassing exon 4 by performing an association study.In addition, was observed a significant association between marker rsl344706 and schizophrenia (P < 1.0 × 10 −5 ) in combined populations.
We genotyped two cSNPs [rs2228261 (Asn470Asn) and rs2292305 (Thr523Ala)] using direct sequencing.In this study, rs2228261 revealed significant association with schizophrenia in both codom-inant and recessive models.Also, rs2292305 was associated with schizophrenia in the recessive model.The results suggest that the THBS1 gene may contribute to the susceptibility of schizophrenia.(46.5% male (773/890), mean age ± SD: 46.9 ± 20.7 years).
Was found significant associations of three SNPs (rs11820062: p0.00011, rs2306365: p0.0031, and rs7119750: p0.0080) with schizophrenia and stronger evidence for association in a multi-marker sliding window haplotype analysis (the lowest p0.00006).It was genotyped six single-nucleotide polymorphisms (SNPs) in this region of FXYD6 in 1142 Han Chinese subjects (576 cases and 566 controls), and performed an association analysis.Significant associations with schizophrenia and the marker rs11544201 (P = 0.0028) and the haplotype rs10790212-rs11544201 (global P = 0.005) were found.
Naz et al.
(2011) [15] Schizophrenia patients and in controls.One hundred patients with schizophrenia and 70 healthy controls participated in the present study.
It was examined the association between the neuregulin 1 and tumor necrosis factor-α (−308) gene polymorphism with schizophrenia.This association was performed on the basis of molecular biology to screen the mutations of neuregulin 1 and tumor necrosis factor-α (−308) gene in schizophrenic patients by polymorphism analysis.Statistical analysis of the observed data shows that there was an association (P = 0.003) between patient's group and controls in terms of genotypes of single-nucleotide polymorphism 1 rather than single-nucleotide polymorphism 2 of neuregulin 1.

ChuanYuan et al. (2011)
[19] 81 Chinese Va schizophrenic patients and 210 Chinese Va healthy subjects The authors found that genetic variants of the FZD3 gene may affect susceptibility to schizophrenia in Chinese Han and Va populations.
(2011) [98] 868 Israeli Jewish participants, of whom 286 were patients with SZ and 582 were healthy controls.SZ patients were included with information regarding ethnicity (Ashkenazi vs. non-Ashkenazi origin).
The results imply a model in which PTPN5 may play a role in normal cognitive functioning and contribute to aspects of the neuropathology of SZ.

Dow et al.
(2011) [6] Whole blood samples were randomly selected from 92 schizophrenia patients from both the Munich and Aberdeen subjects.
Imputation analysis refined association between ADAMTSL3 and schizophrenia, and highlighted additional common variants with similar levels of association.We evaluated the functional consequences of all variants identified by sequencing, or showing direct or imputed association.The strongest evidence for function remained with the originally associated variant, rs950169, suggesting that this variant may be causal of the association.Rare variants were also identified with possible functional impact.
Surface-based analysis of the cortical mantle showed that the dysbindin haplotype was associated with structural differences in the medial orbitofrontal cortex, where the risk carriers showed the highest cortical thickness values and the non-risk carriers the lowest.This study extends previous evidence found on schizophrenic patients to the healthy population, demonstrating the influence of dysbindin risk variants on the neuronal architecture of a specific brain region relevant to the neuropathology of schizophrenia.
The study fails to provide evidence for the contribution of PEA15, ENTPD4, and GAS2L1 genes to the etiology of schizophrenia in the Japanese population.
Sun et al.
Haplotypic association was found for the VNTR-rs6323, VNTR rs1137070, and VNTRrs6323-rs1137070 haplotypes in female subjects.Nevertheless, no significant change of the expression of MAOA mRNA was detected in either female or male patients with paranoid schizophrenia.
The present results suggest that effects on biosynthesis of the neuronal epitope HNK-1, through common B3GAT2 variation, could increase the risk of SCZ, possibly by decreasing cortical area.
Figure 5. Genes that influence the susceptibility of schizophrenia.

Synaptic Regulation and Schizophrenia
Recently, thrombospondin 1 (THBS1) was discovered to be the essential astrocyte-derived synaptogenesis-promoting factor, and has been shown to contribute to the repair of brain injury through its assistance and promotion of neurogenesis [28] [29] [30].
Another work reported a gene PDLIM5 that is localized to the postsynaptic density where it has an important role in limiting the size of dendritic spines-the small synaptic protrusions that serve as the primary sites of excitatory synaptic transmission in the CNS [31] [32].
THBS1 is a 16.39 kb gene located on chromosome 15q15.The chromosome 15q15 region has been reported to be a susceptibility locus on schizophrenia [33] [34] [35], and clinical symptoms of schizophrenia such as periodic catatonia [33] [34].
PDLIM5 was found to be a significant secondary predictor of the paranoid subtype of schizophrenia in this Emirati Arab cohort [18].
Little is known about the biological function of THBS1 in the brain or in schizophrenia.However, given the involvement of astrocytes in the pathogenesis of schizophrenia, and the role of THBS1 in synaptic alteration, it is speculated that THBS1 may be a candidate gene involved in schizophrenia [33] [34] [35].
The finding here of PDLIM5 as a potential marker of schizophrenia subtype expands the role of synaptogenesis in neuropsychiatric disorders more generally [36] [37].

HDAC (Histone Deacetylases)
Several studies suggest the involvement of HDAC dysfunction in major psychotic disorders.First, mRNA expression levels of GAD67 in the prefrontal cortex of patients with schizophrenia were found to be strongly and negatively cor-related with mRNA expression levels of HDACs 1, 3, and 4 [38].In the same way, a common variant in HDAC3 region and epistatic interactions between HDAC9, HDAC10 and HDAC11 genes were found associated with schizophrenia.Second, HDAC activity has been found to be enhanced in the prefrontal cortex of patients with schizophrenia [39].
Given that HDAC proteins interact and form multiprotein complexes [40], and that a HDAC protein seldom operates alone, was tested that epistatic interaction between polymorphisms of HDAC genes may confer an increased risk of schizophrenia [41].

DISC1 Gene
Mutant DISC1 is proposed to contribute to schizophrenia susceptibility by disrupting intracellular transport, neurite modeling, neuronal migration, and proper development of the cerebral cortex [42] [43] [44].A study found three highly polymorphic STR loci in introns 1, 8, and 9 of the DISC1 gene and identified two novel loci, including (ATCC)n1 and (ATCC)n2.The three STRs showed a significant association with schizophrenia.It is noteworthy that the allele distribution of D1S1621 showed a higher risk of schizophrenia [45].
Another study found two SNPs (rs4658971 and rs1538979) to be significantly associated with schizophrenia in terms of both genotypes and allelic distribution.
Additionally, it was found an association of rs2509382 with schizophrenia among male participants [46].

Gene Association Involved with G-Protein and Schizophrenia
In the research of Yamada et al. [47], a gene found is the KCNJ3, also known as GIRK1/KIR3.1.This gene is a novel candidate gene for schizophrenia.The gene encodes a G protein-activated inwardly rectifying potassium channel (GIRK1/KIR3.1)and belongs KIR3.X the subfamily of potassium rectification inside channels.In that study, a single nucleotide polymorphism (SNP) (rs3106653) in the KCNJ3 (potassium inwardly rectifying channel, subfamily J, member 3) gene located at 2q24.1 showed association with schizophrenia in two independent sample sets.Moreover, experiments real-time quantitative RT-PCR showed that KCNJ3 is regulated in the prefrontal cortex schizophrenic patients.Decreased gene expression was observed in the brains from bipolar disorder, as well as schizophrenia patients.This finding may lend support to the ''hypo-NMDA theory of schizophrenia'' [47].
Associations between several polymorphisms of the RGS2 gene and the severity of schizophrenia were identified in patients by data show that the RGS2*G*G*T*C*T haplotype for the rs2746071, rs2746072, rs2746073, rs4606, and rs3767488 polymorphisms of the RGS2 gene is a genetic marker of increased risk of schizophrenia in Russian and Tatar residents of the Republic of Bashkortostan.Furthermore, the RGS2*G/*G genotype and the RGS2*G allele of the polymorphic rs2746071 locus of the RGS2 gene are risk factors for the development of schizophrenia in Russians and Tatars.The RGS2*A/*A genotype, the RGS2*A allele of the polymorphic rs2746071 locus of the RGS2 gene, and the RGS2*A*G*T*C*T haplotype (rs2746071, rs2746072, rs2746073, rs4606, and rs3767488) of the RGS2 gene are protective with regard to the risk of developing paranoid schizophrenia in Russians and Tatars [52].
The conclusion from their study is consistent with results obtained previously, and supports the hypothesis concerning the association of polymorphisms RGS2 gene involvement in the etiology and pathogenesis of schizophrenia [52].

Dopamine Association with Schizophrenia
Dopamine neurons in the substantia nigra are key neurotransmitters in human brain that contribute to the pathogenesis of SCZ [53].
The ''dopamine-serotonin'' imbalance hypothesis of schizophrenia claims that disruption in cortical-subcortical dopamine and serotonin neurotransmission may play an important role in the pathogenesis of schizophrenia [7] [54].
Monoamine oxidase A (MAOA) is the enzyme responsible for degradation of several monoamines, such as dopamine and serotonin that are considered as being two of the most important neurotransmitters involved in the pathophysiology of schizophrenia [7].In a study of Chinese population, all the informative SNPs and the VNTR polymorphism within the MAOA gene were screened.The present results suggest that haplotypes VNTR(L)-rs6323(T), VNTR(L)-rs1137070(C), and VNTR-(L)-rs6323(T)-rs1137070(C) may be associated with the increased risk of paranoid schizophrenia in female subjects [7].Other genes also investigated in susceptibility studies schizophrenia are those encoding proteins the dopaminergic system [4].
One such gene is the rs1800497 single nucleotide polymorphism (SNP) It was identified in exon 8 of repetition and Ankyrin kinase domain containing 1 gene (ANKK1) (MIM 608774) only 10 kb away from D2 (DRD2) dopamine receptor gene (located on 11q23.2,MIM 126450) in the 3'untranslated region.This polymorphism which leads to a substitution glutamic acid to lysine base (E713K) that can high-specificity substrate binding [55] probably modulates the function and expression of DRD2 due to its proximity [56].Despite the fact that the SNP rs1800497 is located ANKK1 the gene appears to be in linkage disequilibrium with several genetic variants DRD2, which could potentially explain a dopaminergic role in the pathogenesis of schizophrenia [4] [57].
Dai et al. [11] demonstrated that the DRD3 gene methylation, subtype of the dopamine receptor family, may have the potential to serve as a gender-specific biomarker to monitor the risk and development of schizophrenia.

Immune-Associated Genes
In recent years, the extended MHC region has been implicated as a main factor in schizophrenia pathogenesis, supported by GWASs of schizophrenia in differ-ent populations [16] [58] [59] [60] [61] [62].Products of the MHC region implicated in the pathogenesis of schizophrenia not only contribute to immune responses but also have general functions in different molecular biological processes [62] [63].Numerous other genes in the extended MHC region have been shown to be involved in the pathogenesis of schizophrenia.MICB, HLA-A, and HLA-B are classic MHC molecules that play central roles in the development of host defense and immunity [62] [64] [65].
In fact, more and more evidence suggested immune-related genes may play important roles in schizophrenia [8].For example, several GWAS revealed that many immune genes are significantly associated with schizophrenia [8] [59] [66] [67].In a replication study, all six SNPs within the NKAPL and PGBD1 genes displayed an association with schizophrenia, but neither of the two SNPs in the ZKSCAN4 gene showed an association [62].
Furthermore, study evidence was found that genetic variants of the RELA gene are associated with the risk for schizophrenia.The v-rel avian reticuloendotheliosis viral oncogene homolog A (RELA) gene encodes the major subunit of the NF-kB protein complex, are abundantly expressed in neurons and glia [68].
The NFKB3 located on chromosome 11q13 showed a suggestive linkage to schizophrenia in a family-based linkage disequilibrium analysis in a Japanese population [69].
A GWAS study in Caucasian population showed significant association with schizophrenia were in a region of LD on chromosome 6p22.1,including several immunity related genes other than the RELA gene [69].However, the biological significance of this gene in susceptibility for schizophrenia might not be large, because 22% of the patients with schizophrenia have homozygous of risk allele in SNP4, but 16% of the controls also are homozygous of risk allele in SNP4.At the same time, the association between the RELA gene and schizophrenia might explain, at least in part, the relation between immune system and schizophrenia [70].
The results not only provide further evidence that schizophrenia is a complex disease involving immune systems [16] [58] [59], but also indicate that the resultant genes identified in the modules are expressed in brain tissues, suggesting these genes may play important roles in brain function [8].

TCF4 and ZNF804A Association with microRNA
The most recent GWAS of schizophrenia identified five new association loci and the strongest new finding was with rs1625579 (P = 1.6 × 10 − 11) within an intron of a putative primary transcript for miRNA137, a known regulator of neuronal development [66] [71].In addition, the TCF4 gene is target of several mi-croRNAs, including the mir-137, which was implicated in SCZ etiology in GWAS studies [72].In particular, the SCZ-associated rs1625579 SNP in miR-137 was correlated with the decreased expression of this microRNA in the dorsolateral prefrontal cortex and a consequent increase in TCF4 levels [73] [74].In another analysis was suggested that ZNF804A may link miRNA137 to the specific molecular pathways or cellular processes implicated in schizophrenia, and experimental evidences were warranted to confirm the direct action between miRNA137 and ZNF804A.It has been suggested that the dysregulation of some miRNAs could be involved in the pathophysiology of schizophrenia [71] [75] [76].
Alternatively, elevated levels of cleaved ECTO-ERBB4 may result in increased NRG1 backward-signaling, and in turn, altered gene expression in schizophrenia.This observation is consistent with their previous findings linking the schizophrenia-associated HAPICE risk haplotype [80] with increased NRG1-Type III mRNA expression [83], and their earlier observations of increased cytoplasmic NRG1-ICD in schizophrenia [77] [84].
Moreover, Chen et al. [85] found strong evidence of association between "delusion factor" and SNPs of Neuregulin3 (NRG3) on chromosome 10q22-q23, a region identified in their previous genome-wide linkage analysis of schizophrenia [86].
Brzustowicz et al. [89] analyzed 21 Canadian families with schizophrenia and found genome-wide significant LOD scores on 13q.Brzustowicz et al. [90], using the same dataset and multipoint analysis and found a maximum LOD score of 3.81 with an empirical p = 0.02 under a recessive-broad model of schizophrenia at D13S793, with an estimated 65% of families linked to this region.Some more recent studies found evidence for a role of chromosome 11p15.4 in bipolar disorder [91], autism [92] and attention deficit hyperactivity disorder [93].This relative lack of previous evidence for chromosome 11p15 contrasted with the strong previous findings available for 13q32 [23].
Thus, although they should not disregard the linkage with 11p15.4,their pri-mary finding is an additional evidence for a potential role for 13q32 region in a small Brazilian sample of families with early onset schizophrenia affected individuals.The allele-sharing model tested is compatible with allelic homogeneity within the families [23].

Genes Involved in Glutamatergic Pathway
Hypoglutamatergic hypothesis, which postulates disturbances in the glutamate neurotransmitter pathway, is one of the neurochemical hypotheses of schizophrenia pathogenesis, which in recent times, has attracted attention from researchers [94].
There is considerable evidence suggesting the involvement of ionotrophic glutamate receptors in SCZ [95] [96].
Changes in affinity of glutamate receptors, transcription of their genes, and expression of their subunits in the prefrontal cortex, hippocampus, and thalamus have been revealed in schizophrenics in post mortem studies [97].
In a study, was examined the association of the protein tyrosine phosphatase non-receptor 5 (PTPN5) gene, which encodes for Striatal-enriched protein tyrosine phosphatase (STEP), with both schizophrenia and cognitive.The variants of protein tyrosine phosphatase non-receptor 5 (PTPN5), the STEP encoding gene involved in glutamate receptor trafficking, are associated with both the diagnosis of SCZ and with neurocognitive function in general [98].
The GRIN2B gene encoding the subunit NR2B NMDA glutamate receptor and its association with schizophrenia was studied in Russians and Tartars, yielding positive results [94].
In genetic studies, previously, they have found some evidence which is components of the glutamate system, namely the GRIN2B gene [99].
In fact, in primary cortical neuronal culture, dysbindin appears to influence exocytotic glutamate release via upregulation of molecules in presynaptic machinery [100] [102].
The GRIA2 ionotropic glutamate receptor gene consists of 16 exons and 15 introns and codes for the GluR2 subunit of the AMPA glutamate receptor [5].
The GRIK2 gene codes for subunit 2 of the kainite glutamate receptor and includes 17 exons.The gene is expressed in the granular layer of the cerebellum, the dentate gyrus of the hippocampus, and the neocortex [5].
The results point to potential involvement of GRIK1 in Caucasian schizophrenia patients [104].
From a functional perspective, FXYD6 is a member of the FXYD protein fam-ily.All members of this family have been shown to modulate Na, K-ATPase and have long-term physiological importance in maintaining cation homeostasis [107].
Interestingly, Na, K-ATPases and glutamate (the major excitatory neurotransmitter in the mammalian brain) are part of the same macromolecular complex and operate as a functional unit to regulate glutamatergic neurotransmission [2] [108].
Zhong et al. [2] reported that the G allele of rs11544201, located in the third exon of FXYD6, is a compelling risk factor for schizophrenia in Han Chinese.

Schizophrenia and Brain Structure
Cognitive deficits are a core feature in SCZ [109], and brain structural changes such as reduced hippocampal and cortical volume [110], the extracellular matrix (ECM) abnormalities [111], as well as reduced cortical thickness [112] have been reported in SCZ patients.
The Human Natural Killer-1 (HNK-1) carbohydrate is among the most characteristic glycoepitopes in the nervous system, involved in neural crest cell migration, neurite outgrowth, neuronal cell adhesion, and synaptic plasticity [113] [114].Kähler et al. [13] suggest that effects on biosynthesis of the neuronal epitope HNK-1, through common variation in B3GAT2, might decrease cortical surface area and thus increase the risk of SCZ.
ADAMTSL3 is a member of the ADAMTS superfamily of proteins, encompassing the 19 human ADAMTS metalloproteases and the seven non-proteolytic ADAMTS-like proteins [6].
ADAMTSL3 itself has been shown to be an important ECM component [115].
Considering the known biology of other superfamily members in the neural ECM, we propose that ADAMTSL3 could plausibly be involved in the PNN abnormalities seen in schizophrenia [111].
Meanwhile, the dysbindin (dystrobrevin-binding protein 1) gene has been indicated as one of the most important schizophrenia susceptibility genes.The dystrobrevin-binding protein 1 (DTNBP1) gene, known as the dysbindin gene, is located in chromosome 6p22.3,encodes 40 -50 kDa protein, which is widely expressed in neurons of the human brain, and plays a significant role in modulating glutamatergic neurotransmission [100] [101].Recently, the C-A-T haplotype (derived from SNPs rs2619539, rs3213207 and rs2619538) reported by Williams et al. [116] has shown a strong association with schizophrenia, reduced dysbindin mRNA expression in human post-mortem analysis of homogenized human brain [117], poorer cognitive performance in working memory [118] and reduced gray matter volume in the prefrontal cortex in patients with schizophrenia

Negative Results
In the study of genes that give the scenario susceptibility to schizophrenia, some genes were tested, and there were no satisfactory results, like PEA15 genes, ENTPD4, Genes GAS2L1 [120] and COMT [12].The PEA15 gene is located at 1q21.1, the gene is ENTPD4 on 8p21.3 and GAS2L1 gene is in 22q12.These three genes are homologous [120].The COMT gene is located at 22q11.2, which is a region that has been implicated in linkage areas of SCZ [12] [121].
This single-nucleotide polymorphism (SNP) has been widely studied for its connection with SCZ susceptibility, but with conflicting results [120].Saito et al. [120] suggests that PEA15 and GAS2L1 do not play an important role in susceptibility to schizophrenia in the Japanese population.Saito et al. [120] left ENTPD4 out of further analysis beyond the screening scan sample subset, since high P values for allelic frequencies suggest that this gene is also unlikely to be related to schizophrenia.Other genes have also been investigated in schizophrenia with predominantly negative results are ZIC2, SLC15A1, and FGF14 [23] (Figure 6).
However, further research on these genes which are located on schizophrenia susceptibility loci, using different SNP's and a large sample set will be required [120].

Strengths
This study involves a search for genetic risk markers for Schizophrenia; Besides, a meta-analysis of studies was carried out, allowing a better understanding and a more precise analysis of the relation of genetic changes and Schizophrenia onset.Most of the analyzed studies have shown an association, in a greater or lesser level, with Schizophrenia.Moreover, there is a need for more meta-analyses of studies on this topic in the literature.

Limitations
Trough this analysis it was not possible to determine what kinds of genetic Figure 6.Genes that were tested negative for susceptibility to schizophrenia.changes presented a higher risk to the schizophrenia development.Besides, the sample size was restricted due to the inclusion criteria.

Conclusion
This study shows numerous genetic targets have been identified related to susceptibility to schizophrenia.Great steps have been taken on this theme, but there is still a real need for more detailed studies like ours, addressing this same topic with a larger sample, in view of that the existence of genetic markers and mutations for SCZ is still controversial.Thus, gathering these markers and mutations already set to develop studies with this approach will enable the prevention, diagnosis and more targeted treatment for patients with schizophrenia, improving their quality of life.

Figure 3 .
Figure 3. Person's correlation coefficient, prepared by the authors.The Pearson correlation was used in order to verify the statistical relationship between the genetic variable and the schizophrenia onset.It is a parameter ranging from −1 to 1, in which −1 means a total negative correlation, 0 means no correlation and means 1 overall positive correlation.According to the analysis from the studies, we found a correlation outcome r = 0.79, which shows quantitatively that the schizophrenia development is associated with a genetic marker.

Figure 4 .
Figure 4. Flow chart showing study selection for the review.Abbreviations: MeSH: Medical Subject Headings.
1091 healthy controls in an attempt to replicate the GWAS results in the Chinese Han population Locus 6p21-p22.1 is significantly associated with schizophrenia in the Chinese Han population.To further investigate the association between polymorphisms at this locus and schizophrenia.Was selected eight other single-nucleotide polymorphisms (SNPs) distributed in or near these genes for a case-control association study in an independent.The examined subjects are represented by 251 patients with paranoid schizophre nia F.20.0xx in accordance with the International Classification of Diseases ICD10 of PS with In the course of the analysis, this study revealed the following: 1) genetic markers of increased risk of developing paranoid schizophrenia in various ethnic groups, including, in Tatars, the GRIN2B*T/*T genotype (p = 0.003; OR = 2.33) and GRIN2B*T allele (p = 0.001; OR = 2.36), rs1805247; in Russians, the GRIN2B*T/*T Continued Gareeva et al. (2013) [94] different ethnicities, including 108 Russians (48 women and 60 men) and 143 Tatars (67 women and 76 men), which were treated in RPH no. 1 at Ufa city.

Table 1 .
Are there schizophrenia genetic markers and mutation?Studies and main findings.
49 females and 59 males) were Russians and 149 (71 females and 78 males) were Tatars.The patients were treated at Republican Mental Hospital no. 1 (Ufa).The mean age was 24.94 ± 8.91 years; the mean age at disease onset, 22.46 ± 7.32 years; and the mean disease history, 3.69 ± 3.85 years.