Phylogeny-wide analysis of social amoeba genomes highlights ancient origins for complex intercellular communication

Multi-omics analysis of aggregative multicellularity

All organisms have to carefully regulate their gene expression, not least during development. mRNA levels are often used as proxy for protein output; however, this approach ignores post-transcriptional effects. In particular, mRNA-protein correlation remains elusive for organisms that exhibit aggregative rather than clonal multicellularity. We addressed this issue by generating a paired transcriptomics and proteomics time series during the transition from uni-to multicellular stage in the social ameba Dictyostelium discoideum. Our data reveals that mRNA and protein levels correlate highly during unicellular growth, but decrease when multicellular development is initiated. This accentuates that transcripts alone cannot accurately predict protein levels. The dataset provides a useful resource to study gene expression during aggregative multicellular development. Additionally, our study provides an example of how to analyze and visualize mRNA and protein levels, which should be broadly applicable to other organisms and conditions.

Comprehensive characterization of extracellular vesicles produced by environmental (Neff) and clinical (T4) strains of Acanthamoeba castellanii

We conducted a comprehensive comparative analysis of extracellular vesicles (EVs) from two Acanthamoeba castellanii strains, Neff (environmental) and T4 (clinical). Morphological analysis via transmission electron microscopy revealed slightly larger Neff EVs (average = 194.5 nm) compared to more polydisperse T4 EVs (average = 168.4 nm). Nanoparticle tracking analysis (NTA) and dynamic light scattering validated these differences. Proteomic analysis of the EVs identified 1,352 proteins, with 1,107 common, 161 exclusive in Neff, and 84 exclusively in T4 EVs. Gene ontology and Kyoto Encyclopedia of Genes and Genomes (KEGG) mapping revealed distinct molecular functions and biological processes and notably, the T4 EVs enrichment in serine proteases, aligned with its pathogenicity. Lipidomic analysis revealed a prevalence of unsaturated lipid species in Neff EVs, particularly triacylglycerols, phosphatidylethanolamines (PEs), and phosphatidylserine, while T4 EVs were enriched in diacylglycerols and diacylglyceryl trimethylhomoserine, phosphatidylcholine and less unsaturated PEs, suggesting differences in lipid metabolism and membrane permeability. Metabolomic analysis indicated Neff EVs enrichment in glycerolipid metabolism, glycolysis, and nucleotide synthesis, while T4 EVs, methionine metabolism. Furthermore, RNA-seq of EVs revealed differential transcript between the strains, with Neff EVs enriched in transcripts related to gluconeogenesis and translation, suggesting gene regulation and metabolic shift, while in the T4 EVs transcripts were associated with signal transduction and protein kinase activity, indicating rapid responses to environmental changes. In this novel study, data integration highlighted the differences in enzyme profiles, metabolic processes, and potential origins of EVs in the two strains shedding light on the diversity and complexity of A. castellanii EVs and having implications for understanding host-pathogen interactions and developing targeted interventions for Acanthamoeba-related diseases.IMPORTANCEA comprehensive and fully comparative analysis of extracellular vesicles (EVs) from two Acanthamoeba castellanii strains of distinct virulence, a Neff (environmental) and T4 (clinical), revealed striking differences in their morphology and protein, lipid, metabolites, and transcripts levels. Data integration highlighted the differences in enzyme profiles, metabolic processes, and potential distinct origin of EVs from both strains, shedding light on the diversity and complexity of A. castellanii EVs, with direct implications for understanding host-pathogen interactions, disease mechanisms, and developing new therapies for the clinical intervention of Acanthamoeba-related diseases.

Evolution of microRNAs in Amoebozoa and implications for the origin of multicellularity

MicroRNAs (miRNAs) are important and ubiquitous regulators of gene expression in both plants and animals. They are thought to have evolved convergently in these lineages and hypothesized to have played a role in the evolution of multicellularity. In line with this hypothesis, miRNAs have so far only been described in few unicellular eukaryotes. Here, we investigate the presence and evolution of miRNAs in Amoebozoa, focusing on species belonging to Acanthamoeba, Physarum and dictyostelid taxonomic groups, representing a range of unicellular and multicellular lifestyles. miRNAs that adhere to both the stringent plant and animal miRNA criteria were identified in all examined amoebae, expanding the total number of protists harbouring miRNAs from 7 to 15. We found conserved miRNAs between closely related species, but the majority of species feature only unique miRNAs. This shows rapid gain and/or loss of miRNAs in Amoebozoa, further illustrated by a detailed comparison between two evolutionary closely related dictyostelids. Additionally, loss of miRNAs in the Dictyostelium discoideum drnB mutant did not seem to affect multicellular development and, hence, demonstrates that the presence of miRNAs does not appear to be a strict requirement for the transition from uni- to multicellular life.

Costs of being a diet generalist for the protist predator Dictyostelium discoideum

Consumers range from specialists that feed on few resources to generalists that feed on many. Generalism has the clear advantage of having more resources to exploit, but the costs that limit generalism are less clear. We explore two understudied costs of generalism in a generalist amoeba predator, Dictyostelium discoideum, feeding on naturally co-occurring bacterial prey. Both involve costs of combining prey that are suitable on their own. First, amoebas exhibit a reduction in growth rate when they switched to one species of prey bacteria from another compared to controls that experience only the second prey. The effect was consistent across all six tested species of bacteria. These switching costs typically disappear within a day, indicating adjustment to new prey bacteria. This suggests that these costs are physiological. Second, amoebas usually grow more slowly on mixtures of prey bacteria compared to the expectation based on their growth on single prey. There were clear mixing costs in three of the six tested prey mixtures, and none showed significant mixing benefits. These results support the idea that, although amoebas can consume a variety of prey, they must use partially different methods and thus must pay costs to handle multiple prey, either sequentially or simultaneously.

Functional Prokaryotic-Like Deoxycytidine Triphosphate Deaminases and Thymidylate Synthase in Eukaryotic Social Amoebae: Vertical, Endosymbiotic, or Horizontal Gene Transfer?

The de novo synthesis of deoxythymidine triphosphate uses several pathways: gram-negative bacteria use deoxycytidine triphosphate deaminase to convert deoxycytidine triphosphate into deoxyuridine triphosphate, whereas eukaryotes and gram-positive bacteria instead use deoxycytidine monophosphate deaminase to transform deoxycytidine monophosphate to deoxyuridine monophosphate. It is then unusual that in addition to deoxycytidine monophosphate deaminases, the eukaryote Dictyostelium discoideum has 2 deoxycytidine triphosphate deaminases (Dcd1Dicty and Dcd2Dicty). Expression of either DcdDicty can fully rescue the slow growth of an Escherichia coli dcd knockout. Both DcdDicty mitigate the hydroxyurea sensitivity of a Schizosaccharomyces pombe deoxycytidine monophosphate deaminase knockout. Phylogenies show that Dcd1Dicty homologs may have entered the common ancestor of the eukaryotic groups of Amoebozoa, Obazoa, Metamonada, and Discoba through an ancient horizontal gene transfer from a prokaryote or an ancient endosymbiotic gene transfer from a mitochondrion, followed by horizontal gene transfer from Amoebozoa to several other unrelated groups of eukaryotes. In contrast, the Dcd2Dicty homologs were a separate horizontal gene transfer from a prokaryote or a virus into either Amoebozoa or Rhizaria, followed by a horizontal gene transfer between them. ThyXDicty, the D. discoideum thymidylate synthase, another enzyme of the deoxythymidine triphosphate biosynthesis pathway, was suggested previously to be acquired from the ancestral mitochondria or by horizontal gene transfer from alpha-proteobacteria. ThyXDicty can fully rescue the E. coli thymidylate synthase knockout, and we establish that it was obtained by the common ancestor of social amoebae not from mitochondria but from a bacterium. We propose horizontal gene transfer and endosymbiotic gene transfer contributed to the enzyme diversity of the deoxythymidine triphosphate synthesis pathway in most social amoebae, many Amoebozoa, and other eukaryotes.

Nonribosomal peptides protect Pseudomonas nunensis 4A2e from amoebal and nematodal predation

