Previous studies reported association between the ZNF804A gene and SZC, SAD and BPD, particularly with the common variant rs1344706. Our finding supports that in Mexican population the rs1344706 is also associated with SCZ, SAD and BPD. It has been suggested that this SNP can alter the expression of the ZNF804A gene at different stages of development, affecting neuronal growth, dendritic spine maintenance, and synaptic function; processes that have been related to the pathophysiology of these disorders (Deans et al. 2016; Hill et al. 2012; Hill & Bray, 2012; Tao et al. 2014). Also, different studies have described that the ZNF804A gene can modify the expression of other genes involved in neurogenesis, synaptogenesis, neurotransmission and inflammation (Chapman et al. 2018; Chen et al. 2015; Deans et al. 2016; Zhou et al. 2017). Therefore, it is possible that rs1344706 SNP, by modifying ZNF804A gene expression may affect the expression of other genes, modifying neuronal processes that would lead to dysfunctional neural circuits and structural alterations reported in patients with SCZ, SAD and BPD. These findings are supported by studies that found an association between rs1344706 and dysfunction in brain connectivity or alterations in gray and white matter (Chang et al. 2017; Wang et al. 2019). However, more studies are still needed to understand how ZNF804A gene is related to the pathophysiology of SZC, SAD and BPD. For example, little is known about the function of the different mRNA isoforms or proteic isoforms encoded by this gene, as well as the change in their expression because of the rs1344706 variants. An allelic specific binding of different transcription and splicing factors has been related to this SNP (Hess et al. 2015; Hess & Glatt, 2014). Nevertheless, several results come from studies in silico, so more experimental studies are required.
It is important to analyze if rs1344706 in LD with other SNPs localized in the same gene and/or the interaction between ZNF804A and other genes may be related with the etiology of these disorders. A bioinformatic study found an interaction between rs1344706 and rs11685645 of Myt1L gene, and between rs1344706 and rs1836796 of NCAM1 gene, concluding that these SNPs may be implicated in the early development of central nervous system (Hess et al. 2015; Riley et al. 2010). Also, interaction study between rs1344706 and a SNP of CACNA1C reported association with increase activation in brain regions of SCZ and BPD patients (Tecelão et al. 2019). In addition, it was reported an interaction effect of rs1344706 and COMT variant on brain activity in healthy individuals (Cui et al. 2019). However, more studies that analyze the interaction between rs1344706 and other SNPs are required. On the other hand, we do not know how these gene variants interact with the environment, in order to understand the biological relationship between SZC, SAD and BPD and the ZNF804A gene.
Genetic studies carried out in Caucasian and Asian populations reported association with the rs1388842; however, we did not find an association in Mexican patients with SCZ, BPD and SAD. Also, we did not find an association between rs3046266 and rs12477430 with these disorders. The present study replicates the previous findings of negative association between rs12477430 and SCZ in Irish population (Riley et al. 2010).
Furthermore, within a targeted 352 bp region located on exon four of ZNF804A gene we did not identify single rare variants in SCZ, BPD and SAD Mexican patients. Similarly, Dwyer S. et al. (2010) in a sequencing study of exons 2, 3 and 4 of ZNF804A gene did not identified rare variants in SZC patients of UK population.
Our results replicate that ZNF804A gene is involved in the genetic risk of SCZ, SAD and BPD. However, much research is still required, as obviously the mechanisms involved in the development of these highly heterogeneous disorders are more complex. The differences between these disorders and even between individuals with the same disorder are probably explained by the numerous combinations of the different genetic and environmental risk factors that interact with each other; coupled with the sum of different additional modifying factors, either genetic and/or environmental, that do not influence the risk of developing a disorder, but influence its clinical presentation; explaining the difficulty to understand the biological role of the ZNF804A gene in this disorders.
In addition, we consider that in the future it would also be important to study the role of ZNF804A gene in other psychiatric disorders, principally neurodevelopmental disorders, or those with psychotic symptoms, such as autism, primary delusional disorder, major depressive disorder with psychotic symptoms and substance-induced psychotic disorders.
The present study has several limitations. First, the small sample size in our analysis resulted in limited power to identify rare variants, principally ultra-rare variants. Second, we only sequenced a targeted region of exon four of the gene, consequently the existence of rare variants in other regions of the ZNF804A gene in SZC, SAD and BPD is possible; therefore, our study did not exclude the involvement of rare variants in the risk of presenting these disorders in Mexican population. It will be important sequencing the complete gene and replicate the association between rs1344706 and SZC, SAD and BPD in a large sample size.