Introduction
Anxiety disorders (AD) are associated with negative family, affective, social and economic impacts on several levels.1 Most adults with AD report that their anxious symptoms had onset early in life, especially during childhood and adolescence.2,3 According to epidemiological data from multiple sources, the prevalence of AD ranges from 3.1 to 17.5% in children and adolescents,4 while the current global prevalence is 7.3%.5 It is accepted that AD are complex disorders linked to multiple genes6 and several different environmental vulnerabilities.7-9
The hypothalamic-pituitary-adrenal (HPA) axis is the central component of the stress system and is responsible for people's reactivity to stress.10,11 This axis mediates the association between early life stress and development of mood disorders in adulthood.12-16 Studies with animal models have described a number of mechanisms underlying the stress-induced neuronal plasticity leading to fear and anxiety,17 but very few studies in humans have been published.
One animal study18 investigated the consequences of imbalance between mRNA levels of the mineralocorticoid receptor (MR) and the glucocorticoid receptor (GR), which are involved with regulation of the HPA axis,19 suggesting that there is significant interaction between these receptors and control of the HPA-axis under stress. This result further intensified the focus on mechanisms involving the HPA axis in behavior.
Serine protease inhibitors (SERPINA) are a family of proteins that includes the corticosteroid-binding globulin (CBG), which is involved in transport of glucocorticoids. Corticosteroid-binding globulin has the ability to bind 80-90% of cortisol in plasma.20
The NR3C1 gene encodes a GR and has been associated with modification of physiological stress response and mood disorders,21 but not in all studies.22 The FKBP5 gene, a co-chaperone of hsp90, is responsible for regulation of GR sensitivity to cortisol. Overexpression of FKBP5 would be expected to result in an attenuation of the negative feedback inhibition.23 Alleles associated with enhanced expression of FKBP5 following GR activation lead to increased GR resistance and decreased HPA axis efficiency in healthy controls. This result prolongs activation of the stress hormone system following exposure to stress.24 Accordingly, genetic polymorphisms of FKBP5 gene have been associated with recurrence of depressive episodes and with incomplete normalization of stress-elicited cortisol secretion and, in interaction with childhood trauma, predict post-traumatic stress disorder (PTSD), aggressive behavior and suicidal behavior,25-28 and are also linked with increased self-reported anxiety after psychosocial stress.29
Genetic variations in the gene for corticotropin releasing hormone (CRH) are associated with behavioral inhibition in children, which is a temperament that has been related to development of AD in adulthood.30 Corticotropin-releasing hormone receptor 1 subtype (CRHR1) is involved in the anxiogenic action of CRH.31 Animal models have shown that mice lacking a functional CRHR1 exhibited less anxious behavior in comparison with wild mice.32 Therefore, our primary objective is to explore, in a sample of adolescents, the role that polymorphisms in candidate HPA-linked genes play in AD and anxiety symptoms. Our hypothesis is that these genetic variants might be associated with early-onset AD in an adolescent sample.
Methods
Sample selection and psychiatric evaluation
Adolescents were selected from a larger community sample enrolled at schools within the area served by the Hospital de Clínicas de Porto Alegre (HCPA), Brazil. The sampling procedures have been described in detail elsewhere.33 The ethical and methodological aspects of this study were approved and all participants provided informed consent in writing.
All subjects who were rated above the 75th percentile in the Screen for Child and Anxiety related Emotional Disorders (SCARED)34 and a random sample of the other 3 quartiles were invited to attend a clinical interview and were also assessed using the Brazilian version of the Schedule for Affective Disorder and Schizophrenia for School-Age Children - Present and Lifetime (K-SADS-PL),35,36 which is a semi-structured interview based on Diagnostic and Statistical Manual of Mental Disorders, fourth edition (DSM-IV) criteria.
Adolescents who were diagnosed with social anxiety disorder, generalized anxiety disorder, separation anxiety disorder or panic disorder according to the K-SADS-PL were classified as cases. Subjects with no AD other than specific phobias were classified as controls. Anxiety severity symptoms were measured using SCARED results expressed as means and standard deviations.
We excluded from this study people who had: 1) significant organic illness; 2) history of bipolar disorder, pervasive developmental disorder or any psychotic disorder; 3) history of alcohol or drug dependency or abuse; 4) clinical suspicion of intellectual disability and 5) subjects whose genotype was not clearly defined.
Participants
The initial sample comprised 240 adolescents who agreed to participate in the psychiatric diagnostic assessment and provided biological material (saliva) for DNA studies. Only unrelated individuals were included (n = 234). Six subjects were excluded because of intellectual disability (3 controls and 3 cases). The final sample comprised 228 adolescents (67.98% females; mean age = 12.81±2.39): 131 cases (74% females; mean age = 13.1±2.36; 65.9% Caucasians) and 97 controls (59.8% females; mean age = 12.54±2.41; 73.4% Caucasians). Ethnicity was determined by self-report on the basis of skin color. Our entire adolescent sample comprised 152 Caucasians, 30 African-Brazilians and 38 mixed individuals with skin color between white and black. Data on ethnicity were missing for 8 participants. There was no difference in ethnicity distribution between cases and controls.
There were a total of 131 cases, 93 (71%) of whom were diagnosed with generalized anxiety disorder, 55 (42%) with social anxiety disorder, 46 (35.1%) with separation anxiety disorder and 9 (6.9%) with panic disorder. Of these, 75 (57.3%) had one AD only, 41 (31.3%) had two AD and 15 (1.5%) had three or four AD, as has been observed elsewhere.37
Selection of polymorphisms of candidate genes
We selected single nucleotide polymorphisms (SNPs) on the basis of the results of searches performed on the MEDLINE, Cochrane and PubMed databases for studies involving dysregulation of the HPA axis. Candidate genes were chosen on the basis of their roles in psychiatric disorders, according to the literature. We found 17 SNPs of candidate HPA linked-genes: FKBP5 (rs3800373, rs9296158, rs1360780, rs9470080, rs4713916); AVPR1B (rs28632197, rs28529127, rs28575468, rs33933482); NR3C1 (rs6198, rs6189); SERPINA6 (rs746530); NR3C2 (rs2070951, rs5522); CRH (rs6159, rs6999100) and CRHR1 (rs878886). Five of these SNPs (rs6189, rs4713916, rs5522, rs6159, rs69999100) were excluded because their minor allele counts (MAF) were below 5%. Our final selection comprised 8 SNPs of candidate genes: FKBP5 (rs3800373, rs9296158, rs1360780, rs9470080 and rs4713916); NR3C1 (rs6198); CRHR1 (rs878886) and SERPINA6 (rs746530).
Genotyping
DNA was extracted from biological samples of saliva using an OrageneTM Kit. All SNPs were detected using real-time polymerase chain reaction (PCR). All of the conditions of the genotyping process have been described in detail elsewhere.38
Statistical analyses
Analyses of associations and the test of Hardy-Weinberg Equilibrium (HWE) were performed using PLINK software.39 The overall level of significance adopted was α = 0.05. Descriptive data were obtained using SPSS version 18 and expressed as means, standard deviations, and n (%).
We performed single marker association tests and we adjusted the p-values to incorporate correction for multiple hypothesis testing with permutations (n = 10,000 permutations) corrected for both gene-wise (p1) and experiment-wise (p2) significance values.
Results
Our total sample was in HWE. Table 1 lists descriptive statistics for the test of HWE, the case-control association study (K-SADS-PL analysis) and the severity of anxiety symptoms according to SCARED. The tests of allelic associations detected no differences between AD cases and controls in the frequencies of our selection of SNPs.
Table 1 - Allelic association between SNPs of HPA-linked genes and K-SADS-PL and SCARED scores (n = 228) in a case-control study of anxiety disorders
K-SADS-PL (n = 228) | SCARED (n = 228) | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|---|
Gene | SNP | MA | FA | FU | χ2 | OR (95%CI) | p | p HWE | β | r² | T | p permuted |
NR3C1 | rs6198 | C | 0.09 | 0.09 | 0.004 | 0.98 (0.51-1.86) | 0.951 | 0.702 | -0.895 | 0.0008 | -0.433 | 0.666 |
FKBP5 | rs3800373 | C | 0.33 | 0.30 | 0.4386 | 1.15 (0.76-1.72) | 0.508 | 0.17 | 0.125 | 4.765e-005 | 0.103 | 0.918 |
rs9296158 | A | 0.38 | 0.35 | 0.414 | 1.14 (0.77-1.68) | 0.520 | 0.253 | 0.658 | 0.001 | 0.563 | 0.574 | |
rs1360780 | T | 0.34 | 0.31 | 0.5865 | 1.17 (0.78-1.74) | 0.444 | 0.178 | 0.128 | 5.098e-005 | 0.107 | 0.915 | |
rs9470080 | T | 0.38 | 0.35 | 0.3483 | 1.12 (0.76-1.66) | 0.555 | 0.318 | 0.575 | 0.001 | 0.489 | 0.625 | |
rs4713916 | A | 0.23 | 0.25 | 0.1193 | 0.93 (0.60-1.43) | 0.730 | 0.858 | -1.923 | 0.009 | -1.431 | 0.154 | |
SERPINA6 | rs746530 | A | 0.37 | 0.44 | 2.231 | 0.75 (0.51-1.09) | 0.135 | 0.489 | -2.155 | 0.014 | -1.766 | 0.079 |
CRHR1 | rs878886 | G | 0.13 | 0.19 | 2.621 | 0.65 (0.39-1.09) | 0.105 | 0.465 | -0.973 | 0.002 | -0.610 | 0.543 |
95%CI = 95% confidence interval; β = regression coefficient; r² = regression r-squared; χ2 = chi-square test; CRHR1 = corticotropin releasing hormone receptor 1; FKBP5 = FK506 binding protein 5; FA = frequency of risk allele in cases; FU = frequency of risk allele in controls; HPA = hypothalamic-pituitary-adrenal; HWE = Hardy-Weinberg equilibrium; K-SADS-PL = The Schedule for Affective Disorder and Schizophrenia for School-Age Children - Present and Lifetime; MA = minor allele; NR3C1 = nuclear receptor subfamily 3, group C, member 1 (glucocorticoid receptor); NR3C2 = nuclear receptor subfamily 3, group C, member 2 (mineralocorticoid receptor); OR = odds ratio; p = Wald-test asymptotic p-value; SCARED = Screen for Child and Anxiety related Emotional Disorders; SERPINA6 = serpin peptidase inhibitor, clade A (alpha-1 antiproteinase, antitrypsin), member 6; SNP = single nucleotide polymorphism; T = Wald test (based on t distribution).
Discussion
Contrary to our hypothesis, we did not find any associations between a selection of genetic variants (FKBP5: rs3800373, rs9296158, rs1360780, rs94700870, rs4713916; NR3C1: rs6198; CRHR1: rs878886 and SERPINA6: rs746530) and AD in our adolescent sample.
The high proportion of psychiatric comorbidities in our sample may have contributed to the absence of evidence of associations between HPA-linked genes and AD. Since our sample is comprised of people from the community, it is possible that adolescents who had externalizing disorder comorbidities could have attended the psychiatric evaluation more frequently.
We observed that there is a paucity of studies investigating genetic variants of the FKBP5 gene in connection with AD,40 although some studies have investigated this gene with relation to other psychiatric disorders such as bipolar disorder and depression and environmental interactions with it in relation to risk of PTSD.28,41,42 Therefore, in an attempt to infer some hypotheses concerning AD, our discussions will be partially based on these other psychiatric disorders.
Minelli et al.40 showed that homozygosis for the T allele of rs1360780 is significantly more frequent in patients affected by major depression with co-morbid AD than in those without co-morbid AD. A study with young adults indicated faster cortisol recovery after acute stress in individuals with the CC genotype who had been exposed to maltreatment during childhood, whereas no association was observed among T allele carriers.43 The T allele and TT genotype of rs1360780 were overrepresented in individuals with depression28,44 and also had a strong trend for association with bipolar disorder in a large family-based sample.41 Additionally, the TT genotype of rs1360780 exhibited incongruent results in relation to faster response to antidepressant treatment with studies observing both associations28 and absence of associations45 to this polymorphism. The FKBP5 polymorphisms (rs4713916, rs3800373 and rs1360780) were also investigated in a sample of affective psychosis cases (depression and bipolar disorder) and controls, but without any associations.46 Moreover, all five of the FKBP5 gene SNPs that we studied were also evaluated in a recent study that suggested an association in depression patients but not with bipolar disorders.47 These results illustrate the challenge of finding associations involving candidate genes and complex disorders.
The NR3C1 gene (GR) has previously been linked to depression and bipolar disorder.21,48 The interaction between rs6198 of the NR3C1 gene and rs9470080 of the FKBP5 gene may modulate the risk for depression by contributing to inhibition of GR signaling,47 although no association between this polymorphism and melancholic features in the course of bipolar disorder was found in another study.22 We also did not find associations between these polymorphism and AD.
Our study also involved the rs878886 polymorphism of the CRHR1 gene. In two German samples49 rs878886 was highly associated with panic disorders, but was not associated with depression in the course of bipolar disorder.22 Among young adults, G-allele carriers showed no acquisition of fear conditioned responses to the threat cue in the uninstructed phase, whereas fear acquisition was present in people homozygous for the C allele.50 Fear acquisition deficit is relevant in the pathogenesis of human AD.51
Finally, we evaluated the rs746530 polymorphism of the SERPINA6 gene in AD. We found only one study describing associations between this specific SNP and multiple somatic symptoms (including symptoms of anxiety).52
Our study has two main limitations. First, our sample size was small and this may have limited our power to detect associations with modest to moderate effect sizes. Second, we did not test HPA axis function by assaying cortisol levels under stress stimuli.
Our study provides additional data about some polymorphisms of HPA-linked genes in AD. Although these HPA-specific SNPs are biologically plausible; there is a paucity of literature concerning their relationships with development of psychiatric disorders. Another important feature of this study was the use of a sample of adolescents. Adolescence is a critical period for development of anxiety and affective disorders in response to stress.53 We believe that these data should encourage future studies investigating the stress axis and functional polymorphisms, taking environmental interactions and epigenetic influences into consideration.
Conclusion
We used a strict methodology to select functional polymorphisms of HPA-linked genes and did not find any associations between these genetic variants and AD in adolescents. We suggest that further studies involving these genes should be designed and conducted using larger adolescent samples, in order to improve our understanding of possible HPA axis mechanisms in AD during this important developmental period.