The rhizosphere is a highly competitive environment forcing bacteria to evolve strategies to oppose their enemies. The production of toxic secondary metabolites allows bacteria to counteract predators. In this study, we describe the anti-predator armamentarium of the soil-derived bacterium Pseudomonas nunensis 4A2e. Based on a genome mining approach, we identified several biosynthetic gene clusters coding for nonribosomal peptide synthetases. Generation of gene deletion mutants of the respective clusters shows a loss of defense capabilities. We isolated the novel lipopeptides keanumycin D and nunapeptins B and C, and fully elucidated their structures by a combination of in-depth mass spectrometry experiments, stable isotope labelling, and chemical synthesis. Additionally, investigation of the quorum sensing-dependent biosynthesis allowed us to elucidate parts of the underlying regulation of the biosynthetic machinery. Ecology-inspired bioassays highlight the role of these peptides as a defence strategy against protozoans and led us to find a previously unknown function against the bacterivorous nematode Oscheius myriophilus.

The protein kinases of Dictyostelia and their incorporation into a signalome

Protein kinases are major regulators of cellular processes, but the roles of most kinases remain unresolved. Dictyostelid social amoebas have been useful in identifying functions for 30% of its kinases in cell migration, cytokinesis, vesicle trafficking, gene regulation and other processes but their upstream regulators and downstream effectors are mostly unknown. Comparative genomics can assist to distinguish between genes involved in deeply conserved core processes and those involved in species-specific innovations, while co-expression of genes as evident from comparative transcriptomics can provide cues to the protein complement of regulatory networks. Genomes and developmental and cell-type specific transcriptomes are available for species that span the 0.5 billion years of evolution of Dictyostelia from their unicellular ancestors. In this work we analysed conservation and change in the abundance, functional domain architecture and developmental regulation of protein kinases across the 4 major taxon groups of Dictyostelia. All data are summarized in annotated phylogenetic trees of the kinase subtypes and accompanied by functional information of all kinases that were experimentally studied. We detected 393 different protein kinase domains across the five studied genomes, of which 212 were fully conserved. Conservation was highest (71%) in the previously defined AGC, CAMK, CK1, CMCG, STE and TKL groups and lowest (26%) in the "other" group of typical protein kinases. This was mostly due to species-specific single gene amplification of "other" kinases. Apart from the AFK and α-kinases, the atypical protein kinases, such as the PIKK and histidine kinases were also almost fully conserved. The phylogeny-wide developmental and cell-type specific expression profiles of the protein kinase genes were combined with profiles from the same transcriptomic experiments for the families of G-protein coupled receptors, small GTPases and their GEFs and GAPs, the transcription factors and for all genes that upon lesion generate a developmental defect. This dataset was subjected to hierarchical clustering to identify clusters of co-expressed genes that potentially act together in a signalling network. The work provides a valuable resource that allows researchers to identify protein kinases and other regulatory proteins that are likely to act as intermediates in a network of interest.

Phylogeny-wide analysis of G-protein coupled receptors in social amoebas and implications for the evolution of multicellularity

G-protein coupled receptors (GPCRs) are seven-transmembrane proteins and constitute the largest group of receptors within eukaryotes. The presence of a large set of GPCRs in the unicellular Amoebozoa was surprising and is indicative of the largely undiscovered environmental sensing capabilities in this group. Evolutionary transitions from unicellular to multicellular lifestyles, like we see in social amoebas, have occurred several times independently in the Amoebozoa, and GPCRs may have been co-opted for new functions in cell-cell communication. Methods We have analysed a set of GPCRs from fully sequenced Amoebozoan genomes by Bayesian inference, compared their phylogenetic distribution and domain composition, and analysed their temporal and spatial expression patterns in five species of dictyostelids. Results We found evidence that most GPCRs are conserved deeply in the Amoebozoa and are probably performing roles in general cell functions and complex environmental sensing. All families of GPCRs (apart from the family 4 fungal pheromone receptors) are present in dictyostelids with family 5 being the largest and family 2 the one with the fewest members. For the first time, we identify the presence of family 1 rhodopsin-like GPCRs in dictyostelids. Some GPCRs have been amplified in the dictyostelids and in specific lineages thereof and through changes in expression patterns may have been repurposed for signalling in multicellular development. Discussion Our phylogenetic analysis suggests that GPCR families 1, 2 and 6 already diverged early in the Amoebozoa, whereas families 3 and 5 expanded later within the dictyostelids. The family 6 cAMP receptors that have experimentally supported roles in multicellular development in dictyostelids ( carA-carD; tasA/B) originated at the root of all dictyostelids and only have weakly associated homologs in Physarum polycephalum. Our analysis identified candidate GPCRs which have evolved in the dictyostelids and could have been co-opted for multicellular development.

Autophagy of the somatic stalk cells likely nurses the propagating spores of Dictyostelid social amoebas

Background:  Autophagy (self-feeding) assists survival of starving cells by partial self-digestion, while dormancy as cysts, spores or seeds enables long-term survival. Starving Dictyostelium amoebas construct multicellular fruiting bodies with spores and stalk cells, with many Dictyostelia still able to encyst individually like their single-celled ancestors. While autophagy mostly occurs in the somatic stalk cells, autophagy gene knock-outs in Dictyostelium discoideum ( D. discoideum) formed no spores and lacked cAMP induction of prespore gene expression. Methods: To investigate whether autophagy also prevents encystation, we knocked-out autophagy genes atg5 and atg7 in the dictyostelid Polysphondylium pallidum, which forms both spores and cysts. We measured spore and cyst differentiation and viability in the knock-out as well as stalk and spore gene expression and its regulation by cAMP. We tested a hypothesis that spores require materials derived from autophagy in stalk cells. Sporulation requires secreted cAMP acting on receptors and intracellular cAMP acting on PKA. We compared the morphology and viability of spores developed in fruiting bodies with spores induced from single cells by stimulation with cAMP and 8Br-cAMP, a membrane-permeant PKA agonist. Results: Loss of autophagy in P. pallidum reduced but did not prevent encystation. Stalk cells still differentiated but stalks were disorganised. However, no spores were formed at all  and cAMP-induced prespore gene expression was lost. D. discoideum spores induced in vitro by cAMP and 8Br-cAMP were smaller and rounder than spores formed multicellularly and while they were not lysed by detergent they germinated not (strain Ax2) or poorly (strain NC4), unlike spores formed in fruiting bodies. Conclusions: The stringent requirement of sporulation on both multicellularity and autophagy, which occurs mostly in stalk cells, suggests that stalk cells nurse the spores through autophagy. This highlights autophagy as a major cause for somatic cell evolution in early multicellularity.

Yellow polyketide pigment suppresses premature hatching in social amoeba

Low-molecular-weight natural products from microbes are indispensable in the development of potent drugs. However, their biological roles within an ecological context often remain elusive. Here, we shed light on natural products from eukaryotic microorganisms that have the ability to transition from single cells to multicellular organisms: the social amoebae. These eukaryotes harbor a large number of polyketide biosynthetic genes in their genomes, yet virtually none of the corresponding products can be isolated or characterized. Using complementary molecular biology approaches, including CRISPR-Cas9, we generated polyketide synthase (pks5) inactivation and overproduction strains of the social amoeba Dictyostelium discoideum. Differential, untargeted metabolomics of wild-type versus mutant fruiting bodies allowed us to pinpoint candidate metabolites derived from the amoebal PKS5. Extrachromosomal expression of the respective gene led to the identification of a yellow polyunsaturated fatty acid. Analysis of the temporospatial production pattern of this compound in conjunction with detailed bioactivity studies revealed the polyketide to be a spore germination suppressor.

Generating polyketide diversity in Dictyostelium: a Steely hybrid polyketide synthase produces alternate products at different developmental stages

