Skip to main content
  • Research article
  • Open access
  • Published:

Combined nanometric and phylogenetic analysis of unique endocytic compartments in Giardia lamblia sheds light on the evolution of endocytosis in Metamonada



Giardia lamblia, a parasitic protist of the Metamonada supergroup, has evolved one of the most diverged endocytic compartment systems investigated so far. Peripheral endocytic compartments, currently known as peripheral vesicles or vacuoles (PVs), perform bulk uptake of fluid phase material which is then digested and sorted either to the cell cytosol or back to the extracellular space.


Here, we present a quantitative morphological characterization of these organelles using volumetric electron microscopy and super-resolution microscopy (SRM). We defined a morphological classification for the heterogenous population of PVs and performed a comparative analysis of PVs and endosome-like organelles in representatives of phylogenetically related taxa, Spironucleus spp. and Tritrichomonas foetus. To investigate the as-yet insufficiently understood connection between PVs and clathrin assemblies in G. lamblia, we further performed an in-depth search for two key elements of the endocytic machinery, clathrin heavy chain (CHC) and clathrin light chain (CLC), across different lineages in Metamonada. Our data point to the loss of a bona fide CLC in the last Fornicata common ancestor (LFCA) with the emergence of a protein analogous to CLC (GlACLC) in the Giardia genus. Finally, the location of clathrin in the various compartments was quantified.


Taken together, this provides the first comprehensive nanometric view of Giardia’s endocytic system architecture and sheds light on the evolution of GlACLC analogues in the Fornicata supergroup and, specific to Giardia, as a possible adaptation to the formation and maintenance of stable clathrin assemblies at PVs.


Endomembrane compartments, while present in a few prokaryotic lineages [1], have evolved and greatly diversified across eukaryotic lineages. A fundamental task performed by some membrane-bounded organelles is endocytosis—the controlled and directed uptake of nutrients and other materials from the extracellular space into the cell by membrane transport. Fluid phase or receptor-bound material at the cell surface is internalized via invaginations and formation of vesicles at the plasma membrane, mediated by clathrin-coated vesicles (CCVs) [2, 3]. In turn, CCVs fuse with early endosomes which mature into late endosomes upon lysosome fusion [4, 5]. Clathrin coats are also involved in protein secretion forming exocytic transport vesicles derived from the trans-Golgi compartment and play a role in Golgi apparatus reassembly after mitotic cell division [6, 7].

Evolutionary adaptations of endocytic pathways to specific environmental niches and nutrient sources are especially relevant to species adopting a fully parasitic or commensal lifestyle [8,9,10,11]. Within the extant Metamonada supergroup [12,13,14], the parasitic protist Giardia lamblia (syn.: intestinalis or duodenalis) evolved a distinct endocytic pathway, which reflects its adaptation to the host intestinal lumen environment. This unicellular parasite is responsible for > 300 million cases annually of water-borne infections causing gastroenteritis—giardiasis—with a higher incidence in low- to middle-income countries [15]. Giardia is the etiological agent for symptomatic gastroenteritis in 15% of children in developing countries, with 1–2% fatality in children with severely compromised health status [16, 17]. There is a strong association of Giardia infections with chronic conditions such as irritable bowel syndrome or inflammatory bowel disease as a result of intestinal barrier function disruption and microbiome dysregulation [18, 19].

The cellular evolution of the Giardia genus as an obligate parasite adapted to the small intestinal niche of vertebrates is characterized by a reduction in subcellular compartment diversity. Peroxisomes, late endosomes and a permanent stacked Golgi complex have not been detected in Giardia [20]. Two nuclei [21], an extensive endoplasmic reticulum (ER) [22], highly reduced mitochondria-derived organelles—the mitosomes [23]—and peripheral vesicles (PVs) [24] are the only membrane-bounded organelles with conserved morphology and function documented in the Giardia trophozoite [25,26,27].

The complex array of PV organelles as the only documented endocytic membrane compartment system in Giardia is responsible for the uptake of fluid-phase and membrane-bound material [27,28,29]. These organelles acidify and presumably serve as digestive compartments with the capability for sorting after processing, similar to early and late endosomes and lysosomes [24]. The static system of PV organelles [27, 30] is mainly restricted to the peripheral cortex below the plasma membrane (PM) of the Giardia trophozoite dorsal side. PV morphology was investigated using high-resolution electron microscopy serial sectioning and three-dimensional reconstruction [27]. These organelles were resolved as tubular structures in close proximity to funnel-shaped invaginations of the PM [27]. Our current working model for bulk fluid-phase uptake of extracellular material into PVs invokes a “kiss and flush” mechanism, whereby acidified PV membranes and the PM transiently form channels at invaginations allowing exchange between PV lumen content and the extracellular space at regular intervals. Endocytosed material is digested and transported towards the cell interior while residual material and waste are flushed to the extracellular space in the next round of membrane fusion, thus completing the PV cycle [26, 27]. The characterized molecular complement of the PV surface includes ESCRT components [11] and focal long-lived albeit non-vesicular accumulations of clathrin heavy chain (CHC) molecules and their main interactors, collectively termed clathrin assemblies [27]. A notable member of these is a putative Giardia clathrin light chain GlCLC (ORF 4259), a scarcely annotated protein presenting strong interaction with GlCHC, similar dynamics as measured by FRAP, matching sub-cellular localization and considerable (ca. 50%) overlap of its predicted 3D structure compared to the corresponding prediction for well-characterized CLCs [27]. The function of these stable focal assemblies, as well as additional components at the interface of the PV membranes and the PM, has proved elusive [27]. However, a transient association of several members of the family of adaptor proteins (AP) suggests a role in dynamic processes linked to the uptake of fluid-phase and receptor-bound material into PVs [26, 27].

In this report, we address open questions concerning G. lamblia’s PV ultrastructure and its associated molecular machinery in a comparative approach with one closely and one more distantly related fornicata and metamonada species, Spironucleus sp. and Tritrichomonas foetus, respectively. Using volumetric electron microscopy and super-resolution light microscopy, we developed a classification of PVs based on organelle morphology. Comparative analysis of Giardia’s PVs with endocytic compartments of fornicata and metamonada species emphasized the genus-specific nature of the Giardia endocytic system architecture. In addition, using a combination of co-immunoprecipitation and phylogeny techniques, we provide evidence that the putative clathrin light chain GlCLC [27] is unique to the Giardia genus and evolved de novo as structurally analogous to CLC after the loss of a bona fide CLC in the last Fornicata common ancestor (LFCA). Taken together, the emergence of a unique and highly polymorphic endocytic system such as the one found in the genus Giardia is linked to the proposed convergent evolution of an independent CLC analogue concomitant with the loss of a mostly conserved CLC orthologue.


Complete FIB-SEM rendering of a G. lamblia trophozoite reveals a novel landscape of vesicular compartments

Volumetric scanning electron microscopy (vSEM) is currently considered the gold standard for the determination of biological ultrastructure [31]. One type of vSEM is focused ion beam electron scanning microscopy (FIB-SEM) where a beam of gallium ions is used to mill and image consecutive layers of an embedded biological sample, resulting in a voxel resolution as low as 1–2 nm [32]. This technique allows for sectioning and imaging of entire cells [33] and was previously implemented for the partial rendering of G. lamblia trophozoite sections [27, 34]. Here, we sectioned for the first time a complete G. lamblia trophozoite at a voxel resolution of 125 nm3 (5 × 5 × 5 nm) after high-pressure freezing (HPF) and embedding. Images representing the sagittal plane adjacent to the cell centre (Fig. 1A, D and Additional file 1: Fig. S1A and S1B) show all the major cell compartments such as the nuclei (Fig. 1A, N), the endoplasmic reticulum (Fig. 1A, ER), mitosomes (Fig. 1A, C, m) and elements of the cytoskeleton: axonemes (Ax), funis (F) and the ventral disc (VD) (Fig. 1D [35];). Two different types of small cytoplasmic organelles are observed: PVs (arrow heads) with heterogenous morphology and smaller and electron-dense membrane vesicles of uniform size and appearance which we termed small vesicles (SVs; asterisks) (Fig. 1B, E).

Fig. 1
figure 1

Complete scanning of a Giardia trophozoite by focused ion beam scanning electron microscopy (FIB-SEM). A, D A whole G. lamblia trophozoite was scanned with an isotropic resolution of 5 nm showing the nucleus (N), the endoplasmic reticulum (ER), cytoskeletal features such as the ventral disc (VD), funis (F) and axonemes (Ax) and mitosomes (m). Peripheral vesicles (PV) are marked by arrowheads. Smaller, electron-dense vesicles (small vesicles (SV)) are also documented (asterisks). B Region of interest of A highlighting PVs of different morphology (arrowheads). C Mitosomes proximal to ER membrane. E Region of interest of D highlighting SVs (asterisks). F Full reconstruction of PVs with ilastik and rendering in Imaris showing at least two PV morphologies: spherical (green) and tubular (violet). Violin plot: 403 spherical and 64 tubular PVs were segmented out. Spherical PVs average a volume of 0.0009243 ± 0.0003322 μm3 with a 95% confidence interval between 9.022 × 10−4 and 9.464 × 10−4 μm3. Tubular PVs average a volume of 0.001067 ± 0.0003322 μm3 with a 95% confidence interval between 0.0009843 and 0.001150 μm3. G Tubular PV lengths are larger than their diameter (p-value < 0.0001) and also larger than spherical PVs (p-value < 0.0001) (values presented in log10 form). H Diameters of spherical PVs average 120.1 ± 9.507 μm (95% CI: [119.2; 121.0]). Tubular PV diameters average 43.01 ± 7.816 μm (95% CI: [41.06; 44.96]). Tubular PV lengths average 771.3 ± 266.2 μm (95% CI: [704.8; 837.8]). I Reconstruction of 269 SVs averages volume at 0.0002525 ± 9.280 × 10-5 μm3 with a 95% confidence interval between 0.0002414 and 0.0002637 μm3 (left violin plot). This equals an average diameter of 77.23 ± 9.666 nm with a 95% confidence interval between 76.07 and 78.40 nm. Spherical PV diameter averages 120.1 ± 9.507 nm with a 95% confidence interval between 119.2 and 121 nm (right violin plot). The difference in the diameter between SVs and spherical PVs is statistically significant (****p-value < 0.0001, t-Student significance test). Scale bars: A, D 2 μm; B, C and E 500 nm. ROI, region of interest

After serial sectioning and alignment with TrakEM [36], we used the supervised machine learning (ML) tool ilastik for pixel-based image segmentation of PVs and SVs [37, 38]. The algorithm collection performs supervised learning and recognition of patterns based on ground truth training provided by the user. Patterns are sorted into classes. Once the algorithm is trained on a subset of image data, it is used to analyse complete datasets and assign features to different classes following a decision tree method [39, 40]. This process enabled the three-dimensional rendering of selected trophozoite features: the complete cytoskeleton, the ER, PVs and mitosomes (Additional file 1: Fig. S1C). In addition, we were able to calculate the volume of the cell at 138 μm3 as well as the average volume of mitosome organelles (N = 14) at 0.001093 ± 0.0005698μm3 with a 95% confidence interval between 0.0007643 and 0.001422 μm3 (Additional file 1: Fig. S1D).

Similarly, supervised ML-assisted pixel segmentation and object clustering analysis allowed the identification of two statistically distinct morphological classes of PVs: spherical and tubular/elongated PVs. Individual PV organelles of both classes (N = 467) were rendered in three dimensions (Fig. 1F and Additional file 2: Video S1). Spherical PVs average a volume of 9.243 × 10−4 ± 3.322 × 10−4 μm3 with a 95% confidence interval between 9.022 × 10−4 and 9.464 × 10−4 μm3 while tubular PVs average a volume of 1.067 × 10−3 ± 3322 × 10−4 μm3 with a 95% confidence interval between 9.843 × 10−4 and 1.150 × 10−3 μm3, a statistically significant difference (t-Student test, p < 0.0001), corroborating PV grouping in these two classes. We also determined the length of tubular PVs to be significantly larger than their diameter and also larger than the diameter of spherical PVs (Fig. 1G). Spherical PVs average a diameter of 120.1 ± 9.507 μm (95% CI: [119.2; 121.0]). The diameter of tubular PVs is calculated to average 43.01 ± 7.816 μm (95% CI: [41.06; 44.96]). Tubular PVs average a length of 771.3 ± 266.2 μm (95% CI: [704.8; 837.8]) (Fig. 1H).

To further investigate the morphological heterogeneity of PVs, we analysed trophozoite ultrastructure using freeze-fracture scanning electron microscopy. We documented PV heterogeneity and the presence of spherical and tubular PV forms (Additional file 3: Fig. S2). Additional ultrastructural studies using transmission electron microscopy (TEM) were consistent with this classification (Additional file 4: Fig. S3A and B).

