Depressive symptoms predict longitudinal changes of chronic inflammation at the transition to adulthood

Background Inflammation is closely related to poor mental and physical health, including depressive symptoms and its specific symptoms. To reveal the linear and nonlinear relationships between depressive symptoms and chronic inflammation levels, and perform further analysis of the associations between symptom-specificity of depressive symptoms and inflammation among young adults by using a prospective design. Methods In this longitudinal study, we examined college students recruited from two universities in China, who were examined at baseline and 2-years follow-up. Depressive symptoms were measured by applying the Patient Health Questionnaire 9 (PHQ-9) at baseline. Plasma levels of four inflammatory biomarkers, including interleukin-6 (IL-6), interleukin-1β (IL-1β), tumor necrosis factor-α (TNF-α), and C reactive protein (CRP) were assayed at baseline and 2-year follow-up. In addition to the conventional generalized linear models, as well as restricted cubic splines were innovatively used to analyze the cross-sectional and longitudinal nonlinear relationships between depressive symptoms and inflammatory biomarkers. Results Generalized linear model analysis revealed that there were no statistical associations between depressive symptoms and any inflammatory biomarker levels. The results of the restricted cubic spline demonstrated a U-shaped nonlinear association between depressive symptoms and ΔIL-1β or ΔTNF-α (changes in baseline and 2-year follow-up), but these associations disappeared after adjusting the confounders. Symptom-specificity of depressive symptoms such as sleeping problems and suicidal ideation were associated with lower IL-1β at baseline or changes in IL-1β levels. Sleeping problems and psychomotor changes at baseline were associated with higher CRP at 2-year follow-up. Suicidal ideation at baseline was associated with changes in TNF-α levels. Conclusion Our findings suggested that symptom-specificity of depressive symptoms was associated with inflammation during a 2-year follow-up at the transition to adulthood. Simultaneously, more research is warranted to seek the directionality of depressive symptoms and chronic inflammation.


Introduction
The inflammatory response is a basic part of the body's innate immune mechanism to defend itself from detrimental stimuli, even though it is often considered a response to physical threats (1). Activation of the inflammatory pathway also affects the endocrine system, including stimulation of the hypothalamic-pituitary-adrenal axis (HPA), thereby causing down-regulation of glucocorticoid receptors, and decreased hormone response to inflammation, and subsequent inflammatory biomarker levels changed (2). Furthermore, common pro-inflammatory inflammatory biomarkers such as IL-6, IL-1b, TNF-a, CRP, etc. (3)(4)(5)(6), were presently recognized. Regardless, systemic inflammation can damage physical health, such as cardiovascular disease and chronic disease. Moreover, there is evidence that it is also strongly related to poor mental health, including major depressive disorders (7). Growing evidence recommends that inflammatory processes are connected with the pathophysiology of depressive symptoms and play a role in the etiology and course of mood disorders. Nonetheless, the precise duty of inflammation in the etiology of depression remains unknown (8, 9). Moreover, inflammation is associated with some but not all depressive symptoms, such as low mood and anhedonia and somatic/neurovegetative symptoms of fatigue, altered sleep, and appetite changes (10). Therefore, a symptom-based treatment may provide insights into the mechanisms of inflammationrelated depression.
In a cross-sectional study, individuals with depression had higher levels of IL-6 than healthy participants and frequently exhibited altered peripheral inflammatory profiles (11). Another cross-sectional study with adult males revealed significantly higher serum TNF-a levels in the depressive symptoms group than in the non-depressive symptoms group (12). It should be pointed out that there is a review that argues that the link between depression and inflammation may be bidirectional (13). A previous prospective study found a weaker bidirectional association between CRP and depressive symptoms in the elderly (14). Data from the UK Biobank manifested that higher inflammatory biomarkers levels were associated with an increased risk of depressive symptoms (15). A recent study reported depression risk score at 12 years old was positively linked to pro-inflammatory biomarkers levels 6 years later in the British (16). The incidence of depressive symptoms in children and adolescents is increasing sharply and the first onset of depression is often in adolescence or early adulthood (17). This is also an important window for identifying modifiable behavioral risk factors to prevent and control depression apparently (18). Nevertheless, most previous studies focus on middle-aged or elderly adults, there is a lack of prospective studies on depressive symptoms that causes elevated levels of inflammation at the transition to adulthood, which is an important phase in the development of depressive symptoms.
College students are in the critical stage of transition from adolescence to adulthood and healthy physical and mental development (19). Owing to most of the previous research participants being adults, there were few related studies on college students and the results were inconsistent obviously. Most research has indicated a linear positive correlation with higher levels of proinflammatory biomarkers in depressed individuals (11,12,20). On the contrary, some research unexpectedly displayed pro-inflammatory biomarkers were negatively correlated with depression (15). Additionally, rare studies have explored the nonlinear relationship between depression and chronic inflammation utilizing restricted cubic splines compared to the more traditional analyses in most studies. To address the indispensable evidence gap, this study was the first to focus on the linear and nonlinear association between depression and chronic inflammation in Chinese college students based on a longitudinal study.

