The genome of Spironucleus salmonicida highlights a fish pathogen adapted to fluctuating environments

Genomics of Preaxostyla Flagellates Illuminates the Path Towards the Loss of Mitochondria

The notion that mitochondria cannot be lost was shattered with the report of an oxymonad Monocercomonoides exilis, the first eukaryote arguably without any mitochondrion. Yet, questions remain about whether this extends beyond the single species and how this transition took place. The Oxymonadida is a group of gut endobionts taxonomically housed in the Preaxostyla which also contains free-living flagellates of the genera Trimastix and Paratrimastix. The latter two taxa harbour conspicuous mitochondrion-related organelles (MROs). Here we report high-quality genome and transcriptome assemblies of two Preaxostyla representatives, the free-living Paratrimastix pyriformis and the oxymonad Blattamonas nauphoetae. We performed thorough comparisons among all available genomic and transcriptomic data of Preaxostyla to further decipher the evolutionary changes towards amitochondriality, endobiosis, and unstacked Golgi. Our results provide insights into the metabolic and endomembrane evolution, but most strikingly the data confirm the complete loss of mitochondria for all three oxymonad species investigated (M. exilis, B. nauphoetae, and Streblomastix strix), suggesting the amitochondriate status is common to a large part if not the whole group of Oxymonadida. This observation moves this unique loss to 100 MYA when oxymonad lineage diversified.

An ancient metalloenzyme evolves through metal preference modulation

Evolution creates functional diversity of proteins, the essential building blocks of all biological systems. However, studies of natural proteins sampled across the tree of life and evaluated in a single experimental system are lacking. Almost half of enzymes require metals, and metalloproteins tend to optimally utilize the physicochemical properties of a specific metal co-factor. Life must adapt to changes in metal bioavailability, including those during the transition from anoxic to oxic Earth or pathogens' exposure to nutritional immunity. These changes can challenge the ability of metalloenzymes to maintain activity, presumptively driving their evolution. Here we studied metal-preference evolution within the natural diversity of the iron/manganese superoxide dismutase (SodFM) family of reactive oxygen species scavengers. We identified and experimentally verified residues with conserved roles in determining metal preference that, when combined with an understanding of the protein's evolutionary history, improved prediction of metal utilization across the five SodFM subfamilies defined herein. By combining phylogenetics, biochemistry and structural biology, we demonstrate that SodFM metal utilization can be evolutionarily fine tuned by sliding along a scale between perfect manganese and iron specificities. Over the history of life, SodFM metal preference has been modulated multiple independent times within different evolutionary and ecological contexts, and can be changed within short evolutionary timeframes.

A chromosome-scale reference genome for Spironucleus salmonicida

Spironucleus salmonicida is a diplomonad causing systemic infection in salmon. The first S. salmonicida genome assembly was published 2014 and has been a valuable reference genome in protist research. However, the genome assembly is fragmented without assignment of the sequences to chromosomes. In our previous Giardia genome study, we have shown how a fragmented genome assembly can be improved with long-read sequencing technology complemented with optical maps. Combining Pacbio long-read sequencing technology and optical maps, we are presenting here this new S. salmonicida genome assembly in nine near-complete chromosomes with only three internal gaps at long repeats. This new genome assembly is not only more complete sequence-wise but also more complete at annotation level, providing more details into gene families, gene organizations and chromosomal structure. This near-complete reference genome will aid comparative genomics at chromosomal level, and serve as a valuable resource for the diplomonad community and protist research.

Measurement(s)	genomic_DNA • sequence_assembly • sequence feature annotation
Technology Type(s)	SMRT Sequencing • sequence assembly process • sequence annotation
Sample Characteristic - Organism	Spironucleus salmonicida

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

Background

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.

Results

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.

Conclusions

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.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12915-022-01402-3.

Precise gene models using long-read sequencing reveal a unique poly(A) signal in Giardia lamblia

During pre-mRNA processing, the poly(A) signal is recognized by a protein complex that ensures precise cleavage and polyadenylation of the nascent transcript. The location of this cleavage event establishes the length and sequence of the 3' UTR of an mRNA, thus determining much of its post-transcriptional fate. Using long-read sequencing, we characterize the polyadenylation signal and related sequences surrounding Giardia lamblia cleavage sites for over 2600 genes. We find that G. lamblia uses an AGURAA poly(A) signal, which differs from the mammalian AAUAAA. We also describe how G. lamblia lacks common auxiliary elements found in other eukaryotes, along with the proteins that recognize them. Further, we identify 133 genes with evidence of alternative polyadenylation. These results suggest that despite pared-down cleavage and polyadenylation machinery, 3' end formation still appears to be an important regulatory step for gene expression in G. lamblia.

Oxygen levels are key to understanding "Anaerobic" protozoan pathogens with micro-aerophilic lifestyles

Publications abound on the physiology, biochemistry and molecular biology of "anaerobic" protozoal parasites as usually grown under "anaerobic" culture conditions. The media routinely used are poised at low redox potentials using techniques that remove O₂ to "undetectable" levels in sealed containers. However there is growing understanding that these culture conditions do not faithfully resemble the O₂ environments these organisms inhabit. Here we review for protists lacking oxidative energy metabolism, the oxygen cascade from atmospheric to intracellular concentrations and relevant methods of measurements of O₂, some well-studied parasitic or symbiotic protozoan lifestyles, their homeodynamic metabolic and redox balances, organism-drug-oxygen interactions, and the present and future prospects for improved drugs and treatment regimes.

Diversity of the Tryptophanase Gene and Its Evolutionary Implications in Living Organisms

Tryptophanase encoded by the gene tnaA is a pyridoxal phosphate-dependent enzyme that catalyses the conversion of tryptophan to indole, which is commonly used as an intra- and interspecies signalling molecule, particularly by microbes. However, the production of indole is rare in eukaryotic organisms. A nucleotide and protein database search revealed tnaA is commonly reported in various Gram-negative bacteria, but that only a few Gram-positive bacteria and archaea possess the gene. The presence of tnaA in eukaryotes, particularly protozoans and marine organisms, demonstrates the importance of this gene in the animal kingdom. Here, we document the distribution of tnaA and its acquisition and expansion among different taxonomic groups, many of which are usually categorized as non-indole producers. This study provides an opportunity to understand the intriguing role played by tnaA, and its distribution among various types of organisms.

Genomic analysis finds no evidence of canonical eukaryotic DNA processing complexes in a free-living protist

Cells replicate and segregate their DNA with precision. Previous studies showed that these regulated cell-cycle processes were present in the last eukaryotic common ancestor and that their core molecular parts are conserved across eukaryotes. However, some metamonad parasites have secondarily lost components of the DNA processing and segregation apparatuses. To clarify the evolutionary history of these systems in these unusual eukaryotes, we generated a genome assembly for the free-living metamonad Carpediemonas membranifera and carried out a comparative genomics analysis. Here, we show that parasitic and free-living metamonads harbor an incomplete set of proteins for processing and segregating DNA. Unexpectedly, Carpediemonas species are further streamlined, lacking the origin recognition complex, Cdc6 and most structural kinetochore subunits. Carpediemonas species are thus the first known eukaryotes that appear to lack this suite of conserved complexes, suggesting that they likely rely on yet-to-be-discovered or alternative mechanisms to carry out these fundamental processes.

Eukaryote-Conserved Methylarginine Is Absent in Diplomonads and Functionally Compensated in Giardia

Methylation is a common posttranslational modification of arginine and lysine in eukaryotic proteins. Methylproteomes are best characterized for higher eukaryotes, where they are functionally expanded and evolved complex regulation. However, this is not the case for protist species evolved from the earliest eukaryotic lineages. Here, we integrated bioinformatic, proteomic, and drug-screening data sets to comprehensively explore the methylproteome of Giardia duodenalis-a deeply branching parasitic protist. We demonstrate that Giardia and related diplomonads lack arginine-methyltransferases and have remodeled conserved RGG/RG motifs targeted by these enzymes. We also provide experimental evidence for methylarginine absence in proteomes of Giardia but readily detect methyllysine. We bioinformatically infer 11 lysine-methyltransferases in Giardia, including highly diverged Su(var)3-9, Enhancer-of-zeste and Trithorax proteins with reduced domain architectures, and novel annotations demonstrating conserved methyllysine regulation of eukaryotic elongation factor 1 alpha. Using mass spectrometry, we identify more than 200 methyllysine sites in Giardia, including in species-specific gene families involved in cytoskeletal regulation, enriched in coiled-coil features. Finally, we use known methylation inhibitors to show that methylation plays key roles in replication and cyst formation in this parasite. This study highlights reduced methylation enzymes, sites, and functions early in eukaryote evolution, including absent methylarginine networks in the Diplomonadida. These results challenge the view that arginine methylation is eukaryote conserved and demonstrate that functional compensation of methylarginine was possible preceding expansion and diversification of these key networks in higher eukaryotes.

