Hostname: page-component-cd9895bd7-gbm5v Total loading time: 0 Render date: 2024-12-13T11:13:22.327Z Has data issue: false hasContentIssue false

Genome-wide association studies in non-anxiety individuals identified novel risk loci for depression

Published online by Cambridge University Press:  22 June 2022

Bolun Cheng
Affiliation:
Key Laboratory of Trace Elements and Endemic Diseases, Collaborative Innovation Center of Endemic Disease and Health Promotion for Silk Road Region, School of Public Health, Health Science Center, Xi’an Jiaotong University, Xi’an 710061, China
Xin Qi
Affiliation:
Precision Medicine Center, The First Affiliated Hospital of Xi’an Jiaotong University, Xi’an, China
Peilin Meng
Affiliation:
Key Laboratory of Trace Elements and Endemic Diseases, Collaborative Innovation Center of Endemic Disease and Health Promotion for Silk Road Region, School of Public Health, Health Science Center, Xi’an Jiaotong University, Xi’an 710061, China
Shiqiang Cheng
Affiliation:
Key Laboratory of Trace Elements and Endemic Diseases, Collaborative Innovation Center of Endemic Disease and Health Promotion for Silk Road Region, School of Public Health, Health Science Center, Xi’an Jiaotong University, Xi’an 710061, China
Xuena Yang
Affiliation:
Key Laboratory of Trace Elements and Endemic Diseases, Collaborative Innovation Center of Endemic Disease and Health Promotion for Silk Road Region, School of Public Health, Health Science Center, Xi’an Jiaotong University, Xi’an 710061, China
Li Liu
Affiliation:
Key Laboratory of Trace Elements and Endemic Diseases, Collaborative Innovation Center of Endemic Disease and Health Promotion for Silk Road Region, School of Public Health, Health Science Center, Xi’an Jiaotong University, Xi’an 710061, China
Yao Yao
Affiliation:
Key Laboratory of Trace Elements and Endemic Diseases, Collaborative Innovation Center of Endemic Disease and Health Promotion for Silk Road Region, School of Public Health, Health Science Center, Xi’an Jiaotong University, Xi’an 710061, China
Yumeng Jia
Affiliation:
Key Laboratory of Trace Elements and Endemic Diseases, Collaborative Innovation Center of Endemic Disease and Health Promotion for Silk Road Region, School of Public Health, Health Science Center, Xi’an Jiaotong University, Xi’an 710061, China
Yan Wen
Affiliation:
Key Laboratory of Trace Elements and Endemic Diseases, Collaborative Innovation Center of Endemic Disease and Health Promotion for Silk Road Region, School of Public Health, Health Science Center, Xi’an Jiaotong University, Xi’an 710061, China
Feng Zhang*
Affiliation:
Key Laboratory of Trace Elements and Endemic Diseases, Collaborative Innovation Center of Endemic Disease and Health Promotion for Silk Road Region, School of Public Health, Health Science Center, Xi’an Jiaotong University, Xi’an 710061, China
*
*Author for correspondence: Feng Zhang, E-mail: fzhxjtu@mail.xjtu.edu.cn

Abstract

Background

Depression is a debilitating mental disorder that often coexists with anxiety. The genetic mechanisms of depression and anxiety have considerable overlap, and studying depression in non-anxiety samples could help to discover novel gene. We assess the genetic variation of depression in non-anxiety samples, using genome-wide association studies (GWAS) and linkage disequilibrium score regression (LDSC).

Methods

The GWAS of depression score and self-reported depression were conducted using the UK Biobank samples, comprising 99,178 non-anxiety participants with anxiety score <5 and 86,503 non-anxiety participants without self-reported anxiety, respectively. Replication analysis was then performed using two large-scale GWAS summary data of depression from Psychiatric Genomics Consortium (PGC). LDSC was finally used to evaluate genetic correlations with 855 health-related traits based on the primary GWAS.

Results

Two genome-wide significant loci for non-anxiety depression were identified: rs139702470 (p = 1.54 × 10−8, OR = 0.29) locate in PIEZO2, and rs6046722 (p = 2.52 × 10−8, OR = 1.09) locate in CFAP61. These associated genes were replicated in two GWAS of depression from PGC, such as rs1040582 (preplication GWAS1 = 0.02, preplication GWAS2 = 2.71 × 10−3) in CFAP61, and rs11661122 (preplication GWAS1 = 8.16 × 10−3, preplication GWAS2 = 8.08 × 10−3) in PIEZO2. LDSC identified 19 traits genetically associated with non-anxiety depression (p < 0.001), such as marital separation/divorce (rg = 0.45, SE = 0.15).

Conclusions

Our findings provide novel clues for understanding of the complex genetic architecture of depression.