Study population
The data came from the ongoing College Student Behavior and Health Cohort Study, which aimed to assess college students' health behavior and mental health in China. Based on the convenient sampling, a medical college and a comprehensive normal college were selected in Hefei City, Anhui Province, and Shangrao City, Jiangxi Province to recruit the freshman students. Then all freshmen from two faculty of each university were selected by cluster random sampling. A prospective cohort study with baseline and 4-wave follow-up at 6-month intervals in 2 years was conducted from freshman year to junior year. The baseline survey was conducted from April to May 2019. A total of 723 questionnaires and blood samples were collected at baseline. Electronic questionnaires were used, which were expected to take 15 to 20 minutes to fill out on a smartphone. Due to uncontrollable factors such as the COVID-19 epidemic and the health status of students, 248 questionnaires and blood samples were collected both at the baseline and 2-year follow-up, and no one was infected by the COVID-19 virus and designated as having experienced definite psychotic symptoms. The authors assert that all procedures used in this study comply with the ethical standards of the relevant national and institutional committees on human experimentation and with the Helsinki Declaration of 1975, as revised in 2008. The research protocol was approved by the Ethics Committee of Anhui Medical University (No. 20170291), and all participating students obtained written informed consent before completing the survey.

Depressive symptoms and symptom-specificity
The Patient Health Questionnaire (PHQ)-9 was applied to evaluate the symptoms of depression at baseline. Total scores for PHQ-9 range from 0 to 27 by summarizing scores of nine screening items with higher scores predicting more severe depressive symptoms. Depressive symptoms were defined as 5 points, 0-4 were divided into the no depressive symptoms group, and ≥ 5 were divided into the depressive symptoms group (21). In the present study, Cronbach's alpha coefficient was 0.917. Compared to clinically diagnosed depression, more depressive symptoms can be captured through self-report, which may lead to an overestimation of the prevalence of depression. However, there was also evidence that self-report measures of depressive symptoms were superior to clinician ratings, since rating scales used routinely in clinical practice may be impractical (22). Otherwise, PHQ-9 is recommended for research and clinical evaluation by the American Psychiatric Association and uses the DSM-IV diagnostic criteria to assess depressive symptomatology, such as sleep, concentration, energy problems, low self-esteem, anhedonia, etc. (23).
Symptom-specificity of depressive symptoms was recorded into "psychological symptoms" and "somatic symptoms" based on the reported PHQ-9 symptoms, including anhedonia, depressed mood, sleeping problems, fatigue, appetite changes, feelings of inadequacy, cognitive problems, psychomotor changes, and suicidal ideation.

