Characteristics of study participants are shown in Table 1. The study population consisted of 5,513 participants whom had all serum markers of the humoral immune system and inflammation measured on the same day. The study population contained more women (63.30%) than men (36.70%) and the mean age of the participants was 51.4 years. Most of the participants had attained ‘Middle’ (39.96%), compared to ‘Low’ (27.30%) and to ‘High’ (23.67%) level of education. Moreover, 83.98% of the individuals had no diseases recorded at the time the measurements were taken, which implies the healthy status of the study population.
Table 1 Descriptive statistics of study population. Participants had both serum markers of the humoral immune system and inflammation measured at the same measurement.
|
N = 5,513 (100%)
|
Sex
|
|
Male
|
2,023 (36.70)
|
Female
|
3,490 (63.30)
|
Age
|
|
Mean (SD)
|
51.4 (17.27)
|
< 40
|
1,523 (27.63)
|
40-50
|
1,165 (21.13)
|
50-65
|
1,460 (26.48)
|
> 65
|
1,365 (24.76)
|
SES*
|
|
Unclassified/Missing
|
1,224 (22.20)
|
Low
|
2,213 (40.14)
|
High
|
2,076 (37.66)
|
Education
|
|
Missing
|
500 (9.07)
|
Low
|
1,505 (27.30)
|
Middle
|
2,203 (39.96)
|
High
|
1,305 (23.67)
|
Charlson Comorbidity Index
|
|
0
|
4630 (83.98)
|
1
|
440 (7.98)
|
2
|
293 (5.31)
|
3+
|
150 (2.72)
|
Biomarkers
|
Mean (SD)
|
IgA (g/L)
|
2.65 (1.39)
|
IgG (g/L)
|
12.65 (3.91)
|
IgM (g/L)
|
1.34 (1.03)
|
CRP (mg/L)
|
6.04 (12.63)
|
Haptoglobin (g/L)
|
1.14 (0.36)
|
Albumin (g/L)
|
42.40 (3.23)
|
White blood cells (109 cells/L)
|
6.82 (3.37)
|
Iron (µmol/L)
|
16.27 (5.54)
|
Total iron-binding capacity (µmol/L)
|
60.27 (8.63)
|
* SES: Socioeconomic status
|
|
Table 2 shows a 9 by 9 scatter plot matrix in which all biomarkers are plotted against each other. All cells with the text in bold indicated correlations above or equal to +/-0.25 that imply interaction and were statistically significant at p-value< 0.0006 corrected for multiple testing. The strongest positive correlation was seen for CRP and haptoglobin (rs=0.50). The rest of the correlation were low as follow: haptoglobin and WBC (rs=0.25) and inverse correlations (haptoglobin and Albumin (rs=-0.31), CRP and Albumin (rs=-0.29), CRP and Iron (rs=-0.28)). Markers of the humoral immune system, IgG and IgA, showed a low correlation (rs=0.22). Moreover, low inverse correlations were observed between Albumin and IgA (rs=-0.26) and IgG (rs=-0.20).
Table 2 Spearman’s ρ rank-order Correlation matrix between all nine serum markers.
|
Biomarker Correlations
|
IgA
|
IgG
|
IgM
|
CRP
|
Albumin
|
Haptoglobin
|
WBC
|
Iron
|
TIBC
|
IgA
|
1.00
|
0.22
P < 0.0001
|
-0.07
P < 0.0001
|
0.11
P < 0.0001
|
-0.26
P < 0.0001
|
0.16
P < 0.0001
|
0.05
P < 0.0001
|
-0.01
P = 0.4639
|
-0.11
P < 0.0001
|
IgG
|
0.22
P < 0.0001
|
1.00
|
0.13
P < 0.0001
|
0.07
P < 0.0001
|
-0.20
P < 0.0001
|
0.02
P = 0.1536
|
-0.03
P = 0.0184
|
-0.01
P = 0.5675
|
-0.07
P < 0.0001
|
IgM
|
-0.07
P < 0.0001
|
0.13
P < 0.0001
|
1.00
|
0.02
P = 0.0972
|
-0.06
P < 0.0001
|
-0.04
P = 0.0034
|
-0.01
P = 0.2745
|
-0.01
P = 0.2809
|
-0.01
P = 0.4373
|
CRP
|
0.11
P < 0.0001
|
0.07
P < 0.0001
|
0.02
P = 0.0972
|
1.00
|
-0.29
P < 0.0001
|
0.50
P < 0.0001
|
0.17
P < 0.0001
|
-0.28
P < 0.0001
|
-0.11
P < 0.0001
|
Albumin
|
-0.26
P < 0.0001
|
-0.20
P < 0.0001
|
-0.06
P < 0.0001
|
-0.29
P < 0.0001
|
1.00
|
-0.31
P < 0.0001
|
-0.08
P < 0.0001
|
0.14
P < 0.0001
|
0.19
P < 0.0001
|
Haptoglobin
|
0.16
P < 0.0001
|
0.02
P = 0.1536
|
-0.04
P = 0.0034
|
0.50
P < 0.0001
|
-0.31
P < 0.0001
|
1.00
|
0.25
P < 0.0001
|
-0.25
P < 0.0001
|
-0.01
P = 0.2710
|
WBC
|
0.05
P < 0.0001
|
-0.03
P =0.0184
|
-0.01
P 0.2745
|
0.17
P < 0.0001
|
-0.08
P < 0.0001
|
0.25
P < 0.0001
|
1.00
|
-0.09
P < 0.0001
|
0.02
P = 0.0717
|
Iron
|
-0.01
P = 0.4639
|
-0.01
P = 0.5675
|
-0.01
P = 0.2809
|
-0.28
P < 0.0001
|
0.14
P < 0.0001
|
-0.25
P < 0.0001
|
-0.09
P < 0.0001
|
1.00
|
-0.06
P < 0.0001
|
TIBC
|
-0.11
P < 0.0001
|
-0.07
P < 0.0001
|
0.01
P = 0.4373
|
-0.11
P < 0.0001
|
0.19
P < 0.0001
|
-0.01
P = 0.2710
|
0.02
P = 0.0717
|
-0.06
P < 0.0001
|
1.00
|
Correlations above or equal to +- 0.25 are in bold. All these correlations are statistically significant (p-value<0.0006). (WBC = White blood cells, TIBC = Total iron-binding capacity)
The hierarchical clustering dendrogram is displayed in Figure 1. Hierarchical clustering formed two large clusters. The first cluster consists of CRP and haptoglobin on the first level. This first level clustered together with WBC on the second level and IgM on the third level. The second cluster consists of two clusters, IgA and IgG, which were clustered together on the second level with iron and albumin and TIBC clustering together.
Figure 1 Hierarchical clustering – Dendrogram displaying results of hierarchical clustering of all nine serum biomarkers. (WBC = White blood cells, TIBC = Total iron-binding capacity)
The scree plot of the principal component analysis displayed the variances of the nine principal components with associated proportion of variance and cumulative proportions (Figure 2). The first four components produced a cumulative proportion of variance of 61% of the data. The first component had a variance of 2 and accounted for 24% of the entire dataset. Component 2 had a variance of 1.4; 15% of the dataset. Component 3 had a variance of 1.1; 12% of the dataset and component 4 had a variance of 1.0, 11% of the dataset. As the first four components comprised the majority of the variance, no further analysis was performed. The loadings of the first four components presented a first component loaded with the acute inflammatory markers CRP and haptoglobin together with albumin, however presenting an inverse correlation. IgA, IgG and WBC were also contributing with an inverse correlation to this component. The second component was driven by the immunoglobulins, mainly IgG and a smaller inverse contribution of the acute inflammation markers. The third component was loaded mainly with IgM (Figure 2).
Figure 2 Principal component analysis – A scree plot showing the variance of the nine components. The table below is describing the loadings for each of the first four components of the principal component analysis which contain most of the variance of the data. The loadings of each component comprise the weight of each biomarker in the equation that define the particular component. For example, the first component comprises negative weights for IgA, IgG, CRP, Haptoglobin and WBC, and positive weights for albumin, iron and TIBC. (WBC = White blood cells, TIBC = Total iron-binding capacity).
The results of the multivariate analysis of variance are displayed in Table 3. Statistical difference in mean was observed for all variables between males and females, except for haptoglobin. The mean values of IgM and iron were not statistically different for the education categories. The mean values of the other biomarkers were statistically different between the Low category and the High education category. The mean values of IgG, IgA, Albumin and Iron markers were statistically different between medium and low education categories, meanwhile the mean values of CRP, Haptoglobin, WBC and TIBC markers were statistically different between Medium and High education. The mean values of haptoglobin and white blood cells were not statistically different from each other in the classes of socio-economic status (SES) and the mean values of IgG, Albumin and Iron were not statistically different between the Low and High category of SES. The mean values of IgA, IgM, CRP and TIBC were statistically different between the Low and High category of SES. The mean values for all the serum biomarkers in the different age categories were statistically different from each other. However, the mean values of IgA, IgM, CRP, Iron and TIBC were not statistically different from each other between the < 40 and 40-50 age categories.
Table 3 Multivariate analysis of variance (MANOVA) - mean values (standard deviation) in each log transformed biomarker for the given factors sex, education, SES and age. Same colour of cells indicates no statistical difference between categories, meaning that cells that share the same colour, independently of the colour, do not present statistical difference between the variance of the values. On the contrary, cells that present different colour are statistically different (Tuckey’s range test). *** P < 0.0001, ** P < 0.05, * P < 0.05 (SES = Socio-economic status, WBC = White blood cells, TIBC = Total iron-binding capacity)
Manova of log transformed biomarkers against exposures
|
|
IgA &
|
IgG &
|
IgM&
|
CRP&
|
Albumin&
|
Haptoglobin&
|
WBC&
|
Iron&
|
TIBC&
|
Sex
|
***
|
**
|
***
|
**
|
***
|
*
|
**
|
***
|
***
|
Male
|
0.94 (0.58)
|
2.51 (0.30)
|
0.01 (0.56)
|
1.17 (1.07)
|
3.76 (0.09)
|
0.09 (0.33)
|
1.88 (0.32)
|
2.76 (0.37)
|
4.07 (0.14)
|
Female
|
0.77 (0.57)
|
2.49 (0.29)
|
0.20 (0.53)
|
1.08 (0.98)
|
3.74 (0.07)
|
0.08 (0.28)
|
1.86 (0.30)
|
2.70 (0.37)
|
4.10 (0.15)
|
|
Education
|
***
|
**
|
*
|
**
|
***
|
***
|
***
|
**
|
***
|
Missing
|
0.94 (0.67)
|
2.54 (0.31)
|
0.08 (0.61)
|
1.25 (1.17)
|
3.71 (0.10)
|
0.11 (0.34)
|
1.92 (0.31)
|
2.67 (0.40)
|
4.06 (0.15)
|
Low
|
0.86 (0.60)
|
2.50 (0.31)
|
0.12 (0.57)
|
1.11 (1.03)
|
3.74 (0.08)
|
0.12 (0.30)
|
1.91 (0.32)
|
2.73 (0.37)
|
4.10 (0.14)
|
Medium
|
0.81 (0.57)
|
2.49 (0.29)
|
0.14 (0.53)
|
1.08 (1.00)
|
3.75 (0.07)
|
0.08 (0.30)
|
1.85 (0.30)
|
2.74 (0.37)
|
4.09 (0.14)
|
High
|
0.79 (0.54)
|
2.48 (0.28)
|
0.16 (0.55)
|
1.12 (0.96)
|
3.75 (0.08)
|
0.04 (0.28)
|
1.82 (0.29)
|
2.72 (0.37)
|
4.08 (0.14)
|
|
SES
|
***
|
***
|
**
|
**
|
***
|
*
|
*
|
***
|
***
|
Unclassified/missing
|
0.89 (0.62)
|
2.53 (0.30)
|
0.10 (0.60)
|
0.17 (1.07)
|
3.73 (0.08)
|
0.09 (0.32)
|
1.88 (0.31)
|
2.68 (0.38)
|
4.08 (0.15)
|
Low
|
0.78 (0.57)
|
2.49 (0.29)
|
0.16 (0.52)
|
1.07 (0.98)
|
3.75 (0.07)
|
0.08 (0.29)
|
1.87 (0.31)
|
2.73 (0.37)
|
4.10 (0.14)
|
High
|
0.85 (0.57)
|
2.48 (0.29)
|
0.12 (0.55)
|
1.13 (1.01)
|
3.75 (0.08)
|
0.07 (0.29)
|
1.86 (0.31)
|
2.74 (0.37)
|
4.08 (0.15)
|
|
Age
|
***
|
***
|
***
|
***
|
***
|
***
|
***
|
**
|
***
|
<40
|
0.74 (0.52)
|
2.53 (0.27)
|
0.22 (0.48)
|
1.02 (0.92)
|
3.78 (0.07)
|
0.01 (0.28)
|
1.83 (0.30)
|
2.74 (0.40)
|
4.11 (0.15)
|
40-50
|
0.78 (0.54)
|
2.46 (0.27)
|
0.19 (0.52)
|
1.00 (0.92)
|
3.76 (0.07)
|
0.07 (0.27)
|
1.88 (0.30)
|
2.73 (0.38)
|
4.10 (0.14)
|
50-65
|
0.84 (0.58)
|
2.47 (0.29)
|
0.08 (0.53)
|
1.16 (1.03)
|
3.74 (0.08)
|
0.11 (0.30)
|
1.86 (0.31)
|
2.74 (0.35)
|
4.09 (0.14)
|
>65
|
0.98 (0.64)
|
2.51 (0.32)
|
0.03 (0.64)
|
1.27 (1.14)
|
3.70 (0.08)
|
0.14 (0.32)
|
1.90 (0.32)
|
2.70 (0.36)
|
4.06 (0.14)
|
& - All the biomarkers have been log transformed.
The interaction analysis presented a statistically significant interaction effect between sex and age in the following markers, IgM, Albumin, WBC and TIBC (Table 4, supplementary figures 1-4 and supplementary table 1).
Table 4 Manova analysis accounting for interaction between external variables. P-value for the manova interaction terms for the independent predictors’ Socioeconomic status (SES), Sex, Age and Education Status. Significant p-values are highlighted in red (Bonferroni α/n = .05/15= .0033)
|
IgA
|
IgG
|
IgM
|
CRP
|
Haptoglobin
|
Albumin
|
WBC
|
Iron
|
TIBC
|
SES*Education
|
.1546
|
.1685
|
.5025
|
.4146
|
.3820
|
.6234
|
.2758
|
.5994
|
.8433
|
SES*Sex
|
.0379
|
.2349
|
.0524
|
.8614
|
.4734
|
.0075
|
.9991
|
.9601
|
.6564
|
Education*Sex
|
.5033
|
.1960
|
.2852
|
.6522
|
.3708
|
.6628
|
.0077
|
.4724
|
.0012
|
SES*Age
|
.1753
|
.7472
|
.7440
|
.5285
|
.2541
|
.0003
|
.0514
|
.1450
|
.8552
|
Education*Age
|
.7581
|
.0185
|
.6127
|
.8787
|
.2086
|
.5353
|
.9512
|
.6132
|
.0316
|
Sex*Age
|
.1171
|
.0076
|
.0015
|
.1113
|
.2690
|
<.0001
|
<.0001
|
.6929
|
.0002
|
Supplementary Figure 1 interactions between Sex and Age for IgM – Sex and Age interactions were statistically significant for IgM in Manova. The plot presents the dependent variable (biomarker) continuous versus the categories resulted of the combinations of the two independent variables studied.
Supplementary Figure 2 interactions between Sex and Age for Albumin – Sex and Age interactions were statistically significant for Albumin in Manova. The plot presents the dependent variable (biomarker) continuous versus the categories resulted of the combinations of the two independent variables studied.
Supplementary Figure 3 interactions between Sex and Age for WBC – Sex and Age interactions were statistically significant for Albumin in WBC. The plot presents the dependent variable (biomarker) continuous versus the categories resulted of the combinations of the two independent variables studied.
Supplementary Figure 4 interactions between Sex and Age for TIBC – Sex and Age interactions were statistically significant for Albumin in TIBC. The plot presents the dependent variable (biomarker) continuous versus the categories resulted of the combinations of the two independent variables studied.
Supplementary Table 1 presents the estimates for interaction between Age and Sex for the four biomarker IgM, Albumin, WBC and TIBC. Please note adelquart is the age variable. We use parametric manova to explore these interactions.