- Open Access
The association between Toxoplasma infection and mortality: the NHANES epidemiologic follow-up study
Parasites & Vectors volume 15, Article number: 284 (2022)
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.
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.
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.
Toxoplasma gondii infection is associated with higher mortality in general population.
Zoonoses continue to be a major public health concern . 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  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.
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.
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.
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.
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/).
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.
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.
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).
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.
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).
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).
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 . Systemic chronic inflammation can lead to several diseases and cause disability and mortality worldwide . 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 . Higher seroprevalence of T. gondii antibody has been found in patients with non-alcoholic fatty liver disease , type 2 diabetes  and cardiovascular disease burden . 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  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 .
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 . 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.
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).
Body mass index
Blood urea nitrogen
Fasting blood glucose
- NHANES III:
The Third National Health and Nutrition Examination Survey
Propensity score matching
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.
Montoya JG, Liesenfeld O. Toxoplasmosis. Lancet. 2004;363:1965–76.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
Sonneville R, Magalhaes E, Meyfroidt G. Central nervous system infections in immunocompromised patients. Curr Opin Crit Care. 2017;23:128–33.
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.
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.
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.
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.
Milne G, Webster JP, Walker M. Toxoplasma gondii: an underestimated threat? Trends Parasitol. 2020;36:959–69.
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.
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.
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.
The authors thank the National Health and Nutrition Examination Survey program for providing the data for this study.
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.
Ethics approval and consent to participate
Consent for publication
The authors have nothing to declare.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
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
- Toxoplasma gondii