Novel association between TGFA, TGFB1, IRF1, PTGS2 and IKBKB single-nucleotide polymorphisms and occurrence, severity and treatment response of major depressive disorder

Background Activation of the immune system might affect the severity of depressive episodes as well as response to the antidepressant treatment. The purpose of this study was to investigate whether the occurrence of variant alleles of analyzed SNPs are involved in prevalence and progression of depression. Moreover, selected genes and SNPs have not been investigated in context of the disease severity and treatment. Therefore, six polymorphisms were selected: g.41354391A>G-TGFB1 (rs1800469), g.132484229C>A-IRF (rs2070729), g.186643058A>G-PTGS2 (rs5275), g.186640617C>T-PTGS2 (rs4648308), g.70677994G>A-TGFA (rs2166975) and g.42140549G>T–IKBKB (rs5029748). Methods A total of 360 (180 patients and 180 controls) DNA samples were genotyped using TaqMan probes. Results We observed that A/G of the rs2166975 TGFA, A/C of rs2070729 IRF1 and G/T of rs5029748 IKBKB were associated with an increased risk of depression development while the T/T of rs5029748 IKBKB, T/T of rs4648308 PTGS2 and G/G of rs2166975 TGFA reduced this risk. We also stratified the study group according to gender and found that genotype A/G and allele G of the rs2166975 TGFA, G/T of rs5029748 IKBKB as well as C allele of rs4648308 PTGS2, homozygote A/A and allele A of rs5275 PTGS2 were associated with increased risk of depression development in men while homozygote G/G of rs5275 PTGS2 decreased this risk. Moreover, C/T of rs4648308 PTGS2 and A/G of rs5275 PTGS2 was positively correlated with the risk of the disease occurrence in women. Furthermore, a gene–gene analysis revealed a link between studied polymorphisms and depression. In addition, A/A of rs1800469 TGFB1 was associated with earlier age of onset of the disease while G/G of this SNP increased severity of the depressive episode. Interestingly, A/C of rs2070729 IRF1 and T/T of rs5029748 IKBKB may modulate the effectiveness of selective serotonin reuptake inhibitors therapy. In conclusion, studied SNPs may modulate the risk of occurrence, age of onset, severity of the disease and response to the antidepressant treatment.