The soil is a rich ecosystem where many ecological interactions are mediated by small molecules, and in which amoebae are low-level predators and also prey. The social amoeba Dictyostelium discoideum has a high genomic potential for producing polyketides to mediate its ecological interactions, including the unique 'Steely' enzymes, consisting of a fusion between a fatty acid synthase and a chalcone synthase. We report here that D. discoideum further increases its polyketide potential by using the StlB Steely enzyme, and a downstream chlorinating enzyme, to make both a chlorinated signal molecule, DIF-1, during its multi-cellular development, and a set of abundant polyketides in terminally differentiated stalk cells. We identify one of these as a chlorinated dibenzofuran with potent anti-bacterial activity. To do this, StlB switches expression from prespore to stalk cells in late development and is cleaved to release the chalcone synthase domain. Expression of this domain alone in StlB null cells allows synthesis of the stalk-associated, chlorinated polyketides. Thus, by altered expression and processing of StlB, cells make first a signal molecule, and then abundant secondary metabolites, which we speculate help to protect the mature spores from bacterial infection.

The polyketide synthase StlA is involved in inducing aggregation in Polysphondylium violaceum

In the social amoeba Dictyostelium discoideum, the polyketide MPBD (4-methyl-5-pentylbenzene-1,3-diol) regulates the gene expressions of cAMP signaling to make cells aggregation-competent and also induces spore maturation. The polyketide synthase StlA is responsible for MPBD biosynthesis in D. discoideum and appears to be conserved throughout the major groups of the social amoeba (Dictyostelia). In this study, we analyzed the function of StlA in Polysphondylium violaceum by identifying the gene sequence and creating the knockout mutants. We found that Pv-stlA- mutants had defects only in cell aggregation but not in spore maturation, indicating that the function of StlA in inducing spore maturation is species-specific. We also found that MPBD could rescue the aggregation defect in Pv-stlA- mutants whereas the mutants normally exhibited chemotaxis to their chemoattractant, glorin. Our data suggest that StlA is involved in inducing aggregation in P. violaceum by acting on signaling pathways other than chemotaxis in P. violaceum.

Novel RNAseq-Informed Cell-type Markers and Their Regulation Alter Paradigms of Dictyostelium Developmental Control

Cell differentiation is traditionally monitored with a few marker genes, which may bias results. To understand the evolution and regulation of the spore, stalk, cup and basal disc cells in Dictyostelia, we previously performed RNAseq on purified cell-types of taxon-group representative dictyostelids. Using promoter-lacZ constructs in D. discoideum, we here investigate the spatio-temporal expression pattern of 29 cell-type specific genes. Genes selected for spore- or cup-specificity in RNAseq were validated as such by lacZ expression, but genes selected for stalk-specificity showed variable additional expression in basal disc, early cup or prestalk populations. We measured responses of 25 genes to 15 single or combined regimes of induction by stimuli known to regulate cell differentiation. The outcomes of these experiments were subjected to hierarchical clustering to identify whether common modes of regulation were correlated with specific expression patterns. The analysis identified a cluster combining the spore and cup genes, which shared upregulation by 8-bromo cyclic AMP and down-regulation by Differentiation Inducing Factor 1 (DIF-1). Most stalk-expressed genes combined into a single cluster and shared strong upregulation by cyclic di-guanylate (c-di-GMP), and synergistic upregulation by combined DIF-1 and c-di-GMP. There was no clustering of genes expressed in other soma besides the stalk, but two genes that were only expressed in the stalk did not respond to any stimuli. In contrast to current models, the study indicates the existence of a stem-cell like soma population in slugs, whose members only acquire ultimate cell fate after progressing to their terminal location during fruiting body morphogenesis.

The Cellular and Developmental Roles of Cullins, Neddylation, and the COP9 Signalosome in Dictyostelium discoideum

Cullins (CULs) are a core component of cullin-RING E3 ubiquitin ligases (CRLs), which regulate the degradation, function, and subcellular trafficking of proteins. CULs are post-translationally regulated through neddylation, a process that conjugates the ubiquitin-like modifier protein neural precursor cell expressed developmentally downregulated protein 8 (NEDD8) to target cullins, as well as non-cullin proteins. Counteracting neddylation is the deneddylase, COP9 signalosome (CSN), which removes NEDD8 from target proteins. Recent comparative genomics studies revealed that CRLs and the CSN are highly conserved in Amoebozoa. A well-studied representative of Amoebozoa, the social amoeba Dictyostelium discoideum, has been used for close to 100 years as a model organism for studying conserved cellular and developmental processes owing to its unique life cycle comprised of unicellular and multicellular phases. The organism is also recognized as an exceptional model system for studying cellular processes impacted by human diseases, including but not limited to, cancer and neurodegeneration. Recent work shows that the neddylation inhibitor, MLN4924 (Pevonedistat), inhibits growth and multicellular development in D. discoideum, which supports previous work that revealed the cullin interactome in D. discoideum and the roles of cullins and the CSN in regulating cellular and developmental processes during the D. discoideum life cycle. Here, we review the roles of cullins, neddylation, and the CSN in D. discoideum to guide future work on using this biomedical model system to further explore the evolutionarily conserved functions of cullins and neddylation.

Engineering the amoeba Dictyostelium discoideum for biosynthesis of a cannabinoid precursor and other polyketides

Aromatic polyketides are natural polyphenolic compounds with a broad spectrum of pharmacological activities. Production of those metabolites in the model organisms Escherichia coli and Saccharomyces cerevisiae has been limited by the extensive cellular engineering needed for the coordinated biosynthesis of polyketides and their precursors. In contrast, the amoeba Dictyostelium discoideum is a native producer of secondary metabolites and harbors a wide, but largely unexplored, repertoire of genes for the biosynthesis of polyketides and terpenoids. Here we present D. discoideum as an advantageous chassis for the production of aromatic polyketides. By expressing its native and cognate plant polyketide synthase genes in D. discoideum, we demonstrate production of phlorocaprophenone, methyl-olivetol, resveratrol and olivetolic acid (OA), which is the central intermediate in the biosynthesis of cannabinoids. To facilitate OA synthesis, we further engineered an amoeba/plant inter-kingdom hybrid enzyme that produced OA from primary metabolites in two enzymatic steps, providing a shortcut in a synthetic cannabinoid pathway using the D. discoideum host system.

Interactome and evolutionary conservation of Dictyostelid small GTPases and their direct regulators

GTP binding proteins known as small GTPases make up one of the largest groups of regulatory proteins and control almost all functions of living cells. Their activity is under, respectively, positive and negative regulation by guanine nucleotide exchange factors (GEFs) and GTPase activating proteins (GAPs), which together with their upstream regulators and the downstream targets of the small GTPases form formidable signalling networks. While genomics has revealed the large size of the GTPase, GEF and GAP repertoires, only a small fraction of their interactions and functions have yet been experimentally explored. Dictyostelid social amoebas have been particularly useful in unravelling the roles of many proteins in the Rac-Rho and Ras-Rap families of GTPases in directional cell migration and regulation of the actin cytoskeleton. Genomes and cell-type specific and developmental transcriptomes are available for Dictyostelium species that span the 0.5 billion years of evolution of the group from their unicellular ancestors. In this work, we identified all GTPases, GEFs and GAPs from genomes representative of the four major taxon groups and investigated their phylogenetic relationships and evolutionary conservation and changes in their functional domain architecture and in their developmental and cell-type specific expression. We performed a hierarchical cluster analysis of the expression profiles of the ~2000 analysed genes to identify putative interacting sets of GTPases, GEFs and GAPs, which highlight sets known to interact experimentally and many novel combinations. This work represents a valuable resource for research into all fields of cellular regulation.

Moving the Research Forward: The Best of British Biology Using the Tractable Model System Dictyostelium discoideum

The social amoeba Dictyostelium discoideum provides an excellent model for research across a broad range of disciplines within biology. The organism diverged from the plant, yeast, fungi and animal kingdoms around 1 billion years ago but retains common aspects found in these kingdoms. Dictyostelium has a low level of genetic complexity and provides a range of molecular, cellular, biochemical and developmental biology experimental techniques, enabling multidisciplinary studies to be carried out in a wide range of areas, leading to research breakthroughs. Numerous laboratories within the United Kingdom employ Dictyostelium as their core research model. This review introduces Dictyostelium and then highlights research from several leading British research laboratories, covering their distinct areas of research, the benefits of using the model, and the breakthroughs that have arisen due to the use of Dictyostelium as a tractable model system.