The reduced ARF regulatory system in Giardia intestinalis pre-dates the transition to parasitism in the lineage Fornicata

Giardia intestinalis is an enteric pathogen with an extremely modified membrane trafficking system, lacking canonical compartments such as the Golgi, endosomes, and intermediate vesicle carriers. By comparison the fornicate relatives of Giardia possess greater endomembrane system complexity. In eukaryotes, the ADP ribosylation factor (ARF) GTPase regulatory system proteins, which consist of the small GTPase ARF1, and its guanine exchange nucleotide factors (GEFs) and GTPase activating proteins (GAPs), coordinate temporal and directional trafficking of cargo vesicles by recognizing and interacting with heterotetrameric coat complexes at pre-Golgi and post-Golgi interfaces. To understand the evolution of this regulatory system across the fornicate lineage, we have performed comparative genomic and phylogenetic analyses of the ARF GTPases, and their regulatory GAPs and GEFs in fornicate genomes and transcriptomes. Prior to our analysis of the fornicates, we first establish that the ARF GAP sub-family ArfGAP with dual PH domains (ADAP) is sparsely distributed but present in at least four eukaryotic supergroups and thus was likely present in the Last Eukaryotic Common Ancestor (LECA). Next, our collective comparative genomic and phylogenetic investigations into the ARF regulatory proteins in fornicates identify a duplication of ARF1 GTPase yielding two paralogues of ARF1F proteins, ancestral to all fornicates and present in all examined isolates of Giardia. However, the ARF GEF and ARF GAP complement is reduced compared with the LECA. This investigation shows that the system was significantly streamlined prior to the fornicate ancestor but was not further reduced concurrent with a transition into a parasitic lifestyle.

Retortamonads from vertebrate hosts share features of anaerobic metabolism and pre-adaptations to parasitism with diplomonads

Although the mitochondria of extant eukaryotes share a single origin, functionally these organelles diversified to a great extent, reflecting lifestyles of the organisms that host them. In anaerobic protists of the group Metamonada, mitochondria are present in reduced forms (also termed hydrogenosomes or mitosomes) and a complete loss of mitochondrion in Monocercomonoides exilis (Metamonada:Preaxostyla) has also been reported. Within metamonads, retortamonads from the gastrointestinal tract of vertebrates form a sister group to parasitic diplomonads (e.g. Giardia and Spironucleus) and have also been hypothesized to completely lack mitochondria. We obtained transcriptomic data from Retortamonas dobelli and R. caviae and searched for enzymes of the core metabolism as well as mitochondrion- and parasitism-related proteins. Our results indicate that retortamonads have a streamlined metabolism lacking pathways for metabolites they are probably capable of obtaining from prey bacteria or their environment, reminiscent of the biochemical arrangement in other metamonads. Retortamonads were surprisingly found do encode homologs of components of Giardia's remarkable ventral disk, as well as homologs of regulatory NEK kinases and secreted lytic enzymes known for involvement in host colonization by Giardia. These can be considered pre-adaptations of these intestinal microorganisms to parasitism. Furthermore, we found traces of the mitochondrial metabolism represented by iron‑sulfur cluster assembly subunits, subunits of mitochondrial translocation and chaperone machinery and, importantly, [FeFe]‑hydrogenases and hydrogenase maturases (HydE, HydF and HydG). Altogether, our results strongly suggest that a remnant mitochondrion is still present.

Metabolic Reconstruction Elucidates the Lifestyle of the Last Diplomonadida Common Ancestor

Diplomonads are a group of microbial eukaryotes found in oxygen-poor environments. There are both parasitic (e.g., Giardia intestinalis) and free-living (e.g., Trepomonas) members in the group.

The identification of ancestral traits is essential to understanding the evolution of any group. In the case of parasitic groups, this helps us understand the adaptation to this lifestyle and a particular host. Most diplomonads are parasites, but there are free-living members of the group nested among the host-associated diplomonads. Furthermore, most of the close relatives within Fornicata are free-living organisms. This leaves the lifestyle of the ancestor unclear. Here, we present metabolic maps of four different diplomonad species. We identified 853 metabolic reactions and 147 pathways present in at least one of the analyzed diplomonads. Our study suggests that diplomonads represent a metabolically diverse group in which differences correlate with different environments (e.g., the detoxification of arsenic). Using a parsimonious analysis, we also provide a description of the putative metabolism of the last Diplomonadida common ancestor. Our results show that the acquisition and loss of reactions have shaped metabolism since this common ancestor. There is a net loss of reaction in all branches leading to parasitic diplomonads, suggesting an ongoing reduction in the metabolic capacity. Important traits present in host-associated diplomonads (e.g., virulence factors and the synthesis of UDP-N-acetyl-d-galactosamine) are shared with free-living relatives. The last Diplomonadida common ancestor most likely already had acquired important enzymes for the salvage of nucleotides and had a reduced capacity to synthesize nucleotides, lipids, and amino acids de novo, suggesting that it was an obligate host-associated organism.

IMPORTANCE Diplomonads are a group of microbial eukaryotes found in oxygen-poor environments. There are both parasitic (e.g., Giardia intestinalis) and free-living (e.g., Trepomonas) members in the group. Diplomonads are well known for their anaerobic metabolism, which has been studied for many years. Here, we reconstructed whole metabolic networks of four extant diplomonad species as well as their ancestors, using a bioinformatics approach. We show that the metabolism within the group is under constant change throughout evolutionary time, in response to the environments that the different lineages explore. Both gene losses and gains are responsible for the adaptation processes. Interestingly, it appears that the last Diplomonadida common ancestor had a metabolism that is more similar to extant parasitic than free-living diplomonads. This suggests that the host-associated lifestyle of parasitic diplomonads, such as the human parasite G. intestinalis, is an old evolutionary adaptation.

The compact genome of Giardia muris reveals important steps in the evolution of intestinal protozoan parasites

Diplomonad parasites of the genus Giardia have adapted to colonizing different hosts, most notably the intestinal tract of mammals. The human-pathogenic Giardia species, Giardia intestinalis, has been extensively studied at the genome and gene expression level, but no such information is available for other Giardia species. Comparative data would be particularly valuable for Giardia muris, which colonizes mice and is commonly used as a prototypic in vivo model for investigating host responses to intestinal parasitic infection. Here we report the draft-genome of G. muris. We discovered a highly streamlined genome, amongst the most densely encoded ever described for a nuclear eukaryotic genome. G. muris and G. intestinalis share many known or predicted virulence factors, including cysteine proteases and a large repertoire of cysteine-rich surface proteins involved in antigenic variation. Different to G. intestinalis, G. muris maintains tandem arrays of pseudogenized surface antigens at the telomeres, whereas intact surface antigens are present centrally in the chromosomes. The two classes of surface antigens engage in genetic exchange. Reconstruction of metabolic pathways from the G. muris genome suggest significant metabolic differences to G. intestinalis. Additionally, G. muris encodes proteins that might be used to modulate the prokaryotic microbiota. The responsible genes have been introduced in the Giardia genus via lateral gene transfer from prokaryotic sources. Our findings point to important evolutionary steps in the Giardia genus as it adapted to different hosts and it provides a powerful foundation for mechanistic exploration of host-pathogen interaction in the G. muris-mouse pathosystem.

Whole-genome sequencing of dog-specific assemblages C and D of Giardia duodenalis from single and pooled cysts indicates host-associated genes