INTRODUCTION
Depression (Major depressive disorder, MDD) is one of the most frequently diagnosed mental diseases. According to World Health Organization, about 350 million people suffer from this disorder all over the world (WHO, 2018). Despite the importance of the problem, pathogenesis of depression is not fully understood. However, there is a growing body of evidence suggesting that immune system impairment and dysregulation is associated with the pathophysiology of MDD. In particular, the "cytokine hypothesis" is widely accepted as one of the mechanisms for the development of depression (Capuron & Miller, 2011). This theory postulates that MDD is a result of elevated expression of pro-inflammatory cytokines, which act as neuromodulators as well as main agents in mediation of the neuroendocrine, neurochemical and behavioral features of the disease (Schiepers, Wichers & Maes, 2005). Some evidence confirmed link between inflammation and depression. Primarily, MDD patients exhibit increased levels of cytokines and other pro-inflammatory markers (Capuron & Miller, 2011). Additionally, medical conditions connected with increased inflammatory response are associated with greater risk of MDD developing (Capuron & Miller, 2011).
One of the cytokine class strongly associated with depression are interferons (IFN), cluster of signaling proteins involved in immune response. More than twenty different IFN proteins have been identified so far and divided into classes. IFN proteins are able to activate immune cells, that is, natural killer cells (NK cells) and macrophages (Pinto & Andrade, 2016). For instance, IFN-a is implicated in modulation of mood, behavior and sleep-wake cycle, partially by its ability to activate the pro-inflammatory cytokine network including, interleukin 1 (IL-1), interleukin 6 (IL-6) and tumor necrosis factor alpha (TNF-a) (Zahiu & Mihai, 2014). IFN and IFN-inducible genes, involved in immunity and inflammation, are transcriptionally regulated by interferon regulatory factor 1 (IRF1) (Tamura et al., 2008). IRF1 was a first identified transcription factor in IFN system and as a member of interferon regulatory factor family, plays important role in controlling expression of aforementioned genes (Kröger et al., 2002). Besides this, IRF1 promotes inflammatory cytokine release and regulates expression of interleukin 12 (IL-12) and interleukin 15 (IL-15), which are involved in MDD (Tamura et al., 2008).
Besides cytokine theory, various inflammatory pathways are thought to be activated in course of depression, including activation of the NF-kB (nuclear factor-kB), what leads to increased levels of pro-inflammatory cytokines (Bierhaus et al., 2003;Pace et al., 2006). NF-kB is a ubiquitous transcriptional factor that regulates expression of genes involved in pleiotropic functions, including pro-inflammatory cytokines and co-stimulatory molecules (Takeda & Akira, 2007;Krakauer, 2008;Zhang, Lenardo & Baltimore, 2017). Inactive NF-kB molecules retain in the cytoplasm by interaction with IkB proteins, allowing to immediate activation in response to adequate impulse (Napetschnig & Wu, 2013). Canonical signaling of NF-kB is activated by IkB kinase (IKK complex), consisting of three subunits, each encoded by separate gene, that is, IKK-a (Inhibitor of nuclear factor kappa-B kinase subunit alpha) encoded by CHUK gene, IKK-B (inhibitor of nuclear factor kappa-B kinase subunit beta) by IKBKB gene and IKK-g (inhibitor of nuclear factor kappa-B kinase subunit gamma) by IKBKG. The activation of IKK is induced by phosphorylation of serine residues in catalytic subunits of kinase complex (Napetschnig & Wu, 2013;Karin & Ben-Neriah, 2000;Cardinez et al., 2018). Therefore, defective expression of NF-kB as the pro-inflammatory transcription factor, caused by alterations in IKBKB gene, may play a role in the development of depression (Napetschnig & Wu, 2013).
Transforming growth factors (TGF) constitute of two classes of polypeptide growth factors, namely TGFA (transforming grow factor a) and TGFB (transforming grow factor β). Important functions of these cytokines are embryonic development and regulation of specific reactions of immune system by their ability to induce T regulatory cells (Treg) (Kissin et al., 2002;Yamagiwa et al., 2001). TGFA is a ligand for epidermal growth factor receptor, which stimulates cell migration and proliferation. These gene and protein have been associated with many types of cancers and other diseases (Ten Dijke & Hill, 2004). Another piece of evidence confirmed that TGFB, an anti-inflammatory cytokine, plays role in brain inflammation as well as in peripheral immune response. Namely, TGFB is mainly involved in regulating inflammatory response by induction of differentiation of CD4 + T cells (Nam et al., 2008;Passos et al., 2010). Another essential function of the protein is cell to cell signaling, and thus controlling of cell growth and differentiation (Ten Dijke & Hill, 2004). In addition, TGFB is able to exert neuroprotective effects in many neurodegenerative disorders (Vivien & Ali, 2006). Information about its role in depression are contradictory. On the one hand, in animal model of depression, the cytokine level is increased and causes imbalance between Treg and Th17 cells (Hong et al., 2013). On the other hand, some studies reported that levels of TGFB in depressed patients are lower than in healthy control group (Musil et al., 2011;Sutcigil et al., 2007). Moreover, TGFB alone is sufficient to stimulate production of pro-inflammatory cytokines for example, IL-1 and TNF-a (Kunzmann et al., 2003). The protein is also able to induce expression of prostaglandin-endoperoxide synthase 2 (PTGS2; cyclooxygenase-2-COX-2) encoded by PTGS2 gene, which is involved in pathogenesis of MDD. PTGS2 besides contribution to processes related to inflammation, also participates in the production of free radicals, which is partly utilized by PTGS2 itself (Aktan, 2004;Hansson, Olsson & Nauseef, 2006). Moreover, COX-2 catalyzes conversion of arachidonic acid (AA) to prostaglandins (PGs), which further intensify inflammation and neurodegenerative processes in central nervous system (CNS) (Minghetti, 2004). In response to growth factors, cytokines and other inflammatory molecules, PTGS2 is immediately expressed and is responsible for the production of prostanoid in both acute and chronic inflammatory conditions (Breyer et al., 2001;Shi et al., 2010). Additionally, in animal model of depression increased expression of PTGS2 was observed in brain regions (Cassano et al., 2006).
The evidence suggests that MDD may be associated with impairment of immune system, caused by defective activity of aforementioned genes. Moreover, genetic factors may play an essential role in development of depression, since genome-wide association studies (GWAS) found several regions significantly associated with MDD (Shyn et al., 2011;Wray et al., 2018). Therefore, the present study examines the prospective relationship between the occurrence, age of onset, severity or antidepressant treatment efficacy of MDD and appearance of single nucleotide polymorphism (SNP) located in inflammatory-related genes, that is, g.132484229C>A of IRF1 (rs2070729, located on 5q31.1), g.186643058A>G of PTGS2 (rs5275, located on 1q31.1), g.186640617C>T of PTGS2 (rs4648308, located on 1q31.1), g.70677994G>A of TGFA (rs2166975, located on 2p13.3), g.41354391A>G of TGFB1 (rs1800469, located on 19q13.2) and g.42140549G>T of IKBKB (rs5029748, located on 8p11.21). Selected SNPs are located within immune genes participating in inflammatory-related signaling pathways. Therefore, they could affect gene expression and protein function and thus contribute to immune disruptions leading to increased risk of MDD.

