Introduction

Subjective well-being(SWB) is a subjective measure of one’s emotion and cognition, consisting of positive and negative affect as well as life satisfaction which stands for happiness1. In contrast, depression can be a serious disorder with significantly increased risk of physical and mental disabilities2,3. The depressive symptoms (DS) are warning signals for one’s psychological health and if not taken seriously, the individual can easily develop into clinical depression4. Happiness and depression are negatively correlated with each other5; they both show heritability at a rate from 30% to 40%, according to many twin and family studies6,7,8. Single nucleotide polymorphisms (SNPs) that are significantly associated with them were discovered as well9,10.

On the other hand, according to previous studies, personality traits had been widely reported as important predictors of subjective well-being (SWB)11,12 and depression10,13,14, which can usually be illustrated as a five-factor model (FFM), including Extraversion (E), Neuroticism (N), Conscientiousness (C), Agreeableness (A), and Openness to Experience (O)15. Twin studies suggested it with a heritability of nearly 40%16 and genome-wide association studies (GWAS) have probed several remarkable associating SNPs17.

In addition, correlation was found between personality and SWB as well as depression10,12,17 in a genetic level. This complex intrinsic correlation suggested a shared genetic inheritance of happiness and depression and we conjecture that personality traits can play an important role as a mediator.

The purpose of this study was to explore the possible shared genetic basis of SWB and DS through the MDD GWAS-supported SNP9 and the mediating role of personality.

Materials and methods

Participants

The participants in this study were 878 Han Chinese college freshmen. The mean age of the participants was 18.7 years (SD = 1.15), and 37.5% of the participants were male. Additional information about participant characteristics was given in Table 1. We also recruited 384 Han Chinese patients with MD (175 males and 209 females; mean age 35.4+/−12.18 years, range18–60 years) in this research to as a validation cohort. The MDD diagnosis were conducted by at least two well-trained psychiatrists using the Structured Clinical Interview for DSM-IV (SCID-I)18. All of them signed the informed consents and the study was appraised and confirmed by the Ethics Committee of the Jining Medical University and Shanghai human genetic resources ethics committee. All methods were carried out in accordance with relevant guidelines and regulation.

Table 1 Demographics and characteristics of the college sample.

The freshmen filled the questionnaire including demographic information, assessment of SWB, DS and personality as well as information feedback to the survey (efficacy, understanding, carefulness, significance). Data cleaning and curation was performed based on physical criterion, clustering analysis, feedback questions and questionnaire response time control (2% omitted in total).

Measurements

Assessment of subjective well-being (SWB) is composed of multiple measuring scales, including Satisfaction With Life Scale (SWLS19), the Positive and Negative Affect Scale (PANAS20) and Oxford Happiness Questionnaire (OHQ21). The final score of subjective well-being was calculated using Principal Component Analysis (PCA).

We measured depressive symptoms (DS) using the revised Chinese version of CES-D (Center for Epidemiologic Studies Depression Scale)22. It can be used as a self-report depressive symptomatology measurement for research in the general population. The CES-D scores range from 0 to 60, with the higher score correspond to more serious depression. Criteria: the total score ≤15 is considered as no depression, and the total score of 16~19 is considered as possible depression, and the total score > = 20 indicates definite depression (moderate to severe)23,24,25.

Personality was assessed by the Big Five Inventory (BFI) measurement on the Big Five dimensions including E, N, A, C and O26.

Supplementary Table 1 shows some details and scoring scheme of above scales.

DNA extraction and genotyping

We collected the peripheral venousblood from each freshmen and MD patient and applied Trizol protocol in DNA extractions. The major depressive disorder GWAS-supported variant rs454214 in TMEM161B-MEF2C9 were genotyped by time of flight (MALDI-TOF) mass spectrometer through MassARRAY Analyzer four platform (AGENA, San Diego, CA). All the primers were designed by My-Sequenom online soft-ware Assay Design Suite v2.0. Each tube in the polymerase chain reaction contained 10 ng genomic DNA dissolved in 5 μl buffer.

Statistical analysis

Descriptive statistics