Giardia duodenalis (syn. Giardia intestinalis or Giardia lamblia) infSAects over 280 million people each year and numerous animals. G. duodenalis can be subdivided into eight assemblages with different host specificity. Unculturable assemblages have so far resisted genome sequencing efforts. In this study, we isolated single and pooled cysts of assemblages C and D from dog faeces by FACS, and sequenced them using multiple displacement amplification and Illumina paired-end sequencing. The genomes of assemblages C and D were compared with genomes of assemblages A and B from humans and assemblage E from ruminants and pigs. The genomes obtained from the pooled cysts and from the single cysts were considered complete (>99 % marker genes observed) and the allelic sequence heterozygosity (ASH) values of assemblages C and D were 0.89 and 0.74 %, respectively. These ASH values were slightly higher than for assemblage B (>0.43 %) and much higher than for assemblages A and E, which ranged from 0.002 to 0.037 %. The flavohaemoglobin and 4Fe-4S binding domain family encoding genes involved in O₂ and NO detoxification were only present in assemblages A, B and E. Cathepsin B orthologs were found in all genomes. Six clades of cathepsin B orthologs contained one gene of each genome, while in three clades not all assemblages were represented. We conclude that whole-genome sequencing from a single Giardia cyst results in complete draft genomes, making the genomes of unculturable Giardia assemblages accessible. Observed differences between the genomes of assemblages C and D on one hand and the assemblages A, B and E on the other hand are possibly associated with host specificity.

Lateral Acquisitions Repeatedly Remodel the Oxygen Detoxification Pathway in Diplomonads and Relatives

Oxygen and reactive oxygen species (ROS) are important stress factors for cells because they can oxidize many large molecules. Fornicata, a group of flagellated protists that includes diplomonads, have anaerobic metabolism but are still able to tolerate fluctuating levels of oxygen. We identified 25 protein families putatively involved in detoxification of oxygen and ROS in this group using a bioinformatics approach and propose how these interact in an oxygen detoxification pathway. These protein families were divided into a central oxygen detoxification pathway and accessory pathways for the synthesis of nonprotein thiols. We then used a phylogenetic approach to investigate the evolutionary origin of the components of this putative pathway in Diplomonadida and other Fornicata species. Our analyses suggested that the diplomonad ancestor was adapted to low-oxygen levels, was able to reduce O₂ to H₂O in a manner similar to extant diplomonads, and was able to synthesize glutathione and l-cysteine. Several genes involved in the pathway have complex evolutionary histories and have apparently been repeatedly acquired through lateral gene transfer and subsequently lost. At least seven genes were acquired independently in different Fornicata lineages, leading to evolutionary convergences. It is likely that acquiring these oxygen detoxification proteins helped anaerobic organisms (like the parasitic Giardia intestinalis) adapt to low-oxygen environments (such as the digestive tract of aerobic hosts).

Molecular Hydrogen Metabolism: a Widespread Trait of Pathogenic Bacteria and Protists

Pathogenic microorganisms use various mechanisms to conserve energy in host tissues and environmental reservoirs. One widespread but often overlooked means of energy conservation is through the consumption or production of molecular hydrogen (H2). Here, we comprehensively review the distribution, biochemistry, and physiology of H2 metabolism in pathogens. Over 200 pathogens and pathobionts carry genes for hydrogenases, the enzymes responsible for H2 oxidation and/or production. Furthermore, at least 46 of these species have been experimentally shown to consume or produce H2 Several major human pathogens use the large amounts of H2 produced by colonic microbiota as an energy source for aerobic or anaerobic respiration. This process has been shown to be critical for growth and virulence of the gastrointestinal bacteria Salmonella enterica serovar Typhimurium, Campylobacter jejuni, Campylobacter concisus, and Helicobacter pylori (including carcinogenic strains). H2 oxidation is generally a facultative trait controlled by central regulators in response to energy and oxidant availability. Other bacterial and protist pathogens produce H2 as a diffusible end product of fermentation processes. These include facultative anaerobes such as Escherichia coli, S Typhimurium, and Giardia intestinalis, which persist by fermentation when limited for respiratory electron acceptors, as well as obligate anaerobes, such as Clostridium perfringens, Clostridioides difficile, and Trichomonas vaginalis, that produce large amounts of H2 during growth. Overall, there is a rich literature on hydrogenases in growth, survival, and virulence in some pathogens. However, we lack a detailed understanding of H2 metabolism in most pathogens, especially obligately anaerobic bacteria, as well as a holistic understanding of gastrointestinal H2 transactions overall. Based on these findings, we also evaluate H2 metabolism as a possible target for drug development or other therapies.

A chromosome-scale reference genome for Giardia intestinalis WB

Giardia intestinalis is a protist causing diarrhea in humans. The first G. intestinalis genome, from the WB isolate, was published more than ten years ago, and has been widely used as the reference genome for Giardia research. However, the genome is fragmented, thus hindering research at the chromosomal level. We re-sequenced the Giardia genome with Pacbio long-read sequencing technology and obtained a new reference genome, which was assembled into near-complete chromosomes with only four internal gaps at long repeats. This new genome is not only more complete but also better annotated at both structural and functional levels, providing more details about gene families, gene organizations and chromosomal structure. This near-complete reference genome will be a valuable resource for the Giardia community and protist research. It also showcases how a fragmented genome can be improved with long-read sequencing technology completed with optical maps.

Measurement(s)	DNA • sequence_assembly • sequence feature annotation
Technology Type(s)	DNA sequencing • sequence assembly process • sequence annotation
Sample Characteristic - Organism	Giardia intestinalis

Machine-accessible metadata file describing the reported data: 10.6084/m9.figshare.11695659

The Oxymonad Genome Displays Canonical Eukaryotic Complexity in the Absence of a Mitochondrion

The discovery that the protist Monocercomonoides exilis completely lacks mitochondria demonstrates that these organelles are not absolutely essential to eukaryotic cells. However, the degree to which the metabolism and cellular systems of this organism have adapted to the loss of mitochondria is unknown. Here, we report an extensive analysis of the M. exilis genome to address this question. Unexpectedly, we find that M. exilis genome structure and content is similar in complexity to other eukaryotes and less "reduced" than genomes of some other protists from the Metamonada group to which it belongs. Furthermore, the predicted cytoskeletal systems, the organization of endomembrane systems, and biosynthetic pathways also display canonical eukaryotic complexity. The only apparent preadaptation that permitted the loss of mitochondria was the acquisition of the SUF system for Fe-S cluster assembly and the loss of glycine cleavage system. Changes in other systems, including in amino acid metabolism and oxidative stress response, were coincident with the loss of mitochondria but are likely adaptations to the microaerophilic and endobiotic niche rather than the mitochondrial loss per se. Apart from the lack of mitochondria and peroxisomes, we show that M. exilis is a fully elaborated eukaryotic cell that is a promising model system in which eukaryotic cell biology can be investigated in the absence of mitochondria.

Patterns of conservation of spliceosomal intron structures and spliceosome divergence in representatives of the diplomonad and parabasalid lineages

Background

Two spliceosomal intron types co-exist in eukaryotic precursor mRNAs and are excised by distinct U2-dependent and U12-dependent spliceosomes. In the diplomonad Giardia lamblia, small nuclear (sn) RNAs show hybrid characteristics of U2- and U12-dependent spliceosomal snRNAs and 5 of 11 identified remaining spliceosomal introns are trans-spliced. It is unknown whether unusual intron and spliceosome features are conserved in other diplomonads.

Results

We have identified spliceosomal introns, snRNAs and proteins from two additional diplomonads for which genome information is currently available, Spironucleus vortens and Spironucleus salmonicida, as well as relatives, including 6 verified cis-spliceosomal introns in S. vortens. Intron splicing signals are mostly conserved between the Spironucleus species and G. lamblia. Similar to ‘long’ G. lamblia introns, RNA secondary structural potential is evident for ‘long’ (> 50 nt) Spironucleus introns as well as introns identified in the parabasalid Trichomonas vaginalis. Base pairing within these introns is predicted to constrain spatial distances between splice junctions to similar distances seen in the shorter and uniformly-sized introns in these organisms. We find that several remaining Spironucleus spliceosomal introns are ancient. We identified a candidate U2 snRNA from S. vortens, and U2 and U5 snRNAs in S. salmonicida; cumulatively, illustrating significant snRNA differences within some diplomonads. Finally, we studied spliceosomal protein complements and find protein sets in Giardia, Spironucleus and Trepomonas sp. PC1 highly- reduced but well conserved across the clade, with between 44 and 62 out of 174 studied spliceosomal proteins detectable. Comparison with more distant relatives revealed a highly nested pattern, with the more intron-rich fornicate Kipferlia bialata retaining 87 total proteins including nearly all those observed in the diplomonad representatives, and the oxymonad Monocercomonoides retaining 115 total proteins including nearly all those observed in K. bialata.