Inflammatory biomarkers
Blood samples were collected at baseline and 2-year follow-up. From 6:00 to 8:00 in the morning, vacuum blood collection tubes (anticoagulation tubes) were used to take 5 mL of fasting venous blood for collection by medical professionals. Participants were instructed to avoid caffeinated foods (chocolate, coffee, tea, cola), alcohol, anti-inflammatory drugs such as aspirin and ibuprofen, and strenuous physical activity the day before venous blood collection. Blood samples were centrifuged at 3 000 rpm for 10 minutes within 2 hours. The upper plasma samples were split into 500 mL cryovials and stored in a -80°C refrigerator. The levels of inflammatory cytokine IL-6, IL-1b, and TNF-a in the supernatant were detected by a liquid-phase protein suspension chip detection instrument. The coefficient of variation of all analytes was less than 5%, and the coefficient of variation between batches of IL-6, IL-1b, and TNF-a was less than 20%, 15%, and 15%, respectively. The detection limit concentration values for IL-6, IL-1b, and TNF-a were 0.11 pg/mL, 0.14 pg/mL, and 0.16 pg/mL, respectively. The level of CRP was measured by immunoturbidimetry using the serum on the day the blood samples were collected.

Covariates
Due to the possible confounding influence of individual and family characteristics, we considered gender (males/females), age, residential area (rural/urban), any siblings (yes/no), selfreported family economy (low/medium or high), self-rated health condition (low/medium or high), a parental education level (primary school and below/middle school/senior high school and above) as covariates. Besides, we also considered some confounding factors that may affect the level of inflammation, including body mass index (BMI) calculated by height and weight, and risk health behaviors such as cigarette use (yes/no) and alcohol use (yes/no), which were classified by answering two modified questions based on the Youth Risk Behaviour Surveillance System questionnaire (24).

Statistical analyses
Nominal data are deliberated as frequency and percentage, and continuous data with Gaussian distribution are expressed as mean and standard deviation (SD). Group differences were tested with Chi-Squared, Kruskal-Wallis tests, or two independent sample t-tests. All statistical tests were two-sided and significance was set at P<0.05. A generalized linear model and restricted cubic spline model were realized to examine the linear and nonlinear associations of PHQ-9 scores with IL-1b, IL-6, TNF-a, and CRP levels. Changes between inflammatory biomarker levels at baseline and 2-year followup were quantified by calculation difference scores (D=follow-up minus baseline). Since the levels of IL-1b, IL-6, TNF-a, and CRP at baseline and follow-up all skewed distribution, logarithmic transformation was performed. Restricted cubic splines were applied to seek a possible nonlinear dependency between the measures. Restricted cubic splines consisted of polynomial functions that fit well with nonlinear relationships. The choice of knots depended on the setting of the research objectives. In the great majority of studies, 3 to 5 knots were generally selected. According to statistical test results, the choice of 3 knots was smoother than 5 knots and can better test the linear deviation. Three knots were prespecified in the 5th, 50th, and 95th percentiles, and one component of the spline function was generated (25)(26)(27). Controlling factors included residential area, any siblings, and parental education level. Further, we also explored the association between the symptomspecificity of depressive symptoms and inflammatory biomarkers by generalized linear models.
Sensitivity analyses were performed by using Chi-Squared tests (categorical variables) or t-tests (continuous variables) to compare the difference in characteristics of participants who were included in the analysis and those who were excluded because of dropping out.
The construction of the cubic spline model used the rms package in the software R version 4.1.2, all other statistical analyzes were performed using SPSS version 23.0.

Characteristics of the study population
As shown in Table 1, 723 college students aged 16 to 26 (M=18.68, SD=0.99) completed the online questionnaires, of which 238 (32.9%) were males. Blood sample collection was completed both at baseline and follow-up. At baseline, the prevalence rate of depressive symptoms was 40.1%.
Students from rural areas suffered from a higher risk of depressive symptoms compared to those from urban areas (P=0.003). Those students who reported lower family economy (P=0.046) or self-rated health conditions (P<0.001) were more likely to report depressive symptoms. The rates of depressive symptoms were higher in respondents who reported lower parental educational levels (P<0.05). Compared with nondrinkers or non-smokers, drinkers or smokers were more likely to report depressive symptoms. Gender, age, BMI, and any siblings did not differ between the participants with depressive symptoms and healthy participants at baseline.

