Skip to main content

Re-testing reported significant SNPs related to suicide in a historical high -risk isolated population from north east India



Genetic diathesis of suicide is supported by family and twin studies. Few candidate gene pathways are known, but does not explain fully the complexity of suicide genetic risk. Recent investigations opting for Genome-Wide Association Studies (GWAS) resulted in finding additional targets, but replication remained a challenge. In this respect small isolated population approach in several complex disease phenotypes is found encouraging. The present study is an attempt to re-test some of the reported significant SNPs for suicide among a small historical high- risk isolated population from Northeast India.


Two hundred ten cases (inclusive of depressed, suicide attempter and depressed + suicide attempter) and 249 controls were considered in the present study which were evaluated for the psychiatric parameters. Sixteen reported significant SNPs for suicide behaviour were re-tested using association approach under various genetic models. Networking by GeneMANIA tool was used for function prediction of the associated genes.


Seven SNPs (of 6 genes) remained significant in different genetic models. On networking genes with significant SNPs IL7, RHEB, CTNN3, KCNIP4, ARFGEF3 are found in interaction with already known candidate gene pathways while SNP rs1109089 (RHEB) gained further support from earlier expression studies. NUGGC gene is in complete isolation.


Small population approach in replicating significant SNPs is useful in complex phenotypes like suicide. This study explored the region-specific demographics of India by identifying vulnerable population for suicide via genetic association analysis in bringing into academic and administrative forum, the importance of suicide as a disease and its biological basis.


Suicide represents a wide range of risk factors encompassing throughout the life period of the subjects. Recently, suicide is considered as a major public health issue and is the second leading cause of death globally among the age groups of 15 to 29 years [1]. Almost 800,000 suicides are completed every year and in this 79% are documented from the low and middle-income countries alone [2]. Suicide death rate in India from 1990 to 2016 was estimated at 17.9 per 100,000 persons, which equates to around 230,000 suicides annually [3]. The complex etiology of suicide involves interaction of genetics with other psychiatric, neurological and environmental conditions. Twin studies demonstrated that suicide has a genetic component and familial basis is evident and is now proven that greater risk of a suicide attempt is involved in the offspring of a candidate with a positive history of completed suicide [4]. Besides, the increased rate of suicide attempt is found in the first-degree relatives of suicide probands [5]. However, the daunting task of complexity lies in the partial genetic contribution of other associated psychiatric disorders such as mood, alcohol/ substance use, schizophrenia etc. along with independent heritable factors for impulsivity, aggression [6]. As of now, few pathways and candidate genes in the central nervous system, serotonergic hypofunction and impaired negative feedback of the hypothalamic-pituitary-adrenal axis are frequently observed in both attempters and those who die by suicide [7] . The contribution of these candidate gene association studies underscores the complexity of suicide genetic risk, leading to several recent investigations opting for hypothesis-free methods like Genome-wide association studies (GWAS) of common variants [8,9,10,11,12,13]. As a result, many known and novel variants have been found [14, 15]. However, replication of the suggestive significant GWAS SNPs in other cohorts remains a challenge [8, 13]. Factors such as sample size [16], phenotype cohort [17], correlation with family history [18], subject follow-ups, lifestyle and cultural differences from population to the population sampled [19], all these together play a vital role in failing to replicate the significance of a GWAS SNP in the pathogenesis of suicide [8, 18]. On the other hand, small isolated populations are known to be useful in replication, as they yield better results even if the sample size is relatively small as compared to large sample size requirements in heterogeneous large populations. Besides, similar co-factors (environment, lifestyle, geography, ethnicity) between case/control holds a benefit along with the advantage of founder effect and reduced genetic heterogeneity in small isolated populations, which is a challenge in general population screening even if sample numbers are large [20]. It is also observed that due to founder effect the frequency of phenotypic traits of a complex disease is likely to be high in small populations [21].