Subjects
The study included a total of 360 participants randomly selected. A group of 180 patients with depression hospitalized at the Department of Adult Psychiatry of the Medical University of Lodz and 180 volunteers without health problems, selected randomly (Table 1). Participants who took part in the experiment were native, not-related Poles. Patients were included based on the criteria set out in ICD-10 (F32.0-7.32.2, F33.0-F33.8). Medical and psychiatric records were obtained in accordance with ICD-10 criteria, using the Standardized Composite International Diagnostic Interview (CIDI). The depression' severity was evaluated using the 21-item Hamilton Depression Rating Scale (HDRS-21). The exclusion criteria included: axis I and II disorders other than MDD, chronic somatic diseases, autoimmune disorders (psoriasis, rheumatoid arthritis, chronic obstructive pulmonary disease, cancer, chronic kidney disease, systemic lupus erythematosus, type 1 diabetes, hepatitis B and C virus and HIV infection), neuroinflammatory and neurodegenerative disorders (including multiple sclerosis, Alzheimer's disease, Parkinson's disease) and central nervous system damage. Furthermore, subjects with familial incidence of mental diseases, other than MDD did not participate in the experiment. Psychiatric examination was conducted by the same psychiatrist, before the subjects were included in the experiment and after 8 weeks of pharmacotherapy with selective serotonin reuptake inhibitor (SSRI). Control group included selected randomly, volunteers with negative history of mental disorders. Participation in the experiment was voluntary. Controls and patients who did not agree to participate in the study were excluded. The purpose of the study was clearly presented, participants were assured that their personal information would be kept confidential. All of the subjects agreed by giving their written consent to participate in the experiment according to the protocol approved by the Bioethics Committee of the Medical University of Lodz (No. RNN/70/14/KE).

SNP selection
Selection of the studied polymorphisms was performed using the public domain of the database for single nucleotide polymorphisms of the National Center for Biotechnology (NCBI dbSNP, www.ncbi.nlm.nih.gov/snp/) (Bethesda, Montgomery County, MD, USA). The criteria used for the SNPs' selection were that the minor allele frequency is greater than 0.05 in the European population, and that they are located in the coding or regulatory region of the genes and may have functional meaning for transcription and protein function. Detailed information about selected polymorphisms are presented in Table 2.

DNA isolation
Genomic DNA was isolated from venous blood in accordance with the manufacturer instructions. Blood samples were collected from control group and patients with MDD. Blood Mini Kit (A&A Biotechnology, Gdynia, Poland) was used to extract nucleic acid. The purity of and concentration of the DNA was measured spectrophotometrically by calculating the ratio between absorbance at 260 nm and 280 nm, using Picodrop TM (Picodrop Limited, Astranet Systems Ltd., Cambridge, UK). Samples were stored at −20 C until use. Responsive (reduction from baseline of ≥ 50% in the total score) 93% Remission (total HRDS1 score ≤7) 66% Table 2 Characteristic of studied polymorphisms.

Statistical analysis
The collected data were analyzed in Statistica 12 (Statsoft, Tulsa, OK, USA), SigmaPlot 14.0 (Systat Software Inc., San Jose, CA, USA), Resampling Stats Add-in for Excel v.4 (Arlington, TX, USA) and StudSize3.02 (CreoStat HB, Florunda, Sweden). The descriptive statistics are shown as medians with interquartile ranges. Normality of the studied group was verified with the Shapiro-Wilk test, homogeneity of variance was checked with Brown-Forsythe test. Accordingly, either the unpaired Student's t test or Mann-Whitney U test was used. To calculate the associations between studied polymorphisms and the occurrence of a disease an unconditional multiple logistic regression model was used.
The results are shown as odds ratio (OR) with 95% confidence interval (95% CI). The OR values were adjusted for the potential confounders, including age and sex. We also stratified results into male and female group and evaluated correlation between case/ control for each polymorphism. In addition, in order to strengthen that the revealed differences were not detected by a pure chance the significant outcomes were further validated with the use of two approaches: the bootstrap-boosted multiple logistic regression (resampling with replacement, 10,000 iterations) and the cross-validated logistic regression (corresponding to the d-jackknife technique), with the patient group acting as the modeled class. This was intended to overcome any possible bias related to relatively low sample sizes. The goodness of fit of logistic regression models showing a significant degree of discrimination between controls and patients was estimated with Hosmer-Lemeshow test. Efficiency of the treatment was calculated using the formula as described before (Czarny et al., 2019): TE-treatment efficiency; HAM-D 0 -score before therapy; HAM-D E -score after therapy.

