Rab11A regulates dense granule transport and secretion during Toxoplasma gondii invasion of host cells and parasite replication

Toxoplasma gondii possesses an armada of secreted virulent factors that enable parasite invasion and survival into host cells. These factors are contained in specific secretory organelles, the rhoptries, micronemes and dense granules that release their content upon host cell recognition. Dense granules are secreted in a constitutive manner during parasite replication and play a crucial role in modulating host metabolic and immune responses. While the molecular mechanisms triggering rhoptry and microneme release upon host cell adhesion have been well studied, constitutive secretion remains a poorly explored aspect of T. gondii vesicular trafficking. Here, we investigated the role of the small GTPase Rab11A, a known regulator of exocytosis in eukaryotic cells. Our data revealed an essential role of Rab11A in promoting the cytoskeleton driven transport of dense granules and the release of their content into the vacuolar space. Rab11A also regulates transmembrane protein trafficking and localization during parasite replication, indicating a broader role of Rab11A in cargo exocytosis at the plasma membrane. Moreover, we found that Rab11A also regulates extracellular parasite motility and adhesion to host cells. In line with these findings, MIC2 secretion was altered in Rab11A-defective parasites, which also exhibited severe morphological defects. Strikingly, by live imaging we observed a polarized accumulation of Rab11A-positive vesicles and dense granules at the apical pole of extracellular motile and invading parasites suggesting that apically polarized Rab11A-dependent delivery of cargo regulates early secretory events during parasite entry into host cells.


Introduction
Toxoplasma gondii (T. gondii) is an obligatory intracellular parasite that belongs to the phylum Apicomplexa, typified by the presence of specific apical secretory organelles called rhoptries and micronemes. Upon contact with the host cell, rhoptry (ROP) and microneme (MIC) proteins are released in order to promote parasite entry by driving the formation of a tight parasite-host cell adhesive membrane structure (termed the moving junction) [1]. ROP proteins also contribute to building the parasitophorous vacuole (PV), within which the parasite rapidly replicates. The molecular mechanisms regulating MIC exocytosis have been well studied leading to the discovery of specific parasite signaling pathways triggering their secretion upon parasite adhesion to host cells [2] [3]. Dense granules (DG) are also parasite secretory organelles critical for its survival, which release effectors modulating host immune and metabolic responses [4]. Dense granule proteins (GRA) also promote the formation of an intravacuolar nanotubular network (IVN), which interconnects parasites during intracellular replication, thereby ensuring the synchronicity of the successive divisions [5] [6] [7]. The IVN also connects the parasite to the PV membrane (PVM), presumably enhancing parasite exchanges with its host, notably for nutrient retrieval and parasite effector release into the host cytosol [8]. In contrast to micronemes and rhoptries, DGs are randomly distributed in the parasite cytosol and the mechanisms regulating their exocytosis at the parasite plasma membrane (PM) have not been elucidated. A recent study reported that in contrast to microneme exocytosis dense granule release is negatively regulated by cytosolic Ca 2+ [9]. In Metazoan, the process of exocytosis implies the active transport of secretory vesicles to the PM, and the secretion of their content into the extracellular environment or their insertion into the PM. In mammalian cells and plants, two different exocytic routes have been described: the "constitutive secretory pathway" supports the sorting of newly synthesized proteins from the endoplasmic reticulum, through the Golgi apparatus to the PM. By contrast, the "recycling pathway" targets to the cell surface internalized material that has been transported to and sorted in the perinuclear or peripheral recycling endosomes [10]. In T. gondii, DGs are considered to be the default constitutive secretory pathway based on the observation that the SAG1-GFP fusion protein (full product or truncated of its GPI anchor (SAG1ΔGPI)) is transported within DGs before being released into the vacuolar space [11]. In addition, proteins whose specific motifs targeting them to other secretory organelles have been deleted are localized in DGs [12] [13]. Yet, there is so far no evidence that transmembrane proteins navigate through the DGs to reach the parasite's PM.
The Rab GTPases belong to the Ras small G protein subfamily and operate as molecular switches that alternate between two conformational states: the GTP-bound "active" form and the GDP-bound "inactive" form [14]. Through their interactions with various effectors, such as coat components, molecular motors, and soluble NSF attachment protein receptors (SNAREs), Rab GTPases serve as multifaceted organizers of almost all membrane trafficking processes, including vesicle budding from the donor compartment, vesicle transport along cytoskeleton tracks and vesicle tethering and fusion at the acceptor membrane [14] [15] [16]. Among the Rab GTPases, Rab11 regulates both constitutive secretory and recycling pathways, thus controlling secretion at the PM [15] [16] [17]. In mammalian cells, Rab11 promotes vesicle transport via their anchoring to both, microtubule [18] and actin-based molecular motors [19]. In addition, Rab11A promotes the tethering of recycling vesicles to the PM in concert with the exocyst complex [20] [21] [22]. Through multiple interactions with effector molecules, Rab11 influences numerous cellular functions including ciliogenesis [23], cytokinesis [24], and cell migration [25] [26]. In contrast to humans, which express over 70 Rabs, T. gondii possesses a limited number of 13 Rabs that include two isoforms of Rab11: Rab11A and Rab11B [27]. In T. gondii as well as in the related apicomplexan genus Plasmodium, Rab11Adefective parasites are unable to complete cytokinesis and show marked defects in the exocytosis-assisted process that leads to proper individualization of daughter cells, otherwise posteriorly connected [28] [29] [30]. In Plasmodium, this process was suggested to be regulated by PI4K-Rab11A-mediated secretion of vesicles from the Trans-Golgi Network (TGN) to the PM [29] [30]. Here, we further explored the functions of Rab11A in T. gondii and demonstrated its key role in regulating DG exocytosis and transmembrane protein delivery at the parasite PM. We also unraveled a novel role for Rab11A in extracellular parasite adhesion and motility, thereby contributing to host cell invasion.