Type
Research Article
Creative Commons
Creative Common License - CCCreative Common License - BY
This is an Open Access article, distributed under the terms of the Creative Commons Attribution licence (http://creativecommons.org/licenses/by/4.0), which permits unrestricted re-use, distribution and reproduction, provided the original article is properly cited.
Copyright
© The Author(s), 2022. Published by Cambridge University Press on behalf of the European Psychiatric Association

Introduction

Depression is one of the complex and common neuropsychiatric disorders affecting approximately 4.4% of the population worldwide, which often presents with low self-esteem, low mood, anhedonia, feeling of worthlessness, fatigue, sense of rejection and guilt, suicidal thoughts, among others [Reference Pandi-Perumal, Monti, Burman, Karthikeyan, BaHammam and Spence1, Reference Villas Boas, Boerngen de Lacerda, Paes, Gubert, Almeida and Rescia2]. The lifetime prevalence of major depression around the world is between 1.0 and 16.9% [Reference Cheng, Chu, Yang, Wen, Jia and Liang3]. Scientists have identified a familial tendency for depression through two decade of family studies [Reference Sullivan, Neale and Kendler4]. However, a study sponsored by the Psychiatric Genomics Consortium (PGC) in 2013 did not find any associated loci of genetic variation in depression [Reference Ripke, Wray, Lewis, Hamilton, Weissman and Breen5], and a subsequent meta-analysis conducted by Hek et al. was similarly inconclusive [Reference Hek, Demirkan, Lahti, Terracciano, Teumer and Cornelis6]. Like other common diseases or traits, depression is thus presumed to be controlled by many genes with minor effects.

Due to the wide range of depression phenotypes, previous studies could not achieve the statistical power required to test these minor effect variations. Therefore, PGC again combined the cohort data from PGC, UK BioBank, FinnGen and 23andMe to find some genetic associations for depression [Reference Wray, Ripke, Mattheisen, Trzaskowski, Byrne and Abdellaoui7]. In the latest study, Levey et al. performed genome-wide association study (GWAS) and meta-analysis using depression data from more than 1.2 million participants in multiple population cohorts, and identified 17 pathogenic genes by fine-mapping GWAS signals combined with transcriptome association analysis [Reference Levey, Stein, Wendt, Pathak, Zhou and Aslan8]. Nevertheless, the existing studies seem to be not account for the influence of comorbidities to the pathogenesis of depression, and lack of further analysis and support of accurate phenotypic cohort data. Despite the previous family studies have shown that about 40% of depression susceptibility is due to genetic effects, little is known about the specific genetic variants involved in depression [Reference Sullivan, Neale and Kendler4].

There is evidence that some mental disorders are typically heritable and share common genetic components [Reference Anttila, Bulik-Sullivan, Finucane, Walters, Bras and Duncan9, Reference Bartlett, Singh and Hunter10]. Anxiety and depression have been demonstrated to be bidirectional risk factors for one another [Reference Jacobson and Newman11]. Previous studies also found considerable overlapped genetic factors between anxiety and depression [Reference Waszczuk, Zavos, Gregory and Eley12, Reference Craske, Rauch, Ursano, Prenoveau, Pine and Zinbarg13]. A recent genetic correlation analysis found approximately 80% genetic correlation between depression and generalized anxiety [Reference Morneau-Vaillancourt, Coleman, Purves, Cheesman, Rayner and Breen14]. However, there is a considerable pathogenesis difference between depression and anxiety [Reference van Dalen, Dierckx, Pasmans, Aendekerk, Mathijssen and Koudstaal15]. For example, the release of some peptides or hormones on the hypothalamic–pituitary–adrenal (HPA) axis is regulated differently in anxiety and depression [Reference Boyer16]. Besides, there are gender differences between the two diseases. Gao et al. found that anxiety was one of the most serious problems among female college students, while depression was relatively more serious among male [Reference Gao, Ping and Liu17]. Therefore, it is reasonable to explore depression specific genetic risk factors without the influence of anxiety. Our aim is to explore depression in the non-anxiety samples, and to eliminate common overlapping risk genes for anxiety and depression.

In this study, four genome-wide association studies were conducted using two depression phenotypes in two non-anxiety cohorts from the UK Biobank, respectively. Then, two GWAS summary statistics from PGC were used to verify the genes corresponding to the candidate loci in our GWAS. Finally, linkage disequilibrium score regression (LDSC) was performed to analyze the genetic correlation between non-anxiety depression and 855 health-related traits.

Methods

Non-anxiety depression samples in the UK Biobank cohort

The phenotypic and genotypic data of this study were derived from UK Biobank health resource (Application 46478), which had recruited 502,656 participants aged between 40 and 69 years [Reference Bycroft, Freeman, Petkova, Band, Elliott and Sharp18]. The present study accessed health-related records of each participant, including age, sex, tobacco and alcohol consumption, and Townsend deprivation index (TDI) from screenshot question or verbal interview within Assessment Center. Anxiety (UK Biobank data fields: 20421 and 20420) and depression (UK Biobank data fields: 20002, 20126 and 20544) were defined based on the general anxiety disorder (GAD-7) and Patient Health Questionnaire (PHQ-9), respectively [Reference Davis, Cullen, Adams, Brailean, Breen and Coleman19]. We used self-reported and mental illness scores to define the phenotype of depression and anxiety, respectively. Ethical approval of UK Biobank was granted by the National Health Service National Research Ethics Service (reference 11/NW/0382). Anxiety and depression score were mean-centered and normalized to one standard deviation (SD) before further analysis. The detailed definitions of mental phenotypes are shown in Supplementary File S1. In this study, individuals with an anxiety score <5 or non-self-reported anxiety were defined as non-anxiety individuals. The samples with anxiety score <5 included 59,334 depression cases (mean ± SD age, 57.22 ± 7.46) and 19,805 controls (mean ± SD age, 55.77 ± 7.40), while the samples with non-self-reported anxiety included 56,603 depression cases (mean ± SD age, 57.21 ± 7.48) and 13,123controls (mean ± SD age, 55.90 ± 7.45). The descriptive characteristics of participants with anxiety score <5 and non-self-reported anxiety are presented in Supplementary Files S2 and S3, respectively.

UK Biobank genotyping, imputation, and quality control

In the UK Biobank, 488,377 participants have genome-wide genotype data. Genome-wide genotyping was conducted using either the Affymetrix UK BiLEVE Axiom or Affymetrix UK Biobank Axiom array. Details of the array design, genotyping, and quality control procedures have been descripted in the published study [Reference Bycroft, Freeman, Petkova, Band, Elliott and Sharp18]. Imputation was conducted by IMPUTE2 against the reference panel of the Haplotype Reference Consortium, 1,000 Genomes and UK10K projects [Reference Bycroft, Freeman, Petkova, Band, Elliott and Sharp18]. Detailed information about these data have been described elsewhere [Reference Canela-Xandri, Rawlik and Tenesa20]. The SNPs with high linkage disequilibrium (r 2 > 0.5) were removed. The participants were restricted to only “White British” according to self-reported ethnicity. The participants who reported inconsistencies between self-reported gender or genetic gender, and were genotyped but not imputed were finally excluded in this study.

Genome-wide association studies of non-anxiety depression

PLINK 2.0 was used to conduct the GWAS of two depression traits in two non-anxiety cohorts, respectively [Reference Purcell, Neale, Todd-Brown, Thomas, Ferreira and Bender21]. For quality control, we removed the SNPs with call rates <90%, Hardy–Weinberg equilibrium (HWE) <0.001, or minor allele frequencies (MAF) <0.01. The kinship coefficients were estimated by KING software (http://people.virginia.edu/~wc9c/KING/) to remove the genetically related subjects [Reference Bycroft, Freeman, Petkova, Band, Elliott and Sharp18]. The GWAS of depression score and self-reported depression were conducted in two non-anxiety cohorts using linear regression and logistic regression assuming an additive model for allelic effects, respectively. The age, sex, TDI, alcohol use frequency/week, smoking frequency/day, and top three principle components of population structure (calculated by UK Biobank) were used as covariates. The SNPs with p < 5.0 × 10−8 were considered to be genome-wide susceptibility significance.

Replication of primary GWAS results

Two large-scale GWAS summary data of depression from the PGC were recruited to verify the accuracy of non-anxiety depression GWAS [Reference Wray, Ripke, Mattheisen, Trzaskowski, Byrne and Abdellaoui7, Reference Howard, Adams, Clarke, Hafferty, Gibson and Shirali22]. Briefly, in replication GWAS 1, a genome-wide association meta-analysis was conducted based on 135,458 depression cases and 344,901 controls using logistic regression [Reference Wray, Ripke, Mattheisen, Trzaskowski, Byrne and Abdellaoui7]. In replication GWAS 2, 807,553 discovery individuals (246,363 depression cases and 561,190 controls) were analyzed using logistic regression [Reference Howard, Adams, Clarke, Hafferty, Gibson and Shirali22]. Detailed descriptions of genotyping, quality control and statistical analysis of these two data sets are available in the published studies [Reference Wray, Ripke, Mattheisen, Trzaskowski, Byrne and Abdellaoui7, Reference Howard, Adams, Clarke, Hafferty, Gibson and Shirali22].

Generating genetic correlations between non-anxiety depression and human traits

We used LDSC to estimate rg for non-anxiety depression with a range of other diseases and health-related traits [Reference Bulik-Sullivan, Finucane, Anttila, Gusev, Day and Loh23]. The purpose of these comparisons was to assess the extent of shared common genetic variants in order to suggest hypotheses about the underlying genetic basis of non-anxiety depression. The overlap of the subjects themselves does not bias rg. These rg are mostly based on independent subject studies and are expected to be unbiased by confounding of genetic and nongenetic effects. rg remains unbiased when GWAS include overlapping samples, but the intercept of LDSC regression is an estimated correlation between association statistics attributable to overlapping samples. In this study, we used the cross-trait LDSC method through the LD Hub v1.9.3 to identify the genetic correlations between non-anxiety depression phenotype and 855 human diseases/traits, including physical and mental diseases, anthropometric markers, living habits and other health-related traits [Reference Zheng, Erzurumluoglu, Elsworth, Kemp, Howe and Haycock24].

Results

Primary analysis of non-anxiety depression

An overview for the GWAS of non-anxiety depression is shown in Figure 1. In the analysis of self-reported depression in anxiety score <5 samples, one SNP reached genome-wide significance: rs6046722, p = 2.52 × 10−8, OR = 1.09 (Supplementary File S4). This SNP is located in an exon of the CFAP61 (cilia and flagella associated protein 61) gene. In the GWAS of depression score in anxiety score <5 samples, rs139702470 reached genome-wide significance: p = 1.54 × 10−8, OR = 0.29 (Supplementary File S4). This is an intronic variant in the PIEZO2 (piezo type mechanosensitive ion channel component 2) gene. In the analysis of depression score in non-self-reported anxiety samples, one SNP reached genome-wide significance: rs139702470, p = 3.66 × 10−8, OR = 0.29 (Supplementary File S4). This SNP is located in an exon of the PIEZO2 gene. Figure 2 shows the LocusZoom plot with data coming directly from our GWAS summary data by querying the corresponding region on chromosomes between 500 and 500 kb region, respectively.

Figure 1. Manhattan plot for the GWAS of depression without anxiety in the UK Biobank cohorts. (A) Linear regression of depression score in anxiety score <5 samples. (B) Linear regression of depression score in non-self-reported anxiety samples. (C) Logistic regression of self-reported depression in anxiety score <5 samples. (D) Logistic regression of self-reported depression in non-self-reported anxiety samples. The red line indicates the p-value threshold for genome-wide significance (p < 5 × 10−8).

Figure 2. LocusZoom plots of depression without anxiety genome-wide significance loci. Association results for SNPs as a function of genomic distance for PIEZO2 and CFAP61 (C20orf26). The top line in each subfigure shows genomic coverage at the locus, with each vertical tick representing the imputed SNPs. Purple diamond indicate SNP at the locus with the strongest association evidence. Each point represents a SNP. Bottom panel shows genes at each locus as annotated in the UCSC Genome Browser Annotation Database. (A) display PIEZO2 in chr18 for GWAS summary of depression score in anxiety score <5 samples. (B) display PIEZO2 in chr18 for GWAS summary of depression score in non-self-reported anxiety samples. (C) display CFAP61 in chr20 for GWAS summary of self-reported depression in anxiety score <5 samples. (D) display CFAP61 in chr20 for GWAS summary of self-reported depression in non-self-reported anxiety samples.

Replication of primary analysis results

The genes corresponding to genome-wide significant loci for non-anxiety depression in the primary analysis were tested in two depression cohorts from the PGC. Both CFAP61 and PIEZO2 were associated with depression in the two replication studies. We observed 9 and 14 candidate SNPs corresponding to CFAP61 and PIEZO2 in two replication datasets (p < 0.05), respectively. For example, rs1040582 (p replication GWAS1 = 0.02, p replication GWAS2 = 2.71 × 10−3) and rs13038510 (p replication GWAS1 = 0.02, p replication GWAS2 = 9.83 × 10−3) were replicated in CFAP61 region (Supplementary File S5), while rs11661122 (p replication GWAS1 = 8.16 × 10−3, p replication GWAS2 = 8.08 × 10−3) and rs11664237 (p replication GWAS1 = 8.58 × 10−4, p replication GWAS2 = 3.46 × 10−3) were replicated in PIEZO2 region (Supplementary File S6).

Linkage disequilibrium score regression

Non-anxiety depression is comorbid with a wide range of other diseases and disorders. To assess the shared genetic architecture between non-anxiety depression and many other traits, genetic correlations (rg) were calculated between our GWAS summary statistics and 855 behavioral and disease traits available via LD Hub. We first focused on the genetic correlations between non-anxiety depression and 12 common mental disorders (Figure 3). Suggestive significant genetic associations were detected for non-anxiety depression with major depressive disorder (rg = 0.27, SE = 0.13) and schizophrenia (rg = 0.15, SE = 0.06). Besides, several weak genetic correlations were observed between non-anxiety depression and many depression-related phenotypes, such as depression (rg = −0.001, SE = 0.09) and number of depression episodes (rg = −0.02, SE = 0.20). Notably, there was no significant genetic correlations between non-anxiety depression and other psychiatric phenotypes such as anxiety (rg = 0.14, SE = 0.18).

Figure 3. Genetic correlations and mental disorders related traits using LD score regression implemented in LD Hub software. The negative rg indicates that an earlier or lower value of a continuous trait was associated with depression without anxiety. The positive rg indicates that a later or higher value of a continuous trait was associated with depression without anxiety.

Of the other behavioral and disease traits, 19 phenotypes were significantly correlated (p < 0.001) with non-anxiety depression (Figure 4), such as overall health rating (rg = 0.21, SE = 0.06), hayfever/allergic rhinitis (rg = −0.30, SE = 0.09), high cholesterol (rg = 0.30, SE = 0.09), current tobacco smoking (rg = 0.21, SE = 0.07), disability or infirmity (rg = 0.18, SE = 0.06), impedance of arm (rg = −0.13, SE = 0.04), and overweight (rg = 0.18, SE = 0.07). Additionally, a novel genetic correlation was observed between non-anxiety depression and marital separation/divorce (rg = 0.45, SE = 0.15) (Figure 4). Detailed results for genetic correlations between non-anxiety depression and other behavioral and disease related traits are summarized in Supplementary File S7.

Figure 4. Significant genetic correlations and other behavioral and disease related traits using LD score regression implemented in LD Hub software. The negative rg indicates that an earlier or lower value of a continuous trait was associated with depression without anxiety. The positive rg indicates that a later or higher value of a continuous trait was associated with depression without anxiety.

Discussion

In this study, we selected approximately 100,000 non-anxiety individuals from the UK Biobank, and conducted GWASs for depression scores and self-reported depression using the two non-anxiety cohorts, respectively. Our GWASs identified two independent SNPs associated with non-anxiety depression. Both of the genes corresponding to the two SNPs were verified by significant association signals across the two replication studies.

The genetic correlations between non-anxiety depression and the general depression phenotype were relatively weak in our study. A recent study using the UK Biobank cohort has found significant overlapping variants between depression and other mental disorders [Reference Brailean, Curtis, Davis, Dregan and Hotopf25]. In contrast, our LDSC results did not find any strong genetic correlation between non-anxiety depression with neuroticism, anorexia nervosa, and anxiety/panic attacks. This supports our hypothesis that the analysis of non-anxiety depression may exclude the effects of comorbidities such as anxiety disorders. Interestingly, a significant genetic correlation was observed between schizophrenia and non-anxiety depression, indicating that schizophrenia may have a strong similarity to the genetic pathogenesis of depression.

Previous studies found that generalized anxiety and depression had a substantial genetic overlap, approximately 80–96% [Reference Morneau-Vaillancourt, Coleman, Purves, Cheesman, Rayner and Breen14, Reference Gustavson, Franz, Panizzon, Reynolds, Xian and Jacobson26]. Our study aimed to explore the depression specific genetic factors by excluding all anxiety individuals from the depression samples. In our LDSC analysis, the GWAS summary data of non-anxiety depression was generated from samples excluded anxiety individuals. On the contrary, the GWAS summary data of depression and number of depression episodes in LD Hub were not excluded anxiety individuals [Reference Zheng, Erzurumluoglu, Elsworth, Kemp, Howe and Haycock24]. In this case, the rg correlation values between non-anxiety depression and depression traits (such as depression [−0.001] and number of depression episodes [−0.02]) were lower than that with anxiety [0.14], illustrating the potentially strong correlation between depression and anxiety. Thus, after the exclusion of the effect of anxiety individuals, we hypothesized that the genetic correlation was inevitably reduced between non-anxiety depression and depressive traits (including anxiety samples).

Examining significant genes that overlap between the current GWAS of non-anxiety depression and the depression studies revealed putative associations with PIEZO2 and CFAP61 [Reference Ripke, Wray, Lewis, Hamilton, Weissman and Breen5, Reference Wray, Ripke, Mattheisen, Trzaskowski, Byrne and Abdellaoui7, Reference Howard, Adams, Clarke, Hafferty, Gibson and Shirali22]. PIEZO2, as the main biological force conduction medium, affects the release of neurotransmitter serotonin (5-HT) [Reference Alcaino, Knutson, Treichel, Yildiz, Strege and Linden27]. The 5-HT plays a central role in brain development, mood regulation, stress response, and the risk of psychiatric disorders, and changes of 5-HT have important implications for behavior and mental health [Reference Brummelte, Mc Glanaghy, Bonnin and Oberlander28]. PIEZO2 is also associated with sensory nerve distribution in the central nervous system of the brainstem and cationic channels activated by brain metastatic cells [Reference Pardo-Pastor, Rubio-Moscardo, Vogel-González, Serra, Afthinos and Mrkonjic29, Reference Umans and Liberles30]. PIEZO2 usually expressed in the cortical and hippocampal pyramidal neurons of the brain and in cerebellar Purkinje cells [Reference Wang and Hamill31]. A whole-exome sequencing analysis in two patients who had unique neuromuscular and skeletal symptoms showed that PIEZO2 was a determinant of mechanosensation in humans [Reference Chesler, Szczot, Bharucha-Goebel, Ceko, Donkervoort and Laubacher32]. Lanier et al. conducted a study of brain injury caused by repeated blast exposure, and indicated that blast may cause PIEZO2 change in sensitivity to mechanical stimuli in the brain and may contribute to cellular injury [Reference Heyburn, Abutarboush, Goodrich, Urioste, Batuure and Statz33].

In our study, CFAP61 was found to be significantly associated with non-anxiety depression and had been identified to play a vital role in primary cilia affecting cerebral cortical development and dysfunction [Reference Liu, Trupiano, Simon, Guo and Anton34]. It was highly expressed in preamygdala, striatum, and hippocampal structures [Reference Rouillard, Gundersen, Fernandez, Wang, Monteiro and McDermott35]. However, the role of CFAP61 in depression has not been well studied, and the role of CFAP61 in the brain or nervous system is also limited. CFAP61 gene identified by our GWAS is a remarkable finding for non-anxiety depression with respect to known biology and points to the potential value of other novel findings from this kind of research. Improving our understanding of the SNPs in non-anxiety depression may help identify the biological differences between depression and anxiety. These results also provide potential evidence for future phenotypic stratification research.

To date, neither PIEZO2 nor CFAP61 has been reported to cause a phenotype related to depression. It is important to note that the closest gene to the GWAS top signal is usually not the causal gene [Reference Zhu, Zhang, Hu, Bakshi, Robinson and Powell36]. Besides, studies have shown that many genetic variants can affect phenotype through distal regulation such as long range enhancer–promoter interactions [Reference Won, de la Torre-Ubieta, Stein, Parikshak, Huang and Opland37], and looping chromatin interactions [Reference Whalen, Truty and Pollard38]. Therefore, the SNPs identified in this study may affect the expression of other depression-related genes nearby. Several genes that are nearly adjacent to PIEZO2 and CFAP61 have been implicated in depression. For example, Dóra et al. carried out whole-exome ultra-high throughput sequencing in brain samples between depression and control subjects, and identified GNAL (near PIEZO2) as one of the genomic region-dependent accumulation of rare variants in depression [Reference Tombácz, Maróti, Kalmár, Csabai, Balázs and Takahashi39]. Moreover, a polymorphism study in the alpha subunit found that GNAL gene was associated with major depression [Reference Zill, Engel, Baghai, Zwanzger, Schüle and Minov40]. As a potassium-dependent Na+/Ca2+ exchanger, NCKX3 (near CFAP61) is mostly abundant in the brain [Reference Tran, Jung, Yoo, Lee and Jeung41]. Behavioral examination in NCKX3 knock-out mice found that depression-related behaviors in Nckx−/− mice were more higher than that in wild type mice [Reference Tran, Jung, Yoo, Lee and Jeung41]. Thus, combined with our findings, the SNPs involved in this study may affect the pathogenesis of depression through regulating these adjacent genes.

The potential links between depression and obesity or smoking have been well speculated and repeatedly examined [Reference Luppino, de Wit, Bouvy, Stijnen, Cuijpers and Penninx42, Reference Mathew, Hogarth, Leventhal, Cook and Hitsman43]. Depression and obesity are common diseases with major public health implications, often cooccurring within individuals [Reference Milaneschi, Simmons, van Rossum and Penninx44]. The current study also demonstrated that there existed significant genetic correlation between non-anxiety depression and obesity-related traits, including overweight, obesity class 2 and impedance of arm. Obesity can be seen as an inflammatory state because weight gain has been shown to activate inflammatory pathways, which in turn has been linked to depression [Reference Shoelson, Herrero and Naaz45Reference Bremmer, Beekman, Deeg, Penninx, Dik and Hack47]. In addition, smoking has been found to phenotypically and genetically correlate with depression [Reference Mathew, Hogarth, Leventhal, Cook and Hitsman43, Reference Fluharty, Taylor, Grabski and Munafò48]. The current study also identified significant genetic correlations between non-anxiety depression and smoking-related traits, including current tobacco smoking and current smoking status. There have been many studies using different methods to demonstrate the complex causal relationship between depression and smoking, such as smoking increasing the risk of depression [Reference Boden, Fergusson and Horwood49], a bidirectional effect [Reference Wootton, Richmond, Stuijfzand, Lawn, Sallis and Taylor50] and no effect reported [Reference Wium-Andersen, Ørsted and Nordestgaard51]. Our results largely support the previous association between depression and obesity or smoking.

Many epidemiological studies convey the same message: divorce threatens mental and physical health [Reference Holt-Lunstad, Smith and Layton52, Reference Sbarra, Law and Portley53]. There was a high genetic correlation between non-anxiety depression and marital separation/divorce (rg = 0.45) in our LDSC. It could be valuable to assess depressive phenotypes from a large cohort associated with access to marital status records. Separated and divorced people have a higher risk of mental illness than married people [Reference Sbarra and Nietert54]. An earlier cross-sectional study focused on immune changes and found that separated/divorced men had poorer immune function and more recent illnesses [Reference Kiecolt-Glaser, Kennedy, Malkoff, Fisher, Speicher and Glaser55]. In fact, the norepinephrine data for newlyweds matched the evidence linking divorce with increased inflammation [Reference Sbarra56]. In addition, those who had symptoms of depression before their divorce were more likely to develop depression after their divorce [Reference Sbarra, Emery, Beam and Ocker57]. Future research should focus on how marriage and divorce can provoke health-relevant immune alterations, especially on immune and hormone-related non-anxiety depression.

Overall health rating is a common self-assessment score in epidemiological studies, and has been widely used as a powerful indicator in public health studies [Reference Wennberg, Rolandsson, van der A, Spijkerman, Kaaks and Boeing58Reference Singh-Manoux, Martikainen, Ferrie, Zins, Marmot and Goldberg60]. We observed moderate genetic correlation between non-anxiety depression and overall health rating (rg = 0.21), which was consistent with a similar research assessing the relationship between overall health rating score and major depressive disorder [Reference Kim and Choi61]. This genetic correlation for overall health rating may offer important new insights into the interrelationship between self-assessment health and depression. Future studies need to assess the pleiotropic confounding factors that could further explain genetic correlation between non-anxiety depression and overall health rating. In these studies, phenotypes should ideally be investigated without composite compositions and based on multiple phenotype indicators.

The principal strength of this study is to conduct a more accurate study on the premise of excluding the influence of anxiety comorbidities of depression in two large cohorts. This has allowed the validation of the effects of variants that have been identified previously to determine whether they maintain an effect on depression. We also recognize limitations in our study. In contrast to Europeans, other ancestors are still less well studied. We hope that the initial results of the UK sample reported here can help to advance the field by encouraging more collaborative research on other ethnic groups. Besides, there are some non-UK ethnic samples in LD Hub, which may have a slight bias on the accuracy of genetic correlation analysis, and future studies could further refine the phenotype.

Conclusion

In summary, we identified two novel loci that should prioritize their further study in the pathology of non-anxiety depression. We examined genetic correlations between non-anxiety depression and 855 health-related phenotypes, largely confirming and strengthening previous observations. This study sheds light on the genetic architecture of depression and provides novel insights for future research of complex comorbid psychiatric traits.

Supplementary Materials

To view supplementary material for this article, please visit http://doi.org/10.1192/j.eurpsy.2022.32.

Data Availability Statement

The data that support the findings of this study are available from the authors.

Author Contributions

Conceptualization: F.Z.; Data curation: B.C.; Formal analysis: X.Q., P.M., S.C.; Investigation: B.C., P.M.; Methodology: S.C., Y.J., F.Z.; Resources: X.Y.; Software: S.C., X.Y.; Supervision: Y.W., F.Z.; Validation: X.Q., P.M., S.C., X.Y., L.L., Y.Y.; Visualization: B.C., L.L., Y.Y.; Writing—original draft: B.C.; Writing—review and editing: B.C., X.Q., Y.J., Y.W., F.Z.

Financial Support

This work was supported by the National Natural Scientific Foundation of China (81922059), the Natural Science Basic Research Plan in Shaanxi Province of China (2021JCW-08), and the Fundamental Research Funds for the Central Universities.

Conflicts of Interest

The authors declare none.

Footnotes

B.C. and X.Q. contributed equally to this work.

References

Pandi-Perumal, SR, Monti, JM, Burman, D, Karthikeyan, R, BaHammam, AS, Spence, DW, et al. Clarifying the role of sleep in depression: a narrative review. Psychiatry Res. 2020;291:113239.CrossRefGoogle ScholarPubMed
Villas Boas, GR, Boerngen de Lacerda, R, Paes, MM, Gubert, P, Almeida, W, Rescia, VC, et al. Molecular aspects of depression: a review from neurobiology to treatment. Eur J Pharmacol. 2019;851:99121.CrossRefGoogle ScholarPubMed
Cheng, B, Chu, X, Yang, X, Wen, Y, Jia, Y, Liang, C, et al. Dietary habit is associated with depression and intelligence: an observational and genome-wide environmental interaction analysis in the UK biobank cohort. Nutrients. 2021;13:1150.CrossRefGoogle ScholarPubMed
Sullivan, PF, Neale, MC, Kendler, KS. Genetic epidemiology of major depression: review and meta-analysis. Am J Psychiatry. 2000;157:1552–62.CrossRefGoogle ScholarPubMed
Ripke, S, Wray, NR, Lewis, CM, Hamilton, SP, Weissman, MM, Breen, G, et al. A mega-analysis of genome-wide association studies for major depressive disorder. Mol Psychiatry. 2013;18:497511.Google ScholarPubMed
Hek, K, Demirkan, A, Lahti, J, Terracciano, A, Teumer, A, Cornelis, MC, et al. A genome-wide association study of depressive symptoms. Biol Psychiatry. 2013;73:667–78.CrossRefGoogle ScholarPubMed
Wray, NR, Ripke, S, Mattheisen, M, Trzaskowski, M, Byrne, EM, Abdellaoui, A, et al. Genome-wide association analyses identify 44 risk variants and refine the genetic architecture of major depression. Nat Genetics. 2018;50:668–81.CrossRefGoogle ScholarPubMed
Levey, DF, Stein, MB, Wendt, FR, Pathak, GA, Zhou, H, Aslan, M. Bi-ancestral depression GWAS in the Million Veteran Program and meta-analysis in >1.2 million individuals highlight new therapeutic directions. Nat Neurosci. 2021;24:954–63.CrossRefGoogle Scholar
Anttila, V, Bulik-Sullivan, B, Finucane, HK, Walters, RK, Bras, J, Duncan, L. Analysis of shared heritability in common disorders of the brain. Science. 2018;360:eaap8757.Google ScholarPubMed
Bartlett, AA, Singh, R, Hunter, RG. Anxiety and Epigenetics. Adv Exp Med Biol. 2017;978:145–66.CrossRefGoogle ScholarPubMed
Jacobson, NC, Newman, MG. Anxiety and depression as bidirectional risk factors for one another: a meta-analysis of longitudinal studies. Psychol Bull. 2017;143:1155–200.CrossRefGoogle ScholarPubMed
Waszczuk, MA, Zavos, HM, Gregory, AM, Eley, TC. The phenotypic and genetic structure of depression and anxiety disorder symptoms in childhood, adolescence, and young adulthood. JAMA Psychiatry. 2014;71:905–16.CrossRefGoogle Scholar
Craske, MG, Rauch, SL, Ursano, R, Prenoveau, J, Pine, DS, Zinbarg, RE. What is an anxiety disorder? Depress Anxiety. 2009;26:1066–85.CrossRefGoogle ScholarPubMed
Morneau-Vaillancourt, G, Coleman, JRI, Purves, KL, Cheesman, R, Rayner, C, Breen, G, et al. The genetic and environmental hierarchical structure of anxiety and depression in the UK Biobank. Depress Anxiety. 2020;37:512–20.CrossRefGoogle ScholarPubMed
van Dalen, M, Dierckx, B, Pasmans, S, Aendekerk, EWC, Mathijssen, IMJ, Koudstaal, MJ, et al. Anxiety and depression in adolescents with a visible difference: a systematic review and meta-analysis. Body Image. 2020;33:3846.CrossRefGoogle ScholarPubMed
Boyer, P. Do anxiety and depression have a common pathophysiological mechanism? Acta Psychiatr Scand Suppl. 2000;102:24–9.CrossRefGoogle Scholar
Gao, W, Ping, S, Liu, X. Gender differences in depression, anxiety, and stress among college students: a longitudinal study from China. J Affect Disord. 2020;263:292300.CrossRefGoogle ScholarPubMed
Bycroft, C, Freeman, C, Petkova, D, Band, G, Elliott, LT, Sharp, K, et al. The UK Biobank resource with deep phenotyping and genomic data. Nature. 2018;562:203–9.CrossRefGoogle ScholarPubMed
Davis, KAS, Cullen, B, Adams, M, Brailean, A, Breen, G, Coleman, JRI, et al.. Indicators of mental disorders in UK Biobank: a comparison of approaches. Int J Methods Psychiatr Res. 2019;28:e1796.CrossRefGoogle ScholarPubMed
Canela-Xandri, O, Rawlik, K, Tenesa, A. An atlas of genetic associations in UK Biobank. Nature Genet. 2018;50:1593–9.CrossRefGoogle ScholarPubMed
Purcell, S, Neale, B, Todd-Brown, K, Thomas, L, Ferreira, MA, Bender, D, et al. PLINK: a tool set for whole-genome association and population-based linkage analyses. Am J Hum Genet. 2007;81:559–75.CrossRefGoogle ScholarPubMed
Howard, DM, Adams, MJ, Clarke, T-K, Hafferty, JD, Gibson, J, Shirali, M, et al. Genome-wide meta-analysis of depression identifies 102 independent variants and highlights the importance of the prefrontal brain regions. Nat Neurosci. 2019;22:343–52.CrossRefGoogle ScholarPubMed
Bulik-Sullivan, B, Finucane, HK, Anttila, V, Gusev, A, Day, FR, Loh, P-R, et al. An atlas of genetic correlations across human diseases and traits. Nat Genet. 2015;47:1236–41.CrossRefGoogle ScholarPubMed
Zheng, J, Erzurumluoglu, AM, Elsworth, BL, Kemp, JP, Howe, L, Haycock, PC, et al. LD Hub: a centralized database and web interface to perform LD score regression that maximizes the potential of summary level GWAS data for SNP heritability and genetic correlation analysis. Bioinformatics. 2017;33:272–9.CrossRefGoogle ScholarPubMed
Brailean, A, Curtis, J, Davis, K, Dregan, A, Hotopf, M. Characteristics, comorbidities, and correlates of atypical depression: evidence from the UK Biobank Mental Health Survey. Psychol Med. 2020;50:1129–38.CrossRefGoogle ScholarPubMed
Gustavson, DE, Franz, CE, Panizzon, MS, Reynolds, CA, Xian, H, Jacobson, KC, et al. Genetic and environmental associations among executive functions, trait anxiety, and depression symptoms in middle age. Clin Psychol Sci. 2019;7:127–42.CrossRefGoogle ScholarPubMed
Alcaino, C, Knutson, KR, Treichel, AJ, Yildiz, G, Strege, PR, Linden, DR, et al. A population of gut epithelial enterochromaffin cells is mechanosensitive and requires Piezo2 to convert force into serotonin release. Proc Natl Acad Sci U S A 2018;115:E7632–e41.CrossRefGoogle ScholarPubMed
Brummelte, S, Mc Glanaghy, E, Bonnin, A, Oberlander, TF. Developmental changes in serotonin signaling: Implications for early brain function, behavior and adaptation. Neuroscience. 2017;342:212–31.CrossRefGoogle ScholarPubMed
Pardo-Pastor, C, Rubio-Moscardo, F, Vogel-González, M, Serra, SA, Afthinos, A, Mrkonjic, S, et al. Piezo2 channel regulates RhoA and actin cytoskeleton to promote cell mechanobiological responses. Proc Natl Acad Sci USA. 2018;115:1925–30.CrossRefGoogle ScholarPubMed
Umans, BD, Liberles, SD. Neural sensing of organ volume. Trends Neurosci. 2018;41:911–24.CrossRefGoogle ScholarPubMed
Wang, J, Hamill, OP. Piezo2, A pressure sensitive channel is expressed in select neurons of the mouse brain: a putative mechanism for synchronizing neural networks by transducing intracranial pressure pulses. Research Square; 2021.CrossRefGoogle Scholar
Chesler, AT, Szczot, M, Bharucha-Goebel, D, Ceko, M, Donkervoort, S, Laubacher, C, et al. The role of PIEZO2 in human mechanosensation. N Engl J Med. 2016;375:1355–64.CrossRefGoogle ScholarPubMed
Heyburn, L, Abutarboush, R, Goodrich, S, Urioste, R, Batuure, A, Statz, J, et al. Repeated low-level blast overpressure leads to endovascular disruption and alterations in TDP-43 and Piezo2 in a rat model of blast TBI. Front Neurol. 2019;10:766.CrossRefGoogle Scholar
Liu, S, Trupiano, MX, Simon, J, Guo, J, Anton, ES. The essential role of primary cilia in cerebral cortical development and disorders. Curr Top Dev Biol. 2021;142:99146.CrossRefGoogle ScholarPubMed
Rouillard, AD, Gundersen, GW, Fernandez, NF, Wang, Z, Monteiro, CD, McDermott, MG, et al. The harmonizome: a collection of processed datasets gathered to serve and mine knowledge about genes and proteins. Database (Oxford) 2016;2016:baw100.CrossRefGoogle ScholarPubMed
Zhu, Z, Zhang, F, Hu, H, Bakshi, A, Robinson, MR, Powell, JE, et al. Integration of summary data from GWAS and eQTL studies predicts complex trait gene targets. Nat Genet. 2016;48:481–7.CrossRefGoogle ScholarPubMed
Won, H, de la Torre-Ubieta, L, Stein, JL, Parikshak, NN, Huang, J, Opland, CK, et al. Chromosome conformation elucidates regulatory relationships in developing human brain. Nature. 2016;538:523–7.CrossRefGoogle ScholarPubMed
Whalen, S, Truty, RM, Pollard, KS. Enhancer-promoter interactions are encoded by complex genomic signatures on looping chromatin. Nat Genet. 2016;48:488–96.CrossRefGoogle ScholarPubMed
Tombácz, D, Maróti, Z, Kalmár, T, Csabai, Z, Balázs, Z, Takahashi, S, et al. High-coverage whole-exome sequencing identifies candidate genes for suicide in victims with major depressive disorder. Sci Rep. 2017;7:7106.CrossRefGoogle ScholarPubMed
Zill, P, Engel, R, Baghai, TC, Zwanzger, P, Schüle, C, Minov, C, et al. Analysis of polymorphisms in the olfactory G-protein Golf in major depression. Psychiatr Genet. 2002;12:1722.CrossRefGoogle ScholarPubMed
Tran, DN, Jung, EM, Yoo, YM, Lee, JH, Jeung, EB. Potassium-dependent sodium/calcium exchanger 3 (Nckx3) depletion leads to abnormal motor function and social behavior in mice. J Physiol Pharmacol. 2020;71:525–36.Google ScholarPubMed
Luppino, FS, de Wit, LM, Bouvy, PF, Stijnen, T, Cuijpers, P, Penninx, BW, et al. Overweight, obesity, and depression: a systematic review and meta-analysis of longitudinal studies. Arch Gen Psychiatry. 2010;67:220–9.CrossRefGoogle ScholarPubMed
Mathew, AR, Hogarth, L, Leventhal, AM, Cook, JW, Hitsman, B. Cigarette smoking and depression comorbidity: systematic review and proposed theoretical model. Addiction. 2017;112:401–12.CrossRefGoogle ScholarPubMed
Milaneschi, Y, Simmons, WK, van Rossum, EFC, Penninx, BW. Depression and obesity: evidence of shared biological mechanisms. Mol Psychiatry. 2019;24:1833.CrossRefGoogle ScholarPubMed
Shoelson, SE, Herrero, L, Naaz, A. Obesity, inflammation, and insulin resistance. Gastroenterology. 2007;132:2169–80.CrossRefGoogle Scholar
Vaccarino, V, Johnson, BD, Sheps, DS, Reis, SE, Kelsey, SF, Bittner, V, et al. Depression, inflammation, and incident cardiovascular disease in women with suspected coronary ischemia: the National Heart, Lung, and Blood Institute-sponsored WISE study. J Am Coll Cardiol. 2007;50:2044–50.CrossRefGoogle ScholarPubMed
Bremmer, MA, Beekman, AT, Deeg, DJ, Penninx, BW, Dik, MG, Hack, CE, et al. Inflammatory markers in late-life depression: results from a population-based study. J Affect Disord. 2008;106:249–55.CrossRefGoogle ScholarPubMed
Fluharty, M, Taylor, AE, Grabski, M, Munafò, MR. The association of cigarette smoking with depression and anxiety: a systematic review. Nicotine Tob Res 2017;19:313.CrossRefGoogle ScholarPubMed
Boden, JM, Fergusson, DM, Horwood, LJ. Cigarette smoking and depression: tests of causal linkages using a longitudinal birth cohort. Br J Psychiatry 2010;196:440–6.CrossRefGoogle ScholarPubMed
Wootton, RE, Richmond, RC, Stuijfzand, BG, Lawn, RB, Sallis, HM, Taylor, GMJ, et al. Evidence for causal effects of lifetime smoking on risk for depression and schizophrenia: a Mendelian randomisation study. Psychol Med. 2020;50:2435–43.CrossRefGoogle ScholarPubMed
Wium-Andersen, MK, Ørsted, DD, Nordestgaard, BG. Tobacco smoking is causally associated with antipsychotic medication use and schizophrenia, but not with antidepressant medication use or depression. Int J Epidemiol. 2015;44:566–77.CrossRefGoogle ScholarPubMed
Holt-Lunstad, J, Smith, TB, Layton, JB. Social relationships and mortality risk: a meta-analytic review. PLoS Med. 2010;7:e1000316.CrossRefGoogle ScholarPubMed
Sbarra, DA, Law, RW, Portley, RM. Divorce and death: a meta-analysis and research agenda for clinical, social, and health psychology. Perspect Psychol Sci. 2011;6:454–74.CrossRefGoogle ScholarPubMed
Sbarra, DA, Nietert, PJ. Divorce and death: forty years of the Charleston heart study. Psychol Sci. 2009;20:107–13.CrossRefGoogle ScholarPubMed
Kiecolt-Glaser, JK, Kennedy, S, Malkoff, S, Fisher, L, Speicher, CE, Glaser, R. Marital discord and immunity in males. Psychosom Med. 1988;50:213–29.CrossRefGoogle ScholarPubMed
Sbarra, DA. Divorce and health: current trends and future directions. Psychosom Med 2015;77:227–36.CrossRefGoogle ScholarPubMed
Sbarra, DA, Emery, RE, Beam, CR, Ocker, BL. Marital dissolution and major depression in Midlife: a propensity score analysis. Clin Psychol Sci. 2014;2:249–57.CrossRefGoogle ScholarPubMed
Wennberg, P, Rolandsson, O, van der A, DL, Spijkerman, AM, Kaaks, R, Boeing, H, et al. Self-rated health and type 2 diabetes risk in the European prospective investigation into cancer and nutrition-interact study: a case-cohort study. BMJ Open. 2013;3:e002436.CrossRefGoogle ScholarPubMed
van der Linde, RM, Mavaddat, N, Luben, R, Brayne, C, Simmons, RK, Khaw, KT, et al. Self-rated health and cardiovascular disease incidence: results from a longitudinal population-based cohort in Norfolk, UK. PLoS One. 2013;8:e65290.CrossRefGoogle ScholarPubMed
Singh-Manoux, A, Martikainen, P, Ferrie, J, Zins, M, Marmot, M, Goldberg, M. What does self rated health measure? Results from the British Whitehall II and French Gazel cohort studies. J Epidemiol Community Health 2006;60:364–72.CrossRefGoogle ScholarPubMed
Kim, J, Choi, Y. Multidimensional factors of determining trajectories in self-rated health and depressive symptoms among Korean children and adolescents. Soc Work Health Care. 2018;57:717–36.CrossRefGoogle ScholarPubMed
Figure 0

Figure 1. Manhattan plot for the GWAS of depression without anxiety in the UK Biobank cohorts. (A) Linear regression of depression score in anxiety score <5 samples. (B) Linear regression of depression score in non-self-reported anxiety samples. (C) Logistic regression of self-reported depression in anxiety score <5 samples. (D) Logistic regression of self-reported depression in non-self-reported anxiety samples. The red line indicates the p-value threshold for genome-wide significance (p < 5 × 10−8).

Figure 1

Figure 2. LocusZoom plots of depression without anxiety genome-wide significance loci. Association results for SNPs as a function of genomic distance for PIEZO2 and CFAP61 (C20orf26). The top line in each subfigure shows genomic coverage at the locus, with each vertical tick representing the imputed SNPs. Purple diamond indicate SNP at the locus with the strongest association evidence. Each point represents a SNP. Bottom panel shows genes at each locus as annotated in the UCSC Genome Browser Annotation Database. (A) display PIEZO2 in chr18 for GWAS summary of depression score in anxiety score <5 samples. (B) display PIEZO2 in chr18 for GWAS summary of depression score in non-self-reported anxiety samples. (C) display CFAP61 in chr20 for GWAS summary of self-reported depression in anxiety score <5 samples. (D) display CFAP61 in chr20 for GWAS summary of self-reported depression in non-self-reported anxiety samples.

Figure 2

Figure 3. Genetic correlations and mental disorders related traits using LD score regression implemented in LD Hub software. The negative rg indicates that an earlier or lower value of a continuous trait was associated with depression without anxiety. The positive rg indicates that a later or higher value of a continuous trait was associated with depression without anxiety.

Figure 3

Figure 4. Significant genetic correlations and other behavioral and disease related traits using LD score regression implemented in LD Hub software. The negative rg indicates that an earlier or lower value of a continuous trait was associated with depression without anxiety. The positive rg indicates that a later or higher value of a continuous trait was associated with depression without anxiety.

Supplementary material: File

Cheng et al. supplementary material

Cheng et al. supplementary material

Download Cheng et al. supplementary material(File)
File 196.4 KB
Submit a response

Comments

No Comments have been published for this article.