Conclusions

Comparisons in diplomonad representatives and species of other closely-related metamonad groups indicates similar patterns of intron structural conservation and spliceosomal protein composition but significant divergence of snRNA structure in genomically-reduced species. Relative to other eukaryotes, loss of evolutionarily-conserved snRNA domains and common sets of spliceosomal proteins point to a more streamlined splicing mechanism, where intron sequences and structures may be functionally compensating for the minimalization of spliceosome components.

Electronic supplementary material

The online version of this article (10.1186/s12862-019-1488-y) contains supplementary material, which is available to authorized users.

Proximity Staining Using Enzymatic Protein Tagging in Diplomonads

The function of many proteins is intrinsically related to their cellular location. Novel methods for ascertainment of the ultrastructural location of proteins have been introduced in recent years, but their implementation in protists has so far not been readily realized. Here, we present an optimized proximity labeling protocol using the APEX system in the salmon pathogen Spironucleus salmonicida. This protocol was also applicable to the human pathogen Giardia intestinalis. Both organisms required extraneous addition of hemin to the growth medium to enable detectable peroxidase activity. Further, we saw no inherent limitation in labeling efficiency coupled to the cellular compartment, as evident with some other proximity labeling systems. We anticipate that the APEX proximity labeling system might offer a great resource to establish the ultrastructural localization of proteins across genetically tractable protists but might require organism-specific labeling conditions.

The diplomonads are a group of understudied eukaryotic flagellates whose most prominent member is the human pathogen Giardia intestinalis. Methods commonly used in other eukaryotic model systems often require special optimization in diplomonads due to the highly derived character of their cell biology. We have optimized a proximity labeling protocol using pea ascorbate peroxidase (APEX) as a reporter for transmission electron microscopy (TEM) to enable the study of ultrastructural cellular details in diplomonads. Currently available TEM-compatible tags require light-induced activation (1, 2) or are inactive in many cellular compartments (3), while ascorbate peroxidase has not been shown to have those limitations. Here, we have optimized the in vivo activities of two versions of pea ascorbate peroxidase (APX^W41F and APEX) using the diplomonad fish parasite Spironucleus salmonicida, a relative of G. intestinalis. We exploited the well-known peroxidase substrates, Amplex UltraRed and 3,3′-diaminobenzidine (DAB), to validate the activity of the two tags and argue that APEX is the most stable version to use in Spironucleus salmonicida. Next, we fused APEX to proteins with established localization to evaluate the activity of APEX in different cellular compartments of the diplomonad cell and used Amplex UltraRed as well as antibodies along with superresolution microscopy to confirm the protein-APEX localization. The ultrastructural details of protein-APEX fusions were determined by TEM, and we observed marker activity in all cellular compartments tested when using the DAB substrate. Finally, we show that the optimized conditions established for S. salmonicida can be used in the related diplomonad G. intestinalis.

IMPORTANCE The function of many proteins is intrinsically related to their cellular location. Novel methods for ascertainment of the ultrastructural location of proteins have been introduced in recent years, but their implementation in protists has so far not been readily realized. Here, we present an optimized proximity labeling protocol using the APEX system in the salmon pathogen Spironucleus salmonicida. This protocol was also applicable to the human pathogen Giardia intestinalis. Both organisms required extraneous addition of hemin to the growth medium to enable detectable peroxidase activity. Further, we saw no inherent limitation in labeling efficiency coupled to the cellular compartment, as evident with some other proximity labeling systems. We anticipate that the APEX proximity labeling system might offer a great resource to establish the ultrastructural localization of proteins across genetically tractable protists but might require organism-specific labeling conditions.

Distribution, Evolution, Catalytic Mechanism, and Physiological Functions of the Flavin-Based Electron-Bifurcating NADH-Dependent Reduced Ferredoxin: NADP+ Oxidoreductase

NADH-dependent reduced ferredoxin:NADP⁺ oxidoreductase (Nfn) is an electron-bifurcating enzyme first discovered in the strict anaerobes Clostridium kluyveri and Moorella thermoacetica. In vivo, Nfn catalyzes the endergonic reduction of NADP⁺ with NADH coupled to the exergonic reduction of NADP⁺ with reduced ferredoxin. Most Nfn homologs consist of two subunits, although in certain species Nfn homologs are fused. In contrast to other electron-bifurcating enzymes, Nfn possess a simpler structure. Therefore, Nfn becomes a perfect model to determine the mechanism of flavin-based electron bifurcation, which is a novel energy coupling mode distributed among anaerobic bacteria and archaea. The crystal structures of Nfn from Thermotoga maritima and Pyrococcus furiosus are known, and studies have shown that the FAD molecule of the NfnB (b-FAD) is the site of electron bifurcation, and other cofactors, including a [2Fe2S] cluster, two [4Fe4S] clusters, and the FAD molecule on the NfnA subunit, contribute to electron transfer. Further, the short-lived anionic flavin semiquinone (ASQ) state of b-FAD is essential for electron bifurcation. Nfn homologs are widely distributed among microbes, including bacteria, archaea, and probably eukaryotes, most of which are anaerobes despite that certain species are facultative microbes and even aerobes. Moreover, potential evidence shows that lateral gene transfer may occur in the evolution of this enzyme. Nfn homologs present four different structural patterns, including the well-characterized NfnAB and three different kinds of fused Nfn homologs whose detailed properties have not been characterized. These findings indicate that gene fusion/fission and gene rearrangement may contribute to the evolution of this enzyme. Under physiological conditions, Nfn catalyzes the reduction of NADP⁺ with NADH and reduced ferredoxin, which is then used in certain NADPH-dependent reactions. Deletion of nfn in several microbes causes low growth and redox unbalance and may influence the distribution of fermentation products. It’s also noteworthy that different Nfn homologs perform different functions according to its circumstance. Physiological functions of Nfn indicate that it can be a potential tool in the metabolic engineering of industrial microorganisms, which can regulate the redox potential in vivo.

Oxygen induces the expression of invasion and stress response genes in the anaerobic salmon parasite Spironucleus salmonicida

Background

Spironucleus salmonicida is an anaerobic parasite that can cause systemic infections in Atlantic salmon. Unlike other diplomonad parasites, such as the human pathogen Giardia intestinalis, Spironucleus species can infiltrate the blood stream of their hosts eventually colonizing organs, skin and gills. How this presumed anaerobe can persist and invade oxygenated tissues, despite having a strictly anaerobic metabolism, remains elusive.

Results

To investigate how S. salmonicida response to oxygen stress, we performed RNAseq transcriptomic analyses of cells grown in the presence of oxygen or antioxidant-free medium. We found that over 20% of the transcriptome is differentially regulated in oxygen (1705 genes) and antioxidant-depleted (2280 genes) conditions. These differentially regulated transcripts encode proteins related to anaerobic metabolism, cysteine and Fe-S cluster biosynthesis, as well as a large number of proteins of unknown function. S. salmonicida does not encode genes involved in the classical elements of oxygen metabolism (e.g., catalases, superoxide dismutase, glutathione biosynthesis, oxidative phosphorylation). Instead, we found that genes encoding bacterial-like oxidoreductases were upregulated in response to oxygen stress. Phylogenetic analysis revealed some of these oxygen-responsive genes (e.g., nadh oxidase, rubrerythrin, superoxide reductase) are rare in eukaryotes and likely derived from lateral gene transfer (LGT) events into diplomonads from prokaryotes. Unexpectedly, we observed that many host evasion- and invasion-related genes were also upregulated under oxidative stress suggesting that oxygen might be an important signal for pathogenesis.

Conclusion

While oxygen is toxic for related organisms, such as G. intestinalis, we find that oxygen is likely a gene induction signal for host invasion- and evasion-related pathways in S. salmonicida. These data provide the first molecular evidence for how S. salmonicida could tolerate oxic host environments and demonstrate how LGT can have a profound impact on the biology of anaerobic parasites.

