Skip to main content

The association between Toxoplasma infection and mortality: the NHANES epidemiologic follow-up study

Abstract

Background

Toxoplasma gondii has been reported to be associated with higher mortality in patients with schizophrenia. This study aimed to explore the relationship between T. gondii infection and 25-year mortality based on data from the Third National Health and Nutrition Examination Survey (NHANES III) database.

Methods

Cases with serum T. gondii antibody test results were included in this study and the corresponding mortality dataset was obtained from the US National Center for Health Statistics (NCHS). Propensity score matching (PSM) was used to match age and sex between groups. The Cox proportional hazards model was used to evaluate the effect of T. gondii infection on mortality.

Results

A total of 14,181 cases were included in the analysis, of which 3831 (27.0%) were seropositive for T. gondii antibody. The median follow-up time of the whole cohort was 22.5 (interquartile range 16.3, 24.5) years. A total of 5082 deaths were observed in this cohort, a mortality rate of 35.8%. All-cause mortality was significantly higher in the seropositive group than in the seronegative group (50.0% vs 30.6%, P < 0.001). Kaplan–Meier analysis showed a significant difference in the survival time between two groups before and after PSM. Multivariate analysis showed that T. gondii infection was independently associated with higher all-cause mortality after adjusting for potential confounders.

Conclusions

Toxoplasma gondii infection is associated with higher mortality in general population.

Graphical Abstract

Background

Zoonoses continue to be a major public health concern [1]. Toxoplasma gondii is a common protozoan parasite that infects approximately 30% of the world’s population [2, 3]. In healthy individuals, T. gondii infection usually results in mild illness most of the time, but it has the potential to cause severe consequences. Toxoplasmosis is a common opportunistic infection in immunocompromised patients and results in high mortality in this patient population [2, 4, 5].

Toxoplasmosis has been linked to a number of diseases, including neuropsychiatric diseases, such as primary neuropathies, behavioral and psychiatric disorders [6, 7], as well as some related metabolic diseases, such as non-alcoholic fatty liver disease [8] and diabetes [9, 10]. Toxoplasma gondii infection is also associated with a worse prognosis of all-cause mortality in individuals with these diseases; for example, it has been shown to be associated with higher all-cause mortality rate in patients with schizophrenia [11, 12]. However, it remains unclear whether T. gondii infection actually affects the all-cause mortality rate in the general population.

The present study aimed to explore the relationship between T. gondii infection and 25-year mortality in the general population based on data retrieved from a nationally representative longitudinal database.

Methods

Study population

The study population was from the Third National Health and Nutrition Examination Survey (NHANES III) 1988–1994, which was a periodic survey conducted by the US Centers for Disease Control and Prevention. The study data are freely available for use by researchers from all over the world. We selected cases with serum T. gondii antibody test results for the present study. All individuals were followed up until December 2015, and the data on survival status and death dates were collected. The corresponding mortality dataset was obtained from the US National Center for Health Statistics (NCHS), which is also freely accessible. All NHANES protocols had been approved by the NCHS ethics review board, and participants or their proxies had provided informed consent prior to participation. As the present study only utilized anonymized data available to the public, no additional ethics approval was required.

Definitions

