Direct genotyping of Toxoplasma gondii from amniotic fluids based on B1 gene polymorphism using minisequencing analysis
- Jean-Marc Costa†1Email author,
- Alexandre Alanio†2, 3,
- Sandrine Moukoury1,
- Vincent Clairet1,
- Monique Debruyne1,
- Jean-Dominique Poveda1 and
- Stéphane Bretagne†2, 3
© Costa et al.; licensee BioMed Central Ltd. 2013
Received: 23 July 2013
Accepted: 12 November 2013
Published: 19 November 2013
Because some Toxoplasma gondii genotypes may be more virulent in pregnant women, discriminating between them appears valuable. Currently, the main genotyping method is based on single copy microsatellite markers, which limit direct genotyping from amniotic fluids (AFs) to samples with a high parasitic load. We investigated whether the multicopy gene B1 could type the parasite with a higher sensitivity. To estimate the amplifiable DNA present in AFs, we first compared three different PCR assays used for Toxoplasma infection diagnosis: the P30-PCR, targeting the single copy gene P30; the B1-PCR, targeting the repeated B1 gene; and RE-PCR, targeting the repeated element.
Of the 1792 AFs analyzed between 2008 and 2011, 73 were RE-PCR positive. Of those, 49 (67.1%) were P30-PCR and B1-PCR positive, and 14 (19.2%) additional AFs were B1-PCR positive only.
All 63 BI-positive AFs (France n = 49; overseas n = 14) could be genotyped based on an analysis of eight nucleotide polymorphisms (SNPs) located within the B1 gene. Following high-resolution melting (HRM) analysis, minisequencing was carried out for each of the eight SNPs. DNA from six reference strains was included in the study, and AFs were assigned to one of the three major lineages (Types I, II, and III). In total, 26 genotypes were observed, and the hierarchical clustering distinguished two clades in lineages II (IIa, n = 30 and IIb, n = 4) and III (IIIa n = 23 and IIIb n = 6). There was an overrepresentation of overseas isolates in Clade IIb (4/4, 100%) and Clade IIIa (8/22; 36.4%) (p <0.0001), whereas medical interruption and fetal death were overrepresented in Clade IIb (2/4, 50%) and Clade IIIa (4/23, 17.4%) (p = 0.049).
Although the current genotyping system cannot pretend to replace multilocus typing, we clearly show that targeting the multicopy B1 gene yields a genotyping capacity of AFs around 20% better than when single copy targets are used. The present genotyping method also allows clear identification of genotypes of potential higher virulence.
KeywordsToxoplasma gondii P30 gene B1 gene AF487550 repeated element Genotyping Minisequencing analysis
Infection by the protozoan parasite Toxoplasma gondii can yield symptoms ranging from asymptomatic infections to miscarriages or meningoencephalitis in immunocompromised patients. In prenatal infections, fetal damages are usually linked to the time of infection during the pregnancy. The earlier on, the more likely miscarriages and severe debilitating infections become . However, genotyping has revealed genetic diversity in this worldwide parasite, and three main lineages (Types I, II, and III) , now grouped into 15 haplogroups in six major clades, have been described . Setting aside epidemiological and taxonomic studies, genotype determination is clinically relevant, as some genotypes harbor a specific virulence [4–6]. The need for a simple genotyping system, effective in clinical specimens, is clear.
Microsatellite markers have emerged as a very convenient means for direct genotyping since first reported . Up to 15 markers in a multiplex PCR format are now suggested for genotyping . However, only 62% (15/24) of the AFs sent to the French National Center for toxoplasmosis could be genotyped (http://cnrtoxoplasmose.chu-reims.fr/wp-content/uploads/2012/10/CNR-TOXOPLASMOSE-RAPPORT-ACTIVITES-2011.pdf). The small amount of DNA extracted from clinical samples for single locus typing of microsatellite markers probably explains this dearth of typable isolates. Testing microsatellite markers on samples with low T. gondii DNA content can lead to the absence of detectable peaks or peaks of low intensity not easily distinghishable from non specific PCR products . Consequently, if for any reason the non-typable samples correspond to specific genotypes, correlations drawn between symptoms and genotypes or conclusions on epidemiologic trends could be biased.
Several options for improving the rate of typability exist. One possibility is to cultivate the parasite before DNA extraction. However, cultures introduce bias, as not every isolate will grow with the same efficiency. The risk of false positives, especially in the context of a routine laboratory dealing with pregnant women, precludes a nested PCR format . Hence, our choice was to use a multicopy gene as a target. We recently investigated the possibility of genotyping T. gondii using the polymorphism of the repeated B1 gene with high-resolution melting (HRM) analysis and minisequencing (SNaPshot) using reference strains of different types . Our typing system using B1 was based on a single PCR amplification whose yield is easy to manage with an internal control  (as opposed to a multiplex PCR, where failure of one of the primer sets can have multiple causes).
We present here the results of analysis of the amniotic fluids (AFs) sent to our laboratory under suspicion of congenital toxoplasmosis. We added a formal comparison of the three targets used for Toxoplasma infection diagnosis: the single copy gene P30 , the repeated B1 gene , and the repeated element (RE) . Besides this previously unreported comparison, our aim was to find clues indicating the presence of amplifiable DNA in the sample.
Between 2008 and 2011, 1792 AFs were analyzed for the presence of T. gondii DNA. Tests were required because of proven or suspected maternal infection following systematic serological screening for anti-Toxoplasma antibodies (n = 1309; 73%), as required by French regulations , or for investigation of abnormalities revealed upon ultrasound examination (n = 483; 27%). Samples were sent from 465 different medical centers in France and overseas (French Caribbean, French Guyana, La Réunion Island, New Caledonia, and North Africa). After analysis, the remaining material was stored at -80°C, Informed consent for AF sampling, biochemical and microbiological testing, and data analysis was obtained from all women. In compliance with French regulations of prenatal diagnoses, patients provided informed consent for sample collection and analysis. Local ethics committee approval was not required, as the procedures were part of routine care. The French National Commission on Computing and Freedom was notified of the database (registration number: 1406332).
Toxoplasma gondii DNA detection
Total DNA was extracted from one ml of each AF using the Total Nucleic Acid Isolation kit (Roche Diagnostics, Meylan, France) on a MagNA Pure Compact apparatus, in accordance with the manufacturer's instructions. DNA was eluted with 100 μl of elution buffer, of which 5 μl were used for each PCR assay. Each AF was tested with three quantitative PCR (qPCR) assays targeting the P30 gene (P30-PCR; GenBank access number AY187278), the B1 gene (B1-PCR; GenBank access number AF179871), and the repeated element (RE-PCR; GenBank access number AF487550) , as previously described .
PCR reactions were carried out in a LightCycler LC480 Instrument (Roche Diagnostics, Meylan, France) in a final volume of 20 μl, using the LC480 Probes Master Kit (Roche Diagnostics, Meylan, France) with 0.5 μM of each primer and 0.25 μM of each probe. Following an initial incubation step of 1 minute at 50°C and a denaturation step of 10 minutes at 95°C, amplification was performed for 50 cycles of denaturation (95°C, 10 seconds, and ramping rate 4.4°C/second), annealing (56°C, 10 seconds, and ramping rate 2.2°C/second), and extension (72°C, 15 seconds, and ramping rate 4.4°C/second).
Toxoplasma gondii genotyping
If detectable for the B1 gene, T. gondii DNA was genotyped based on an analysis of eight nucleotide polymorphisms (SNPs) located within the B1 gene, as described elsewhere . All SNPs were studied during a single independent PCR in a final volume of 20 μl with the High Resolution Melting Master Kit (Roche Diagnostics, Meylan, France), 3 mM MgCl2, and each primer (Sigma, Paris, France) at a concentration of 0.5 μM. The reaction mixture was initially incubated for 10-min steps at 95°C. Amplification was performed for 50 cycles of denaturation (95°C for 10 sec; ramp rate, and 4.4°C/s), annealing (60°C for 10 sec; ramp rate, and 2.2°C/s), and extension (72°C for 15 sec; ramp rate, and 4.4°C/s) in a LightCycler 480 Instrument (Roche Diagnostics, Meylan, France). HRM of PCR products was then performed at 95°C for 1 min, at 40°C for 1 min, and with increasing temperature from 65°C to 95°C at a rate of 1°C/s with 25 acquisitions per °C. Following purification of the PCR products with ExoSAP-IT (USB Europe, Staufen, Germany), SNaPshot analysis was carried out using the SNaPshot Multiplex kit (Applied Biosystems, Courtaboeuf, France) for each of the eight examined SNPs, and the reactions were run on an ABI3130XL genetic analyzer and analyzed using the Genescan software. DNA from reference strains for each of the three main lineages, Type I (RH strain), Type II (B7 strain), and Type III (C5 strain), were kindly provided by Asis Khan (David Sibley laboratory), while strains from Africa (Africa 1 [RMS-2003-DJO] and Africa 2 [CCH-2004-NIA]), and the Caribbean (CCH-2005-REN) were purchased from the French National Center of Reference for toxoplasmosis.
Toxoplasma gondii genotype analysis and genotype-clinical data association
Genotypes were determined for the B1-PCR positive AFs and six reference strains based on determinations of the transition/transversion matrix  for each locus tested with SNaPshot analysis. The predominant allele (>50% of the samples) was considered the wild-type allele, and was used as the reference sequence. The polymorphic allele was then classified according to the transition/transversion matrix. In the case of heterozygosity (a double peak at the polymorphic locus), indicating a variable proportion of polymorphic repeated sequence, 50% of the corresponding transition/transversion matrix value was arbitrary attributed.
Based on the combination of the eight markers, we determined the final genotype using an arbitrary implementation. Relatedness between the different genotypes was evaluated based on hierarchical clustering (complete linkage; Euclidian distance metric) using the open-source genomic analysis software MeV v4.6.1 (The TM4 Development Group), obtained from http://mev.tm4.org. The different clades were arbitrarily determined using a threshold of 1.225 (Euclidian).
Statistical analysis and calculations
The graphs, Chi-squared tests, and ANOVA tests were performed using Prism v4.0 (GraphPAD Software, San Diego, CA). Fold changes from Ct values were calculated following this formula : Fold change (FC) = 2-ΔCt.
Results and discussion
According to the RE-PCR assay, which is used for routine testing in our diagnostic laboratory , T. gondii DNA was detected in 82/1792 (4.6%) AFs. Of the 82 infected AFs, 73 had enough sample left to be stored at -80°C, and therefore were available for simultaneous comparison of the three DNA targets. The nine AFs with no DNA left for further analysis were from different trimesters (seroconversion in first trimester n = 2, second trimester n = 5, and unknown n = 2) and of different origin (France n = 8; overseas n = 1); therefore, it is improbable that their exclusion impacted the subsequent analysis. The 73 RE-PCR positive AFs came from 65 different cities and 43 different medical centers.
As the current microsatellite typing method is based on single copy targets , the use of a multicopy gene such as B1 was expected to improve the number of typable samples. Using the B1 gene, we increased the number of typable infecting T. gondii from 67% (49/73)—expected given the ability to amplify a single copy gene—to 86% (63/73). However, the B1 gene’s sensitivity does not compete to that of RE. Of the 73 AFs available for genotyping studies, minisequencing analysis was not possible in 10/73 (13.7%) samples. These 10 non-genotypable AFs exhibited a RE-PCR Cq ranging from 32.4 to 36.3, i.e., a low quantity of T. gondii DNA. As a consequence, we expected 10-15% of non-typable samples to return positive only when using the RE-PCR assay. The 10 non-genotypable AFs were not associated with a given trimester (seroconversion in second trimester n = 3, third trimester n = 6, and unknown n = 1; p value = 0.26). These 10 non-typable AFs were all from France. This possible association of low amount of T. gondii DNA with French AFs might be due to comparatively early diagnoses due to the monthly mandatory screening of seroconversion during pregnacy instaured in France for several decades (see below).
The repartition of the French and overseas samples was significantly different in Clades IIa, IIb, IIIa and IIIb (Figure 2), with an overrepresentation of overseas isolates (p < .0001) in Clade IIb (4/4, 100%) and Clade IIIa (9/22, 36,4%). Among the 56 samples for which data was available, the repartition of medical interruption and fetal death versus no abnormality significantly differed in the four clades, with an overrepresentation (p = 0.049) of medical interruption and fetal death in Clade IIb (2/4, 50%) and Clade IIIa (4/23, 17.4%). As already reported, the European infections were mainly due to Type II strains [7, 8], whereas strains from French Carribean, French Guyana, and Reunion Island were confirmed to be more grave [4, 6].
However, seriousness depends not only on the genotype, but also on numerous other factors, such as the time of infection [1, 7]. In our study, we were able to confirm the association between medical interruption and fetal death and the first trimester of seroconversion (p = 0.0032). One might suspect that T. gondii infections during overseas pregnancies were diagnosed later, resulting in poorer outcomes. However, there was no association between the four clades and the date of seroconversion (p = 0.51). Similarly, seriousness was not associated (p = 0.33) with a geographical origin. Therefore, lower standards of care in overseas pregnancies are unlikely to bias the greater seriousness of some clades.
Our study highlights the value of using the polymorphism of the repeated B1 gene for directly genotyping T. gondii from clinical specimens without the use of any isolation or culture. Our typing system clearly individualized four clades, including the three main lineages of the parasite. Using the homogenous Type II genotype for background, we were able to individualize overseas genotypes with a higher virulence. However, we fully recognize that taxonomic studies with a single locus are limited, and therefore we did not attempt to further develop our markers in this context.
We would like to thank the clinicians for their cooperation in collecting samples. We also want to acknowledge Dr Gabrielle Cremer for the language revision of the manuscript. Our study received no external funding.
- Wallon M, Peyron F, Cornu C, Vinault S, Abrahamowicz M, Kopp CB, Binquet C: Congenital toxoplasma infection: monthly prenatal screening decreases transmission rate and improves clinical outcome at age 3 years. Clin Infect Dis. 2013, 56: 1223-1231. 10.1093/cid/cit032.View ArticlePubMedGoogle Scholar
- Saeij JP, Boyle JP, Boothroyd JC: Differences among the three major strains of Toxoplasma gondii and their specific interactions with the infected host. Trends Parasitol. 2005, 21: 476-481. 10.1016/j.pt.2005.08.001.View ArticlePubMedGoogle Scholar
- Su C, Khan A, Zhou P, Majumdar D, Ajzenberg D, Darde ML, Zhu XQ, Ajioka JW, Rosenthal BM, Dubey JP, Sibley LD: Globally diverse Toxoplasma gondii isolates comprise six major clades originating from a small number of distinct ancestral lineages. Proc Natl Acad Sci U S A. 2012, 109: 5844-5849. 10.1073/pnas.1203190109.View ArticlePubMedPubMed CentralGoogle Scholar
- Delhaes L, Ajzenberg D, Sicot B, Bourgeot P, Darde ML, Dei-Cas E, Houfflin-Debarge V: Severe congenital toxoplasmosis due to a Toxoplasma gondii strain with an atypical genotype: case report and review. Prenat Diagn. 2010, 30: 902-905. 10.1002/pd.2563.View ArticlePubMedGoogle Scholar
- Switaj K, Master A, Borkowski PK, Skrzypczak M, Wojciechowicz J, Zaborowski P: Association of ocular toxoplasmosis with type I Toxoplasma gondii strains: direct genotyping from peripheral blood samples. J Clin Microbiol. 2006, 44: 4262-4264. 10.1128/JCM.01786-06.View ArticlePubMedPubMed CentralGoogle Scholar
- Mercier A, Ajzenberg D, Devillard S, Demar MP, de Thoisy B, Bonnabau H, Collinet F, Boukhari R, Blanchet D, Simon S, et al: Human impact on genetic diversity of Toxoplasma gondii: example of the anthropized environment from French Guiana. Infect Genet Evol. 2011, 11: 1378-1387. 10.1016/j.meegid.2011.05.003.View ArticlePubMedGoogle Scholar
- Costa JM, Darde ML, Assouline B, Vidaud M, Bretagne S: Microsatellite in the beta-tubulin gene of Toxoplasma gondii as a new genetic marker for use in direct screening of amniotic fluids. J Clin Microbiol. 1997, 35: 2542-2545.PubMedPubMed CentralGoogle Scholar
- Ajzenberg D, Collinet F, Mercier A, Vignoles P, Darde ML: Genotyping of Toxoplasma gondii isolates with 15 microsatellite markers in a single multiplex PCR assay. J Clin Microbiol. 2010, 48: 4641-4645. 10.1128/JCM.01152-10.View ArticlePubMedPubMed CentralGoogle Scholar
- Bretagne S, Costa JM: Towards a nucleic acid-based diagnosis in clinical parasitology and mycology. Clin Chim Acta. 2006, 363: 221-228. 10.1016/j.cccn.2005.05.051.View ArticlePubMedGoogle Scholar
- Costa JM, Cabaret O, Moukoury S, Bretagne S: Genotyping of the protozoan pathogen Toxoplasma gondii using high-resolution melting analysis of the repeated B1 gene. J Microbiol Methods. 2011, 86: 357-363. 10.1016/j.mimet.2011.06.017.View ArticlePubMedGoogle Scholar
- Costa JM, Ernault P, Gautier E, Bretagne S: Prenatal diagnosis of congenital toxoplasmosis by duplex real-time PCR using fluorescence resonance energy transfer hybridization probes. Prenat Diagn. 2001, 21: 85-88. 10.1002/1097-0223(200102)21:2<85::AID-PD18>3.0.CO;2-1.View ArticlePubMedGoogle Scholar
- Savva D, Morris JC, Johnson JD, Holliman RE: Polymerase chain reaction for detection of Toxoplasma gondii. J Med Microbiol. 1990, 32: 25-31. 10.1099/00222615-32-1-25.View ArticlePubMedGoogle Scholar
- Burg JL, Grover CM, Pouletty P, Boothroyd JC: Direct and sensitive detection of a pathogenic protozoan, Toxoplasma gondii, by polymerase chain reaction. J Clin Microbiol. 1989, 27: 1787-1792.PubMedPubMed CentralGoogle Scholar
- Homan WL, Vercammen M, De Braekeleer J, Verschueren H: Identification of a 200- to 300-fold repetitive 529 bp DNA fragment in Toxoplasma gondii, and its use for diagnostic and quantitative PCR. Int J Parasitol. 2000, 30: 69-75. 10.1016/S0020-7519(99)00170-8.View ArticlePubMedGoogle Scholar
- Berger F, Goulet V, Le Strat Y, Desenclos JC: Toxoplasmosis among pregnant women in France: risk factors and change of prevalence between 1995 and 2003. Rev Epidemiol Sante Publique. 2009, 57: 241-248. 10.1016/j.respe.2009.03.006.View ArticlePubMedGoogle Scholar
- Reischl U, Bretagne S, Kruger D, Ernault P, Costa JM: Comparison of two DNA targets for the diagnosis of Toxoplasmosis by real-time PCR using fluorescence resonance energy transfer hybridization probes. BMC Infect Dis. 2003, 3: 7-10.1186/1471-2334-3-7.View ArticlePubMedPubMed CentralGoogle Scholar
- Costa JM, Bretagne S: Variation of B1 gene and AF146527 repeat element copy numbers according to Toxoplasma gondii strains assessed using real-time quantitative PCR. J Clin Microbiol. 2012, 50: 1452-1454. 10.1128/JCM.06514-11.View ArticlePubMedPubMed CentralGoogle Scholar
- Sung W-K: Algorithms in bioinformatics. A Practical Introduction. Edited by: Chapman & Hall. 2011, Boca Raton, FL: CRC press Taylor & Francis GroupGoogle Scholar
- Saeed AI, Sharov V, White J, Li J, Liang W, Bhagabati N, Braisted J, Klapa M, Currier T, Thiagarajan M, et al: TM4: a free, open-source system for microarray data management and analysis. Biotechniques. 2003, 34: 374-378.PubMedGoogle Scholar
- Jones CD, Okhravi N, Adamson P, Tasker S, Lightman S: Comparison of PCR detection methods for B1, P30, and 18S rDNA genes of T. gondii in aqueous humor. Invest Ophthalmol Vis Sci. 2000, 41: 634-644.PubMedGoogle Scholar
- Wahab T, Edvinsson B, Palm D, Lindh J: Comparison of the AF146527 and B1 repeated elements, two real-time PCR targets used for detection of Toxoplasma gondii. J Clin Microbiol. 2010, 48: 591-592. 10.1128/JCM.01113-09.View ArticlePubMedGoogle Scholar
- The pre-publication history for this paper can be accessed here:http://www.biomedcentral.com/1471-2334/13/552/prepub
This article is published under license to BioMed Central Ltd. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.