RESULTS
Single nucleotide polymorphisms of genes encoding IRF1, IKBKB, TGFA, TGFB and PTGS2 as a risk of MDD The distribution of genotypes and alleles in both depressed and control groups was in agreement with Hardy-Weinberg equilibrium. Results are presented in Table 3. The results demonstrated that the A/G genotype of the g.70677994G>A (rs2166975) Single-nucleotide polymorphisms of genes encoding IRF1, IKBKB, TGFA, TGFB and PTGS2 and MDD occurrence in male and female population Since women show two-times higher risk of MDD occurrence compared to men, we decided to investigated the association between prevalence of the disease in stratified male/female population and all studied SNPs. Results are presented in Table 4. The results demonstrated that in the case of g.70677994G>A (rs2166975) polymorphism of the TGFA, the A/G genotype increased the risk of MDD in men, but not in women. Moreover, allele A of this SNP was associated with decreased chance of the disease, while allele G was strongly correlated with higher risk of MDD. Furthermore, in male population allele G and G/G homozygote of the g.186643058A>G (rs5275) of PTGS2 decreased risk of depression while, allele A and A/A homozygote of the same polymorphism was associated with increased risk of the occurrence of the disease. Additionally, it was found that A/G genotype of this SNP was correlated with higher risk of MDD in the female group. Another SNP of PTGS2 gene, g.186640617C>T (rs4648308) was associated with MDD risk in both studied groups. Precisely, C/T genotype was positively correlated with the risk of the occurrence of MDD in women. Similarly, allele C of the mentioned polymorphism  increased prevalence of the disease among men, while allele T decreased this risk. We also found that genotypes of g.42140549G>T (rs5029748) polymorphism of IKBKB gene were related with appearance of MDD in male population. Particularly, G/T genotype was connected with increased risk of depression, while T/T genotype of the same SNP decreased this risk.