Electronic supplementary material

The online version of this article (10.1186/s12915-019-0634-8) contains supplementary material, which is available to authorized users.

Single cell transcriptomics: moving towards multi-omics

Single-cell multi-omics analysis helps characterize multiple layers of molecular features at a single-cell scale to provide insights into cellular processes and functions.

Diversity and Evolution of Sensor Histidine Kinases in Eukaryotes

Histidine kinases (HKs) are primary sensor proteins that act in cell signaling pathways generically referred to as "two-component systems" (TCSs). TCSs are among the most widely distributed transduction systems used by both prokaryotic and eukaryotic organisms to detect and respond to a broad range of environmental cues. The structure and distribution of HK proteins are now well documented in prokaryotes, but information is still fragmentary for eukaryotes. Here, we have taken advantage of recent genomic resources to explore the structural diversity and the phylogenetic distribution of HKs in the prominent eukaryotic supergroups. Searches of the genomes of 67 eukaryotic species spread evenly throughout the phylogenetic tree of life identified 748 predicted HK proteins. Independent phylogenetic analyses of predicted HK proteins were carried out for each of the major eukaryotic supergroups. This allowed most of the compiled sequences to be categorized into previously described HK groups. Beyond the phylogenetic analysis of eukaryotic HKs, this study revealed some interesting findings: 1) characterization of some previously undescribed eukaryotic HK groups with predicted functions putatively related to physiological traits; 2) discovery of HK groups that were previously believed to be restricted to a single kingdom in additional supergroups, and 3) indications that some evolutionary paths have led to the appearance, transfer, duplication, and loss of HK genes in some phylogenetic lineages. This study provides an unprecedented overview of the structure and distribution of HKs in the Eukaryota and represents a first step toward deciphering the evolution of TCS signaling in living organisms.

Aquatic Parasite Cultures and Their Applications

In this era of unprecedented growth in aquaculture and trade, aquatic parasite cultures are essential to better understand emerging diseases and their implications for human and animal health. Yet culturing parasites presents multiple challenges, arising from their complex, often multihost life cycles, multiple developmental stages, variable generation times and reproductive modes. Furthermore, the essential environmental requirements of most parasites remain enigmatic. Despite these inherent difficulties, in vivo and in vitro cultures are being developed for a small but growing number of aquatic pathogens. Expanding this resource will facilitate diagnostic capabilities and treatment trials, thus supporting the growth of sustainable aquatic commodities and communities.

Fe-S cluster assembly in the supergroup Excavata

The majority of established model organisms belong to the supergroup Opisthokonta, which includes yeasts and animals. While enlightening, this focus has neglected protists,  organisms that represent the bulk of eukaryotic diversity and are often regarded as primitive eukaryotes. One of these is the “supergroup” Excavata, which comprises unicellular flagellates of diverse lifestyles and contains species of medical importance, such as Trichomonas, Giardia, Naegleria, Trypanosoma and Leishmania. Excavata exhibits a continuum in mitochondrial forms, ranging from classical aerobic, cristae-bearing mitochondria to mitochondria-related organelles, such as hydrogenosomes and mitosomes, to the extreme case of a complete absence of the organelle. All forms of mitochondria house a machinery for the assembly of Fe–S clusters, ancient cofactors required in various biochemical activities needed to sustain every extant cell. In this review, we survey what is known about the Fe–S cluster assembly in the supergroup Excavata. We aim to bring attention to the diversity found in this group, reflected in gene losses and gains that have shaped the Fe–S cluster biogenesis pathways.

The draft genome of Kipferlia bialata reveals reductive genome evolution in fornicate parasites

The fornicata (fornicates) is a eukaryotic group known to consist of free-living and parasitic organisms. Genome datasets of two model fornicate parasites Giardia intestinalis and Spironucleus salmonicida are well annotated, so far. The nuclear genomes of G. intestinalis assemblages and S. salmonicida are small in terms of the genome size and simple in genome structure. However, an ancestral genomic structure and gene contents, from which genomes of the fornicate parasites have evolved, remains to be clarified. In order to understand genome evolution in fornicates, here, we present the draft genome sequence of a free-living fornicate, Kipferlia bialata, the divergence of which is earlier than those of the fornicate parasites, and compare it to the genomes of G. intestinalis and S. salmonicida. Our data show that the number of protein genes and introns in K. bialata genome are the most abundant in the genomes of three fornicates, reflecting an ancestral state of fornicate genome evolution. Evasion mechanisms of host immunity found in G. intestinalis and S. salmonicida are absent in the K. bialata genome, suggesting that the two parasites acquired the complex membrane surface proteins on the line leading to the common ancestor of G. intestinalis and S. salmonicida after the divergence from K. bialata. Furthermore, the mitochondrion related organelles (MROs) of K. bialata possess more complex suites of metabolic pathways than those in Giardia and in Spironucleus. In sum, our results unveil the process of reductive evolution which shaped the current genomes in two model fornicate parasites G. intestinalis and S. salmonicida.

Diversity within diversity: Parasite species richness in poison frogs assessed by transcriptomics

Symbionts (e.g., endoparasites and commensals) play an integral role in their host's ecology, yet in many cases their diversity is likely underestimated. Although endoparasites are traditionally characterized using morphology, sequences of conserved genes, and shotgun metagenomics, host transcriptomes constitute an underused resource to identify these organisms' diversity. By isolating non-host transcripts from host transcriptomes, individual host tissues can now simultaneously reveal their endoparasite species richness (i.e., number of different taxa) and provide insights into parasite gene expression. These approaches can be used in host taxa whose endoparasites are mostly unknown, such as those of tropical amphibians. Here, we focus on the poison frogs (Dendrobatidae) as hosts, which are a Neotropical clade known for their bright coloration and defensive alkaloids. These toxins are an effective protection against vertebrate predators (e.g., snakes and birds), bacteria, and skin-biting ectoparasites (e.g., mosquitoes); however, little is known about their deterrence against eukaryotic endoparasites. With de novo transcriptomes of dendrobatids, we developed a bioinformatics pipeline for endoparasite identification that uses host annotated RNA-seq data and set of a priori parasite taxonomic terms, which are used to mine for specific endoparasites. We found a large community of helminths and protozoans that were mostly restricted to the digestive tract and a few systemic parasites (e.g., Trypanosoma). Contrary to our expectations, all dendrobatid frogs regardless of the presence of alkaloid defenses have endoparasites, with their highest species richness located in the frog digestive tract. Some of these organisms (e.g., roundworms) might prove to be generalists, as they were not found to be co-diversifying with their frog hosts. We propose that endoparasites may escape poison frogs' chemical defenses by colonizing tissues with fewer alkaloids than the frog's skin, where most toxins are stored.

14-3-3 Regulates Actin Filament Formation in the Deep-Branching Eukaryote Giardia lamblia

The phosphoserine/phosphothreonine-binding protein 14-3-3 is known to regulate actin; this function has been previously attributed to sequestration of phosphorylated cofilin. 14-3-3 was identified as an actin-associated protein in the deep-branching eukaryote Giardia lamblia; however, Giardia lacks cofilin and all other canonical actin-binding proteins (ABPs). Thus, the role of G. lamblia 14-3-3 (Gl-14-3-3) in actin regulation was unknown. Gl-14-3-3 depletion resulted in an overall disruption of actin organization characterized by ectopically distributed short actin filaments. Using phosphatase and kinase inhibitors, we demonstrated that actin phosphorylation correlated with destabilization of the actin network and increased complex formation with 14-3-3, while blocking actin phosphorylation stabilized actin filaments and attenuated complex formation. Giardia's sole Rho family GTPase, Gl-Rac, modulates Gl-14-3-3's association with actin, providing the first connection between Gl-Rac and the actin cytoskeleton in Giardia. Giardia actin (Gl-actin) contains two putative 14-3-3 binding motifs, one of which (S330) is conserved in mammalian actin. Mutation of these sites reduced, but did not completely disrupt, the association with 14-3-3. Native gels and overlay assays indicate that intermediate proteins are required to support complex formation between 14-3-3 and actin. Overall, our results support a role for 14-3-3 as a regulator of actin; however, the presence of multiple 14-3-3-actin complexes suggests a more complex regulatory relationship than might be expected for a minimalistic parasite. IMPORTANCEGiardia lacks canonical actin-binding proteins. Gl-14-3-3 was identified as an actin interactor, but the significance of this interaction was unknown. Loss of Gl-14-3-3 results in ectopic short actin filaments, indicating that Gl-14-3-3 is an important regulator of the actin cytoskeleton in Giardia. Drug studies indicate that Gl-14-3-3 complex formation is in part phospho-regulated. We demonstrate that complex formation is downstream of Giardia's sole Rho family GTPase, Gl-Rac. This result provides the first mechanistic connection between Gl-Rac and Gl-actin in Giardia. Native gels and overlay assays indicate intermediate proteins are required to support the interaction between Gl-14-3-3 and Gl-actin, suggesting that Gl-14-3-3 is regulating multiple Gl-actin complexes.