Rab11A localizes to dynamic cytoplasmic vesicles
To investigate T. gondii Rab11A localization, we raised a polyclonal antibody in mice, which recognized a unique protein of the expected size of 25kDa in a total extract of type I RHΔKu80 parasites ( Fig 1A). Next, we performed immunofluorescence assays (IFA) in fixed RHΔKu80 tachyzoites. Rab11A displayed distinct localizations depending on the cell cycle stage. During the G1 phase, Rab11A was localized in cytoplasmic vesicles and as previously described [28], a signal was also detected at the Golgi/Endosome-Like Compartment (ELC) region ( Fig 1B). IFA confirmed the co-distribution of Rab11A with the TGN marker TgSortilin-like Receptor (TgSORTLR) (S1 Fig). Consistent with this observation, Rab11A was found to be mostly localized adjacent to the Rab5A signal defining the early ELC, previously shown to be tightly associated with the TGN [31] (S1 Fig). During cytokinesis, the Golgi localization of Rab11A was also detected in emerging daughter cells, together with a strong enrichment of the protein at the apical tip of the growing buds, reflecting a possible Rab11A-dependent transport of newly synthesized material between these two locations ( Fig 1B and S1B Fig). Rab11A also accumulated at the basal pole of the parasite at the end of cytokinesis (Fig 1B).
In order to get further insights into the dynamic localization of Rab11A, we used the previously established transgenic ddFKBP-myc-mCherryRab11A-RHΔKu80 parasites (from here designated as mcherryRab11A-WT parasites) [28] [32]. In this strain, the expression of Rab11A fused to a mCherry tag is under the control of an N-terminal ddFKBP tag, which allows regulation of recombinant protein levels by the inducer Shield-1. Using super-resolution live imaging of parasites expressing the Inner Membrane Complex protein IMC3-YFP and mCherryRab11A-WT, we clearly observed bi-directional trajectories of Rab11A-positive vesicles between the basal and apical poles of the parasite both within the parasite cytosol (Fig 1C and S1 Movie) and along the parasite cortex delineated by the IMC3-YFP staining (Fig 1C and S2 Movie). During cytokinesis, videomicroscopy highlights the presence of Rab11A at the Golgi area of daughter cells and the transport of Rab11A-positive vesicles along the newly formed daughter bud scaffold (S3 Movie). In addition, consistent with our IFA imaging, we detected a dynamic localization of Rab11A at the basal pole of replicating parasites (S4 Movie). Interestingly, we also noticed Rab11A-positive vesicles and tubular-like structures within the residual body region (Fig 1C, RB). This region has been recently described to harbor a dense actino-myosin network that interconnects intracellular dividing tachyzoites [6] [7], suggesting that Rab11A may regulate actin-dependent material exchanges between parasites, or the dynamics of this cell-to-cell connecting network. In line with this observation, after transient expression of actin chromobodies coupled to Emerald GFP (Cb-E) that specifically label filamentous actin [7], we visualized Rab11A-positive vesicles moving along actin-positive structures at the parasite cortex (Fig 1D, upper panel and S5 Movie) or anchored to dynamic F-actin structures within the parasite cytosol ( Fig 1D, lower panel and S6 Movie). As previously observed [7], we also detected dynamic F-actin structures at the Golgi/ ELC area that co-distribute with the Rab11A signal (S5 Movie). To investigate whether Rab11A-positive vesicle movements depend on the actin cytoskeleton, we treated IMC3-YFP/ mcherryRab11A-WT tachyzoites with cytochalasin D (CD) for 30 min before recording parasites by live imaging. Depolymerizing actin filaments by CD led to the formation of cytosolic and cortical Rab11A-positive clusters that in contrast to non-treated parasites, displayed confined trajectories as illustrated by the tracking of their displacement ( Fig 1E "tracking", S7 Movie).
Collectively, these data demonstrated that Rab11A-positive vesicle movement is dependent on actin cytoskeleton activity and that Rab11A might participate in (i) vesicle budding from the TGN/ELC, (ii) cargo transport between the apical and basal poles of the parasite and (iii) material exchange between the replicating parasites via release of vesicles at the basal pole.

Rab11A-positive vesicles dynamically co-distribute with DGs
The DG-mediated secretory pathway is considered in T. gondii to be the default constitutive secretory pathway based on the observation that soluble SAG1 protein truncated of its GPI anchor (SAG1ΔGPI) is transported within DGs before being released into the vacuolar space [10] [33]. Interestingly, the dynamic motion of Rab11A-positive vesicles was similar to recently described actin and myosin F-dependent movements of DGs [33] and Rab11A is known as a regulator of exocytosis in other eukaryotic systems [15].
In order to explore dense granule dynamics in relation to Rab11A, we expressed SAG1ΔG-PI-GFP in mcherryRab11A-WT parasites. Using live imaging, we confirmed that DG content A-Western blot analysis with specific anti-Rab11A antibodies detects a unique band at 25kDa in a RHΔKU80 parasite lysate. B-Analysis of Rab11A localization in fixed RHΔKU80 parasites using antibodies recognizing Rab11A, IMC3 and ROP2/3, as indicated. Bars: 1 μm. C-Sequences of images extracted from S1 Movie and S2 Movie (left images, white frames) showing the dynamic bi-directional movement of Rab11A-positive vesicles in the cytosol (upper sequence) and along the parasite cortex (lower sequence) of mcherryRab11A-WT and IMC3-YFP expressing parasites. Tracking of vesicle trajectory is also shown. Images on the right show a zoom of the residual body (RB) region indicated by a yellow frame in the corresponding vacuoles. The arrow indicates Rab11A-positive tubular-like structures connecting the parasites Bars: 2 μm. D-Sequences of images extracted from S5 Movie and S6 Movie (left images, white frames) showing the dynamic movement of Rab11A-positive vesicles along the parasite cortex (upper sequence) and in the cytosol (lower sequence) of mcherryRab11A-WT and Cb-Emerald GFP (Cb-E) expressing parasites. Bars: 2 μm. E-Images extracted from movie S7 Movie. mcherryRab11A-WT and IMC3-YFP expressing parasites were treated with cytochalasin D (CD) for 30 min before being recorded. Rab11A-positive vesicles are detected in clusters displaying confined trajectories (right image: "tracking"). Bar: 2 μm. https://doi.org/10.1371/journal.ppat.1008106.g001

