Skip to main content

Investigating the relationship between depression and breast cancer: observational and genetic analyses

A Commentary to this article was published on 06 June 2023

Abstract

Background

Both depression and breast cancer (BC) contribute to a substantial global burden of morbidity and mortality among women, and previous studies have observed a potential depression-BC link. We aimed to comprehensively characterize the phenotypic and genetic relationships between depression and BC.

Methods

We first evaluated phenotypic association using longitudinal follow-up data from the UK Biobank (N = 250,294). We then investigated genetic relationships leveraging summary statistics from the hitherto largest genome-wide association study of European individuals conducted for depression (N = 500,199), BC (N = 247,173), and its subtypes based on the status of estrogen receptor (ER + : N = 175,475; ER − : N = 127,442).

Results

Observational analysis suggested an increased hazard of BC in depression patients (HR = 1.10, 95%CIs = 0.95–1.26). A positive genetic correlation between depression and overall BC was observed (\({r}_{g}\) = 0.08, P = 3.00 × 10–4), consistent across ER + (\({r}_{g}\) = 0.06, P = 6.30 × 10–3) and ER − subtypes (\({r}_{g}\) = 0.08, P = 7.20 × 10–3). Several specific genomic regions showed evidence of local genetic correlation, including one locus at 9q31.2, and four loci at, or close, to 6p22.1. Cross-trait meta-analysis identified 17 pleiotropic loci shared between depression and BC. TWAS analysis revealed five shared genes. Bi-directional Mendelian randomization suggested risk of depression was causally associated with risk of overall BC (OR = 1.12, 95%Cis = 1.04–1.19), but risk of BC was not causally associated with risk of depression.

Conclusions

Our work demonstrates a shared genetic basis, pleiotropic loci, and a putative causal relationship between depression and BC, highlighting a biological link underlying the observed phenotypic relationship; these findings may provide important implications for future studies aimed reducing BC risk.

Peer Review reports

Background

Individuals with mental health disorders are often at an increased risk for an array of subsequent medical conditions including cancer [1]. One example of this phenotypic association is depression and breast cancer (BC), both of which disproportionately or primarily affect women and contribute to a substantial global burden of morbidity and mortality [2, 3]. Hypothesized underlying biological mechanisms include enhanced inflammation and oxidative stress, inhibited immune surveillance, and dysfunctional activation of the autonomic nervous system and the hypothalamic–pituitary–adrenal axis [1]. Population-based evidence regarding the depression-BC relationship, however, remains inconsistent. While one meta-analysis aggregating data from 11 cohort studies involving 182,241 individuals identified an increased (though non-significant) risk of BC among individuals with depression (pooled relative risk (RR) = 1.13, 95% confidence intervals (95%CIs) = 0.94–1.36) [4], another meta-analysis restricted to cohort studies with follow-up periods longer than ten years reported a significant association (pooled RR = 2.50, 95%CIs = 1.06–5.91) [5]. Nevertheless, phenotypic correlations derived from observational studies can be subject to bias, confounding, and reverse causality [6].

One way to disentangle these conflicting findings is to investigate the potential genetic underpinnings of comorbid disorders. Twin studies have established that both depression and BC are under genetic influence, with heritability estimates of 37% and 31%, respectively [7, 8]. More recently, large-scale genome-wide analysis of common variants have found significant genetic correlations between depression and multiple female reproductive phenotypes [9,10,11] including age at menarche and age at natural menopause (well-established risk factors for BC [12]). Multiple loci have further been identified (assessed in GWAS catalog [13] on Jan 15, 2023) as influencing both traits (i.e., TENM2, BTN2A1, ESR1, ASTN2, SLC6A15). These results suggest that depression and BC may be linked by shared biology, though the extent and nature of such links remains unclear.

Recent advances in statistical genetics have yielded a range of methods to enable comprehensive genome-wide cross-trait analyses to characterize the shared and distinct genetic influences across traits, driving forward epidemiologic associations with novel insights into the underlying biological mechanisms [6]. Here we apply these methods to dissect the genetic and phenotypic relationships between depression and BC. Specifically, leveraging the hitherto largest observational and genetic data, we quantified phenotypic association, global and local genetic correlations, pleiotropic loci, and potential causal relationships. Given the role estrogen receptor (ER) plays in the pathogenesis of depression [14], we further investigated the shared genetic architecture with depression across BC subtypes characterized by distinct ER status. Figure 1 illustrates the overall design of the study, with Additional file 1: Figure S1 providing an accompanying graphical abstract.

Fig. 1
figure 1

Flowchart of the overall study design. We first used longitudinal follow-up data from the UK Biobank to evaluate the phenotypic association between depression and breast cancer. We then leveraged summary statistics from the hitherto largest genome-wide association study conducted for each trait to characterize the shared genetic architecture through a genome-wide cross-trait analysis. BC: breast cancer; ER: estrogen receptor; GWAS: genome-wide association study; GTEx: Genotype-Tissue Expression project; MR: Mendelian randomization

Methods

We first used data from UK Biobank (UKBB) to evaluate the phenotypic association. We then leveraged summary statistics from genome-wide association study (GWAS) conducted for each trait to characterize the shared genetic architecture through a genome-wide cross-trait analysis. The Strengthening the Reporting of Observational Studies in Epidemiology (STROBE) checklists for both observational and genetic studies are detailed in Additional file 2.

Data sources

UK Biobank data

UKBB is a large-scale prospective cohort study with more than 500,000 UK participants (55% women) aged 40–69 years when recruited in 2006–2010 [15]. At recruitment, all participants gave informed consent for participation and follow-up. Overall, 503,317 participants consented to join the study cohort and visited an assessment center, among which we only considered women of European descent (N = 257,150). We defined a diagnosis of depression as the ICD-10 (international classification of diseases, 10th revision) codes F32, F33, F34, F38, and F39, and a diagnosis of BC as the ICD-10 code C50 and the ICD-9 code 174, using linked records of hospital admission in UKBB (Data-Field: 41,202 and Data-Field: 41,204). After excluding 6,856 participants with a history of BC at baseline, 250,294 participants were finally included.

GWAS datasets

GWAS summary data for depression was obtained from a meta-analysis of 807,553 individuals (246,363 cases and 561,190 controls, all of European ancestry) combining the three largest existing genetic studies of depression (UKBB, 23andMe and the Psychiatric Genomics Consortium (PGC)) [9]. The phenotypes ranged from self-reported help-seeking for problems with “nerves, anxiety, tension or depression” (termed “broad depression”) (51.8%), self-reported clinical diagnosis of major depression (MD) (30.7%), and clinically ascertained diagnosis of MD (17.5%). Independent trait-associated genome-wide significant single nucleotide polymorphisms (SNPs) were identified at a P-threshold of 5 × 10–8, after clumping SNPs in linkage disequilibrium (LD) (\({r}^{2}\) > 0.10 across a 3.0 Mb window) with the top-associated SNPs.

We extracted the effect size and relevant information of the 102 GWAS-identified depression-associated significant SNPs and used those SNPs as instrumental variables (IVs) in Mendelian randomization (MR) analyses [16] (Additional file 3: Table S1). We also retrieved the full set of summary statistics for depression, in which data from 23andMe were excluded (due to limited data availability).