The Macronuclear Genome of Stentor coeruleus Reveals Tiny Introns in a Giant Cell

The giant, single-celled organism Stentor coeruleus has a long history as a model system for studying pattern formation and regeneration in single cells. Stentor [1, 2] is a heterotrichous ciliate distantly related to familiar ciliate models, such as Tetrahymena or Paramecium. The primary distinguishing feature of Stentor is its incredible size: a single cell is 1 mm long. Early developmental biologists, including T.H. Morgan [3], were attracted to the system because of its regenerative abilities-if large portions of a cell are surgically removed, the remnant reorganizes into a normal-looking but smaller cell with correct proportionality [2, 3]. These biologists were also drawn to Stentor because it exhibits a rich repertoire of behaviors, including light avoidance, mechanosensitive contraction, food selection, and even the ability to habituate to touch, a simple form of learning usually seen in higher organisms [4]. While early microsurgical approaches demonstrated a startling array of regenerative and morphogenetic processes in this single-celled organism, Stentor was never developed as a molecular model system. We report the sequencing of the Stentor coeruleus macronuclear genome and reveal key features of the genome. First, we find that Stentor uses the standard genetic code, suggesting that ciliate-specific genetic codes arose after Stentor branched from other ciliates. We also discover that ploidy correlates with Stentor's cell size. Finally, in the Stentor genome, we discover the smallest spliceosomal introns reported for any species. The sequenced genome opens the door to molecular analysis of single-cell regeneration in Stentor.

Transcriptomic analysis of diplomonad parasites reveals a trans-spliced intron in a helicase gene in Giardia

Background

The mechanisms by which DNA sequences are expressed is the central preoccupation of molecular genetics. Recently, ourselves and others reported that in the diplomonad protist Giardia lamblia, the coding regions of several mRNAs are produced by ligation of independent RNA species expressed from distinct genomic loci. Such trans-splicing of introns was found to affect nearly as many genes in this organism as does classical cis-splicing of introns. These findings raised questions about the incidence of intron trans-splicing both across the G. lambliatranscriptome and across diplomonad diversity in general, however a dearth of transcriptomic data at the time prohibited systematic study of these questions.

Methods

I leverage newly available transcriptomic data from G. lamblia and the related diplomonad Spironucleus salmonicidato search for trans-spliced introns. My computational pipeline recovers all four previously reported trans-spliced introns in G. lamblia, suggesting good sensitivity.

Results

Scrutiny of thousands of potential cases revealed only a single additional trans-spliced intron in G. lamblia, in the p68 helicase gene, and no cases in S. salmonicida. The p68 intron differs from the previously reported trans-spliced introns in its high degree of streamlining: the core features of G. lamblia trans-spliced introns are closely packed together, revealing striking economy in the implementation of a seemingly inherently uneconomical molecular mechanism.

Discussion

These results serve to circumscribe the role of trans-splicing in diplomonads both in terms of the number of genes effected and taxonomically. Future work should focus on the molecular mechanisms, evolutionary origins and phenotypic implications of this intriguing phenomenon.

An Interactome-Centered Protein Discovery Approach Reveals Novel Components Involved in Mitosome Function and Homeostasis in Giardia lamblia

Protozoan parasites of the genus Giardia are highly prevalent globally, and infect a wide range of vertebrate hosts including humans, with proliferation and pathology restricted to the small intestine. This narrow ecological specialization entailed extensive structural and functional adaptations during host-parasite co-evolution. An example is the streamlined mitosomal proteome with iron-sulphur protein maturation as the only biochemical pathway clearly associated with this organelle. Here, we applied techniques in microscopy and protein biochemistry to investigate the mitosomal membrane proteome in association to mitosome homeostasis. Live cell imaging revealed a highly immobilized array of 30–40 physically distinct mitosome organelles in trophozoites. We provide direct evidence for the single giardial dynamin-related protein as a contributor to mitosomal morphogenesis and homeostasis. To overcome inherent limitations that have hitherto severely hampered the characterization of these unique organelles we applied a novel interaction-based proteome discovery strategy using forward and reverse protein co-immunoprecipitation. This allowed generation of organelle proteome data strictly in a protein-protein interaction context. We built an initial Tom40-centered outer membrane interactome by co-immunoprecipitation experiments, identifying small GTPases, factors with dual mitosome and endoplasmic reticulum (ER) distribution, as well as novel matrix proteins. Through iterative expansion of this protein-protein interaction network, we were able to i) significantly extend this interaction-based mitosomal proteome to include other membrane-associated proteins with possible roles in mitosome morphogenesis and connection to other subcellular compartments, and ii) identify novel matrix proteins which may shed light on mitosome-associated metabolic functions other than Fe-S cluster biogenesis. Functional analysis also revealed conceptual conservation of protein translocation despite the massive divergence and reduction of protein import machinery in Giardia mitosomes.

Organelles with endosymbiotic origin are present in virtually all extant eukaryotes and have undergone considerable remodeling during > 1 billion years of evolution. Highly diverged organelles such as mitosomes or plastids in some parasitic protozoa are the product of extensive secondary reduction. They are sufficiently unique to generate interest as targets for pharmacological intervention, in addition to providing a rich ground for evolutionary cell biologists. The so-called mitochondria-related organelles (MROs) comprise mitosomes and hydrogenosomes, with the former having lost any role in energy metabolism along with the organelle genome. The mitosomes of the intestinal pathogen Giardia lamblia are the most highly reduced MROs known and have proven difficult to investigate because of their extreme divergence and their unique biophysical properties. Here, we implemented a novel strategy aimed at systematic analysis of the organelle proteome by iterative expansion of a protein-protein interaction network. We combined serial forward and reverse co-immunoprecipitations with mass spectrometry analysis, data mining, and validation by subcellular localization and/or functional analysis to generate an interactome network centered on a giardial Tom40 homolog. This iterative ab initio proteome reconstruction provided protein-protein interaction data in addition to identifying novel organelle proteins and functions. Building on this data we generated information on organelle replication, mitosome morphogenesis and organelle dynamics in living cells.

Fine-Tuning Motile Cilia and Flagella: Evolution of the Dynein Motor Proteins from Plants to Humans at High Resolution

The flagellum is a key innovation linked to eukaryogenesis. It provides motility by regulated cycles of bending and bend propagation, which are thought to be controlled by a complex arrangement of seven distinct dyneins in repeated patterns of outer- (OAD) and inner-arm dynein (IAD) complexes. Electron tomography showed high similarity of this axonemal repeat pattern across ciliates, algae, and animals, but the diversity of dynein sequences across the eukaryotes has not yet comprehensively been resolved and correlated with structural data. To shed light on the evolution of the axoneme I performed an exhaustive analysis of dyneins using the available sequenced genome data. Evidence from motor domain phylogeny allowed expanding the current set of nine dynein subtypes by eight additional isoforms with, however, restricted taxonomic distributions. I confirmed the presence of the nine dyneins in all eukaryotic super-groups indicating their origin predating the last eukaryotic common ancestor. The comparison of the N-terminal tail domains revealed a most likely axonemal dynein origin of the new classes, a group of chimeric dyneins in plants/algae and Stramenopiles, and the unique domain architecture and origin of the outermost OADs present in green algae and ciliates but not animals. The correlation of sequence and structural data suggests the single-headed class-8 and class-9 dyneins to localize to the distal end of the axonemal repeat and the class-7 dyneins filling the region up to the proximal heterodimeric IAD. Tracing dynein gene duplications across the eukaryotes indicated ongoing diversification and fine-tuning of flagellar functions in extant taxa and species.