Genome-Wide Analysis of Four Pathotypes of Wheat Rust Pathogen (Puccinia graminis) Reveals Structural Variations and Diversifying Selection

Diseases caused by Puccinia graminis are some of the most devastating diseases of wheat. Extensive genomic understanding of the pathogen has proven helpful not only in understanding host- pathogen interaction but also in finding appropriate control measures. In the present study, whole-genome sequencing of four diverse P. graminis pathotypes was performed to understand the genetic variation and evolution. An average of 63.5 Gb of data per pathotype with about 100× average genomic coverage was achieved with 100-base paired-end sequencing performed with Illumina Hiseq 1000. Genome structural annotations collectively predicted 9273 functional proteins including ~583 extracellular secreted proteins. Approximately 7.4% of the genes showed similarity with the PHI database which is suggestive of their significance in pathogenesis. Genome-wide analysis demonstrated pathotype 117-6 as likely distinct and descended through a different lineage. The 3-6% more SNPs in the regulatory regions and 154 genes under positive selection with their orthologs and under negative selection in the other three pathotypes further supported pathotype 117-6 to be highly diverse in nature. The genomic information generated in the present study could serve as an important source for comparative genomic studies across the genus Puccinia and lead to better rust management in wheat.

Evolution of Multicellular Complexity in The Dictyostelid Social Amoebas

Multicellularity evolved repeatedly in the history of life, but how it unfolded varies greatly between different lineages. Dictyostelid social amoebas offer a good system to study the evolution of multicellular complexity, with a well-resolved phylogeny and molecular genetic tools being available. We compare the life cycles of the Dictyostelids with closely related amoebozoans to show that complex life cycles were already present in the unicellular common ancestor of Dictyostelids. We propose frost resistance as an early driver of multicellular evolution in Dictyostelids and show that the cell signalling pathways for differentiating spore and stalk cells evolved from that for encystation. The stalk cell differentiation program was further modified, possibly through gene duplication, to evolve a new cell type, cup cells, in Group 4 Dictyostelids. Studies in various multicellular organisms, including Dictyostelids, volvocine algae, and metazoans, suggest as a common principle in the evolution of multicellular complexity that unicellular regulatory programs for adapting to environmental change serve as "proto-cell types" for subsequent evolution of multicellular organisms. Later, new cell types could further evolve by duplicating and diversifying the "proto-cell type" gene regulatory networks.

Inhibiting Neddylation with MLN4924 Suppresses Growth and Delays Multicellular Development in Dictyostelium discoideum

Neddylation is a post-translational modification that is essential for a variety of cellular processes and is linked to many human diseases including cancer, neurodegeneration, and autoimmune disorders. Neddylation involves the conjugation of the ubiquitin-like modifier neural precursor cell expressed developmentally downregulated protein 8 (NEDD8) to target proteins, and has been studied extensively in various eukaryotes including fungi, plants, and metazoans. Here, we examine the biological processes influenced by neddylation in the social amoeba, Dictyostelium discoideum, using a well-established inhibitor of neddylation, MLN4924 (pevonedistat). NEDD8, and the target of MLN4924 inhibition, NEDD8-activating enzyme E1 (NAE1), are highly conserved in D. discoideum (Nedd8 and Nae1, respectively). Treatment of D. discoideum cells with MLN4924 increased the amount of free Nedd8, suggesting that MLN4924 inhibited neddylation. During growth, MLN4924 suppressed cell proliferation and folic acid-mediated chemotaxis. During multicellular development, MLN4924 inhibited cyclic adenosine monophosphate (cAMP)-mediated chemotaxis, delayed aggregation, and suppressed fruiting body formation. Together, these findings indicate that neddylation plays an important role in regulating cellular and developmental events during the D. discoideum life cycle and that this organism can be used as a model system to better understand the essential roles of neddylation in eukaryotes, and consequently, its involvement in human disease.

Abundantly expressed class of noncoding RNAs conserved through the multicellular evolution of dictyostelid social amoebas

Aggregative multicellularity has evolved multiple times in diverse groups of eukaryotes, exemplified by the well-studied development of dictyostelid social amoebas, for example, Dictyostelium discoideum However, it is still poorly understood why multicellularity emerged in these amoebas while the majority of other members of Amoebozoa are unicellular. Previously, a novel type of noncoding RNA, Class I RNAs, was identified in D. discoideum and shown to be important for normal multicellular development. Here, we investigated Class I RNA evolution and its connection to multicellular development. We identified a large number of new Class I RNA genes by constructing a covariance model combined with a scoring system based on conserved upstream sequences. Multiple genes were predicted in representatives of each major group of Dictyostelia and expression analysis confirmed that our search approach identifies expressed Class I RNA genes with high accuracy and sensitivity and that the RNAs are developmentally regulated. Further studies showed that Class I RNAs are ubiquitous in Dictyostelia and share highly conserved structure and sequence motifs. In addition, Class I RNA genes appear to be unique to dictyostelid social amoebas because they could not be identified in outgroup genomes, including their closest known relatives. Our results show that Class I RNA is an ancient class of ncRNAs, likely to have been present in the last common ancestor of Dictyostelia dating back at least 600 million years. Based on previous functional analyses and the presented evolutionary investigation, we hypothesize that Class I RNAs were involved in evolution of multicellularity in Dictyostelia.

Cancer and the breakdown of multicellularity: What Dictyostelium discoideum, a social amoeba, can teach us

Ancient pathways promoting unicellularity and multicellularity are associated with cancer, the former being pro-oncogenic and the latter acting to suppress oncogenesis. However, there are only a limited number of non-vertebrate models for studying these pathways. Here, we review Dictyostelium discoideum and describe how it can be used to understand these gene networks. D. discoideum has a unicellular and multicellular life cycle, making it possible to study orthologs of cancer-associated genes in both phases. During development, differentiated amoebae form a fruiting body composed of a mass of spores that are supported atop a stalk. A portion of the cells sacrifice themselves to become non-reproductive stalk cells. Cheating disrupts the principles of multicellularity, as cheater cells alter their cell fate to preferentially become spores. Importantly, D. discoideum has gene networks and several strategies for maintaining multicellularity. Therefore, D. discoideum can help us better understand how conserved genes and pathways involved in multicellularity also influence cancer development, potentially identifying new therapeutic avenues.

Unusual Occurrence of Two Bona-Fide CCA-Adding Enzymes in Dictyostelium discoideum

Dictyostelium discoideum, the model organism for the evolutionary supergroup of Amoebozoa, is a social amoeba that, upon starvation, undergoes transition from a unicellular to a multicellular organism. In its genome, we identified two genes encoding for tRNA nucleotidyltransferases. Such pairs of tRNA nucleotidyltransferases usually represent collaborating partial activities catalyzing CC- and A-addition to the tRNA 3'-end, respectively. In D. discoideum, however, both enzymes exhibit identical activities, representing bona-fide CCA-adding enzymes. Detailed characterization of the corresponding activities revealed that both enzymes seem to be essential and are regulated inversely during different developmental stages of D. discoideum. Intriguingly, this is the first description of two functionally equivalent CCA-adding enzymes using the same set of tRNAs and showing a similar distribution within the cell. This situation seems to be a common feature in Dictyostelia, as other members of this phylum carry similar pairs of tRNA nucleotidyltransferase genes in their genome.

Evolution of multicellularity in Dictyostelia

The well-orchestrated multicellular life cycle of Dictyostelium discoideum has fascinated biologists for over a century. Self-organisation of its amoebas into aggregates, migrating slugs and fruiting structures by pulsatile cAMP signalling and their ability to follow separate differentiation pathways in well-regulated proportions continue to be topics under investigation. A striking aspect of D. discoideum development is the recurrent use of cAMP as chemoattractant, differentiation inducing signal and second messenger for other signals that control the developmental programme. D. discoideum is one of >150 species of Dictyostelia and aggregative life styles similar to those of Dictyostelia evolved many times in eukaryotes. Here we review experimental studies investigating how phenotypic complexity and cAMP signalling co-evolved in Dictyostelia. In addition, we summarize comparative genomic studies of multicellular Dictyostelia and unicellular Amoebozoa aimed to identify evolutionary conservation and change in all genes known to be essential for D. discoideum development.