We proceeded with the rendering of 269 SVs—small spherical vesicles, with distinctly higher electron density than PVs and what could be a coat on the cytoplasmic side of the delimiting membrane (Fig. 1I). SVs were also identified by TEM (Additional file 4: Fig. S3), proximal to the PM. SVs average a volume of 2.525 × 10−4 ± 9.280 × 10−5 μm3 with a 95% confidence interval between 2.414 × 10−4 and 2.637 × 10−4 μm3 (Fig. 1I, left violin plot). This equals an average diameter of 77.23 ± 9.666 nm with a 95% confidence interval between 76.07 and 78.40 nm, differing significantly from spherical PVs which average 120.1 ± 9.507 nm with a 95% confidence interval between 119.2 and 121 nm (p < 0.0001) (Fig. 1I, right violin plot). Thus, there is statistical support for SVs as a distinct category of membrane-bounded vesicles (Additional file 4: Fig. S3A and S3C).

Taken together, these findings lead us to hypothesize that, unlike previously thought, PVs are morphologically heterogenous and may comprise different functional categories [41,42,43]. However, these data are currently insufficient to determine whether distinct morphologies correlate with distinct functions.

Combining super-resolution microscopy with ML-assisted image analysis identifies three classes of endocytic compartments in G. lamblia trophozoites

FIB-SEM as a technique is not well-suited to the investigation of large cell numbers, and TEM cannot readily provide 3D volumetric information on subcellular compartments. Hence, to address PV heterogeneity in more detail, we continued our investigation of Giardia endocytic compartments by super-resolution light microscopy (SRM) techniques and ML-assisted image analysis of compartment shapes.

The dimensions of Giardia endocytic compartments are well below the diffraction limit of conventional light microscopy [44]. To overcome the Abbe diffraction barrier, we used stimulated emission depletion microscopy (STED), potentially achieving a lateral (x, y) and axial (z) resolution of 25–50 nm and 60–100 nm, respectively. This technique overall decreases the point spread function signal from the illuminated region [45, 46] and allows for accurate imaging of trophozoite PV lumina loaded with a highly photostable fluid phase marker (10-kDa Dextran-Alexa Fluor 594) which is readily taken up into PVs via the fluid phase endocytic pathway (Fig. 2A and Additional file 5: Video S2). In addition to spherical and tubular PVs documented in FIB-SEM, using STED, we also determined the presence of polymorphic dextran-labelled organelles, i.e. spherical PVs with elongated rods (Fig. 2A, representative image of N = 15 cells). All labelled PVs were further analysed using the ML-assisted algorithm of the ilastik program suite. We first performed a supervised pixel segmentation followed by a supervised object classification. In this second step, we defined and trained the classifier in three organelle morphologies: spherical, tubular and polymorphic. The latter comprised characteristics of both vesicular and tubular classes, generally with spherical centres with attached tubular protrusions (Fig. 2B). After organelle classification, we measured their projected areas. Spherical organelles (N = 1684) have an average projected area of 0.0205 ± 0.0169 μm2 with a 95% confidence interval between 0.0197 and 0.0213 μm2. Tubular endocytic organelles (N = 835) present an average projected area of 0.0453 ± 0.0278 μm2 with a 95% confidence interval between 0.0435 and 0.0472 μm2. Polymorphic organelles (N = 400) have an average projected area of 0.0981 ± 0.0429 μm2 with a 95% confidence interval between 0.0939 and 0.102 μm2. ANOVA analysis reveals that each of the three categories is indeed significantly distinct (p < 0.0001) based on the projected surface area (Fig. 2C). This lends further support to the possibility that PV morphological heterogeneity may have functional implications.

Fig. 2
figure 2

Super-resolution imaging of Giardia lamblia peripheral vesicles with stimulated emission depletion (STED). A Giardia trophozoites loaded with 10-kDa Dextran-AlexaFluor 594 were imaged using confocal and STED microscopy. Dorsal (upper row) and ventral (lower row) regions are represented. In contrast to confocal imaging, STED microscopy allows to separate individual organelles and to visualize different endocytic compartment morphologies. ROI, region of interest. B Organelle segmentation with ilastik distinguishes three dextran-labelled PV categories. C PV areas were calculated post-segmentation on maximum projections of the dorsal regions of 15 cells, using ilastik. Spherical PVs (green, N = 1684) have an average area of 0.0205 ± 0.0169 μm2 with a 95% confidence interval between 0.0197 and 0.0213, tubular PVs (blue, N = 835) have an average area of 0.0453 ± 0.0278 μm2 with a 95% confidence interval between 0.0435 and 0.0472 and polymorphic PVs (magenta, N = 400) have an average area of 0.0981 ± 0.0429 μm2 with a 95% confidence interval between 0.0939 and 0.102. The differences in the area are statistically significant (ANOVA; p-value < 0.0001). Scale bars: A 5 μm and 1 μm (ROIs), B 5 μm

Although a STED microscopy-based approach clearly allows the resolution of individual organelles as small as PVs, the distinctly lower axial resolution remains limiting for the three-dimensional rendering of organelles. Therefore, to push the boundaries of resolution and to further characterize PV morphology, we employed single-molecule localization microscopy (SMLM) [47,48,49].