Minimal cytosolic iron-sulfur cluster assembly machinery of Giardia intestinalis is partially associated with mitosomes

Iron-sulfur (Fe-S) clusters are essential cofactors that enable proteins to transport electrons, sense signals, or catalyze chemical reactions. The maturation of dozens of Fe-S proteins in various compartments of every eukaryotic cell is driven by several assembly pathways. The ubiquitous cytosolic Fe-S cluster assembly (CIA) pathway, typically composed of eight highly conserved proteins, depends on mitochondrial Fe-S cluster assembly (ISC) machinery. Giardia intestinalis contains one of the smallest eukaryotic genomes and the mitosome, an extremely reduced mitochondrion. Because the only pathway known to be retained within this organelle is the synthesis of Fe-S clusters mediated by ISC machinery, a likely function of the mitosome is to cooperate with the CIA pathway. We investigated the cellular localization of CIA components in G. intestinalis and the origin and distribution of CIA-related components and Tah18-like proteins in other Metamonada. We show that orthologs of Tah18 and Dre2 are missing in these eukaryotes. In Giardia, all CIA components are exclusively cytosolic, with the important exception of Cia2 and two Nbp35 paralogs, which are present in the mitosomes. We propose that the dual localization of Cia2 and Nbp35 proteins in Giardia might represent a novel connection between the ISC and the CIA pathways.

On the reversibility of parasitism: adaptation to a free-living lifestyle via gene acquisitions in the diplomonad Trepomonas sp. PC1

Background

It is generally thought that the evolutionary transition to parasitism is irreversible because it is associated with the loss of functions needed for a free-living lifestyle. Nevertheless, free-living taxa are sometimes nested within parasite clades in phylogenetic trees, which could indicate that they are secondarily free-living. Herein, we test this hypothesis by studying the genomic basis for evolutionary transitions between lifestyles in diplomonads, a group of anaerobic eukaryotes. Most described diplomonads are intestinal parasites or commensals of various animals, but there are also free-living diplomonads found in oxygen-poor environments such as marine and freshwater sediments. All these nest well within groups of parasitic diplomonads in phylogenetic trees, suggesting that they could be secondarily free-living.

Results

We present a transcriptome study of Trepomonas sp. PC1, a diplomonad isolated from marine sediment. Analysis of the metabolic genes revealed a number of proteins involved in degradation of the bacterial membrane and cell wall, as well as an extended set of enzymes involved in carbohydrate degradation and nucleotide metabolism. Phylogenetic analyses showed that most of the differences in metabolic capacity between free-living Trepomonas and the parasitic diplomonads are due to recent acquisitions of bacterial genes via gene transfer. Interestingly, one of the acquired genes encodes a ribonucleotide reductase, which frees Trepomonas from the need to scavenge deoxyribonucleosides. The transcriptome included a gene encoding squalene-tetrahymanol cyclase. This enzyme synthesizes the sterol substitute tetrahymanol in the absence of oxygen, potentially allowing Trepomonas to thrive under anaerobic conditions as a free-living bacterivore, without depending on sterols from other eukaryotes.

Conclusions

Our findings are consistent with the phylogenetic evidence that the last common ancestor of diplomonads was dependent on a host and that Trepomonas has adapted secondarily to a free-living lifestyle. We believe that similar studies of other groups where free-living taxa are nested within parasites could reveal more examples of secondarily free-living eukaryotes.

Electronic supplementary material

The online version of this article (doi:10.1186/s12915-016-0284-z) contains supplementary material, which is available to authorized users.

Static Clathrin Assemblies at the Peripheral Vacuole-Plasma Membrane Interface of the Parasitic Protozoan Giardia lamblia

Giardia lamblia is a parasitic protozoan that infects a wide range of vertebrate hosts including humans. Trophozoites are non-invasive but associate tightly with the enterocyte surface of the small intestine. This narrow ecological specialization entailed extensive morphological and functional adaptations during host-parasite co-evolution, including a distinctly polarized array of endocytic organelles termed peripheral vacuoles (PVs), which are confined to the dorsal cortical region exposed to the gut lumen and are in close proximity to the plasma membrane (PM). Here, we investigated the molecular consequences of these adaptations on the Giardia endocytic machinery and membrane coat complexes. Despite the absence of canonical clathrin coated vesicles in electron microscopy, Giardia possesses conserved PV-associated clathrin heavy chain (GlCHC), dynamin-related protein (GlDRP), and assembly polypeptide complex 2 (AP2) subunits, suggesting a novel function for GlCHC and its adaptors. We found that, in contrast to GFP-tagged AP2 subunits and DRP, CHC::GFP reporters have no detectable turnover in living cells, indicating fundamental differences in recruitment to the membrane and disassembly compared to previously characterized clathrin coats. Histochemical localization in electron tomography showed that these long-lived GlCHC assemblies localized at distinctive approximations between the plasma and PV membrane. A detailed protein interactome of GlCHC revealed all of the conserved factors in addition to novel or highly diverged proteins, including a putative clathrin light chain and lipid-binding proteins. Taken together, our data provide strong evidence for giardial CHC as a component of highly stable assemblies at PV-PM junctions that likely have a central role in organizing continuities between the PM and PV membranes for controlled sampling of the fluid environment. This suggests a novel function for CHC in Giardia and the extent of molecular remodeling of endocytosis in this species.

In canonical clathrin mediated endocytosis (CME) models, the concerted action of ca. 50 proteins mediates the uptake of extracellular components. The key player in this process is clathrin which coats transport intermediates called clathrin coated vesicles (CCV). The intestinal parasite Giardia lamblia has undergone extensive remodeling during colonization of the mammalian duodenum. Here, we report on unique features of this parasite’s endocytic system, consisting of fixed peripheral vacuoles (PV) in close proximity to the exposed plasma membrane (PM), with no discernible CCVs. Using state-of-the-art imaging strategies, we show that the surface of Giardia trophozoites is pock-marked with PM invaginations reaching to the underlying PV membrane. Co-immunoprecipitation and analysis of protein dynamics reveal that, in line with the absence of CCVs, giardial clathrin assemblies have no dynamic behavior. CHC still remains associated to AP2 and dynamin, both conserved dynamic CME components, and to a newly identified putative clathrin light chain. The emerging model calls for giardial clathrin organized into static cores surrounded by dynamic interaction partners, and most likely involved in the regulation of fusion between the PM and the PVs in a “kiss-and-flush”-like mechanism. This suggests that Giardia harbors a conceptually novel function for clathrin in endocytosis, which might be a consequence of host-parasite co-evolution.

A Eukaryote without a Mitochondrial Organelle

The presence of mitochondria and related organelles in every studied eukaryote supports the view that mitochondria are essential cellular components. Here, we report the genome sequence of a microbial eukaryote, the oxymonad Monocercomonoides sp., which revealed that this organism lacks all hallmark mitochondrial proteins. Crucially, the mitochondrial iron-sulfur cluster assembly pathway, thought to be conserved in virtually all eukaryotic cells, has been replaced by a cytosolic sulfur mobilization system (SUF) acquired by lateral gene transfer from bacteria. In the context of eukaryotic phylogeny, our data suggest that Monocercomonoides is not primitively amitochondrial but has lost the mitochondrion secondarily. This is the first example of a eukaryote lacking any form of a mitochondrion, demonstrating that this organelle is not absolutely essential for the viability of a eukaryotic cell.

Comparative Cell Biology and Evolution of Annexins in Diplomonads

Annexins are proteins that associate with phospholipids in a Ca²⁺-dependent fashion. These proteins have been intensely studied in animals and plants because of their importance in diverse cellular processes, yet very little is known about annexins in single-celled eukaryotes, which represent the largest diversity of organisms. The human intestinal parasite Giardia intestinalis is known to have more annexins than humans, and they contribute to its pathogenic potential. In this study, we investigated the annexin complement in the salmon pathogen Spironucleus salmonicida, a relative of G. intestinalis. We found that S. salmonicida has a large repertoire of annexins and that the gene family has expanded separately across diplomonads, with members showing sequence diversity similar to that seen across kingdom-level groups such as plants and animals. S. salmonicida annexins are prominent components of the cytoskeleton and membrane. Two annexins are associated with a previously unrecognized structure in the anterior of the cell.