Cold climate adaptation is a plausible cause for evolution of multicellular sporulation in Dictyostelia

Unicellular protozoa that encyst individually upon starvation evolved at least eight times into organisms that instead form multicellular fruiting bodies with spores. The Dictyostelia are the largest and most complex group of such organisms. They can be subdivided into 4 major groups, with many species in groups 1–3 having additionally retained encystment. To understand fitness differences between spores and cysts, we measured long-term survival of spores and cysts under climate-mimicking conditions, investigated spore and cyst ultrastructure, and related fitness characteristics to species ecology. We found that spores and cysts survived 22 °C equally well, but that spores survived wet and dry frost better than cysts, with group 4 spores being most resilient. Spore walls consist of three layers and those of cysts of maximally two, while spores were also more compacted than cysts, with group 4 spores being the most compacted. Group 4 species were frequently isolated from arctic and alpine zones, which was rarely the case for group 1–3 species. We inferred a fossil-calibrated phylogeny of Dictyostelia, which showed that its two major branches diverged 0.52 billion years ago, following several global glaciations. Our results suggest that Dictyostelium multicellular sporulation was a likely adaptation to a cold climate.

Loss of the Polyketide Synthase StlB Results in Stalk Cell Overproduction in Polysphondylium violaceum

Major phenotypic innovations in social amoeba evolution occurred at the transition between the Polysphondylia and group 4 Dictyostelia, which comprise the model organism Dictyostelium discoideum, such as the formation of a new structure, the basal disk. Basal disk differentiation and robust stalk formation require the morphogen DIF-1, synthesized by the polyketide synthase StlB, the des-methyl-DIF-1 methyltransferase DmtA, and the chlorinase ChlA, which are conserved throughout Dictyostelia. To understand how the basal disk and other innovations evolved in group 4, we sequenced and annotated the Polysphondylium violaceum (Pvio) genome, performed cell type-specific transcriptomics to identify cell-type marker genes, and developed transformation and gene knock-out procedures for Pvio. We used the novel methods to delete the Pvio stlB gene. The Pvio stlB- mutants formed misshapen curly sorogens with thick and irregular stalks. As fruiting body formation continued, the upper stalks became more regular, but structures contained 40% less spores. The stlB- sorogens overexpressed a stalk gene and underexpressed a (pre)spore gene. Normal fruiting body formation and sporulation were restored in Pvio stlB- by including DIF-1 in the supporting agar. These data indicate that, although conserved, stlB and its product(s) acquired both a novel role in the group 4 Dictyostelia and a role opposite to that in its sister group.

Characterization of Mollivirus kamchatka, the First Modern Representative of the Proposed Molliviridae Family of Giant Viruses

Microbes trapped in permanently frozen paleosoils (permafrost) are the focus of increasing research in the context of global warming. Our previous investigations led to the discovery and reactivation of two Acanthamoeba-infecting giant viruses, Mollivirus sibericum and Pithovirus sibericum, from a 30,000-year old permafrost layer. While several modern pithovirus strains have since been isolated, no contemporary mollivirus relative was found. We now describe Mollivirus kamchatka, a close relative to M. sibericum, isolated from surface soil sampled on the bank of the Kronotsky River in Kamchatka, Russian Federation. This discovery confirms that molliviruses have not gone extinct and are at least present in a distant subarctic continental location. This modern isolate exhibits a nucleocytoplasmic replication cycle identical to that of M. sibericum Its spherical particle (0.6 μm in diameter) encloses a 648-kb GC-rich double-stranded DNA genome coding for 480 proteins, of which 61% are unique to these two molliviruses. The 461 homologous proteins are highly conserved (92% identical residues, on average), despite the presumed stasis of M. sibericum for the last 30,000 years. Selection pressure analyses show that most of these proteins contribute to virus fitness. The comparison of these first two molliviruses clarify their evolutionary relationship with the pandoraviruses, supporting their provisional classification in a distinct family, the Molliviridae, pending the eventual discovery of intermediary missing links better demonstrating their common ancestry.IMPORTANCE Virology has long been viewed through the prism of human, cattle, or plant diseases, leading to a largely incomplete picture of the viral world. The serendipitous discovery of the first giant virus visible under a light microscope (i.e., >0.3 μm in diameter), mimivirus, opened a new era of environmental virology, now incorporating protozoan-infecting viruses. Planet-wide isolation studies and metagenome analyses have shown the presence of giant viruses in most terrestrial and aquatic environments, including upper Pleistocene frozen soils. Those systematic surveys have led authors to propose several new distinct families, including the Mimiviridae, Marseilleviridae, Faustoviridae, Pandoraviridae, and Pithoviridae We now propose to introduce one additional family, the Molliviridae, following the description of M. kamchatka, the first modern relative of M. sibericum, previously isolated from 30,000-year-old arctic permafrost.

Transcriptional Responses of Dictyostelium discoideum Exposed to Different Classes of Bacteria

Dictyostelium discoideum amoebae feed by ingesting bacteria, then killing them in phagosomes. Ingestion and killing of different bacteria have been shown to rely on largely different molecular mechanisms. One would thus expect that D. discoideum adapts its ingestion and killing machinery when encountering different bacteria. In this study, we investigated by RNA sequencing if and how D. discoideum amoebae respond to the presence of different bacteria by modifying their gene expression patterns. Each bacterial species analyzed induced a specific modification of the transcriptome. Bacteria such as Bacillus subtilis, Klebsiella pneumoniae, or Mycobacterium marinum induced a specific and different transcriptional response, while Micrococcus luteus did not trigger a significant gene regulation. Although folate has been proposed to be one of the key molecules secreted by bacteria and recognized by hunting amoebae, it elicited a very specific and restricted transcriptional signature, distinct from that triggered by any bacteria analyzed here. Our results indicate that D. discoideum amoebae respond in a highly specific, almost non-overlapping manner to different species of bacteria. We additionally identify specific sets of genes that can be used as reporters of the response of D. discoideum to different bacteria.

Phylogeny-wide conservation and change in developmental expression, cell-type specificity and functional domains of the transcriptional regulators of social amoebas

Background

Dictyostelid social amoebas self-organize into fruiting bodies, consisting of spores and up to four supporting cell types in the phenotypically most complex taxon group 4. High quality genomes and stage- and cell-type specific transcriptomes are available for representative species of each of the four taxon groups. To understand how evolution of gene regulation in Dictyostelia contributed to evolution of phenotypic complexity, we analysed conservation and change in abundance, functional domain architecture and developmental regulation of their transcription factors (TFs).

Results

We detected 440 sequence-specific TFs across 33 families, of which 68% were upregulated in multicellular development and about half conserved throughout Dictyostelia. Prespore cells expressed two times more TFs than prestalk cells, but stalk cells expressed more TFs than spores, suggesting that gene expression events that define spores occur earlier than those that define stalk cells. Changes in TF developmental expression, but not in TF abundance or functional domains occurred more frequently between group 4 and groups 1–3, than between the more distant branches formed by groups 1 + 2 and 3 + 4.

Conclusions

Phenotypic innovation is correlated with changes in TF regulation, rather than functional domain- or TF acquisition. The function of only 34 TFs is known. Of 12 TFs essential for cell differentiation, 9 are expressed in the cell type for which they are required. The information acquired here on conserved cell type specifity of 120 additional TFs can effectively guide further functional analysis, while observed evolutionary change in TF developmental expression may highlight how genotypic change caused phenotypic innovation.

Structure of Dictyostelium discoideum telomeres. Analysis of possible replication mechanisms