Gene-gene interactions of IRF1, IKBKB, TGFA, TGFB and PTGS2 and the risk of MDD
In this research, we also studied whether the combined genotypes of investigated polymorphism are associated with appearance of MDD. Results are presented in Table 5.
In reference to effect of combined genotypes, it was found that G/G-T/T genotypes of Notes: * 'Adjusted OR' means OR adjusted for sex and age; for significant comparisons the superscript b means the bootstrap-boosted OR (resampling with replacement, 10,000 iterations); all OR values without bootstrap analysis were calculated using cross-validation algorithm. Statistical power (1 − β) (calculated at a = 0.05) for significant comparisons given in superscripts. p < 0.05 along with corresponding ORs are in bold. Notes: * 'Adjusted OR' means OR adjusted for sex and age; for significant comparisons the superscript b means the bootstrap-boosted OR (resampling with replacement, 10,000 iterations); all OR values without bootstrap analysis were calculated using cross-validation algorithm. Statistical power (1 − β) (calculated at a = 0.05) for significant comparisons given in superscripts. p < 0.05 along with corresponding ORs are in bold. g.70677994G>A (rs2166975)-TGFA and g.186640617C>T (rs4648308)-PTGS2 was associated with decreased risk of depression occurrence, while A/G-C/T genotypes increased this risk. The A/G-A/C genotypes of g.70677994G>A (rs2166975)-TGFA and g.132484229C>A (rs2070729)-IRF1 as well as A/G-A/G genotypes of g.70677994G>A (rs2166975) TGFA and g.41354391A>G (rs1800469)-TGFB also increased the risk of the disease. Furthermore, higher risk of MDD occurrence was associated with the G/T-A/G genotypes of g.42140549G>T (rs5029748)-IKBKB and g.186643058A>G (rs5275)-PTGS2, however the T/T-G/G genotypes reduced this risk. In the case of linked genotypes of g.70677994G>A (rs2166975)-TGFA and g.186643058A>G (rs5275)-PTGS2, we found that link between A/G-G/G of this genes was associated with higher risk of appearance of the MDD, while G/G-G/G as well as A/A-G/G genotypes decreased this chance. Similarly, A/G-G/T combined genotypes of g.70677994G>A (rs2166975)-TGFA and g.42140549G>T (rs5029748)-IKBKB increased risk of MDD but G/G-T/T genotypes of the same SNP were associated with lower risk of disease incidence. Moreover, carriers of A/C-A/G combined genotypes of g.132484229C>A (rs2070729)-IRF1 and g.186643058A>G (rs5275)-PTGS, A/C-G/T of g.132484229C>A (rs2070729)-IRF1 and g.42140549G>T (rs5029748)-IKBKB as well as G/T-C/T genotypes of g.42140549G>T (rs5029748)-IKBKB and g.186640617C>T (rs4648308)-PTGS2 had a greater risk of MDD appearance.
Single-nucleotide polymorphisms of genes encoding IRF1, IKBKB, TGFA, TGFB, PTGS2 and the age of the first episode of MDD and the severity classification on the hamilton depression rating scale To estimate whether the investigated polymorphisms may had an impact on the age of the first episode of MDD, patients were stratified in accordance to genotype and their age of onset was compared (Fig. 1). A significant difference was found between A/A and A/G genotypes as well as A/A and G/G genotypes of g.41354391A>G (rs1800469)-TGFB1. Carriers of A/A genotype had their first episode significantly earlier compared to other genotypes.
In the case of the impact of genotypes of the investigated SNPs on the episode severity measured using the Hamilton Depression Rating Scale (HDRS) (Fig. 2), significant differences was found between carriers of A/A and G/G genotypes of g.41354391A>G (rs1800469)-TGFB1.
Single-nucleotide polymorphisms of genes encoding IRF1, IKBKB, TGFA, TGFB, PTGS2 and effectiveness of depression treatment We also evaluated impact of the studied polymorphisms on the effectiveness of antidepressant treatment with selective serotonin reuptake inhibitor (SSRI) (Fig. 3). Regarding the effect of investigated SNPs on treatment efficiency, differences was found between A/A and A/C genotypes of g.132484229C>A (rs2070729)-IRF1 as well as G/G and T/T genotypes of g.42140549G>T (rs5029748)-IKBKB.