Annexins are multifunctional, calcium-binding proteins found in organisms across all kingdoms. Most studies of annexins from single-celled eukaryotes have focused on the alpha-giardins, proteins assigned to the group E annexins, expressed by the diplomonad Giardia intestinalis. We have characterized the annexin gene family in another diplomonad parasite, Spironucleus salmonicida, by phylogenetic and experimental approaches. We constructed a comprehensive phylogeny of the diplomonad group E annexins and found that they are abundant across the group with frequent gene duplications and losses. The annexins of S. salmonicida were found to be related to alpha-giardins but with better-preserved type II Ca²⁺ coordination sites. Two annexins were confirmed to bind phospholipids in a Ca²⁺-dependent fashion but with different specificities. Superresolution and confocal microscopy of epitope-tagged S. salmonicida annexins revealed localization to distinct parts of the cytoskeleton and membrane. The ultrastructural details of the localization of several annexins were determined by proximity labeling and transmission electron microscopy. Two annexins localize to a novel cytoskeletal structure in the anterior of the cell. Our results show that the annexin gene family is expanded in diplomonads and that these group E annexins are associated mostly with cytoskeletal and membrane structures.

IMPORTANCE Annexins are proteins that associate with phospholipids in a Ca²⁺-dependent fashion. These proteins have been intensely studied in animals and plants because of their importance in diverse cellular processes, yet very little is known about annexins in single-celled eukaryotes, which represent the largest diversity of organisms. The human intestinal parasite Giardia intestinalis is known to have more annexins than humans, and they contribute to its pathogenic potential. In this study, we investigated the annexin complement in the salmon pathogen Spironucleus salmonicida, a relative of G. intestinalis. We found that S. salmonicida has a large repertoire of annexins and that the gene family has expanded separately across diplomonads, with members showing sequence diversity similar to that seen across kingdom-level groups such as plants and animals. S. salmonicida annexins are prominent components of the cytoskeleton and membrane. Two annexins are associated with a previously unrecognized structure in the anterior of the cell.

Time-Dependent Transcriptional Changes in Axenic Giardia duodenalis Trophozoites

Giardia duodenalis is the most common gastrointestinal protozoan parasite of humans and a significant contributor to the global burden of both diarrheal disease and post-infectious chronic disorders. Although G. duodenalis can be cultured axenically, significant gaps exist in our understanding of the molecular biology and metabolism of this pathogen. The present study employed RNA sequencing to characterize the mRNA transcriptome of G. duodenalis trophozoites in axenic culture, at log (48 h of growth), stationary (60 h), and declining (96 h) growth phases. Using ~400-times coverage of the transcriptome, we identified 754 differentially transcribed genes (DTGs), mainly representing two large DTG groups: 438 that were down-regulated in the declining phase relative to log and stationary phases, and 281 that were up-regulated. Differential transcription of prominent antioxidant and glycolytic enzymes implicated oxygen tension as a key factor influencing the transcriptional program of axenic trophozoites. Systematic bioinformatic characterization of numerous DTGs encoding hypothetical proteins of unknown function was achieved using structural homology searching. This powerful approach greatly informed the differential transcription analysis and revealed putative novel antioxidant-coding genes, and the presence of a near-complete two-component-like signaling system that may link cytosolic redox or metabolite sensing to the observed transcriptional changes. Motif searching applied to promoter regions of the two large DTG groups identified different putative transcription factor-binding motifs that may underpin global transcriptional regulation. This study provides new insights into the drivers and potential mediators of transcriptional variation in axenic G. duodenalis and provides context for static transcriptional studies.

Diversity and origins of anaerobic metabolism in mitochondria and related organelles

Across the diversity of life, organisms have evolved different strategies to thrive in hypoxic environments, and microbial eukaryotes (protists) are no exception. Protists that experience hypoxia often possess metabolically distinct mitochondria called mitochondrion-related organelles (MROs). While there are some common metabolic features shared between the MROs of distantly related protists, these organelles have evolved independently multiple times across the breadth of eukaryotic diversity. Until recently, much of our knowledge regarding the metabolic potential of different MROs was limited to studies in parasitic lineages. Over the past decade, deep-sequencing studies of free-living anaerobic protists have revealed novel configurations of metabolic pathways that have been co-opted for life in low oxygen environments. Here, we provide recent examples of anaerobic metabolism in the MROs of free-living protists and their parasitic relatives. Additionally, we outline evolutionary scenarios to explain the origins of these anaerobic pathways in eukaryotes.

Probing the Biology of Giardia intestinalis Mitosomes Using In Vivo Enzymatic Tagging

Giardia intestinalis parasites contain mitosomes, one of the simplest mitochondrion-related organelles. Strategies to identify the functions of mitosomes have been limited mainly to homology detection, which is not suitable for identifying species-specific proteins and their functions. An in vivo enzymatic tagging technique based on the Escherichia coli biotin ligase (BirA) has been introduced to G. intestinalis; this method allows for the compartment-specific biotinylation of a protein of interest. Known proteins involved in the mitosomal protein import were in vivo tagged, cross-linked, and used to copurify complexes from the outer and inner mitosomal membranes in a single step. New proteins were then identified by mass spectrometry. This approach enabled the identification of highly diverged mitosomal Tim44 (GiTim44), the first known component of the mitosomal inner membrane translocase (TIM). In addition, our subsequent bioinformatics searches returned novel diverged Tim44 paralogs, which mediate the translation and mitosomal insertion of mitochondrially encoded proteins in other eukaryotes. However, most of the identified proteins are specific to G. intestinalis and even absent from the related diplomonad parasite Spironucleus salmonicida, thus reflecting the unique character of the mitosomal metabolism. The in vivo enzymatic tagging also showed that proteins enter the mitosome posttranslationally in an unfolded state and without vesicular transport.

Prevalent structural disorder carries signature of prokaryotic adaptation to oxic atmosphere

Microbes have adopted efficient mechanisms to contend with environmental changes. The emergence of oxygen was a major event that led to an abrupt change in Earth's atmosphere. To adjust with this shift in environmental condition ancient microbes must have undergone several modifications. Although some proteomic and genomic attributes were proposed to facilitate survival of microorganisms in the presence of oxygen, the process of adaptation still remains elusive. Recent studies have focused that intrinsically disordered proteins play crucial roles in adaptation to a wide range of ecological conditions. Therefore, it is likely that disordered proteins could also play indispensable roles in microbial adaptation to the aerobic environment. To test this hypothesis we measured the disorder content of 679 prokaryotes from four oxygen requirement groups. Our result revealed that aerobic proteomes are endowed with the highest protein disorder followed by facultative microbes. Minimal disorder was observed in anaerobic and microaerophilic microbes with no significant difference in their disorder content. Considering all the potential confounding factors that can modulate protein disorder, here we established that the high protein disorder in aerobic microbe is not a by-product of adaptation to any other selective pressure. On the functional level, we found that the high disorder in aerobic proteomes has been utilized for processes that are important for their aerobic lifestyle. Moreover, aerobic proteomes were found to be enriched with disordered binding sites and to contain transcription factors with high disorder propensity. Based on our results, here we proposed that the high protein disorder is an adaptive opportunity for aerobic microbes to fit with the genomic and functional complexities of the aerobic lifestyle.

Identification of obscure yet conserved actin-associated proteins in Giardia lamblia

Consistent with its proposed status as an early branching eukaryote, Giardia has the most divergent actin of any eukaryote and lacks core actin regulators. Although conserved actin-binding proteins are missing from Giardia, its actin is utilized similarly to that of other eukaryotes and functions in core cellular processes such as cellular organization, endocytosis, and cytokinesis. We set out to identify actin-binding proteins in Giardia using affinity purification coupled with mass spectroscopy (multidimensional protein identification technology [MudPIT]) and have identified >80 putative actin-binding proteins. Several of these have homology to conserved proteins known to complex with actin for functions in the nucleus and flagella. We validated localization and interaction for seven of these proteins, including 14-3-3, a known cytoskeletal regulator with a controversial relationship to actin. Our results indicate that although Giardia lacks canonical actin-binding proteins, there is a conserved set of actin-interacting proteins that are evolutionarily indispensable and perhaps represent some of the earliest functions of the actin cytoskeleton.