Telomeres are nucleo-protein structures that protect the ends of eukaryotic chromosomes. They are not completely synthesized during DNA replication and are elongated by specific mechanisms. The structure of the telomeres and the elongation mechanism have not been determined in Dictyostelium discoideum. This organism presents extrachromosomal palindromic elements containing two copies of the rDNA, also present at the end of the chromosomes. In this article the structure of the terminal region of the rDNA is shown to consist of repetitions of the A(G)n sequence where the number of Gs is variable. These repeats extend as a 3’ single stranded region. The G-rich region is preceded by four tandem repetitions of two different DNA motifs. D. discoideum telomere reverse transcriptase homologous protein (TERTHP) presented RNase-sensitive enzymatic activity and was required to maintain telomere structure since terthp-mutant strains presented reorganizations of the DNA terminal regions. These modifications were different in several terthp-mutants and changed with their prolonged culture and subcloning. However, the terthp gene is not essential for D. discoideum proliferation. Telomeres could be maintained in terthp-mutant strains by homologous recombination mechanisms such as ALT (Alternative Lengthening of Telomeres) or HAATI (heterochromatin amplification-mediated and telomerase-independent). In agreement with this hypothesis, the expression of mRNAs coding for several proteins involved in homologous recombination was induced in terthp-mutant strains. Extrachromosomal rDNA could serve as substrate in these DNA homologous recombination reactions.

A telomerase with novel non-canonical roles: TERT controls cellular aggregation and tissue size in Dictyostelium

Telomerase, particularly its main subunit, the reverse transcriptase, TERT, prevents DNA erosion during eukaryotic chromosomal replication, but also has poorly understood non-canonical functions. Here, in the model social amoeba Dictyostelium discoideum, we show that the protein encoded by tert has telomerase-like motifs, and regulates, non-canonically, important developmental processes. Expression levels of wild-type (WT) tert were biphasic, peaking at 8 and 12 h post-starvation, aligning with developmental events, such as the initiation of streaming (~7 h) and mound formation (~10 h). In tert KO mutants, however, aggregation was delayed until 16 h. Large, irregular streams formed, then broke up, forming small mounds. The mound-size defect was not induced when a KO mutant of countin (a master size-regulating gene) was treated with TERT inhibitors, but anti-countin antibodies did rescue size in the tert KO. Although, conditioned medium (CM) from countin mutants failed to rescue size in the tert KO, tert KO CM rescued the countin KO phenotype. These and additional observations indicate that TERT acts upstream of smlA/countin: (i) the observed expression levels of smlA and countin, being respectively lower and higher (than WT) in the tert KO; (ii) the levels of known size-regulation intermediates, glucose (low) and adenosine (high), in the tert mutant, and the size defect's rescue by supplemented glucose or the adenosine-antagonist, caffeine; (iii) the induction of the size defect in the WT by tert KO CM and TERT inhibitors. The tert KO's other defects (delayed aggregation, irregular streaming) were associated with changes to cAMP-regulated processes (e.g. chemotaxis, cAMP pulsing) and their regulatory factors (e.g. cAMP; acaA, carA expression). Overexpression of WT tert in the tert KO rescued these defects (and size), and restored a single cAMP signaling centre. Our results indicate that TERT acts in novel, non-canonical and upstream ways, regulating key developmental events in Dictyostelium.

Discrimination Experiments in Entamoeba and Evidence from Other Protists Suggest Pathogenic Amebas Cooperate with Kin to Colonize Hosts and Deter Rivals