PLOS PATHOGENS
Rab11A, a novel regulator of dense granule secretion in Toxoplasma gondii was efficiently released, as illustrated by the localization of the GFP signal in the vacuolar space (Fig 2A). GFP-positive DGs detected in the parasite cytosol displayed a significant and dynamic co-distribution with mcherryRab11A-WT positive vesicles (Fig 2B and 2C). In replicating parasites, 33,7% of the DG population co-distributed over time with Rab11A-positive vesicles, while 26,1% of Rab11A-positive vesicles co-distributed with DGs. This shows that Rab11A-positive vesicles and DGs are distinct intracellular compartments that transiently interact with each other. Consistent with this notion, fluorescent signal intensity profiles indicated that GFP-positive DGs and mcherryRab11A-positive vesicles are closely apposed (Fig 2A  and 2B). This is also clearly visualized in S8 Movie (Fig 2D), in which a DG is observed docked onto a Rab11A-positive vesicle, the latter being anchored at the periphery of the parasite, and both compartments are simultaneously transported along the parasite cortex ( Fig 2D). We tracked this GFP-positive DG motion (Fig 2E and 2F; S9 Movie and S10 Movie) and fitted the recorded xy positions over time using mathematical models of "directed" or "diffusive" motion (see M&M) [34]. We confirmed that the DG trajectory 2 is consistent with "directed" motion (fitted curve, Fig 2F) characteristic of a vesicle moving along cytoskeleton tracks, in contrast to the trajectories 1 and 3, characteristic of "confined" diffusive motions [34]. This together with the observed inhibition of Rab11A-positive vesicle and DG movements upon CD treatment ( Fig 1D) [33], suggests that Rab11A promotes DG transport by mediating their anchoring along actin filaments, at least at the parasite cortex.

Rab11A promotes DG exocytosis
To assess whether Rab11A regulates DG transport, docking or the later step of fusion at the PM, we used a previously established parasite strain that over-expresses in a rapidly inducible manner an inactive GDP locked version of Rab11A fused to the mCherry fluorescent reporter (DDmCherrycmycRab11A-DN-RHΔKu80; from hereon called mCherryRab11A-DN and distinguished from mCherryRab11A-WT) [28] [32]. By WB, we confirmed that both Rab11A-WT and Rab11A-DN proteins were expressed in similar amounts after 4 h induction with Shield-1 ( Fig 3A). First, we monitored DG release in fixed Rab11A-WT and Rab11A-DN intracellular tachyzoites following gentle saponin permeabilization, which improved detection of secreted GRA proteins localized in the vacuolar space and at the PVM. To rule out any indirect effect of the previously described cytokinesis defect on DG secretion in Rab11A-DN parasites [28], we pre-treated freshly egressed extracellular tachyzoites for 1 h with Shield-1 before seeding them on a fibroblast monolayer and analyzed DG secretion 2h and 4h after parasite invasion ( Fig 3B). We observed a drastic block of GRA1 and GRA3 secretion in Rab11A-DN parasites in contrast to Rab11A-WT in which both proteins were typically released in the vacuolar space or decorated the PVM (Fig 3B and 3C). To further analyze the role of Rab11A in DG secretion, we also expressed SAG1ΔGPI-GFP in mcherryRab11A-DN parasites. In contrast to Rab11A-WT parasites, Rab11A-DN parasites were impaired in their ability to release SAGΔGPI-GFP into the PV space ( Fig 3D). Consequently, DGs were densely packed in the cytosol, which impaired reliable automatic tracking of all vesicles and therefore the quantification of the percentage of "directed" versus "diffusive" or "confined" trajectories in the total DG population. Nonetheless, DGs appeared to mostly display diffusive and confined motions (Fig 3E: tracking of 4 DG and their resulting trajectories shown in S11 Movie, which were further analyzed in Fig 3F). In particular, the accumulation of DGs observed at the altered interface between the two segregating daughter cells accounted for a local quasi static behavior as illustrated by their confined trajectories ( Fig  3F: trajectories 2, 3 and SM11). Rare longer trajectories could be detected along the cortex of the parasites (such as illustrated for trajectory 1), however they never fitted, with good probability, a model of directed motion. In agreement, analysis of cortical DG trajectories in

PLOS PATHOGENS
Rab11A, a novel regulator of dense granule secretion in Toxoplasma gondii Rab11A-WT and Rab11A-DN parasites revealed a significant increase in the coefficient of diffusion of Rab11A-DN trajectories, suggesting a role for Rab11A in regulating DG directed transport along the parasite cytoskeleton ( Fig 3G). Finally, we performed an experiment in which we washed out 0.5 μM-(S12 Movie) or 1 μM-(S13 Movie) Shield-1 pre-induced Rab11A-DN parasites in order to arrest the expression of the Rab11A-DN protein. We clearly observed, 4h after Shield-1 removal, a strong accumulation of GFP-positive DGs at the parasite PM together with the re-initialization of their content release (S12 Movie), and a pronounced signal at the defective interface between dividing parasites (S13 Movie). This observation suggests that Rab11A may be required for DG docking/tethering at the PM.
Collectively, these data indicate that Rab11A regulates both the directed transport of DG along cytoskeleton tracks ( Fig 1D and Fig 2D, 2E and 2F) and their exocytosis into the PV space.