The largest available GWAS summary data for overall BC was obtained from a meta-analysis of 247,173 individuals (133,384 clinically ascertained BC cases and 113,789 controls, all of European ancestry), combining data from 82 participating studies of the Breast Cancer Association Consortium (BCAC) and 11 other BC genetic studies [17]. Top-associated SNPs in the combined meta-analysis reaching a P-threshold of 5 × 10–8 were reported. This GWAS identified 32 novel susceptibility loci in addition to 153 previously detected loci (for a total of 185 genome-wide significant loci). For subtype-specific BC, the hitherto largest GWAS summary data was obtained from a meta-analysis of 122,977 BC cases (69,501 ER + BC and 21,468 ER − BC) and 105,974 controls (all of European ancestry) combining data from BCAC and DRIVE (Discovery, Biology and Risk of Inherited Variants in Breast Cancer Consortium) [18]. No overlapping participating studies were shared between the BCs and depression GWASs.

We extracted the effect size and relevant information of the 185 GWAS-identified BC-associated significant SNPs for reverse-direction MR analysis. As trait-associated SNPs were not reported for subtype-specific BC, we determined IVs as the lead SNPs reaching genome-wide significance (P < 5 × 10–8) after removing SNPs in LD (\({r}^{2}\) > 0.001 across a 10.0 Mb window). We also downloaded full summary statistics for all BCs (overall and subtypes).

Details on the characteristics of each included data set are presented in Additional file 3: Table S2.

Statistical analysis

Observational analysis

We constructed a Cox proportional hazards regression model with exposure to depression modeled as a time-dependent variable. We used three sets of adjustments to minimize the role of confounding. Estimates in model 1 were adjusted only for age, assessment center, and the top 40 genetic principal components. Estimates in model 2 were adjusted additionally for income, Townsend deprivation index, body mass index (BMI), smoking, drinking, physical activity, sleep duration, and education. Estimates in model 3 were adjusted, on top of model 2, for family history of BC, reproductive factors (parity, age at menarche, menopausal status, use of oral contraceptives, and hormone replacement therapy), diagnosis of other mental health disorders, and treatment/medicine for antidepressants or antipsychotics. All analyses were conducted using SAS version 9.4 (SAS Institute, Cary, NC).

Global genetic correlation analysis

To evaluate an overall shared genetic basis between depression and BC, we estimated their global genetic correlation using cross-trait LD Score Regression (LDSC) [19]. LDSC requires only GWAS summary statistics as input and quantifies the average sharing of genetic effect between pairs of traits. The genetic correlation (\({r}_{g}\)) ranges from –1 to 1, with –1 indicating a perfect negative correlation and 1 indicating a perfect positive correlation. A Bonferroni corrected P-value of 0.017 was used as significant threshold (P < 0.05/3, number of overall BC and subtypes).

Local genetic correlation analysis

To identify genomic regions that contribute disproportionately to the global genetic correlation between depression and BC, we further estimated local genetic correlation in 1,703 pre-defined LD-independent regions using both Heritability Estimator from Summary Statistics (ρ-HESS) and Pairwise-GWAS (GWAS-PW) [20, 21]. While ρ-HESS provides a precise quantification of each genomic region, GWAS-PW uses a Bayesian framework to calculate the posterior probabilities of association (PPA) of a genomic region under four models [21]. Significant local signal was determined if the P-value from ρ-HESS survived multiple corrections (Pρ-HESS < 0.05/1703), or if PPA of model 3 (PPA _3, the probability that a genomic region associated with both traits) from GWAS-PW was larger than 0.5.

Cross-trait meta-analysis

Significant genetic correlation either suggest genetic variants having independent effect on both traits (pleiotropy) or genetic variants influencing one trait via its effect on the other (causality). We next performed a cross-trait meta-analysis using Cross-Phenotype Association (CPASSOC) to identify pleiotropic loci [22]. Based on the fixed-effect model, CPASSOC provides two estimates, SHom and SHet, to combine summary statistics across traits while controlling population structure and cryptic relatedness. SHom is the most powerful when genetic effect sizes are homogenous, which is unlikely to be true when meta-analyzing multiple traits. As an extension of SHom, SHet maintains statistical power even under the presence of heterogeneity by assigning more weights to the larger trait-specific effect sizes, and was thus adopted for the analysis herein.

We applied PLINK’s “clumping” function to obtain independent loci (parameters: –clump-p1 5e-8 –clump-p2 1e-5 –clump-r2 0.2 –clump-kb 500) [23]. Within each locus, the variant with the lowest P-value was kept as index SNP. Significant pleiotropic SNPs were defined as index variants satisfying PCPASSOC < 5 × 10–8 and Psingle-trait < 1 × 10–3 (for both traits). Novel pleiotropic SNPs were defined as significant pleiotropic SNPs that did not reach genome-wide significance in single-trait GWASs (5 × 10–8 < Psingle-trait < 1 × 10–3), were independent (\({r}^{2}\)< 0.20) of previously reported genome-wide significant SNPs (of depression and BC), and none of their neighboring SNPs (± 500 kb) reached P < 5 × 10–8 in single-trait GWASs.

We used Ensembl Variant Effect Predictor (VEP) [24] and 3DSNP [25] for detailed functional annotation of the identified pleiotropic SNPs. While VEP selects candidate genes based on simple physical proximity [24], 3DSNP annotates the regulatory function of SNPs by exploring their 3D interactions with genes mediated by chromatin loops [25].

Fine mapping credible set analysis