Giardia PVs in trophozoites were loaded with a 10-kDa Dextran-Alexa Fluor 647 fluid phase marker with a high degree of photostability to survive repeated cycles of photoactivation and excitation in SMLM experiments [50, 51]. After the acquisition, images were reconstructed using the ImageJ plugin ThunderStorm which performs signal centroid calculation, image reconstruction and output [52, 53]. Dextran uptake in PVs was confirmed using conventional widefield microscopy (Fig. 3A). STORM image reconstruction shows the subcellular distribution of the fluorescent marker and defines individual organelle lumina (Fig. 3B, representative image of N = 10 cells). A closer inspection revealed the presence of morphologically distinct endocytic organelles as previously observed in our FIB-SEM and STED datasets (Fig. 3B, ROI and Additional file 6: Video S3). We again used the supervised ML-assisted algorithm in ilastik to classify the different morphologies. After a pixel segmentation routine, we performed object classification using supervised ground truth training on subsets of organelle images. Three categories of PVs were defined: spherical, tubular and polymorphic (Fig. 3C). To test whether the morphological categorization was consistent with categorization based on organelle volume, we calculated the average lumina volumes of > 4000 organelles from the three PV categories. ANOVA testing of organelle volumes for vesicular (0.00507 ± 0.00336 μm3, N = 1989, 95% confidence interval: [0.00492, 0.00522] μm3), tubular (0.0103 ± 0.00925 μm3, N = 838, 95% confidence interval: [0.00967, 0.0109] μm3) and polymorphic (0.0227 ± 0.0214 μm3, N = 1494, 95% confidence interval: [0.0216, 0.0238] μm3 organelles confirmed statistically significant (p < 0.0001) morphological differences (Fig. 3D and summarized in Additional file 7: Table S1).

Fig. 3
figure 3

Super-resolution imaging of Giardia lamblia PVs by stochastic object reconstruction microscopy (STORM). A Widefield microscopy-based imaging of a Giardia trophozoite loaded with 10-kDa Dextran-Alexa Fluor 647. B Reconstruction of single-molecule events using the Fiji plugin Thunderstorm. As with STED imaging, different PV morphologies are observed. C PVs were segmented (N = 10 cells) with the help of ilastik, and volumes were calculated and plotted in (D). Spherical PVs (green, N = 1989) present an average volume of 0.00507 ± 0.00336 μm3 with a 95% confidence interval between 0.00492 and 0.00522 μm3, tubular PVs (blue, N = 838) present an average volume of 0.0103 ± 0.00925 μm3 with a 95% confidence interval between 0.00967 and 0.0109 μm3 and polymorphic PVs (magenta, N = 1494) present an average volume of 0.0227 ± 0.0214 μm3 with a 95% confidence interval between 0.0216 and 0.0238 μm3. The differences in the area are statistically significant (ANOVA; p-value < 0.0001). Based on these and previously shown data, the renaming of PVs to peripheral endocytic compartments (PECs) is proposed. Scale bars: 5 μm and 1 μm (ROI)

Taken together, the data generated using three distinct imaging techniques clearly demonstrate PV heterogeneity which may be linked to distinct functions and/or maturation states in this unique endocytic system. To reflect this novel finding and taking into account that these endocytic and peripherally localized organelles are neither proper vesicles nor canonical vacuoles, we propose renaming PVs to peripheral endocytic compartments (PECs).

Comparative analysis of endocytic and secretory organelles in Giardia, Spironucleus sp. and T. foetus

Giardia spp. have evolved a unique cell architecture including a dedicated organelle for attachment to the small intestinal epithelium—the ventral disc (VD) [35, 54]. In turn, this innovation defines a distinct dorsal-ventral as well as antero-posterior polarization of the flagellated trophozoite, marked by swimming directionality. PVs/PECs localize exclusively to the dome-shaped dorsal parasite PM except for a small circular patch at the centre of the VD called the bare zone [26, 27]. The result is a maximally decentralized architecture of the Giardia endocytic system forming a single-layer interface of what we now appreciate as 3 morphologically distinct organelle classes between the cell exterior and the cytoplasm/ER [27, 30]. We asked whether this type of decentralized sub-PM localization and polymorphic morphology of endocytic compartments was also represented in other tractable, phylogenetically related members of the Diplomonadida as well as more distant metamonada lineages which do not have a VD, i.e. Spironucleus vortens and Spironucleus salmonicida and the parabasalid Tritrichomonas foetus.

The diplomonads S. vortens and S.salmonicida are amongst the closest tractable relatives of G. lamblia that can be grown axenically under similar conditions [55,56,57,58,59]. Their endocytic compartments and machineries are partially characterized, with some reports of large vacuolar structures detected by electron microscopy in trophozoites [60, 61]. Unlike Giardia, both species lack dorso-ventral polarization but display a distinct antero-posterior axis. Putative endocytic organelles in S. vortens have been detected by fluorescence microscopy of live and fixed cells after incubation with fluorophore-coupled dextran [27]. To further investigate these endocytic compartments, we incubated S. vortens and S. salmonicida trophozoites with a 10-kDa Dextran-TexasRed fluid phase marker (Fig. 4). In stark contrast to the distinctly arrayed PV/PEC labelling seen in Giardia lamblia (Fig. 4A), confocal microscopy revealed the presence of several dispersed labelled organelles in both S. vortens (Fig. 4B) and S. salmonicida (Fig. 4C). Spironucleus spp. display several relatively large globular membrane compartments, similar to those observed in well-characterized model organisms lacking a fixed subcellular localization [5, 62]. While S. salmonicida endocytic compartments localize mostly at the cell periphery (Fig. 4C), S. vortens organelles present both peripheral and central localizations (Fig. 4B and Additional file 8: Video S4). We also assessed endosome morphology in T. foetus using the same labelled dextran-based approach. Similar to Spironucleus species, T. foetus presents an antero-posterior axis but no attachment organelle nor dorso-ventral polarization. Similar to Spironucleus spp., T. foetus accumulated the endocytosed fluid phase maker in several globular endocytic compartments (Fig. 4D) consistent with previous reports on vacuolar structures identified by electron microscopy [63].

Fig. 4
figure 4

Uptake of fluorescently labelled dextran in metamonada and discoba members: G. lamblia, S. vortens, S. salmonicida and T. foetus after 30 min. A G. lamblia cells present endocytic compartments spread in the cell periphery, unresolved by conventional light microscopy. On the other hand, B S. vortens, C S. salmonicida and D T. foetus present vesicular endocytic compartments. Scale bars: 5 μm (full cells) and 1 μm (ROI)

Taken together, these data show how, in closely related protozoa lacking dorsal-ventral polarization and a dedicated attachment organelle, endocytic organelles appear to have no specific localization. This lends support to the notion that PV/PEC organelle architecture is intimately associated to the emergence of the VD, both structures as adaptations to the mammalian small intestine niche [27].

To visualize and measure the morphological parameters of Spironucleus and T. foetus endocytic compartments, we performed 2D-STED imaging and transmission electron microscopy (TEM). S. vortens cells loaded with 10-kDa Dextran-Alexa Fluor 594 showed accumulation of the fluid phase marker in roughly spherical organelles (Fig. 5A). Labelled endocytic vacuoles have an average diameter of 468 ± 206 nm (95% confidence interval [421; 515] nm, N = 10 cells) (Fig. 5A, violin plot). Volumetric rendering of 3D reconstructed optical sections documents the uniformly globular morphology of these organelles (Additional file 9: Video S5). TEM analysis revealed an ellipsoid shape of endocytic organelles in S. vortens with an average maximal diameter of 844 ± 335 nm with a 95% confidence interval between [763;905] nm (Fig. 5B, violin-plot). The dimensions measured in TEM represent those of the membrane-delimited organelle. In contrast, the dimensions measured by STED represent a projection of the fluid phase marker distribution within the available organelle lumen. The fact that the former (844 ± 335 nm) is larger than the latter (468 ± 206 nm) indicates that these organelles may contain additional cargo which prevents the endocytosed fluid phase marker to distribute in the complete compartment volume delimited by the organelle membranes. TEM investigation in S. salmonicida cells (Additional file 10: Fig. S4A) showed the presence of small globular vacuoles (V) close to the PM (Additional file 10: Fig. S4B) with an average diameter of 205 ± 62.6 nm (N = 114) with a 95% confidence interval of 193 and 217 nm (Additional file 10: Fig. S4E). These vacuoles are smaller than the ones found in S. vortens (Additional file 10: Fig. S4C, D, F; p-value < 0.0001). In these conditions, neither coated vesicles nor a stacked Golgi apparatus could be documented in S. vortens or S. salmonicida.

Fig. 5
figure 5

Super-resolution imaging of S. vortens endocytic compartments with STED and transmission electron microscopy (TEM). A Following incubation with dextran-TexasRed for 30 min, S. vortens display elongated endocytic compartments with an average diameter of 468 ± 206 nm with a 95% confidence interval between 421 and 515 nm (upper violin plot, N = 10 cells). B TEM imaging detects several endosome-like vacuoles (V) throughout the cell, with an average diameter of 844 ± 335 nm with a 95% confidence interval between 763 and 905 nm (lower violin plot). Measurements were done manually. The endoplasmic reticulum (ER) and a nucleus (N) are also highlighted in the images. Scale bars: 5 μm (full field of view) and 1 μm (ROIs)

2D-STED analysis of T. foetus cells incubated with 10-kDa Dextran-Alexa Fluor 594 revealed a roughly circular distribution of the marker within endocytic vacuoles (Fig. 6A and Additional file 11: Video S6) with an average maximal diameter of 517 ± 251 nm (95% confidence interval [455; 580] nm, N = 10 cells) (Fig. 6A, violin plot). TEM imaging revealed the presence of two distinct classes of endosome-like vesicles (Fig. 6B) based on the electron density of the lumen. Low-density vesicles were identified both at the cell periphery and in central areas termed vacuoles (V); vesicles of higher electron density were previously identified as digestive vacuoles (DVs) [63] and contain structured material and membranes. Analysis of TEM micrographs showed that DVs are significantly larger than vacuoles, with an average diameter of 764 ± 203 nm (N = 50) (95% confidence interval [707; 822] nm). Vacuoles in turn have an average diameter of 246 ± 100 nm (N = 153) in (95% confidence interval [230; 262] nm) (Fig. 6B, violin plot). Stacked Golgi organelles are abundant in TEM micrographs of T. foetus trophozoites, as documented previously [64] (Additional file 12: Fig. S5A). Consistent with a more canonical architecture of the membrane trafficking system in T. foetus, coated vesicles were observed in the cytosol particularly in the vicinity of Golgi stacks (Additional file 12: Fig. S5B) [63, 65, 66]. These vesicles averaged a diameter of 58.4 ± 13.1 nm (N = 128) (95% confidence interval [56.1; 60.7] nm) corresponding to the size range of clathrin-coated vesicles (CCVs) [67]. In our ultrastructure observations, we did not detect multivesicular bodies nor vacuoles containing intra-luminal vesicles, neither in Spironucleus spp. nor T. foetus.

Fig. 6
figure 6

Super-resolution imaging of T. foetus endocytic compartments with STED and TEM. A 2D-STED analysis of T. foetus cells loaded with 10-kDa Dextran-Alexa Fluor 594 illuminates endosomes with globular structures at an average diameter of 517 ± 251 nm with a 95% confidence interval between 455 and 580 nm (upper violin plot, N = 10 cells). B TEM investigation of T. foetus cells reveals two kinds of endosome-like vesicles. Low electron-density vesicles found both at the cell periphery and centre are deemed vacuoles (V) with an average diameter of 246 ± 100 nm (N = 153) with a 95% confidence interval of 230 and 262 nm. Vesicles of higher electron density and with noticeable content are digestive vacuoles (DVs). DVs are larger than vacuoles, with an average diameter of 764 ± 203 nm (N = 50) with a 95% confidence interval of 707 and 822 nm (lower violin plot). Scale bars: 5 μm (full field of view) and 1 μm (ROIs)

Finally, to probe the dynamics of endocytic compartments in G. lamblia, S. vortens, S. salmonicida and T. foetus, cells were exposed to 10-kDa Dextran-TexasRed for 5, 10, 20 or 30 min, fixed chemically, and imaged by confocal microscopy (Fig. 7). The number of G. lamblia PV/PECs labelled with the fluid-phase marker increased over time, with the label accumulating strictly at the cell periphery (Fig. 7A). In contrast, endocytic compartments in S. vortens were first visualized at the PM and were then observed at more central locations of the cell at later time points. Given the overall increase in fluorescent intensity and the motile nature of these organelles, it appears there is a constant uptake of dextran over the analysed period (Fig. 7B). In these conditions, S. salmonicida and T. foetus vacuoles both appear diffused within the cell cytoplasm although S. salmonicida organelles show a decrease and then marked increase in dextran content (Fig. 7C) while T. foetus organelles present irregular fluctuations in dextran content (Fig. 7D).

Fig. 7
figure 7

Time course on fluorescent dextran uptake in selected metamonada and discoba specimens: G. lamblia, S. vortens, S. salmonicida and T. foetus after 30 min. A PV/PECs in G. lamblia start acquiring external material right after 5 min of incubation with dextran. Over 30 min of dextran incubation, the endocytic marker does not leave PV/PECs. In B S. vortens, C S. salmonicida and D T. foetus, dextran is up taken in the small vesicles. In S. vortens, these small vesicles tend to agglomerate at the centre of the cell while in S. salmonicida, these vesicles seem to stay peripheral. In T. foetus, vesicles bearing dextran from the periphery of the cell also migrate to the interior of the cell. All scale bars: 5 μm

Pan-eukaryotic searches for endocytic markers CHC and CLC reveal the loss of a bona fide CLC within the fornicata lineage and the emergence of putative CLC analogues

Previously, we established that Giardia clathrin heavy chain (GlCHC) associates to discrete static foci at the dorsal PM of trophozoites, in close proximity to PV/PECs. Furthermore, GlCHC strongly interacts with a putative albeit highly diverged Giardia clathrin light chain homologue previously named GlCLC. Native co-IP experiments demonstrate how CHC is invariably associated to CLC, consistent with the lack of a detectable cytoplasmic pool of CHC and suggesting that virtually all known CHCs and CLCs exist in a complex [26, 27]. However, whether GlCLC is truly a CLC direct homologue (i.e. orthologous) is unclear, as is the prevalence of this protein in other members of the lineage Fornicata. To investigate the occurrence of both CHC and CLC orthologues in selected eukaryotic lineages, we employed protein homology searches based on hidden Markov models (HMM) [68] using as query an alignment of canonical and documented CHC or CLC sequences from several protozoa and metazoan species (Additional file 13: Tables S2 and S5) [69,70,71,72,73,74,75]. In this search, we considered assembled read data from RNA-seq experiments (transcriptomics) as reliable as genomic sequence data [76]. In this case, we used the reference CHC or CLC sequences and performed tblastn searches. Nucleotide sequences from each reliable hit (lowest e-value) were translated and subjected to a reciprocal blast-p analysis to validate protein identity. We found CHC homologues in all selected genomes and transcriptomes we searched, highlighting the likely essential nature of CHC (Fig. 8A and Additional file 13: Tables S4 and S5).

Fig. 8
figure 8

Homology search for bona fide CHC and CLC reveals the key importance of CHC and patchy conservation of CLC. A CHC and CLC homology searches demonstrate the conservation of CHC orthologs in all analysed lineages. However, several lineages appeared to have lost a readily detectable bona fide CLC including all selected Fornicata, certain Chromerids [77], some ciliates [78] and members of the Hemimastigophora. B Ab initio in silico protein modelling of G. lamblia ACLC (formerly GlCLC/Gl4529), T. brucei CLC, T. foetus CLC, M. exilis CLC, D. discoideum and H. sapiens CLC sequences using AlphaFold, the current standard in in silico protein structure modelling. TM-align and RMSD values were calculated showing a close structure analogy between the predicted structures (table)

GlCHC is a clearly divergent ortholog compared to its counterpart in eukaryotic model organisms, with only 24% amino acid identity to human CHC [25]. A domain analysis of selected CHC sequences (Additional file 14: Fig. S6 and Additional file 15: Table S10) reveals that GlCHC contains fewer α-helical domains than other analysed CHC sequences, further highlighting its divergence. We also performed an in-depth search for the CHC triskelion uncoating “QLMLT” motif which we documented previously to be missing in Giardia [27, 79, 80]. Notably, this motif appears to be only present in Metazoa and in the closely related Filastera and Choanoflagellata [81,82,83]. In Fungi, only a partial “L(M)TL” motif was identified, and we were unable to detect a conserved uncoating motif in CHC sequences of members of the Archaeplastida, Amoebozoa or SAR supergroups (Additional file 16: Fig. S7).

In stark contrast to CHC, the search for bona fide CLC sequences did not retrieve any reliable predictions in available genomes and transcriptomes from species of the Fornicata lineage, including the lineages Hexamitidae, Retortamonas and Carpediemonas-like organisms [59, 75, 84,85,86]. Importantly, this search did not return the putative, highly diverged GlCLC [27]. There are documented CLC orthologues in members of the Discoba, such as Trypanosoma brucei CLC (Tb927.10.14760) [87] and in the parabasalid Trichomonas vaginalis (TVAG_29749) [88, 89]. Furthermore, we readily identified a CLC homologue in T. foetus (gene accession OHT14195.1, forward HMMer e-value of 1.00E−26 and reverse Blastp e-value of 2.00E−11, returning the human CLC homologue) (Fig. 9A). Therefore, while bona fide CLC orthologues can be readily identified in Discoba and Euglenozoa members and Preaxostyla—as in the metamonad Monocercomonoides exilis [73]—no sequence could be found amongst the members of the Fornicata lineage. Furthermore, we were unable to identify a bona fide CLC sequence within the newly documented transcriptome of Hemimastigophora [90]. Beyond Fornicata, we failed to identify bona fide CLC in chromerids such as Cryptosporidium parvum and Babesia bovis [77] and in some ciliate lineages, such as Tetrahymena thermophila and Paramecium caudatum [78].

Fig. 9
figure 9

GlCHC and GlACLC foci associate to Giardia PV/PECs with different stoichiometries. A An epitope-tagged GlCHC reporter was used to localize the foci of GlCHC deposition at the cell periphery beneath the PM, by STED microscopy. B These foci are well-resolved with STED microscopy (left graph). C The same was done with an epitope-tagged GlACLC reporter, and D the foci resolved as in B (right graph). E The foci were segmented with ilastik. Areas were determined and diameters calculated in an automatic procedure assuming spherical geometry. GlCHC foci present an average diameter of 134 ± 36.6 nm (N = 4524) with a 95% confidence interval between 132 and 135 nm. GlACLC foci average a diameter of 159 ± 48.8 nm (N = 984) with a 95% confidence interval between 156 and 162 nm. GlACLC foci are larger than GlCHC in significant statistical manner (p-value < 0.0001, t-Student test). F GlCHC foci associate in different stoichiometry to different classes of PV/PECs. Spherical PV/PECs are either not associated to clathrin or with just one focus: a mean of 0.488 ± 0.159 foci per spherical with a 95% confidence interval between 0.403 and 0.573 foci per spherical organelle. Tubular PV/PECs associate with one focus of GlCHC: a mean of 1.15 ± 0.287 with a 95% confidence interval between 0.994 and 1.3 foci per tubular organelle. Polymorphic PV/PECs associate with 3 or more GlCHC foci with an average of 3.85 ± 1.14 with a 95% confidence interval between 3.25 and 4.46 foci per polymorphic organelle. These distributions are statistically significant with a p-value < 0.0001 (ANOVA analysis)

Given that GlCLC’s predicted 3D structure is reminiscent of CLCs (Zumthor et al. [27]) but it could not be retrieved as related to a bona fide CLC, with its only known orthologue found in Giardia muris (Fig. 8A), GlCLC was further analysed using the HHPred suit, in the attempt to find distantly related non-Giardia sequences [91]. This search retrieved no robust prediction for a non-Giardia sequence (Additional file 13: Table S9). Given that the degree of divergence is such that no reliable claim to orthology can currently be supported and no orthologue for GlCLC can be found outside the Giardia genus, we propose the renaming of GlCLC to Giardia lamblia analogous to clathrin light chain—GlACLC—as a CLC structural analogue acquired and retained in the last Giardia common ancestor (LGCA). This appears to correlate with the loss of a bona fide CLC with the last Fornicata common ancestor (LFCA). To test the extent of environmental pressure on this protein family’s evolution, we calculated synonymous vs non-synonymous mutation ratios (ω = ks/kn) for GlACLC homologues (Additional file 17: Fig. S8). Interestingly, known sequences for all Giardia isolates present a ω < 1 which indicates that current sequences are not under selective pressure to evolve. To further investigate the structural analogy of GlACLC to canonical CLCs, we performed in silico modelling of its C-terminal domain using the new standard in ab initio protein structure modelling—AlphaFold—based on deep-learning neural networks [92, 93] (Fig. 8B and Additional file 20: Table S11). Template modelling score independent of sequence (Tm-align) and root mean square deviation (RMSD) [94, 95] values provide substantial evidence for the structural analogy of GlACLC and canonical CLCs, in line with previous observations [27]. The newly predicted structures for GlACLC have a stronger resemblance to the predicted structure of a mammalian clathrin light chain [96]. Altogether, the presented in silico data strongly suggest GlACLC to be a structural analogue of CLC.

Finally, we asked whether de novo acquisition of a CLC analogue with a divergent sequence, but the preservation of structural features had occurred independently in other Diplomonadida lineages. In a first approach, we selected S. salmonicida the closest genetically tractable and sequenced relative to Giardia [59, 97] in which a bona fide CLC can not be detected. An epitope-tagged variant of the ca. 210 kDa S. salmonicida CHC orthologue (SsCHC-3xHA, ORF Ss50377_14164) distributes in a punctate pattern throughout the trophozoite cytosol (Additional file 18: Fig. S9A), reminiscent of GlCHC focal assemblies (Additional file 18: Fig. S9B and Additional file 19: Video S7). SsCHC-3xHA was utilized as an affinity handle to define a putative SsCHC interactome in a single native co-IP and protein identification experiment (Additional file 18: Fig. S9C, Additional file 20: Table S11). Amongst several endocytosis-related proteins (Ss-dynamin, Ss-β-adaptin, Ss-calmodulin and SsSec7), one ORF namely, Ss50377_11905, was found to be prominently pulled down and contains several coil-coil domains, at a predicted weight of 39 kDa. The 150-amino acid C-terminus of the protein was modelled in AlphaFold and superimposed with CLC structures (Additional file 18: Fig. S9E). TM-align values within structure similarity (0.5 or above), and RMSD values of circa 5-6 Å suggest Ss50377_11905 may be a S. salmonicida structural CLC analogue (Additional file 21: Fig. S10). Using the Ss11905 sequence, we retrieved a putative orthologue only in the available transcriptome of the related diplomonad Trepomonas sp. [56, 57], namely, TPC1_16039 (forward tblastn e-value of 1E−5 and reverse blastp e-value of 4E−12), and in no other selected fornicate lineage (Additional file 22: Table S12).

GlCHC foci associate with different classes of giardial PV/PECs with varying stoichiometry

Having established the presence of at least three different PV/PEC morphologies, presumably corresponding to different cryptic organelles, and having established that the GlACLC is a unique Giardia protein, we wanted to better understand the association of the clathrin complexes with the various PV/PECs and if there is any correlation between PV/PEC morphology and clathrin assemblies. We first wanted to confirm whether GlCHC and GlACLC are consistently found together in foci. To address this question, we used STED microscopy to investigate epitope-tagged GlCHC (GlCHC-HA) deposition at distinct foci at the dorsal PM and the cell periphery, consistent with PV/PECs location (Fig. 9A). Segmentation of foci using a ML-assisted ilastik tool allowed to determine the dimensions of GlCHC at an average diameter of 134 ± 36.6 nm (N = 4524) (95% confidence interval [132; 135] nm) (Fig. 9B). Similar to GlCHC, the subcellular distribution of epitope-tagged GlACLC-HA showed an identical pattern consistent with its demonstrated direct interaction with GlCHC (Fig. 9C) [27]. Segmentation of foci using a ML-assisted ilastik tool determined the dimensions of GlACLC foci at an average diameter of 159 ± 48.8 nm (N = 984) (95% confidence interval [156; 162] nm) (Fig. 9D). Notably, the average size of GlACLC foci is larger than that of GlCHC foci (p < 0.0001, t-Student test) (Fig. 9E).

Using STED microscopy, we further determined the number of GlCHC foci showing signal overlap with the three classes of dextran-Texas Red loaded PV/PECs (Fig. 9F). By calculating the degree of the signal overlap between GlCHC foci and PV/PEC lumina, we determined that spherical PV/PECs are associated with at most one GlCHC focus with an average of 0.488 ± 0.159 foci per spherical PV/PEC (95% confidence interval [0.403; 0.573]). Tubular PV/PECs associated with at least one GlCHC focus with an average of 1.15 ± 0.287 foci per tubular PV/PEC (95% confidence interval [0.994; 1.3]). Polymorphic PV/PECs associated with three or more GlCHC foci with an average of 3.85 ± 1.14 foci per PV/PEC (95% confidence interval [3.25; 4.46]). Taken together, we find a directly proportional and statistically significant ratio of clathrin foci to PV/PEC size and type (ANOVA; p-value < 0.0001). This is in line with the possibility that PV/PEC morphological heterogeneity is correlated with organelle functional diversity, as measured by association to clathrin assemblies consisting of GlCHC and GlACLC.


The Giardia endocytic organelle system consists of three classes of membrane compartments

Subsequent to ingestion and excystation, Giardia trophozoites attach to the intestinal lumen, proliferating and encysting on localized foci throughout the mucosa of the small intestine [98]. Nutrients required for this propagation are taken up from the environment through PV/PEC-mediated endocytosis of fluid phase and membrane-bound material [24, 27, 30, 99,100,101,102]. Despite the essential nature of these endocytic organelles, complete resolution of the ultrastructure of the Giardia endocytic pathway remains unsolved. To address this, we performed an ultrastructural investigation of G. lamblia endocytic compartments to obtain a nanometric view of their morphology as defined by their membrane as well as the lumen accessible to fluid phase markers in labelling experiments [27, 30].

We began by dissecting an entire G. lamblia trophozoite using scanning electron microscopy and focused our analysis on PVs. These structures were segmented and rendered in three dimensions. Using this method unambiguously detected at least two distinct classes of PV morphologies, with some being obviously globular in shape while others presenting a more tubular nature. After expanding our analysis of PVs to super-resolution light microscopy methods STED and STORM [103], we determined that PVs are present in three discernible morphologies: spherical, tubular and polymorphic. Thus, we proposed the renaming of these organelles into peripheral endocytic compartments (PECs).

Compared to endosome-like vacuoles in Carpediemonas-like organisms (CLOs) [104,105,106] and large vesicular endosome-like structures observed in S. salmonicida and S. vortens and the more distantly related Parabasalia member, T. foetus, specific and complete remodelling of endosomes has occurred in the Giardia genus. T. foetus, except for the presence of endosome-like vesicles, presents digestive vacuoles and a stacked Golgi apparatus (Fig. 10). Coated vesicles, likely CCVs, are observed near the T. foetus Golgi apparatus and the PM. In our analysis, we could not confirm fluid phase material uptake through the cytostome present in Spironucleus sp. [60] and dextran accumulated in spherical vesicles of different dimensions and unknown origin, similar to endosomes. Figure 10A, B summarizes the results of our comparative analysis and highlights the unique endocytic system in Giardia where, unlike related species and other excavates, PV/PEC-mediated uptake is restricted to the dorsal side of the cell [27, 107] while the ventral side is deputed to attachment to host structures. Interestingly, endosome and lysosome tubulation has been documented in macrophages [42, 43] and are linked with different physiological states of the organelles and subsequent function in the cell—such as prompting the cell for phagocytosis. This allows for the hypothesis that different kinds of PV/PECs correspond to different stages in organelle maturation/function. In line with this, we provide evidence for different stochiometric associations of CHC foci with different kinds of PV/PECs, although it is possible polymorphic PV/PECs are associated to more foci simply due to their larger area (Fig. 2C).

Fig. 10
figure 10

Endosome-like organelle models in Giardia, Spironucleus and Tritrichomonas and proposed evolution of CLC. A Simplified cartoons of the endocytic and secretory pathway in Giardia lamblia (Giardiinae), Spironucleus sp. (Hexamitidae) and Tritrichomonas foetus (Parabasalia). PV/PECs are a hallmark of the Giardia lineage, while more canonical vesicular endosomes are present in both the Spironucleus lineage and Parabasalia. B Simplified evolutionary model for bona fide CLC. The last eukaryotic common ancestor possessed a bona fide CLC which was lost at the last Fornicata common ancestor. In at least two derived lineages—Giardinae and Spironucleus spp.—de novo protein analogue to CLC was acquired independently

FIB-SEM ML-based analysis of a full trophozoite also yielded unprecedented views of large numbers of ca. 80 nm SVs which had been previously albeit anecdotally reported [27, 108,109,110]. Based on their electron-dense surface, SVs appear coated and are not related to CHC foci at the PV/PECs-PM interface [26, 27]. Given the absence of a specific marker for these compartments, their identity remains elusive although their coated appearance suggests protein trafficking routes (COPI or COPII) may be involved [25, 111, 112]. However, recent reports in Giardia of immuno-EM detection of peroxisome-like proteins in small dense ca. 100 nm vesicles [113] provide an open playing field for the definitive identification of these compartments.

G. lamblia possesses a highly divergent clathrin heavy chain and a newly acquired clathrin light chain analogue

We performed an in-depth search for CHC homologues within excavates and other key eukaryotic groups. We found that CHC is conserved in all of these organisms, underlining the vital role of CHC in eukaryotic organisms. The sequence divergence of the giardial CHC protein is reflected in an overall decrease in the number of α-helical domains which are essential for the formation of the triskelion leg, and hence necessary for coat assembly [69]. Thus, the reduction in α-helical domains during GlCHC evolution may have led to a lower propensity of GlCHC forming triskelion assemblies and membrane coats. So far, none of the many attempted methods to detect GlCHC in association with small vesicles has been able to show anything other than an exclusive focal localization at PVs/PEC membrane interfaces [27]. Also, the GlCHC protein does not contain the C-terminal uncoating motif “QLMLT” nor is this motif present in the CHC homologues of any diplomonad. In fact, this motif appears to be only present in Metazoa and in the closely related Filastera and Choanoflagellata [81,82,83] despite the documented ability to form and uncoat bona fide CCVs in some protozoa [114]. In Fungi, only a partial “L(M)TL” motif was identified. We could not detect a conserved uncoating motif in CHC sequences of members of the Archaeplastida, Amoebozoa or SAR supergroups (Additional file 17: Fig. 7). Taken together, this data indicates the uncoating QLMLT motif is apparently specific to and likely and invention of the Holozoa lineage. This observation points to as yet uncharacterized uncoating mechanisms are present in other lineages. For example, clathrin-mediated endocytosis is essential in the parasitic protist Trypanosoma brucei and CCVs have been documented in this organism [70, 71, 114, 115]. Clathrin and other coat proteins associated with CCVs need to be recycled. While HSC70 is documented in T. brucei and likely involved in clathrin uncoating [80], no bona fide uncoating motif has been documented [71, 87, 114].

In contrast to GlCHC, the evolution of the previously identified putative GlCLC/Gl4259 protein presents a different and surprising natural history. This protein was identified as the strongest interactor of GlCHC [27] and is present in all sequenced Giardia lineages. GlCLC/Gl4259 has no measurable sequence conservation but a high degree of structural similarity to bona fide CLCs, warranting its proposed renaming to GlACLC. Aside from the Giardia genus, we were unable to identify homologues for GlACLC in any other eukaryotic taxa, nor could we find any orthologues of CLC in any available Fornicata genome/transcriptome sequence, suggesting that the last Fornicata common ancestor (LFCA) lacked a canonical CLC. Taken together, the available data is currently insufficient to decide between two mutually exclusive evolutionary scenarios: (a) secondary loss of a canonical CLC in the last fornicate common ancestor, with acquisition of a structurally and functionally related GlACLC, or (b) massive sequence divergence of the original, rendering it undetectable by even our most sensitive methods. This latter scenario could have been driven by significant changes in function particularly in Giardia where once-dynamic membrane coating machinery has evolved to become a static structural element supporting interfaces between the plasma membrane and the endocytic system. The discovery of a strong interactor of CHC in the closely related S. salmonicida Ss11905, with structural similarity to GlACLC as well as to bona fide CLCs is consistent with both scenarios. Notably, this protein neither retrieves GlACLC nor CLCs in BLAST searches, leaving no evidence of direct homology. By contrast, robust predictions for CLC homologues were made for members of the Preaxostyla, Discoba and Parabasalia lineages. Although scenario 2 is possible there is no evidence supporting it, and until some emerges, we must favour scenario one, loss and emergence of a convergent structural analogue, which is supported by the evidence in hand. Other lineages appear also to have lost a bona fide CLC, like C. parvum and T. thermophila [77, 78], but perhaps similar investigations to ours of CHC may identify CLC analogues/divergent homologues. Taken together, this data suggests that the constraints on the CHC primary structure are higher than on CLC even after massive changes in clathrin coat function with demonstrated complete losses in some protists. Members of the Giardia genus as well as S. salmonicida have no identifiable bona fide CLC, yet at least the giardial GlACLC has retained its function as a CHC interacting partner.


Our data provide a robust understanding of Giardia, Hexamitidae members and Tritrichomonas foetus endocytic pathway organellar ultrastructure. Contrary to Spironucleus or Tritrichomonas and other excavates, Giardia underwent a complete remodelling of its endocytic machinery. Our investigation revealed its organelles to be polymorphic in nature, justifying the proposed name change to peripheral endocytic compartments. Furthermore, the analysis of GlCHC sequences highlights its divergence which is likely due to a massive reorganization of the endocytic pathway in these species, while the origin and evolution of CLC structural and to some extent functional homologs in Giardia (GlACLC) and in certain Hexamitidae members (S. salmonicida and Trepomonas sp. PC1) remains uncertain.


Cell culture and transfection

Giardia intestinalis strain WB (clone C6; ATCC catalogue number 50803) trophozoites were grown using standard methods as described in Morf et. al. [116]. Episomally transfected parasites were obtained via electroporation of the circular pPacV-Integ-based plasmid prepared in E. coli as described in Zumthor et al. [27] Transfectants were selected using Puromycin (final conc. 50 μg ml−1; InvivoGen). S. vortens and S. salmonicida were cultured as described before [58, 59]. S. salmonicida was transfected using a modified PAC vector and selected with Puromycin (final conc. 50 μg ml−1; InvivoGen) [97]. T. foetus was axenically grown also as described [63].

Construction of expression vectors

S. salmonicida CHC sequence (SS50377_14164) was amplified with the primers ATATTTAATTAAGGCGGATCTATAGTTTCTTGGAATACTAAAATAGGA (forward) and TATGCGGCCGCCACCAGTTATCAGCGGGTGCC (reverse) containing a MluI and a NotI rectrictyion site, respectively. The genomic sequence amplified contained a 5′ UTR region of 179 bp which encodes a putative promoter. The genomic fragment was inserted in the previously described vector pSpiro-PAC-3xHA-C [97].

Focused ion beam scanning electron microscopy (FIB-SEM) of a full Giardia trophozoite and image analysis

Wild-type Giardia lamblia trophozoites were subject to high-pressure freezing and processed as established in [27]. Ion milling and imaging were performed in a Auriga 40 Crossbeam system (Zeiss, Oberkochen, Germany) using the FIBICS Nanopatterning engine (Fibics Inc., Ottawa, Canada) following the aforementioned established protocol. The pixel size was set to 5 nm, obtaining isotropic imaging. Alignment of the dataset was performed resorting to the ImageJ plugin Sift [53]. Image segmentation was done using the semi-autonomous algorithm ilastik [38]. The routine of pixel and object classification are used. Algorithm training was performed in a small representative region of the dataset which was then applied to the complete dataset. Imaris (Bitplane AG) was used for three-dimensional rendering and volume measuring. Diameters (d) of spherical PV/PECs were measured based on radius (r) calculation after volume (Vspherical) determination in Imaris. Thus, \({V}_{\mathrm{spherical}}=\frac{4}{3}\pi {r}^3\leftrightarrow r=\sqrt[3]{\frac{3{V}_{\mathrm{spherical}}}{4\pi }}\leftrightarrow d=2r\). Diameter (d) and length (L) of tubular PVs/PECs (and assuming a cylindric shape) were calculated after lateral area (A) and volume (V) determination in Imaris. Thus, A = 2πrL + 2πr2. As L r, it can be simplified to A = 2πrL. With \(V=\pi {r}^2L\leftrightarrow r=\frac{2V}{A}\) and \(L=\frac{A}{2\pi r}\). Finally, d = 2r. Graphical representations were shown in their log10 for easier comparison.

Transmission electron microscopy analysis of Giardia lamblia, Spironucleus spp. and Tritrichomonas foetus cells and analysis

G. lamblia, S. vortens, S. salmonicida and Tritrichomonas foetus samples were subject to high-pressure freezing and processed as we previously established [27, 117]. Samples were imaged in a FEI CM100 transmission electron microscope. The pixel size was assigned to 0.8 nm. Tiles were obtained automatically after the determination of the focal point. Tiles were aligned with TrakEM2 [36].

Immunofluorescence assays

Chemically fixed cells for subcellular recombinant protein localization were prepared as previously described [118]. HA-epitope tagged recombinant proteins were detected using a rat-derived monoclonal anti-HA antibody (dilution 1:200, Roche) followed by a secondary anti-rat antibody coupled to AlexaFluor 488 fluorophores (dilution 1:200, Invitrogen). Samples were embedded in Vectashield (VectorLabs) or Prolong Diamond Mounting medium (Invitrogen) containing 4′,6-diamidino-2-phenylindole (DAPI) for nuclear staining.

Fluid-phase marker uptake

Dextran uptake assays were performed as described in [27, 117] using dextran 10 kDa at 2 mg/mL (Invitrogen). Coupled fluorophore was chosen based on the image technique chosen. Immunostaining was performed as described above with the exception of using only 0.02% Triton-X100 (Sigma) in 2% BSA (Sigma) for permeabilization, to prevent leakage and loss of dextran signal. Intensities were calculated with a costume-developed macro in Fiji/ImageJ [53], resorting to WEKA algorithms for segmentation [119].

Laser scan confocal microscopy (LSCM)

Imaging was performed in an inverted confocal laser scanning microscope Leica SP8 using appropriate parameters. Confocal images were subsequently deconvolved using Huygens Professional ( and analysed using Fiji/ImageJ [53].

Stimulated emission depletion (STED) microscopy

Sample preparation was done as described for LSCM. For imaging, samples were mounted in ProLong Diamond antifade reagent (Thermo Fisher Scientific). Super-resolution microscopy was performed on a LSCM SP8 gSTED 3X Leica (Leica Microsystems) using appropriate gating settings. Nuclear labelling was omitted due to possible interference with the STED laser. A pulse depletion laser of 775 nm at 100% strength was used to deplete the signal coming from samples using the fluorophore Alexa Fluor 594. The signal from samples containing Alexa Fluor 488 was depleted with the depletion laser line 592 nm at 50% strength. Pinhole was kept at 1 AU. Images were deconvolved using Huygens Professional ( After deconvolution, the signal was segmented following a pixel and object classification routine in ilastik. Thresholding was processed in Fiji/ImageJ [53] with the respective calculation of organelle area.

Single-molecule localization microscopy (SMLM)

Cells were fixed onto a coverslip using a cytospin (6 min, 600 g). Samples were then embedded in Vectashield-based imaging medium [50]. Excess buffer was dried up, and samples were sealed. Single-molecule imaging was performed on a Leica SR-GSD 3D microscope (Leica Microsystems) as described in [120] with cylindrical lenses, in order to image the apical cell region, giving a z-depth of about 800 nm. A minimum of 100,000 events were recorded. Image reconstruction was performed with the ImageJ plugin Thunderstorm [52]. Reconstructed images were segmented following a pixel and object classification routine in ilastik [37, 38]. Thresholding and volume calculation was performed in Imaris (Bitplane AG).

Native co-immunoprecipitation of S. salmonicida CHC

Co-immunoprecipitation assays on control wild-type and transgenic S. salmonicida bearing the HA-tagged CHC were processed as previously established [27] in non-cross-linking condition agent.

Protein analysis and sample preparation for mass spectrometry (MS)-based protein identification

SDS-PAGE analysis was performed on 4–10% polyacrylamide gels under reducing conditions. Blotting was done as described in [118] using a primary rat-derived anti-HA antibody (dilution 1:500, Roche) followed by an anti-rat (dilution 1:2000; Southern Biotech) antibody coupled to horseradish peroxidase. Gels for mass spectroscopy (MS) analysis were stained with Instant blue (Expedeon) and de-stained with ultrapure water. MS-based protein identification was performed as previously reported [27].

In silico co-immunoprecipitation dataset analysis

The co-IP datasets derived from transgenic cells expressing epitope-tagged “baits” as affinity handles were filtered using dedicated control co-IP datasets generated from non-transgenic wild-type parasites to identify candidate interaction partners unique to bait-specific datasets. This was done using Scaffold4 ( Unless otherwise indicated, bait-derived co-IP data was filtered using high-stringency parameters (exclusive spectrum counts at 95-2-95, 0% FDR) and manually curated to rank putative interaction partners in a semi-quantitative fashion using ESCs as a proxy for relative abundance. Only proteins with more than 10 hits were considered. Proteins in both datasets were only considered if present 3-fold in the transgenic line versus the control. In silico analysis of hypothetical proteins was mainly carried out using BLASTp for protein homology detection ( and HHPred ( for protein homology detection based on hidden Markov model (HMM-HMM) comparisons and a cut-off at e-value < 0.05 was implemented to assign in silico annotation to otherwise non-annotated proteins of unknown function [91].

Protein structure was modelled with the ab initio modelling tool AlphaFold ( from Alphabet, powered by Google DeepMind ( deep learning neural network algorithms [92, 93]. Modelling was done via Google Colab in a Jupyter notebook environment(//

The TM-align calculation was performed online on the server Pymol (the PyMOL Molecular Graphics System, version 2.0 Schrödinger, LLC.) was used for protein structure prediction visualization, superimposing and RMSD calculation using the cealign command.

Homologue search and phylogenetic analysis and tree construction

CHC and CLC sequences were probed amongst several available genomes and transcriptomes with special focus within the fornicata members. Query protein sequences for CHC and CLC from several pan-eukaryotic representatives were obtained and aligned using MUSCLE v.3.8.31 [121] (Additional file 13: Table S2). Resulting alignments were used to generate hidden Markov models using the hmmbuild option and HHMer searches were made on all available genomes with an e-value cut-off of 0.01 [68]. Hits were considered valid if reciprocal BLASTp returned a Homo sapiens homologue with a e-value < 0.05. Transcriptome searches were carried out resorting to tBLASTn searches using the Homo sapiens and Monocercomonoides exilis respective sequences for CHC or CLC. Once a hit was found, it was translated into an amino acid sequence and was considered valid if it pulled a Homo sapiens homologue with an e-value < 0.05. All found sequences can be found in Additional file 13: Tables S2 to S9. Protein domain searches were performed at the Conservate Domain Database (CDD), through the Pfam database [122, 123]. The interPro and SMART platforms were also used for domain classification [124, 125]. Synonymous vs non-synonymous mutation ratio was calculated with an available online software ( following maximum likelihood parameters.

Statistical analysis and further used software

All data was analysed for statistical significance and plotted using the Prism 9 (Graphpad, software. Images were composed using the Affinity Designer software ( Video processing was made using Da Vinci Resolve v17.3.

Availability of data and materials

All data and materials are available. Access to raw mass spectrometry data is provided through the ProteomeXchange Consortium on the PRIDE platform [126]. Data is freely available using the project accession number and project DOI PXD020201 (


  1. Heimerl T, Flechsler J, Pickl C, Heinz V, Salecker B, Zweck J, et al. A complex endomembrane system in the archaeon Ignicoccus hospitalis tapped by Nanoarchaeum equitans. Front Microbiol. 2017;8:1–13.

    Article  Google Scholar 

  2. Kaksonen M, Roux AAA. Mechanisms of clathrin-mediated endocytosis. Nat Rev Mol Cell Biol. 2018;19:313–26.

    Article  CAS  PubMed  Google Scholar 

  3. Robinson MS. Forty years of clathrin-coated vesicles. Traffic. 2015;16:1210–38.

    Article  CAS  PubMed  Google Scholar 

  4. Naslavsky N, Caplan S. The enigmatic endosome–sorting the ins and outs of endocytic trafficking. J Cell Sci. 2018;131(13):jcs216499.

  5. Huotari J, Helenius A. Endosome maturation. EMBO J. 2011;30:3481–500.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  6. Jaiswal JK, Rivera VM, Simon SM. Exocytosis of post-Golgi vesicles is regulated by components of the endocytic machinery. Cell. 2009;137:1308–19.

    Article  PubMed  PubMed Central  Google Scholar 

  7. Radulescu AE, Siddhanta A, Shields D. A role for clathrin in reassembly of the Golgi apparatus. Mol Biol Cell. 2007;18:94–105.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  8. Poulin R, Randhawa HS. Evolution of parasitism along convergent lines: from ecology to genomics. Parasitology. 2015;142:S6–15.

    Article  PubMed  Google Scholar 

  9. Dacks JB, Field MC. Evolutionary origins and specialisation of membrane transport. Curr Opin Cell Biol. 2018;53:70–6.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  10. Jackson AP, Otto TD, Aslett M, Armstrong SD, Bringaud F, Schlacht A, et al. Kinetoplastid hylogenomics reveals the evolutionary innovations associated with the origins of parasitism. Curr Biol. 2016;26:161–72.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  11. Pipaliya S, Santos R, Salas-Leiva D, Balmer EA, Wirdnam CD, Roger AJ, et al. Unexpected organellar locations of ESCRT machinery in Giardia intestinalis and complex evolutionary dynamics spanning the transition to parasitism in the lineage Fornicata. BMC Biol. 2021;19:1–23.

    Article  CAS  Google Scholar 

  12. Hampl V, Hug L, Leigh JW, Dacks JB, Lang BF, Simpson AGB, et al. Phylogenomic analyses support the monophyly of Excavata and resolve relationships among eukaryotic “supergroups.”. PNAS. 2009;106:3859–64.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  13. Burki F, Roger AJ, Brown MW, Simpson AGB. The new tree of eukaryotes. Trends Ecol Evol. 2020;35:43–55.

    Article  PubMed  Google Scholar 

  14. Hug LA, Baker BJ, Anantharaman K, Brown CT, Probst AJ, Castelle CJ, Butterfield CN, Hernsdorf AW, Amano Y, Ise K, Suzuki Y. A new view of the tree of life. Nature Microbiol. 2016;1(5):1–6.

  15. Caccì OSM, Ryan U. Molecular epidemiology of giardiasis. Mol Biochem Parasitol. 2008;160:75–80.

    Article  CAS  Google Scholar 

  16. Kotloff KL, Nataro JP, Blackwelder WC, Nasrin D, Farag TH, Panchalingam S, et al. Burden and aetiology of diarrhoeal disease in infants and young children in developing countries (the Global Enteric Multicenter Study, GEMS): a prospective, case-control study. Lancet. 2013;382:209–22.

    Article  PubMed  Google Scholar 

  17. Lanata CF, Fischer-Walker CL, Olascoaga AC, Torres CX, Aryee MJ, Black RE, et al. Global causes of diarrheal disease mortality in children <5 years of age: a systematic review. PLoS One. 2013;8:e72788.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  18. Allain T, Amat CB, Motta JP, Manko A, Buret AG. Interactions of Giardia sp. with the intestinal barrier: epithelium, mucus, and microbiota. Tissue Barriers. 2017;5:1–16.

    Article  CAS  Google Scholar 

  19. Fekete E, Allain T, Siddiq A, Sosnowski O, Buret AG. Giardia spp. and the gut microbiota: dangerous liaisons. Front Microbiol. 2021;11:618106.

  20. Faso C, Hehl AB. Membrane trafficking and organelle biogenesis in Giardia lamblia: use it or lose it. Int J Parasitol. 2011;41:471–80.

    Article  CAS  PubMed  Google Scholar 

  21. Benchimol M. The nuclei of Giardia lamblia - new ultrastructural observations. Arch Microbiol. 2005;183:160–8.

    Article  CAS  PubMed  Google Scholar 

  22. Soltys BJ, Falah M, Gupta RS. Identification of endoplasmic reticulum in the primitive eukaryote Giardia lamblia using cryoelectron microscopy and antibody to BiP. J Cell Sci. 1996;109:1909–17.

    Article  CAS  PubMed  Google Scholar 

  23. Tovar J, León-Avila G, Sánchez LB, Sutak R, Tachezy J, Van Der Giezen M, et al. Mitochondrial remnant organelles of Giardia function in iron-sulphur protein maturation. Nature. 2003;426:172–6.

    Article  CAS  PubMed  Google Scholar 

  24. Lanfredi-Rangel A, Attias M, de Carvalho TM, Kattenbach WM, de Souza W. The peripheral vesicles of trophozoites of the primitive protozoan Giardia lamblia may correspond to early and late endosomes and to lysosomes. J Struct Biol. 1998;123:225–35.

    Article  CAS  PubMed  Google Scholar 

  25. Marti M, Regös A, Li Y, Schraner EM, Wild P, Müller N, et al. An ancestaral secretory apparatus in the protozoan parasite Giardia intestinalis. J Biol Chem. 2003;278:24837–48.

    Article  CAS  PubMed  Google Scholar 

  26. Cernikova L, Faso C, Hehl AB. Phosphoinositide-binding proteins mark, shape and functionally modulate highly-diverged endocytic compartments in the parasitic protist Giardia lamblia; 2020.

    Book  Google Scholar 

  27. Zumthor JP, Cernikova L, Rout S, Kaech A, Faso C, Hehl AB. Static clathrin assemblies at the peripheral vacuole plasma membrane interface of the parasitic protozoan Giardia lamblia. PLoS Pathog. 2016;12:1–33.

    Article  CAS  Google Scholar 

  28. Rivero MR, Miras SL, Quiroga R, Rópolo AS, Touz MC. Giardia lamblia low-density lipoprotein receptor-related protein is involved in selective lipoprotein endocytosis and parasite replication. Mol Microbiol. 2011;79:1204–19.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  29. Frontera LS, Moyano S, Quassollo G, Lanfredi-Rangel A, Rópolo AS, Touz MC. Lactoferrin and lactoferricin endocytosis halt Giardia cell growth and prevent infective cyst production. Sci Rep. 2018;8:18020.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  30. Abodeely M, DuBois KN, Hehl A, Stefanic S, Sajid M, DeSouza W, et al. A contiguous compartment functions as endoplasmic reticulum and endosome/lysosome in Giardia lamblia. Eukaryot Cell. 2009;8:1665–76.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  31. Titze B, Genoud C. Volume scanning electron microscopy for imaging biological ultrastructure. Biol Cell. 2016;108:307–23.

    Article  CAS  PubMed  Google Scholar 

  32. Kizilyaprak C, Daraspe J, Humbel BM. Focused ion beam scanning electron microscopy in biology. J Microsc. 2014;254:109–14.

    Article  CAS  PubMed  Google Scholar 

  33. Wei D, Jacobs S, Modla S, Zhang S, Young CL, Cirino R, et al. High-resolution three-dimensional reconstruction of a whole yeast cell using focused-ion beam scanning electron microscopy. BioTechniques. 2012;53:41–8.

    Article  CAS  PubMed  Google Scholar 

  34. Tůmová P, Nohýnková E, Klingl A, Wanner G. A rapid workflow for the characterization of small numbers of unicellular eukaryotes by using correlative light and electron microscopy. J Microbiol Methods. 2020;172:105888.

  35. Dawson SC. An insider’s guide to the microtubule cytoskeleton of Giardia. Cell Microbiol. 2010;12:588–98.

    Article  CAS  PubMed  Google Scholar 

  36. Cardona A, Saalfeld S, Schindelin J, Arganda-Carreras I, Preibisch S, Longair M, Tomancak P, Hartenstein V, Douglas RJ. TrakEM2 software for neural circuit reconstruction. PloS One. 2012;7(6):e38011.

  37. Berg S, Kutra D, Kroeger T, Straehle CN, Kausler BX, Haubold C, et al. ilastik: interactive machine learning for (bio)image analysis. Nat Methods. 2019;16:1226–32.

    Article  CAS  PubMed  Google Scholar 

  38. Sommer C, Straehle C, Koethe U, Hamprecht FA. Ilastik: Interactive learning and segmentation toolkit. In2011 IEEE international symposium on biomedical imaging: From nano to macro 2011 Mar 30 (pp. 230-233). IEEE.

  39. Kan A. Machine learning applications in cell image analysis. Immunol Cell Biol. 2017;95:525–30.

    Article  PubMed  Google Scholar 

  40. Sommer C, Gerlich DW. Machine learning in cell biology – teaching computers to recognize phenotypes. J Cell Sci. 2013;126(126):5529–39.

    CAS  PubMed  Google Scholar 

  41. Poteryaev D, Datta S, Ackema K, Zerial M, Spang A. Identification of the switch in early-to-late endosome transition. Cell. 2010;141:497–508.

    Article  CAS  PubMed  Google Scholar 

  42. Suresh B, Saminathan A, Chakraborty K, Cui C, Krishnan Y. Tubular lysosomes harbor active ion gradients and poise macrophages for phagocytosis. PNAS. 2020;18:2020.12.05.413229.

    Google Scholar 

  43. Hipolito VEB, Ospina-Escobar E, Botelho RJ. Lysosome remodelling and adaptation during phagocyte activation. Cell Microbiol. 2018;20:1–8.

    Article  CAS  Google Scholar 

  44. Combs CA, Shroff H. Fluorescence microscopy: a concise guide to current imaging methods. Curr Protoc Neurosci. 2017;2017:2.1.1–2.1.25.

    Google Scholar 

  45. Klar TA, Jakobs S, Dyba M, Egner A, Hell SW. Fluorescence microscopy with diffraction resolution barrier broken by stimulated emission. Proc Natl Acad Sci U S A. 2000;97:8206–10.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  46. Willig KI, Rizzoli SO, Westphal V, Jahn R, Hell SW. STED microscopy reveals that synaptotagmin remains clustered after synaptic vesicle exocytosis. Nature. 2006;440:935–9.

    Article  CAS  PubMed  Google Scholar 

  47. Kao HP, Kao HP, Verkman a S, Verkman a S. Tracking of single fluorescent particles in three dimensions: use of cylindrical optics to encode particle position. Biophys J. 1994;67:1291–300.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  48. Jones SA, Shim S-H, He J, Zhuang X. Fast, three-dimensional super-resolution imaging of live cells. Nat Methods. 2011;8:499–505.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  49. Huang B, Wang W, Bates M, Zhuang X. Three-dimensional super-resolution reconstruction microscopy. Science (1979). 2008;319:810–3.

    CAS  Google Scholar 

  50. Olivier N, Keller D, Rajan VS, Gönczy P, Manley S. Simple buffers for 3D STORM microscopy. Biomed Opt Express. 2013;4:885.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  51. Dempsey GT, Vaughan JC, Chen KH, Bates M, Zhuang X. Evaluation of fluorophores for optimal performance in localization-based super-resolution imaging. Nat Methods. 2011;8:1027–36.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  52. Ovesný M, Křížek P, Borkovec J, Švindrych Z, Hagen GM. ThunderSTORM: a comprehensive ImageJ plug-in for PALM and STORM data analysis and super-resolution imaging. Bioinformatics. 2014;30:2389–90.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  53. Schindelin J, Arganda-Carreras I, Frise E, Kaynig V, Longair M, Pietzsch T, et al. Fiji: an open-source platform for biological-image analysis. Nat Methods. 2012;9:676–82.

    Article  CAS  PubMed  Google Scholar 

  54. Brown JR, Schwartz CL, Heumann JM, Dawson SC, Hoenger A. A detailed look at the cytoskeletal architecture of the Giardia lamblia ventral disc. J Struct Biol. 2016;194:38–48.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  55. Jørgensen A, Sterud E. Phylogeny of Spironucleus (Eopharyngia: Diplomonadida: Hexamitinae). Protist. 2007;158:247–54.

    Article  PubMed  CAS  Google Scholar 

  56. Kolisko M, Cepicka I, Hampl V, Leigh J, Roger AJ, Kulda J, et al. Molecular phylogeny of diplomonads and enteromonads based on SSU rRNA, alpha-tubulin and HSP90 genes: implications for the evolutionary history of the double karyomastigont of diplomonads. BMC Evol Biol. 2008;8:1–14.

    Article  CAS  Google Scholar 

  57. Xu F, Jerlström-Hultqvist J, Kolisko M, Simpson AGB, Roger AJ, Svärd SG, et al. On the reversibility of parasitism: adaptation to a free-living lifestyle via gene acquisitions in the diplomonad Trepomonas sp. PC1. BMC Biol. 2016;14:1–15.

    PubMed  PubMed Central  CAS  Google Scholar 

  58. Paull GC, Matthews RA. Spironucleus vortens, a possible cause of hole-in-the-head disease in cichlids. Dis Aquat Organ. 2001;45:197–202.

    Article  CAS  PubMed  Google Scholar 

  59. Xu F, Jerlström-Hultqvist J, Einarsson E, Astvaldsson A, Svärd SG, Andersson JO. The genome of Spironucleus salmonicida highlights a fish pathogen adapted to fluctuating environments. PLoS Genet. 2014;10(2):e1004053.

  60. Sterud E, Poynton SL. Spironucleus vortens (Diplomonadida) in the Ide, Leuciscus idus (L.)(Cyprinidae): a warm water hexamitid flagellate found in northern Europe. J Eukaryot Microbiol. 2002;49(2):137–45.

  61. Ástvaldsson Á, Hultenby K, Svärd SG, Jerlström-Hultqvist J. Proximity staining using enzymatic protein tagging in diplomonads. mSphere. 2019;4:1–15.

    Article  Google Scholar 

  62. Day KJ, Casler JC, Glick BS. Budding yeast has a minimal endomembrane system. Dev Cell. 2018;44:56–72.e4.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  63. Lealda N, Silva C, Elias CA. Tritrichomonas foetus: ultrastructure of endocytosis and cytochemistry. Exp Parasitol. 1986;62:405–15.

    Article  Google Scholar 

  64. Rosa IDA, Caruso MB, Rodrigues SP, Geraldo RB, Kist LW, Bogo MR, et al. New insights on the Golgi complex of Tritrichomonas foetus. Parasitology. 2014;141:241–53.

    Article  CAS  Google Scholar 

  65. Schlacht A, Herman EK, Klute MJ, Field MC, Dacks JB. Missing pieces of an ancient puzzle: evolution of the eukaryotic membrane-trafficking system. Cold Spring Harbor Perspect Biol. 2014;6(10):a016048.

  66. Midlej V, Pereira-Neves A, Kist LW, Bogo MR, Benchimol M. Ultrastructural features of Tritrichomonas mobilensis and comparison with Tritrichomonas foetus. Vet Parasitol. 2011;182:171–80.

    Article  PubMed  Google Scholar 

  67. Traub LM. Regarding the amazing choreography of clathrin coats. PLoS Biol. 2011;9:3–7.

    Article  CAS  Google Scholar 

  68. Eddy SR. Accelerated profile HMM searches. PLoS Comput Biol. 2011;7(10):e1002195.

  69. Kirchhausen T, Owen D, Harrison SC. Molecular structure, function, and dynamics of clathrin-mediated membrane traffic. Cold Spring Harbor Perspect Biol. 2014;6(5):a016725.

  70. Morgan GW, Allen CL, Jeffries TR, Hollinshead M, Field MC. Developmental and morphological regulation of clathrin-mediated endocytosis in Trypanosoma brucei. J Cell Sci. 2001;114:2605–15.

    Article  CAS  PubMed  Google Scholar 

  71. Adung’a VO, Gadelha C, Field MC. Proteomic analysis of clathrin interactions in trypanosomes reveals dynamic evolution of endocytosis. Traffic. 2013;14:440–57.

    Article  PubMed  CAS  Google Scholar 

  72. Kaksonen M, Toret CP, Drubin DG. A modular design for the clathrin- and actin-mediated endocytosis machinery. Cell. 2005;123:305–20.

    Article  CAS  PubMed  Google Scholar 

  73. Karnkowska A, Vacek V, Zubáčová Z, Treitli SC, Petrželková R, Eme L, Novák L, Žárský V, Barlow LD, Herman EK, Soukal P. A eukaryote without a mitochondrial organelle. Curr Biol. 2016;26(10):1274–84.

  74. Karnkowska A, Treitli SC, Brzoň O, Novák L, Vacek V, Soukal P, et al. The oxymonad genome displays canonical eukaryotic complexity in the absence of a mitochondrion. Mol Biol Evol. 2019.

  75. Füssy Z, Vinopalová M, Treitli SC, Pánek T, Smejkalová P, Čepička I, Doležal P, Hampl V. Retortamonads from vertebrate hosts share features of anaerobic metabolism and pre-adaptations to parasitism with diplomonads. Parasitol Int. 2021;82:102308.

  76. Cheon S, Zhang J, Park C. Is phylotranscriptomics as reliable as phylogenomics? Mol Biol Evol. 2020;37:3672–83.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  77. Woo YH, Ansari H, Otto TD, Linger CMK, Olisko MK, Michálek J, et al. Chromerid genomes reveal the evolutionary path from photosynthetic algae to obligate intracellular parasites. Elife. 2015;4:1–41.

    Article  Google Scholar 

  78. Richardson E, Dacks JB. Distribution of membrane trafficking system components across ciliate diversity highlights heterogenous organelle-associated machinery. Traffic. 2022.

  79. Fotin A, Cheng Y, Grigorieff N, Walz T, Harrison SC, Kirchhausen T. Structure of an auxilin-bound clathrin coat and its implications for the mechanism of uncoating. Nature. 2004;429:2004.

    Google Scholar 

  80. Rapoport I, Boll W, Yu A, Bocking T, Kirchhausen T. A motif in the clathrin heavy chain required for the Hsc70/auxilin uncoating reaction. Mol Biol Cell. 2008;19:3250–63.

    Article  Google Scholar 

  81. Suga H, Chen Z, De Mendoza A, Sebé-Pedrós A, Brown MW, Kramer E, et al. The Capsaspora genome reveals a complex unicellular prehistory of animals. Nat Commun. 2013;4:1–9.

    Article  CAS  Google Scholar 

  82. King N, Westbrook MJ, Young SL, Kuo A, Abedin M, Chapman J, et al. The genome of the choanoflagellate Monosiga brevicollis and the origin of metazoans. Nature. 2008;451:783–8.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  83. Fairclough SR, Chen Z, Kramer E, Zeng Q, Young S, Robertson HM, et al. Premetazoan genome evolution and the regulation of cell differentiation in the choanoflagellate Salpingoeca rosetta. Genome Biol. 2013;14:1–15.

    Article  CAS  Google Scholar 

  84. Tanifuji G, Takabayashi S, Kume K, Takagi M, Nakayama T, Kamikawa R, et al. The draft genome of Kipferlia bialata reveals reductive genome evolution in fornicate parasites. PLoS One. 2018.

  85. Leger MM, Kolisko M, Kamikawa R, Stairs CW, Kume K, Čepička I, Silberman JD, Andersson JO, Xu F, Yabuki A, Eme L. Organelles that illuminate the origins of Trichomonas hydrogenosomes and Giardia mitosomes. Nat Ecol Evol. 2017;1(4):1–7.

  86. Salas-Leiva DE, Tromer EC, Curtis BA, Jerlström-Hultqvist J, Kolisko M, Yi Z, et al. Genomic analysis finds no evidence of canonical eukaryotic DNA processing complexes in a free-living protist. Nat Commun. 2021;12:1–13.

    CAS  Google Scholar 

  87. Manna PT, Obado SO, Boehm C, Gadelha C, Sali A, Chait BT, et al. Lineage-specific proteins essential for endocytosis in trypanosomes. J Cell Sci. 2017;130:1379–92.

    CAS  PubMed  PubMed Central  Google Scholar 

  88. Aurrecoechea C, Brestelli J, Brunk BP, Carlton JM, Dommer J, Fischer S, et al. GiardiaDB and TrichDB: integrated genomic resources for the eukaryotic protist pathogens Giardia lamblia and Trichomonas vaginalis. Nucleic Acids Res. 2009;37(SUPPL. 1):526–30.

    Article  CAS  Google Scholar 

  89. Carlton JM, Hirt RP, Silva JC, Delcher AL, Schatz M, Zhao Q, et al. Draft genome sequence of the sexually transmitted pathogen Trichomonas vaginalis. Science (1979). 2007;315:207–13.

    Google Scholar 

  90. Lax G, Eglit Y, Eme L, Bertrand EM, Roger AJ, Simpson AGB. Hemimastigophora is a novel supra-kingdom-level lineage of eukaryotes. Nature. 2018;564:410–4.

    Article  CAS  PubMed  Google Scholar 

  91. Zimmermann L, Stephens A, Nam SZ, Rau D, Kübler J, Lozajic M, et al. A completely reimplemented MPI bioinformatics toolkit with a new HHpred server at its core. J Mol Biol. 2017;430:1–7.

    Google Scholar 

  92. Tunyasuvunakool K, Adler J, Wu Z, Green T, Zielinski M, Žídek A, et al. Highly accurate protein structure prediction for the human proteome. Nature. 2021;596:590–6.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  93. Jumper J, Evans R, Pritzel A, Green T, Figurnov M, Ronneberger O, et al. Highly accurate protein structure prediction with AlphaFold. Nature. 2021;596:583–9.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  94. Zhang Y, Skolnick J. TM-align: a protein structure alignment algorithm based on the TM-score. Nucleic Acids Res. 2005;33:2302–9.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  95. Kufareva I, Abagyan R. Methods of protein structure comparison. Methods Mol Biol. 2012;857:231–57.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  96. Wilbur JD, Hwang PK, Ybe JA, Lane M, Sellers BD, Jacobson MP, et al. Conformation switching of clathrin light chain regulates clathrin lattice assembly. Dev Cell. 2010;18:841–8.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  97. Jerlström-Hultqvist J, Einarsson E, Svärd SG. Stable transfection of the diplomonad parasite Spironucleus salmonicida. Eukaryot Cell. 2012;11:1353–61.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  98. Barash NR, Nosala C, Pham JK, Mcinally SG, Gourguechon S, Dawson SC. Giardia colonizes and encysts in high-density foci in the murine small intestine. mSphere. 2017;2:1–20.

    Article  Google Scholar 

  99. Adam RD. Biology of Giardia lamblia. Clin Microbiol Rev. 2001;14:447–69.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  100. Carranza PG, Lujan HD. New insights regarding the biology of Giardia lamblia. Microbes Infect. 2009;12:71–80.

    Article  PubMed  CAS  Google Scholar 

  101. Cotton JA, Beatty JK, Buret AG. Host parasite interactions and pathophysiology in Giardia infections. Int J Parasitol. 2011;41:925–33.

    Article  CAS  PubMed  Google Scholar 

  102. Touz M, Feliziani C, Rópolo A. Membrane-associated proteins in Giardia lamblia. Genes (Basel). 2018;9:404.

    Article  CAS  Google Scholar 

  103. Jacquemet G, Carisey AF, Hamidi H, Henriques R, Leterrier C. The cell biologist's guide to super-resolution microscopy. J Cell Sci. 2020;133(11):jcs240713.

  104. Hamann E, Tegetmeyer HE, Di R, Littmann S, Ahmerkamp S, Chen J, et al. Syntrophic linkage between predatory Carpediemonas and specific prokaryotic populations. ISME J. 2017;11:1205–17.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  105. Yubuki N, Simpson AGB, Leander BS. Comprehensive ultrastructure of Kipferlia bialata provides evidence for character evolution within the Fornicata (Excavata). Protist. 2013;164:423–39.

    Article  CAS  PubMed  Google Scholar 

  106. Yubuki N, Huang SSCC, Leander BS. Comparative ultrastructure of fornicate excavates, including a novel free-living relative of diplomonads: Aduncisulcus paluster gen. et sp. nov. Protist. 2016;167:584–96.

    Article  PubMed  Google Scholar 

  107. Ebneter JA, Hehl AB. The single epsin homolog in Giardia lamblia localizes to the ventral disk of trophozoites and is not associated with clathrin membrane coats. Mol Biochem Parasitol. 2014;197:24–7.

    Article  CAS  PubMed  Google Scholar 

  108. Benchimol M. A new set of vesicles in Giardia lamblia. Exp Parasitol. 2002;102:30–7.

    Article  PubMed  Google Scholar 

  109. McCaffery JM, Gillin FD. Giardia lamblia: ultrastructural basis of protein transport during growth and encystation. Exp Parasitol. 1994;79:220–35.

    Article  CAS  PubMed  Google Scholar 

  110. McCaffery JM, Faubert GM, Gillin FD. Traffic of a trophozoite variant surface protein and a major cyst wall epitope during growth, encystation, and antigenic switching. Exp Parasitol. 1994;79:236–49.

    Article  CAS  PubMed  Google Scholar 

  111. Stefanic S, Palm D, Svärd SG, Hehl AB. Organelle proteomics reveals cargo maturation mechanisms associated with Golgi-like encystation vesicles in the early-diverged protozoan Giardia lamblia *. J Biol Chem. 2006;281:7595–604.

    Article  CAS  PubMed  Google Scholar 

  112. Stefanic S, Morf L, Kulangara C, Regös A, Sonda S, Schraner E, et al. Neogenesis and maturation of transient Golgi-like cisternae in a simple eukaryote. J Cell Sci. 2009;122:2846–56.

    Article  CAS  PubMed  Google Scholar 

  113. Acosta-Virgen K, Chávez-Munguía B, Talamás-Lara D, Lagunes-Guillén A, Martínez-Higuera A, Lazcano A, et al. Giardia lamblia: identification of peroxisomal-like proteins. Exp Parasitol. 2018;191:36–43.

    Article  CAS  PubMed  Google Scholar 

  114. Link F, Borges AR, Jones NG, Engstler M. To the surface and back: exo- and endocytic pathways in Trypanosoma brucei. Front Cell Dev Biol. 2021;9:1–15.

    Article  Google Scholar 

  115. Allen CL, Goulding D, Field MC. Clathrin-mediated endocytosis is essential in Trypanosoma brucei. EMBO J. 2003;22:4991–5002.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  116. Morf L, Spycher C, Rehrauer H, Fournier CA, Morrison HG, Hehl AB. The transcriptional response to encystation stimuli in Giardia lamblia is restricted to a small set of genes. Eukaryot Cell. 2010;9:1566–76.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  117. Gaechter V, Schraner E, Wild P, Hehl AB. The single dynamin family protein in the primitive protozoan giardia lamblia is essential for stage conversion and endocytic transport. Traffic. 2008;9:57–71.

    Article  CAS  PubMed  Google Scholar 

  118. Konrad C, Spycher C, Hehl AB. Selective condensation drives partitioning and sequential secretion of cyst wall proteins in differentiating Giardia lamblia. PLoS Pathog. 2010;6:e1000835.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  119. Arganda-Carreras I, Kaynig V, Rueden C, Eliceiri KW, Schindelin J, Cardona A, et al. Trainable Weka segmentation: a machine learning tool for microscopy pixel classification. Bioinformatics. 2017;33:2424–6.

    Article  CAS  PubMed  Google Scholar 

  120. Mateos JM, Guhl B, Doehner J, Barmettler G, Kaech A, Ziegler U. Topographic contrast of ultrathin cryo-sections for correlative super-resolution light and electron microscopy. Sci Rep. 2016;6:34062.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  121. Edgar RC. MUSCLE: Multiple sequence alignment with high accuracy and high throughput. Nucleic Acids Res. 2004.

  122. Lu S, Wang J, Chitsaz F, Derbyshire MK, Geer RC, Gonzales NR, et al. CDD/SPARCLE: the conserved domain database in 2020. Nucleic Acids Res. 2020;48:D265–8.

    Article  CAS  PubMed  Google Scholar 

  123. Mistry J, Chuguransky S, Williams L, Qureshi M, Salazar GA, Sonnhammer ELL, et al. Pfam: the protein families database in 2021. Nucleic Acids Res. 2021;49:D412–9.

    Article  CAS  PubMed  Google Scholar 

  124. Mitchell AL, Attwood TK, Babbitt PC, Blum M, Bork P, Bridge A, et al. InterPro in 2019: improving coverage, classification and access to protein sequence annotations. Nucleic Acids Res. 2019;47:D351–60.

    Article  CAS  PubMed  Google Scholar 

  125. Letunic I, Bork P. 20 years of the SMART protein domain annotation resource. Nucleic Acids Res. 2017;46:493–6.

    Article  CAS  Google Scholar 

  126. Perez-Riverol Y, Csordas A, Bai J, Bernal-Llinares M, Hewapathirana S, Kundu DJ, et al. The PRIDE database and related tools and resources in 2019: improving support for quantification data. Nucleic Acids Res. 2019;47:D442–50.

    Article  CAS  PubMed  Google Scholar 

Download references


Imaging and image analysis were performed with equipment from the Centre of Microscopy and Image Analysis (ZMB) of the University of Zurich. We thank the following members of the ZMB for the technical and scientific support: Dr. Jana Döhner, Dr. Moritz Kirchmann, Dr. Dominik Hänni, Dr. José Mateos and Dr. Urs Ziegler.


ABH and CF are funded by the Swiss National Foundation grants 31003A-166437 and PR00P3_179813, respectively.

Author information

Authors and Affiliations



RS, ABH and CF designed and curated the study. RS performed all the experiments and analysed all the ensuing experimental data with the exception of S. salmonicida culturing, transfection and uptake experiments were performed by AA and SS, and SEM experiments were performed by JPZ. SVP, RS and JBD performed the molecular phylogeny analyses. RS, ABH, JBD and CF wrote and revised the manuscript. All authors read and approved the final manuscript prior to submission.

Authors’ information

Twitter handles: @Filipe_RFRS (Rui Santos); @AsgeirAstvalds (Ásgeir Ástvaldsson); @pipa_sw (Shweta V. Pipaliya); @DacksLab1 (Joel B. Dacks); @LabFaso (Carmen Faso).

Corresponding author

Correspondence to Carmen Faso.

Ethics declarations

Ethics approval and consent to participate

No animal experiments were performed nor any patient-derived material or data was used for this work.

Consent for publication

All authors consent to the publication of this manuscript and its contents.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Additional file 1: Fig. S1.

Rendering of a G. lamblia trophozoite scanned with FIB-SEM reveals the cell’s inner ultrastructure. (A) 3D view of acquired FIB-SEM trophozoite data. (B) Single slice showing inner cellular structures such as cytoskeleton elements at the median body (MB), the ventral disk (VD), the endoplasmic reticulum (ER), mitosomes (m) and peripheral vacuoles (PV), highlighted in the region of interest (ROI). (C) Segmentation of different categories of the dataset: cell volume (138 μm3), cytoskeleton, endoplasmic reticulum, peripheral vacuoles, small vesicles and mitosomes. (D) Mitosome volume (N = 14, violin-plot) was determined post segmentation at an average volume of 0.001093±0.0005698μm3 in a 95% confidence interval between [0.0007643, 0.001422] μm3.

Additional file 2: Video S1. Three-dimensional rendering of endocytic compartments in G. lamblia derived from FIB-SEM sectioning and imaging. Scale bar 1 μm.

Additional file 3: Fig. S2.

Cryo-SEM of freeze-fractured trophozoites reveals varying vacuolar morphology in Giardia lamblia. (A) Overview of cryo-preserved Giardia trophozoites subjected to freeze-fracture and SEM imaging. Nuclei (N), Endoplasmic Reticulum (ER), Ventral Disk (VD) and peripheral endocytic compartments (PEC) and plasma membrane (PM) are clearly identifiable. (B and C) Insets showing different PV/PEC morphology: vesicular (asterisk) and tubular (hashtag). Scale bar: (A) 2 μm and (B and C) 500 nm.

Additional file 4: Fig. S3.

TEM investigation of Giardia lamblia endocytic and secretory pathway. (A) Overview of a trophozoite. Different PV/PEC structures, vesicular and tubular are observed, together with small vesicles (SV). The N (nucleus) and ER are also highlighted. (B) Close up on tubular PV/PECs (hashtag). (C) Close up on vesicular PV/PECs (asterisk) and SVs (arrowhead). Scale bars: (A) 2 μm, (B) 1 μm and (C) 500 nm.

Additional file 5: Video S2. Comparison between confocal and STED imaging of Giardia PV/PECs. Scale bar: 3 μm.

Additional file 6: Video S3. Tri-dimensional reconstruction of PV/PECs from STORM data.

Additional file 7: Table S1.

PV/PECs volume comparison as calculated in FIB-SEM and STORM experiments.

Additional file 8: Video S4. Tri-dimensional confocal imaging of S. vortens with Dextran-Texas Red. Both peripheral and near-nuclear endosome-like vacuoles are observed.

Additional file 9: Video S5. Tri-dimensional STED imaging of S. vortens with Dextran Alexa Fluor 594 reveals endosome-like vacuoles in greater detail.

Additional file 10: Fig. S4.

TEM investigation of S. salmonicida endocytic and secretory pathway. (A) S. salmonicida presents vacuolar formations close to the plasma membrane. Cells also present a prominent endoplasmic reticulum (ER; blue-framed inset). (B) Highlight of vacuolar formations (V) and ER. (C) Second cell displaying an abundance of PV close to its plasma membrane. (D) Highlight of vacuoles (V) and the prominent ER that connects to the plasma membrane (asterisk). (E) S. salmonicida PVs average a diameter of 205±62.6 nm (N=114) in a 95% confidence interval of [193;217]. (F) S. vortens peripheral vacuoles are larger than S. salmonicida vacuoles in a statistically significant manner (p-value < 0.0001). Diameters were manually determined. Scale bars: (A and C) 2 μm and (B and D) 500 nm.

Additional file 11: Video S6. Tri-dimensional STED imaging of T. foetus with Dextran Alexa Fluor 594 reveals endosome-like vacuoles in greater detail.

Additional file 12: Fig. S5.

TEM investigation of T. foetus Golgi vesicles. (A) More than one Golgi apparatus (G) can be found per cell. These organelles resemble canonical stacked Golgi releasing small coated vesicles. (B) These vesicles average a diameter of 58.4±13.1 nm (N=128) in a 95% confidence interval of [56.1;60.7] nm. Scale bar: (A) 500 nm.

Additional file 13: Table S2.

Queries used for CHC HHM profile building. Table S3. Results from Pan-Eukaryotic search of CHC homologues in available Proteomes. Table S5. Queries used for CLC HHM profile building. Table S6. Results for GlCLC search in available Giardia genomes. Table S7. Results for bona fide CLC present in other genomes/transcriptomes. Table S8. Comparison of GlCLC with bona fide CLC. Table S9. Ten best hits from HHPred.

Additional file 14: Fig. S6.

Pan-Eukaryotic prediction of clathrin heavy chain protein domains. Pfam analysis of predicted protein domains for several clathrin heavy chain proteins sequences from the following species: Giardia lamblia, Spironucleus vortens, Spironucleus salmonicida, Trepomonas sp., Hexamita inflata, Dysnectes brevis, Kipferlia bialata, Carpediomonas membranifera, Aduncisulcus paluster, Chilomastix cuspidata, Trypanosoma brucei, Naegleria gruberi, Tritrichomonas foetus, Monocercomonoides exilis, Tetrahymena thermophila, Hemimastix kukwestjiik, Chlamydomonas reinhardtii, Dyctiostilium discoideum, Saccharomyces cerevisiae, Caenorhabditis elegans, Homo sapiens, Salpingoeca rosetta, Capsospora owczarzaki and Monosiga brevicollis. A general decrease in domain complexity is observed in excavates compared with higher eukaryotes. CLOs: Carpediomonas-like organisms. Diplom: Diplomonada.

Additional file 15: Table S10.

Search for domains from Clathrin heavy chain super family repeats.

Additional file 16: Fig. S7.

The QLMLT motif is exclusive to Holozoa. Alignment of the C-terminii of CHC sequences from selected Opisthokonta, Archaeplastida, Amoebozoa and SAR species highlights the present of the QLMLT uncoating motif only in Holozoa supergroup. The positioning of the QLMLT is highlighted in blue.

Additional file 17: Fig. S8.

Calculation of Giardia ACLC synonymous vs non- synonymous mutation ratio (ω = ks/kn). (A) Phylogenetic tree resulting of maximum likelihood analysis of the Giardia ACLC sequences. Each node is represented by a number. (B) Overall ω < 1 indicating there is no selective pressure on Giardia ACLC.

Additional file 18: Fig. S9.

SsCHC is distributed in the cytosol and interacts with a putative light chain structural analogue. (A) SsCHC was tagged C-terminally with three HA tags and immune-localized to the cellcytosol. Signal was observed in 88% of the analyzed cells (N = 171). (B) High resolution imaging of SsCHC using confocal imaging reveals CHC foci. (C) Single native co-IP analysis of HA-tagged SsCHC reporter (predicted at ca. 210 kDa) including distribution of the 171 proteins found in higher abundance with respect to a control co-IP experiment using extracts of non-transgenic Ss cells. (D) Qualitative immunoblot analysis of samples from the single native co-IP of HA-tagged SsCHC reporter. I: soluble native co-IP input; P: insoluble cell debris post cell lysis; F: native co-IP flow-through; B: anti-HA beads. WT: non-transgenic Ss cells. MW: molecular weight. (E) Ab initio in silico protein modelling with AlphaFold of Ss11905, GlACLC, TbCLC and HsCLC. TM-align and RMSD scores for predicted structures of Giardia ACLC, Trypanosoma brucei CLC and Ss11905 with respect to Homo sapiens CLC show overall structural conservation with respect to a bona fide CLC. Scale bars: (A) 20 μm. (B) 5 μm.

Additional file 19: Video S7. Tri-dimensional high resolution confocal imaging and representation of SsCHC-3xHA foci in the cell cytoplasm.

Additional file 20: Table S11.

SsCHC co-IP results.

Additional file 21: Fig. S10.

Trepomonas sp. PC1 also harbours a putative CLC analogue. Ab initio protein modelling of TPC1_16039, orthologous to Ss11905 in combination with Ss11905, GlACLC, TbCLC and HsCLC. RMSD and TM-align cores show overall structural conservation with respect to a bona fide CLC.

Additional file 22: Table S12.

Selected CLC-related sequences used for predictive in silico modelling analyses.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Santos, R., Ástvaldsson, Á., Pipaliya, S.V. et al. Combined nanometric and phylogenetic analysis of unique endocytic compartments in Giardia lamblia sheds light on the evolution of endocytosis in Metamonada. BMC Biol 20, 206 (2022).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: