Endosomal vesicle fusion machinery is involved with the contractile vacuole in Dictyostelium discoideum.

Contractile vacuoles (CVs), enigmatic osmoregulatory organelles, share common characteristics, such as a requirement for RAB11 and high levels of V-ATPase. These commonalities suggest a conserved evolutionary origin for the CVs with implications for understanding of the last common ancestor of eukaryotes and eukaryotic diversification more broadly. A taxonomically broader sampling of CV-associated machinery is required to address this question further. We used a transcriptomics-based approach to identify CV-associated gene products in Dictyostelium discoideum. This approach was first validated by assessing a set of known CV-associated gene products, which were significantly upregulated following hypo-osmotic exposure. Moreover, endosomal and vacuolar gene products were enriched in the upregulated gene set. An upregulated SNARE protein (NPSNB) was predominantly plasma membrane localised and enriched in the vicinity of CVs, supporting the association with this organelle found in the transcriptomic analysis. We therefore confirm that transcriptomic approaches can identify known and novel players in CV function, in our case emphasizing the role of endosomal vesicle fusion machinery in the D. discoideum CV and facilitating future work to address questions regarding the deep evolution of eukaryotic organelles.

Proksee: in-depth characterization and visualization of bacterial genomes.

Proksee (https://proksee.ca) provides users with a powerful, easy-to-use, and feature-rich system for assembling, annotating, analysing, and visualizing bacterial genomes. Proksee accepts Illumina sequence reads as compressed FASTQ files or pre-assembled contigs in raw, FASTA, or GenBank format. Alternatively, users can supply a GenBank accession or a previously generated Proksee map in JSON format. Proksee then performs assembly (for raw sequence data), generates a graphical map, and provides an interface for customizing the map and launching further analysis jobs. Notable features of Proksee include unique and informative assembly metrics provided via a custom reference database of assemblies; a deeply integrated high-performance genome browser for viewing and comparing analysis results at individual base resolution (developed specifically for Proksee); an ever-growing list of embedded analysis tools whose results can be seamlessly added to the map or searched and explored in other formats; and the option to export graphical maps, analysis results, and log files for data sharing and research reproducibility. All these features are provided via a carefully designed multi-server cloud-based system that can easily scale to meet user demand and that ensures the web server is robust and responsive.

Evolution of factors shaping the endoplasmic reticulum.

Endomembrane system compartments are significant elements in virtually all eukaryotic cells, supporting functions including protein synthesis, post-translational modifications and protein/lipid targeting. In terms of membrane area the endoplasmic reticulum (ER) is the largest intracellular organelle, but the origins of proteins defining the organelle and the nature of lineage-specific modifications remain poorly studied. To understand the evolution of factors mediating ER morphology and function we report a comparative genomics analysis of experimentally characterized ER-associated proteins involved in maintaining ER structure. We find that reticulons, REEPs, atlastins, Ufe1p, Use1p, Dsl1p, TBC1D20, Yip3p and VAPs are highly conserved, suggesting an origin at least as early as the last eukaryotic common ancestor (LECA), although many of these proteins possess additional non-ER functions in modern eukaryotes. Secondary losses are common in individual species and in certain lineages, for example lunapark is missing from the Stramenopiles and the Alveolata. Lineage-specific innovations include protrudin, Caspr1, Arl6IP1, p180, NogoR, kinectin and CLIMP-63, which are restricted to the Opisthokonta. Hence, much of the machinery required to build and maintain the ER predates the LECA, but alternative strategies for the maintenance and elaboration of ER shape and function are present in modern eukaryotes. Moreover, experimental investigations for ER maintenance factors in diverse eukaryotes are expected to uncover novel mechanisms.

Characterisation of eight cattle with Swyer syndrome by whole-genome sequencing.

Swyer syndrome is where an individual has the karyotype of a typical male yet is phenotypically a female. The lack of a (functional) SRY gene located on the Y-chromosome is implicated in some cases of the Swyer syndrome, although many Swyer individuals with an apparently fully functional SRY gene have also been documented. The present study undertook whole genome sequence analyses of eight cattle with suspected Swyer syndrome and compared their genome to that of both a control male and female. Sequence analyses coupled with female phenotypes confirmed that all eight individuals had the 60,XY sex reversal Swyer syndrome. Seven of the eight Swyer syndrome individuals had a deletion on the Y chromosome encompassing the SRY gene (i.e., SRY-). The eighth individual had no obvious mutation in the SRY gene (SRY+) or indeed in any reported gene associated with sex reversal in mammals; a necropsy was performed on this individual. No testicles were detected during the necropsy. Histological examination of the reproductive tract revealed an immature uterine body and horns with inactive glandular tissue of normal histological appearance; both gonads were elongated, a characteristic of most reported cases of Swyer in mammals. The flanking sequence of 11 single nucleotide polymorphisms within 10 kb of the SRY gene are provided to help diagnose some cases of Swyer syndrome. These single nucleotide polymorphisms will not, however, detect all cases of Swyer syndrome since, as evidenced from the present study (and other studies), some individuals with the Swyer condition still contain the SRY gene (i.e., SRY+).

Evaluating the potential of third generation metagenomic sequencing for the detection of BRD pathogens and genetic determinants of antimicrobial resistance in chronically ill feedlot cattle.

BACKGROUND: Bovine respiratory disease (BRD) is an important cause of morbidity and mortality and is responsible for most of the injectable antimicrobial use in the feedlot industry. Traditional bacterial culture can be used to diagnose BRD by confirming the presence of causative pathogens and to support antimicrobial selection. However, given that bacterial culture takes up to a week and early intervention is critical for treatment success, culture has limited utility for informing rapid therapeutic decision-making. In contrast, metagenomic sequencing has the potential to quickly resolve all nucleic acid in a sample, including pathogen biomarkers and antimicrobial resistance genes. In particular, third-generation Oxford Nanopore Technology sequencing platforms provide long reads and access to raw sequencing data in real-time as it is produced, thereby reducing the time from sample collection to diagnostic answer. The purpose of this study was to compare the performance of nanopore metagenomic sequencing to traditional culture and sensitivity methods as applied to nasopharyngeal samples from segregated groups of chronically ill feedlot cattle, previously treated with antimicrobials for nonresponsive pneumonia or lameness. RESULTS: BRD pathogens were isolated from most samples and a variety of different resistance profiles were observed across isolates. The sequencing data indicated the samples were dominated by Moraxella bovoculi, Mannheimia haemolytica, Mycoplasma dispar, and Pasteurella multocida, and included a wide range of antimicrobial resistance genes (ARGs), encoding resistance for up to seven classes of antimicrobials. Genes conferring resistance to beta-lactams were the most commonly detected, while the tetH gene was detected in the most samples overall. Metagenomic sequencing detected the BRD pathogens of interest more often than did culture, but there was limited concordance between phenotypic resistance to antimicrobials and the presence of relevant ARGs. CONCLUSIONS: Metagenomic sequencing can reduce the time from sampling to results, detect pathogens missed by bacterial culture, and identify genetically encoded determinants of resistance. Increasing sequencing coverage of target organisms will be an essential component of improving the reliability of this technology, such that it can be better used for the surveillance of pathogens of interest, genetic determinants of resistance, and to inform diagnostic decisions.

Genomics and transcriptomics yields a system-level view of the biology of the pathogen Naegleria fowleri.

BACKGROUND: The opportunistic pathogen Naegleria fowleri establishes infection in the human brain, killing almost invariably within 2 weeks. The amoeba performs piece-meal ingestion, or trogocytosis, of brain material causing direct tissue damage and massive inflammation. The cellular basis distinguishing N. fowleri from other Naegleria species, which are all non-pathogenic, is not known. Yet, with the geographic range of N. fowleri advancing, potentially due to climate change, understanding how this pathogen invades and kills is both important and timely. RESULTS: Here, we report an -omics approach to understanding N. fowleri biology and infection at the system level. We sequenced two new strains of N. fowleri and performed a transcriptomic analysis of low- versus high-pathogenicity N. fowleri cultured in a mouse infection model. Comparative analysis provides an in-depth assessment of encoded protein complement between strains, finding high conservation. Molecular evolutionary analyses of multiple diverse cellular systems demonstrate that the N. fowleri genome encodes a similarly complete cellular repertoire to that found in free-living N. gruberi. From transcriptomics, neither stress responses nor traits conferred from lateral gene transfer are suggested as critical for pathogenicity. By contrast, cellular systems such as proteases, lysosomal machinery, and motility, together with metabolic reprogramming and novel N. fowleri proteins, are all implicated in facilitating pathogenicity within the host. Upregulation in mouse-passaged N. fowleri of genes associated with glutamate metabolism and ammonia transport suggests adaptation to available carbon sources in the central nervous system. CONCLUSIONS: In-depth analysis of Naegleria genomes and transcriptomes provides a model of cellular systems involved in opportunistic pathogenicity, uncovering new angles to understanding the biology of a rare but highly fatal pathogen.

Regulation of early endosomes across eukaryotes: Evolution and functional homology of Vps9 proteins.

Endocytosis is a crucial process in eukaryotic cells. The GTPases Rab 5, 21 and 22 that mediate endocytosis are ancient eukaryotic features and all available evidence suggests retained conserved function. In animals and fungi, these GTPases are regulated in part by proteins possessing Vps9 domains. However, the diversity, evolution and functions of Vps9 proteins beyond animals or fungi are poorly explored. Here we report a comprehensive analysis of the Vps9 family of GTPase regulators, combining molecular evolutionary data with functional characterization in the non-opisthokont model organism Trypanosoma brucei. At least 3 subfamilies, Alsin, Varp and Rabex5 + GAPVD1, are found across eukaryotes, suggesting that all are ancient features of regulation of endocytic Rab protein function. There are examples of lineage-specific Vps9 subfamily member expansions and novel domain combinations, suggesting diversity in precise regulatory mechanisms between individual lineages. Characterization of the Rabex5 + GAPVD1 and Alsin orthologues in T. brucei demonstrates that both proteins are involved in endocytosis, and that simultaneous knockdown prevents membrane recruitment of Rab5 and Rab21, indicating conservation of function. These data demonstrate that, for the Vps9-domain family at least, modulation of Rab function is mediated by evolutionarily conserved protein-protein interactions.

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.

Identification and characterisation of a cryptic Golgi complex in Naegleria gruberi.

Although the Golgi complex has a conserved morphology of flattened stacked cisternae in most eukaryotes, it has lost the stacked organisation in several lineages, raising the question of what range of morphologies is possible for the Golgi. In order to understand this diversity, it is necessary to characterise the Golgi in many different lineages. Here, we identify the Golgi complex in Naegleria, one of the first descriptions of an unstacked Golgi organelle in a non-parasitic eukaryote, other than fungi. We provide a comprehensive list of Golgi-associated membrane trafficking genes encoded in two species of Naegleria and show that nearly all are expressed in mouse-passaged N. fowleri cells. We then study distribution of the Golgi marker (Ng)CopB by fluorescence in Naegleria gruberi, identifying membranous structures that are disrupted by Brefeldin A treatment, consistent with Golgi localisation. Confocal and immunoelectron microscopy reveals that NgCOPB localises to tubular membranous structures. Our data identify the Golgi organelle for the first time in this major eukaryotic lineage, and provide the rare example of a tubular morphology, representing an important sampling point for the comparative understanding of Golgi organellar diversity.This article has an associated First Person interview with the first author of the paper.

Extreme genome diversity in the hyper-prevalent parasitic eukaryote Blastocystis.

Blastocystis is the most prevalent eukaryotic microbe colonizing the human gut, infecting approximately 1 billion individuals worldwide. Although Blastocystis has been linked to intestinal disorders, its pathogenicity remains controversial because most carriers are asymptomatic. Here, the genome sequence of Blastocystis subtype (ST) 1 is presented and compared to previously published sequences for ST4 and ST7. Despite a conserved core of genes, there is unexpected diversity between these STs in terms of their genome sizes, guanine-cytosine (GC) content, intron numbers, and gene content. ST1 has 6,544 protein-coding genes, which is several hundred more than reported for ST4 and ST7. The percentage of proteins unique to each ST ranges from 6.2% to 20.5%, greatly exceeding the differences observed within parasite genera. Orthologous proteins also display extreme divergence in amino acid sequence identity between STs (i.e., 59%-61% median identity), on par with observations of the most distantly related species pairs of parasite genera. The STs also display substantial variation in gene family distributions and sizes, especially for protein kinase and protease gene families, which could reflect differences in virulence. It remains to be seen to what extent these inter-ST differences persist at the intra-ST level. A full 26% of genes in ST1 have stop codons that are created on the mRNA level by a novel polyadenylation mechanism found only in Blastocystis. Reconstructions of pathways and organellar systems revealed that ST1 has a relatively complete membrane-trafficking system and a near-complete meiotic toolkit, possibly indicating a sexual cycle. Unlike some intestinal protistan parasites, Blastocystis ST1 has near-complete de novo pyrimidine, purine, and thiamine biosynthesis pathways and is unique amongst studied stramenopiles in being able to metabolize α-glucans rather than ß-glucans. It lacks all genes encoding heme-containing cytochrome P450 proteins. Predictions of the mitochondrion-related organelle (MRO) proteome reveal an expanded repertoire of functions, including lipid, cofactor, and vitamin biosynthesis, as well as proteins that may be involved in regulating mitochondrial morphology and MRO/endoplasmic reticulum (ER) interactions. In sharp contrast, genes for peroxisome-associated functions are absent, suggesting Blastocystis STs lack this organelle. Overall, this study provides an important window into the biology of Blastocystis, showcasing significant differences between STs that can guide future experimental investigations into differences in their virulence and clarifying the roles of these organisms in gut health and disease.

Pan genome of the phytoplankton Emiliania underpins its global distribution.

Coccolithophores have influenced the global climate for over 200 million years. These marine phytoplankton can account for 20 per cent of total carbon fixation in some systems. They form blooms that can occupy hundreds of thousands of square kilometres and are distinguished by their elegantly sculpted calcium carbonate exoskeletons (coccoliths), rendering them visible from space. Although coccolithophores export carbon in the form of organic matter and calcite to the sea floor, they also release CO2 in the calcification process. Hence, they have a complex influence on the carbon cycle, driving either CO2 production or uptake, sequestration and export to the deep ocean. Here we report the first haptophyte reference genome, from the coccolithophore Emiliania huxleyi strain CCMP1516, and sequences from 13 additional isolates. Our analyses reveal a pan genome (core genes plus genes distributed variably between strains) probably supported by an atypical complement of repetitive sequence in the genome. Comparisons across strains demonstrate that E. huxleyi, which has long been considered a single species, harbours extensive genome variability reflected in different metabolic repertoires. Genome variability within this species complex seems to underpin its capacity both to thrive in habitats ranging from the equator to the subarctic and to form large-scale episodic blooms under a wide variety of environmental conditions.

Algal genomes reveal evolutionary mosaicism and the fate of nucleomorphs.

Cryptophyte and chlorarachniophyte algae are transitional forms in the widespread secondary endosymbiotic acquisition of photosynthesis by engulfment of eukaryotic algae. Unlike most secondary plastid-bearing algae, miniaturized versions of the endosymbiont nuclei (nucleomorphs) persist in cryptophytes and chlorarachniophytes. To determine why, and to address other fundamental questions about eukaryote-eukaryote endosymbiosis, we sequenced the nuclear genomes of the cryptophyte Guillardia theta and the chlorarachniophyte Bigelowiella natans. Both genomes have >21,000 protein genes and are intron rich, and B. natans exhibits unprecedented alternative splicing for a single-celled organism. Phylogenomic analyses and subcellular targeting predictions reveal extensive genetic and biochemical mosaicism, with both host- and endosymbiont-derived genes servicing the mitochondrion, the host cell cytosol, the plastid and the remnant endosymbiont cytosol of both algae. Mitochondrion-to-nucleus gene transfer still occurs in both organisms but plastid-to-nucleus and nucleomorph-to-nucleus transfers do not, which explains why a small residue of essential genes remains locked in each nucleomorph.

An ER-directed transcriptional response to unfolded protein stress in the absence of conserved sensor-transducer proteins in Giardia lamblia.

The protozoan Giardia lamblia has a minimized organelle repertoire, and most strikingly lacks a classical stacked Golgi apparatus. Nevertheless, Giardia trophozoites constitutively secrete variant surface proteins, and dramatically increase the volume of protein secretion during differentiation to cysts. Eukaryotic cells have evolved an elaborate system for quality control (QC) of protein folding and capacity in the endoplasmic reticulum (ER). Upon ER-overload, an unfolded protein response (UPR) is triggered on transcriptional/translational level aiming at alleviating ER stress. In Giardia, a minimized secretory machinery and absence of glycan-dependent QC suggests that a genetically conserved UPR (or functional equivalent) to cope with insults to the secretory system has been eliminated. We tested this hypothesis of UPR elimination by profiling the transcriptional response during induced ER-folding stress. We show that on the contrary, ER-folding stress triggers a stressor-specific, ER-directed response with upregulation of only ~ 30 genes, with different kinetics and scope compared with the UPR of other eukaryotes. Computational genomics revealed conserved cis-acting motifs in upstream regions of responder genes capable of stressor-specific gene regulation in transfected cells. Interestingly, the sensors/transducers of folding stress, well conserved in model eukaryotes, are absent in Giardia suggesting the presence of a novel version of this essential eukaryotic function.

Bacterial enrichment prior to third-generation metagenomic sequencing improves detection of BRD pathogens and genetic determinants of antimicrobial resistance in feedlot cattle.

Introduction: Bovine respiratory disease (BRD) is one of the most important animal health problems in the beef industry. While bacterial culture and antimicrobial susceptibility testing have been used for diagnostic testing, the common practice of examining one isolate per species does not fully reflect the bacterial population in the sample. In contrast, a recent study with metagenomic sequencing of nasal swabs from feedlot cattle is promising in terms of bacterial pathogen identification and detection of antimicrobial resistance genes (ARGs). However, the sensitivity of metagenomic sequencing was impeded by the high proportion of host biomass in the nasal swab samples. Methods: This pilot study employed a non-selective bacterial enrichment step before nucleic acid extraction to increase the relative proportion of bacterial DNA for sequencing. Results: Non-selective bacterial enrichment increased the proportion of bacteria relative to host sequence data, allowing increased detection of BRD pathogens compared with unenriched samples. This process also allowed for enhanced detection of ARGs with species-level resolution, including detection of ARGs for bacterial species of interest that were not targeted for culture and susceptibility testing. The long-read sequencing approach enabled ARG detection on individual bacterial reads without the need for assembly. Metagenomics following non-selective bacterial enrichment resulted in substantial agreement for four of six comparisons with culture for respiratory bacteria and substantial or better correlation with qPCR. Comparison between isolate susceptibility results and detection of ARGs was best for macrolide ARGs in Mannheimia haemolytica reads but was also substantial for sulfonamide ARGs within M. haemolytica and Pasteurella multocida reads and tetracycline ARGs in Histophilus somni reads. Discussion: By increasing the proportion of bacterial DNA relative to host DNA through non-selective enrichment, we demonstrated a corresponding increase in the proportion of sequencing data identifying BRD-associated pathogens and ARGs in deep nasopharyngeal swabs from feedlot cattle using long-read metagenomic sequencing. This method shows promise as a detection strategy for BRD pathogens and ARGs and strikes a balance between processing time, input costs, and generation of on-target data. This approach could serve as a valuable tool to inform antimicrobial management for BRD and support antimicrobial stewardship.

Multivesicular bodies in the enigmatic amoeboflagellate Breviata anathema and the evolution of ESCRT 0.

Endosomal sorting complexes required for transport (ESCRTs) are heteromeric protein complexes required for multivesicular body (MVB) morphogenesis. ESCRTs I, II, III and III-associated are ubiquitous in eukaryotes and presumably ancient in origin. ESCRT 0 recruits cargo to the MVB and appears to be opisthokont-specific, bringing into question aspects of the current model of ESCRT mechanism. One caveat to the restricted distribution of ESCRT 0 was the previous limited availability of amoebozoan genomes, the supergroup closest to opisthokonts. Here, we significantly expand the sampling of ESCRTs in Amoebozoa. Our electron micrographic and bioinformatics evidence confirm the presence of MVBs in the amoeboflagellate Breviata anathema. Searches of genomic databases of amoebozoans confirm the ubiquitous nature of ESCRTs I-III-associated and the restriction of ESCRT 0 to opisthokonts. Recently, an alternate ESCRT 0 complex, centering on Tom1 proteins, has been proposed. We determine the distribution of Tom1 family proteins across eukaryotes and show that the Tom1, Tom1L1 and Tom1L2 proteins are a vertebrate-specific expansion of the single Tom1 family ancestor, which has indeed been identified in at least one member of each of the major eukaryotic supergroups. This implies a more widely conserved and ancient role for the Tom1 family in endocytosis than previously suspected.

Evolutionary reconstruction of the retromer complex and its function in Trypanosoma brucei.

Intracellular trafficking and protein sorting are mediated by various protein complexes, with the retromer complex being primarily involved in retrograde traffic from the endosome or lysosome to the Golgi complex. Here, comparative genomics, cell biology and phylogenetics were used to probe the early evolution of retromer and its function. Retromer subunits Vps26, Vps29 and Vps35 are near universal, and, by inference, the complex was an ancient feature of eukaryotic cells. Surprisingly, we found DSCR3, a Vps26 paralogue in humans associated with Down's syndrome, in at least four eukaryotic supergroups, implying a more ancient origin than previously suspected. By contrast, retromer cargo proteins showed considerable interlineage variability, with lineage-specific and broadly conserved examples found. Vps10 trafficking probably represents an ancestral role for the complex. Vps5, the BAR-domain-containing membrane-deformation subunit, was found in diverse eukaryotes, including in the divergent eukaryote Trypanosoma brucei, where it is the first example of a BAR-domain protein. To determine functional conservation, an initial characterisation of retromer was performed in T. brucei; the endosomal localisation and its role in endosomal targeting are conserved. Therefore retromer is identified as a further feature of the sophisticated intracellular trafficking machinery of the last eukaryotic common ancestor, with BAR domains representing a possible third independent mechanism of membrane-deformation arising in early eukaryotes.

The mitochondrial genome and a 60-kb nuclear DNA segment from Naegleria fowleri, the causative agent of primary amoebic meningoencephalitis.

Naegleria fowleri is a unicellular eukaryote causing primary amoebic meningoencephalitis, a neuropathic disease killing 99% of those infected, usually within 7-14 days. Naegleria fowleri is found globally in regions including the US and Australia. The genome of the related nonpathogenic species Naegleria gruberi has been sequenced, but the genetic basis for N. fowleri pathogenicity is unclear. To generate such insight, we sequenced and assembled the mitochondrial genome and a 60-kb segment of nuclear genome from N. fowleri. The mitochondrial genome is highly similar to its counterpart in N. gruberi in gene complement and organization, while distinct lack of synteny is observed for the nuclear segments. Even in this short (60-kb) segment, we identified examples of potential factors for pathogenesis, including ten novel N. fowleri-specific genes. We also identified a homolog of cathepsin B; proteases proposed to be involved in the pathogenesis of diverse eukaryotic pathogens, including N. fowleri. Finally, we demonstrate a likely case of horizontal gene transfer between N. fowleri and two unrelated amoebae, one of which causes granulomatous amoebic encephalitis. This initial look into the N. fowleri nuclear genome has revealed several examples of potential pathogenesis factors, improving our understanding of a neglected pathogen of increasing global importance.

Evolutionary analysis of cellular reduction and anaerobicity in the hyper-prevalent gut microbe Blastocystis.

Blastocystis is the most prevalent microbial eukaryote in the human and animal gut, yet its role as commensal or parasite is still under debate. Blastocystis has clearly undergone evolutionary adaptation to the gut environment and possesses minimal cellular compartmentalization, reduced anaerobic mitochondria, no flagella, and no reported peroxisomes. To address this poorly understood evolutionary transition, we have taken a multi-disciplinary approach to characterize Proteromonas lacertae, the closest canonical stramenopile relative of Blastocystis. Genomic data reveal an abundance of unique genes in P. lacertae but also reductive evolution of the genomic complement in Blastocystis. Comparative genomic analysis sheds light on flagellar evolution, including 37 new candidate components implicated with mastigonemes, the stramenopile morphological hallmark. The P. lacertae membrane-trafficking system (MTS) complement is only slightly more canonical than that of Blastocystis, but notably, we identified that both organisms encode the complete enigmatic endocytic TSET complex, a first for the entire stramenopile lineage. Investigation also details the modulation of mitochondrial composition and metabolism in both P. lacertae and Blastocystis. Unexpectedly, we identify in P. lacertae the most reduced peroxisome-derived organelle reported to date, which leads us to speculate on a mechanism of constraint guiding the dynamics of peroxisome-mitochondrion reductive evolution on the path to anaerobiosis. Overall, these analyses provide a launching point to investigate organellar evolution and reveal in detail the evolutionary path that Blastocystis has taken from a canonical flagellated protist to the hyper-divergent and hyper-prevalent animal and human gut microbe.

Mechanism and evolution of the Zn-fingernail required for interaction of VARP with VPS29.

VARP and TBC1D5 are accessory/regulatory proteins of retromer-mediated retrograde trafficking from endosomes. Using an NMR/X-ray approach, we determined the structure of the complex between retromer subunit VPS29 and a 12 residue, four-cysteine/Zn++ microdomain, which we term a Zn-fingernail, two of which are present in VARP. Mutations that abolish VPS29:VARP binding inhibit trafficking from endosomes to the cell surface. We show that VARP and TBC1D5 bind the same site on VPS29 and can compete for binding VPS29 in vivo. The relative disposition of VPS29s in hetero-hexameric, membrane-attached, retromer arches indicates that VARP will prefer binding to assembled retromer coats through simultaneous binding of two VPS29s. The TBC1D5:VPS29 interaction is over one billion years old but the Zn-fingernail appears only in VARP homologues in the lineage directly giving rise to animals at which point the retromer/VARP/TBC1D5 regulatory network became fully established.