Entamoeba histolytica is one of the least understood protists in terms of taxa, clone, and kin discrimination/recognition ability. However, the capacity to tell apart same or self (clone/kin) from different or nonself (nonclone/nonkin) has long been demonstrated in pathogenic eukaryotes like Trypanosoma and Plasmodium, free-living social amebas (Dictyostelium, Polysphondylium), budding yeast (Saccharomyces), and in numerous bacteria and archaea (prokaryotes). Kin discrimination/recognition is explained under inclusive fitness theory; that is, the reproductive advantage that genetically closely related organisms (kin) can gain by cooperating preferably with one another (rather than with distantly related or unrelated individuals), minimizing antagonism and competition with kin, and excluding genetic strangers (or cheaters = noncooperators that benefit from others' investments in altruistic cooperation). In this review, we rely on the outcomes of in vitro pairwise discrimination/recognition encounters between seven Entamoeba lineages to discuss the biological significance of taxa, clone, and kin discrimination/recognition in a range of generalist and specialist species (close or distantly related phylogenetically). We then focus our discussion on the importance of these laboratory observations for E. histolytica's life cycle, host infestation, and implications of these features of the amebas' natural history for human health (including mitigation of amebiasis).

The putative G protein-coupled receptor GrlD mediates extracellular polyphosphate sensing in Dictyostelium discoideum

Five or more orthophosphates bound together by high-energy phosphoanhydride bonds are highly ubiquitous inorganic molecules called polyphosphate. Polyphosphate acts as a signaling molecule eliciting a number of responses in eukaryotic cells, but the mechanisms mediating these effects are poorly understood. Proliferating Dictyostelium discoideum cells accumulate extracellular polyphosphate. At extracellular concentrations similar to those observed in stationary phase cells, polyphosphate inhibits proteasome activity and proliferation, and induces aggregation. Here we identify GrlD as a putative G protein–coupled receptor that mediates binding of extracellular polyphosphate to the cell surface. Cells lacking GrlD do not respond to polyphosphate-induced proteasome inhibition, aggregation, or proliferation inhibition. Polyphosphate also elicits differential effects on cell-substratum adhesion and cytoskeletal F-actin levels based on nutrient availability, and these effects were also mediated by GrlD. Starving cells also accumulate extracellular polyphosphate. Starved cells treated with exopolyphosphatase failed to aggregate effectively, suggesting that polyphosphate also acts as a signaling molecule during starvation-induced development of Dictyostelium. Together, these results suggest that a eukaryotic cell uses a G protein–coupled receptor to mediate the sensing and response to extracellular polyphosphate.

A well supported multi gene phylogeny of 52 dictyostelia

The Dictyostelid social amoebas are a popular model system for cell- and developmental biology and for evolution of sociality. Small subunit (SSU) ribosomal DNA-based phylogenies subdivide the known 150 species into four major and some minor groups, but lack resolution within groups, particularly group 4, and, as shown by genome-based phylogenies of 11 species, showed errors in the position of the root and nodes separating major clades. We are interested in the evolution of cell-type specialization, which particularly expanded in group 4. To construct a more robust phylogeny, we first included 7 recently sequenced genomes in the genome-based phylogeny of 47 functionally divergent proteins and next selected 6 proteins (Agl, AmdA, PurD, PurL, RpaA, SmdA) that independently or in sets of two fully reproduced the core-phylogeny. We amplified their coding regions from 34 Dictyostelium species and combined their concatenated sequences with those identified in the 18 genomes to generate a fully resolved phylogeny. The new AAPPRS based phylogeny (after the acronym of the 6 proteins) subdivides group 4 into 2 branches. These branches further resolve into 5 clades, rather than the progressively nested group 4 topology of the SSU rDNA tree, and also re-orders taxa in the other major groups. Ancestral state reconstruction of 25 phenotypic traits returned higher "goodness of fit" metrics for evolution of 19 of those traits over the AAPPRS tree, than over the SSU rDNA tree. The novel tree provides a solid framework for studying the evolution of cell-type specialization, signalling and other cellular processes in particularly group 4, which contains the model Dictyostelid D. discoideum.

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.

Diversity and Functional Evolution of Terpene Synthases in Dictyostelid Social Amoebae

Dictyostelids, or social amoebae, have a unique life style in forming multicellular fruiting bodies from unicellular amoeboids upon starvation. Recently, dictyostelids were found to contain terpene synthase (TPS) genes, a gene type of secondary metabolism previously known to occur only in plants, fungi and bacteria. Here we report an evolutionary functional study of dictyostelid TPS genes. The number of TPS genes in six species of dictyostelids examined ranges from 1 to 19; and the model species Dictyostelium purpureum contains 12 genes. Using in vitro enzyme assays, the 12 TPS genes from D. purpureum were shown to encode functional enzymes with distinct product profiles. The expression of the 12 TPS genes in D. purpureum is developmentally regulated. During multicellular development, D. purpureum releases a mixture of volatile terpenes dominated by sesquiterpenes that are the in vitro products of a subset of the 12 TPS genes. The quality and quantity of the terpenes released from D. purpureum, however, bear little resemblance to those of D. discoideum, a closely related dictyostelid. Despite these variations, the conserved clade of dictyostelid TPSs, which have an evolutionary distance of more than 600 million years, has the same biochemical function, catalyzing the formation of a sesquiterpene protoillud-7-ene. Taken together, our results indicate that the dynamic evolution of dictyostelid TPS genes includes both purifying selection of an orthologous group and species-specific expansion with functional divergence. Consequently, the terpenes produced by these TPSs most likely have conserved as well as species-adaptive biological functions as chemical languages in dictyostelids.

The Chemoattractant Glorin Is Inactivated by Ester Cleavage during Early Multicellular Development of Polysphondylium pallidum

Among the amoebozoan species capable of forming fruiting bodies, the dictyostelid social amoebae stand out since they form true multicellular organisms by means of single cell aggregation. Upon food depletion, cells migrate across gradients of extracellular signals initiated by cells in aggregation centers. The model species that is widely used to study multicellular development of social amoebae, Dictyostelium discoideum, uses cyclic adenosine monophosphate (cAMP) as a chemoattractant to coordinate aggregation. Molecular phylogeny studies suggested that social amoebae evolved in four major groups, of which groups 1 and 2 are paraphyletic to groups 3 and 4. During early development, intercellular communication with cAMP appears to be restricted to group 4 species. Cells of group 1 and 2 taxa do not respond chemotactically to extracellular cAMP and likely use a dipeptide chemoattractant known as glorin ( N-propionyl-γ-L-glutamyl-L-ornithin-δ-lactam-ethylester) to regulate aggregation. Directional migration of glorin-responsive cells requires the periodic breakdown of the chemoattractant. Here, we identified an extracellular enzymatic activity (glorinase) in the glorin-responsive group 2 taxon Polysphondylium pallidum leading to the inactivation of glorin. We determined the inactivation mechanism to proceed via hydrolytic ethyl ester cleavage of the γ-glutamyl moiety of glorin. Synthetic glorinamide, in which the ethyl ester group was substituted by an ethyl amide group, had glorin-like biological activity but was resistant to degradation by glorinase. Our observations pave the way for future investigations toward an ancient eukaryotic chemotaxis system.

The transcription factor Spores Absent A is a PKA dependent inducer of Dictyostelium sporulation

Sporulation in Dictyostelium fruiting bodies evolved from amoebozoan encystation with both being induced by cAMP acting on PKA, but with downstream components still being unknown. Using tagged mutagenesis to find missing pathway components, we identified a sporeless mutant defective in a nuclear protein, SpaA. Expression of prespore genes was strongly reduced in spaA- cells, while expression of many spore stage genes was absent. Chromatin immunoprecipitation (ChIP) of a SpaA-YFP gene fusion showed that (pre)spore gene promoters bind directly to SpaA, identifying SpaA as a transcriptional regulator. SpaA dependent spore gene expression required PKA in vivo and was stimulated in vitro by the membrane-permeant PKA agonist 8Br-cAMP. The PKA agonist also promoted SpaA binding to (pre)spore promoters, placing SpaA downstream of PKA. Sequencing of SpaA-YFP ChIPed DNA fragments revealed that SpaA binds at least 117 (pre)spore promoters, including those of other transcription factors that activate some spore genes. These factors are not in turn required for spaA expression, identifying SpaA as the major trancriptional inducer of sporulation.

An Autocrine Proliferation Repressor Regulates Dictyostelium discoideum Proliferation and Chemorepulsion Using the G Protein-Coupled Receptor GrlH

In eukaryotic microbes, little is known about signals that inhibit the proliferation of the cells that secrete the signal, and little is known about signals (chemorepellents) that cause cells to move away from the source of the signal. Autocrine proliferation repressor protein A (AprA) is a protein secreted by the eukaryotic microbe Dictyostelium discoideum. AprA is a chemorepellent for and inhibits the proliferation of D. discoideum. We previously found that cells sense AprA using G proteins, suggesting the existence of a G protein-coupled AprA receptor. To identify the AprA receptor, we screened mutants lacking putative G protein-coupled receptors. We found that, compared to the wild-type strain, cells lacking putative receptor GrlH (grlH¯ cells) show rapid proliferation, do not have large numbers of cells moving away from the edges of colonies, are insensitive to AprA-induced proliferation inhibition and chemorepulsion, and have decreased AprA binding. Expression of GrlH in grlH¯ cells (grlH¯/grlH^OE) rescues the phenotypes described above. These data indicate that AprA signaling may be mediated by GrlH in D. discoideum.

Little is known about how eukaryotic cells can count themselves and thus regulate the size of a tissue or density of cells. In addition, little is known about how eukaryotic cells can sense a repellant signal and move away from the source of the repellant, for instance, to organize the movement of cells in a developing embryo or to move immune cells out of a tissue. In this study, we found that a eukaryotic microbe uses G protein-coupled receptors to mediate both cell density sensing and chemorepulsion.

Multiple Roots of Fruiting Body Formation in Amoebozoa

Establishment of multicellularity represents a major transition in eukaryote evolution. A subgroup of Amoebozoa, the dictyosteliids, has evolved a relatively simple aggregative multicellular stage resulting in a fruiting body supported by a stalk. Protosteloid amoeba, which are scattered throughout the amoebozoan tree, differ by producing only one or few single stalked spores. Thus, one obvious difference in the developmental cycle of protosteliids and dictyosteliids seems to be the establishment of multicellularity. To separate spore development from multicellular interactions, we compared the genome and transcriptome of a Protostelium species (Protostelium aurantium var. fungivorum) with those of social and solitary members of the Amoebozoa. During fruiting body formation nearly 4,000 genes, corresponding to specific pathways required for differentiation processes, are upregulated. A comparison with genes involved in the development of dictyosteliids revealed conservation of >500 genes, but most of them are also present in Acanthamoeba castellanii for which fruiting bodies have not been documented. Moreover, expression regulation of those genes differs between P. aurantium and Dictyostelium discoideum. Within Amoebozoa differentiation to fruiting bodies is common, but our current genome analysis suggests that protosteliids and dictyosteliids used different routes to achieve this. Most remarkable is both the large repertoire and diversity between species in genes that mediate environmental sensing and signal processing. This likely reflects an immense adaptability of the single cell stage to varying environmental conditions. We surmise that this signaling repertoire provided sufficient building blocks to accommodate the relatively simple demands for cell-cell communication in the early multicellular forms.

A Deep Hidden Diversity of Dictyostelia

Dictyostelia is a monophyletic group of transiently multicellular (sorocarpic) amoebae, whose study is currently limited to laboratory culture. This tends to favour faster growing species with robust sorocarps, while species with smaller more delicate sorocarps constitute most of the group's taxonomic breadth. The number of known species is also small (∼150) given Dictyostelia's molecular depth and apparent antiquity (>600 myr). Nonetheless, dictyostelid sequences are rarely recovered in culture independent sampling (ciPCR) surveys. We developed ciPCR primers to specifically target dictyostelid small subunit (SSU or 18S) rDNA and tested them on total DNAs extracted from a wide range of soils from five continents. The resulting clone libraries show mostly dictyostelid sequences (∼90%), and phylogenetic analyses of these sequences indicate novel lineages in all four dictyostelid families and most genera. This is especially true for the species-rich Heterostelium and Dictyosteliaceae but also the less species-rich Raperosteliaceae. However, the most novel deep branches are found in two very species-poor taxa, including the deepest branch yet seen in the highly divergent Cavenderiaceae. These results confirm a deep hidden diversity of Dictyostelia, potentially including novel morphologies and developmental schemes. The primers and protocols presented here should also enable more comprehensive studies of dictyostelid ecology.

The emerging picture of the mitochondrial protein import complexes of Amoebozoa supergroup

BACKGROUND

The existence of mitochondria-related organelles (MROs) is proposed for eukaryotic organisms. The Amoebozoa includes some organisms that are known to have mitosomes but also organisms that have aerobic mitochondria. However, the mitochondrial protein apparatus of this supergroup remains largely unsampled, except for the mitochondrial outer membrane import complexes studied recently. Therefore, in this study we investigated the mitochondrial inner membrane and intermembrane space complexes, using the available genome and transcriptome sequences.

RESULTS

When compared with the canonical cognate complexes described for the yeast Saccharomyces cerevisiae, amoebozoans with aerobic mitochondria, display lower differences in the number of subunits predicted for these complexes than the mitochondrial outer membrane complexes, although the predicted subunits appear to display different levels of diversity in regard to phylogenetic position and isoform numbers. For the putative mitosome-bearing amoebozoans, the number of predicted subunits suggests the complex elimination distinctly more pronounced than in the case of the outer membrane ones.

CONCLUSION

The results concern the problem of mitochondrial and mitosome protein import machinery structural variability and the reduction of their complexity within the currently defined supergroup of Amoebozoa. This results are crucial for better understanding of the Amoebozoa taxa of both biomedical and evolutionary importance.

ProtozoaDB 2.0: A Trypanosoma Brucei Case Study

Over the last decade new species of Protozoa have been sequenced and deposited in GenBank. Analyzing large amounts of genomic data, especially using Next Generation Sequencing (NGS), is not a trivial task, considering that researchers used to deal or focus their studies on few genes or gene families or even small genomes. To facilitate the information extraction process from genomic data, we developed a database system called ProtozoaDB that included five genomes of Protozoa in its first version. In the present study, we present a new version of ProtozoaDB called ProtozoaDB 2.0, now with the genomes of 22 pathogenic Protozoa. The system has been fully remodeled to allow for new tools and a more expanded view of data, and now includes a number of analyses such as: (i) similarities with other databases (model organisms, the Conserved Domains Database, and the Protein Data Bank); (ii) visualization of KEGG metabolic pathways; (iii) the protein structure from PDB; (iv) homology inferences; (v) the search for related publications in PubMed; (vi) superfamily classification; and (vii) phenotype inferences based on comparisons with model organisms. ProtozoaDB 2.0 supports RESTful Web Services to make data access easier. Those services were written in Ruby language using Ruby on Rails (RoR). This new version also allows a more detailed analysis of the object of study, as well as expanding the number of genomes and proteomes available to the scientific community. In our case study, a group of prenyltransferase proteinsalready described in the literature was found to be a good drug target for Trypanosomatids.

A tripeptidyl peptidase 1 is a binding partner of the Golgi pH regulator (GPHR) in Dictyostelium

Mutations in tripeptidyl peptidase 1 (TPP1) have been associated with late infantile neuronal ceroid lipofuscinosis (NCL), a neurodegenerative disorder. TPP1 is a lysosomal serine protease, which removes tripeptides from the N-terminus of proteins and is composed of an N-terminal prodomain and a catalytic domain. It is conserved in mammals, amphibians, fish and the amoeba Dictyostelium discoideum. D. discoideum harbors at least six genes encoding TPP1, tpp1A to tpp1F. We identified TPP1F as binding partner of Dictyostelium GPHR (Golgi pH regulator), which is an evolutionarily highly conserved intracellular transmembrane protein. A region encompassing the DUF3735 (GPHR_N) domain of GPHR was responsible for the interaction. In TPP1F, the binding site is located in the prodomain of the protein. The tpp1F gene is transcribed throughout development and translated into a polypeptide of ∼65 kDa. TPP1 activity was demonstrated for TPP1F-GFP immunoprecipitated from D. discoideum cells. Its activity could be inhibited by addition of the recombinant DUF3735 domain of GPHR. Knockout tpp1F mutants did not display any particular phenotype, and TPP1 activity was not abrogated, presumably because tpp1B compensates as it has the highest expression level of all the TPP1 genes during growth. The GPHR interaction was not restricted to TPP1F but occurred also with TPP1B. As previous reports show that the majority of the TPP1 mutations in NCL resulted in reduction or loss of enzyme activity, we suggest that Dicyostelium could be used as a model system in which to test new reagents that could affect the activity of the protein and ameliorate the disease.

Summary: Interaction of Dictyostelium tripeptidyl peptidase 1 with GPHR could be relevant for studies of the human enzyme, which is associated with a neurodegenerative disorder.

Multiple Dictyostelid Species Destroy Biofilms of Klebsiella oxytoca and Other Gram Negative Species

Dictyostelids are free-living phagocytes that feed on bacteria in diverse habitats. When bacterial prey is in short supply or depleted, they undergo multicellular development culminating in the formation of dormant spores. In this work, we tested isolates representing four dictyostelid species from two genera (Dictyostelium and Polysphondylium) for the potential to feed on biofilms preformed on glass and polycarbonate surfaces. The abilities of dictyostelids were monitored for three hallmarks of activity: 1) spore germination on biofilms, 2) predation on biofilm enmeshed bacteria by phagocytic cells and 3) characteristic stages of multicellular development (streaming and fructification). We found that all dictyostelid isolates tested could feed on biofilm enmeshed bacteria produced by human and plant pathogens: Klebsiella oxytoca, Pseudomonas aeruginosa, Pseudomonas syringae, Erwinia amylovora 1189 (biofilm former) and E. amylovora 1189 Δams (biofilm deficient mutant). However, when dictyostelids were fed planktonic E. amylovora Δams the bacterial cells exhibited an increased susceptibility to predation by one of the two dictyostelid strains they were tested against. Taken together, the qualitative and quantitative data presented here suggest that dictyostelids have preferences in bacterial prey which affects their efficiency of feeding on bacterial biofilms.