DISCUSSION
There is strong amount of evidence that inflammation is undeniably associated with major depressive disorder. Moreover, it was confirmed that some inflammatory genes and presence of their genetics variants play important role in MDD development. Additionally, several loci/chromosomal regions connected with MDD were mapped by genome-wide linkage analysis, that is, 1q32. 1, 2p25.1, 3p21.1, 3p26.1, 3q26.1, 6p22.3, 8q22.2, 8q22.3, 8q12.1, 8q23.3, 11p14.2-p14.3, 13q31.1-q31.3, 15q25.2 and 19q12 (McGuffin et al., 2005;Shyn et al., 2011;Sullivan et al., 2013). Selected candidate genes in current study are located in proximity to the above mentioned regions of chromosomes. In this research, we genotyped six polymorphic variants of TGFA, TGFB1, IRF1 and PTGS2 genes; and to our knowledge, none of this SNPs have been studied in the context of severity and treatment response in depression before. However, these SNPs were included in GWAS but only one of them, that is, rs2070729, had p value below 0.05.
The first of investigated polymorphisms in this study was g.70677994G>A (rs2166975)-TGFA. The SNP is localized on 2p13.3 and it is responsible for synonymous change Val159Val. This terminal amino acid is present in the precursor protein and is necessary for glycosylation during protein maturation as well as protein localization to the cell surface (Briley et al., 1997). In our study, we were the first to show a link between rs2166975 polymorphism of TGFA and depression. The results confirmed that A/G genotype of rs2166975 is more frequently distributed in patients suffering from depression. Interestingly, the same genotype increased the risk of MDD only in man population. In the case of the gene-gene interactions between polymorphism of TGFA and other SNPs, analysis confirmed that A/G-A/C combined genotypes of rs2166975-TGFA and rs207072-IRF1 are associated with higher chance to develop MDD. In addition, A/G-G/G genotypes of rs2166975-TGFA and rs5275-PTGS2 is associated with higher risk of MDD, while G/G-G/G homozygotes decreased this chance. It was indicated that rs2166975 in TGFA gene, showed association with the risk of cleft palate (Morkūniené et al., 2007). Furthermore, another study confirmed, using transmission disequilibrium test, that minor allele of rs2166975 was over-transmitted to cleft-palate cases (Carter et al., 2010). Although there are no studies investigating role of rs2166975 polymorphism in depression or any other psychiatric disorders, our results suggest important role of investigated polymorphism in pathophysiology and course of depression. The second studied SNP, g.41354391A>G (rs1800469) of TGFB1, is located on 19q13.2 in the proximal negative regulatory region of the gene. The human TGFB1 protein is considered to be one of the immunosuppressive cytokines, which plays crucial role in CNS development (Sousa Vde et al., 2004). It is responsible for such functions as astrocyte differentiation, synaptogenesis and neuronal migration (De Sampaio e Spohr et al., 2002;Sousa Vde et al., 2004;Feng & Ko, 2008;Siegenthaler & Miller, 2004). Our results show that rs1800469 polymorphism is associated with both severity of depressive episodes and age of the onset of the disease. Precisely, carriers of G/G genotype are characterized by more severe episodes than A/A genotype carriers, which may correlate with increased concentrations of TGFB1. Moreover, a significant difference in the age of the first episode of MDD was found between A/A and A/G genotypes, as well as A/A and G/G genotypes of rs1800469-TGFB1. In accordance to our findings, TGFB levels were found to be increased in people suffering from MDD (Davami et al., 2016;Kim et al., 2007;Kim et al., 2008) as well as in Chronic HBV-Infected Patients (CHB) with mild depression symptoms (Bahramabadi et al., 2017). It has been reported that rs1800469, is not associated with neither Alzheimer's disease risk (Chang et al., 2013) nor Schizophrenia (Kapelski et al., 2015). However rs1800469 of TGFB1 is associated with altered plasma levels of TGFB1, which may modulate a susceptibility to MDD (Shah et al., 2006;Wang et al., 2008). Data suggest that, allele G is associated with lower expression of TGFB1 (Shah et al., 2006). On the other hand, another study confirmed that genotypes A/G and G/G was correlated with increased plasma TGFB1 concentrations, indicating that G allele is associated with higher production of the protein (Wang et al., 2008). It was found that other SNPs of TGFB1 could be associated with MDD. In the case of rs1800470 (codon 10), genotype T/T is significantly more frequently distributed in depressed patients (Mihailova et al., 2016). Moreover, another study revealed that C/C genotype of the same SNP is positively correlated with higher risk of depression development and more severe episodes of the disease (Caraci et al., 2012). Although TGFB1 is considered to play important role in psychoneuroimmunology, there is only few research about its association with mental disorders, and interestingly there is no other studies investigated role of mentioned rs1800469 in MDD.
In our study we also investigated whether SNPs in PTGS2 gene are involved in MDD development. As mentioned in Introduction, PTGS2 participates in inflammatory processes partly related with neurodegeneration in CNS (Minghetti, 2004). There is evidence demonstrating that rs20417 polymorphism of PTGS2 may play a role in MDD. Precisely, presence of G allele is strongly associated with increased risk of depression development (Gałecki et al., 2010). However, we have not included this polymorphism in our study. Instead, we explored g.186640617C>T (rs4648308) polymorphism located on 1q31.1. There are evidence of its involvement in depression. Precisely, allele T and C/T genotype (in positive strand allele A and G/A genotype) of mentioned SNP are associated with significantly increased risk of IFN-a-induced depression (Su et al., 2010). Part of our result are consistent with this findings, namely, we found that C/T heterozygote increased risk of MDD in woman, as well as the C allele increased this chance in man group. On the contrary, we also reported that T/T genotype carriers of this SNP are less likely to develop depression in general population. Similarly, in man group allele T was also negatively correlated with depression prevalence.
Second polymorphism of PTGS2 gene, g.186643058A>G (rs5275) located on 1q31.1, is a functional SNP, which modulates expression of PTGS2. We were first to found that allele G is connected with higher chance of MDD occurrence. Additionally, it is confirmed that this SNP is associated with severe pain in lung cancer patients. Namely, A/A and A/G (in forward strand T/T and T/C) carriers experience more severe pain than G/G carriers (Reyes-Gibby et al., 2009;Reyes-Gibby et al., 2013). However, Mendlewicz et al. (2012) found no association between PTGS2 rs5275 polymorphism and treatment response and remission of MDD. Still, there are no other studies investigated aforementioned SNPs in PTGS2 gene in context of MDD.