- Open Access
Stimulation and quantification of Babesia divergens gametocytogenesis
© The Author(s). 2016
- Received: 21 March 2016
- Accepted: 27 July 2016
- Published: 8 August 2016
Babesia divergens is the most common blood parasite in Europe causing babesiosis, a tick-borne malaria-like disease. Despite an increasing focus on B. divergens, especially regarding veterinary and human medicine, the sexual development of Babesia is poorly understood. Development of Babesia sexual stages in the host blood (gametocytes) plays a decisive role in parasite acquisition by the tick vector. However, the exact mechanism of gametocytogenesis is still unexplained.
Babesia divergens gametocytes are characterized by expression of bdccp1, bdccp2 and bdccp3 genes. Using previously described sequences of bdccp1, bdccp2 and bdccp3, we have established a quantitative real-time PCR (qRT-PCR) assay for detection and assessment of the efficiency of B. divergens gametocytes production in bovine blood. We analysed fluctuations in expression of bdccp genes during cultivation in vitro, as well as in cultures treated with different drugs and stimuli.
We demonstrated that all B. divergens clonal lines tested, originally derived from naturally infected cows, exhibited sexual stages. Furthermore, sexual commitment was stimulated during continuous growth of the cultures, by addition of specific stress-inducing drugs or by alternating cultivation conditions. Expression of bdccp genes was greatly reduced or even lost after long-term cultivation, suggesting possible problems in the artificial infections of ticks in feeding assays in vitro.
Our research provides insight into sexual development of B. divergens and may facilitate the development of transmission models in vitro, enabling a more detailed understanding of Babesia-tick interactions.
- Babesia divergens
- bdccp genes
Babesia are protozoan intracellular parasites infecting various vertebrates including humans. All representatives of the genus are cosmopolitan, tick-transmitted pathogens that belong to the most common blood parasites of mammals . Babesia forms a sister clade to Theileria and together they form a group referred to as Piroplasmida [1, 2]. Babesiosis caused by Babesia divergens, the most common blood parasite in Europe, is a disease in human and veterinary medicine that is occurring with increasing incidence . Babesia is evolutionarily related to Plasmodium , the agent of malaria, and both protists share many features in parasite development, such as asexual multiplication in the red blood cells (RBCs) of the vertebrate host and sexual development in the internal organs of the arthropod vector [4, 5].
Gametocytes represent essential developmental sexual stages of apicomplexan life-cycles and, in the case of Babesia, they determine the ability to infect the tick [6, 7]. The commitment from asexual growth to sexual maturation already occurs in the blood stream of the vertebrate host [7, 8]. Unlike Plasmodium, Babesia gametocytes are barely distinguishable from other asexual stages. For this reason, only laborious electron microscopy has reliably described gametocytogenesis in cultures of Babesia bigemina  or in the blood of hamsters infected with Babesia microti . The only case of gametocyte detection by light microscopy was described after stimulation of B. bigemina in vitro by addition of xanthurenic acid (XA)  or a gut homogenate from fully engorged Rhipicephalus (Boophilus) microplus ticks .
Babesia gametocytes are also poorly characterized at the molecular level. Several genes, such as heat shock protein 20 and rhoptry-associated protein 1a were believed to be transcribed in B. bigemina sexual stages. However, transcription of these genes was later found not to be exclusive for gametocytes and was also detected in other parasite stages . To date, the only molecular assay enabling specific recognition of Babesia sexual stages is based on the analysis of a highly conserved family of proteins named CCp . CCp proteins are, in general, characterized by the presence of at least one Limulus coagulation factor C (LCCL) domain [15, 16] and are often involved in cell adhesion . Gene orthologs from the highly-conserved CCp family have been identified in numerous apicomplexan parasites [16–18], including Babesia and Theileria species [14, 19]. Transcription of ccp genes was found to be restricted to gametocytes in vertebrate blood, while translation occurs in the arthropod vector to mediate gamete fertilization [14, 16, 20–24]. Based on post-genomic bioinformatic analyses of Babesia and Plasmodium genomes, three bdccp genes (bdccp1, bdccp2 and bdccp3) were thoroughly characterized and described as markers of B. divergens sexual stages . The transcripts of bdccp1, bdccp2 and bdccp3 genes were also detected in gametocytes appearing in cultures of B. divergens, B. bigemina, Babesia bovis and Theileria equi [14, 19]. Moreover, antibody targeted to BdCCp2 protein enabled visualization of B. divergens sexual stages exclusively in the midgut of Ixodes ricinus .
Here, we have established qRT-PCR conditions for the assessment of the efficiency of B. divergens gametocytes production in cultures in vitro by measuring the expression of bdccp genes. This technique is a unique tool to monitor the kinetics of B. divergens sexual stages. We analysed changes in expression of bdccp genes following variations in cultivation conditions and identified stimuli that significantly increased gametocytemia. Practical applications of our results have the potential to facilitate further detailed research in the field of Babesia-tick interactions.
Strains of B. divergens were isolated from bovine blood during the acute phases of babesioses as described earlier . 11 isolates of B. divergens from different geographical locations within France were cultivated and cloned by limited dilution . The first two digits in the description of each clone (Additional file 1: Table S1) refer to the French county of origin. Isolate Rouen 87 originated from human blood . Babesia divergens isolates were cultivated in vitro in a suspension of bovine erythrocytes obtained from a parasite-free cow (serologically negative and culture tested) as described [25, 26]. Parasitemia was monitored using the commercial Diff-Quik Stain Set (Siemens) and RBC smears.
Selection of target and reference genes, primer design, DNA extraction and PCR
List of oligonucleotides
Gene name and sequence reference
Sequence 5′-3′ and amplicon length (bp)
Sequence 5′-3′ and amplicon length (bp)
Polymorphisms in the selected reference and target genes were evaluated in 11 B. divergens clonal lines from different locations in France (listed in Additional file 1: Table S1) and compared to the B. divergens genome and other available sequences of bdccp genes. Partial gene sequences and amplicons for all target and reference genes were amplified with the sequencing primers, purified with ExoSAP-IT® (USB) and sequenced. Sequences were analysed by BioEdit v7.2.5 software.
Quantitative analysis of expression of bdccp genes
Total RNA was extracted by a combination of TRIzol® Reagent (Ambion) and NucleoSpin® RNA extraction kit (Macherey-Nagel). Briefly, 50 μl of pelleted RBCs were mixed with 200 μl of TRIzol and supplemented with 40 μl of chloroform (Sigma-Aldrich), thoroughly vortexed and centrifuged (12,000× g, 15 min, 4 °C). The aqueous phase (about 100 μl) was mixed with the same volume of 70 % ethanol and loaded onto the NucleoSpin® RNA extraction kit column. Subsequent procedures were carried out according to the manufacturer’s instructions. Evaluation of quantity and quality of RNA was performed using NanoDrop (Thermo scientific) and Experion (Bio-Rad) analyses. Residual gDNA was removed by DNase digestion with TURBO DNA-free™ Kit (Ambion) according to the manufacturer’s protocol. The absence of residual DNA was verified by lack of amplicon by PCR using qRT-PCR primers for gapdh.
Reverse transcription was performed by the SuperScript™III First-Strand Synthesis System for RT-PCR (Invitrogen) using a combination of Oligo(dT) and random hexamers according to the manufacturer’s instructions. qRT-PCR assay was performed using HOT FIREPol® EvaGreen® qPCR Mix Plus (Rox) (Solis BioDyne) in the 7300 Real-Time PCR System (Applied Biosystems). For each biological sample, three technical replicates were performed. Each assay included a standard curve generated from triplicate reactions of a 10-fold serial dilution of template. Based on standard curves, reaction efficiency and specificity were verified for each assay and each gene separately; the value of R2 > 0.98 (the correlating coefficient obtained for the standard curve) and slopes between -3.58 (reaction efficiency 90 %) and -3.10 (110 %) were accepted . For all genes, the dissociation curve analysis was performed to exclude the formation of primer-dimers and to confirm the specificity of primers. The stability of reference gene expression was tested and evaluated by comparisons of all reference genes.
The qRT-PCR results were analysed using Applied Biosystems 7300 Real-Time PCR instrument software. For analysis of the results, a comparative Ct (cycle threshold) (2-ΔΔCt) method was used [32, 33]. Mean values from technical replicates were assessed and only standard deviation (SD) values ≤ 0.5 were accepted. Target gene expression was normalized using gapdh and actin and compared using the Student t-test with Welsh’s corrections.
Analysis of expression of bdccp genes in cultures in vitro
Overview of experimental conditions and resulting effects on expression of bdccp genes
B. divergens clones
Expression of bdccp genes
Continuous culture growth
2210A G2 1802A G8
The initial parasitemia was set up at 0.1 % and expression of bdccp genes was analyzed daily for all five days post culture initiation (DPI); cultivation was performed without medium replacement.
2210A G2 6903C E2 Rouen F5
B. divergens clones 2210A G2 and 6903C E2, were continuously propagated in vitro for ≈ 1 year. Samples before and after long-term cultivation were analyzed; parasitemia was equal for all analysed samples to minimize variations in the expression of bdccp genes. Expression of bdccp genes by B. divergens clone Rouen F5 was analyzed by PCR using gDNA and cDNA.
2210A G2 Rouen G11
The range of efficient doses of both drugs was determined following parasite growth monitoring in vitro for 48 h  to select effective concentrations of drugs (imidocarbe 179.5 nM, 359 nM and 718 nM; atovaquone 10 nM, 40 nM and 75 nM). The culture without drug treatment was used as a control. The effect of drug treatment was measured 2 DPI; starting parasitemia was 2 %.
increased or decreased*
Altered cultivation temperature and air environment
XA was added at 100 μM concentration and its effect was tested after 24 h of parasites cultivation either under standard (37 °C, 5 % CO2) or altered conditions (28 °C, air). As a control, cultures without XA were used. A starting parasitemia was set up 6 % in order to reach > 10 % parasitemia level (experiment design setting taken from ).
Combination of altered cultivation and XA addition
2210A G2 Rouen G11 7101A D11
Different clonal lines were mixed in the same ratio and expression of bdccp genes was analysed in cultures cultivated for 24 h and 48 h. As a control, clones were cultivated independently; starting parasitemia was 2 %.
RBCs lysate addition
2210A G2 Rouen G11
Lysate of uninfected RBCs was added into the culture to simulate cultivation medium corresponding with 10 % parasitemia. Analyses were performed after 24 and 48 h of cultivation; the control was represented by a culture without lysate addition; starting parasitemia was 2 %.
2210A G2 Rouen G11
Hematocrit increase was simulated by doubling the quantity of RBCs in the medium and analyses were performed after 24 h and 48 h of cultivation; standard in vitro culture was used as a control; starting parasitemia was 2 %.
High parasitemia maintenance
2210A G2 Rouen G11
Analyses were performed at the starting point (0 DPI), where parasitemia was starting at 10 %, and 1 and 2 DPI. Media were changed daily.
Cultivation without FCS
2210A G2 Rouen G11
Altered cultivation conditions (cultivation in medium without FCS) were maintained for 24 h in culture with 10 % parasitemia. Analyses were performed 0 and 1 DPI; starting parasitemia was 2 %.
Statistical analyses were performed in R (version 3.2.2), a software environment for statistical computing (https://www.r-project.org/), using the Student t-test with Welsh’s correction or ANOVA followed by Tukey’s multiple comparisons test, assuming that the Bartlett test of homogeneity of variances was passed. Graphs were designed in GraphPad Prism (version 6). For graphical representations of the results and statistical analyses, mean values (± standard deviation, SD) from three biological replicates (independent experiments) were assessed.
Analysis of gene polymorphisms
Previously, polymorphisms had not been detected in 18S rDNA sequences from several B. divergens isolates . Hence sequence FJ944825 (GenBank) was used as an 18S rDNA reference gene. For other reference genes, we did not detect any polymorphisms in the b-tubulin gene and only two synonymous substitutions in gapdh and actin genes. bdccp1 was found to be highly conserved (one synonymous substitution) compared to bdccp2 (6 substitutions in the coding regions, 2 non-synonymous, resulting in 5 different sequences) and bdccp3 (7 substitutions, 2 localized in introns, all synonymous, resulting in 7 different nucleotide sequences) genes (Additional file 2: Figure S1). The qRT-PCR primers were designed only in the conserved regions.
Optimization of expression of bdccp genes
qRT-PCR was optimized as recommended by MIQE  for reference (18S, gapdh, actin, b-tubulin and 18S), as well as for target (bdccp1, bdccp2, bdccp3) genes. Standard curves of reference, target genes and qRT-PCR parameters are summarized in Additional file 3: Figure S2. Comparisons between reference genes using Ct values showed that gapdh and actin were the most stably expressed (Additional file 4: Figure S3) and these genes were therefore selected as references for further analyses. Sample normalization to gapdh or actin were consistent and no significant differences were recorded.
Expression of bdccp genes under standard cultivation conditions
Expression of bdccp genes under stress conditions
Atovaquone, another effective anti-babesial drug that induces cellular oxidative stress and is commonly used in malaria and human babesiosis treatments , caused a significant reduction in growth of B. divergens at concentrations of 40 nM (moderate inhibitory effect, t (4) = 17.88, P = 0.002) and 75 nM (complete inhibitory effect, t (4) = 22.02, P = 0.002). At a concentration of 40 nM, drug treatment resulted in significantly reduced expression of bdccp1 and bdccp2 genes: 2.2 (t (4) = 4.19, P = 0.017) and 2.0 times (t (4) = 4.55, P = 0.032), respectively, whereas at 75 nM, atovaquone significantly increased bdccp3 transcript levels (1.8 times, t (4) = −4.28, P = 0.016) (Fig. 3c, d).
A reduction in cultivation temperature from 37 to 28 °C in combination with a change in environmental conditions from 5 % CO2 to an air atmosphere, resulted in significant inhibition (t (4) = 9.84, P = 0.004) of parasite division as well as in a significant increase in expression of bdccp1 (t (4) = -5.57, P = 0.005) and bdccp2 (t (4) = -8.32, P = 0.002) genes (Fig. 4a, b). Treatment with XA, a metabolic intermediate of tryptophan degradation, has been proposed to increase the development of sexual stages in B. bigemina in vitro . In our experiments with XA treatment and cultivation at 37 °C and 5 % CO2 we identified conditions that significantly increased expression of all bdccp1, bdccp2 and bdcpp3 genes: 1.9 (t (4) = -3.97, P = 0.017), 2.4 (t (4) = -11.97, P < 0.001) and 1.6 (t (4) = -5.27, P = 0.006) times, respectively (Fig. 4a) without any inhibitory effect on culture growth (Fig. 4b). Combining XA treatment with altered cultivation conditions (28 °C, air atmosphere) resulted in significantly increased expression of bdccp1 (2.0 times, t (4) = -3.54, P = 0.029) and bdccp2 (2.4 times, t (4) = -8.39, P = 0.001) genes but culture growth was significantly inhibited (t (4) = 10.80, P = 0.002) (Fig. 4b). All other stress factors tested did not result in a significant increase in expression of bdccp genes (Table 2).
The production of gametocytes in the host blood is a prerequisite for successful parasite transmission to the arthropod vector. Plasmodium gametocytemia, which could be quantified by simple light microscopy [39, 40], was demonstrated to closely correlate with mosquito infection [41–44]. However, such a simple morphological identification is not possible for Babesia gametocytes, preventing controllable infections of ticks. Based on similarities between these two parasites , we presumed that similarly to Plasmodium, changes in the expression of Babesia sexual stage-specific bdccp genes would correlate with actual numbers of gametocytes in the total intra-erythrocytic parasite population. Using previously described sequences of bdccp1, bdccp2 and bdccp3 genes , we have developed a qRT-PCR assay to detect and quantify gametocyte densities in B. divergens cultures in vitro. Based on comparisons between the reference genes we chose actin and gapdh as references for our assays (Additional file 4: Figure S3). The 18S rDNA exhibited lower stability than actin, gapdh or b-tubulin. This result differs considerably from the generally accepted view that 18S rDNA is one of the most stably expressed genes [45, 46].
The selection of specific target and reference gene primers, universal for most of the B. divergens strains, was absolutely critical for further reliable assessment of the gametocytes production efficiency by qRT-PCR. Despite the fact that CCp proteins are presumed to be conserved among the apicomplexan parasites [16, 19, 20], no data were available about single nucleotide polymorphisms of ccp genes among various strains within one species. We demonstrated that between the 11 B. divergens clonal lines, nucleotide sequences of ccp genes varied, especially for bdccp2 and bdccp3 genes (Additional file 2: Figure S1). On the contrary, the bdccp1 gene was highly conserved. The sequences of reference genes seemed to be highly conserved. Some studies questioned the suitability of actin and gapdh reference genes because of their variabilities , but our results did not support this (Additional file 2: Figure S1) and confirmed their suitability.
The appearance of gametocytes in the blood is a crucial event that it is still not fully understood. Referring to the recent knowledge on Plasmodium, commitment towards the sexual development occurs randomly, asynchronously and is governed by the genetic and environmental factors , as demonstrated by detailed studies performed on Plasmodium (see reviews [8, 36, 37, 48–50]). To date, only one study has been dedicated to this subject in Babesia (B. bigemina) .
We tested the effects of various factors and conditions on gametocytogenesis in B. divergens cultures. Our results demonstrated the ability of B. divergens to produce gametocytes (measured by expression of bdccp genes in several bovine strains; Additional file 5: Figure S4) after a short term cultivation. On the contrary, long-term cultivation led to a significant decrease or even absence of expression of bdccp genes (Fig. 1), suggesting that these cultures had halted production of gametocytes and were probably no longer infectious for ticks. Similarly to Babesia, the disappearance of gametocytes from long-term maintained Plasmodium falciparum cultures has also been described (reviewed in ), therefore only fresh cultures with low passage numbers should be used for tick or mosquito infection studies.
The enhancement of Babesia sexual commitment was observed after several days of cultivation without medium changes, but minor variations were recorded in the bdccp genes expression of various B. divergens strains (Fig. 2, Additional file 6: Figure S5). Such phenomena could be explained by the stochastic differentiation mechanism, that was previously reported for Theileria . A rapid expansion of a Plasmodium population (intensive multiplication of asexual stages) also resulted in an increase in gametocytogenesis [34, 35]. A possible explanation of this phenomenon is the accumulation of metabolites under stress conditions, as high parasitemia or regular medium exchanges did not alter levels of bdccp transcripts (Table 2). This change is probably induced by the accumulation of metabolic waste in the blood, as an addition of a lysis solution of healthy RBCs had no significant effect (Table 2). Nevertheless, hemolysis products of both infected and healthy RBCs influenced production of gametocytes of P. falciparum in vitro as well as Plasmodium chabaudi in vivo [52, 53]. As previously shown, mixed population of Plasmodium species could result in an increase of gametocytemia and promoted more successful transmission into the vector [36, 54, 55] despite some contradictory results . We did not observe this phenomenon for B. divergens isolates (Table 2), however the choice of strains could greatly influence results, depending on their modes of interaction (neutral or synergistic instead of antagonistic).
Addition of inhibitory drugs certainly represents a stressful condition for the parasite. Numerous experiments performed on Plasmodium proved that treatment with anti-malarial drugs had an effect on the recruitment of gametocytes, both in vivo and in vitro (reviewed in [8, 48, 50, 57]). We have tested the effects of imidocarbe and atovaquone, the widely used anti-babesial drugs. Imidocarbe has been used for over 20 years as the drug of choice for the treatment and prophylaxis of animal babesiosis . The mode of action of imidocarbe still remains unclear, although disruption of polyamine metabolism or a blockage of inositol influx into parasitized cells was proposed . In our experiment, imidocarbe treatment significantly increased all three bdccp transcripts, while overall parasitemia was greatly decreased (Fig. 3). This implies that this drug either stimulated sexual commitment to the sexual pathway or has a lower impact on gametocytes as they are metabolically less active compared to asexual stages.
Atovaquone is widely used to treat babesiosis (and malaria) in humans  and causes oxidative stress in the parasite by inhibition of the mitochondrial electron transfer . This drug displayed remarkable activity against asexual stages. In gametocytes, only bdccp3 gene transcription was significantly increased (Fig. 3). It was previously demonstrated that atovaquone treatment had different effects on the various maturation stages of P. falciparum gametocytes [60, 61]. Therefore, we speculate that differences in gene expression of bdccp1 and bdccp2 compared to bdccp3 upon atovaquone application could be also related to the age of Babesia gametocytes.
Physical or chemical alterations of the parasite environment that mimic transition from the blood stream to the vector gut (temperature decrease from 37 to 28 °C, CO2 decrease from 5 % to air environment and addition of a gut homogenate from fully engorged ticks or XA) have been shown to have an effect on the Babesia sexual development [11, 12, 62]. We observed a similar stimulation of B. divergens sexual commitment after changes in the cultivation environment and/or XA addition using analysis of bdccp genes transcription. However, no apparent cumulative effect was observed when combining these stimuli. We demonstrated that XA addition into the culture under standard cultivation conditions (37 °C and 5 % CO2) significantly stimulated B. divergens sexual commitment without inhibiting parasite growth. This is in contrast to previously published results for B. bigemina, where no change in gametocyte development occurred upon XA treatment of a culture propagated under the same conditions . In mosquitoes, XA naturally produced inside the gut is able to induce gamete formation and exflagellation of Plasmodium parasites [63, 64]. As exflagellation does not occur in the Babesia life cycle, the exact effect of XA on Babesia sexual development remains to be elucidated. It is not verified yet whether XA is produced inside the tick gut. If so, one can speculate that Babesia gametocytes might be stimulated in the host blood by the tick gut contents regurgitated during the week-long feeding of the adult tick female. This hypothesis could be supported by the studies demonstrating that gametocyte development was stimulated after addition of tick gut homogenate [12, 62]. Further investigation is needed to provide an unequivocal answer.
Compared to Plasmodium, sexual development of Babesia is poorly understood. Our research provided insight into sexual development of B. divergens during either standard cultivation conditions in vitro or cultivation under stress by different stimuli. Using our newly introduced quantification assay of bdccp genes transcripts by qRT-PCR we have shown that levels of gametocytes fluctuate during B. divergens culture in vitro and identified conditions that significantly increased the transcription of bdccp genes (and thus gametocytemia). By setting these conditions we should be able to perform studies focusing on the transmission and persistence of Babesia in the tick vector using an artificial membrane feeding system of ticks . Research aimed to identify and characterize molecular mechanisms of interaction between the parasite and the tick vector could accelerate discovery of effective therapies or vaccines blocking Babesia transmission.
BLAST, basic local alignment search tool; bp, base pairs; cDNA, complementary DNA; Ct, cycle threshold; DPI, days post (culture) initiation; FCS, fetal calf serum; gapdh, glyceraldehyde 3-phosphate dehydrogenase; gDNA, genomic DNA; LCCL, Limulus coagulation factor C; qRT-PCR, quantitative real-time PCR; RBCs, red blood cells; SD, standard deviation; XA, xanthurenic acid
This work was supported by the grants No. 13-11043S, 13-27630P and 15-12006Y from the Czech Science Foundation (GA CR) and by funds from the Institut National de la Recherche Agronomique (INRA). M.J. is supported by the Dual Czech-French PhD. program (Doctorat en co-tutelle) provided by the French Institute in Prague, Czech Republic and by GAJU (Grant Agency of University of South Bohemia) 155/2013/P.
The funding agencies played no role in the design or implementation of the study, analysis or interpretation of the data, or the preparation and submission of the manuscript.
Availability of data and materials
The datasets supporting the conclusions of this article are included within the article and its Additional files.
Conceived and designed the experiments: MJ, CB, LM. Performed the experiments: MJ. Analysed the data: MJ, CB, OH, PK, LM. Contributed reagents/materials/analysis tools: CB, OH, PK, LM. Wrote the paper: MJ, OH, PK, LM. All authors read and approved the final manuscript.
The authors declare that they have no competing interests.
Consent for publication
Ethics approval and consent to participate
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
- Schnittger L, Rodriguez AE, Florin-Christensen M, Morrison DA. Babesia: a world emerging. Infect Genet Evol. 2012;12(8):1788–809.View ArticlePubMedGoogle Scholar
- Arisue N, Hashimoto T. Phylogeny and evolution of apicoplasts and apicomplexan parasites. Parasitol Int. 2014;64(3):254-9.Google Scholar
- Zintl A, Mulcahy G, Skerrett HE, Taylor SM, Gray JS. Babesia divergens, a bovine blood parasite of veterinary and zoonotic importance. Clin Microbiol Rev. 2003;16(4):622–36.View ArticlePubMedPubMed CentralGoogle Scholar
- Mehlhorn H, Shein E. The piroplasms: life cycle and sexual stages. Adv Parasitol. 1984;23:37–103.View ArticlePubMedGoogle Scholar
- Smith TG, Walliker D, Ranford-Cartwright LC. Sexual differentiation and sex determination in the Apicomplexa. Trends Parasitol. 2002;18(7):315–23.View ArticlePubMedGoogle Scholar
- Homer MJ, Aguilar-Delfin I, Telford SR, Krause PJ, Persing DH. Babesiosis. Clin Microbiol Rev. 2000;13(3):451–69.View ArticlePubMedPubMed CentralGoogle Scholar
- Chauvin A, Moreau E, Bonnet S, Plantard O, Malandrin L. Babesia and its hosts: adaptation to long-lasting interactions as a way to achieve efficient transmission. Vet Res. 2009;40(2):37.View ArticlePubMedPubMed CentralGoogle Scholar
- Dyer M, Day KP. Commitment to gametocytogenesis in Plasmodium falciparum. Parasitol Today. 2000;16(3):102–7.View ArticlePubMedGoogle Scholar
- Ribeiro MF, Patarroyo JH. Ultrastructure of Babesia bigemina gametes obtained in “in vitro” erythrocyte cultures. Vet Parasitol. 1998;76(1–2):19–25.View ArticlePubMedGoogle Scholar
- Rudzinska MA, Spielman A, Riek RF, Lewengrub SJ, Piesman J. Intraerythrocytic ‘gametocytes’ of Babesia microti and their maturation in ticks. Can J Zool. 1979;57(2):424–34.View ArticlePubMedGoogle Scholar
- Mosqueda J, Falcon A, Antonio Alvarez J, Alberto Ramos J, Oropeza-Hernandez LF, Figueroa JV. Babesia bigemina sexual stages are induced in vitro and are specifically recognized by antibodies in the midgut of infected Boophilus microplus ticks. Int J Parasitol. 2004;34(11):1229–36.View ArticlePubMedGoogle Scholar
- Gough JM, Jorgensen WK, Kemp DH. Development of tick gut forms of Babesia bigemina in vitro. J Eukaryot Microbiol. 1998;45(3):298–306.View ArticlePubMedGoogle Scholar
- Vichido R, Falcon A, Ramos JA, Alvarez A, Figueroa JV, Norimine J, et al. Expression analysis of heat shock protein 20 and rhoptry-associated protein 1a in sexual stages and kinetes of Babesia bigemina. Ann N Y Acad Sci. 2008;1149:136–40.View ArticlePubMedGoogle Scholar
- Becker CA, Malandrin L, Depoix D, Larcher T, David PH, Chauvin A, et al. Identification of three CCp genes in Babesia divergens: novel markers for sexual stages parasites. Mol Biochem Parasitol. 2010;174(1):36–43.View ArticlePubMedGoogle Scholar
- Trexler M, Bányai L, Patthy L. The LCCL module. Eur J Biochem. 2000;267(18):5751–7.View ArticlePubMedGoogle Scholar
- Dessens JT, Sinden RE, Claudianos C. LCCL proteins of apicomplexan parasites. Trends Parasitol. 2004;20(3):102–8.View ArticlePubMedGoogle Scholar
- Tosini F, Trasarti E, Pozio E. Apicomplexa genes involved in the host cell invasion: the Cpa135 protein family. Parassitologia. 2006;48(1–2):105–7.PubMedGoogle Scholar
- Templeton TJ, Enomoto S, Chen WJ, Huang CG, Lancto CA, Abrahamsen MS, Zhu G. A genome-sequence survey for Ascogregarina taiwanensis supports evolutionary affiliation but metabolic diversity between a Gregarine and Cryptosporidium. Mol Biol Evol. 2010;27(2):235–48.View ArticlePubMedGoogle Scholar
- Bastos RG, Suarez CE, Laughery JM, Johnson WC, Ueti MW, Knowles DP. Differential expression of three members of the multidomain adhesion CCp family in Babesia bigemina, Babesia bovis and Theileria equi. PLoS One. 2013;8(7):e67765.View ArticlePubMedPubMed CentralGoogle Scholar
- Pradel G, Hayton K, Aravind L, Iyer LM, Abrahamsen MS, Bonawitz A, et al. A multidomain adhesion protein family expressed in Plasmodium falciparum is essential for transmission to the mosquito. J Exp Med. 2004;199(11):1533–44.View ArticlePubMedPubMed CentralGoogle Scholar
- Lavazec C, Moreira CK, Mair GR, Waters AP, Janse CJ, Templeton TJ. Analysis of mutant Plasmodium berghei parasites lacking expression of multiple PbCCp genes. Mol Biochem Parasitol. 2009;163(1):1–7.View ArticlePubMedGoogle Scholar
- Becker CA, Malandrin L, Larcher T, Chauvin A, Bischoff E, Bonnet SI. Validation of BdCCp2 as a marker for Babesia divergens sexual stages in ticks. Exp Parasitol. 2013;133(1):51–6.View ArticlePubMedGoogle Scholar
- Scholz SM, Simon N, Lavazec C, Dude MA, Templeton TJ, Pradel G. PfCCp proteins of Plasmodium falciparum: gametocyte-specific expression and role in complement-mediated inhibition of exflagellation. Int J Parasitol. 2008;38(3–4):327–40.View ArticlePubMedGoogle Scholar
- Simon N, Scholz SM, Moreira CK, Templeton TJ, Kuehn A, Dude MA, Pradel G. Sexual stage adhesion proteins form multi-protein complexes in the malaria parasite Plasmodium falciparum. J Biol Chem. 2009;284(21):14537–46.View ArticlePubMedPubMed CentralGoogle Scholar
- Malandrin L, L’Hostis M, Chauvin A. Isolation of Babesia divergens from carrier cattle blood using in vitro culture. Vet Res. 2004;35(1):131–9.View ArticlePubMedGoogle Scholar
- Malandrin L, Jouglin M, Moreau E, Chauvin A. Individual heterogeneity in erythrocyte susceptibility to Babesia divergens is a critical factor for the outcome of experimental spleen-intact sheep infections. Vet Res. 2009;40(4):25.View ArticlePubMedPubMed CentralGoogle Scholar
- Chauvin A, Valentin A, Malandrin L, L’Hostis M. Sheep as a new experimental host for Babesia divergens. Vet Res. 2002;33(4):429–33.View ArticlePubMedGoogle Scholar
- Jackson AP, Otto TD, Darby A, Ramaprasad A, Xia D, Echaide IE, et al. The evolutionary dynamics of variant antigen genes in Babesia reveal a history of genomic innovation underlying host-parasite interaction. Nucleic Acids Res. 2014;42(11):7113–31.View ArticlePubMedPubMed CentralGoogle Scholar
- Altschul SF, Gish W, Miller W, Myers EW, Lipman DJ. Basic local alignment search tool. J Mol Biol. 1990;215(3):403–10.View ArticlePubMedGoogle Scholar
- Malandrin L, Jouglin M, Sun Y, Brisseau N, Chauvin A. Redescription of Babesia capreoli (Enigk and Friedhoff, 1962) from roe deer (Capreolus capreolus): isolation, cultivation, host specificity, molecular characterisation and differentiation from Babesia divergens. Int J Parasitol. 2010;40(3):277–84.View ArticlePubMedGoogle Scholar
- Bustin SA, Benes V, Garson JA, Hellemans J, Huggett J, Kubista M, et al. The MIQE guidelines: minimum information for publication of quantitative real-time PCR experiments. Clin Chem. 2009;55(4):611–22.View ArticlePubMedGoogle Scholar
- Livak KJ, Schmittgen TD. Analysis of relative gene expression data using real-time quantitative PCR and the 2(-Delta Delta C(T)) Method. Methods. 2001;25(4):402–8.View ArticlePubMedGoogle Scholar
- Schmittgen TD, Livak KJ. Analyzing real-time PCR data by the comparative C(T) method. Nat Protoc. 2008;3(6):1101–8.View ArticlePubMedGoogle Scholar
- Alano P, Carter R. Sexual differentiation in malaria parasites. Annu Rev Microbiol. 1990;44:429–49.View ArticlePubMedGoogle Scholar
- Lobo CA, Kumar N. Sexual differentiation and development in the malaria parasite. Parasitol Today. 1998;14(4):146–50.View ArticlePubMedGoogle Scholar
- Babiker HA, Schneider P, Reece SE. Gametocytes: insights gained during a decade of molecular monitoring. Trends Parasitol. 2008;24(11):525–30.View ArticlePubMedPubMed CentralGoogle Scholar
- Talman AM, Domarle O, McKenzie FE, Ariey F, Robert V. Gametocytogenesis: the puberty of Plasmodium falciparum. Malar J. 2004;3:24.View ArticlePubMedPubMed CentralGoogle Scholar
- Vannier E, Krause PJ. Human babesiosis. N Engl J Med. 2012;366(25):2397–407.View ArticlePubMedGoogle Scholar
- Wampfler R, Mwingira F, Javati S, Robinson L, Betuela I, Siba P, et al. Strategies for detection of Plasmodium species gametocytes. PLoS One. 2013;8(9):e76316.View ArticlePubMedPubMed CentralGoogle Scholar
- Schneider P, Reece SE, van Schaijk BC, Bousema T, Lanke KH, Meaden CS, et al. Quantification of female and male Plasmodium falciparum gametocytes by reverse transcriptase quantitative PCR. Mol Biochem Parasitol. 2015;199(1–2):29–33.PubMedGoogle Scholar
- Churcher TS, Bousema T, Walker M, Drakeley C, Schneider P, Ouédraogo AL, Basáñez MG. Predicting mosquito infection from Plasmodium falciparum gametocyte density and estimating the reservoir of infection. Elife. 2013;2:e00626.View ArticlePubMedPubMed CentralGoogle Scholar
- Da DF, Churcher TS, Yerbanga RS, Yaméogo B, Sangaré I, Ouedraogo JB, et al. Experimental study of the relationship between Plasmodium gametocyte density and infection success in mosquitoes; implications for the evaluation of malaria transmission-reducing interventions. Exp Parasitol. 2015;149:74–83.View ArticlePubMedGoogle Scholar
- Morlais I, Nsango SE, Toussile W, Abate L, Annan Z, Tchioffo MT, et al. Plasmodium falciparum mating patterns and mosquito infectivity of natural isolates of gametocytes. PLoS One. 2015;10(4):e0123777.View ArticlePubMedPubMed CentralGoogle Scholar
- Stone W, Gonçalves BP, Bousema T, Drakeley C. Assessing the infectious reservoir of falciparum malaria: past and future. Trends Parasitol. 2015;31(7):287–96.View ArticlePubMedGoogle Scholar
- Bustin SA. Quantification of mRNA using real-time reverse transcription PCR (RT-PCR): trends and problems. J Mol Endocrinol. 2002;29(1):23–39.View ArticlePubMedGoogle Scholar
- Bas A, Forsberg G, Hammarström S, Hammarström ML. Utility of the housekeeping genes 18S rRNA, beta-actin and glyceraldehyde-3-phosphate-dehydrogenase for normalization in real-time quantitative reverse transcriptase-polymerase chain reaction analysis of gene expression in human T lymphocytes. Scand J Immunol. 2004;59(6):566–73.View ArticlePubMedGoogle Scholar
- Josling GA, Llinas M. Sexual development in Plasmodium parasites: knowing when it’s time to commit. Nat Rev Microbiol. 2015;13(9):573–87.View ArticlePubMedGoogle Scholar
- Dixon MW, Thompson J, Gardiner DL, Trenholme KR. Sex in Plasmodium: a sign of commitment. Trends Parasitol. 2008;24(4):168–75.View ArticlePubMedGoogle Scholar
- Liu Z, Miao J, Cui L. Gametocytogenesis in malaria parasite: commitment, development and regulation. Future Microbiol. 2011;6(11):1351–69.View ArticlePubMedGoogle Scholar
- Bousema T, Drakeley C. Epidemiology and infectivity of Plasmodium falciparum and Plasmodium vivax gametocytes in relation to malaria control and elimination. Clin Microbiol Rev. 2011;24(2):377–410.View ArticlePubMedPubMed CentralGoogle Scholar
- Pieszko M, Weir W, Goodhead I, Kinnaird J, Shiels B. ApiAP2 factors as candidate regulators of stochastic commitment to merozoite production in Theileria annulata. PLoS Negl Trop Dis. 2015;9(8):e0003933.View ArticlePubMedPubMed CentralGoogle Scholar
- Schneweis S, Maier WA, Seitz HM. Haemolysis of infected erythrocytes--a trigger for formation of Plasmodium falciparum gametocytes? Parasitol Res. 1991;77(5):458–60.View ArticlePubMedGoogle Scholar
- Carter LM, Schneider P, Reece SE. Information use and plasticity in the reproductive decisions of malaria parasites. Malar J. 2014;13:115.View ArticlePubMedPubMed CentralGoogle Scholar
- Taylor LH, Walliker D, Read AF. Mixed-genotype infections of the rodent malaria Plasmodium chabaudi are more infectious to mosquitoes than single-genotype infections. Parasitology. 1997;115(Pt 2):121–32.View ArticlePubMedGoogle Scholar
- Bousema JT, Drakeley CJ, Mens PF, Arens T, Houben R, Omar SA, Gouagna LC, Schallig H, Sauerwein RW. Increased Plasmodium falciparum gametocyte production in mixed infections with P. malariae. Am J Trop Med Hyg. 2008;78(3):442–8.PubMedGoogle Scholar
- Nsango SE, Abate L, Thoma M, Pompon J, Fraiture M, Rademacher A, Berry A, Awono-Ambene PH, Levashina EA, Morlais I. Genetic clonality of Plasmodium falciparum affects the outcome of infection in Anopheles gambiae. Int J Parasitol. 2012;42(6):589–95.View ArticlePubMedGoogle Scholar
- Butterworth AS, Skinner-Adams TS, Gardiner DL, Trenholme KR. Plasmodium falciparum gametocytes: with a view to a kill. Parasitology. 2013;140(14):1718–34.View ArticlePubMedGoogle Scholar
- Mosqueda J, Olvera-Ramirez A, Aguilar-Tipacamu G, Canto GJ. Current advances in detection and treatment of babesiosis. Curr Med Chem. 2012;19(10):1504–18.View ArticlePubMedPubMed CentralGoogle Scholar
- Müller IB, Hyde JE. Antimalarial drugs: modes of action and mechanisms of parasite resistance. Future Microbiol. 2010;5(12):1857–73.View ArticlePubMedGoogle Scholar
- Adjalley SH, Johnston GL, Li T, Eastman RT, Ekland EH, Eappen AG, et al. Quantitative assessment of Plasmodium falciparum sexual development reveals potent transmission-blocking activity by methylene blue. Proc Natl Acad Sci USA. 2011;108(47):E1214–23.View ArticlePubMedPubMed CentralGoogle Scholar
- Fleck SL, Pudney M, Sinden RE. The effect of atovaquone (566C80) on the maturation and viability of Plasmodium falciparum gametocytes in vitro. Trans R Soc Trop Med Hyg. 1996;90(3):309–12.View ArticlePubMedGoogle Scholar
- Weber G, Friedhoff KT. Preliminary observations on the ultrastructure of suppossed sexual stages of Babesia bigemina (Piroplasmea). Z Parasitenkd. 1977;53(1):83–92.View ArticlePubMedGoogle Scholar
- Billker O, Lindo V, Panico M, Etienne AE, Paxton T, Dell A, Rogers M, Sinden RE, Morris HR. Identification of xanthurenic acid as the putative inducer of malaria development in the mosquito. Nature. 1998;392(6673):289–92.View ArticlePubMedGoogle Scholar
- Garcia GE, Wirtz RA, Barr JR, Woolfitt A, Rosenberg R. Xanthurenic acid induces gametogenesis in Plasmodium, the malaria parasite. J Biol Chem. 1998;273(20):12003–5.View ArticlePubMedGoogle Scholar
- Kröber T, Guerin PM. In vitro feeding assays for hard ticks. Trends Parasitol. 2007;23(9):445–9.View ArticlePubMedGoogle Scholar
- Malandrin L, Marchand AM, Chauvin A. Development of a microtitre-based spectrophotometric method to monitor Babesia divergens growth in vitro. J Microbiol Methods. 2004;58(3):303–12.View ArticlePubMedGoogle Scholar