Introduction
Toxoplasma gondii is a zoonotic parasite that causes opportunistic infection globally. Although roughly one-third of the world's population is affected by T. gondii infection, seroprevalence of anti-T gondii immunoglobulin G (IgG) and IgM antibodies varies greatly by geographic regions and countries, with lower seroprevalence in European countries and the United States (Pappas et al., Reference Pappas, Roussos and Falagas2009). Regional differences in seroprevalence can range from 10 to 80% (Robert-Gangneux and Darde, Reference Robert-Gangneux and Darde2012), but the lack of systematic reporting of seropositivity makes it difficult to understand associated risk factors (Pappas et al., Reference Pappas, Roussos and Falagas2009). Although T. gondii infection is found in many warm-blooded animals, domestic cats are the definitive host (Chiang et al., Reference Chiang, Kuo, Chen, Yang, Kao, Ji and Fang2014; Achaw et al., Reference Achaw, Tesfa, Zeleke, Worku, Addisu, Yigzaw and Tegegne2019; Ybanez et al., Reference Ybanez, Busmeon, Viernes, Langbid, Nuevarez, Ybanez and Nishikawa2019). Numerous studies demonstrated that all cats in the world can shed T. gondii oocytes in feces, leading to environmental contamination and public health threats (Dabritz and Conrad, Reference Dabritz and Conrad2010; Berger-Schoch et al., Reference Berger-Schoch, Herrmann, Schares, Müller, Bernet, Gottstein and Frey2011; Lilly and Wortham, Reference Lilly and Wortham2013; Zulpo et al., Reference Zulpo, Sammi, Dos Santos, Sasse, Martins, Minutti, Cardim, de Barros, Navarro and Garcia2018). Human are likely to be infected via ingestion of T. gondii oocyte-contaminated soil and water, and/or consumption of raw or undercooked meat containing T. gondii tissue cysts (Elmore et al., Reference Elmore, Jones, Conrad, Patton, Lindsay and Dubey2010; Berger-Schoch et al., Reference Berger-Schoch, Herrmann, Schares, Müller, Bernet, Gottstein and Frey2011; Chiang et al., Reference Chiang, Kuo, Chen, Yang, Kao, Ji and Fang2014). Stray cats with increased exposure to pathogens are more likely to be infected by T. gondii than domestic cats (Hartmann et al., Reference Hartmann, Addie, Belak, Boucraut-Baralon, Egberink, Frymus, Gruffydd-Jones, Hosie, Lloret, Lutz, Marsilio, Mostl, Pennisi, Radford, Thiry, Truyen and Horzinek2013; Chalkowski et al., Reference Chalkowski, Wilson, Lepczyk and Zohdy2019; Khodaverdi and Razmi, Reference Khodaverdi and Razmi2019; Palerme et al., Reference Palerme, Lamperelli, Gagne, Cazlan, Zhang and Olds2019). Domestic cats, however, can transmit T. gondii to their owners, so having a cat is an independent risk factor for T. gondii seropositivity in Ethiopia and Taiwan (Chiang et al., Reference Chiang, Kuo, Chen, Yang, Kao, Ji and Fang2014; Achaw et al., Reference Achaw, Tesfa, Zeleke, Worku, Addisu, Yigzaw and Tegegne2019).
A survey revealed that among 1783 Taiwanese healthy blood donors, 161 (9.0%) donors were seropositive for anti-T. gondii IgG only and 5 (0.28%) donors were seropositive for both IgM and IgG, suggesting previous T. gondii infection (Chiang et al., Reference Chiang, Hsieh, Kuo, Chiu, Lin, Fan, Fang and Ji2012). In healthy individuals, T. gondii infection is typically mild with flu-like symptoms and lymphadenopathy (Hill et al., Reference Hill, Chirukandoth and Dubey2005), whereas central nervous system (CNS) disorders and ocular disease are also noted in congenital T. gondii infection (Chiang et al., Reference Chiang, Kuo, Chen, Yang, Kao, Ji and Fang2014). However, immunocompromised individuals such as those with HIV/AIDS (Hung et al., Reference Hung, Chen, Hsieh, Hsiao, Sheng and Chang2005), those undergoing chemotherapy or transplant recipients (Dubey and Jones, Reference Dubey and Jones2008) and neonates (Hu et al., Reference Hu, Chen, Su, Hsieh, Jeng, Liao, Su, Lin and Hsieh2006), may develop serious complications such as seizures, encephalitis or lung disease resembling tuberculosis. Notably, studies from around the world have reported associations between toxoplasma seropositivity and psychiatric disorders, such as psychosis, schizophrenia, bipolar disorder, major depression and anxiety (Lindgren et al., Reference Lindgren, Torniainen-Holm, Harkanen, Dickerson, Yolken and Suvisaari2018; Achaw et al., Reference Achaw, Tesfa, Zeleke, Worku, Addisu, Yigzaw and Tegegne2019; Palerme et al., Reference Palerme, Lamperelli, Gagne, Cazlan, Zhang and Olds2019; Stepanova et al., Reference Stepanova, Kondrashin, Sergiev, Morozova, Turbabina, Maksimova, Romanov, Kinkulkina, Lazareva and Morozov2019; Ybanez et al., Reference Ybanez, Busmeon, Viernes, Langbid, Nuevarez, Ybanez and Nishikawa2019). However, inconsistent and even controversial results were reported.
A recent cross-sectional survey found significant associations between T. gondii seropositivity and psychotic-like symptoms, but not schizophrenia, in Finland (Lindgren et al., Reference Lindgren, Torniainen-Holm, Harkanen, Dickerson, Yolken and Suvisaari2018). Although, patients diagnosed with schizophrenia had a higher incidence of latent toxoplasmosis than healthy controls in Russia (Stepanova et al., Reference Stepanova, Kondrashin, Sergiev, Morozova, Turbabina, Maksimova, Romanov, Kinkulkina, Lazareva and Morozov2019). Results of a review and meta-analysis (Sutterland et al., Reference Sutterland, Fond, Kuin, Koeter, Lutter, van Gool, Yolken, Szoke, Leboyer and de Haan2015) concluded higher odds of anti-toxoplasma IgG in patients with schizophrenia and bipolar disorder, but not in those with major depression. Clearly, there is still much to be learned about associations between T. gondii and specific mental disorders. We hypothesized that infection with T. gondii may increase the probability of developing psychiatric disorders. The comprehensive data from the National Health Insurance Research Database (NHIRD) of Taiwan may help to test our hypothesis. Therefore, the present nationwide population-based, longitudinal cohort study aimed to evaluate associations between infection with T. gondii and specific psychiatric disorders in a large population-based sample from the NHIRD.
Patients and methods
Data sources
In the present longitudinal cohort study, the data analysed were derived from the NHIRD of Taiwan, which covers 99.9% of the Taiwanese population (Wu et al., Reference Wu, Majeed and Kuo2010). The NHIRD includes all demographic and clinical information on outpatient and inpatient claims data. All clinical diagnoses and procedures are recorded based on the International Classification of Diseases, 9th Revision, Clinical Modification (ICD-9-CM).
Ethical considerations
The present study protocol was approved by the Institutional Review Board (IRB) of Tri-Service General Hospital, Taiwan (TSGHIRB No.: B-109-23). Because all patient data were anonymized in the NHIRD, informed consent was waived.
Study population
The case group consisted of patients newly diagnosed with toxoplasmosis (T gondii infection), who were identified from the NHIRD (1997–2013) based on ICD-9-CM 130. The date when toxoplasmosis was diagnosed was defined as the index date, and patients diagnosed with psychiatric disorders before the index date were excluded to ensure the temporal order between T. gondii infection and psychiatric disorders. In addition, patients younger than 18 years were excluded. In addition, the non-toxoplasmosis control group was established by propensity score matching for age, gender and index year with a 4-fold ratio of the case group (ratio 4:1), thereby reducing selection bias.
Primary outcome
In the case group, the patients with psychiatric disorders were tracked from the index date to the onset of all recorded psychiatric disorders, whereas the patients without psychiatric disorders were followed from the index date to the end of the study period. The primary outcome was the incidence of psychiatric disorders, consisting of schizophrenia (ICD-9-CM 295), bipolar disorder (ICD-9-CM 296.0, 296.4–296.8), depression (ICD-9-CM 296.2–296.3, 300.4, 311), anxiety (ICD-9-CM 300) and other psychiatric disorders (ICD-9-CM 295–312), while excluding overlapping codes. Eligible patients with the following other psychiatric disorders were identified in the present study, such as unspecified affective psychoses (296.90), other specified affective psychoses (296.99), unspecified psychosis (298.9), acute alcoholic intoxication, continuous (303.01), other and unspecified alcohol dependence (303.90), other and unspecified alcohol dependence, continuous (303.91), opioid type dependence, unspecified (304.00), alcohol abuse (305.0), alcohol abuse, unspecified (305.00), tobacco use disorder (305.1), psychogenic respiratory malfunction (306.1), tension headache (307.81), unspecified adjustment reaction (309.9) and postconcussion syndrome (310.2).
Covariates
The sociodemographic variables of the case and control groups, included age, insurance premium, season of diagnosis, place of residence, urbanization level and hospital level. The eligible patients were divided into three age groups: 18–44 years, 45–64 years and ⩾65 years. Insurance premiums in New Taiwan Dollars (NTD) were classified into three groups: <18 000, 18 000–34 999 and ⩾35 000. Regarding place of residence, Taiwan was divided into five regions, such as northern, middle, southern and eastern Taiwan, as well as outlying islands. Urbanization levels of residences were categorized into four levels: from the highest urbanization (level 1) to the lowest urbanization (level 4). The levels of hospital care were classified into three groups: local hospitals, regional hospitals and medical centres.
Statistical analysis
Categorical variables are expressed as counts and percentages; continuous variables are expressed as mean and standard deviation (mean ± s.d.). Differences in categorical variables between case and control groups were examined by χ 2 test or Fisher's exact test, whereas continuous variables were examined by t-test. Differences in the incidence of psychiatric disorders between the two groups were assessed by the Kaplan–Meier method with a log-rank test, and the result was presented as a survival curve. Univariate and multivariate Cox proportional hazards regression analyses were performed to examine the risk of developing psychiatric disorders. In addition, stratified analyses were performed based on the demographics and psychiatric disorders to explore the hazard ratios (HRs) between patients with and without toxoplasmosis. All P values were two-sided and P < 0.05 was considered statistically significant. All statistical analyses were conducted using the statistical software package SAS software version 9.4 (SAS Institute Inc., Cary, NC, USA).
Results
In the present longitudinal cohort study, 259 patients with toxoplasmosis were included in the case group, and 1036 patients without toxoplasmosis were included in the control group based on propensity score matching at a 1:4 ratio (Fig. 1). The demographic and clinical characteristics of the total population, case group and control group are summarized in Table 1. No significant differences were found in age, age group, gender, insurance premium and season of diagnosis between the case and control groups. In contrast, level of hospital care, urbanization level and place of residence were significantly different between patients with and without toxoplasmosis (all P < 0.001). Furthermore, patients with toxoplasmosis had significantly higher Charlson comorbidity index (CCI) scores than those without toxoplasmosis (P < 0.001) (Table 1).
Fig. 1. Flowchart of patient selection.
Table 1. Demographic and clinical characteristics of the total population, case group and control group
Categorical variables are presented as counts (percentages); continuous variables are presented as mean ± s.d.
Significant P values are shown in bold.
Patients with toxoplasmosis had a significantly higher incidence of psychiatric disorders than those without toxoplasmosis (P = 0.016). Thirty out of 259 (11.6%) patients in the case group were diagnosed with psychiatric disorders, and 73 out of 1036 (7.0%) patients in the control group were diagnosed with psychiatric disorders (Table 2). Particularly, the incidence of anxiety was significantly higher in patients with toxoplasmosis than in those without toxoplasmosis (P = 0.002) (Table 2). In addition, Kaplan–Meier analysis of the cumulative incidence of psychiatric disorders over 14 years of follow-up demonstrated a significant difference in incidence of psychiatric disorders between the case group and control group (log-rank P < 0.001), and a significant difference was observed as early as the first year of follow-up (Fig. 2).
Fig. 2. Kaplan–Meier for cumulative incidence of psychiatric disorders stratified by Toxoplasmosis using log-rank test.
Table 2. Incidence of psychiatric disorders in the total population, case group and control group
Values are presented as counts (percentage).
Significant P values are shown in bold.
The risk factors for psychiatric disorders in the total population were identified by univariate and multivariate Cox regression analyses (Table 3). After adjusting for age, gender, insurance premium, season of diagnosis, level of hospital care, urbanization level and CCI, multivariate Cox proportional hazards regression analyses revealed that in the total population, patients with toxoplasmosis had a significantly higher adjusted HR (aHR) for psychiatric disorders [aHR = 3.79, 95% confidence interval (CI) = 2.37, 6.07] compared to those without toxoplasmosis. Younger patients also had significantly higher HRs for psychiatric disorders (45–64 vs ⩾65: aHR = 2.23, 95% CI = 1.14, 4.36; 18–44 vs ⩾65: aHR = 5.50, 95% CI = 2.89, 10.44). Male patients had a significantly higher risk for psychiatric disorders than female patients (aHR = 2.44, 95% CI = 1.52, 3.89). Medical visits in the autumn were significantly associated with a lower risk for psychiatric disorders compared to medical visits in the spring (aHR = 0.47, 95% CI = 0.26, 0.86). Compared to seeking medical help in local hospitals, seeking medical assistance in regional hospitals and hospital centres were significantly associated with higher HRs for psychiatric disorders (regional hospital vs local hospital: aHR = 1.67, 95% CI = 1.37, 2.22; hospital centre vs local hospital: aHR = 1.82, 95% CI = 1.49, 2.37). Patients with more severe comorbidities had a significantly higher risk for psychiatric disorders (aHR = 1.08, 95% CI = 1.02, 1.15) (Table 3).
Table 3. Risks of psychiatric disorders in the total population evaluated by Cox regression analyses
aHR, adjusted hazard ratio; CI, confidence interval; NA, not applicable.
Significant P values are shown in bold.
Table 4 summarizes the risk of developing psychiatric disorders in patients with or without toxoplasmosis stratified by demographic variables. After adjusting for significant variables identified in univariate Cox regression analyses shown in Table 3, multivariate Cox regression analyses revealed that the risk of developing psychiatric disorders was significantly higher in patients with toxoplasmosis than in those without toxoplasmosis, regardless of being stratified by almost all demographic and clinical variables. Although patients who had insurance premiums less than NTD18 000 were significantly associated with an elevated risk of psychiatric disorders (aHR = 3.78, 95% CI = 2.36, 6.10) (Table 4), only one event occurred in patients with an insurance premium higher than NTD18 000.
Table 4. Multivariate Cox regression analysis for psychiatric disorders in patients with or without toxoplasmosis stratified by demographics
aHR, adjusted hazard ratio; CI, confidence interval; NA, not applicable.
Significant P values are shown in bold.
Furthermore, Cox proportional hazards regression model was applied to estimate the risk for developing each psychiatric disorder in patients with toxoplasmosis in relative to those without toxoplasmosis (Table 5). After adjusting for significant variables identified in univariate Cox regression analyses shown in Table 3, multivariate Cox regression analyses revealed that compared to the control group, the case group had significantly higher risk of bipolar disorder (aHR = 3.60, 95% CI = 2.07, 7.26), depression (aHR = 4.94, 95% CI = 2.15, 11.80) and anxiety (aHR = 5.36, 95% CI = 2.98, 25.88). In contrast, no significant differences in the HRs of developing schizophrenia and other psychiatric disorders between the two groups were observed (Table 5).
Table 5. The risk for specific psychiatric disorders between the case and control groups evaluated by Cox regression analyses
aHR, adjusted hazard ratio; CI, confidence interval.
Significant P values are shown in bold.
Discussion
In the present Taiwanese population-based cohort study, patients with toxoplasmosis had a significantly higher risk for psychiatric disorders than the non-toxoplasmosis control population, demonstrated by consistently higher aHR, regardless of stratification by demographic variables. In particular, compared to patients without toxoplasmosis, those with confirmed toxoplasmosis had significantly higher risk of bipolar disorder, depression and anxiety. However, no significant differences were found in the risk of schizophrenia and other psychiatric disorders between the patients with and without toxoplasmosis.
To our best knowledge, the present study is the first to explore associations between toxoplasmosis and several psychiatric disorders using a national population sample in Taiwan. Notably, we found no significant differences in the risk of schizophrenia between patients with and without T. gondii infection. In consistent with our findings, no significant association between T. gondii seropositivity and risk of schizophrenia was reported by two population-based studies conducted in Finland and New Zealand (Sugden et al., Reference Sugden, Moffitt, Pinto, Poulton, Williams and Caspi2016; Lindgren et al., Reference Lindgren, Torniainen-Holm, Harkanen, Dickerson, Yolken and Suvisaari2018) and a case-control study carried out in the Netherlands (de Witte et al., Reference de Witte, van Mierlo, Litjens, Klein, Bahn and Osterhaus2015).
Nevertheless, overwhelming evidence suggests otherwise. A recent large internet-based cohort study revealed that toxoplasmosis was significantly associated with schizophrenia in Czechia (Flegr and Horáček, Reference Flegr and Horáček2020), and such association has been previously documented via systemic review and meta-analysis (Torrey et al., Reference Torrey, Bartko and Yolken2012; Sutterland et al., Reference Sutterland, Fond, Kuin, Koeter, Lutter, van Gool, Yolken, Szoke, Leboyer and de Haan2015). Latent toxoplasmosis incidence, evidenced by anti-T. gondii IgM and IgG concentrations, was significantly higher in patients with schizophrenia than healthy controls in Russia (Stepanova et al., Reference Stepanova, Kondrashin, Sergiev, Morozova, Turbabina, Maksimova, Romanov, Kinkulkina, Lazareva and Morozov2019). Comparably, elevated serum levels of anti-T. gondii IgG and IgM in patients with schizophrenia in relative to those of controls were found in China (Chen et al., Reference Chen, Chen, Hou, Zheng, Yang, Ke, Hu and Tan2019). Furthermore, an analysis of the Danish Blood Donor Study suggested that T. gondii infection was a possible causative factor for schizophrenia due to the temporality of pathogen exposure (Burgdorf et al., Reference Burgdorf, Trabjerg, Pedersen, Nissen, Banasik, Pedersen, Sorensen, Nielsen, Larsen, Erikstrup, Bruun-Rasmussen, Westergaard, Thorner, Hjalgrim, Paarup, Brunak, Pedersen, Torrey, Werge, Mortensen, Yolken and Ullum2019).
Regarding the discrepancies among previous studies examining the relationship between T. gondii infection and schizophrenia, Yolken et al. (Reference Yolken, Torrey and Dickerson2017) hypothesized that the timing of patient evaluation might be one of possible explanations in a nested case-control study. They found that patients, who were diagnosed with recent onset psychosis and took a serologic test soon after psychosis manifested, had elevated odds of seropositivity for IgG antibodies against T. gondii; however, patients with established schizophrenia did not have significantly higher odds of T. gondii IgG seropositivity (Yolken et al., Reference Yolken, Torrey and Dickerson2017). Several schizophrenia medications were demonstrated to possess anti-T. gondii activity in cell culture (Jones-Brando et al., Reference Jones-Brando, Torrey and Yolken2003), which may partially explain non-reactivity of anti-T. gondii IgG antibodies in established schizophrenia. In addition, serum levels of anti-T. gondii antibodies were suggested to decline in the absence of persistent exposure to T. gondii (Rougier et al., Reference Rougier, Montoya and Peyron2017). In the present population-based study, toxoplasmosis was defined based on ICD diagnostic code without serological evidence. The ICD-9-CM coding system does not contain a specific code for recent onset psychosis. Hence, additional prospective longitudinal cohort studies are required to further investigate T. gondii infection in schizophrenia patients relying on antibody titres and T. gondii exposure history.
In the present study, patients with confirmed toxoplasmosis had significantly higher risk of bipolar disorder, and other recent studies have also reported associations between T. gondii infection and bipolar disorder (Wu et al., Reference Wu, Majeed and Kuo2010; Sutterland et al., Reference Sutterland, Fond, Kuin, Koeter, Lutter, van Gool, Yolken, Szoke, Leboyer and de Haan2015). In addition, Del Grande et al. (Reference Del Grande, Galli, Schiavi, Dell'Osso and Bruschi2017) demonstrated the link between T. gondii and bipolar disorder, emphasizing the human immune response, the presence of latent toxoplasmosis and an increased risk of bipolar disorder and suicidal/aggressive behaviour. Findings of this report suggest that re-exposure to the parasite could modify the individual's immune response and reactivate latent toxoplasmosis (Del Grande et al., Reference Del Grande, Galli, Schiavi, Dell'Osso and Bruschi2017), as has been suggested for schizophrenia (Sutterland et al., Reference Sutterland, Fond, Kuin, Koeter, Lutter, van Gool, Yolken, Szoke, Leboyer and de Haan2015). Furthermore, the potential of latent toxoplasmosis to provoke the development of bipolar disorder has been suggested (Afifi et al., Reference Afifi, Jiman-Fatani, Al-Rabia, Al-Hussainy, El Saadany and Mayah2018).
The present study showed that both depression and anxiety were more likely to occur in patients with toxoplasmosis than in those without, whereas varying results have been found in the literature. The elevated levels of T. gondii IgG antibodies were correlated with depression and anxiety in pregnant women (Groër et al., Reference Groër, Yolken, Xiao, Beckstead, Fuchs, Mohapatra, Seyfang and Postolache2011), but another study reported no association between T. gondii seropositivity and prenatal depression (Alvarado-Esquivel et al., Reference Alvarado-Esquivel, Martínez-Martínez, Sánchez-Anguiano, Hernández-Tinoco, Castillo-Orona, Salas-Martínez, Sifuentes-Álvarez, Sandoval-Carrillo, Salas-Pacheco, Liesenfeld and Antuna-Salcido2017). A community-based study found that T. gondii seropositivity was associated with anxiety, but not depression (Markovitz et al., Reference Markovitz, Simanek, Yolken, Galea, Koenen, Chen and Aiello2015), and a cohort study indicated that toxoplasmosis was associated with anxiety (Flegr and Horáček, Reference Flegr and Horáček2020). Two meta-analyses concluded that T. gondii infection was not associated with depression (Wang et al., Reference Wang, Zhang, Lei, Liu, Zhou, Liu, Wang, Yang, Zhang, Fan and Xie2014; Nayeri Chegeni et al., Reference Nayeri Chegeni, Sharif, Sarvi, Moosazadeh, Montazeri, Aghayan, Balalami, Gholami, Hosseininejad, Saberi, Anvari, Gohardehi and Daryani2019). Notably, a Danish study found that T. gondii seropositivity was associated with moderate–severe depression, but not mild–severe (Bay-Richter et al., Reference Bay-Richter, Buttenschøn, Mors, Eskelund, Budac, Kærlev and Wegener2019a; Nayeri Chegeni et al., Reference Nayeri Chegeni, Sharif, Sarvi, Moosazadeh, Montazeri, Aghayan, Balalami, Gholami, Hosseininejad, Saberi, Anvari, Gohardehi and Daryani2019); therefore, severity of depression may be one possible reason for conflicting relationships between T. gondii infection and depression reported in the literature.
Vector-borne parasitic diseases have been noted previously for associations with mental disorders. Our previous nationwide population-based cohort study (Lin et al., Reference Lin, Huang, Chung, Lin, Chen, Tsao, Chien and Chiueh2019), demonstrated a higher risk for psychiatric disorders among women with Trichomonas vaginalis infection, specifically a higher likelihood of trichomoniasis in psychiatric patients. Although cats act as a reservoir for T. gondii infections (Lilly and Wortham, Reference Lilly and Wortham2013; Zulop et al., Reference Zulpo, Sammi, Dos Santos, Sasse, Martins, Minutti, Cardim, de Barros, Navarro and Garcia2018), the vectors only pose a risk in the initial cystic stage when cats shed oocytes (Hartmann et al., Reference Hartmann, Addie, Belak, Boucraut-Baralon, Egberink, Frymus, Gruffydd-Jones, Hosie, Lloret, Lutz, Marsilio, Mostl, Pennisi, Radford, Thiry, Truyen and Horzinek2013). If sheep, cattle and pigs ingest T. gondii oocyte-contaminated soil and water, tissue cysts will develop, and humans are likely to be infected if they consume undercooked meat (Berger-Schoch et al., Reference Berger-Schoch, Herrmann, Schares, Müller, Bernet, Gottstein and Frey2011). Actual clinical disease does not occur often in T. gondii positive cats, manifesting mainly in immunosuppressed cats that may display CNS effects if disease symptoms finally do appear (Hartmann et al., Reference Hartmann, Addie, Belak, Boucraut-Baralon, Egberink, Frymus, Gruffydd-Jones, Hosie, Lloret, Lutz, Marsilio, Mostl, Pennisi, Radford, Thiry, Truyen and Horzinek2013). This pattern may present clues about the mental effects of T. gondii infection in humans.
Toxoplasma gondii has a tropism for muscle and brain tissue, where it may establish chronic infection through cyst formation (Weiss and Dubey, Reference Weiss and Dubey2009). The intracellular parasite T. gondii can alter dopamine metabolism within neurons; such metabolic fluctuation may cause behavioural changes in the host (Parlog et al., Reference Parlog, Schluter and Dunay2015). In a mouse model of T. gondii infection, degenerating neurons shown upregulated complement proteins and were surrounded by activated microglia, which are both signs of neurodegeneration typical of psychiatric disorders (Li et al., Reference Li, Severance, Viscidi, Yolken and Xiao2019). Changes in neuronal connectivity and synaptic plasticity also have been suggested as part of the mechanism responsible for psychiatric disorders upon T. gondii infection (Parlog et al., Reference Parlog, Schluter and Dunay2015).
Genetic background is also suggested to be an important influence on the response to latent toxoplasmosis (Bay-Richter et al., Reference Bay-Richter, Petersen, Liebenberg, Elfving and Wegener2019b; Torrey and Yolken, Reference Torrey and Yolken2019). Supportive of the role of gene-mediated immune response to T. gondii, Wang et al. (Reference Wang, Avramopoulos, Lori, Mulle, Conneely, Powers, Duncan, Almli, Massa, McGrath, Schwartz, Goes, Weng, Wang, Yolken, Ruczinski, Gillespie, Jovanovic, Ressler, Pulver and Pearce2019) identified two schizophrenia-related genes with significant toxoplasmosis-associated variants. Kano et al. (Reference Kano, Hodgkinson, Jones-Brando, Eastwood, Ishizuka, Niwa, Choi, Chang, Chen, Velivela, Leister, Wood, Chowdari, Ducci, Caycedo, Heinz, Newman, Cascella, Mortensen, Zandi, Dickerson, Nimgaonkar, Goldman, Harrison, Yolken and Sawa2018) found that the specific mental health-relate gene (disrupted in schizophrenia, or DISC1) activated transcription factors against T. gondii infection, and that the genotype DISC1 Phe/Phe was associated with higher titres of T. gondii antibodies in sera. Furthermore, via a comprehensive system analysis of T. gondii-infected human brains, Ngô et al. (Reference Ngô, Zhou, Lorenzi, Wang, Kim, Zhou, El Bissati, Mui, Fraczek, Rajagopala, Roberts, Henriquez, Montpetit, Blackwell, Jamieson, Wheeler, Begeman, Naranjo-Galvis, Alliey-Rodriguez, Davis, Soroceanu, Cobbs, Steindler, Boyer, Noble, Swisher, Heydemann, Rabiah, Withers, Soteropoulos, Hood and McLeod2017) identified several susceptibility genes for congenital toxoplasmosis and suggested that T. gondii modulates critical signalling pathways, which contributes to the development of epilepsy, neurodegeneration and cancer.
Strengths and limitations
One of the main strengths of the present study is the utilization of the NHIRD, which contains nationwide longitudinal data of 99.9% of Taiwan's population, and minimizes discrepancies and biases. Hence, the NHIRD is an excellent data source for population-level analysis. However, because only data from Taiwanese patients are used, the current findings may not be generalizable to other populations or other countries. The use of secondary retrospective data does not support inferences of causation and only indicates associations. The NHIRD does not provide data of individual patients' lifestyle, family history, the date when T. gondii infection occurred, and serological data, which limits the analysis of potential influencing factors. The diagnoses of toxoplasmosis and psychiatric disorders were identified in the present study based on ICD-9 codes from the NHIRD database, and it is possible that coding errors and misclassification may exist, as noted in a previous study (Kern et al., Reference Kern, Maney, Miller, Tseng, Tiwari, Rajan, Aron and Pogach2006). Although one-third of the world population was estimated to be infected by T. gondii (Montoya and Liesenfeld, Reference Montoya and Liesenfeld2004), latent T. gondii infection was not indicated in the NHIRD.
Conclusions
The present nationwide population-based analysis revealed that T. gondii infection in Taiwan significantly increases risk for developing bipolar disorder, depression and anxiety, but not schizophrenia or other psychiatric disorders. The epidemiologic role of T. gondii in the development of psychiatric disorders is an urgent public health issue that may benefit from increasing the awareness of this association and its prominent risk factors among the general public, healthcare professionals and public health agencies.
Acknowledgements
We would like to thank the National Health Insurance programme of Taiwan for making the comprehensive research database, the NHIRD, accessible. Moreover, we would like to thank the National Defense Medical Center team for support.
Financial support
This work was supported by the Tri-Service General Hospital Songshan Branch, Taiwan (TSGHSB-C108-10, TSGH-SS-E-109035) to Hsin-An Lin, and the Tri-Service General Hospital, Taiwan (TSGH-B-109010) to Wu-Chien Chien.
Conflict of interest
None.
Ethical standards
The current study protocol was reviewed and approved by the Institutional Review Board of Tri-Service General Hospital, Taiwan. Because all patient data were deidentified in the National Health Insurance Research Database (NHIRD) of Taiwan, informed consent was waived.