The Idu-Mishmi is an isolated small endogamous Tibeto-Burman speaking tribe inhabiting Lower and Upper Dibang Valleys in Arunachal Pradesh, India and numbered around 15,000 souls in 2001 census [22]. Mene [23,24,25] reported that, between 1971 and 2010, 218 cases of suicide occurred in this tribe in the 10–29 age group with an estimated suicide rate of 58 per 100,000 individuals surpassing the national average. In a systematic sampling and phenotypic annotation, our studies on Idu- Mishmi reported a high rate of attempted suicide (14.2% compared to the general urban population frequency of 0.4–4.2%) and significant association of depression and endo-phenotypes (impulsivity and aggression) [26,27,28]. The key objective of the present study is whether suggestive significant SNPs reported in GWAS studies conducted in 10 years period (2004 to 2015) and published in a comprehensive review [8], targeting suicide behavior can be replicated in the historical isolated small Idu-Mishmi population with a high rate of suicide attempt.


The study subjects recruited were based on their psychiatric trait assessment of suicide attempt, depression and family history. The controls were negative for both psychiatric traits screened and family history. Idu-Mishmi is a close-knit community distributed over a few settlements and suicide occurrences are a shared memory of 3 to 4 generations accordingly, families can be identified. To control the relatedness in a small endogamous population sampling was done on critical screening where primary and secondary degree relatives were not included in the study. The psychiatric assessment of suicide behavior was done based on the Columbia Suicide Severity Rating Scale (C-SSRS) and PHQ-9 used for depression assessment. The reliability and validity of these scales used were described in our study [26, 27]. All the studied individuals with a suicide attempt, depression or both with or without positive family history have been considered as cases (N = 210) and 249 healthy controls were recruited without any of the given traits and negative family history. The sample numbers for genotyping vary between SNPs due to degradation in storage and transportation of blood samples from remote areas.

The suggestive SNPs with significant threshold at p < 0.001 [8] were selected and included in the present study. Genotyping in the study population was done with ARMS PCR Technique and SNP location traced based on coordinate position (GRCh38.p7) followed by sequence retrieval from ensemble genome browser []. Formatting of the retrieved sequences was done using ApE (v2.0.55). []. Forward and reverse (normal and mutant) primers were designed using Primer3Plus []. To increase the specificity, weak secondary mismatches were incorporated at the penultimate position manually for a strong primary mismatch and a strong secondary mismatch was introduced for a weak primary mismatch for prevention of any false positives as two mismatches PCR will not continue and in case of one mismatch, the PCR will be initiated [29]. All the 16 primer sequences used are shown in Table S1. Each primer was standardized for specific Tm and specificity using gradient PCR then validation PCR was carried out on 20 samples to check for specificity and Tm. To ensure the selectivity and specificity re-validation of same samples were done twice with similar conditions to ensure similar genotypes in each run. In-silico PCR was also carried out for all the primers and expected band sizes were carefully observed for the exclusion of any false positives. Extreme care was taken in designing and validation, after a complete testing only validated primers were used for genotyping by ARMS technique for all the cases and controls. Genotyping was performed using agarose gel electrophoresis on a 2% gel stained with ethidium bromide and visualized in T-Genius, Syngene gel documentation system.

