Comparative virulence of Caribbean, Brazilian and European isolates of Toxoplasma gondii
Parasites & Vectors volume 12, Article number: 104 (2019)
Toxoplasma gondii is a zoonotic parasite of global importance. The outcome of infection in humans can depend on a number of factors including the infecting stage of the parasite, inoculating dose and virulence of the infecting strain. Molecular epidemiological studies have demonstrated an abundance of atypical strains of T. gondii in South America, many of which have been associated with more severe sequelae of infection. The aim of this study was to compare the virulence of T. gondii strains isolated in the Caribbean to a virulent Brazilian strain and an avirulent European strain.
One hundred and twenty Swiss CD-1 mice were split into 8 groups of 15 mice and each group was inoculated with 200 tachyzoites of one of 8 isolates, comprising ToxoDB genotypes #1, #141, #265, #13, #3 and #6. Five mice per group were euthanized at day 8 post-inoculation (p.i.) and parasite burden was determined in heart, lungs and eyes using quantitative PCR. Lungs and brain were also examined by histopathology and immunohistochemistry. The remaining 10 mice per group were part of a survival experiment to assess virulence. DNA was extracted from tachyzoites of each of the 8 T. gondii isolates and genotyped at four ROP gene loci, including ROP5, ROP16, ROP17 and ROP18 to look for association with markers of virulence.
Infection with ToxoDB genotype #13 from the Caribbean resulted in 100% of mice being euthanized which was comparative to infection with the virulent Brazilian strain (ToxoDB genotype #6). Significantly higher parasite burdens were recorded in the lungs and eyes of mice infected with ToxoDB genotypes #13 and #6. Genotyping of ROP loci revealed that the virulent Caribbean isolates had a different ROP18/ROP5 allelic profile (3/1) to the virulent Brazilian isolate (1/3); however, the avirulent Caribbean isolate (ToxoDB genotype #1) had the same ROP18/ROP5 profile as the avirulent European isolate (ToxoDB #3) (both 2/2). Caribbean isolates of intermediate virulence (ToxoDB #141 and #265) all had the same ROP18/ROP5 allelic profile (2/2).
Isolates from the Caribbean with ToxoDB genotype #13 were acutely virulent for mice and comparable to a known virulent Brazilian isolate. The ROP protein allelic profile of the virulent Caribbean and Brazilian isolates differed indicating that perhaps other factors are involved in predicting virulence. Understanding virulence is important for predicting disease outcome in humans and may also aid vaccine design as well as drug discovery.
Toxoplasma gondii is a ubiquitous intracellular parasite capable of infecting all warm-blooded animals, including humans . Felids are the only definitive hosts shedding large numbers of oocysts into the environment in their faeces. Transmission routes of the parasite to humans include ingestion of infectious oocysts directly from the environment or in contaminated food or water, ingestion of tissue cysts in undercooked/raw infected meat, and vertical transmission from mother to foetus during a primary infection. Infection in immune competent people is usually mild or asymptomatic; however, immune compromised people and congenitally infected infants can suffer severe and life-threatening clinical signs, such as encephalitis, blindness or even death .
The population structure of T. gondii was originally thought to be clonal, consisting of three dominant lineages (designated Types I, II and III) ; however, it has become apparent that the parasite has much greater genetic variability and some genotypes are distributed in distinct geographical patterns . In Europe and North America, clonal Types II and III dominate (as well as Type 12 in North America), whilst in South and Central America there is an abundance of atypical (non-clonal) strains with no clear genotypes dominating . Virulence of T. gondii has been assessed in mice and varies depending on the infecting strain, with Type I genotypes causing 100% mortality irrespective of the inoculating dose and Types II and III causing intermediate or no mortality, depending on dose . A recent study correlating virulence markers with published virulence data for over 200 isolates from all over the world noted that isolates from North America, Europe, North Africa and Asia were non-lethal to mice at low infectious doses whereas a large proportion of isolates from South America were lethal to mice .
Polymorphic rhoptry proteins (ROP) are secreted by the parasite upon invasion and have been shown to be involved in host-parasite interactions and immune evasion . Recently, ROP5 and ROP18 alleles were identified as key determinants of virulence in mice . Studies on the population genetics and epidemiology of T. gondii have suggested a link between the geographical distribution of isolates and more severe clinical signs in immune competent humans. In South America (specifically, French Guiana and Suriname), there have been reports of severe systemic toxoplasmosis resulting in death of immune competent people infected with an atypical strain of the parasite . Similarly, the prevalence, severity and risk of ocular toxoplasmosis is much greater in Brazil than in North America and Europe .
The virulence of T. gondii strains in outbred mice appears to correlate with disease manifestations in humans making the mouse a useful tool for assessing virulence and predicting possible outcomes of human infection [10, 11]. The aim of this study was to determine the virulence of previously identified atypical strains of T. gondii isolated in the Caribbean  and compare them with a known virulent Brazilian strain  and an avirulent European strain .
In vitro culture of tachyzoites
Toxoplasma gondii tachyzoites of the strains listed in Table 1 were grown in pre-cultured Vero cells (ATCC® CCL-81™) in IMDM (supplemented with 2% FCS, 200 IU/ml penicillin and 200 μg/ml streptomycin) at a ratio of one Vero cell to two T. gondii tachyzoites. On the day of inoculation, tachyzoites were enumerated using a haemocytometer and one extra vial per group (Table 1) was prepared and stored at -20 °C for DNA extraction and confirmation of genotype by PCR-RFLP.
Mouse virulence study
One hundred and twenty female outbred Swiss CD-1 mice were housed in groups of five with access to food and water ad libitum. The mice were split into 8 groups of 15 mice. Within each group, 5 mice were euthanised on day 8 post-inoculation (p.i.) to allow for comparison of parasite burden and distribution at one time-point, and 10 mice were part of a survival experiment to assess virulence of the isolates. All mice were inoculated, intraperitoneally, with 200 tachyzoites (Table 1). Although only a single dose was chosen; this was based on previous experiments demonstrating that virulent isolates resulted in 100% mortality irrespective of dose [15,16,17,18] and that avirulent and intermediately virulent isolates may also be identified at a lower dose [19, 20]. Mice were monitored and scored twice daily (at least) for signs of clinical toxoplasmosis. Any mice that reached a maximum permissible score for 2 consecutive days and were deemed unlikely to recover were humanely euthanised in accordance with the Animals (Scientific Procedures) Act 1986. Remaining mice were euthanised when they reached the end of the experiment (day 8 p.i. or day 28 p.i.). Immediately following euthanasia, mice were bled by cardiac puncture and sera isolated as previously described . One cerebral hemisphere, one lung and both eyes from each mouse were transferred to clean 2 ml tubes and stored at -80 °C for DNA extraction, while one cerebral hemisphere and one lung were stored in 10% buffered formalin for histopathology examination (described below). All sera from mice in the survival experiment (10 per group) were tested for antibodies to T. gondii using an indirect ELISA (ID Screen® Toxoplasmosis Indirect Multi-species, IDvet, Montpellier, France) according to the manufacturer’s instructions. Sera from the 5 mice per group which were euthanised at day 8 p.i. were not tested as the ELISA detects IgG and these samples would have been negative (too early in infection).
DNA extraction and quantitative PCR
One cerebral hemisphere from each mouse was homogenised in 1 ml PBS using a syringe and 18G needle followed by a 21G needle. One lung from each mouse was homogenised in 1 ml Nuclei Lysis Solution (Promega Corporation, Southampton, UK) in Precellys tubes containing CK28 ceramic beads. Both eyes from each mouse were homogenised, together, in 1 ml Nuclei Lysis Solution using an 18G needle. DNA was extracted from 400 μl of each homogenised tissue using the Wizard® genomic DNA purification protocol (Promega Corporation, Southampton, UK) . The final pellet of DNA was resuspended in 200 μl nuclease-free water and stored at -20 °C until required for PCR.
Quantitative PCR, targeting the 529 bp repeat element, was carried out in triplicate according to a previously described method  with slight modifications. The 20 μl reaction mixture consisted of 10 μl 2× Lightcycler® 480 Probes Master mix (Roche), 0.7 μM of each primer (Tox-9Fand Tox-11R), 0.1 μM of Tox-TP1, 0.2 μM of CIAC-probe, 0.02 fg of CIAC and 250 ng of template DNA in 8 μl. Sequences for primers and probes, as well as PCR conditions have been described previously [12, 22].
Genetic characterisation of T. gondii by PCR-RFLP
DNA was extracted from aliquots of mouse inocula (tachyzoites representing each isolate) using proteinase K (Promega). Briefly, 100 μg/ml proteinase K was added to each sample prior to incubation at 50 °C for 90 min, followed by 95 °C for 15 min. Samples were centrifuged at 13,000× g for 5 min at 4 °C and the supernatant (containing the DNA) was transferred to a clean 2 ml tube and stored at -20 °C until required. The ToxoDB genotype of each isolate was confirmed using a multiplex nested PCR-RFLP targeting 10 genetic markers, including SAG1, SAG2 (5′- and 3′-SAG2 and Alt.SAG2), SAG3, BTUB, GRA6, c22-8, c29-2, L358, PK1 and Apico as previously described [23, 24]. Genotypes were determined using RFLP banding profiles of reference strains RH (Type I), M4 (Type II variant; used for all markers except Apico where Type II strain Me49 was used instead) and NED (Type III) .
Genotyping of loci known to be associated with virulence (ROP5, ROP16 and ROP18) was carried out using PCR-RFLP methods as previously described [6, 24]. For ROP18, two sets of primers were used: one to amplify a repetitive sequence (DEL) in the promoter region of archetypal Type I and II alleles; and one to amplify the upstream promoter insertion sequence (UPS) of the archetypal Type III allele . As the UPS fragment is only amplified for Type III alleles no restriction digest was required. A double digest was performed on the DEL fragment and banding patterns were compared to reference strains RH and M4 .
Histopathology and immunohistochemistry
Formalin-fixed samples of brain and lung were processed routinely and embedded in paraffin wax. A formalin fixed paraffin-wax embedded liver from a mouse known to be infected with T. gondii was used as a positive control. Sections were cut at 5 µm and stained with haematoxylin and eosin (HE) for histological examination. For immunohistochemistry (IHC), serial sections were cut at 5 µm and placed on treated glass slides (Superfrost Plus, Menzel-Glaser, Braunschweig, Germany). For IHC, tissues were deparaffinised and hydrated before endogenous peroxidase was blocked with 3% hydrogen peroxide and methanol solution for 20 min at room temperature. For heat-induced epitope retrieval (HIER), slides were placed into citrate buffer (0.01M citric acid, pH6) and then autoclaved for ten minutes at 121 °C. Normal goat serum (25%) was applied to all slides for at least 30 min at room temperature. Toxoplasma gondii rabbit polyclonal antibody (Thermo-Fisher Scientific, Illinois, USA), diluted 1:600, was applied to all samples overnight at 4 °C. For the negative control, normal rabbit serum (HKV1) was used at 1:500. After 18–24 hours, goat anti-rabbit HRP-labelled polymer (DAKO EnVision+, Carpinteria, California, USA) was applied to the slides for 30 min at room temperature. AEC (3-amino-9-ethylcarbazole) chromogen solution (Vector Laboratories, Peterborough, UK) was made as instructed in the manufacturers guidelines and applied to slides for 30 min at room temperature. All tissues were counterstained with hematoxylin, then mounted with Marienfeld microscope cover slips (Lauda-Königshofen, Germany) using Shandon Consul-Mount™ Histology Formulation (Xylene based). Sections of organ samples collected from the eight groups of mice were assessed microscopically and the severity of pathology was based on the number and size of lesions observed. The presence of T. gondii was assessed by immunohistochemistry.
The virulence of the isolates was assessed by estimating survival curves using the Kaplan-Meier method based on 10 mice per group. The overall difference between survival curves was tested for statistical significance using the asymptotic log-rank test for multiple groups  based on a chi-square probability distribution with 7 degrees of freedom. This was followed by post-hoc pairwise comparisons between curves using 2-group log-rank tests, with the resulting P-values being corrected for multiple comparisons using the Benjamini-Hochberg’s procedure  to control for false discovery rate (FDR).
Parasite burden in lung, brain and eyes was measured using a quantitative PCR and the data were analysed using a linear mixed model (LMM) fitted by restricted maximum likelihood (REML) to log-transformed data, with “Group” and “Organ”, and the interaction between them, as fixed effects and “Animal” as a random effect. Post-hoc pairwise comparisons between organs and groups were conducted from the LMM estimates, with the corresponding P-values being adjusted for FDR.
Statistical test significance was assessed at the usual 5% significance level. The statistical analyses were conducted on the R system for statistical computing v3.4 .
Mouse survival experiment
During the 28-day survival experiment, only 45% (36/80) of mice survived until the end of the experiment (Fig. 1). Forty-four mice (55%) had to be euthanised due to signs of clinical toxoplasmosis. Of the 44 mice euthanised, 20 had been infected with ToxoDB genotype #13 (St. Kitts isolates TgCkStK9 and TgCkStK11; Groups 5 and 6, respectively), 10 with ToxoDB genotype #6 (Brazilian BrI isolate TgCatBr71; Group 8), 9 with ToxoDB genotype #141 (5 with St. Kitts isolate TgCkStK13 and 4 mice with St. Kitts isolate TgCkStK10; Groups 2 and 3, respectively), 3 with ToxoDB genotype #265 (St. Kitts isolate TgCkStK2; Group 4) and 2 with ToxoDB genotype #1 (M4; Group 7). Kaplan-Meier survival curves were plotted to present the survival rate of mice for each group (Fig. 1). A log-rank test demonstrated an overall statistically significant difference in survival between all groups of mice (χ2 = 54.07, df = 7, P < 0.001). Pairwise comparisons revealed that mice infected with ToxoDB genotype #13 (Groups 5 and 6) and ToxoDB genotype #6 (Group 8) had a significantly lower survival rate in comparison to all other groups (P < 0.001). Mice infected with ToxoDB genotype #1 (Group 1) had a statistically significantly higher survival rate compared to all other groups of mice (P < 0.05) except Group 7 (European M4 isolate; ToxoDB genotype #3) and Group 4 (ToxoDB genotype #265) (P = 0.019 and P = 0.094, respectively).
Distribution and burden of parasites
Parasite distribution and burden were investigated in five mice from each group euthanised at day 8 post-infection. Using the 529-bp qPCR, T. gondii DNA was detected in the lungs of 100% (40/40) of mice, in the brain tissue of 77.5% (31/40) of mice and in the eyes of 80% (32/40) of mice. At day 8 p.i., there was a statistically significantly higher mean parasite burden in the lungs compared to brain and eyes (F(2,63) = 126.71, P < 0.001; Fig. 2). The highest burden in the lungs (Fig. 2a) was evident in mice in Group 8 inoculated with the Brazilian BrI isolate (average T. gondii DNA concentration of 173.6 pg) followed by mice in Group 5 inoculated with the St. Kitts ToxoDB genotype #13 (average T. gondii DNA concentration of 93.1 pg). Mice in Group 1, inoculated with the Type II St. Kitts isolate TgCkStK12 (ToxoDB genotype #1), demonstrated the lowest lung burden with an average T. gondii DNA concentration of only 0.4 pg. Similar to the lung, the highest parasite burden in the brain (Fig. 2b) was observed in mice in Group 8 (average T. gondii DNA concentration of 7.7 pg), followed by Group 2 (average T. gondii DNA concentration of 4.2 pg) and then Group 5 (average T. gondii DNA concentration of 3.5 pg). The parasite burden in the eyes (Fig. 2c) was low compared to lung and brain. Again, the highest burden was observed in mice in Group 8 (11.3 pg), followed by Group 5 (1.5 pg) and then Group 7 (1.2 pg). When comparing parasite burden as a whole, the overall levels were generally higher in Groups 5 and 8, representing mice infected with ToxoDB genotype #13 (St. Kitts isolate TgCkStK9) and ToxoDB genotype #6 (Brazilian isolate TgCatBr71), respectively. For Group 8, the overall parasite burden was statistically significantly higher that all other groups (P < 0.001), except for Group 5 (t(32) = -1.48, P = 0.173) and marginally Groups 2 and 7 (t(32) = -2.28, P = 0.054; t(32) = 2.22, P = 0.055). The pattern of differences in mean parasite burden between groups did not change statistically significantly across organs (F(14,63) = 1.55, P = 0.120).
Genetic characterisation and markers for virulence
Analysis of each of the 8 isolates by PCR-RFLP targeting 10 genetic markers (SAG1, SAG2 (5’, 3’ and Alt.), SAG3, BTUB, GRA6, c29-2, c22-8, L358, PK1 and Apico) revealed their expected genotypes  confirming that they had not changed during in vitro culture of tachyzoites for inoculation (data not shown). Results of PCR-RFLP analysis of polymorphic loci ROP5, ROP16, ROP17 and ROP18 are given in Table 2. For ROP5 and ROP18, three different alleles were identified, with allele 3 predominating for both loci, particularly ROP18. Loci ROP16 and ROP17 were less polymorphic with just two alleles being identified, although allele 1 predominated for both. Both the virulent isolates from St. Kitts (ToxoDB genotype #13; Groups 5 and 6) displayed a different profile to the virulent isolate from Brazil (ToxoDB genotype #6; Group 8). Of the 3 isolates of intermediate virulence, two isolates (St. Kitts isolates TgCkStK13 and TgCkStK10) displayed the same allelic profile (alleles 3, 3, 1, 1 for ROP5, ROP18, ROP16 and ROP17, respectively) and the other (St. Kitts isolate TgCkStK2) differed at the ROP17 locus (allele 2 rather than allele 1).
All mice in the survival experiment which were euthanised later than day 16 p.i. (44/80) had detectable anti-T. gondii IgG antibodies in their sera (data not shown). Mice in the survival experiment which were euthanised earlier than this time point were negative by ELISA (data not shown).
Lung and brain tissue from five mice from each group euthanised at day 8 p.i. were examined. Inflammatory infiltrate and necrotic foci were present in the lungs of all mice examined but these were more severe in the mice in Group 5 (St. Kitts isolate TgCkStK9, ToxoDB genotype #13). Mice from Groups 5, 6 and 8 all showed large numbers of parasite in the lungs. Although the lungs from mice in Group 8 displayed the highest burden of parasites, the pathology observed was not as severe as in Group 5. Very low numbers of parasites and mild/no pathology were seen in the brains of infected mice, which likely reflects the early time point of sample collection (day 8 p.i.).
In the present study, two isolates from St. Kitts, which were both ToxoDB genotype #13, proved to be highly virulent for mice and showed a level of pathogenicity similar to that of the known virulent Brazilian isolate BrI (ToxoDB genotype #6). All mice infected with these isolates (Groups 5, 6 and 8) showed severe signs of clinical toxoplasmosis, with maximum permissible scores for 2 consecutive days, and 100% were euthanised by days 12, 14 and 16, respectively for Groups 5, 8 and 6. Mice in Groups 5 and 8 also had higher parasite burdens in their lungs compared to other mice and mice in Group 5 had higher inflammatory infiltrate and necrotic foci in comparison to mice from the other groups. We have previously reported the potential virulence of ToxoDB genotype #13 in a bioassay where mice were inoculated with T. gondii-infected chicken tissues and 8 out 9 mice which received this genotype had to be euthanised due to clinical toxoplasmosis . It is noteworthy that the findings in the present study are in contrast to previous studies involving ToxoDB genotype #13 which have all described the isolate as avirulent for mice [20, 29,30,31,32,33,34,35]. In six of the eight studies, virulence of the isolates was determined following a bioassay where tissues from dogs , chickens [31, 33], cats , goats [20, 35] and a howler monkey  were homogenised and inoculated into mice with no prior determination of parasite burden in the tissues. Therefore, it is possible that the inoculating dose may have been lower than that used in the present study and resulted in a less severe infection. In the present study, two different isolates of genotype #13 were used and both resulted in 100% mortality confirming the virulence of this isolate in this study. Toxoplasma gondii isolates of the same ToxoDB genotype but different levels of virulence have been described previously. Rego et al.  described two isolates (TgPgBrPI 8 and TgPgBrPI 1), both designated ToxoDB genotype #163, which resulted in 100% mortality and 0% mortality, respectively, in mice infected with 1, 10, 100 or 1000 tachyzoites. Also, isolates with the same ToxoDB genotype (#7) have been described as 100% virulent in one study  and of intermediate virulence in another  despite both studies using a titrated dose of tachyzoites.
ToxoDB genotype #13 has been isolated from a patient from the Caribbean with toxoplasmic lymphadenopathy and a history of toxoplasmic encephalitis  and phylogenetic analysis revealed that the strain clustered more closely with Type I strains than other strains . This isolate (PSP-2003-ERO) had originally been genotyped using microsatellite markers and was defined as “Caribbean 1” genotype . Caribbean 1 genotype (genotyped using microsatellite markers) has also been isolated at post mortem from an AIDS patient with toxoplasmic encephalitis and a transplant patient who died from pulmonary toxoplasmosis  and although these isolates remain to be genotyped using the panel of PCR-RFLP markers, it is highly likely that they are also ToxoDB genotype #13 . This further highlights the potential virulence of this genotype. Indeed, it has been demonstrated that T. gondii isolates containing many Type I alleles are more pathogenic and likely to cause more severe disease [11, 40]. ToxoDB genotype #13 has Type I alleles at 5 out of 10 of the PCR-RFLP markers used to define it. It is also amongst the top ten most frequently isolated genotypes in South America .
Of note in the present study, was the burden of T. gondii in the eyes of mice infected with the Brazilian BrI isolate (Group 8). These mice had a statistically significantly higher burden in comparison to most of the other groups. BrI (ToxoDB genotype #6) is the second most frequently occurring isolate in Central and South America , and severe ocular toxoplasmosis has been shown to be diagnosed more often in this part of the world. For example, ocular involvement in T. gondii-infected patients in Europe and the USA has been estimated at 2% in comparison to almost 18% in Brazil [11, 41, 42]. Furthermore, congenitally infected children in Brazil are five times more likely to develop ocular lesions than European children, and the lesions are more severe and likely to lead to visual impairment .
ToxoDB genotype #141 proved to be of intermediate virulence in the present study, resulting in 50% and 40% mortality of mice in Groups 2 and 3, respectively. This genotype was previously isolated from a fox in the USA and was shown to be acutely virulent for mice even at low doses . In contrast, in another study, this genotype was isolated from cats in St. Kitts (Caribbean) and was shown to be avirulent for mice . In the latter study, the determination of virulence was based on a bioassay involving the inoculation of mice with digested tissues from seropositive cats and not a controlled study with a defined dose of parasites (as in the former study). This could explain the discrepancy between the results and highlights the potential pitfalls in assigning virulence based on bioassays in which the dose is not controlled, as outlined above for ToxoDB genotype #13. Differences in virulence may also be a reflection of the methods used to assign genotypes. Whilst PCR-RFLP is a widely used method, it is based on restriction enzymes detecting single nucleotide polymorphisms at a limited number of loci so it is possible that isolates which appear to have the same ToxoDB genotype may display undetected polymorphisms at other loci. Perhaps multi-locus sequence typing of virulence markers would offer a more robust approach, although this is more costly and time-consuming .
Rhoptry proteins (ROP) have been shown to play an important role in host-parasite interactions. They are secreted by the rhoptry bulb of T. gondii upon invasion of the host cell and act at the parasitophorous vacuole to inhibit the localisation of immunity-regulated GTPases (IRGs) and guanylate-binding proteins (GBPs), both of which have been shown to be host resistance factors to pathogens [44, 45]. Genetic crosses between the archetypal lineages of T. gondii and the subsequent analysis of their progeny resulted in the identification of ROP18 and ROP5 as key determinants of virulence in mice [46, 47]. Two other proteins, ROP16 and ROP17, have also been shown to play a role in the modulation of host responses to T. gondii invasion and the avoidance of parasite clearance [48, 49]. In a recent study, 240 T. gondii strains from South America and Asia were analysed using PCR-RFLP targeting ROP18, ROP5, ROP16 and ROP17 and the authors concluded that allele 1 of ROP18 in combination with allele 3 of ROP5 were key determinants of virulence in mice . In the present study, ToxoDB genotype #13 (Groups 5 and 6) and ToxoDB genotype #6 (Group 8) were both associated with 100% mortality in mice; therefore, it was expected that they would have the virulent ROP profile. However, only ToxoDB genotype #6 had allele 1 of ROP18 and allele 3 of ROP5, whereas ToxoDB genotype #13 had allele 3 of ROP18 and allele 1 of ROP5. It may be the case that other genetic factors are involved in virulence or that although the virulent genotype profile applies to most virulent strains, it does not apply to all. For example, amongst the strains analysed by Shwab et al. , one was associated with 75% mortality in mice and had the same ROP18 and ROP5 alleles as ToxoDB genotype #13 in this study. Similarly, 39 out of 106 strains, accounting for 37% of those classified as lethal for mice, had allele 4 of ROP18 and allele 3 of ROP5. Also, out of 134 non-virulent isolates, only 9 had the ROP18/5 allelic combination of 3/1. Whole genome sequencing of virulent and avirulent genotype #13 isolates may help identify differences between the isolates and potentially other markers for virulence.
The isolates of intermediate virulence in this study had the same ROP profile except for one isolate which differed at the ROP17 locus. It is of note that this isolate (ToxoDB genotype #265), which had the same ROP17 allele as the avirulent isolates (ToxoDB genotypes #1 and #3), was associated with 30% mortality which is the threshold for designating an isolate as being of intermediate virulence (≥ 30% mortality) or avirulent (< 30% mortality). Perhaps there is a role for ROP17 in determining virulence, although this has not been described before. The ROP18/ROP5 allelic profile of 3/3 displayed by the isolates of intermediate virulence in this study is in accordance with previous studies where this allelic profile has been described for isolates causing varying levels of mortality . The avirulent isolates had the same ROP profile with allele 2 at all loci, which is expected in Type II/Type II variant strains and has been described before .
Whilst it appears that ROP proteins play a role in determining virulence, it is likely to also involve a number of other parasite and host genetic factors. The majority of investigative studies into the mechanisms behind the conferrence of virulence by ROP proteins have been focused on mice but recent studies have shown that neither ROP18 nor ROP5 affected the survival of T. gondii in human cells . In mice, ROP18 has been shown to phosphorylate and inactivate IFN-γ-induced IRGs thus preventing disruption of the parasitophorous vacuole and killing of the parasite; however, humans only have two IRG genes so it is unlikely that the same mechanism confers virulence in the human host .
The results of this study have demonstrated that the ToxoDB genotype #13 isolates from the Caribbean displayed the same level of virulence as the ToxoDB genotype #6 from Brazil, resulting in 100% mortality of mice. The ToxoDB genotype #1 isolate from the Caribbean was similar to the avirulent European isolate (ToxoDB genotype #3). Whilst the isolates used in this study originated from animals, they represent the strains circulating in the environments of the Caribbean, Europe and South America which may cause disease in humans. Given the more severe sequelae associated with atypical strains it is imperative to understand their level of virulence so that disease outcome may potentially be predicted. This becomes even more important with the knowledge that it is possible for previously infected humans (and animals) to become re-infected with an atypical strain [17, 51]. A greater understanding of parasite virulence not only allows for the potential to predict disease outcome, it may also aid vaccine design, a more specific treatment regime as well as the development of new drug compounds . Saraf et al.  recently called for a standardized methodology for determining T. gondii virulence in mice to allow for direct comparison between studies. Whilst this is much needed, it must also be recognised that ethical guidelines for working with laboratory animals differ between regions (particularly Europe and the Americas) and have to be taken into consideration when conducting experiments.
Dubey JP. The history of Toxoplasma gondii - the first 100 years. J Eukaryot Microbiol. 2008;55:467–75.
Tenter AM, Heckeroth AR, Weiss LM. Toxoplasma gondii: from animals to humans. Int J Parasitol. 2000;30:1217–58.
Howe DK, Sibley LD. Toxoplasma gondii comprises three clonal lineages: correlation of parasite genotype with human disease. J Infect Dis. 1995;172:1561–6.
Shwab EK, Zhu XQ, Majumdar D, Pena HF, Gennari SM, Dubey JP, Su C. Geographical patterns of Toxoplasma gondii genetic diversity revealed by multilocus PCR-RFLP genotyping. Parasitology. 2014;141:453–61.
Robert-Gangneux F, Darde ML. Epidemiology of and diagnostic strategies for toxoplasmosis. Clin Microbiol Rev. 2012;25:264–96.
Shwab EK, Jiang T, Pena HF, Gennari SM, Dubey JP, Su C. The ROP18 and ROP5 gene allele types are highly predictive of virulence in mice across globally distributed strains of Toxoplasma gondii. Int J Parasitol. 2016;46:141–6.
Hunter CA, Sibley LD. Modulation of innate immunity by Toxoplasma gondii virulence effectors. Nat Rev Microbiol. 2012;10:766–78.
Carme B, Demar M, Ajzenberg D, Darde ML. Severe acquired toxoplasmosis caused by wild cycle of Toxoplasma gondii, French Guiana. Emerg Infect Dis. 2009;15:656–8.
Grigg ME, Dubey JP, Nussenblatt RB. Ocular toxoplasmosis: lessons from Brazil. Am J Ophthalmol. 2015;159:999–1001.
Saraf P, Shwab EK, Dubey JP, Su C. On the determination of Toxoplasma gondii virulence in mice. Exp Parasitol. 2017;174:25–30.
Xiao J, Yolken RH. Strain hypothesis of Toxoplasma gondii infection on the outcome of human diseases. Acta Physiol. 2015;213:828–45.
Hamilton CM, Kelly PJ, Boey K, Corey TM, Huynh H, Metzler D, et al. Predominance of atypical genotypes of Toxoplasma gondii in free-roaming chickens in St. Kitts. West Indies. Parasit Vectors. 2017;10:104.
Pena HF, Gennari SM, Dubey JP, Su C. Population structure and mouse-virulence of Toxoplasma gondii in Brazil. Int J Parasitol. 2008;38:561–9.
Benavides J, Maley S, Pang Y, Palarea J, Eaton S, Katzer F, et al. Development of lesions and tissue distribution of parasite in lambs orally infected with sporulated oocysts of Toxoplasma gondii. Vet Parasitol. 2011;179:209–15.
Dubey JP, Verma SK, Calero-Bernal R, Cassinelli AB, Kwok OC, Van Why K, et al. Isolation and genetic characterization of Toxoplasma gondii from black bears (Ursus americanus), bobcats (Lynx rufus), and feral cats (Felis catus) from Pennsylvania. J Eukaryot Microbiol. 2015;62:410–5.
Carneiro AC, Andrade GM, Costa JG, Pinheiro BV, Vasconcelos-Santos DV, Ferreira AM, et al. Genetic characterization of Toxoplasma gondii revealed highly diverse genotypes for isolates from newborns with congenital toxoplasmosis in southeastern Brazil. J Clin Microbiol. 2013;51:901–7.
Edwards JF, Dubey JP. Toxoplasma gondii abortion storm in sheep on a Texas farm and isolation of mouse virulent atypical genotype T. gondii from an aborted lamb from a chronically infected ewe. Vet Parasitol. 2013;192:129–36.
Dubey JP, Van Why K, Verma SK, Choudhary S, Kwok OC, Khan A, et al. Genotyping Toxoplasma gondii from wildlife in Pennsylvania and identification of natural recombinants virulent to mice. Vet Parasitol. 2014;200:74–84.
Sanchez-Sanchez R, Ferre I, Regidor-Cerrillo J, Gutierrez-Exposito D, Ferrer LM, Arteche-Villasol N, et al. Virulence in mice of a Toxoplasma gondii Type II isolate does not correlate with the outcome of experimental infection in pregnant sheep. Front Cell Infect Microbiol. 2018;8:436.
Rego WMF, Costa JGL, Baraviera RCA, Pinto LV, Bessa GL, Lopes REN, Vitor RWA. Association of ROP18 and ROP5 was efficient as a marker of virulence in atypical isolates of Toxoplasma gondii obtained from pigs and goats in Piaui, Brazil. Vet Parasitol. 2017;247:19–25.
Katzer F, Canton G, Burrells A, Palarea-Albaladejo J, Horton B, Bartley PM, et al. Immunization of lambs with the S48 strain of Toxoplasma gondii reduces tissue cyst burden following oral challenge with a complete strain of the parasite. Vet Parasitol. 2014;205:46–56.
Opsteegh M, Langelaar M, Sprong H, den Hartog L, De Craeye S, Bokken G, et al. Direct detection and genotyping of Toxoplasma gondii in meat samples using magnetic capture and PCR. Int J Food Microbiol. 2010;139:193–201.
Hamilton CM, Kelly PJ, Bartley PM, Burrells A, Porco A, Metzler D, et al. Toxoplasma gondii in livestock in St. Kitts and Nevis, West Indies. Parasit Vectors. 2015;8:166.
Su C, Shwab EK, Zhou P, Zhu XQ, Dubey JP. Moving towards an integrated approach to molecular detection and identification of Toxoplasma gondii. Parasitology. 2010;137:1–11.
Burrells A, Bartley PM, Zimmer IA, Roy S, Kitchener AC, Meredith A, et al. Evidence of the three main clonal Toxoplasma gondii lineages from wild mammalian carnivores in the UK. Parasitology. 2013;140:1768–76.
Peto R, Peto J. Asymptotically efficient rank invariant test procedures. J R Stat Soc Ser A. 1972;135:185–207.
Benjamini Y, Hochberg Y. Controlling the false discovery rate: a practical and powerful approach to multiple testing. J R Stat Soc Ser B. 1995;57:289–300.
R Development Core Team. R: A language and environment for statistical computing. Vienna: R Foundation for Statistical Computing; 2016. https://www.R-project.org/
de Almeida JC, de Melo RPB, de Morais Pedrosa C, da Silva Santos M, de Barros LD, Garcia JL, et al. First isolation and RFLP genotyping of Toxoplasma gondii from crab-eating fox (Cerdocyon thous Linnaeus, 1766). Acta Trop. 2017;169:26–9.
Dubey JP, Tiwari K, Chikweto A, Deallie C, Sharma R, Thomas D, et al. Isolation and RFLP genotyping of Toxoplasma gondii from the domestic dogs (Canis familiaris) from Grenada, West Indies revealed high genetic variability. Vet Parasitol. 2013;197:623–6.
Chikweto A, Sharma RN, Tiwari KP, Verma SK, Calero-Bernal R, Jiang T, et al. Isolation and RFLP genotyping of Toxoplasma gondii in free-range chickens (Gallus domesticus) in Grenada, West Indies, revealed widespread and dominance of clonal type III parasites. J Parasitol. 2017;103:52–5.
Dubey JP, Moura L, Majumdar D, Sundar N, Velmurugan GV, Kwok OC, et al. Isolation and characterization of viable Toxoplasma gondii isolates revealed possible high frequency of mixed infection in feral cats (Felis domesticus) from St Kitts, West Indies. Parasitology. 2009;136:589–94.
de Oliveira LN, Costa Junior LM, de Melo CF, Ramos Silva JC, Bevilaqua CM, Azevedo SS, et al. Toxoplasma gondii isolates from free-range chickens from the northeast region of Brazil. J Parasitol. 2009;95:235–7.
Pena HF, Marvulo MF, Horta MC, Silva MA, Silva JC, Siqueira DB, et al. Isolation and genetic characterisation of Toxoplasma gondii from a red-handed howler monkey (Alouatta belzebul), a jaguarundi (Puma yagouaroundi), and a black-eared opossum (Didelphis aurita) from Brazil. Vet Parasitol. 2011;175:377–81.
Ragozo AM, Yai LE, Oliveira LN, Dias RA, Goncalves HC, Azevedo SS, et al. Isolation of Toxoplasma gondii from goats from Brazil. J Parasitol. 2009;95:323–6.
Howe DK, Summers BC, Sibley LD. Acute virulence in mice is associated with markers on chromosome VIII in Toxoplasma gondii. Infect Immun. 1996;64:5193–8.
Ghosn J, Paris L, Ajzenberg D, Carcelain G, Darde ML, Tubiana R, et al. Atypical toxoplasmic manifestation after discontinuation of maintenance therapy in a human immunodeficiency virus type 1-infected patient with immune recovery. Clin Infect Dis. 2003;37:e112–4.
Ajzenberg D, Banuls AL, Su C, Dumetre A, Demar M, Carme B, Darde ML. Genetic diversity, clonality and sexuality in Toxoplasma gondii. Int J Parasitol. 2004;34:1185–96.
Ajzenberg D, Yera H, Marty P, Paris L, Dalle F, Menotti J, et al. Genotype of 88 Toxoplasma gondii isolates associated with toxoplasmosis in immunocompromised patients and correlation with clinical findings. J Infect Dis. 2009;199:1155–67.
Mercier A, Devillard S, Ngoubangoye B, Bonnabau H, Banuls AL, Durand P, et al. Additional haplogroups of Toxoplasma gondii out of Africa: population structure and mouse-virulence of strains from Gabon. PLoS Negl Trop Dis. 2010;4:e876.
Glasner PD, Silveira C, Kruszon-Moran D, Martins MC, Burnier Junior M, Silveira S, et al. An unusually high prevalence of ocular toxoplasmosis in southern Brazil. Am J Ophthalmol. 1992;114:136–44.
Holland GN. Ocular toxoplasmosis: a global reassessment. Part I: epidemiology and course of disease. Am J Opthalmol. 2003;136:973–88.
Gilbert RE, Freeman K, Lago EG, Bahia-Oliveira LM, Tan HK, Wallon M, et al. Ocular sequelae of congenital toxoplasmosis in Brazil compared with Europe. PLoS Negl Trop Dis. 2008;2:e277.
Jensen KD, Camejo A, Melo MB, Cordeiro C, Julien L, Grotenbreg GM, et al. Toxoplasma gondii superinfection and virulence during secondary infection correlate with the exact ROP5/ROP18 allelic combination. mBio. 2015;6:e02280.
Kim BH, Shenoy AR, Kumar P, Bradfield CJ, MacMicking JD. IFN-inducible GTPases in host cell defense. Cell Host Microbe. 2012;12:432–44.
Saeij JP, Boyle JP, Coller S, Taylor S, Sibley LD, Brooke-Powell ET, et al. Polymorphic secreted kinases are key virulence factors in toxoplasmosis. Science. 2006;314:1780–3.
Behnke MS, Fentress SJ, Mashayekhi M, Li LX, Taylor GA, Sibley LD. The polymorphic pseudokinase ROP5 controls virulence in Toxoplasma gondii by regulating the active kinase ROP18. PLoS Pathog. 2012;8:e1002992.
Yamamoto M, Standley DM, Takashima S, Saiga H, Okuyama M, Kayama H, et al. A single polymorphic amino acid on Toxoplasma gondii kinase ROP16 determines the direct and strain-specific activation of Stat3. J Exp Med. 2009;206:2747–60.
Etheridge RD, Alaganan A, Tang K, Lou HJ, Turk BE, Sibley LD. The Toxoplasma pseudokinase ROP5 forms complexes with ROP18 and ROP17 kinases that synergize to control acute virulence in mice. Cell Host Microbe. 2014;15:537–50.
Niedelman W, Gold DA, Rosowski EE, Sprokholt JK, Lim D, Farid Arenas A, et al. The rhoptry proteins ROP18 and ROP5 mediate Toxoplasma gondii evasion of the murine, but not the human, interferon-gamma response. PLoS Pathog. 2012;8:e1002784.
Elbez-Rubinstein A, Ajzenberg D, Darde ML, Cohen R, Dumetre A, Yera H, et al. Congenital toxoplasmosis and reinfection during pregnancy: case report, strain characterization, experimental model of reinfection, and review. J Infect Dis. 2009;199:280–5.
Alday PH, Doggett JS. Drugs in development for toxoplasmosis: advances, challenges, and current status. Drug Des Devel Ther. 2017;11:273–93.
The authors would like to thank the staff of the High Security Unit at the Moredun Research Institute for performing the mouse inoculations and assisting with mouse monitoring. The authors also thank Jackie Thomson for assistance with mouse inoculations.
Funding for this study was provided by Ross University School of Veterinary Medicine, the Moredun Research Institute and the Scottish Government (Rural and Environment Science and Analytical Services Division). SO received a scholarship from the Coordination for the Improvement of Higher Education Personnel (CAPES), Brazil (PDSE, project no. 88881.131561/2016-01).
Availability of data and materials
All data generated or analysed during this study are included in this published article.
CMH designed the study and prepared the manuscript. CMH cultured the T. gondii isolates, monitored the mice and performed mouse post mortems, DNA extractions, PCRs and genotyping. LB and FC performed the histopathology and immunohistochemistry on mouse tissues and wrote relevant sections of the manuscript. SO performed DNA extractions and PCRs on some of the mouse tissues. AB and PMB assisted with mouse inoculations, mouse monitoring and post mortems. RPBM performed the ROP genotyping. JP performed the statistical analysis and contributed to the relevant sections of the manuscript. EAI, PJK and FK contributed to study design and provided comments on the manuscript. All authors read and approved the final manuscript.
Ethics approval and consent to participate
All animal procedures complied with the Animals (Scientific Procedures) Act 1986 and were approved by the Moredun Research Institute Experimental and Ethical Review Committee (E60/15).
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
Hamilton, C.M., Black, L., Oliveira, S. et al. Comparative virulence of Caribbean, Brazilian and European isolates of Toxoplasma gondii. Parasites Vectors 12, 104 (2019). https://doi.org/10.1186/s13071-019-3372-4
- Toxoplasma gondii