Associations between baseline PHQ-9 scores and inflammatory biomarkers
Tables 2-4 depict the linear relationships between baseline PHQ-9 scores and inflammatory biomarker levels at baseline or 2-year follow-up, as well as the changes in inflammatory biomarker levels between follow-up and baseline. In the generalized linear models, individuals with higher PHQ-9 total scores at baseline were more likely to have lower CRP at follow-up and changes in CRP between follow-up and baseline. After controlling for residential area, self-reported family economy, self-rated health condition, father's education level, mother's education level, cigarette use, and alcohol use, the results remained no statistical significance. Figures 1-3 illustrate the nonlinear relationship between baseline PHQ-9 scores and inflammatory biomarker levels at baseline or 2year follow-up, as well as the changes in inflammatory biomarker levels between follow-up and baseline. As illustrated in Figure 1, baseline PHQ-9 scores were not found non-linearly associated with baseline inflammatory biomarker levels.
No nonlinear associations were found between PHQ-9 scores at baseline and inflammatory biomarker levels at 2-year follow-up. After adjusting for residential area, self-reported family economy, self-rated health condition, father's education level, mother's education level, cigarette use, and alcohol use, the results did not change significantly ( Figure 2). Figure 3 shows the U-shaped nonlinear associations of PHQ-9 scores with changes in IL-1b, IL-6, TNF-a, and CRP levels between 2-year follow-up and baseline. As PHQ-9 scores increased, changes in IL-1b decreased until the PHQ-9 scores were around 4 and increased eventually (P nonlinear = 0.013), However, the effects were attenuated and no longer significant after adjusting for covariates (P nonlinear = 0.059). Similarly, Ushaped associations were found between baseline PHQ-9 scores and changes in TNF-a levels (P nonlinear = 0.009), but no nonlinear associations persisted after controlling for residential area, self-reported family economy, self-rated health condition, father's education level, mother's education level, cigarette use, and alcohol use. No associations were found for change in IL-6 and CRP levels between 2-year follow-up and baseline.

Associations between the symptomspecificity of depressive symptoms at baseline and inflammatory biomarkers
As shown in Figures 4-6, after adjustment for residential area, self-reported family economy, self-rated health condition, father's    Inflammatory cytokines were log-transformed before analysis; the crude model was not adjusted by any variables, the adjusted model was adjusted by residential area, self-reported family economy, self-rated health condition, father's education level, mother's education level, cigarette use and alcohol use. B, regression coefficient; CI, confidence interval; IL-1b, interleukin-1b; IL-6, interleukin-6; TNF-a, tumor necrosis factor-a; CRP, C reactive protein. Frontiers in Immunology frontiersin.org education level, mother's education level, cigarette use and alcohol use, sleeping problems, and suicidal ideation were associated with lower IL-1b at baseline. Sleeping problems and psychomotor changes at baseline were associated with higher CRP at 2-year follow-up. Similarly, we also found sleeping problems and suicidal ideation at baseline were associated with changes in IL-1b levels. Suicidal ideation at baseline was associated with changes in TNFa levels. Other individual depressive symptoms were not associated with these inflammatory biomarker levels.

Sensitivity analysis
In a sensitivity analysis, we compared demographic differences between the dropout and the included samples, and the only statistical difference found was in gender. Compared to participants included in the study sample, participants who were not included were more likely to be females. Supplementary Table S1 shows detailed general characteristic information.