Rab11A regulates transmembrane protein localization at the PM
Based on our previous study [28], we proposed that Rab11A is required for the delivery of vesicles containing SAG1 and probably other surface proteins, from the endosomal network to the plasmalemma of daughter cells, where new PM is synthesized, similar to its function described in other eukaryotes. This prompted us to investigate whether during replication Rab11A might regulate the localization of other surface proteins in T. gondii. We transiently transfected Rab11A-WT and -DN parasites with plasmids encoding the transmembrane HA-tagged Glucose transporter 1 (GT1) [36], or the Ty-tagged rhomboïd protease 4 (ROM4) [37]. In contrast to the rhomboid protease ROM1 that localizes to micronemes, ROM4 was found to be targeted to the tachyzoite PM, suggesting that it is transported through the constitutive pathway [37] [38]. Similar to DGs, GT1 and ROM4 proteins were retained in intracellular vesicles and were no longer delivered to the parasite PM ( Fig 3H). In addition, we took advantage of the impaired exocytosis activity in Rab11A-DN parasites to study whether different populations of secretory vesicles may co-exist during parasite replication. Co-localization studies in fixed Rab11A-DN parasites showed that ROM4 and GRA3 partially co-localize, but were also detected in distinct vesicular compartments ( Fig 3I). This may reflect a distinct timing of protein synthesis and vesicle release from the Golgi to the PM. However, this observation also suggests the existence of different regulatory pathways for the trafficking of protein localized at the PM vs proteins secreted into the vacuolar space. In particular, transmembrane proteins may be actively recycled during parasite division, as suggested in a previous study on the retromer subunit TgVPS35 [39], and more recently during extracellular parasite motility [40]. Thus, Rab11A may not only play a role in the regulation of DG protein release into the vacuolar space but also in the trafficking of proteins localized at the PM during parasite replication.
Importantly, unlike GRA protein secretion, DG biogenesis was not impaired in Rab11A-DN parasites as assessed by transmission electron microscopy (Fig 4). In addition, supporting a major disturbance in DG exocytosis, the IVN could not be detected in the drastically reduced vacuolar space characterized by the PVM being closely apposed to the parasite DGs and their resulting trajectories, which were analyzed in F-. F-Tracking of DGs in Rab11A-DN expressing parasites indicates mostly confined (as exemplified for DG trajectories 2, 3) and diffusive (trajectories 1, 4) motions. G-Mean diffusion coefficient (D) calculated from 10 cortical trajectories manually tracked in Shield-1 induced Rab11A-WT and Rab11A-DN parasites. Data show mean ± SD (unpaired Student's t-test �� p<0,01). H-IFA showing the glucose transporter GT1 and Romboïd protein ROM4 (green) retained in intra-cytosolic vesicles in Shield-1 induced Rab11A-DN parasites, while being efficiently delivered at the plasma membrane in induced Rab11A-WT parasites. The parasite cortex is delineated by GAP45 (red). Bars: 2 μm. I-IFA showing the localization of the proteins GRA3 (red) and ROM4 (green) in distinct vesicles in Shield-1 (+ S) induced Rab11A-DN expressing parasites during parasite replication. In Rab11A-WT expressing parasites, GRA3 and ROM4 localized at the vacuolar membrane and at the parasite plasma membrane, respectively. Bars: 2 μm. https://doi.org/10.1371/journal.ppat.1008106.g003

PLOS PATHOGENS
Rab11A, a novel regulator of dense granule secretion in Toxoplasma gondii PM (Fig 4B and 4C). We also detected the previously described defect in daughter cell segregation [28] (Fig 4C, arrows). Presumably, in addition to the contribution of the mother cell plasma membrane, delivery of new membrane is required to complete daughter cell segregation at the end of cytokinesis and this process may be regulated by Rab11A. In line, the requirement for de novo lipid synthesis to complete daughter cell segregation has been already demonstrated in other studies [41].

Rab11A regulates adhesion and motility of extracellular parasites
A role for Rab11A in parasite invasion has been previously demonstrated [32]. To explore which steps of parasite entry (e.g. adhesion, motility, and invasion) were altered, we treated extracellular Rab11A-WT and -DN parasites with Shield-1 for 2 h before monitoring their ability to adhere to host cells. We found that Rab11A-DN tachyzoites were severely impaired in surface attachment to human fibroblast (HFF) monolayers compared to Rab11A-WT parasites ( Fig 5A). Furthermore, parasites that successfully adhered exhibited a strong defect in motility, as quantified by the percentage of parasites displaying a SAG1-positive trail deposit ( Fig 5B). Importantly, compared to Rab11A-WT parasites, the morphology of adherent motile Rab11A-DN parasites was altered, the latter being wider and shorter, losing their typical arc shape ( Fig 5C). Analysis of individual parasites imaged by Scanning EM (n = 70 for WT and DN) confirmed that Rab11A-DN parasites display a significant increase in circularity and accordingly, a decrease in the aspect ratio (AR: major axis/minor axis) (Fig 5D). Of note, conoid extrusion was only slightly decreased in Shield-1 induced Rab11A-DN compared to Rab11A-WT ( Fig 5E). Along this line, changes in morphology were not correlated with major perturbations in formation and organization of the sub-pellicular microtubule network in Rab11A-DN parasites compared to Rab11A-WT (Fig 5F).
An impaired recruitment of late glideosome components at daughter cell buds has been previously reported in dividing Rab11A-DN parasites [28] and could account for the motility defect. However, we induced Rab11A-DN protein expression in non-dividing extracellular parasites and accordingly we did not observe any significant defect in the localization of GAP45 and Myosin Light Chain 1 (MLC1) at the cortex of extracellular parasites (S3 Fig). This indicates that the morphological defect observed in Rab11A-DN parasites is not correlated with a significant perturbation of glideosome component localization.
The microneme protein MIC2, a transmembrane protein released at the PM of the parasite, promotes parasite motility and adhesion [42] [43]. First, we confirmed by IFA that MIC2-positive micronemes were detected at the apical pole of extracellular induced Rab11A-DN parasites, indicating no defect in their formation and localization (S4A Fig). Secretion of microneme proteins by extracellular parasites can be triggered by ethanol, a step followed by their release from the parasite PM after cleavage by proteases. Notably, ROM4 has been shown to promote MIC2 trimming at the parasite PM [38]. Since ROM4 was no longer present at the PM of replicating Rab11A-DN parasites, we investigated whether a similar defect could be observed in 2h Shield-1 induced extracellular Rab11A-DN parasites. As previously observed for glideosome components, ROM4 localization at the PM was not perturbed in extracellular induced Rab11A-DN parasites (S4B Fig). Next, we performed excretion/secretion assays to assess the transport of the MIC2 protein to the parasite PM and its subsequent release into the culture medium. Western blot quantification of the Excreted-Secreted Antigen (ESA) fractions demonstrated a significant reduction in MIC2 release upon induction of microneme exocytosis by ethanol ( Fig 5G). Accordingly, a slight increase in MIC2 protein levels was observed in the pellet fraction, also indicating that the decrease in MIC2 secretion is not due to a defect in protein synthesis. As observed by IFA, a reduced level of constitutive GRA1 secretion was also

PLOS PATHOGENS
Rab11A, a novel regulator of dense granule secretion in Toxoplasma gondii

PLOS PATHOGENS
Rab11A, a novel regulator of dense granule secretion in Toxoplasma gondii detected by WB, which correlated with GRA1 accumulation in the parasite pellet fraction ( Fig  5G). Together, these data suggest that the defect of extracellular Rab11A-DN parasites in host cell adhesion and motility may be due to impaired MIC2 delivery to the PM.
Lastly, Rab11A-DN parasites that successfully adhered to the surface of host cells, displayed only a mild defect in host cell invasion (Fig 5H). This was supported by the observation of a correctly formed RON4-positive junction by invading Rab11A-DN parasites (Fig 5H).
Collectively, these data demonstrate that Rab11A promotes parasite invasion by regulating parasite motility and adhesion to host cells, but not the formation of the moving junction.

Rab11A-positive vesicles accumulate at the apical pole during parasite motility and host cell invasion
The active role of Rab11A in parasite motility and adhesion led us to explore the localization of Rab11A in motile extracellular and invading parasites. Live imaging of mcherryRab11A-WT revealed an unexpected polarized accumulation of Rab11A-positive vesicles at two main foci localized at the apical tip of extracellular adhering and motile parasites (Fig 6A, S14 Movie). Apical accumulation of Rab11A appeared to be prolonged during host cell invasion (Fig 6B,  S15 Movie). Apically polarized localization of Rab11A in invading parasites was further confirmed in fixed parasites after labeling of the moving junction with RON4 ( Fig 6C).

Rab11A regulates polarized secretion of DG content during parasite motility and host cell invasion
Next, we assessed whether Rab11A regulates the transport of DGs, not only during parasite replication (Fig 3) but also during parasite motility and invasion. Similarly to the live imaging data (Fig 6), we found Rab11A at two foci localized at the apex of extracellular parasites that had been allowed to move on coverslips prior to fixation (Fig 7A). These Rab11A foci co-localized with the DG protein GRA1 suggesting that Rab11A may regulate apical transport and/or anchoring of DGs at the apical pole of motile extracellular parasites. A similar co-recruitment of Rab11A and DGs at two apical foci was observed during host cell invasion (Fig 7B, white  arrows and S5A Fig). Most importantly, we observed a complete inhibition of this polarized DG apical localization in extracellular motile Rab11A-DN parasites ( Fig 7A) and during host cell invasion (Fig 7B and S5A Fig). This demonstrates that Rab11A regulates the apical accumulation of DGs during the early steps of parasite motility and entry into host cells.

Discussion
In this study, we unraveled an essential role of Rab11A in the delivery of transmembrane proteins to the parasite PM and the release of DG proteins into the vacuolar space during intracellular replication.
In other eukaryotic systems, Rab11A localizes to the endocytic recycling compartment (ERC) and has been implicated in trafficking of internalized receptors from the ERC to the for each condition; unpaired Student's t-test ��� p<0,0001). E-Histogram showing the percentage of Shield-1 induced Rab11A-WT and Rab11A-DN extracellular parasites displaying an extruded conoid (unpaired Student's t-test � p = 0.011). F-IFA showing the localization of the sub-pellicular microtubule network (green) in Shield-1 induced Rab11A-WT and Rab11A-DN (red) parasites. Bars: 2 μm. G-Western blot analysis of excreted-secreted antigen assays (ESA) performed with Shield-1 induced (+S) extracellular RHΔKU80, Rab11A-WT and Rab11A-DN expressing parasites revealed a defect in MIC2 and GRA1 protein secretion. Eno2 was used as a loading control. Secreted MIC2 proteins (ESA fraction) and intracellular GRA1 proteins (pellet fraction) were quantified from 3 independent experiments and expressed as fold-change compared to induced RHΔKU80 parasites. H-Quantification of the percentage of Shield-1 induced extracellular Rab11A-WT and Rab11A-DN expressing parasites, which have invaded host cells. Data show mean ± SEM of three independent experiments (unpaired Student's t-test � p = 0,0488). Fluorescence images show Shield-1 induced mcherryRab11A-DN (red) invading host cells, as illustrated by the moving junction positive for RON4 (green). Bars: 1 μm. https://doi.org/10.1371/journal.ppat.1008106.g005

PLOS PATHOGENS
Rab11A, a novel regulator of dense granule secretion in Toxoplasma gondii PM [8]. Rab11A also localizes to the TGN compartment, where it regulates transport of material from this compartment to the ERC, or to the PM [17]. Similarly, during T. gondii cytokinesis, Rab11A mostly localizes at the Golgi of daughter cells, and at the tip of growing buds, suggesting a polarized transport of de novo synthetized material between these two locations during daughter cell emergence. Interestingly, a similar apically polarized localization of Rab11A was also evident during extracellular parasite motility. Thus, one may envision that components of the apical complex, a microtubule-rich structure from which emanates subpellicular microtubules [44], may control Rab11A-dependent recruitment of specific cargos at the apical pole of the parasite. In particular, RING2, a component of the apical polar ring, was

PLOS PATHOGENS
Rab11A, a novel regulator of dense granule secretion in Toxoplasma gondii shown to function in constitutive and cGMP-stimulated secretion of microneme proteins [45]. Recently, two other components of the apical polar ring, APR1 and the Kinesin A, have also been reported to regulate MIC2 secretion [44]. Hence, it will be of interest to investigate whether Rab11A interacts with components of the apical polar ring to promote recruitment and/or exocytosis of micronemes and DG during extracellular motility and invasion.
Moreover, videomicroscopy recordings of mCherryRab11A-WT in intracellular parasites revealed highly dynamic Rab11A-positive vesicles displaying bidirectional trajectories between the apical and the basal poles, with an accumulation at the basal pole of replicating parasites. This suggests that Rab11A may contribute to the formation of the residual body, an organized structure that interconnects parasites during replication, or in the regulation of parasite exchanges, the two processes being likely tightly correlated. In this context, we observed Rab11A-positive vesicles and tubular-like structures in the region of the residual body. This region was recently reported to harbor a dense actin-myosin network that connects the

PLOS PATHOGENS
Rab11A, a novel regulator of dense granule secretion in Toxoplasma gondii parasites within the PV ensuring synchronous divisions [6] [7]. Thus, Rab11A may also contribute to the regulation of this actin network function and dynamics. Indeed, in plants dysregulated Rab11A activity affects actin organization in the apical region of growing pollen tubes [46]. Supporting the hypothesis of a specific interaction between Rab11A and the actino-myosin cytoskeleton, depolymerizing actin filaments alters Rab11A-positive vesicle displacements. A role for the complex Myosin Vb-FIP2-Rab11A in promoting actin-mediated transport of vesicles has been previously observed in mammalian cells [47] [48] [49]. So far, no homologues of Rab11-family interacting proteins (FIPs) have been identified in T. gondii and Plasmodium. Nonetheless, P. falciparum Rab11A was reported to directly interact with the myosin light chain 1 (MLC1/MTIP), which therefore links Rab11A-mediated vesicular transport to unconventional myosins and the actin cytoskeleton [28]. As actin depolymerization resulted in the formation of both cytosolic and peripheral Rab11A-positive static clusters, it's possible that distinct myosins regulate different steps of Rab11A/DG transport e.g. MyoF in the cytosol and from the TGN [33], MyoA at the parasite cortex where the glideosome is located [28], and MyoJ in the cell-to-cell connecting network [6]. Further studies using parasite strains deleted for these molecular motors will address this question.
Co-distribution studies indicated that Rab11A-positive vesicles associate with dense granules in a dynamic manner. However, we did not observe Rab11A at the limiting membrane of DG. Rather, these two compartments appear to transiently dock one with each other enabling joint transient motions that were particularly evident at the cortex of the parasite. Indeed, tracking of the trajectories of both Rab11A-positive vesicles and DG revealed that Rab11Apositive vesicles promoted DG anchoring at the parasite cortex and their rapid "directed" transport. This mode of transport called "hitchhiking" has been recently described in different cell types and has emerged as a novel mechanism to control organelle movement [50]. During this process, the "hitchhiker" benefits from distinct molecular motors present at the surface of the "vehicle". In addition, endosomes represent multifunctional platforms that receive specific signals and could drive transport of hitchhiker cargo to particular regions of the cell. Notably, co-movement of cargo may facilitate interactions at membrane contact sites important for organelle maturation, fusion and/or material exchange. Related to this last aspect, we found that over-expression of Rab11A-DN led to a complete block in DG secretion. We observed that restoration of Rab11A functions by washing out Shield-1 correlated with an accumulation of Rab11A-positive vesicles at the parasite plasma membrane suggesting a role for Rab11A in vesicle docking/tethering at the PM, which remains to be formally demonstrated. In other eukaryotic systems, Rab11A is known to promote vesicle docking and fusion at the PM via its interaction with the exocyst complex and SNARE proteins, respectively [15]. However, homologues of the different exocyst complex subunits could not be identified in T. gondii [51]. Thus, unexplored mechanisms of Rab11A-mediated vesicle docking at the PM may exist in T. gondii and Rab11A-interacting SNAREs remain to be identified. One may envision that Rab11A drives DGs to sites that favor exocytosis by promoting interactions with regulatory factors involved in vesicle fusion.
Benefiting from the fast and efficient induction of the Rab11A-DN protein expression in extracellular parasites, we confirmed the previously described defect in host cell invasion [30]. Of note, our numerous attempts to generate parasites expressing C-terminal tagged Rab11A failed, and notably, our attempts to apply the rapidly inducible AID knock-down system also failed [52]. This is likely due to the fact that the C-terminal domain of the Rabs contains one or two cysteines recognized by geranylgeranyl-transferases to induce their isoprenylation, a modification required for their association with membranes. The impaired cell invasion of Rab11A-DN expressing parasites results from a strong defect in parasite adhesion to host cells. Indeed, parasites that successfully adhered to host cells were only mildly perturbed in host cell

PLOS PATHOGENS
Rab11A, a novel regulator of dense granule secretion in Toxoplasma gondii entry. Moreover, secretion of MIC2, an adhesin essential for parasite adhesion and motility was reduced upon dysregulation of Rab11A activity. Secretion of the GPI-anchored protein SAG1 is also altered in Rab11A-DN expressing parasites [28]. Thus, it's likely that the altered secretion of these two host cell adhesins contributes to the decrease in adhesion and motility of Rab11A-DN parasites. Consistent with a role of Rab11A in the regulation of surface protein trafficking, we also found a strong defect in the localization of the romboïd protease ROM4 and the glucose transporter GT1 at the PM, indicating a broader role of Rab11A in the regulation of surface protein trafficking. Presumably, distinct exocytic pathways exist in T. gondii, such as described in other organisms. In particular, whether a distinct endosome recycling compartment is present in T. gondii requires further exploration. Previous studies highlighted that T. gondi has functionally repurposed its endocytic system to serve as secretory pathway of this fast replicating intracellular parasite [51] [53]. In this context, the TGN appears to be a hybrid compartment to which the endosomal markers (Rab5 and Rab7) are tightly associated [31]. Therefore, one might envision that material internalized from the PM reaches this hydrid TGN/ELC compartment before being re-directed to other target membranes, such as the rhoptries, the PM, and the degradative vacuole (VAC). Such a recycling process has been recently observed during extracellular parasite motility [38]. Recycling of mother material during daughter cell emergence may also follow this indirect secretory pathway, while de novo synthetized proteins may traffic directly from the TGN to the PM.
Finally, during extracellular parasite motility and invasion, imaging of both live and fixed parasites revealed an unexpected polarized accumulation of Rab11A-positive vesicles towards two main foci located just beneath the conoid. In mammalian cells, Rab11A-dependent polarized secretion towards the leading edge of motile cells is essential to promote persistent migration [54]. This process not only provides additional membrane ensuring the extension of the leading edge, but also contributes to the translocation of regulatory factors involved in actin and microtubule cytoskeleton activity. In T. gondii, apical delivery of some effectors may regulate actin and microtubule cytoskeleton activity and thereby parasite motility. Such regulatory mechanisms have been demonstrated for the lysine methyltransferase, AKMT (Apical complex lysine (K) methyltransferase) localized at the conoid [55]. It has been also recently shown that the DG protein GRA8 contributes to parasite motility by regulating conoid extrusion and organization of the microtubule network [56]. Thus, future research will aim to identify the cargos that are apically delivered in a Rab11A-dependent manner and their putative role in regulating parasite motility. Interestingly, the apical accumulation of DGs that we observed in extracellular motile parasites has been previously described during parasite invasion [57]. Thus, an alternative explanation would be that the parasite "prepares its arrival" at the host cell, anticipating the burst of DG secretion that occurs during invasion by promoting their anchoring at the apical pole. In such a scenario, a second signal (vacuole closure?) would then trigger their fusion and content release into the vacuolar space.
Importantly, at present we cannot explain the defect in morphology we observed in extracellular motile Rab11A-DN parasites. This may be linked to a dysregulation of the actin cytoskeleton activity, or be related to a defect in the dynamics of an endo-exocytosis activity required for parasite forward movement, both potentially leading to shape deformation of the moving parasite. Along with the first hypothesis, actin staining (Cb-E transfected parasites) in extracellular parasites that were allowed to move on coated coverslips before fixation revealed a strong accumulation of actin at the basal pole of the parasite (S5B Fig), such as recently described in invading parasites [58].
Therefore, identifying Rab11A interactors will be an important future goal, as it will improve our understanding of the mechanisms regulating the distinct exocytic pathways in T. gondii. In particular, it will be important to characterize the molecular mechanisms involved

PLOS PATHOGENS
Rab11A, a novel regulator of dense granule secretion in Toxoplasma gondii in anchoring Rab11A-positive vesicle to actin or microtubule molecular motors, and of a possible process of vesicle docking/tethering at the PM, both during parasite motility and intracellular replication. Finally, exploring a putative functional interaction between Rab11A and the apical complex may lead to the discovery of novel regulated secretory mechanisms essential to ensure parasite virulence.

Parasite culture and transfection
Toxoplasma gondii Type I RHΔKU80ΔHXGPRT parasites were grown on confluent Human Foreskin Fibroblast (HFF) cells (CCD-1112Sk (ATCC)) which were cultured in complete DMEM (gibcoLife Technologies) supplemented with 10% Fetal Bovine Serum (GibcoLife Technologies) and 1% Penicillin/Streptavidin (gibcoLife Technologies). To obtain the transgenic parasites, the RHΔKU80ΔHXGPRT parental strain was transfected by electroporation following standard procedures with 50μg of the plasmids listed in Table 1.
Following transfection, parasites were subjected to drug selection and verified for the transfection efficiency by immunofluorescence analysis. Subsequently the parasites were subjected to cloning by serial dilution.

Production of the anti-Rab11A antibodies
Recombinant purified GST-Rab11A protein was used to raise a TgRab11A specific mouse polyclonal antibody. The cleavage site present between the GST tag and Rab11A was digested with Precision protease (GE life science). GST-Rab11A bound to agarose beads was washed with 10 bed volumes of Cleavage buffer (50mM Tris HCl, pH7.0, 150mM NaCl, 1mM EDTA, 1mM DTT) at 4˚C. Precision protease (40 units) was added to the cleavage buffer and incubated with the beads at 4˚C overnight. The purified Rab11A was collected in the supernatant. 50μg of the purified recombinant protein suspended in Freund's Adjuvant were injected intraperitoneally into mice over a series of 4 boosts. Following the third boost, a sample of serum was collected and tested by western blot for antibody reactivity using a total protein extract of parasites. Once specific antibody activity was detected mice were sacrificed and serum collected and stored at -20˚C.

Protein sample preparation and Western Blot
Parasites were lysed in a TritonX100 1%, NaCl 150mM, Hepes 10mM buffer and the lysate was loaded onto a 10% SDS polyacrylamide gel using. Proteins were blotted onto a nitrocellulose membrane and probed with the indicated primary antibodies followed by species-specific secondary antibodies conjugated with HRP. Antibodies were diluted in 5% milk dissolved in TNT buffer (0.1M Tris HCl pH7.6, 0.15M NaCl and 0.1% v/v Tween20). The probed nitrocellulose membranes were visualized using the ECL Western blotting substrate (Pierce).

PLOS PATHOGENS
Rab11A, a novel regulator of dense granule secretion in Toxoplasma gondii

Immunofluorescence assay
Confluent HFF monolayers were grown on coverslips and infected with parasites prior to fixing with 4% PFA for 15 min. After quenching with 50mM NH 4 Cl, the coverslips were permeabilized with 0.2% triton dissolved in 5% FBS-PBS for 30 min. Coverslips were incubated with primary antibodies in 0.1% triton dissolved in 2% FBS-PBS and then washed thrice with 1X PBS. Alternatively, the coverslips were incubated with primary antibodies in 0.01% Saponin diluted in 2%FBS-PBS for 1 h. Incubation with secondary antibodies was performed in 0.1% triton or 0.01% Saponin dissolved in 2%FBS-PBS for 30 min. To label invading parasites, freshly egressed extracellular parasites expressing Rab11A-WT and Rab11A-DN were induced with Shield-1 for 2 h and seeded onto HFF monolayers in a 24-well plate at a concentration of 2 � 10 6 parasites (Rab11A-WT) and 4 � 10 6 (Rab11A-DN) /500μl complete medium containing Shield-1 per coverslip. The plate was centrifuged for 2 min at 1000rpm at room temperature to trigger adhesion and synchronized invasion events. The plate was immediately shifted to a water bath at 37˚C and the parasites were fixed with 4% PFA-sucrose at the following time points-0, 2 and 5 min. Coverslips were washed with PBS and adherent or invading parasites labeled without permeabilization with the anti-SAG1 antibody and a secondary anti-mouse AlexaFluor405 antibody. After washing with PBS, parasites were permeabilized with 0.05% saponin for 10 min, followed by a blocking step with 5% FBS-PBS for 30 min. Next, coverslips were incubated with rabbit anti-RON4 antibodies and secondary anti-rabbit AlexaFluor488 to label the moving junction. Depending on the experiment, additional primary antibodies were added to detect GRA1, GRA3 and Rab11A during parasite invasion. Images were acquired using a Zeiss LSM880 confocal microscope equipped with an airyscan module. The antibodies used in this study are listed in Table 2.

PLOS PATHOGENS
Rab11A, a novel regulator of dense granule secretion in Toxoplasma gondii

Invasion assay
Freshly egressed extracellular parasites expressing Rab11A-WT and Rab11A-DN were harvested and treated for 2 h with 1 μM of Shield-1. Induced parasites were counted and seeded onto HFF monolayers in a 24-well plate at a concentration of 2 � 10 6 parasites (Rab11A-WT) or 4 � 10 6 parasites (Rab11A-DN) / 500μl complete medium containing Shield-1 / coverslip. The plate was centrifuged for 2 min at 1000rpm at RT to trigger immediate adhesion and synchronized invasion events. Parasites were then shifted to 37˚C for 1 h. The coverslips were washed with PBS-three times prior to fixation. Cells were fixed in 4% PFA for 10 min and subjected to a red/green invasion assay. Briefly, adherent external parasites were labeled without permeabilization with mouse anti-TgSAG1 antibodies, followed by secondary anti-mouse antibodies coupled to Alexa488. After cell permeabilization with TritonX100 0.1%, invaded intracellular parasites were detected using rabbit anti-TgGAP45 antibodies followed with a secondary antirabbit antibodies coupled to Alexa594. All parasites labeled green-red were considered as extracellular, while parasites exclusively red (positive for GAP45) were considered intracellular. At least 300 parasites were counted for each condition performed in triplicate. Data represent mean values ± SEM from three independent biological experiments.

Motility (Trail deposition) Assay
Glass slides were coated with 100μg/ml BSA-PBS and incubated at 37˚C for 1 h. The slides were washed three times with PBS and allowed to dry. Freshly egressed extracellular Rab11A-WT and Rab11A-DN expressing parasites were harvested and treated for 2 h with 1 μM of Shield-1. Induced parasites were counted and suspended in HHE buffer (HBSS, 10mM HEPES, 1mM EGTA) containing 1 μM of Shield-1. 1 � 10 6 (Rab11A-WT) or 2 � 10 6 (Rab11A-DN) parasites were seeded per well and incubated for 15 min at 37˚C. Parasites were then fixed with 4% PFA in PBS for 10 min at RT. A standard IFA protocol was followed wherein primary mouse anti-SAG1 antibodies were used followed by goat anti-mouse secondary antibodies conjugated to Alexa Fluor 488. 200 parasites per coverslip were counted for the presence or absence of a SAG1-positive trail. With internal triplicates, the experiment was performed 3-times. Mean values ± SEM was calculated.

Adhesion assay
Freshly egressed extracellular Rab11A-WT and Rab11A-DN parasites were harvested and treated for 2 h with 1 μM of Shield-1. Parasites were then counted and resuspended in Endo buffer (44.7mM K 2 SO 4 , 10mM Mg 2 SO 4 , 100mM sucrose, 5mM glucose, 20mM Tris, 0.35% wt/vol BSA-pH 8.2) containing 1 μM cytochalasin D and 1 μM of Shield-1. 2 � 10 6 parasites were then seeded onto confluent HFF cells grown on glass coverslips, spun down for 2 min at 1000rpm and incubated for 15 min at 37˚C in the presence of 1 μM cytochalasin D and Shield-1. The coverslips were washed with PBS before fixation with PFA 4% for 10 min. The Red/ Green assay was performed (see "Invasion assay"). At least 200 parasites were counted for each condition performed in triplicate. Data represent mean values ± SEM from three independent biological experiments.

Transmission electron microscopy (TEM)
After infection of a confluent HFF monolayer, cells containing replicating shield-1 induced Rab11A-WT and Rab11A-DN expressing parasites were detached with a scraper, spun down and fixed with 1% glutaraldehyde in 0.1 M sodium cacodylate pH 6.8 overnight at 4˚C. Cells were post-fixed with 1% osmium tetroxide and 1.5% potassium ferricyanide for 1 h, then with 1% uranyl acetate for 45 min, both in distilled water at RT in the dark. After washing, cells were dehydrated in graded ethanol solutions then finally infiltrated with epoxy resin and cured for 48 hs at 60˚C. Sections of 70-80 nm thickness on formvar-coated grids were observed with a Hitachi H7500 TEM (Elexience, France), and images were acquired with a 1 Mpixel digital camera from AMT (Elexience, France).

Scanning Electron microscopy (SEM)
Parasites were seeded on BSA coated-glass coverslips for 15 min at 37˚C before being fixed with 2.5% glutaraldehyde in 0.1 M sodium cacodylate for 30 min. After washing, cells were treated with 1% osmium tetroxide in water, in the dark for 1 hour. Cells were dehydrated with increasing ethanol concentration baths. After two pure ethanol baths, cells were air-dried with hexamethyldisilazan (HMDS). Finally, dry coverslips were mounted on stubs and coated with 5 nm platinum (Quorum Technologies Q150T, Milexia, France) and cells were imaged at 2 kV by a secondary electron detector with a Zeiss Merlin Compact VP SEM (Zeiss, France). The circularity and aspect ratio (AR) parameters were calculated using the ImageJ software. First, each parasite contour was manually delineated on the SEM images and the pluggin "Analyze"> "Set measurements"> Shape descriptor" was applied for all defined ROI (individual parasites, n = 70) to extract the circularity and aspect ratio parameters. Circularity describes how close an object is to a true circle and is calculated using the formula: circularity = 4D � (area/perimeter 2 ). A circularity value of 1 indicates a perfect circle. As the value approaches 0, it indicates an increasingly elongated shape. The aspect ratio describes the proportional relationship of an object's width to its height and is calculated using the formula: AR = major axis/minor axis.

PLOS PATHOGENS
Rab11A, a novel regulator of dense granule secretion in Toxoplasma gondii Videomicroscopy Time-lapse video microscopy was conducted in LabTek chambers installed on an Eclipse Ti inverted confocal microscope (Nikon France Instruments, Champigny sur Marne, France) with a temperature and CO 2 -controlled stage and chamber (Okolab), equipped with two Prime 95B Scientific Cameras (Photometrics, UK) and a CSU W1 spinning disk (Yokogawa, Roper Scientific, France). The microscope was piloted using MetaMorph software (Universal Imaging Corporation, Roper Scientific, France). A live-SR module (Gataca Systems, France) was added to the system to improve the obtained resolutions. Exposure time of 500 ms was used for the simultaneous acquisition of the GFP and mCherry channels, in dual camera mode (with band pass filters 525/50 nm and 578/105 nm, dichroic mirror at 560 nm, and laser excitation at 488 nm and 561 nm). Videos were captured at 2 frame/second.

Automatic tracking and vesicle co-distribution using the Imaris software
Automatic tracking of vesicles using the Imaris software (Bitplane, Oxford Instruments) was applied on the recorded videos retrieved from the GFP and mcherry channels of SAGΔG-PI-GFP / mcherryRab11A-WT expressing parasites. We first used the tool "Spot detector" for selecting-filtering spot size and intensity values for each channel. Next, we manually removed detection of false GFP-positive spots (notably detected in the vacuolar space due to the secretion of the SAGΔGPI protein). The tool "Track Manager" was used to manually correct the obtained tracks when required and to extract the xy positions of a given spot over time enabling to calculate the Mean Square Displacement (MSD) using MATLAB (see below). The tool "spot co-localization" was used to calculate the percentage of co-distribution between DG and Rab11A-postive vesicles. A distance of 300 nm between the spots was selected corresponding to the average size of the vesicles. At a given time point and for the entire vacuole, the number of all detected green spots, as well as the number of green spots co-distributing with the red spots were extracted to calculate the co-distribution percentage. This was repeated over 5 consecutive time points every 2 s for the first 10 s of recording to avoid bleaching of the fluorescent signals. The mean co-distribution percentage over these 5 time points was calculated per vacuole. The mean +/-SD of 10 vacuoles was then calculated.

Manual tracking and mathematical modeling with MATLAB
When indicated, the manual tracking plugin from the ImageJ software (https://imagej-nihgov/ij/) was applied on the images obtained with the MetaMorph software to extract in time the spatial xy positions of the fluorescent vesicles. In order to track and model the type of motion of the vesicle, images were processed in MATLAB (www.mathworks.com) by applying fit function ('poly1' or 'poly2' options).
MSD was calculated thanks to a MATLAB script according to the formula:

PLOS PATHOGENS
Rab11A, a novel regulator of dense granule secretion in Toxoplasma gondii

Statistics
Means and SEM / SD were calculated in GraphPad (Prism). P-values were calculated using the Student's t-test assuming equal variance, unpaired samples and using two-tailed distribution.

PLOS PATHOGENS
Rab11A, a novel regulator of dense granule secretion in Toxoplasma gondii