INTRODUCTION
Over the past few decades, persistent viral pathogens such as herpesviruses and chronic bacterial infections such as Helicobacter pylori have been hypothesized to be involved in the aetiology of cardiovascular disease (CVD) as well as other chronic diseases of ageing because of the ability of such pathogens to cause damage to various tissues and organs of the body via both direct and indirect pathways [Reference O'Connor, Taylor and Hughes1]. For example, herpesviruses as well as H. pylori have been detected in tissues of the vascular [Reference Izadi2, Reference Melnick3], central nervous [Reference Fraser4], and gastrointestinal systems [Reference Marshall and Warren5] where these pathogens may cause tissue damage, induce localized inflammation and/or directly influence cellular processes that contribute to the aetiology of CVD [Reference Epstein6], cognitive impairment [Reference Kountouras7, Reference Tarter8], and cancer [Reference Ferreira9, Reference Soroceanu and Cobbs10]. However, even in the absence of direct tissue invasion, these pathogens may cause damage to tissues and organs of the body indirectly via the process of molecular mimicry whereby antibodies targeted against such pathogens may cross-react with and attack host tissues expressing proteins homologous to those contained by the pathogen [Reference Epstein11]. Furthermore, infection with and reactivation of persistent pathogens may trigger the release of pro-inflammatory cytokines such as interleukin-6 (IL-6) and the acute phase protein C-reactive protein (CRP) [Reference Nazmi12, Reference Jackson13], elevated levels of which have been shown to be predictive of CVD as well as other chronic diseases with inflammation-related aetiology [14].
Findings from studies examining the role of individual pathogens in the aetiology of outcomes such as CVD [Reference Simanek, Dowd and Aiello15, Reference Sealy-Jefferson16], physical and cognitive impairment [Reference Tarter8, Reference Aiello17, Reference Aiello18], and metabolic disorders [Reference Jeon19, Reference Xia20] have, however, been mixed. Researchers have hypothesized that more important than the contribution of any single pathogen to the aetiology of various chronic diseases may be the cumulative damage resulting from exposure to multiple pathogens over one's lifetime [Reference Epstein11]. Indeed, in recent years numerous studies have aimed to examine whether there is a dose–response relationship between the total number of pathogens for which individuals are seropositive (i.e. total pathogen burden level) and a wide array of chronic disease outcomes [Reference Katan21–Reference Zhu24]. However, despite these hypothesized relationships few studies have examined the association between total pathogen burden level and all-cause mortality and those that exist have been conducted in selected older populations such as those with coronary artery disease (CAD) and/or from specific geographical regions [Reference Zhu25, Reference Elkind26]. In addition, in the majority of previous studies it has been assumed that there is a monotonic relationship between the total number of pathogens for which individuals are seropositive and mortality, irrespective of pathogen type. For this reason, a ‘black box’ paradigm has prevailed in which the role that seropositivity for specific combinations of pathogens may play in the aetiology of chronic disease and mortality has been largely ignored. Overall, whether a graded relationship between total pathogen burden and all-cause mortality exists and the role that seropositivity for specific combinations of pathogens may play in predicting all-cause mortality in the general US population remains unknown.
The objectives of this study are therefore: (1) to examine the association between total pathogen burden with four persistent pathogens [cytomegalovirus (CMV), herpes simplex virus (HSV)-1, HSV-2 and H. pylori] and all-cause mortality, (2) to examine the association between seropositivity for each specific combination of these pathogens and all-cause mortality and (3) to assess whether the total number or specific combination of pathogens for which individuals are seropositive better predict all-cause mortality in a US representative population of individuals aged ⩾25 years.
METHODS
Study population
Data were obtained from the National Health and Nutrition Examination Survey (NHANES) III (1988–1994), a population-based, multistage stratified probability survey which collected information on the health and nutrition of the US civilian non-institutionalized population. A total of 33 994 subjects were interviewed in NHANES III. Our study sample was limited to subjects that were aged ⩾25 years at the time of examination (N = 15 242), were tested for serostatus for CMV, HSV-1, HSV-2 and H. pylori (N = 6525), and had mortality follow-up to 31 December 2006 (N = 6522). All participants in NHANES III provided informed consent at the time of enrolment.
Laboratory analyses
CMV-specific immunoglobulin G (IgG) was measured by a commercially available enzyme linked immunosorbent assay (ELISA) (Quest International Inc., USA) with sera values near the ELISA cut-off confirmed with a second ELISA assay (bioMérieux Inc., USA). If there was disagreement between the results from the first two tests, an Immunofluorescence Assay (Bion International Inc., USA) was performed and results from this test were considered the final result [Reference Staras27]. HSV-1 and HSV-2 seropositivity was assessed by solid-phase enzymatic immunodot assays using purified glycoprotein gG-1 of HSV-1 and gG-2 of HSV-2 as the antigen [Reference Schillinger28]. Subjects aged ⩾20 years from phase 1 of NHANES III (1988–1991) were tested for H. pylori antibody using H. pylori IgG ELISA (Wampole Laboratories, USA) [29].
Measures
Total pathogen burden was constructed as a five-level categorical variable in which individuals were categorized according to the total number of the four pathogens (CMV, HSV-1, HSV-1 and/or H. pylori) for which they were seropositive (range 0–4, referent category 0). In addition, a 16-level categorical variable was constructed in which individuals were classified according to which of the 16 possible combinations of pathogens they were seropositive for: CMV+, HSV-1+, HSV-2+, H. pylori+, CMV+/HSV-1+, CMV+/HSV-2+, CMV+/H. pylori+, HSV-1+/HSV-2+, HSV-1+/H. pylori+, HSV-2/H. pylori+, CMV+/HSV-1+/HSV-2+, CMV+/HSV-1+/H. pylori+, CMV+/HSV-2+/H. pylori+, HSV-1+/HSV-2+/H. pylori+, CMV+/HSV-1+/HSV-2+/H. pylori+ or seronegative to all four pathogens (referent category). Mortality status as of 31 December 2006, of participants was obtained by NCHS primarily from the National Death Index; however, other sources of mortality status included indication of deceased status from the Social Security Administration, the Centers for Medicare and Medicaid Services, or death certificate review [30]. Cause of death was coded using the International Classification of Diseases (ICD), Ninth Revision up to 1998 and ICD, Tenth Revision for 1999–2006. All deaths before 1999 were re-coded by the NCHS into comparable ICD, Tenth Revision codes [30]. Follow-up time for each person was calculated as the difference in months between the NHANES III examination date and the last known date alive or censored [30]. Persons who survived the entire follow-up period were administratively censored on 31 December 2006.
Covariates
Covariates of interest included age, gender, race/ethnicity, education level, body mass index (BMI; kg/m2) and smoking status. Age in years at examination was self-reported and treated as a continuous variable. Gender was dichotomized as female and male. Race/ethnicity was self-reported as non-Hispanic white, non-Hispanic black, Mexican-American, or Other and individuals were categorized as non-Hispanic white, non-Hispanic black or Other (Mexican-American or Other were combined into Other due to small sample sizes). Years of education was self-reported and treated as continuous. BMI was computed by NCHS from weight and standing height, measured at examination, and treated as continuous. Smoking status was self-reported and categorized into never/former smoker (did not smoke ⩾100 cigarettes in one's lifetime or smoked ⩾100 cigarettes in one's lifetime but do not currently smoke) and current smoker (smoked ⩾100 cigarettes in one's lifetime and currently smoke).
Statistical analyses
Statistical analyses were performed using SAS v. 9.4, with SAS-callable SUDAAN, v. 10.0.1 (SAS Institute Inc., USA) [31]. All analyses used appropriate weights and adjustments for strata and clustering used in the complex study design in NHANES III. We first estimated the weighted population proportion of individuals seropositive for each combination of pathogens within each total pathogen burden level. Next, bivariate associations between total pathogen burden level, all-cause mortality and covariates of interest were assessed using Pearson's χ 2 tests of independence for difference in proportions and Cochran–Mantel–Haenzel tests for linear trend for categorical variables and independent sample t tests for difference in means across mortality status or analysis of variance for comparison of means across total pathogen burden levels, for continuous variables. Covariates associated with total pathogen burden level and mortality in bivariate analyses and/or hypothesized a priori to be relevant confounders in the association between total pathogen burden and mortality were included in the confounder-adjusted models.
Next, we used Cox proportional hazard models to estimate the crude and confounder-adjusted hazard ratio (HR) and 95% confidence interval (CI) for the association between total pathogen burden level and all-cause mortality treating seronegativity to all four pathogens as the referent category. Models were first unadjusted, then adjusted for sociodemographic characteristics (age, gender, race/ethnicity, education level) and last for clinical factors (BMI, smoking status). These analyses were repeated with models in which total pathogen burden level was treated as a continuous ordinal variable and the adjusted Wald F statistic was estimated to test whether there was a statistically significant linear trend for all-cause mortality with increasing total pathogen burden level. Last, we estimated the HR and 95% CI for all-cause mortality using Cox proportional hazard models in which the 16-level pathogen burden variable was included (referent category was seronegativity for all four pathogens) to estimate the independent and joint effects of seropositivity for each combination of pathogens and all-cause mortality. In addition, we conducted sensitivity analyses in which we used proc multtest in SAS v. 9·4 to adjust for multiple comparisons between each pathogen combination and mortality, using the false discovery methods described by Benjamini & Hochberg [Reference Benjamini and Hochberg32].
RESULTS
The number and population-weighted proportion of individuals seropositive for each combination of pathogens is shown in Table 1. Nearly 22% of individuals were seropositive for one pathogen; 29·5% for two pathogens; 29·1% for three pathogens; and 8·7% for all four pathogens (see Table 1). Overall, the greatest weighted proportion of the US population were CMV+/HSV-1+/H. pylori+, CMV+/HSV-1+ or HSV-1+ only (18·8%, 18·5% and 11·5%, respectively) (see Table 1). Table 2 shows population-weighted estimates of demographic and clinical characteristics of those included in our analyses and the unadjusted bivariate associations between covariates of interest and total pathogen burden level. The mean [±standard error (s.e.)] age of those included in our analyses was 47·7 (±0·4) years, 52·4% were female and the majority were non-Hispanic white (79·6%). Individuals had a mean (±s.e.) of 12·3 (±0·1) years of education, a mean (±s.e.) BMI of 26·5 (±0·1) kg/m2, and ~30% were current smokers. Older age, female gender, non-Hispanic black and Other race/ethnicity, lower education level and higher BMI were statistically significantly associated with increasing total pathogen burden level (see Table 2). The population-weighted proportion of individuals who died during the mean 15·0 years of follow-up was 21% (see Table 3). Those who died were more likely to be of older age, lower education level, higher BMI, less likely to be never smokers and less likely to be categorized as Other race/ethnicity (see Table 3).
Population-weighted proportion of individuals seropositive for each combination of pathogens by total pathogen burden level in subjects aged ⩾25 years in the National Health and Nutrition Examination Survey III (1988–1994)
CMV, Cytomegalovirus; HSV-1, herpes simplex virus-1; HSV-2, herpes simplex virus-2.
Population-weighted demographic and clinical characteristics by total pathogen burden level in subjects aged ⩾25 years assessed for mortality status in the National Health and Nutrition Examination Survey III (1988–1994)
s.e., Standard error.
† Two-sided P values for Cochran–Mantel–Haenszel tests for trend for categorical variables and analysis of variance for tests of linear trend for continuous variables.
‡ N = 6480 due to 42 subjects missing data on education level.
§ N = 6499 due to 23 subjects missing data on body mass index.
* P < 0·05.
Population-weighted demographic and clinical characteristics by mortality status in individuals aged ⩾25 years tested for pathogens in the National Health and Nutrition Examination Survey III (1988–1994)
s.e., Standard error.
† Two-sided P values for Pearson's χ 2 tests for difference in proportions or Cochran–Mantel–Haenszel tests for trend for categorical variables and t tests for difference in means for continuous variables.
‡ N = 6480 due to 42 subjects missing data on education level.
§ N = 6499 due to 23 subjects missing data on body mass index.
* P < 0·05.
Table 4 shows the HR (95% CI) from crude- and covariate-adjusted Cox proportional hazard models examining the association between total pathogen burden level and all-cause mortality, adjusting first for age, gender, race/ethnicity, and education level (see model 2) and second additionally controlling for BMI and smoking status (see model 3). In the crude model, we observed a statistically significantly graded relationship between increasing total pathogen burden level and all-cause mortality (P < 0·001 for linear trend). However, the rate of mortality associated with each pathogen burden level compared to those seronegative for all four pathogens was attenuated and no longer statistically significant in the fully adjusted model (see Table 4). In addition, there was no overall association (P = 0·528) nor a statistically significant linear trend (P = 0·733) for the association between total pathogen burden level and all-cause mortality (see model 3, Table 4).
Association between total pathogen burden level and all-cause mortality in subjects aged ⩾25 years in the National Health and Nutrition Examination Survey III (1988–1994)
† Model 1 unadjusted.
‡ Model 2 adjusted for age, gender, race/ethnicity, and education level (N = 6480).
§ Model 3 additionally adjusted for body mass index and smoking status (N = 6457).
|| Two-tailed P value for test for linear triend estimated from model with total pathogen burden level treated as a continuous ordinal variable.
* P < 0·05.
The independent and joint effects of seropositivity for each possible combination of pathogens on all-cause mortality, compared to those seronegative for all four pathogens (referent category) is shown in Table 5. In the crude model, there were statistically significant independent effects of CMV (HR 3·63, 95% CI 2·17–6·09), HSV-1 (HR 2·02, 95% CI 1·22–3·35) and H. pylori (HR 3·45, 95% CI 1·45–8·20) seropositivity on all-cause mortality and combined seropositivity for all four pathogens was associated with the highest rate of all-cause mortality (HR 7·85, 95% CI 5·27–11·70). In the fully adjusted model, however, the independent effects of each individual pathogen on all-cause mortality were no longer statistically significant (see model 3, Table 5). Moreover, seropositivity for all four pathogens no longer yielded the highest rate of all-cause mortality (HR 1·38, 95% CI 0·92–2·06). By contrast, after covariate adjustment, combined seropositivity for CMV and HSV-2 (HR 1·95, 95% CI 1·13–3·35) and for HSV-2 and H. pylori (HR 2·07, 95% CI 0·59–7·21) were associated with greater all-cause mortality rates than that yielded for all other combinations of pathogens (see model 3, Table 5), although the latter association was not statistically significant. After adjustment for multiple comparisons, the joint effect of CMV and HSV-2 on all-cause mortality was no longer statistically significant (P = 0·270).
Association between each pathogen combination and all-cause mortality in subjects aged ⩾25 years in the National Health and Nutrition Examination Survey III (1988–1994)
CMV, Cytomegalovirus; HSV-1, herpes simplex virus-1; HSV-2, herpes simplex virus-2.
† Model 1 unadjusted (N = 6522).
‡ Model 2 adjusted for age, gender, race/ethnicity and education level (N = 6480).
§ Model 3 additionally adjusted for body mass index and smoking status (N = 6457).
* P < 0·05.
DISCUSSION
To the best of our knowledge, this study is the first to examine the association between both the total number and specific combinations of pathogens for which individuals are seropositive and all-cause mortality in a large, nationally representative sample of US adults. After adjustment for age, gender, race/ethnicity, education level, BMI and smoking status, we did not observe a statistically significantly graded relationship between increased total pathogen burden level with CMV, HSV-1, HSV-2 and/or H. pylori and all-cause mortality. By contrast, we found that compared to those seronegative for all four pathogens, individuals who were CMV+/HSV-2+ and HSV-2/H. pylori+ had the greatest rates of all-cause mortality. Our results suggest that future studies should think outside the ‘black box’ paradigm in which the number but not type of pathogens for which individuals are seropositive has been the focus, in order to identify the specific combinations of pathogens for which seropositivity is most predictive of chronic disease and consequently mortality in the United States.
We are only aware of two previous studies that have examined the association between total pathogen burden and all-cause mortality [Reference Zhu25, Reference Elkind26]. Zhu et al. examined the association between total pathogen burden with CMV, hepatitis A virus, HSV-1, and/or HSV-2 and all-cause mortality in a mostly male cohort of Utah residents with existing CAD (mean age, 65·3 years) [Reference Zhu25]. The authors identified a statistically significant trend for all-cause mortality with increasing total pathogen burden level (P = 0·020); however, consistent with our findings, there were no statistically significant associations between each level of total pathogen burden and all-cause mortality in models adjusting for age, gender, and various CVD risk factors [Reference Zhu25]. Importantly, Zhu et al. did not control for socioeconomic status and the study population was older, from a single region of the United States and limited to individuals with a history of CAD [Reference Zhu25]. For these reasons, it is possible that the statistically significant linear trend observed by Zhu et al. may have resulted from residual confounding by socioeconomic factors and/or may reflect demographic or clinical differences between those included in this previous study and the general US population [Reference Zhu25].
In a study by Elkind et al., instead of categorizing individuals according to the total number of pathogens for which they were seropositive, the authors constructed an infectious burden (IB) index based on the strength of the association between the pathogens for which individuals were seropositive and rate of stroke in 1625 individuals (mean age 68·4 years) in Northern Manhattan [Reference Elkind26]. Specifically, the authors estimated the hazard rate of stroke associated with seropositivity for each of five pathogens (Chlamydia pneumoniae, CMV, HSV-1, HSV-2 and/or H. pylori) and then summed the β coefficients corresponding to the log hazard rate yielded for each pathogen for which individuals were seropositive [Reference Elkind26]. Elkind et al. then estimated the HR (95% CI) for all-cause mortality per 1 s.d. (0·33) increase in the IB index, finding that after adjusting for demographic and cardiovascular risk factors, there was no statistically significant association between increasing IB index and all-cause mortality (HR 1·12, 95% CI 0·99–1·26) [Reference Elkind26].
The results from the study by Elkind et al. suggest, similarly to our findings, that there is no statistically significant linear association between increasing total pathogen burden and all-cause mortality [Reference Elkind26]. Importantly, however, because a 1 s.d. increase in the IB index is not necessarily equivalent to seropositivity for an additional pathogen, it is difficult to directly compare the results from Elkind et al. [Reference Elkind26] to those observed in our study as well as those from the study by Zhu et al. [Reference Zhu25]. Moreover, given that the IB index constructed by Elkind et al. was based on the strength of the association between seropositivity for individual pathogens and rate of stroke, the associations observed by the authors may more closely reflect the rate of death associated specifically with increased risk of infection-related stroke than between infectious burden and mortality, per se [Reference Elkind26]. Last, although Elkind et al. constructed the IB index based upon the pathogens for which individuals were seropositive, because the authors did not report the rates of all-cause mortality associated with specific IB index scores, whether seropositivity for particular combinations of pathogens were most predictive of mortality in this previous study, cannot be ascertained [Reference Elkind26].
When we examined the independent and joint effects of each combination of pathogens on all-cause mortality, compared to those seronegative for all four pathogens, we found that there was not a strictly monotonic relationship between increasing total pathogen burden level and all-cause mortality, regardless of pathogen combination. For example, in the fully adjusted model, the combined effects of seropositivity for CMV and HSV-2 as well as HSV-2 and H. pylori (although not statistically significant) yielded higher rates of all-cause mortality than that yielded for all combinations of three pathogens as well as for seropositivity for all four pathogens. While CMV, HSV-2 and H. pylori have all been linked to numerous chronic disease outcomes including CVD [Reference Simanek, Dowd and Aiello15, Reference Sun33, Reference de Luis34] cancer [Reference Ferreira9, Reference Soroceanu and Cobbs10, Reference Thomas35] and cognitive impairment [Reference Tarter8, Reference Beydoun36], there is a paucity of studies examining the biological interaction between these pathogens. Indeed, we are unaware of any studies which have examined the interactive effects of HSV-2 and H. pylori in the aetiology of chronic disease outcomes or mortality. In vitro studies suggest, however, that HSV-2 replication may be induced by CMV [37, Reference Colberg-Poley, Isom and Rapp38], thus it is possible that HSV-2 and CMV may indeed biologically interact to increase rates of chronic disease and thus mortality, beyond that of other pathogens included in our study. For this reason, although the combined effect of CMV and HSV-2 seropositivity was no longer statistically significant after adjustment for multiple comparisons, further examination of the effect of seropositivity for this combination of pathogens on mortality may be warranted. Our findings also suggest that there may be negative biological interaction between and/or competing antagonism among pathogens acting through similar biological pathways such that there will not always be a strictly monotonic relationship between increasing total pathogen burden level and mortality. Importantly, violations of this assumption will go undetected in models which only consider the total number and not types of pathogens for which individuals are seropositive. For these reasons, future studies should consider thinking outside the total pathogen burden ‘black box’ to also consider the role that seropositivity for and/or biological interactions between specific combinations of pathogens may play in the aetiology of various chronic diseases and consequently mortality.
There are a few limitations to our study. First, we were unable to include other persistent pathogens that have also been associated with chronic diseases such as CVD as well as mortality including C. pneumoniae [Reference Strachan39] in our measure of total pathogen burden, because only selected pathogens were serologically tested for in NHANES III. For this reason, we cannot rule out that the associations between specific combinations of pathogens and all-cause mortality observed in our study are due to unmeasured co-infection with other persistent pathogens not measured in our study. Second, the number of individuals seropositive for certain combinations of pathogens was small (i.e. N = 13 for those seropositive for HSV-2+ and H. pylori+ combined), thus power to detect associations between all combinations of pathogens and mortality may have been limited. For these reasons, future studies which include a comprehensive array of pathogens, and possibly even microbiome data, are needed to fully understand which pathogens, alone or in combination, are most strongly associated with all-cause mortality in the general US population. Nonetheless, NHANES III provided us with the largest and only population-based study available for addressing questions related to total pathogen burden and all-cause mortality in the general US population.
Overall, our findings suggest that interventions targeting the prevention or treatment of particular combinations of persistent pathogens may be more effective for reducing all-cause mortality than those focused solely on reducing overall pathogen burden. For this reason, further investigations into the total pathogen burden ‘black box’ in which researchers focus on the number but not type of pathogens for which individuals are seropositive, may ultimately misinform preventive measures aimed at reducing infection-related death. Future studies should therefore aim to examine the role of specific combinations of pathogens in the aetiology of chronic disease and to identify which pathogens, in addition to those identified in our study, individually or in combination, are most predictive of mortality in the general US population.
ACKNOWLEDGEMENTS
This work was supported by National Institutes of Health (1-R21-NR-011181-01) (to J.B.D. and A.E.A.) and National Institutes of Aging, and National Institutes of Health (1-R01-AG-040115) (to J.B.D. and A.E.A.).
DECLARATION OF INTEREST
None.