Discussion
To our knowledge, this is the first study on simultaneous analysis of linear and nonlinear relationships between depressive symptoms and inflammation at the transition to adulthood by a prospective study design, and further analysis of the associations between symptom-specificity of depressive symptoms and inflammation was performed. Our findings indicated a significant U-shaped nonlinear link between depressive symptoms at baseline and changes in IL-1b and TNF-a levels between 2-year follow-up and baseline, but these associations disappeared after adjusting the confounders. Symptom-specificity of depressive symptoms such as sleeping problems, psychomotor changes, and suicidal ideation were associated with inflammation. This highlights the need for further research into the biological mechanisms of depression in adolescence and the possibility of early intervention in this special population, which may help reduce the risk of adverse outcomes.
The transition from adolescence to adulthood has been seen as a decisive phase where mental health can undergo dramatic changes by changing roles and situations (19). For the reason, of the differences in assessment and research design, the results of the systematic review indicated that the prevalence of depressive symptoms in college students was obviously higher than that of the general population, which reported prevalence ranged from 10% to 85%, with a weighted average prevalence of 30.6%, which was also basically consistent with the 40.1% detection rate of depressive symptoms in our survey (28). In a longitudinal prospective study from childhood to adulthood, up to 72% of patients with major depressive disorder relapsed at intervals of 3-5 years after the first episode (29). In addition, individuals with the onset and recurrent episodes of depression in adolescence represented a particularly severe group and were strongly associated with impairment in multiple psychosocial domains (30). Moreover, early adult life is considered a specific risky age, and research on young adults with depressive symptoms across the transition period was of concern.
We identified a significant nonlinear association between depressive symptoms and changes in inflammatory cytokine TNF-a level between 2 years of follow-up and baseline at the transition to adulthood, in line with a previous systematic review supporting a link between chronic inflammation with depressive symptoms (31). Genome-wide association studies (GWAS) reported so far on depression only found a genetic association of TNF-a with depression (32). Furthermore, animal experimental Inflammatory cytokines were log-transformed before analysis; the crude model was not adjusted by any variables, the adjusted model was adjusted by residential area, self-reported family economy, self-rated health condition, father's education level, mother's education level, cigarette use and alcohol use. B, regression coefficient; CI, confidence interval; IL-1b, interleukin-1b; IL-6, interleukin-6; TNF-a, tumor necrosis factor-a; CRP, C reactive protein. a D represented changes in inflammatory biomarkers levels between 2-year follow-up and baseline.

