|Year : 2014 | Volume
| Issue : 6 | Page : 736-743
5-Hydroxytryptamine (serotonin) 2A receptor gene polymorphism is associated with schizophrenia
Subash Padmajeya Sujitha1, Asha Nair1, Moinak Banerjee2, Srinivasan Lakshmanan3, Sampth Harshavaradhan4, Soosiah Gunasekaran1, Anilkumar Gopinathan1
1 School of Biosciences & Technology, VIT University, Vellore, India
2 Rajiv Gandhi Centre for Biotechnology, Thiruvananthapuram, India
3 Government Vellore Medical College, Vellore, India
4 Sri Narayani Hospital & Research Centre, Vellore, India
|Date of Submission||27-Aug-2013|
|Date of Web Publication||3-Mar-2015|
Professor, School of Biological Sciences, VIT University, Vellore 632 014, Tamil Nadu
Source of Support: None, Conflict of Interest: None
| Abstract|| |
Background & objectives: Schizophrenia, the debilitating neuropsychiatric disorder, is known to be heritable, involving complex genetic mechanisms. Several chromosomal regions associated with schizophrenia have been identified during the past; putative gene (s) in question, to be called the global signature for the pathophysiology of the disease, however, seems to evade us. The results obtained from the several population-wise association-non association studies have been diverse. w0 e therefore, undertook the present study on Tamil speaking population in south India to examine the association between the single nucleotide polymorphisms (SNPs) at the serotonin receptor gene (5HT2A) and the occurrence of the disease.
Methods: Blood samples collected from 266 cases and 272 controls were subjected to genotyping (PCR amplification of candidate SNPs, RFLP and sequencing). The data on the SNPs were subjected to statistical analysis for assessing the gene frequencies in both the cases and the controls.
Results: The study revealed significant association between the genotypic frequencies of the serotonin receptor polymorphism and schizophrenia. SNP analysis revealed that the frequencies of GG (30%, rs6311) and CC genotypes (32%, rs6313), were higher in patients (P<0.05) than in controls. The study also showed presence of G and C alleles in patients. s0 ignificant levels of linkage disequilibrium (LD) were found to exist between the genotype frequencies of rs6311 and rs6313.
Interpretation & conclusions: This study indicated an association between the SNPs (rs6311 and rs6313) of the serotonin receptor 5HT2A and schizophrenia. HapMap analysis revealed that in its genotype distribution, the Tamil speaking population was different from several other populations across the world, signifying the importance of such ethnicity-based studies to improve our understanding of this complex disease.
Keywords: GC/CC homozygocity - 5HT2A receptor - linkage disequilibrium - rs6313 - schizophrenia
|How to cite this article:|
Sujitha SP, Nair A, Banerjee M, Lakshmanan S, Harshavaradhan S, Gunasekaran S, Gopinathan A. 5-Hydroxytryptamine (serotonin) 2A receptor gene polymorphism is associated with schizophrenia. Indian J Med Res 2014;140:736-43
|How to cite this URL:|
Sujitha SP, Nair A, Banerjee M, Lakshmanan S, Harshavaradhan S, Gunasekaran S, Gopinathan A. 5-Hydroxytryptamine (serotonin) 2A receptor gene polymorphism is associated with schizophrenia. Indian J Med Res [serial online] 2014 [cited 2020 Feb 26];140:736-43. Available from: http://www.ijmr.org.in/text.asp?2014/140/6/736/152447
Schizophrenia is a debilitating mental disease, the aetiology of which is still enigmatic. Investigations involving a large body of data obtained from families and twins reveal the involvement of complex genetic mechanisms leading to the high degree of heritability of the disease , . Researchers during the past two decades were able to identify several chromosomal regions associated with schizophrenia  . But so far, no single candidate gene locus or genetic variant that can be considered as a global signature for the pathophysiology of the disease, has been identified.
s0 chizophrenia has been implicated to be caused by abnormalities in neurotransmitters and its pathways. Previous studies have resulted in several lines of evidences suggesting the involvement of serotoninergic system in the onset of the disease , . Whole - genome linkage analyses have further paved the way for the identification of chromosomal regions that harbour the genes related to schizophrenia  . Lin and coworkers  in their pioneering work involving a mixed population of UK and Japanese families, have suggested a linkage between the 13q14 - q21 and schizophrenia. From then on, this part of the chromosome that houses the 5-hydroxytryptamine (5HT) receptors became the focus of attention of several investigators and this locus was suggested to enable the serotonergic system to regulate diverse neuroendocrine functions as well as depression, anxiety, and fatigue , . Subsequent attempts have provided clues on the possible involvement of two single nucleotide polymorphisms (SNPs rs6311 and rs6313) from the same gene locus (encoding 5-HT receptor, 5HT2A), in the pathogenesis of multiple neuropsychiatric disorders  . The SNP at the promoter region (rs6311) has been shown to affect the promoter activity of the gene, leading to functional consequences. The inadequacy of the promoter activity (at rs6311) has been suggested to cause neuropsychiatric disorders including schizophrenia in several populations , . However, in spite of many studies published including meta-analyses, the results obtained so far have been diverse. A combined analysis involving the Russians, Tatars and Bhashkirs, revealed a significant association to exist between the genotype frequencies at rs6311 and schizophrenia in the Russian population; the Tatars and the Bashkirs, however, did not show any significant association between the SNP and schizophrenia  . In another association study involving rs6313, a significant allelic heterogeneity was found to exist between the East Asian and European populations  . In the East Asian population, including the Chinese, Japanese, Singaporean and Korean populations, a prevalence of T allele could be noticed associated with the SNP (rs6313) and schizophrenia , . Several other studies, however, found no association between the expression of either C or T allele (rs6313) and the disease , . In view of the central role played by the serotoninergic system in neural transmission and maintaining the neuropsychological balance, and the ethnicity-based discrepancies observed in various association studies, the present study was carried out to understand the association (if any) in between the serotonin receptor (5HT2A) gene variants and schizophrenia in Tamil speaking population in south India.
| Material & Methods|| |
This study of Tamil speaking population residing in and around Vellore (Tamil Nadu), included 266 unrelated case subjects consecutively selected, diagnosed of schizophrenia from the Psychiatry out-patient department, Government Vellore Medical College, and the Sri Narayani Hospital and Research Center (Vellore) from June 2011 to December 2012, with the following inclusion and exclusion criteria.
Inclusion criteria: those who were diagnosed as schizophrenia patients as per the Diagnostic and Statistical Manual of Mental Disorders - Fourth Edition  , the case subjects, irrespective of drug naïve or medicated, and age ranging between 20 and 60 yr, of Tamil Nadu origin (Vellore district, India). The study population comprised ancestrally Tamil - speaking population only; information on their ethnicity was obtained through interviews and profiling (demographic data) forms.
Exclusion criteria: it was ensured that the matched control subjects 272 numbers with respect to age and ethnicity from among healthy volunteers working in VIT university chosen had neither any family history of any psychotic problems or symptoms, nor had they taken any neuroleptics in the past. Written informed consent was obtained from each subject prior to study. Necessary Ethical Clearances were obtained from the respective institutions [Institutional Review Board (IRB), Government Vellore Medical College, Vellore, Institutional Ethics Committee (IEC)/Institutional Review Board (IRB), Sri Narayani Hospital and Research Centre, Vellore, and University Human Ethical Committee, VIT University, Vellore]. The demographic details collected from patients and controls included the age, gender, marital status, mother tongue, status of current medication, age of onset of the symptoms, family history (if any) of schizophrenia and substance abuse (if applicable) ([Table 1]).
b0 lood samples (10 ml) were collected from the subjects in EDTA-coated tubes and kept at 4 o C and were transported immediately to the Molecular e0 ndocrinology laboratory (VIT University, Vellore) and subsequently processed for DNA isolation.
Genotyping: Genomic DNA was isolated from the leucocytes using phenol-chloroform method  . The polymorphisms (rs6311 and rs6313) of 5HT2A receptor gene were amplified in Applied Biosystems 9902 thermal cycler, USA in 15 μl reaction volume [0.5 μl of 100 μM primer (IDT, USA), 0.5 μl of 10 mM of dNTP mix (Merck), 1.0 μl of 10X Taq buffer A (Merck), 1.0 μl of Taq polymerase (Merck) and 50 ng of template DNA (Merck)]; primers , used for the amplification are shown in [Table 2]. Identical amplification protocols were followed for both the polymorphisms: i.e. initial denaturation at 94°C for 8 min, followed by 40 cycles of denaturation at 94°C for 30 sec, annealing at 55° C for 30 sec, chain lengthening at 72°C for 30 sec, and final extension at 72° C for 10 min.
|Table 2: Primer sequences and the annealing temperatures adopted for genotyping of 5HT2A variants (rs6311 and rs6313) |
Click here to view
The PCR amplicons were visualized by resolving in 1 per cent agarose gel ([Figure 1]). The PCR products were also subjected to restriction fragment length polymorphism (RFLP) using MSPI restriction enzyme (NEB, England), as per manufacturer's protocol. The digested products were resolved in 2.5 per cent agarose gel.
|Figure 1: PCR amplicons of 5HT2A SNPs. M, 100 bp marker; Lanes 1 - 3, PCR amplicons of rs6311; Lanes 4-6, PCR amplicons of rs6313.|
Click here to view
PCR amplification resulted in the appearance of distinct bands with amplicon size of 365 bp and 446 bp for rs6311 and rs6313, respectively ([Figure 1]). RFLP and subsequent electrophoresis, generated bands of 238 and 127 bp for subjects with SNP rs6311, indicating the presence of G/G genotype. For patients with SNP rs6313, bands of 255 and 191 bp were generated signifying the presence of C/C genotypes ([Figure 2]). Further, for A/G genotype (rs6311), the restriction digestion generated three bands at 365, 238 and 127 bp. For C/T genotype (rs6313), the digestion generated three bands at 446, 255 and 191 bp ([Figure 2]). Representative samples were sequenced using a DNA sequencer (ABI 3730XL Genetic Analyzer, USA), and the sequences were deposited in the GenBank and accession numbers were obtained (JX682563 for CC, KC603613 for AA, KC609429 for GG, KC609430 for AG, KF158420 for TC, JX656724 for TT genotype).
|Figure 2: RFLP products of rs6311 and rs6313. M, 100bp ladder; Lane 1, AG genotype; Lane 2, GG genotype; Lane 3, AA genotype; Lane 4, TC genotype; Lane 5, CC genotype; Lane 6, TT genotype|
Click here to view
Statistical analysis: The allele and the genotype frequencies for the test and the controls were assessed for its possible association with schizophrenia using Chi square test (GraphPad Prism, version 5.0, USA) 95 per cent confidence interval was considered to be the requisite for rejecting null hypothesis.
Genotypic and allelic frequencies of all the HapMap populations (involving Africans, Asians and Europeans) were obtained from HapMap data  for comparison among ethnicities (http://hapmap.ncbi.nlm.nih.gov/ cgi-perl/gbrowse/hapmap28-B36/), and were subjected to Chi square analysis to compare with the Tamil population (present study). Linkage disequilibrium patterns were analyzed and viewed with Haploview (Version 4.2)  . (www.broadinstitute.org/scientific-community/science/programs/medical-and-populationgenetics/ haploview/haploview).
In silico analysis: FAST SNP analysis was done to make predictions on the risk factor. The output score (obtainable from the database) is based on the risk levels in a scale of 0 - 5 (such as, no risk, very low, low, medium, high, and very high risks, respectively)  . Functional SNP database  was used to analyze both the SNPs for its deleterious effects, if any, with respect to the functional categories such as protein coding, splicing regulation, transcriptional regulation, and post-translation, assessed through the respective algorithms.
| Results|| |
Our analysis revealed that 78 of 262 (30%) patients possessed G/G genotypes, 140 patients had A/G genotypes (53%), and 44 patients (17%) had A/A genotypes ([Table 3]). o0 f the 272 controls, 60 (22%) had G /G genotype, 148 and 64 individuals had A/G (54%) and A/A (24%) genotypes, respectively. o0 f the 266 patients, 85 had C/C genotypes (32%), while 141 patients belonged to C/T genotypes (53%); 40 patients were of T/T genotype (15%).
|Table 3: Genotype and allele frequencies of rs6311 and rs6313 in Tamil population |
Click here to view
Application of the Hardy-Weinberg equilibrium and subsequent Chi-square analysis revealed the existence of perceptible differences between the genotype distributions and allelic frequencies of the patients and the controls. Frequencies of the G/G (30%) and the C/C genotypes (32%) (rs6311 and rs6313, respectively), were found to be significantly higher in patients than in the controls; the G allele (56%) and the C allele (58%) frequencies were also found to be higher in patients ([Table 3]). Chi-square analysis on the rs6311 and rs6313 revealed significantly higher frequencies of A/A (rs6311) (24%) and T/T (rs6313) (22%) genotypes in controls compared with cases ([Table 3]). The power analysis revealed a power of 91 per cent for both the SNPs rs6311 and rs6313, with a relative risk of 1.5 (α = 0.05), under a multiplicative model.
In silico analysis: The linkage disequilibrium analyses with Haploview (Version 4.2) generated an LOD value (Logarithm of Odds) as high as 82.37, evincing the occurrence of LD in our study population.
Comparison of genotype frequencies of both SNPs of the Tamil speaking population with those of other populations, obtained from HapMap Database revealed the existence of considerable discrepancy in terms of the relative distribution of the genotypes. For example, GG/AA ratios of rs6311 and CC/TT ratios of rs6313 appeared to differ between the Tamil population (present study) and other populations such as LWK, ASW, GIH, MEX, YRI and MKK. c0 omparison with some other populations such as CHB, CHD, JPT and TSI however, did not show any significant difference in the genotype distribution ([Figure 3]A, B).
|Figure 3: Genotype frequencies of (A) rs6311 and (B) rs6313 in Tamil speaking population and the HapMap populations. P*<0.05, **<0.001; ASW, African ancestry in Southwest USA; CEU, Utah residents with Northern and Western European ancestry from the CEPH collection; CHB, Han Chinese in Beijing, China; CHD, Chinese in Metropolitan Denver, Colorado; GIH, Gujarati Indians in Houston, Texas; JPT, Japanese in Tokyo, Japan; LWK, Luhya in Webuye, Kenya; MEX, Mexican ancestry in Los Angeles, California; MKK, Maasai in Kinyawa, Kenya; TSI, Tuscan in Italy; YRI, Yoruban in Ibadan, Nigeria; TAM, Tamil speaking population, India.|
Click here to view
F-SNP provided the functional information on the SNPs in coding region with respect to splicing and transcriptional regulations  . F-SNP analysis of our data revealed that rs6311 had a significant role in the transcriptional regulation with the FS score of 0.5 and rs6313 has an evident role in the splicing regulation with the FS score of 0.33.
Analysis of the sequence data generated from the amplicons of rs6311 and rs6313 revealed the synonymous nature of these SNPs. Fast-SNP database was used as a tool to foresee the role, if any, of these synonymous SNPs in the regulatory and/or the coding regions of the candidate gene. Fast-SNP predicted a risk ranking of low to medium i.e. 1-3 (in 0-5 scale) for rs6311. Presumably, this situation might lead to alteration in the promoter/regulatory region; those for rs6313, however, was found to have a very low - very low risk ranking i.e. 1-1 (in 0-5 scale).
| Discussion|| |
Our study showed a significant association of the homozygocities of G/G (rs6311) and C/C (rs6313) genotypes with schizophrenia patients and significantly higher frequencies of A/A (rs6311) and T/T (rs6313) in the control population compared with cases.
Although the Haploview analysis revealed the existence of LD between rs6311 and rs6313 in the present study population, data on other populations showed varying results. While the Jewish, Asturian (Northern Spain) and Caucasian populations, for instance, demonstrate the existence of LD between these polymorphisms  , LD does not seem to occur in the ASW population (African ancestry in South West USA), wherein these candidate SNPs are known to be placed in different blocks (Block 9 and 10 respectively, Haploview Version 4.2)  .
The present findings on the schizophrenia patients with G/G and C/C genotypes at rs6311 and rs6313, respectively, are in consonance with those of several other populations. In a meta-analysis performed on the Chinese population, an association was found to exist between rs6311 and schizophrenia in most of the sub-populations assessed  . In another study involving the European population, a significant association between rs6313 and schizophrenia was reported  ; the C allele was associated with diagnosis of schizophrenia in this population. Another meta-analysis revealed the occurrence of strong correlation between rs6313 and schizophrenia in Caucasians  . The allele frequency analyses in the Han population in China, have revealed the predominance of T allele over C allele (rs6313) in schizophrenia patients  . In another study, involving 93 patients from Korea, no significant difference was found between the genotype frequencies (rs6313) and schizophrenia  . In a study on the Turkish population involving 102 patients and 107 controls, no correlation was found between schizophrenia and the genotype frequencies of either of the SNPs (rs6311 and rs6313)  . A positive correlation existed between the genotype frequencies at rs6313 and schizophrenia in a cohort of Japanese population  . However, no association (between rs6313 and schizophrenia) was found in a similar study, but using a larger population from Japan (Kyushu Islands)  .
The SNPs rs6311, rs6313 are the results of synonymous mutations that may lead to functional inadequacies at mRNA levels and disorder in transcriptional activities. In humans, synonymous mutations are suggested to result in deficient protein synthesis caused by functionally defective RNA  . The allelic variant rs6311 is suggested to result in inadequate promoter activity, which in turn may affect the transcription of the gene  . The two variants (rs6311 and rs6313) are in complete LD, the C allele of rs6313 is shown to be associated with lower mRNA and lower protein expression than the T variant  , which may possibly explain the prevailing occurrence of C allele (and CC genotype) in schizophrenia patients in the present study.
Fast SNP analysis suggested that rs6311 was functionally significant with the risk ranking from low - medium (1-3), leading to altering the promoter/regulatory region. Both the SNPs with FS score of 0.5 and 0.33, respectively for rs6311 and rs6313 have been suggested to play a significant role in the transcriptional regulation , .
To conclude, the present study showed an association between the SNPs occurring in 5HT2A and schizophrenia in Tamil speaking population in India. That these two SNPs exhibited considerable degrees of LD and had their origin in the same receptor (5HT2A). This association (between the SNPs and the disease) poses a question if these point mutations (at 5HT2A) could be considered as genetic marker for the disease. Due to the proposed polygenic mode of inheritance of the disease and the variable distribution of these genotypes in several populations distributed across the different continents, these SNPs cannot be considered as exclusive genetic markers. Our study had two major limitations: (i) The study relied on the data obtained exclusively from the Tamil speaking population residing in and around Vellore at Tamil Nadu. s0 amples from other parts of Tamil Nadu State were not collected. (ii) The population sub-types (based on caste) could not be considered.
However, that these two SNPs could potentially contribute to the genesis of the disease in Tamil speaking population should be addressed through future research involving expression studies.
| Acknowledgment|| |
Financial support from the Department of Biotechnology (DBT), Government of India, and the supports from the VIT Management (Vellore) are gratefully acknowledged. Authors thank the Government Vellore Medical College and the Sri Narayani Hospital and Research Centre, for providing the blood samples, and the Directorate of Medical Education for ethical clearance to facilitate the study.
| References|| |
Sullivan PF. The genetics of schizophrenia. PLoS Med
Kim Y, Zerwas S, Trace SE, Sullivan PF. Schizophrenia genetics: Where next? Schizophr Bull
Bertram L. Genetic research in schizophrenia: new tools and future perspectives. Schizophr Bull
Bleich A, Brown SL, van Praag, HM. A serotonergic theory of schizophrenia. In: Brown SL, van Praag HM, editors. The role of serotonin in psychiatric disorders
. New York: Brunner/Mazel; 1991. p. 183-214.
Iqbal N, van Praag HM. The role of serotonin in schizophrenia.EurNeuropsychopharmacol
(Suppl) : 11-23.
Lin MW, Curtis D, Williams N, Arranz M, Nanko S, Collier D, et al
. Suggestive evidence for linkage of schizophrenia to markers on chromosome 13q14.1 - q32. Psychiatr Genet
Norton N, Owen MJ. HTR2A: Association and expression studies in neuropsychiatric genetics. Ann Med
Smith AK, Dimulescu I, Falkenberg VR, Narasimhan S, Heim C, Vernon SD, et al
. Genetic evaluation of the serotonergic system in chronic fatigue syndrome. Psychoneuroendocrinology
Myers RL, Airey DC, Manier DH, Shelton RC, Sanders-Bush E. Polymorphisms in the regulatory region of the human serotonin 5-HT2A receptor gene (HTR2A) influence gene expression. Biol Psychiatry
Abdolmaleky HM, Yaqubi S, Papageorgis P, Lambert AW, Ozturk S, Sivaraman V, et al
. Epigenetic dysregulation of HTR2A in the brain of patients with schizophrenia and bipolar disorder. Schizophr Res
Gaysina D, Zainullina A, Khusnutdinova E. Inherited neuropsychiatric disorders in Russia. In: Kumar D, editor. Genomics and health in the developing world
. Oxford New York: University Press; 2012. p. 1356-67.
Abdolmaleky HM, Faraone SV, Glatt SJ, Tsuang MT. Meta-analysis of association between the T102C polymorphism of the 5HT2A receptor gene and schizophrenia. Schizophr Res
Polesskaya OO, Sokolov BP. Differential expression of the''C'' and ''T'' alleles of the 5-HT2A receptor gene in the temporal cortex of normal individuals and schizophrenics. J Neurosci Res
Khait VD, Huang YY, Zalsman G, Oquendo MA, Brent DA, Harkavy-Friedman JM, et al
. Association of serotonin 5-HT2A receptor binding and the T102C polymorphism in depressed and healthy Caucasian subjects. Neuropsychopharmacology
Yildiz SH, Akilli A, Bagcioglu E, Ozdemir Erdogan M, Coskun KS, Alpaslan AH, et al.
Association of schizophrenia with T102C (rs6313) and 1438 A/G (rs6311) polymorphisms of HTR2A
gene. Acta Neuropsychiatr
American Psychiatric Association (APA). Diagnostic and statistical manual of mental disorders (DSM-IV).
Arlington, Va, USA: APA; 1994.
Sambrook, J, Russell DW. Molecular cloning: a laboratory manual
, 3 rd
ed. New York: Cold Spring Harbor Laboratory; 2001.
Nöthen MM, Rietschel M, Erdmann J, Oberländer H, Möller HJ, Nader D, I. Genetic variation of the 5-HT2A receptor and response to clozapine. Lancet
Higuchi R. Simple and rapid preparation of samples for PCR. In: Erlich HA, editor. PCT technologies
. New York: Oxford University Press, 1992.p. 31-8.
Available from: http://www.hapmap.ncbi.nlm.nih.gov/ cgi-perl/gbrowse/hapmap 28-B36
, accessed on June 20, 2013.
Barrett JC, Fry B, Maller J, Daly MJ. Haploview: analysis and visualization of LD and haplotype maps. Bioinformatics
Yuan HY, Chiou JJ, Tseng WH, Liu CH, Liu CK, Lin YJ, et al
. FASTSNP: an always up-to-date and extendable service for SNP function analysis and prioritization. Nucleic Acids Res
Lee PH, Shatkay H. F-SNP: Computationally predicted functional SNPs for disease association studies. Nucleic Acids Res
Saiz PA, Garcia-Portilla MP, Arango C, Morales B, Alvarez VV, Coto E , et al
. Association study of serotonin 2A receptor (5-HT2A) and serotonin transporter (5-HTT) gene polymorphisms with schizophrenia. Prog Neuropsychopharmacol Biol Psychiatry
Gu L, Long J, Yan Y, Chen Q, Pan R, Xie X, et al
. HTR2A-1438A/G polymorphism influences the risk of schizophrenia but not bipolar disorder or major depressive disorder: a meta-analysis. J Neurosci Res
Williams J, Spurlock G, McGuffin P, Mallet J, Nothen MM, Gill M, et al
. Association between schizophrenia and T102C polymorphism of the 5-hydroxytryptamine type 2a-receptor gene. European Multicentre Association Study of Schizophrenia (EMASS) Group. Lancet
Chen CH, Lee YR, Wei FC, Koong FJ, Hwu HG, Hsiao KJ. Lack of allelic association between 102T/C polymorphism of serotonin receptor type 2A gene and schizophrenia in Chinese. Psychiatr Genet
Lee SK, Wang SK, Chee IS, Lee YH, Shin SC. No association between 102T/C and 452His/Tyr polymorphisms of 5-HT2A receptor gene and schizophrenia in Korean population. J Korean Neuropsychiatr Assoc
Inayama Y, Yoneda H, Ishida T, Nonomura Y, Kono Y, Koh J, et al
. An association between schizophrenia and a serotonin receptor DNA marker (5HTR2). Neuropsychopharmacology
Shinkai T, Ohmori O, Kojima H, Terao T, Suzuki T, Abe K. Negative association between T102C polymorphism of the 5-HT2A receptor gene and schizophrenia in Japan. Hum Hered
[Figure 1], [Figure 2], [Figure 3]
[Table 1], [Table 2], [Table 3]