The presence of T. gondii antibody was tested using the indirect enzyme immunoassay method. The diagnostic threshold for T. gondii infection was 7 IU/ml as described in the original report. Samples with results ≥ 7 IU/ml were considered to indicate T. gondii infection. The entire dataset and details on the methods used can be retrieved from the NHANES website (https://www.cdc.gov/nchs/nhanes/about_nhanes.htm).

Hypertension was defined as a systolic blood pressure ≥ 135 mmHg or diastolic blood pressure  ≥ 85 mmHg. Patients with a significant history of hypertension and/or undertreatment of hypertension were also included.

People with type 2 diabetes were defined as having a history of diabetes, currently using insulin or oral hypoglycemic agents, with a fasting blood glucose (FBG)  ≥ 7.0 mmol/l or glycated hemoglobin (HbA1c) ≥ 6.5% or 2-h postprandial glucose  ≥ 11.0 mmol/l.

Body mass index (BMI, kg/m2) was defined as the weight in kilograms divided by the square of the height in meters.

Income was measured based on the poverty income ratio, an index reflecting the ratio of household income to the household poverty level determined by area of residence and household size. Low family income was defined as a ratio below 1.

In NHANES III, data on education was collected as a numerical variable (number of years of education), ranging from zero to 17 years. Low educational level was defined as an educational level of lower than the 12th grade.

Laboratory examinations

Routine laboratory examinations, including blood routine tests, C-reactive protein (CRP), FBG, HbA1c, alanine aminotransferase (ALT), aspartate aminotransferase (AST), albumin, serum, creatinine, blood urea nitrogen (BUN) and uric acid, among others, were all collected from the dataset.

Statistical analysis

Categorical variables were expressed as percentages. Continuous variables were expressed as means ± standard deviation. The Student t-test (for variables normally distributed), the Mann–Whitney U-test (for variables non-normally distributed) and the Chi-square test (for categorical variables) were used to compare the differences between the seropositive and seronegative groups. We used propensity score matching (PSM) to select age- and sex-matched participants between the two groups. Patients were matched 1:1 based on their propensity scores, and the match tolerance value was 0.0001. The Cox proportional hazards model was used to evaluate the effect of T. gondii infection on all-cause mortality. All tests were two-tailed and a P-value < 0.05 was considered to be statistically significant. All analysis was conducted in R 3.6.2 (https://www.r-project.org/).

Results

Baseline characteristics of the study population

A total of 19,598 participants with follow-up data were screened for this study (Fig. 1). After excluding participants without T. gondii antibody test (N = 4,085), 1332 without other key data, 14,181 cases were included in final analysis. Among them 6,672 (47.0%) were male and the average age was 46.3 ± 19.5 years old. A total of 2128 (15.0%) had diabetes and 5859 (41.3%) had hypertension. A total of 3831 (27.0%) cases were seropositive in T. gondii antibody. More details were shown in Table 1.

Fig. 1
figure 1

Flow chart of case selection. NHANES, Third National Health and Nutrition Examination Survey

Table 1 Comparison of Toxoplasma gondii-seropositive and -seronegative groups before propensity score matching

Comparison between the T. gondii-seropositive group and the T. gondii-seronegative group before and after PSM

According to the results of the T. gondii antibody test, participants were divided into a T. gondii-seropositive group and a T. gondii-seronegative group. The comparison of the baseline characteristics before PSM are shown in Table 1. Patients in the seropositive group were significantly older than their seronegative counterparts (mean ± SD, 53.7 ± 19.7 vs. 43.6 ± 18.7; P < 0.001). There was a higher proportion of males in the seropositive group than in the seronegative group (49.9% vs. 46.0%; P < 0.001). Diabetes, hypertension and metabolic derangements were more common in the seropositive group.

As age and sex are factors closely associated with life expectancy, and a significant difference in these factors was found between the seropositive and seronegative group, we used PSM to match age and sex between two groups. A total of 3806 pairs of cases were selected in which age and sex were comparable after PSM (Table 2). It was noteworthy that the severity of metabolic derangements, the proportion of patients with diabetes and BMI level remained significantly different between two groups even after matching for age and sex.

Table 2 Comparison between the T. gondii-seropositive and -seronegative groups after propensity score matching

All-cause mortality in T. gondii-seropositive and -seronegative groups

The median follow-up time of the whole cohort was 22.5 (interquartile range 16.3, 24.5) years. A total of 5082 deaths were observed in this cohort, which is a mortality rate of 35.8%. All-cause mortality was significantly higher in the seropositive group than in the seronegative group (50.0% vs. 30.6%; P < 0.001). Kaplan–Meier analysis showed a significant difference in the survival time between two groups (Fig. 2a).

Fig. 2
figure 2

Kaplan–Meier analysis between the Toxoplasma gondii-seropositive group and the T. gondii-seronegative group before PSM (a) and after PSM (b). The P-values were calculated by the log-rank test. c Forest plot of the hazard ratio for mortality based on Cox regression analysis. CI Confidence interval, HR hazard ratio, PSM Propensity score matching

After PSM, the difference in the mortality rate between the seropositive group and seronegative group was attenuated but remained statistically significant (49.7% vs. 47.3; P = 0.037). The same result was obtained in the Kaplan–Meier analysis (P = 0.014; Fig. 2b).

Cox multivariate analysis for the impact of T. gondii infection on all-cause mortality after PSM

We used Cox analysis to explore the correlation between T. gondii infection and all-cause mortality in the PSM population. The results are shown in Table 3. In the univariate analysis, T. gondii antibody seropositivity, race, sex, age, low educational level, hypertension, diabetes, waist circumference, related metabolic parameters and liver function and renal function test results were significantly different between the survival group and mortality group. Variables with a P-value < 0.05 in Table 3 were considered to be candidates for inclusion in the multivariate regression. Collinearity analysis showed a strong collinearity between BMI and waist, FBG and HbA1c, AST and ALT, and creatinine and BUN. For the variables with collinearity, we used the receiver operating characteristic (ROC) analysis to select the variable with a higher area under the ROC curve. The results showed that waist, HbA1c, AST and BUN were the most suitable for inclusion in the multivariate Cox regression analysis.

Table 3 Comparison between the survival and all-cause mortality groups after propensity score matching

In order to explore the association between T. gondii infection and mortality, we established three models with the combination of different variables to adjust for different confounding factors. Model 1 was adjusted for race, sex and age; model 2 was adjusted for race, sex, age, low educational level, diabetes and hypertension; and model 3 was adjusted for additional metabolic profiles and biomarkers of liver and renal functions based on model 2. The results were sustained as being statistically significant in the different models (model 1: hazard ratio [HR] 1.102, 95% confidence interval [CI] 1.033–1.176, P = 0.003; model 2: HR 1.081, 95% CI 1.013–1.154, P = 0.018; model 3: HR 1.144, 95% CI 1.070–1.224, P < 0.001), indicating the association of T. gondii infection and all-cause mortality was independent of other confounders (Table 4; Fig. 2c).

Table 4 Multivariate COX regression analysis after propensity score matching

The difference in risk factors for mortality in the T. gondii-seropositive and T. gondii-seronegative groups

The results of the Cox multivariate analysis in the T. gondii-seropositive and T. gondii-seronegative groups are shown in Additional file 1: Table S1. The effect of most factors was similar in the two groups, with the exception of CRP. CRP was significantly associated with the risk of mortality in the T. gondii-seropositive group (HR 1.070, 95% CI 1.020–1.121, P = 0.005), while no statistical significance was found in the seronegative group (P = 0.121).

Discussion

The present study used data from NHANES III to assess the longitudinal relationship between T. gondii infection at baseline and subsequent mortality after adjusting for relevant confounding factors. The main finding of this prospective study is that T. gondii antibody positivity was associated with increased all-cause mortality.

Toxoplasmosis is generally regarded to be a threat to immunocompromised patients or those with congenital toxoplasmosis [2, 13]. Several studies have also shown a higher risk for mortality from natural causes in persons with schizophrenia infected with T. gondii [4, 11, 12, 14]. Our study is the first to show that T. gondii infection is associated with higher mortality in the general population.

Several possible hypotheses might explain this finding. First, T. gondii infections are associated with elevated biomarkers of chronic inflammation, such as CRP [15]. Systemic chronic inflammation can lead to several diseases and cause disability and mortality worldwide [16]. The results of the present study show that the CRP level was associated with higher mortality in the T. gondii-seropositive group, indicating the possible role of chronic inflammation in an increased mortality. Second, T. gondii can influence host metabolism of fatty acids, lipids and energy in the liver [17]. Higher seroprevalence of T. gondii antibody has been found in patients with non-alcoholic fatty liver disease [8], type 2 diabetes [9] and cardiovascular disease burden [3]. The higher risk of mortality in the toxoplasmosis population might be partially attributed to the increased metabolic morbidity following T. gondii infection. Third, toxoplasmosis has been showed to be linked to some neuropsychiatric disorders [18] and a higher risk of traffic accidents and suicide attempts [19, 20]. The mechanism underlying the impact of toxoplasmosis on emotional and mental problems includes proinflammatory immune response to infection and endocrine and neurotransmitter dysregulation [19].

This study was based on a general population survey database and had a larger sample size and longer follow-up duration than other similar studies. Although the data are strong due to the number and quality of the data, as well as the multiple confounders we have tried to control, there are many possible explanations of this statistically significant association. Toxoplasma gondii infection has been shown to be associated with deregulated metabolic functions [21, 22] and cardiovascular disease burden [3]. The factors that increase metabolic diseases can also increase the risk for toxoplasmosis acquisition. We have attempted to control various known confounding factors, including metabolic profiles, socioeconomic status and the medical history, but even with very good multivariate models, using this type of data it is impossible to separate which factor is really influencing the outcome. Notably, the effect of toxoplasmosis on mortality was not as strong as that of diabetes, hypertension and sex. Moreover, as shown in the Fig. 2, the entire effect of toxoplasmosis attenuated when age and sex were controlled. Therefore, we are unable to rule out the possibility that the rest of the effect was mediated by other, yet unknown confounding factor(s). A well-designed prospective study involving people with metabolic risk may possibly answer this question.

Conclusions

In conclusion, in this epidemiologic follow-up study, toxoplasmosis was found to be associated with higher all-cause mortality in the general population. The results suggest that toxoplasmosis remains a neglected parasitic infection which requires more public health action.

Availability of data and materials

Publicly available datasets were analyzed in this study. The raw data used in the article are available on National Health and Nutrition Examination Survey website (https://www.cdc.gov/nchs/nhanes/index.htm).

Abbreviations

ALT:

Alanine aminotransferase

AST:

Aspartate aminotransferase

BMI:

Body mass index

BUN:

Blood urea nitrogen

CRP:

C-reactive protein

FBG:

Fasting blood glucose

HbA1c:

Glycosylated hemoglobin

HR:

Hazard ratio

NHANES III:

The Third National Health and Nutrition Examination Survey

PSM:

Propensity score matching

References

  1. Liu Q, Cao L, Zhu XQ. Major emerging and re-emerging zoonoses in China: a matter of global health and socioeconomic development for 1.3 billion. Int J Infect Dis. 2014;25:65–72.

    Article  Google Scholar 

  2. Montoya JG, Liesenfeld O. Toxoplasmosis. Lancet. 2004;363:1965–76.

    Article  CAS  Google Scholar 

  3. Flegr J, Prandota J, Sovičková M, Israili ZH. Toxoplasmosis–a global threat, correlation of latent toxoplasmosis with specific disease burden in a set of 88 countries. PLoS ONE. 2014;9:e90203.

    Article  Google Scholar 

  4. Matta SK, Rinkenberger N, Dunay IR, Sibley LD. Toxoplasma gondii infection and its implications within the central nervous system. Nat Rev Microbiol. 2021;19:467–80.

    Article  CAS  Google Scholar 

  5. Rojas-Pirela M, Medina L, Rojas MV, Liempi AI, Castillo C, Pérez-Pérez E, et al. Congenital transmission of Apicomplexan Parasites: a review. Front Microbiol. 2021;12:751648.

    Article  Google Scholar 

  6. Daher D, Shaghlil A, Sobh E, Hamie M, Hassan ME, Moumneh MB, et al. Comprehensive overview of Toxoplasma gondii-induced and associated diseases. Pathogens. 2021;10:1351.

    Article  CAS  Google Scholar 

  7. Xiao J, Prandovszky E, Kannan G, Pletnikov MV, Dickerson F, Severance EG, et al. Toxoplasma gondii: biological parameters of the connection to schizophrenia. Schizophr Bull. 2018;44:983–92.

    Article  Google Scholar 

  8. Huang J, Zhang H, Liu S, Wang M, Wan B, Velani B, et al. Is Toxoplasma gondii infection correlated with nonalcoholic fatty liver disease?—a population-based study. BMC Infect Dis. 2018;18:629.

    Article  CAS  Google Scholar 

  9. Majidiani H, Dalvand S, Daryani A, Galvan-Ramirez ML, Foroutan-Rad M. Is chronic toxoplasmosis a risk factor for diabetes mellitus? A systematic review and meta-analysis of case-control studies. Braz J Infect Dis. 2016;20:605–9.

    Article  Google Scholar 

  10. Han Y, Nie L, Ye X, Zhou Z, Huang S, Zeng C, et al. The association between Toxoplasma gondii infection and hypertensive disorders in T2DM patients: a case-control study in the Han Chinese population. Parasitol Res. 2018;117:689–95.

    Article  Google Scholar 

  11. Dickerson F, Boronow J, Stallings C, Origoni A, Yolken R. Toxoplasma gondii in individuals with schizophrenia: association with clinical and demographic factors and with mortality. Schizophr Bull. 2007;33:737–40.

    Article  Google Scholar 

  12. Roe K. The link between Toxoplasma gondii infections and higher mortality in COVID-19 patients having schizophrenia. Eur Arch Psychiatry Clin Neurosci. 2022;272:167–8.

    Article  Google Scholar 

  13. Zhou Z, Ortiz Lopez HIA, Pérez GE, Burgos LM, Farina JM, Saldarriaga C, et al. Toxoplasmosis and the heart. Curr Probl Cardiol. 2021;46:100741.

    Article  Google Scholar 

  14. Sonneville R, Magalhaes E, Meyfroidt G. Central nervous system infections in immunocompromised patients. Curr Opin Crit Care. 2017;23:128–33.

    Article  Google Scholar 

  15. Egorov AI, Converse RR, Griffin SM, Styles JN, Sams E, Hudgens E, et al. Latent Toxoplasma gondii infections are associated with elevated biomarkers of inflammation and vascular injury. BMC Infect Dis. 2021;21:188.

    Article  CAS  Google Scholar 

  16. Furman D, Campisi J, Verdin E, Carrera-Bastos P, Targ S, Franceschi C, et al. Chronic inflammation in the etiology of disease across the life span. Nat Med. 2019;25:1822–32.

    Article  CAS  Google Scholar 

  17. He JJ, Ma J, Elsheikha HM, Song HQ, Zhou DH, Zhu XQ. Proteomic profiling of mouse liver following acute Toxoplasma gondii infection. PLoS ONE. 2016;11:e0152022.

    Article  Google Scholar 

  18. Flegr J, Escudero DQ. Impaired health status and increased incidence of diseases in Toxoplasma-seropositive subjects—an explorative cross-sectional study. Parasitology. 2016;143:1974–89.

    Article  CAS  Google Scholar 

  19. Milne G, Webster JP, Walker M. Toxoplasma gondii: an underestimated threat? Trends Parasitol. 2020;36:959–69.

    Article  Google Scholar 

  20. Sutterland AL, Kuin A, Kuiper B, van Gool T, Leboyer M, Fond G, et al. Driving us mad: the association of Toxoplasma gondii with suicide attempts and traffic accidents—a systematic review and meta-analysis. Psychol Med. 2019;49:1608–23.

    Article  Google Scholar 

  21. Kankova S, Flegr J, Calda P. An elevated blood glucose level and increased incidence of gestational diabetes mellitus in pregnant women with latent toxoplasmosis. Folia Parasit (Praha). 2015;62:2015.056. https://doi.org/10.14411/fp.2015.056.

  22. Salem DA, Salem NA, Hendawy SR. Association between Toxoplasma gondii infection and metabolic syndrome in obese adolescents: a possible immune-metabolic link. Parasitol Int. 2021;83:102343.

    Article  CAS  Google Scholar 

Download references

Acknowledgements

The authors thank the National Health and Nutrition Examination Survey program for providing the data for this study.

Funding

This study was supported by the Startup Fund for scientific research of Fujian Medical University ((No. 2019QH1285), and Natural Science Foundation of Fujian Province (No. 2021J011263). All authors read and approved the final manuscript.

Author information

Authors and Affiliations

Authors

Contributions

JH, JZ and DL designed the study and wrote the article. JH and JZ designed the study and analyzed the data. BL, LL and HW were involved in manuscript preparation and data collection. DL and SL revised the final article for important intellectual content. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Dongliang Li.

Ethics declarations

Ethics approval and consent to participate

Not applicable.

Consent for publication

Not applicable.

Competing interests

The authors have nothing to declare.

Additional information

Publisher's Note

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

Supplementary Information

Additional file 1:

Table S1. Multivariate COX regression analysis grouped by the T. gondii antibody

Rights and permissions

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

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Huang, J., Zheng, J., Liu, B. et al. The association between Toxoplasma infection and mortality: the NHANES epidemiologic follow-up study. Parasites Vectors 15, 284 (2022). https://doi.org/10.1186/s13071-022-05398-1

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s13071-022-05398-1

Keywords