FIGURE 1
Restricted cubic spline models of the associations between PHQ-9 total scores and inflammatory biomarkers at baseline (n=723). IL-1b, interleukin-1b; IL-6, interleukin-6; TNF-a, tumor necrosis factor-a; CRP, C reactive protein. a The blue plot represented the crude model which was not adjusted by any variables, and the red plot represented the adjusted model which was adjusted by residential area, self-reported family economy, self-rated health condition, father's education level, mother's education level, cigarette use and alcohol use.
results showed that TNF-a produced a depression-like state in mice reinforcing the idea that inflammatory components may play a crucial role in the pathophysiology of depressive symptoms (33).
The current study sheds light on a dose-response relationship between depressive symptoms at baseline and changes in inflammatory cytokine IL-1b levels after 2 years of follow-up at the transition to adulthood. A recent study showed the involvement of the IL-1 family in the early stage of depression, especially for pro-inflammatory cytokine IL-1b (34). Another study found that focused interventions for individuals with acute coronary syndrome with high IL-1b levels may reduce the risk of future depression (35). Our findings confirm and extend previous epidemiological studies suggesting that depression could be related to the development of chronic inflammation in young adults.
However, we found no statistical cross-sectional or longitudinal link between depressive symptoms and the inflammatory biomarkers IL-6 and CRP, which may be due to random measurement errors due to its diurnal variation (36). CRP was involved in the acute phase response and was sensitive to short-term effects such as infection, so there may be a mix of chronic and acute effects (16). Consistent with our findings, a recent longitudinal survey also reported no association between depressive symptoms and IL-6 (37). Whether there was an association between depression and inflammation, previous research was also conflicting (38,39). Another longitudinal study found that more severe depression at baseline predicted higher levels of IL-6 at follow-up (40). Settled levels of inflammatory cytokines such as IL-6 were often interpreted as a chronic, lowgrade inflammatory state. It should be noted, however, that IL-6 was not only a pro-inflammatory factor, but they were also antiinflammatory and play a role in tissue maintenance and repair as well as in immune preparation (6). The inflammatory system has been the focus of treatment trials for depression (41). Nevertheless, FIGURE 3 Restricted cubic spline models of the associations between PHQ-9 scores at baseline and changes in inflammatory biomarkers (n=248). IL-1b, interleukin-1b; IL-6, interleukin-6; TNF-a, tumor necrosis factor-a; CRP, C reactive protein. a The blue plot represented the crude model which was not adjusted by any variables, and the red plot represented the adjusted model which was adjusted by residential area, self-reported family economy, self-rated health condition, father's education level, mother's education level, cigarette use and alcohol use. b D represented changes between 2-year follow-up with baseline. Restricted cubic spline models of the associations between PHQ-9 scores at baseline and inflammatory biomarkers at follow-up (n=248). IL-1b, interleukin-1b; IL-6, interleukin-6; TNF-a, tumor necrosis factor-a; CRP, C reactive protein. a The blue plot represented the crude model which was not adjusted by any variables, and the red plot represented the adjusted model which was adjusted by residential area, self-reported family economy, self-rated health condition, father's education level, mother's education level, cigarette use and alcohol use. b F represented 2-year follow-up.
not all human beings with depression show increased inflammation. For example, there is great heterogeneity in the presentation of depression symptoms (42).
Our findings support the hypothesis that inflammation may be associated with certain specific symptoms such as sleep problems, suicidal ideation, and psychomotor changes in individuals with depression, which did not change statistically after adjustment for key controlling variables. A recent metaanalytic of 15 population-based cohort studies confirmed the symptom-specificity of the relationship between depression and inflammation (43). A prior study using data from UK Biobank and NESDA cohorts indicated that chronic inflammation was associated with somatic/neurovegetative symptoms of sleeping problems (10). Another research demonstrated chronic inflammation was not strongly linked to overall depression scores or severity, but rather related to specific features of depression marked by eating, appetite, and tiredness (44). The profile of symptoms identified can be used to define the subgroups of depressed individuals most likely to benefit from anti-inflammatory therapy. Our results may have vital implications for future research as it suggests a more targeted, symptom-focused approach to exploring the link between depression and chronic inflammation.
In addition to the causative pathway, depressive symptoms and inflammation may be associated, at least in part, by sharing pathophysiological processes. This was suggested by previous research manifesting that both depression and inflammation are associated with the same genetic variants (45,46). For example, inflammatory cytokine gene polymorphisms are associated with depression (47). There may be multiple interactions between depression and inflammatory cytokines. Traditionally, the biological mechanisms of depression regulation were complex. The impact of depression on chronic inflammation may be a multifactorial process involving the interaction of multiple mechanisms. More exploration is vital to further figure out these complex connections and to try to clarify their underlying mechanisms. We should focus attention on the need for further investigation and evaluation into the biological mechanisms FIGURE 4 Generalized linear models of the associations between symptom-specificity of depressive symptoms and inflammatory biomarkers at baseline (n=723). IL-1b, interleukin-1b; IL-6, interleukin-6; TNF-a, tumor necrosis factor-a; CRP, C reactive protein. a The blue plot represented the crude model which was not adjusted by any variables, and the red plot represented the adjusted model which was adjusted by residential area, self-reported family economy, self-rated health condition, father's education level, mother's education level, cigarette use and alcohol use. Generalized linear models of the associations between symptom-specificity of depressive symptoms at baseline and inflammatory biomarkers at follow-up (n=248). IL-1b, interleukin-1b; IL-6, interleukin-6; TNF-a, tumor necrosis factor-a; CRP, C reactive protein. a The blue plot represented the crude model which was not adjusted by any variables, and the red plot represented the adjusted model which was adjusted by residential area, self-reported family economy, self-rated health condition, father's education level, mother's education level, cigarette use and alcohol use. b F represented 2-year follow-up.
underlying depressive symptoms among adolescence and propose effective intervention methods. A growing body of research has recently linked an unhealthy lifestyle to depression and inflammation. Hence, lifestyle changes may be effective, low-cost preventive interventions (48, 49).
The sensitivity analysis, comparing the included and excluded samples, illustrate that the included participants were more likely to be males, which may lead to extra caution in the extrapolation of the results. Evidence indicated that depression and IL-6 were more closely correlated in women than in men, which may also lead to bias in our findings (50). Nevertheless, a systematic review suggested that gender differences in the association between depressive symptoms and systemic inflammation were not significant (51). More relevant studies comparing men and women would be required. Besides, because of the potentially confounding effects of the menstrual cycle and estrogen on inflammatory biomarkers, systematic studies targeting women are highly warranted (52,53).
Several strengths of this study should be addressed. Firstly, this is one of the few studies that comprehensively explore the causal associations linking depressive symptoms and multiple inflammatory biomarkers longitudinally. Most previous research was cross-sectional, our design was a step forward in evaluating the temporal association between depressive symptoms and chronic inflammation. Secondly, this research applied an innovative analytic approach to restricted cubic splines to analyze the nonlinear association of depressive symptoms with inflammation. Conversely, the contributions of this study are required to be considering its limitations in mind. Blood samples were only collected at two points to measure inflammatory biomarker levels, and the duration of follow-up may be a factor in the difference in results. However, no studies have shown how long the effects of depressive symptoms on chronic inflammation will last. Secondly, a large sample size lost to follow-up may lead to biased results, but we conducted a sensitivity analysis of participants who were included and excluded, except for gender, and general characteristics were statistically comparable. Thirdly, depressive symptoms were assessed by self-reporting rather than clinical diagnoses, which may not avoid recall bias and reporting bias. However, studies have found that convenient self-report measures of depression may be more helpful in predicting adverse outcomes than clinical scores (22, 23). Future research would preferably apply the clinical standard for the diagnosis of depressive symptoms.
In summary, our study identified that there was a significant nonlinear relationship between baseline depressive symptoms and changes in IL-1b and TNF-a levels during a 2-year follow-up, but these associations disappeared after adjusting the confounders. Symptom-specificity of depressive symptoms such as sleeping problems, psychomotor changes, and suicidal ideation were associated with inflammation during a 2-year follow-up at the transition to adulthood. Collectively, our results support the effect of depressive symptoms on chronic inflammation, and future research with larger samples and more well-designed waves of follow-up is imperative to verify our findings. Nonetheless, this information could be useful for the early identification of depressive symptoms and may provide clues for new symptombased psychological and pharmacological treatments for individuals with elevated levels of inflammation.

Data availability statement
The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation. Generalized linear models of the associations between symptom-specificity of depressive symptoms at baseline and changes in inflammatory biomarkers between follow-up and baseline (n=248). IL-1b, interleukin-1b; IL-6, interleukin-6; TNF-a, tumor necrosis factor-a; CRP, C reactive protein. a The blue plot represented the crude model which was not adjusted by any variables, the red plot represented adjusted model which was adjusted by residential area, self-reported family economy, self-rated health condition, father's education level, mother's education level, cigarette use and alcohol use. b D represented changes between 2-year follow-up with baseline.

Ethics statement
The studies involving human participants were reviewed and approved by Ethics Committee of Anhui Medical University (No. 20170291). Written informed consent to participate in this study was provided by the participants' legal guardian/next of kin.

Author contributions
SZ was responsible for study conceptualization, formal analysis, and writing the original draft of the manuscript. YQ, DZ, TL, YX, LZ, and YY contributed to data collection. XW and ST supervised the study, were responsible for funding acquisition, and reviewed and edited the manuscript. FT contributed to resource acquisition. All authors contributed to the article and approved the submitted version.