Allele frequencies were calculated by gene counting method and HW equilibrium test by on line portal ( Differences in age and sex among cases and controls were presented as percentages with Z test and p- values. Inheritance models (dominant, recessive and additive) were employed for the categorization and interpretation of the data. Age, sex-adjusted bivariate analyses were performed using the SPSS package (v16.0, SPSS Chicago). Multiple testing by Bonferroni correction was performed on the significant SNPs to remove false positives, the critical value was set to 0.05 and the number of tests to 25 and the corrected critical value was kept 0.002. Network analyses were done by GeneMANIA (, a user- friendly flexible web site that uses a wealth of genomics and proteomics data and finds functionally similar genes with a gene list query [30]. The default parameters of the network with a percentage contribution of each are as given in Table S4. This study was approved by the Institutional Ethical Board of the Department of Anthropology, Delhi University and all subjects participated in the study gave their written informed consent.


Table 1 shows age, sex distribution in cases and controls of the study. Overall females are at slightly higher risk (p = 0.08186), however between age groups > 19 years show higher risk in both sexes, with significant difference in males (p = 0.0004). Information on each of the selected SNPs, their genomic context within the respective genes, associated psychiatry traits reported and minor allele frequencies (MAF) of the present population in comparison to other continental populations available from European, South East Asian and South Asian populations are given in Table S2. The variation of MAF frequency in controls is within the range reported for East Asian and South Asian populations. Genotypes and allele frequencies of the 16 selected SNPs with HW test of significance along with a percentage of the difference between observed and expected heterozygotes are shown in Table S3. Ten of 16 SNPs are with significant HW deviation among controls with low observed heterozygote frequencies. Bivariate analyses with age-sex adjusted ODDs ratios in dominant, recessive and additive genetic models were computed and 7 SNPs (of 6 genes) remained significant with Bonferroni corrected p-values (Table 2). Significant SNPs of genes RHEB, NUGCC are at risk in the dominant model, IL7 in the recessive model, ARFGEE3, KCNIP4 in the additive model and CTNN3 in dominant and recessive models. However, simultaneous protection in alternative models is also seen. All are at a higher threshold of significance (p < 2.50E− 08), whereas ARFGEE3 gene SNP at p < 0.035. Figure S1 shows the networking of 6 genes with significant found SNPs in our population along with candidate genes of serotonergic system (TPH1, TPH2, HTR1A, HTR1B, HTR2A, SLC6A4), dopaminergic and adrenergic system (DRD2, AKT1, AKTIP, ADRA2), catabolism of monoamines system (COMT, MAOA), HPA axis pathway (CRHR1, CRHR2, FKBP5, CRHBP, NR3C1, AVPR1B) and neurotrophic processors (BDNF, TRKB, CCKBR, NGF, NTRKR2, HOMER1, NPTX2). Five out of 6 genes were observationally found interacting/closer with specific candidate gene pathways, i.e. KCNIP4 (Neurotrophic processors), RHEB (HPA axis), IL7& CTNNA3 (Serotonergic), ARFGEF3 (Dopaminergic and Adrenergic). However, NUGGC gene is having no interaction with any of the genes in the network and is lying in distant isolation.

Table 1 Age-sex distribution among cases and controls along with Z-test values
Table 2 Association analysis of SNPs in different genetic models (age and sex adjusted) between cases and controls with Bonferroni corrected p -values


The present study is born out of the fact that significant SNPs in several GWAS studies conducted in various populations, failed to replicate in other populations/cohorts because of various factors like large sample size requirement, phenotypic heterogeneity and lifestyle cultural differences etc. This remained a major challenge in understanding the pathogenesis of suicide, in-spite of the fact that these SNPs may be important and may provide an opportunity for discovery of novel bio-pathological pathways or strengthening known pathways in biomarker discovery for suicide risk and therapeutic intervention. Founder effect events in historically isolated small populations are useful in finding genes not only in Mendelian disorders but also in polygenes of several complex disease phenotypes with the advantage of narrowing down on the sub-phenotype heterogeneity [31,32,33,34]. To address this issue, we designed the present study in a historical small isolated endogamous Idu Mishmi population having the highest rate of suicide attempt (14.2%) compared to general urban population (0.4–4.2%) with depression as a significant covariate described in our earlier studies [26, 27]. Genetic variants identified in a small population are not restricted only by founder effects, but they can also be mapped in larger populations that help in identification of new or strengthening known pathways underlying the effect of these SNPs in other complex diseases [20]. The Human genome has been explored to have around 10 million variants/SNPs varying individually [35], therefore it becomes a strong point to consider GWAS SNPs in the pathobiology of suicide. Some of the reported suggestive significant SNPs for suicide were chosen for re-testing in our small high-risk isolated population association study on the subjects recruited after careful evaluation of psychiatric traits, suicide attempt and depression, and we found several SNPs were highly significant.

In pathway analysis by GeneMANIA ( the suggestive significant SNP genes are associated with specific biological networks (pathway, co-expression, shared protein domains, physical interaction, co-localization, genetic interaction) of the candidate genes, observationally. RHEB, CTNNA3, KCNIP4, IL7, ARFGEF3 were interacting with the candidate genes via connecting genes/pathways. RHEB gene product is a GTP-binding protein called RAS homolog enriched in the brain. The main function of this gene is involved in mTOR pathway and regulation of the cell cycle. The SNP of this gene (rs1109089) which has been re-tested and found significant in the present study, has already been implicated in suicide [36]. mTOR pathway is well known for its antidepressant drug response [37,38,39]. Emerging quick and effective antidepressants are now available based on RHEB mediated mTOR pathway which is showing an effect on treatment-resistant subjects [40, 41]. CTNNA3 gene has roles implicated in the formation of stretch-resistant cell-cell adhesion complexes and is reported for causing mental ailments such as schizophrenia [42], besides the gene is reported to be expressed in the cerebellum [43]. KCNIP4 gene is found significant for various drug targets being a Kv channel-interacting protein 4 family [44]. IL7 gene in the immune system plays a role in psychiatric disorders is a fact and a lot of research claims this point, along with other cell types, it is also produced in neurons and it has close interconnections with serotonergic pathways via IL9 as shown in figure S1. A study reported high expression levels of IL7 in affected males and low levels in affected females [45]. ARFGEF3 gene showed connection via a connecting gene to the candidate gene DRD2 of a dopaminergic and adrenergic pathway. The genes with significant variants found suggestive of their role in the pathogenicity of suicide attempt and depression in the present study. However, more insights are required to explore functional validation of the role of these Genes/SNPs in causing the phenotypes considered.

The most significant observation of the study is NUGGC gene (Nuclear GTPase SLIP- GC) with significant risk allele in our bivariate analysis, is an outlier in the network (Figure S1). NUGGC inhibits function of the activation-induced cytidine deaminase AICDA [46] and helps in maintenance of genome stability by reduction of somatic hyper-mutations in B-Cells [47]. The functional validation of this gene may lead to discovery of entirely new bio-pathway in suicide research.

Conclusion and future prospects

Worldwide suicide is considered as the leading cause of death. This study explores the region-specific demographic assessment of the suicidal genetic risk. Re-testing of reported significant SNPs in our study population and functional correlation with known candidate gene bio-pathways contributed to the advantage of historically isolated populations in deciphering genetics of complex phenotypes like suicide. Having such a small population with high prevalence of the targeted phenotype is an added advantage. Besides, the finding of involvement of a gene related to genome stability may be important in finding entirely novel bio-pathways in future suicide research. However, the findings of this study need to be considered as exploratory and further functional validations are definitely in need.

Limitations: Subject selection is a challenge in a small population genetics association study. Significant HW deviation is a concern, low observed heterozygote frequencies are as expected in local small endogamous population. However, 3 genes with suggestive significant SNPs re-tested in our study were in HW equilibrium for controls. In studies aiming re-testing or replication of previous studies, one of the serious limitations is different methods of phenotypic annotation used, which also applies to the present study.

Availability of data and materials

All the analysis data will be made available upon request.



Genome Wide Association Study


Single Nucleotide Polymorphism


Columbia Suicide Severity Rating Scale


Patient Health Quesstionaire-9


Amplification Refractory Mutation System


Polymerase Chain Reaction


Melting Temperature




Minor Allele Frequency


  1. 1.

    Bachmann S. Epidemiology of suicide and the psychiatric perspective. Int J Environ Res Public Health. 2018;15(7):1425 Multidisciplinary Digital Publishing Institute (MDPI); [cited 2019 Feb 22]. Available from:

    Article  Google Scholar 

  2. 2.

    Suicide: one person dies every 40 seconds. [cited 2019 Oct 30]. Available from:

  3. 3.

    Armstrong G, Vijayakumar L. Suicide in India: a complex public health tragedy in need of a plan; 2018. [Cited 2019 Feb 22]; Available from:

    Google Scholar 

  4. 4.

    Brent DA, Oquendo M, Birmaher B, Greenhill L, Kolko D, Stanley B, et al. Peripubertal suicide attempts in offspring of suicide attempters with siblings concordant for suicidal behaviour. Am J Psychiatry. 2003;160:1486–93 [Cited 2020 Apr 12]. Available from:

    Article  Google Scholar 

  5. 5.

    Melhem NM, Brent DA, Ziegler M, Iyengar S, Kolko D, Oquendo M, et al. Familial pathways to early-onset suicidal behaviour: familial and individual antecedents of suicidal behaviour. Am J Psychiatry. 2007;164:1364–70.

    Article  Google Scholar 

  6. 6.

    Mann JJ, Arango VA, Avenevoli S, Brent DA, Champagne FA, Clayton P, et al. Candidate Endophenotypes for genetic studies of suicidal behaviour. Biol Psychiatry. 2009;65:556–63 [Cited 2019 July 22]; Available from:

    CAS  Article  Google Scholar 

  7. 7.

    Zai CC, de Luca V, Strauss J, Tong RP, Sakinofsky I, Kennedy JL. Genetic factors and suicidal behaviour. Neurobiol Basis Suicide. 2012; CRC Press/Taylor & Francis; [Cited 2019 Mar 25]. Available from:

  8. 8.

    Mirkovic B, Laurent C, Podlipski M-A, Frebourg T, Cohen D, Gerardin P. Genetic association studies of suicidal behaviour: a review of the past 10 years, Progress, limitations, and future directions. Front Psychiatry Front Media SA. 2016;7:158 [Cited 2019 Mar 25]. Available from:

    Google Scholar 

  9. 9.

    Perlis RH, Huang J, Purcell S, Fava M, Rush AJ, Sullivan PF, et al. Genome-wide association study of suicide attempts in mood disorder patients. Am J Psychiatry. 2010;167:1499–507 [Cited 2019 Jul 22]. Available from:

    Article  Google Scholar 

  10. 10.

    Schosser A, Butler AW, Ising M, Perroud N, Uher R, Ng MY, et al. Genomewide association scan of suicidal thoughts and behaviour in major depression. PLoS One. 2011;6:e20690.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  11. 11.

    Willour VL, Seifuddin F, Mahon PB, Jancic D, Pirooznia M, Steele J, et al. A genome-wide association study of attempted suicide. Mol Psychiatry. 2012;17:433–44 [Cited 2019 Jul 22]. Available from:

    CAS  Article  Google Scholar 

  12. 12.

    Mullins N, Perroud N, Uher R, Butler AW, Cohen-Woods S, Rivera M, et al. Genetic relationships between suicide attempts, suicidal ideation and major psychiatric disorders: a genome-wide association and polygenic scoring study. Am J Med Genet Part B Neuropsychiatr Genet. 2014;165:428–37 [Cited 2019 May 8]. Available from:

    Article  Google Scholar 

  13. 13.

    Galfalvy H, Haghighi F, Hodgkinson C, Goldman D, Oquendo MA, Burke A, et al. A genome-wide association study of suicidal behaviour. Am J Med Genet Part B Neuropsychiatr Genet. 2015;168:557–63 [Cited 2019 Mar 25]. Available from:

    CAS  Article  Google Scholar 

  14. 14.

    Mullins N, Bigdeli TB, Børglum AD, Coleman JRI, Demontis D, Mehta D, et al. GWAS of suicide attempt in psychiatric disorders and association with major depression polygenic risk scores. Am J PsychiatryAmerican Psychiatric Association. 2019;176:651–60.

    Article  Google Scholar 

  15. 15.

    Levey DF, Polimanti R, Cheng Z, Zhou H, Nuñez YZ, Jain S, et al. Genetic associations with suicide attempt severity and genetic overlap with major depression. Transl Psychiatry. 2019;9:22 [Cited 2019 Mar 25]. Available from:

    Article  Google Scholar 

  16. 16.

    Ordaz SJ, Goyer MS, Ho TC, Singh MK, Gotlib IH. Network basis of suicidal ideation in depressed adolescents. J Affect Disord. 2018;226:92–9 [Cited 2019 Mar 25]. Available from:

    Article  Google Scholar 

  17. 17.

    Sher L, Grunebaum MF, Burke AK, Chaudhury S, Mann JJ, Oquendo MA. Depressed multiple-suicide­attempters – a high-risk phenotype. Crisis. 2017;38:367–75 [Cited 2019 Mar 25]. Available from:

    Article  Google Scholar 

  18. 18.

    Lizardi D, Sher L, Sullivan GM, Stanley B, Burke A, Oquendo MA. Association between familial suicidal behaviour and frequency of attempts among depressed suicide attempters. Acta Psychiatr Scand. 2009;119:406–10 NIH Public Access; [Cited 2019 Mar 25]. Available from:

    CAS  Article  Google Scholar 

  19. 19.

    Perez-Rodriguez MM, Baca-Garcia E, Oquendo MA, Blanco C. Ethnic differences in suicidal ideation and attempts. Prim Psychiatry. 2008;15:44–53 NIH Public Access; [Cited 2019 Mar 25]. Available from:

    PubMed  PubMed Central  Google Scholar 

  20. 20.

    Kristiansson K, Naukkarinen J, Peltonen L. Isolated populations and complex disease gene identification. Genome Biol BioMed Central. 2008;9:109 [Cited 2019 Jul 17]. Available from:

    Article  Google Scholar 

  21. 21.

    Moltke I, Grarup N, Jørgensen ME, Bjerregaard P, Treebak JT, Fumagalli M, et al. A common Greenlandic TBC1D4 variant confers muscle insulin resistance and type 2 diabetes. Nature. 2014;512:190–3 Nature Publishing Group; [cited 2019 Jul 22]. Available from:

    CAS  Article  Google Scholar 

  22. 22.

    Baruah TKM. The Idu Mishmi. India: Government of Assam; 1960.

    Google Scholar 

  23. 23.

    Mene T. Suicide among the Idu Mishmi tribe of Arunachal Pradesh [doctoral dissertation]. Itanagar: Rajiv Gandhi University; 2011.

    Google Scholar 

  24. 24.

    Mene T. Underestimation of suicide deaths: a study of the Idu Mishmi tribe of Arunachal Pradesh. Econ Polit Wkly. 2013;52:129–33 Available from:

    Google Scholar 

  25. 25.

    Mene T. Suicide: a study of the Idu Mishmi tribe of Arunachal Pradesh. Resarun J Directorate ResGovernment of Arunachal Pradesh. India. 2013;37:20–6.

    Google Scholar 

  26. 26.

    Singh PK, Singh RK, Biswas A, Rao VR. High rate of suicide attempt and associated psychological traits in an isolated tribal population of north-East India. J Affect Disord. 2013;151:673–8 [Cited 2019 Mar 25]. Available from:

    Article  Google Scholar 

  27. 27.

    Singh PK, Rao VR. Explaining suicide attempt with personality traits of aggression and impulsivity in a high risk tribal population of India. Pendyala G, editor. PLoS One. 2018;13:e0192969Public Library of Science; [Cited 2019 Mar 18]. Available from.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  28. 28.

    Singh PK. Molecular anthropological study of depression and suicide in Idu-Mishmi tribe of Arunachal Pradesh [doctoral dissertation]. Delhi: [Department of Anthropology]: University of Delhi; 2015.

    Google Scholar 

  29. 29.

    Little S. Amplification-refractory mutation system (ARMS) analysis of point mutations. Curr Protoc hum genet. Hoboken: Wiley; 2001. [Cited 2019 Mar 25]. p. Unit 9.8. Available from:

    Google Scholar 

  30. 30.

    GeneMANIA update 2018. [Cited 2020 May 13]. Available from:

  31. 31.

    Heutink P, Oostra BA. Gene finding in genetically isolated populations. Hum Mol Genet. 2002;11:2507–15 Narnia; [Cited 2019 Jul 17]. Available from:

    CAS  Article  Google Scholar 

  32. 32.

    Liu F, Arias-Vásquez A, Sleegers K, Aulchenko YS, Kayser M, Sanchez-Juan P, et al. A Genomewide screen for late-onset Alzheimer disease in a genetically isolated Dutch population. Am J Hum Genet. 2007;81:17–31 [Cited 2019 Jul 17]. Available from:

    CAS  Article  Google Scholar 

  33. 33.

    Lowe JK, Maller JB, Pe’er I, Neale BM, Salit J, Kenny EE, et al. Genome-wide association studies in an isolated founder population from the Pacific Island of Kosrae. Gibson G, editor. PLoS Genet. 2009;5:e1000365 [Cited 2019 Jul 17]. Available from:

    Article  Google Scholar 

  34. 34.

    Hatzikotoulas K, Gilly A, Zeggini E. Using population isolates in genetic association studies. Brief Funct Genomics. 2014;13:371–7 [Cited 2019 Jul 17]. Available from:

    Article  Google Scholar 

  35. 35.

    1000 Genomes Project Consortium T 1000 GP, Abecasis GR, Altshuler D, Auton A, Brooks LD, Durbin RM, et al. A map of human genome variation from population-scale sequencing. Nature. 2010;467:1061–73 Europe PMC Funders; [Cited 2019 Mar 25]. Available from:

    Article  Google Scholar 

  36. 36.

    Menke A, Domschke K, Czamara D, Klengel T, Hennings J, Lucae S, et al. Genome-wide association study of antidepressant treatment-emergent suicidal ideation. Neuropsychopharmacology. Nature Publishing Group; 2012;37:797–807. [Cited 2019 Apr 22]. Available from:

  37. 37.

    Rodrigues AL, Réus G, Quevedo J. mTOR signaling in the neuropathophysiology of depression: current evidence. J Receptor Ligand Channel Res. 2015;8:65–74 Dove Press; [Cited 2019 May 21]. Available from:

    Article  Google Scholar 

  38. 38.

    Harraz MM, Tyagi R, Cortés P, Snyder SH. Antidepressant action of ketamine via mTOR is mediated by inhibition of nitrergic Rheb degradation. Mol Psychiatry. 2016;21:313–9 NIH Public Access; [Cited 2019 May 21]. Available from:

    CAS  Article  Google Scholar 

  39. 39.

    Flory JD, Donohue D, Muhie S, Yang R, Miller SA, Hammamieh R, et al. Gene expression associated with suicide attempts in US veterans. Transl Psychiatry. 2017;7:e1226 Nature Publishing Group; [Cited 2019 May 21]. Available from:

    CAS  Article  Google Scholar 

  40. 40.

    Duman RS. Ketamine and rapid-acting antidepressants: a new era in the battle against depression and suicide. F1000Research. 2018;7:1 Faculty of 1000 Ltd; [Cited 2019 May 21]. Available from:

    Article  Google Scholar 

  41. 41.

    Niculescu AB, Le-Niculescu H, Levey DF, Phalen PL, Dainton HL, Roseberry K, et al. Precision medicine for suicidality: from universality to subtypes and personalization. Mol Psychiatry. 2017;22:1250–73 Nature Publishing Group.

    CAS  Article  Google Scholar 

  42. 42.

    Uher R. Gene-environment interactions in severe mental illness. Front Psychiatry. 2014;5:48 Frontiers Research Foundation; Available from:

    Article  Google Scholar 

  43. 43.

    Liu W, Yan H, Zhou D, et al. The depression GWAS risk allele predicts smaller cerebellar gray matter volume and reduced SIRT1 mRNA expression in Chinese population. Transl Psychiatry. 2019;9:333.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  44. 44.

    Perroud N, Uher R, Ng MYM, Guipponi M, Hauser J, Henigsberg N, et al. Genome-wide association study of increasing suicidal ideation during antidepressant treatment in the GENDEP project. Pharm J. 2012;12:68–77 Nature Publishing Group; [Cited 2020 Apr 7]. Available from:

    CAS  Google Scholar 

  45. 45.

    Hall JR, Wiechmann A, Edwards M, Johnson LA, O’Bryant SE. IL-7 and depression: the importance of gender and blood fraction. Behav Brain Res Elsevier BV. 2016;315:147–9.

    CAS  Article  Google Scholar 

  46. 46.

    Richter K, Brar S, Ray M, Pisitkun P, Bolland S, Verkoczy L, et al. Speckled-like pattern in the germinal center (SLIP-GC), a nuclear GTPase expressed in activation-induced deaminase-expressing lymphomas and germinal center B cells. J Biol Chem. 2009;284:30652–61 [Cited 2019 Jun 9]. Available from:

    CAS  Article  Google Scholar 

  47. 47.

    Richter K, Burch L, Chao F, Henke D, Jiang C, Daly J, et al. Altered pattern of immunoglobulin Hypermutation in mice deficient in Slip-GC protein. J Biol Chem. 2012;287:31856–65 [Cited 2019 Jun 9]. Available from:

    CAS  Article  Google Scholar 

Download references


The authors acknowledge Genome Foundation, for facilitation of Lab and equipment’s for execution of experiments and patients and their family for participating in the study.

Informed consent

All subjects (human) participated in the study gave their written consent and before taking the consent the details of the study were explained to participants in their native language.


This work was supported by Delhi University Faculty Grant and ICMR Emeritus Medical Scientist Fellowship (ICMR/74/1/2016- Pers, EMS) to VRR and Project Assistant support to GG. The funding source has no involvement in study design; in the collection, analysis and interpretation of data; in the writing of the report; and in the decision to submit an article for publication.

Author information




Authors Name: Prof. V. R. Rao, Contribution: Project design, Funding acquisition & Implementation, manuscript finalization. Dr. Piyoosh Kumar Singh, Contribution: Participated in project design and paper writing, involved in field operations and data collection, curation, analysis. Mr. Gaurav Gupta, Contribution: Laboratory experiments, data curation, analysis and paper writing. Dr. Ravi Deval, Contribution: Data curation, analysis. Mr. Shashank Upadhyay, Contribution: Software applications, analysis. Dr. Anshuman Mishra, Contribution: Manuscript Finalization and Updation. All authors read and approved the final manuscript

Corresponding author

Correspondence to V. R. Rao.

Ethics declarations

Ethics approval and consent to participate

All procedures performed in the study involving human participants were in accordance with the ethical standards and was approved by the Institutional Ethical Board, Department of Anthropology, Delhi University.

Consent for publication

All subjects gave their consent for conducting research on their samples and publishing the same. The personal identity of any subject is not disclosed in the manuscript.

Competing interests

No Competing Interest.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary information

Additional file 1: Table S1.

Primer sequences used for genotyping along with melting temperature and amplicon base pairs. Table S2. List of suggestive significant SNPs chosen from a review of 10 year study for re-testing in the present study. Table S3. Distribution of genotypes, allele frequencies and HW test of significance among cases and controls of the 16 selected SNPs in the present study. Table S4. Network pathway percentage contribution of different default parameters of the figure S1. Figure S1. Pathway network analysis of 6 significant genes of the present study along with candidate genes of different known pathways of suicide behaviour. Default settings used for creation of network are as described in Table S4.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Gupta, G., Deval, R., Mishra, A. et al. Re-testing reported significant SNPs related to suicide in a historical high -risk isolated population from north east India. Hereditas 157, 31 (2020).

Download citation


  • Suicide
  • GWAS
  • SNPs
  • Replication