Because an index SNP is not likely to be the causal SNP, we further identified a credible set of variants that were 99% likely to contain causal variants at each of the pleiotropic loci through FM-summary (https://github.com/hailianghuang/FM-summary) [26]. Briefly, FM-summary is a Bayesian fine-mapping algorithm that maps only the primary signal and uses a flat prior with steepest descent approximation, assuming at least one causal variant exists within a given region.

Colocalization analysis

To investigate whether the same variants are responsible for two GWAS signals as opposed to distinct genetic variants in proximity, we next performed a colocalization analysis using Coloc [27]. This method provides posterior probabilities for five mutually exclusive hypotheses regarding the sharing of causal variants in a genomic region. We extracted summary statistics for variants within 500 kb of each shared index SNP, and calculated the posterior probability for H4 (PPH4, the probability that both traits associated through sharing a single causal variant). A locus was considered colocalized if PPH4 was greater than 0.5.

Transcriptome-wide association studies

Many genetic variants influence complex traits by modulating gene expression; thus, investigating overlapped genes may help clarify causal mechanisms. To identify associations between depression and BC with regard to gene expression in specific tissues, we conducted a transcriptome-wide association study (TWAS) using FUSION [28] by integrating GWAS summary data with expression weights across 49 tissues from GTEx (Genotype-Tissue Expression, version 8). We first performed 49 TWASs for each trait, one tissue-trait pair at a time, to identify an independent set of gene-tissue pairs. We then intersected single-trait TWAS results to examine if there were shared gene-tissue pairs across traits. The Bonferroni correction was used within each tissue to account for multiple comparisons.

Mendelian randomization analysis

We finally performed a comprehensive two-sample MR analysis to detect potential causal relationships. We used the random-effect inverse-variance weighted (IVW) approach as the primary approach [29]. This method pools the Wald ratio estimate of each SNP and obtains the average casual effect estimate between two traits. To reduce biased estimate due to pleiotropic effects of genetic instruments, we adopted two complementary methods: MR-Egger regression (detects and corrects for bias due to directional pleiotropy) [30], and weighted median approach (provides a consistent causal estimation with even ≥ 50% invalid IVs) [31]. As these two approaches were less powerful than IVW in detecting true causal effects, we defined a significant causal estimate as significant in IVW (P-value < 0.05/3, number of overall BC and subtypes) and showing directional consistency in MR-Egger regression and weighted median approach.

Additional sensitivity analyses were conducted to validate MR model assumptions (i.e., relevance, independence, and exclusion restriction) [32], including: (i) exclusion of palindromic IVs with strand ambiguity; (ii) exclusion of pleiotropic IVs associated with potential confounding traits (confirmed by GWAS Catalog [13]); (iii) leave-one-out analysis where one SNP was removed at a time and IVW was conducted based on the remaining SNPs; and (iv) MR-Pleiotropy Residual Sum and Outlier (MR-PRESSO) to evaluate the presence of horizontal pleiotropy and to re-calculate causal effects after removing the detected outliers [33]. A multi-variable MR (MVMR) approach [34] was further employed to account for the effect from major confounders, including BMI [35], smoking initiation [36], alcohol consumption [36], physical activity [37], sleep duration [38], and educational attainment [39], where confounders were incorporated together with depression, one at a time as well as simultaneously. To evaluate whether genetically predicted BC exerts a causal effect on depression, a reverse-direction MR was performed where BC-associated independent SNPs were used as IVs.

MR analyses were conducted using packages “TwoSampleMR” (version 0.5.4), “MRPRESSO” (version 1.0), and “MendelianRandomization” (version 0.7.0) in software R (version 4.1.0).

Results

Phenotypic association

Baseline characteristics of UKBB participants included in the observational analysis are presented in Additional file 3: Table S3. In total, participants were followed for 3,014,168 person-years (11.41 ± 2.95 years), during which 529 depression patients and 9,516 depression-free individuals developed BC (Table 1). After adjusting for age, assessment center, and the top 40 genetic principal components, we observed a positive association between depression and risk of BC (hazard ratio (HR) = 1.075, 95%CIs = 0.963–1.199). The effect was strengthened when additionally adjusted for income, Townsend deprivation index, and lifestyle-related factors (HR = 1.090, 95%CIs = 0.951–1.250). In the fully adjusted model, the effect stabilized to 1.096 (95%CIs = 0.951–1.264).

Table 1 Observational association between depression and breast cancer

Global and local genetic correlation

A limited but significant positive genetic correlation was found for depression with overall BC (\({r}_{g}\) = 0.08, P = 3.00 × 10–4), which was also observed in both ER + (\({r}_{g}\) = 0.06, P = 6.30 × 10–3) and ER − subtypes (\({r}_{g}\) = 0.08, P = 7.20 × 10–3), all withstanding Bonferroni correction (Table 2).

Table 2 Global genetic correlation between depression and breast cancer

When the whole genome was partitioned, one genomic region at 6p22.1 (harboring a previous-reported depression locus ZSCAN12) was identified as contributing a significant local genetic correlation to depression and overall BC (Pρ-HESS = 1.83 × 10–5). When examining cancer subtypes, a region at 6p22.2 (harboring a previous-reported depression locus ABT1) was identified as contributing to both depression and ER + BC (Pρ-HESS = 7.30 × 10–6), which was further replicated by GWAS-PW (PPA_3 = 0.73). Additionally, GWAS-PW identified three regions showing high probability of shared associations between depression and ER + BC, including two regions near 6p22.2 (6p22.1, PPA_3 = 0.71; 6p22.3–22.2, PPA_3 = 0.89), and one region at 9q31.2 (PPA_3 = 0.52, harboring a previous-reported BC locus KLF4) (Fig. 2). No significant shared region was found for depression with ER − BC.

Fig. 2
figure 2

Local genetic correlation between depression and breast cancer. Red bars in Manhattan plots and red dots in QQ plots represent loci showing significant local genetic correlation with ρ-HESS after multiple testing adjustments (P < 0.05/1,703). Blue bars in Manhattan plots and blue dots in QQ plots represent additional loci showing a high probability of shared associations using GWAS-PW (PPA _3 > 0.5). a Manhattan plot showing the local genetic correlation, genetic covariance, and SNP-heritability estimated from ρ-HESS between depression and overall breast cancer. b QQ plot presenting region-specific P-values for ρ-HESS-estimated local genetic correlation between depression and overall breast cancer. c Manhattan plot showing the local genetic correlation, genetic covariance, and SNP-heritability estimated from ρ-HESS between depression and ER + breast cancer. d QQ plot presenting region-specific P-values for ρ-HESS-estimated local genetic correlation between depression and ER + breast cancer. DEP: depression; BC: breast cancer; ER: estrogen receptor

Cross-trait meta-analysis and pleiotropic loci

Motivated by the significant genetic overlap observed for depression and BC, we continued to explore at the individual variant level. In total, 17 independent loci reached genome-wide significance in CPASSOC (PCPASSOC < 5 × 10–8 and Pdepression < 1 × 10–3 and PBC < 1 × 10–3), including 12 loci shared between depression and overall BC, and nine loci shared between depression and ER + BC (Table 3). No significant shared locus was found for depression with ER– BC.

Table 3 Genome-wide significant loci shared between depression and breast cancer in cross-trait meta-analysis (PCPASSOC < 5 × 10–8 and Psingle-trait < 1 × 10–3)

Of these 17 pleiotropic SNPs, interestingly, SNP rs2403907 represented the strongest shared signal across both overall (PCPASSOC = 1.55 × 10–33) and ER + BC (PCPASSOC = 1.91 × 10–32), located at an intergenic region at 21q21.1. SNPs rs67981811 (overall BC: PCPASSOC = 6.25 × 10–24; ER + BC: PCPASSOC = 3.14 × 10–23) and rs17693963 (overall BC: PCPASSOC = 8.96 × 10–21; ER + BC: PCPASSOC = 8.85 × 10–21) were the second and the third strongest shared signals, both located at the 6p22.1 region that harbors a depression-associated risk locus [9].

In addition to “known” pleiotropic SNPs (SNPs significantly associated with both traits in previous GWASs, N = 6) or “single-trait-driven” pleiotropic SNPs (SNPs significantly associated with one of the two traits in previous GWASs, N = 10), we identified one novel shared SNP (re56101042) for depression and overall BC (PCPASSOC = 6.53 × 10–9). This SNP is located at 14q32.32 near snoU13 and TRAF3. The former gene represents a small nucleolar RNA recently implicated as a biomarker in vulnerability to brain injury [40], while the latter gene functions in immune modulation and inflammation via its regulation of nuclear factor-κB signaling [41], which has been implicated in the pathogenesis of both psychiatric disorders and tumorigenesis [1]. Four genes (AMN, ANKRD9, MIR4309, and RCOR1) were additionally identified to interact with rs56101042 through three-dimensional chromatin looping [25]. Detailed annotations for each SNP discovered by CPASSOC are shown in Additional file 3: Tables S4-S5.

Identification of causal variants and colocalization

For each of the 17 pleiotropic SNPs, we further determined a 99% credible set of causal SNPs, providing targets for downstream experimental analysis (Additional file 3: Tables S6-S7). We found 413 candidate causal SNPs across all loci shared by depression and overall BC, and 282 candidate causal SNPs across all loci shared by depression and ER + BC. In particular, SNPs rs5757946 (22q13.1) and rs11049366 (12p11.22) were identified as having a posterior probability of 1.00 in the 99% credible set.

We also performed colocalization analysis to determine whether pleiotropic SNPs driving the associations in two traits were the same (Additional file 3: Table S8). A substantial proportion of shared loci (five out of the 12 loci shared by depression and overall BC; and five out of the nine loci shared by depression and ER + BC) colocalized at the same candidate SNPs (eight of which showed a PPH4 > 0.90), reinforcing shared causal associations.

Transcriptome-wide association studies and shared genes

Results from tissue-specific TWAS revealed gene-level genetic overlap (single-trait results shown in Additional file 3: Table S9 and Additional file 1: Figure S2). After multiple corrections, a total of five genes were shared by depression and overall BC of which two were also shared by depression and ER + BC, enriched in six tissues including brain, blood, uterus, colon, artery and heart (Table 4). Four genes were located at pleiotropic loci identified in cross-trait meta-analysis, including FLOT1 and HLA-S at 6p21.33, ENSG00000247934.4 at 12p11.22, and GRAP2 at 22q13.1.

Table 4 Significant shared transcriptome-wide association analysis results between depression and breast cancer

While data on ENSG00000247934.4 was limited, it represents an antisense to CCDC91 previously implicated in BC and brain morphology estimates (compared with data from GWAS Catalog [13]). Among other TWAS-significant shared genes, FLOT1 was previously implicated in depression and aerodigestive squamous cell carcinoma. We also identified three novel genes, HLA-S, ZSWIM1, and GRAP2. As a pseudogene located at the major histocompatibility complex (MHC) region, HLA-S was reported to be associated with a broad range of conditions including schizophrenia, gastric, and hepatocellular cancers. By regulating T helper cells development, ZSWIM1 plays a critical role in immune system [42], and signals a poor prognosis of ovarian cancer [43]. Similarly, GRAP2 was implicated in T cell-mediated immune responses [44], as well as in medullary thyroid cancer and lung adenocarcinoma prognosis [45, 46].

Mendelian randomization

Finally, we conducted a two-sample MR analysis combining 102 GWAS-identified depression-associated SNPs (F-statistic ranging from 91.5 to 329.6 (Additional file 3: Table S1)) as an IV to examine a potential causal effect of depression on BC. Using IVW, genetic liability to depression was significantly associated with an increased risk of overall BC (odds ratio (OR) = 1.12, 95%CIs = 1.04–1.19, P = 1.40 × 10–3). The estimates remained directionally consistent in MR-Egger regression and the weighted median approach, despite larger statistical uncertainties (Fig. 3). No indication of horizontal pleiotropy was observed (PMR-Egger intercept = 0.86). Sensitivity analyses excluding palindromic SNPs or pleiotropic SNPs supported the robustness of these results. Leave-one-out analysis demonstrated the absence of outlying variants (Additional file 1: Figure S3) which was validated by MR-PRESSO (OR = 1.13, 95%CIs = 1.06–1.21, P = 5.00 × 10–4). MVMR accounting for potential confounders generated similar results with even more pronounced magnitude and significance, suggesting a depression-BC causal relationship independent of common confounders (Additional file 1: Figure S4).

Fig. 3
figure 3

Bi-directional causal relationship underlying depression and breast cancer. a Estimates of causal effects for genetically predicted depression with overall breast cancer, ER + breast cancer, and ER– breast cancer, respectively. b Estimates of causal effects for genetically predicted overall breast cancer, genetically predicted ER + breast cancer, and genetically predicted ER– breast cancer with depression. Boxes represent the point estimates of causal effects, and error bars represent 95% confidence intervals. Inverse-variance weighted approach was adopted as the primary analysis. MR-Egger regression, weighted median, and MR-PRESSO approaches were adopted as sensitivity analysis. BC: breast cancer; ER: estrogen receptor

In subgroup analysis defined by ER status, we found a non-significant association of depression with ER + BC (IVW OR = 1.08, 95%CIs = 0.99–1.18, P = 6.98 × 10–2; PMR-Egger intercept = 0.86), and a suggestive significant association of depression with ER– BC (IVW OR = 1.12, 95%CIs = 1.01–1.24, P = 3.77 × 10–2; PMR-Egger intercept = 0.82), though neither withstood multiple testing corrections (Fig. 3). In MVMR, while a putative causal effect of depression on ER + BC was observed, the causal effect of depression with ER– BC attenuated to null after adjustment for each confounder, and none survived multiple corrections (Additional file 1: Figure S4).

No evidence of reverse causality was found such that genetically predicted BC did not seem to influence depression risk (IVW OR = 1.01, P = 0.49; MR-Egger OR = 0.98, P = 0.20; weighted median OR = 1.01, P = 0.34). Similar null effects were identified for ER + and ER– subtypes (Fig. 3).

Discussion

To our knowledge, this is the largest observational and genetic analysis to systematically investigate the phenotypic association as well as the shared genetic architecture underlying depression and BC. Our observational study demonstrated a positive association of depression with BC risk, which was corroborated by findings from genetic analyses. We found evidence supporting a significant shared genetic basis, both globally and regionally, indicating shared biology between the two traits. This genetic overlap was further decomposed into both pleiotropy and causality, reflected by the pleiotropic loci identified in CPASSOC, the shared genes revealed by TWAS, and the putative causal relationship demonstrated by MR.

Overall, the genetic correlation between depression and BC was limited (\({r}_{g}\) = 0.08) but statistically significant. In addition, significant local genetic correlations were found for several specific genomic regions including a locus at 9q31.2, and 4 loci at, or close to 6p22.1. Chr6q22.1 represents an important locus at MHC, a region well-known for containing multiple genes essential for the adaptive immune system as well as the anti-tumoral immune response [47]. The MHC region has been associated with depression in some prior studies [9,10,11, 48], with the strongest association located in the classical or extended class I region; moreover, a study has also identified a 66% down-regulation of MHC class I expression in BC tissues [49], collectively highlighting its role in the etiological connection underlying these two traits. Leveraging genotype imputation, a recent (also the largest) study on depression and classical MHC variants indicated an overall non-significant association [50], calling for additional efforts to dissect signals within the extended MHC.

Significant genetic overlap can be the result of biological or horizontal pleiotropy (in which a genetic variant affects multiple traits) and/or mediated or vertical pleiotropy (in which a genetic variant affects a trait via its effect on an intermediate trait) [51]. In our downstream analysis performed to explore these alternatives, 17 pleiotropic loci between depression and BC were identified, harboring genes which were previously implicated in neural development and brain functions (i.e., HCN1, ABT1, ZSCAN12, ZNF184, KLC1, MKL1, SGSM3, CCDC91), or biological processes related to tumor growth (i.e., ABT1, ZSCAN12, ZNF184, KLC1, TTC28, MKL1, CCDC91, OLFM4, DDR1). Notably, several shared loci were located at genes that play essential roles in immune, inflammatory, and stress responses (i.e., PSORS1C1, TRAF3, APOPT1, PRSS16, BTN3A2). Multiple genes showed strong evidence of colocalization, especially among those located near the MHC (i.e., ABT1, BTN3A2, ZNF184, PRSS16, ZSCAN12, SLC17A3).

One advantage of meta-analyzing GWASs of different traits is that combining association evidence across multiple studies can reveal signals which might not have reached genome-wide significance in a single-trait analysis. Indeed, we found a novel locus (lead SNP rs56101042) shared between depression and BC, only reported previously as associated with allergic disease [52]. While its direct involvement in the development of depression or BC is not known, this SNP interacts with several genes that have already been reported (i.e., TRAF3 and RCOR1 [53, 54]), or that could play a role in the pathogenesis of the two traits (i.e., AMN, ANKRD9, and MIR4309). Encoding a transmembrane protein named amnionless, AMN is involved in the plasma membrane transport of Cubilin, and may influence cancer progression through the Cubilin-FGF8 (fibroblast growth factor 8) interaction [55]. AMN may also function in psychiatric illness given its role in the vitamin B12 transportation into nervous system [56]. ANKRD9 encodes an ankyrin repeat domain protein which was implicated in drug-induced mental disorder [57], and recognized as a tumor suppressor in gastric cancer among Asians [58]. MIR4309 has been found to be significantly elevated among patients with disordered sleep (a core symptom of depression) [59], as well as up-regulated in tissues of gastric cancer [60]. Further experimental studies would be needed to provide more detailed functional annotation of rs56101042, particularly in relation to the onset of depression and BC.

By integrating data from GWAS and GTEx tissue expression, TWAS suggested shared mechanistic hypotheses between depression and BC on a gene-tissue pair level. The four loci identified in both CPASSOC and TWAS analysis implicate common biological mechanisms in depression and BC regulation, involving cell proliferation [61], brain structure [62, 63], and immune response [44]. In addition to tissues of brain, uterus, and blood (well-recognized as relevant to depression and/or BC [64,65,66]), our TWAS reported shared regulatory features in the digestive (i.e., colon transverse) and cardiovascular systems (i.e., artery coronary and hear left ventricle), suggesting the possibility of shared pathways extending to a wider range of organs. A growing body of literature has shown that functional disorders of the digestive system may play a critical role in a variety of psychiatric disorders including depression. One of the speculated mechanisms is the communication of gut microbiome with the brain through the brain-gut axis [67]. Interestingly, gut microbiome may also affect breast carcinogenesis by promoting anti-tumor immunity and/or modulating systemic estrogen levels [68]. As another major cause of mortality and morbidity, cardiovascular disease shares a number of common risk factors with both conditions, and represents one of the most prevalent comorbidities for depression and BC in clinical settings [69, 70]. More studies are needed to clarify the role of these hypothesized mechanisms.

To date, epidemiological evidence regarding the depression-BC relationship has been inconclusive, but our observational and MR analyses showed consistent evidence for a putative causal role of depression in BC among European populations. Our MR results extended existing MR analyses [71, 72] with substantially improved accuracy and precision of causal estimates by applying: (i) the most updated GWAS of overall BC; (ii) a broader range of complementary analysis to validate the MR assumptions; and (iii) an extensive MVMR design to evaluate the independent effect of depression. The estimated causal effects remained directionally consistent across a variety of sensitivity analyses, as well as after adjustment of multiple confounders, supporting the robustness of our MR findings. When extended to BC subtypes, we found consistently positive associations with greater statistical uncertainties. Future studies with larger sample sizes are warranted to establish or rule out a potential causal link in a subtype-specific manner.

Several limitations need to be acknowledged. Using data of European ancestry populations restricted the generalizability of our results, and future studies including more diverse ancestries are needed. The broad depression phenotype used in the current study is regarded as a relatively non-specific index of vulnerability to psychological distress, rather than highly specific for MD [73]. However, we performed sensitivity analysis restricting to clinically ascertained MD cases and identified consistent results (fully adjusted observational analysis using only ICD-10 codes F32 and F33 for MD diagnosis: HR = 1.091, 95%CIs = 0.946–1.258; genetic correlation and MR analyses applying an older MD GWAS [10] involving only clinically defined MD cases: \({r}_{g}\) = 0.07, P = 3.00 × 10–4; ORMR = 1.08, P = 4.30 × 10–2). Additionally, due to limited data availability at the time of conducting the analysis, we were unable to utilize sex-specific GWAS data of depression to match with the female-specific outcome BC. However, sex-heterogeneity did not seem to play a significant role when using female-specific depression GWAS summary statistics which now become available [74]. Although underpowered, female-specific genetic correlation analysis and MR analysis yielded statistically significant and directionally consistent results (\({r}_{g}\) = 0.09, P = 1.60 × 10−3; ORIVW = 1.11, P = 0.04). Using sufficiently powered female-specific data would be valuable for future studies. Although depression is highly heterogeneous, we did not investigate depression subtypes (i.e., typical vs. atypical) or additional BC subtypes due to limited well-powered GWAS available for these phenotypes. Finally, validation of the study findings in additional follow-up cohorts and GWAS is needed. Future in-depth experimental works are also warranted to better understand pathophysiological mechanisms that may underlie genetic overlap between depression and BC.

Conclusions

In sum, the current study advances our understanding of previously reported associations between depression and BC by providing evidence of genetic correlation, pleiotropic loci, and a possible causal effect of depression on BC. Our findings further suggest potential biological mechanisms linking these disorders and may have implications for future studies aimed reducing BC risk.

Availability of data and materials

Data from UK Biobank is available for open access to scientific researcher (www.biobank.ac.uk). Depression GWAS summary statistics (without 23andMe) is publicly available from https://datashare.ed.ac.uk/handle/10283/3203. GWAS summary statistics for overall and subtype-specific BCs are publicly available from https://bcac.ccge.medschl.cam.ac.uk/bcacdata/oncoarray/oncoarray-and-combined-summary-result/.

Abbreviations

BC:

Breast cancer

BMI:

Body mass index

CIs:

Confidence intervals

CPASSOC:

Cross-Phenotype Association

ER:

Estrogen receptor

GTEx:

Genotype-Tissue Expression

GWAS:

Genome-wide association study

GWAS-PW:

Pairwise-GWAS

HESS:

Heritability Estimator from Summary Statistics

HR:

Hazard ratio

IV:

Instrumental variable

IVW:

Inverse-variance weighted

LD:

Linkage disequilibrium

LDSC:

Linkage disequilibrium Score Regression

MD:

Major depression

MHC:

Major histocompatibility complex

MR:

Mendelian randomization

MR-PRESSO:

MR-Pleiotropy Residual Sum and Outlier

MVMR:

Multi-variable Mendelian randomization

OR:

Odds ratio

PPA:

Posterior probabilities of association

r g :

Genetic correlation

RR:

Relative risk

SNP:

Single nucleotide polymorphism

TWAS:

Transcriptome-wide association study

UKBB:

UK Biobank

VEP:

Ensembl Variant Effect Predictor

References:

  1. Bortolato B, Hyphantis TN, Valpione S, et al. Depression in cancer: The many biobehavioral pathways driving tumor progression. Cancer Treat Rev. 2017;52:58–70. https://doi.org/10.1016/J.CTRV.2016.11.004.

    Article  PubMed  Google Scholar 

  2. James SL, Abate D, Abate KH, et al. Global, regional, and national incidence, prevalence, and years lived with disability for 354 Diseases and Injuries for 195 countries and territories, 1990–2017: A systematic analysis for the Global Burden of Disease Study 2017. Lancet. 2018;392(10159):1789–858. https://doi.org/10.1016/S0140-6736(18)32279-7/ATTACHMENT/DB6E3413-74DC-43AE-B7CC-CA155C28589E/MMC2.PDF.

    Article  Google Scholar 

  3. Sung H, Ferlay J, Siegel RL, et al. Global Cancer Statistics 2020: GLOBOCAN Estimates of Incidence and Mortality Worldwide for 36 Cancers in 185 Countries. CA Cancer J Clin. 2021;71(3):209–49. https://doi.org/10.3322/CAAC.21660.

    Article  PubMed  Google Scholar 

  4. Sun HL, Dong XX, Cong YJ, et al. Depression and the risk of breast cancer: a meta-analysis of cohort studies. Asian Pac J Cancer Prev. 2015;16(8):3233–9. https://doi.org/10.7314/APJCP.2015.16.8.3233.

    Article  PubMed  Google Scholar 

  5. Oerlemans MEJ, van den Akker M, Schuurman AG, Kellen E, Buntinx F. A meta-analysis on depression and subsequent cancer risk. Clin Pract Epidemiol Ment Health. 2007;3:29. https://doi.org/10.1186/1745-0179-3-29.

    Article  PubMed  PubMed Central  Google Scholar 

  6. Zhu Z, Hasegawa K, Camargo CA, Liang L. Investigating asthma heterogeneity through shared and distinct genetics: Insights from genome-wide cross-trait analysis. J Allergy Clin Immunol. 2021;147(3):796–807. https://doi.org/10.1016/j.jaci.2020.07.004.

    Article  CAS  PubMed  Google Scholar 

  7. Sullivan PF, Neale MC, Kendler KS. Genetic epidemiology of major depression: review and meta-analysis. Am J Psychiatry. 2000;157(10):1552–62. https://doi.org/10.1176/APPI.AJP.157.10.1552.

    Article  CAS  PubMed  Google Scholar 

  8. Möller S, Mucci LA, Harris JR, et al. The heritability of breast cancer among women in the nordic twin study of cancer. Cancer Epidemiol Biomarkers Prev. 2016;25(1):145–50. https://doi.org/10.1158/1055-9965.EPI-15-0913/68674/AM/THE-HERITABILITY-OF-BREAST-CANCER-AMONG-WOMEN-IN.

    Article  PubMed  Google Scholar 

  9. Howard DM, Adams MJ, Clarke TK, et al. Genome-wide meta-analysis of depression identifies 102 independentvariants and highlights the importance of the prefrontal brainregions. Nat Neurosci. 2019;22(3):343. https://doi.org/10.1038/S41593-018-0326-7.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  10. Wray NR, Ripke S, Mattheisen M, et al. Genome-wide association analyses identify 44 risk variants and refine the genetic architecture of major depression. Nat Genet. 2018;50(5):668–81. https://doi.org/10.1038/S41588-018-0090-3.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  11. Howard DM, Adams MJ, Shirali M, et al. Genome-wide association study of depression phenotypes in UK Biobank identifies variants in excitatory synaptic pathways. Nat Commun. 2018;9(1). https://doi.org/10.1038/S41467-018-03819-3

  12. Shah R, Rosso K, Nathanson SD. Pathogenesis, prevention, diagnosis and treatment of breast cancer. World J Clin Oncol. 2014;5(3):283. https://doi.org/10.5306/WJCO.V5.I3.283.

    Article  PubMed  PubMed Central  Google Scholar 

  13. MacArthur J, Bowler E, Cerezo M, et al. The new NHGRI-EBI Catalog of published genome-wide association studies (GWAS Catalog). Nucleic Acids Res. 2017;45(Database issue):D896. https://doi.org/10.1093/NAR/GKW1133.

    Article  CAS  PubMed  Google Scholar 

  14. Hwang WJ, Lee TY, Kim NS, Kwon JS. The Role of Estrogen Receptors and Their Signaling across Psychiatric Disorders. Int J Mol Sci. 2021;22(1):373. https://doi.org/10.3390/IJMS22010373.

    Article  CAS  Google Scholar 

  15. Sudlow C, Gallacher J, Allen N, et al. UK biobank: an open access resource for identifying the causes of a wide range of complex diseases of middle and old age. PLoS Med. 2015;12(3). https://doi.org/10.1371/JOURNAL.PMED.1001779

  16. Lawlor DA, Harbord RM, Sterne JAC, Timpson N, Davey Smith G. Mendelian randomization: using genes as instruments for making causal inferences in epidemiology. Stat Med. 2008;27(8):1133–63. https://doi.org/10.1002/SIM.3034.

    Article  PubMed  Google Scholar 

  17. Zhang H, Ahearn TU, Lecarpentier J, et al. Genome-wide association study identifies 32 novel breast cancer susceptibility loci from overall and subtype-specific analyses. Nat Genet. 2020;52(6):572–81. https://doi.org/10.1038/s41588-020-0609-2.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  18. Michailidou K, Lindström S, Dennis J, et al. Association analysis identifies 65 new breast cancer risk loci. Nat. 2017;551(7678):92–4. https://doi.org/10.1038/nature24284.

    Article  CAS  Google Scholar 

  19. Bulik-Sullivan B, Finucane HK, Anttila V, et al. An atlas of genetic correlations across human diseases and traits. Nat Genet. 2015;47(11):1236–41. https://doi.org/10.1038/ng.3406.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  20. Shi H, Mancuso N, Spendlove S, Pasaniuc B. Local Genetic Correlation Gives Insights into the Shared Genetic Architecture of Complex Traits. Am J Hum Genet. 2017;101(5):737–51. https://doi.org/10.1016/J.AJHG.2017.09.022.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  21. Pickrell JK, Berisa T, Liu JZ, Ségurel L, Tung JY, Hinds DA. Detection and interpretation of shared genetic influences on 42 human traits. Nat Genet. 2016;48(7):709. https://doi.org/10.1038/NG.3570.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  22. Zhu X, Feng T, Tayo BO, et al. Meta-analysis of correlated traits via summary statistics from GWASs with an application in hypertension. Am J Hum Genet. 2015;96(1):21–36. https://doi.org/10.1016/J.AJHG.2014.11.011.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  23. Purcell S, Neale B, Todd-Brown K, et al. PLINK: a tool set for whole-genome association and population-based linkage analyses. Am J Hum Genet. 2007;81(3):559–75. https://doi.org/10.1086/519795.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  24. Zerbino D, Achuthan P, Akanni W, et al. Ensembl 2018. Nucleic Acids Res. 2018;46(D1):D754–61. https://doi.org/10.1093/NAR/GKX1098.

    Article  CAS  PubMed  Google Scholar 

  25. Lu Y, Quan C, Chen H, Bo X, Zhang C. 3DSNP: a database for linking human noncoding SNPs to their three-dimensional interacting genes. Nucleic Acids Res. 2017;45(D1):D643–9. https://doi.org/10.1093/NAR/GKW1022.

    Article  CAS  PubMed  Google Scholar 

  26. Farh KK-H, Marson A, Zhu J, et al. Genetic and epigenetic fine mapping of causal autoimmune disease variants. Nat. 2014;518(7539):337–43. https://doi.org/10.1038/nature13835.

    Article  CAS  Google Scholar 

  27. Giambartolomei C, Vukcevic D, Schadt EE, et al. Bayesian Test for Colocalisation between Pairs of Genetic Association Studies Using Summary Statistics. PLOS Genet. 2014;10(5):e1004383. https://doi.org/10.1371/JOURNAL.PGEN.1004383.

    Article  PubMed  PubMed Central  Google Scholar 

  28. Gusev A, Ko A, Shi H, et al. Integrative approaches for large-scale transcriptome-wide association studies. Nat Genet. 2016;48(3):245–52. https://doi.org/10.1038/NG.3506.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  29. Burgess S, Scott RA, Timpson NJ, Smith GD, Thompson SG. Using published data in Mendelian randomization: a blueprint for efficient identification of causal risk factors. Eur J Epidemiol. 2015;30(7):543–52. https://doi.org/10.1007/S10654-015-0011-Z.

    Article  PubMed  PubMed Central  Google Scholar 

  30. Bowden J, Smith GD, Burgess S. Mendelian randomization with invalid instruments: effect estimation and bias detection through Egger regression. Int J Epidemiol. 2015;44(2):512–25. https://doi.org/10.1093/IJE/DYV080.

    Article  PubMed  PubMed Central  Google Scholar 

  31. Bowden J, Smith GD, Haycock PC, Burgess S. Consistent Estimation in Mendelian Randomization with Some Invalid Instruments Using a Weighted Median Estimator. Genet Epidemiol. 2016;40(4):304–14. https://doi.org/10.1002/GEPI.21965.

    Article  PubMed  PubMed Central  Google Scholar 

  32. Davey Smith G, Ebrahim S. ‘Mendelian randomization’: can genetic epidemiology contribute to understanding environmental determinants of disease? Int J Epidemiol. 2003;32(1):1–22. https://doi.org/10.1093/IJE/DYG070.

    Article  Google Scholar 

  33. Verbanck M, Chen C-Y, Neale B, Do R. Detection of widespread horizontal pleiotropy in causal relationships inferred from Mendelian randomization between complex traits and diseases. Nat Genet. 2018;50(5):693–8. https://doi.org/10.1038/s41588-018-0099-7.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  34. Burgess S, Thompson SG. Multivariable Mendelian Randomization: The Use of Pleiotropic Genetic Variants to Estimate Causal Effects. Am J Epidemiol. 2015;181(4):251–60. https://doi.org/10.1093/AJE/KWU283.

    Article  PubMed  PubMed Central  Google Scholar 

  35. Pulit SL, Stoneman C, Morris AP, et al. Meta-analysis of genome-wide association studies for body fat distribution in 694 649 individuals of European ancestry. Hum Mol Genet. 2019;28(1):166–74. https://doi.org/10.1093/HMG/DDY327.

    Article  CAS  PubMed  Google Scholar 

  36. Liu M, Jiang Y, Wedow R, et al. Association studies of up to 1.2 million individuals yield new insights into the genetic etiology of tobacco and alcohol use. Nat Genet. 2019;51(2):237–44. https://doi.org/10.1038/S41588-018-0307-5.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  37. Klimentidis YC, Raichlen DA, Bea J, et al. Genome-wide association study of habitual physical activity in over 377,000 UK Biobank participants identifies multiple variants including CADM2 and APOE. Int J Obes. 2018;42(6):1161–76. https://doi.org/10.1038/s41366-018-0120-3.

    Article  CAS  Google Scholar 

  38. Dashti HS, Jones SE, Wood AR, et al. Genome-wide association study identifies genetic loci for self-reported habitual sleep duration supported by accelerometer-derived estimates. Nat Commun. 2019;10(1):1–12. https://doi.org/10.1038/s41467-019-08917-4.

    Article  CAS  Google Scholar 

  39. Lee JJ, Wedow R, Okbay A, et al. Gene discovery and polygenic prediction from a genome-wide association study of educational attainment in 1.1 million individuals. Nat Genet. 2018;50(8):1112–21. https://doi.org/10.1038/s41588-018-0147-3.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  40. Di Pietro V, O’Halloran P, Watson CN, et al. Unique diagnostic signatures of concussion in the saliva of male athletes: the Study of Concussion in Rugby Union through MicroRNAs (SCRUM). Br J Sports Med. 2021;55(24):1395. https://doi.org/10.1136/BJSPORTS-2020-103274.

    Article  PubMed  Google Scholar 

  41. Häcker H, Tseng PH, Karin M. Expanding TRAF function: TRAF3 as a tri-faced immune regulator. Nat Rev Immunol. 2011;11(7):457–68. https://doi.org/10.1038/nri2998.

    Article  CAS  PubMed  Google Scholar 

  42. Ko KK, Powell MS, Hogarth PM. ZSWIM1: a novel biomarker in T helper cell differentiation. Immunol Lett. 2014;160(2):133–8. https://doi.org/10.1016/J.IMLET.2014.01.016.

    Article  CAS  PubMed  Google Scholar 

  43. Pawar A, Chowdhury OR, Chauhan R, Talole S, Bhattacharjee A. Identification of key gene signatures for the overall survival of ovarian cancer. J Ovarian Res. 2022;15(1).https://doi.org/10.1186/S13048-022-00942-0

  44. Yablonski D. Bridging the Gap: Modulatory Roles of the Grb2-Family Adaptor, Gads, in Cellular and Allergic Immune Responses. Front Immunol. 2019;10:1704. https://doi.org/10.3389/FIMMU.2019.01704/FULL.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  45. Ludwig L, Oswald F, Hoang-Vu C, et al. Expression of the Grb2-related RET adapter protein Grap-2 in human medullary thyroid carcinoma. Cancer Lett. 2009;275(2):194–7. https://doi.org/10.1016/J.CANLET.2008.10.010.

    Article  CAS  PubMed  Google Scholar 

  46. S S, X D, J W, et al. GRAP2 is a Prognostic Biomarker and Correlated with Immune Infiltration in Lung Adenocarcinoma. Published online 10 Jan, 2022. https://doi.org/10.21203/RS.3.RS-1227511/V1

  47. Rock KL, Reits E, Neefjes J. Present Yourself! By MHC Class I and MHC Class II Molecules. Trends Immunol. 2016;37(11):724–37. https://doi.org/10.1016/J.IT.2016.08.010.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  48. Li H, Chang H, Song X, et al. Integrative analyses of major histocompatibility complex loci in the genome-wide association studies of major depressive disorder. Neuropsychopharmacology. 2019;44(9):1552. https://doi.org/10.1038/S41386-019-0346-3.

    Article  PubMed  PubMed Central  Google Scholar 

  49. Kaneko K, Ishigami S, Kijima Y, et al. Clinical implication of HLA class I expression in breast cancer. BMC Cancer. 2011;11:454. https://doi.org/10.1186/1471-2407-11-454.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  50. Glanville KP, Coleman JRI, Hanscombe KB, et al. Classical Human Leukocyte Antigen Alleles and C4 Haplotypes Are Not Significantly Associated With Depression. Biol Psychiatry. 2020;87(5):419–30. https://doi.org/10.1016/J.BIOPSYCH.2019.06.031.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  51. Solovieff N, Cotsapas C, Lee PH, Purcell SM, Smoller JW. Pleiotropy in complex traits: challenges and strategies. Nat Rev Genet. 2013;14(7):483–95. https://doi.org/10.1038/NRG3461.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  52. Ferreira MAR, Vonk JM, Baurecht H, et al. Age-of-onset information helps identify 76 genetic variants associated with allergic disease. PLOS Genet. 2020;16(6):e1008725. https://doi.org/10.1371/JOURNAL.PGEN.1008725.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  53. Lv H, Pan G, Zheng G, et al. Expression and functions of the repressor element 1 (RE-1)-silencing transcription factor (REST) in breast cancer. J Cell Biochem. 2010;110(4):968–74. https://doi.org/10.1002/JCB.22610.

    Article  CAS  PubMed  Google Scholar 

  54. Monaghan CE, Nechiporuk T, Jeng S, et al. REST corepressors RCOR1 and RCOR2 and the repressor INSM1 regulate the proliferation-differentiation balance in the developing brain. Proc Natl Acad Sci U S A. 2017;114(3):E406–15. https://doi.org/10.1073/PNAS.1620230114.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  55. Kozyraki R, Cases O. Cubilin, the Intrinsic Factor-Vitamin B12 Receptor in Development and Disease. Curr Med Chem. 2020;27(19):3123–50. https://doi.org/10.2174/0929867325666181008143945.

    Article  CAS  PubMed  Google Scholar 

  56. Luder AS, Tanner SM, de la Chapelle A, Walter JH. Amnionless (AMN) mutations in Imerslund-Gräsbeck syndrome may be associated with disturbed vitamin B12 transport into the CNS. J Inherit Metab Dis. 2008;31:493–6. https://doi.org/10.1007/S10545-007-0760-2.

    Article  PubMed  Google Scholar 

  57. de la Fuente RM, Ibi D, Saunders JM, et al. HDAC2-dependent Antipsychotic-like Effects of Chronic Treatment with the HDAC Inhibitor SAHA in Mice. Neuroscience. 2018;388:102–17. https://doi.org/10.1016/J.NEUROSCIENCE.2018.07.010.

    Article  Google Scholar 

  58. Lee Y, Lim B, Lee SW, et al. ANKRD9 is associated with tumor suppression as a substrate receptor subunit of ubiquitin ligase. Biochim Biophys Acta - Mol Basis Dis. 2018;1864(10):3145–53. https://doi.org/10.1016/J.BBADIS.2018.07.001.

    Article  CAS  PubMed  Google Scholar 

  59. Liu Z, Yang L, Zhao Y, et al. Reproducibility of quantitative real-time PCR assay in microRNA expression profiling and comparison with microarray analysis in narcolepsy. Springerplus. 2015;4(1):1–6. https://doi.org/10.1186/S40064-015-1613-3/FIGURES/3.

    Article  Google Scholar 

  60. Liu D, Hu X, Zhou H, Shi G, Wu J. Identification of aberrantly expressed miRNAs in gastric cancer. Gastroenterol Res Pract. 2014. https://doi.org/10.1155/2014/473817

  61. Lin C, Wu Z, Lin X, et al. Knockdown of FLOT1 impairs cell proliferation and tumorigenicity in breast cancer through upregulation of FOXO3a. Clin Cancer Res. 2011;17(10):3089–99. https://doi.org/10.1158/1078-0432.CCR-10-3068/84191/AM/KNOCKDOWN-OF-FLOT1-IMPAIRS-CELL-PROLIFERATION-AND.

    Article  CAS  PubMed  Google Scholar 

  62. Sprooten E, Knowles EE, McKay DR, et al. Common genetic variants and gene expression associated with white matter microstructure in the human brain. Neuroimage. 2014;97:252. https://doi.org/10.1016/J.NEUROIMAGE.2014.04.021.

    Article  CAS  PubMed  Google Scholar 

  63. Van Der Meer D, Kaufmann T, Shadrin AA, et al. The genetic architecture of human cortical folding. Sci Adv. 2021;7(51):9446. https://doi.org/10.1126/SCIADV.ABJ9446/SUPPL_FILE/SCIADV.ABJ9446_DATA_S1_TO_S9.ZIP.

    Article  Google Scholar 

  64. Phillips JR, Hewedi DH, Eissa AM, Moustafa AA. The Cerebellum and Psychiatric Disorders. Front Public Heal. 2015;3:66. https://doi.org/10.3389/FPUBH.2015.00066.

    Article  Google Scholar 

  65. Shen T-C, Hsia TC, Hsiao CL, et al. Patients with uterine leiomyoma exhibit a high incidence but low mortality rate for breast cancer. Oncotarget. 2017;8(20):33014–23. https://doi.org/10.18632/ONCOTARGET.16520.

    Article  PubMed  PubMed Central  Google Scholar 

  66. Chaussabel D, Pascual V, Banchereau J. Assessing the human immune system through blood transcriptomics. BMC Biol. 2010;8(1):1–14. https://doi.org/10.1186/1741-7007-8-84/FIGURES/3.

    Article  Google Scholar 

  67. Bastiaanssen TFS, Cussotto S, Claesson MJ, Clarke G, Dinan TG, Cryan JF. Gutted! Unraveling the Role of the Microbiome in Major Depressive Disorder. Harv Rev Psychiatry. 2020;28(1):26–39. https://doi.org/10.1097/HRP.0000000000000243.

    Article  PubMed  PubMed Central  Google Scholar 

  68. Fernández MF, Reina-Pérez I, Astorga JM, Rodríguez-Carrillo A, Plaza-Díaz J, Fontana L. Breast Cancer and Its Relationship with the Microbiota. Int J Environ Res Public Health. 2018;15(8). https://doi.org/10.3390/IJERPH15081747

  69. Ardekani AM, Naeini MM. The Role of MicroRNAs in Human Diseases. Avicenna J Med Biotechnol. 2010;2(4):161 (Accessed March 1, 2022. /pmc/articles/PMC3558168/).

    CAS  PubMed  PubMed Central  Google Scholar 

  70. Hare DL, Toukhsati SR, Johansson P, Jaarsma T. Depression and cardiovascular disease: a clinical review. Eur Heart J. 2014;35(21):1365–72. https://doi.org/10.1093/EURHEARTJ/EHT462.

    Article  PubMed  Google Scholar 

  71. Zhu G-L, Xu C, Yang K, et al. Causal relationship between genetically predicted depression and cancer risk: a two-sample bi-directional mendelian randomization. BMC Cancer. 2022;22(1):1–10. https://doi.org/10.1186/S12885-022-09457-9.

    Article  Google Scholar 

  72. Ren Q, Luo F, Ge S, Chen P. Major depression disorder may causally associate with the increased breast cancer risk: Evidence from two-sample mendelian randomization analyses. Cancer Med. 2022;00:1–13. https://doi.org/10.1002/CAM4.5043.

    Article  Google Scholar 

  73. Cai N, Revez JA, Adams MJ, et al. Minimal phenotyping yields genome-wide association signals of low specificity for major depression. Nat Genet. 2020;52(4):437–47. https://doi.org/10.1038/s41588-020-0594-5.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  74. Silveira PP, Pokhvisneva I, Howard DM, Meaney MJ. A sex-specific genome-wide association study of depression phenotypes in UK Biobank. Mol Psychiatry 2023. Published online February 7, 2023:1–11. https://doi.org/10.1038/s41380-023-01960-0

Download references

Acknowledgements

Not applicable.

Funding

Open access funding provided by Karolinska Institute. This study was supported by funds from the National Key R&D Program of China (2022YFC3600600, 2022YFC3600604), the National Natural Science Foundation of China (82204170, U22A20359, 81874283, 81673255), the Recruitment Program for Young Professionals of China, the Promotion Plan for Basic Medical Sciences and the Development Plan for Cutting-Edge Disciplines, Sichuan University, and other Projects from West China School of Public Health and West China Fourth Hospital, Sichuan University. The funders of the study had no role in study design, data collection, data analysis, data interpretation, or writing of the report.

Author information

Authors and Affiliations

Authors

Contributions

X.J. and X.W. conceived and designed the study. X.W., W.Z., X.Z., L.Z., M.X., Y.H., and J.X contributed to analysis of the data. X.W. and X.J. drafted the article. J.W.S., P.K., J.L., and B.Z. contributed to the interpretations of the findings and the critical revision of the article. X.J. supervised the study. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Xia Jiang.

Ethics declarations

Ethics approval and consent to participate

UK Biobank had obtained ethics approval from the North West Multi-Centre Research Ethics Committee which covers the UK (approval number: 11/NW/0382) and had obtained informed consent from all participants. The current study was approved by the UK Biobank access management board. Participants who had withdrawn consent by the time of the analysis were excluded from the analysis dataset. The GWAS summary statistics used in the present study are aggregated level of data which do not contain any personal identifiers. The original GWAS have obtained ethical approval from relevant ethics review committees.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher’s Note

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

Supplementary Information

12916_2023_2876_MOESM1_ESM.pdf

Additional file 1: Figure S1. Graphical abstract. Figure S2. Number of TWAS significant genes for depression, overall BC, and subtype-specific BC. Figure S3. Box plot of betas in leave-one-out analysis. Figure S4. Results of multi-variable Mendelian randomization analysis.

Additional file 2. STROBE and STROBE-MR checklists.

12916_2023_2876_MOESM3_ESM.xlsx

Additional file 3: Table S1. Depression IVs. Table S2. Summary of GWAS data. Table S3. Baseline characteristics of UK Biobank participants. Tables S4-S5. Annotation of CPASSOC-identified SNPs. Tables S6-S7. 99% credible sets of CPASSOC-identified SNPs. Table S8. Results from colocalization analysis for CPASSOC-identified SNPs. Table S9. TWAS significant gene-tissue pairs for depression, overall BC, and subtype-specific BC.

Rights and permissions

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

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Wu, X., Zhang, W., Zhao, X. et al. Investigating the relationship between depression and breast cancer: observational and genetic analyses. BMC Med 21, 170 (2023). https://doi.org/10.1186/s12916-023-02876-w

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12916-023-02876-w

Keywords