A polycystin-type transient receptor potential (Trp) channel that is activated by ATP

ATP and ADP are ancient extra-cellular signalling molecules that in Dictyostelium amoebae cause rapid, transient increases in cytosolic calcium due to an influx through the plasma membrane. This response is independent of hetero-trimeric G-proteins, the putative IP3 receptor IplA and all P2X channels. We show, unexpectedly, that it is abolished in mutants of the polycystin-type transient receptor potential channel, TrpP. Responses to the chemoattractants cyclic-AMP and folic acid are unaffected in TrpP mutants. We report that the DIF morphogens, cyclic-di-GMP, GABA, glutamate and adenosine all induce strong cytoplasmic calcium responses, likewise independently of TrpP. Thus, TrpP is dedicated to purinergic signalling. ATP treatment causes cell blebbing within seconds but this does not require TrpP, implicating a separate purinergic receptor. We could detect no effect of ATP on chemotaxis and TrpP mutants grow, chemotax and develop almost normally in standard conditions. No gating ligand is known for the human homologue of TrpP, polycystin-2, which causes polycystic kidney disease. Our results now show that TrpP mediates purinergic signalling in Dictyostelium and is directly or indirectly gated by ATP.

Summary: We show that a Trp channel related to the mammalian polycystin channel, rather than a P2X receptor, is responsible for the purinergic stimulation of cytosolic calcium levels in Dictyostelium cells.

Improved annotation with de novo transcriptome assembly in four social amoeba species

Background

Annotation of gene models and transcripts is a fundamental step in genome sequencing projects. Often this is performed with automated prediction pipelines, which can miss complex and atypical genes or transcripts. RNA sequencing (RNA-seq) data can aid the annotation with empirical data. Here we present de novo transcriptome assemblies generated from RNA-seq data in four Dictyostelid species: D. discoideum, P. pallidum, D. fasciculatum and D. lacteum. The assemblies were incorporated with existing gene models to determine corrections and improvement on a whole-genome scale. This is the first time this has been performed in these eukaryotic species.

Results

An initial de novo transcriptome assembly was generated by Trinity for each species and then refined with Program to Assemble Spliced Alignments (PASA). The completeness and quality were assessed with the Benchmarking Universal Single-Copy Orthologs (BUSCO) and Transrate tools at each stage of the assemblies. The final datasets of 11,315-12,849 transcripts contained 5,610-7,712 updates and corrections to >50% of existing gene models including changes to hundreds or thousands of protein products. Putative novel genes are also identified and alternative splice isoforms were observed for the first time in P. pallidum, D. lacteum and D. fasciculatum.

Conclusions

In taking a whole transcriptome approach to genome annotation with empirical data we have been able to enrich the annotations of four existing genome sequencing projects. In doing so we have identified updates to the majority of the gene annotations across all four species under study and found putative novel genes and transcripts which could be worthy for follow-up. The new transcriptome data we present here will be a valuable resource for genome curators in the Dictyostelia and we propose this effective methodology for use in other genome annotation projects.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-017-3505-0) contains supplementary material, which is available to authorized users.