For all analyses, we set statistical significance at P value <0.05. Supplementary Table 2 summarizes the R packages we use. Descriptive statistics were calculated using statistical packages in R Studio including the basic information of the sample, the mean and standard deviation of SWB, depression symptoms and personality. T-test was used to compare the groups of different genders.

Correlation analysis

In this study, correlation analysis and its heat map were achieved using “Corrplot” package in R Studio. Association analyses were conducted by R package “SNPassoc” with five different genetic models (codominant, dominant, recessive, over-dominant and log-additive models, respectively).

To test whether DS and SWB share the same affecting genetic variant, rs454214, we first analyzed the correlation between DS and SWB, and the correlation between rs454214 and the scores of both SWB and DS. Moreover, we performed another analysis to validate the association between rs454214 and clinical depression (387 cases with MD, 878 controls). Information of allele and genotype frequency were obtained by using the SHEsis webserver (http://analysis.bio-x.cn/myAnalysis.php) and the “SNPassoc” package. The result of Hardy-Weinberg equilibrium test was also generated this way.

Mediation analysis

We conjecture that personality play a mediating role in the path from rs454214 genotype to SWB/DS phenotype. To test this hypothesis, we performed another correlation analysis for 5 dimensions of personality with SWB/DS. We also conducted the same analysis for rs454214 and the 5 dimensions of personality.

Furthermore, we tried to verify the mediational models in which personality mediate the relationships between rs454214 and SWB/DS; the R package “semMediation” was employed in this analysis. And we used the “mediation” R package to validate the result of each mediator respectively with bootstrapping of the indirect effect. Bootstrapping was used over alternative tests to avoid Type 1 errors that may arise from non-normal distributions of an indirect effect27. It could also provide additional information about percentage of mediating variables explaining the association between X and Y, i.e., the genotype and phenotype.

Results

Degree of depressive symptoms and subjective well-being

Descriptions of statistical scales are shown in Table 1. The mean score of CESD is 13.77, belonging to “no depressive symptoms” level (<=15). The degree of satisfaction with life is “Neutral” (20). And the average score of OHQ reach the level of “Rather happy” (4–5). In this study, there shows no significant difference in terms of DS, SWB, and personality scores between males and females (P > 0.05).

After Cattell Scree Test28, Parallel Analysis29 combined with Kaiser-Harris Criterion30 of Principal Component Analysis (PCA), the first principal component was taken as the final measurement score of Subjective Well-being.

Rs454214 is a shared variant of SWB/DS

Correlation analysis showed that SWB was strongly negatively correlated to DS (β = −0.85, P < 0.001, Fig. 1). We tested five genetic models and found significant association between the candidate SNP and SWB/DS in three of them, i.e., codominant, recessive, and overdominant, (P < 0.05, Table 2). Codominant model was then excluded due to 95% confidence interval constraint (contains zero). In the overdominant model, TC genotype referred to a higher state of SWB and lower score of DS (compared to TT and CC), while the recessive model indicated that TC and TT genotypes play a more positive role in SWB comparing to the CC genotype.

Figure 1
figure 1

Heat map illustrating correlations between personality traits and depressive symptoms and subjective well-being. The values in the lower triangular matrix is correlation coefficients. In the upper triangular matrix, big circles in dark colors indicates strong correlations. Positive correlations are showed in red while negative correlations are showed in blue. DS: depressive symptoms, SWB: Subjective well-being, BFI-E: Extraversion, BFI-A: Agreeableness, BFI-C: Conscientiousness, BFI-O: Openness to experience, BFI-N: Neuroticism. All correlation coefficients were significant at 0.001 level. *p < 0.05; **p < 0.01; ***p < 0.001.

Table 2 Association between rs454214 and subjective well-being as well as depressive symptoms.

For validation purpose, we recruited 384 Han Chinese patients with MD to examine whether rs454214 is also associated with clinical depression. It turned out that rs454214 passed the Hardy−Weinberg equilibrium test (HWE, P > 0.05) both in MD (case) and college samples (control). The genotype frequency of rs454214 between cases and controls was consistent with the recessive model (P < 0.05) according to five genetic models analyzed by R package “SNPassoc”.

Correlation between personality, SWB and DS

Correlation analysis indicates that the phenotyes of the five personality traits are correlated with SWB/DS. The Heat map (Fig. 1) illustrates that Subjective Well-being is positively correlated with E, A, C and O, and is negatively associated with N in a significant level of 0.001. On the other side, E, A, C and O are negatively associated with DS, while N is positively associated with DS. Likewise, E, A, C, O, and N have significant inter-correlation relationships with one another (P < 0.001).

Association analysis of personality

After the association study between personality and rs454214 (Supplementary Table 3), different genotypes of rs454214 are connected to all the five aspects of personality, fitting the overdominant model. N and A were also significant in codominant model (P < 0.05) but failed to fitting the 95% confidence interval constraint (contains zero).

It showed that the CT genotype of rs454214 is associated with higher score of E, A, O, C, and lower score of N. In the previous two analyses hereinabove, we had found the correlation between personality and SWB/DS as well as rs454214. Altogether, the inter-correlations among the variables provided initial support to the hypothetical conjectured indirect effects.

Mediation analysis

The multi-mediation analysis and the results in the plot generated by “semMediation” (Fig. 2) revealed the indirect effects of rs454214 on subjective well-being through E, A, O, N (P < 0.001) and those on DS through E, A, and N. Heterozygotes (TC genotype) can act as a protective factor, reducing the level of N and increasing the level of A and E to inhibit depression, and promoted happiness by raising the level of N and suppressing the degree of E, A, and O.

Figure 2
figure 2

Multi-mediation plot shows indirect effect between rs454214, personality traits, subjective well-being and depressive symptoms. Full lines indicate significant mediated paths in the full model (p < 0.05). Dashed lines indicate insignificant ones (p > 0.05). Estimates of indirect effect on subjective well-being are shown in red. Estimates of indirect effect on depressive symptoms are shown in blue. Estimates of total effect on them are shown in black. *p < 0.05; **p < 0.01; ***p < 0.001.

We then tested whether aspects of personality fully or partially mediate the relationship between rs454214 and SWB/DS, by examining whether the direct effect of rs454214 on SWB/DS (controlling for personality) was statistically significant. As demonstrated in Supplementary Fig. 1, the direct effects of rs454214 on subjective well-being and DS are not significant after controlling the impacts of mediators (β = 0.02, 0.09 respectively, P > 0.05), supporting an inference of a full mediation.

Supplementary Table 4 shows the verifying result of individual mediation model test by the “Mediation” R package. It is shown that the total effect of each model is significant (P < 0.05). We also found that rs454214 facilitated or suppressed SWB/DS through personality (E, A, O, N/E, A, N) with significant P value of indirect effect. As zero was not contained in the 95% confidence interval, the conjecture of indirect effect was supported. Direct effect and proportion information can also be found in the Supplementary Table 4. None of the direct effects were significant indicating the full mediation model. Among the five personality aspects, the mediating effect of N explained the highest significant percentage of the total effect (71% for both SWB and DS, P < 0.05), followed by E (67% for SWB and 61% for DS, P < 0.05) and A (63% for SWB and 59% for DS P < 0.05). O could explain 46% of inheritance of SWB due to rs454214 (P < 0.05).

Discussion

The result demonstrated a relatively good situation of psychological health state in the college participants. Whether gender affects SWB and DS was controversial, and in this study, we did not find significant difference between males and females in terms of SWB/DS and personality traits. But it could vary from center to center.

Different from previous studies which measure SWB from a single aspect, multiple scales were adopted to measure SWB for each individual and PCA (principal component analysis) was additionally adopted, according to the SWB model proposed by Diener1, making our measurement of SWB more thorough and robust.

We found rs454214 in TMEM161B-MEF2C9 to be a shared genetic loci of DS and SWB in health college participants. Rs454214 was also linked to clinical depression in patients with MDD. It supports our conjecture that happiness and depression partially share the same genetic mechanism. Adolescent depression is difficult to diagnose and easy to be neglected, making depressive symptoms important warning signals25. Our direct replication of the genome-wide significant associations with depressive symptoms in an independent depression sample provided further confirmation for this knowledge. Therefore, cares for the mental health of college students and depression prevention are critically important; previous studies also share this point of view31,32.

We found that rs454214 might affect DS and SWB in three genetic models: codominant, recessive, and overdominant models. It showed that the TT genotype was the risk factor of DS while the TC genotype was the protective factor. The superiority of heterozygotes was evident in the result. This inverted U-shaped relation was similar to the findings in sickle-cell disease and cognitive control processes. It was recently found that homozygotes of βS sickle mutation can lead to sickle-cell disease while heterozygotes contribute to higher protection33. The heterozygotes of CHRNA4 polymorphism (rs1044396) showed higher activity during cognitive engagement than the homozygotes in the cingulo-opercular (CO) network34.

Rs454214 is located in the upstream of MEF2C gene, which is a transcription factor regulating neuron maturation processes35,36,37. It can promote cognition (learning and memory) through negative feedback regulation that inhibits the over-generation of excitatory synapses and rebuilds the neural network continuously. It is highly expressed in the frontal cortex, entorhinal cortex, dentate gyrus, and amygdala of CNS38. These brain regions play essential roles in reward, disgust, stress, and emotional control, which may lead to the phenotype of depression and happiness39,40. Previous studies also found strong correlation between personality traits and these brain regions41,42, so it is no surprise for us to identify that rs454214 is associated with personality traits and that it affects happiness and depression.

As we found a pairwise causal relationship between rs454214, personality traits and DS/SWB combining with results of the previous studies, mediation analysis was suitable for this study to explain the effect of rs454214, personality traits on DS/SWB. Recent studies reported that Bayesian Structural Equation Model (BSEM) can become an equivalent model for the general Mediated Model we use43. Different from standard practice relies on frequentist methods (Bootstrapping and Sobel method) we use, alternative Bayesian approach employs Monte Carlo method and is easier to work with complex models44.

Besides suggesting the mediating role of personality, the result of mediation analysis also revealed the differences between different personality traits. In previous research on the relationship between personality, subjective well-being and depression, N and E were regarded as the most prominent personality trait45,46. Neuroticism is a dimension of personality which often presents negative emotions such as depression and is negatively correlated with happiness12,13. E on the other hand, showed opposite pattern as N14, which were again confirmed in our result of correlation analysis and mediation analysis. Furthermore, N and E accounted for the highest proportion of mediating effect between genetic variant (rs454214) and SWB, as well as DS, indicating the importance of E and N in the inheritance of happiness and depression.

In contrast, most of the studies found that A, C, and O had little or no effect on happiness and depression12,13,47. In our study however, we found that A played an important role in the genetic of subjective well-being and depressive symptoms. This discovery was supported by another previous study with Han Chinese sample48. It illustrated the complexity of happiness and depression that the influence of personality on them can be affected by other factors such as ethnic background, cultural differences, life events. In addition, O also mediated the inheritance of happiness to certain extent. Previous study indicated that this may be achieved by affecting the efficiency of information processing within the brain49.

One of the limitations of this study is the relatively small size and that only one SNP was intensively exploited for limited candidate SNP pool. More genetic candidates of depression and happiness will be included with larger sample size in the future work. In addition, the Mediated Model we used was relative simple and only included personality so far. Therefore, even the present model partly demonstrates the important role of personality in the inherited process of DS and SWB, it inevitably excluded effect from some other underlying factors and may not accurately estimating the whole complex genetic and psychological network. More factors such as social support, which is believed to play a more complex role in the existing multi-mediation model50,51, learning and memory,which can be influenced by gene MEF2C through regulating synaptic transmission35,52, will also be included in future work. With the increase of considering factors, machine learning models like BSEM (Bayesian network for continuous variable43) and mixed graphical models (both continuous and discrete variables53) may be fit for the construction of complex network. We will also conduct longitudinal research to track how heredity, personality and other factors influence happiness and depression as participants progress through college study.

In conclusion, this study found a relatively comforting situation of the mental health condition of college participants and found the major depressive disorder GWAS-supported variant rs454214 was a shared genetic variant of the subjective well-being and depressive symptoms, which is highly possibly mediated by personality traits. These findings add supporting evidences to the genetic mechanisms of happiness and depression and implies that more attentions need to be paid to college students’ mental health from the perspective of cognition and personality so that early prevention of depression can be achieved.