Comparative genomics of the social amoebae Dictyostelium discoideum and Dictyostelium purpureum

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.

Chromosome-level genome assembly and annotation of the social amoeba Dictyostelium firmibasis

Dicytostelium firmibasis is a member of Dictyostelia, a group of social amoebae that upon starvation display aggregative multicellularity where the amoebae transition from uni- to multicellular life. The D. firmibasis genome assembly that is currently available is of limited use due to its low contiguity, large number of undetermined bases, and lack of annotations. Here we used Nanopore long read sequencing, complemented with Illumina sequencing, and developmental transcriptomics as well as small RNA-sequencing, to present a new, fully annotated, chromosome-level D. firmibasis genome assembly. The new assembly contains no undetermined bases, and consists mainly of six large contigs representing the chromosomes, as well as a complete mitochondrial genome. This new genome assembly will be a valuable tool, allowing comprehensive comparison to Dictyostelium discoideum, the dictyostelid genetically tractable model. Further, the new genome will be important for studies of evolutionary processes governing the transition from unicellular to multicellular organisms and will aid in the sequencing and annotation of other dictyostelids genomes, many of which are currently of poor quality.

Isolation and Structure Determination of New Pyrones from Dictyostelium spp. Cellular Slime Molds Coincubated with Pseudomonas spp

Cellular slime molds are excellent model organisms in the field of cell and developmental biology because of their simple developmental patterns. During our studies on the identification of bioactive molecules from secondary metabolites of cellular slime molds toward the development of novel pharmaceuticals, we revealed the structural diversity of secondary metabolites. Cellular slime molds grow by feeding on bacteria, such as Klebsiella aerogenes and Escherichia coli, without using medium components. Although changing the feeding bacteria is expected to affect dramatically the secondary metabolite production, the effect of the feeding bacteria on the production of secondary metabolites is not known. Herein, we report the isolation and structure elucidation of clavapyrone (1) from Dictyostelium clavatum, intermedipyrone (2) from D. magnum, and magnumiol (3) from D. intermedium. These compounds are not obtained from usual cultural conditions with Klebsiella aerogenes but obtained from coincubated conditions with Pseudomonas spp. The results demonstrate the diversity of the secondary metabolites of cellular slime molds and suggest that widening the range of feeding bacteria for cellular slime molds would increase their application potential in drug discovery.

Paraburkholderia symbionts isolated from Dictyostelium discoideum induce bacterial carriage in other Dictyostelium species

The social amoeba Dictyostelium discoideum engages in a complex relationship with bacterial endosymbionts in the genus Paraburkholderia, which can benefit their host by imbuing it with the ability to carry prey bacteria throughout its life cycle. The relationship between D. discoideum and Paraburkholderia has been shown to take place across many strains and a large geographical area, but little is known about Paraburkholderia's potential interaction with other dictyostelid species. We explore the ability of three Paraburkholderia species to stably infect and induce bacterial carriage in other dictyostelid hosts. We found that all three Paraburkholderia species successfully infected and induced carriage in seven species of Dictyostelium hosts. While the overall behaviour was qualitatively similar to that previously observed in infections of D. discoideum, differences in the outcomes of different host/symbiont combinations suggest a degree of specialization between partners. Paraburkholderia was unable to maintain a stable association with the more distantly related host Polysphondylium violaceum. Our results suggest that the mechanisms and evolutionary history of Paraburkholderia's symbiotic relationships may be general within Dictyostelium hosts, but not so general that it can associate with hosts of other genera. Our work further develops an emerging model system for the study of symbiosis in microbes.

Single-cell phenotypic plasticity modulates social behavior in Dictyostelium discoideum

In Dictyostelium chimeras, "cheaters" are strains that positively bias their contribution to the pool of spores, i.e., the reproductive cells resulting from development. On evolutionary time scales, the selective advantage; thus, gained by cheaters is predicted to undermine collective functions whenever social behaviors are genetically determined. Genotypes; however, are not the sole determinant of spore bias, but the relative role of genetic and plastic differences in evolutionary success is unclear. Here, we study chimeras composed of cells harvested in different phases of population growth. We show that such heterogeneity induces frequency-dependent, plastic variation in spore bias. In genetic chimeras, the magnitude of such variation is not negligible and can even reverse the classification of a strain's social behavior. Our results suggest that differential cell mechanical properties can underpin, through biases emerging during aggregation, a "lottery" in strains' reproductive success that may counter the evolution of cheating.

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.

Spindly is a nucleocytosolic O-fucosyltransferase in Dictyostelium and related proteins are widespread in protists and bacteria

O-GlcNAcylation is a prominent modification of nuclear and cytoplasmic proteins in animals and plants and is mediated by a single O-GlcNAc transferase (OGT). Spindly (Spy), a paralog of OGT first discovered in higher plants, has an ortholog in the apicomplexan parasite Toxoplasma gondii, and both enzymes are now recognized as O-fucosyltransferases (OFTs). Here we investigate the evolution of spy-like genes and experimentally confirm OFT activity in the social amoeba Dictyostelium-a protist that is more related to fungi and metazoa. Immunofluorescence probing with the fucose-specific Aleuria aurantia lectin (AAL) and biochemical cell fractionation combined with western blotting suggested the occurrence of nucleocytoplasmic fucosylation. The absence of reactivity in mutants deleted in spy or gmd (unable to synthesize GDP-Fuc) suggested monofucosylation mediated by Spy. Genetic ablation of the modE locus, previously predicted to encode a GDP-fucose transporter, confirmed its necessity for fucosylation in the secretory pathway but not for the nucleocytoplasmic proteins. Affinity capture of these proteins combined with mass spectrometry confirmed monofucosylation of Ser and Thr residues of several known nucleocytoplasmic proteins. As in Toxoplasma, the Spy OFT was required for optimal proliferation of Dictyostelium under laboratory conditions. These findings support a new phylogenetic analysis of OGT and OFT evolution that indicates their occurrence in the last eukaryotic common ancestor but mostly complementary presence in its eukaryotic descendants with the notable exception that both occur in red algae and plants. Their generally exclusive expression, high degree of conservation, and shared monoglycosylation targets suggest overlapping roles in physiological regulation.

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.

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.

Role of an FNIP Repeat Domain-Containing Protein Encoded by Megavirus Baoshan during Viral Infection

FNIP repeat domain-containing protein (FNIP protein) is a little-studied atypical leucine-rich repeat domain-containing protein found in social amoebae and mimiviruses. Here, a recently reported mimivirus of lineage C, Megavirus baoshan, was analyzed for FNIP protein genes. A total of 82 FNIP protein genes were identified, each containing up to 26 copies of the FNIP repeat, and mostly having an F-box domain at the N terminus. Both nucleotide and amino acid sequences of FNIP repeat were highly conserved. Most of the FNIP protein genes clustered together tandemly in groups of two to 14 genes. Nearly all FNIP protein genes shared similar expression patterns and were expressed 4 to 9 h postinfection. A typical viral FNIP protein, Mb0983, was selected for functional analysis. Protein interactome analysis identified two small GTPases, Rap1B and Rab7A, that interacted with Mb0983 in cytoplasm. The overexpression of Mb0983 in Acanthamoeba castellanii accelerated the degradation of Rap1B and Rab7A during viral infection. Mb0983 also interacted with host SKP1 and cullin-1, which were conserved components of the SKP1-cullin-1-F-box protein (SCF)-type ubiquitin E3 ligase complex. Deletion of the F-box domain of Mb0983 not only abolished its interaction with SKP1 and cullin-1 but also returned the speed of Rap1B and Rab7A degradation to normal in infected A. castellanii. These results suggested that Mb0983 is a part of the SCF-type ubiquitin E3 ligase complex and plays a role in the degradation of Rap1B and Rab7A. They also implied that other viral F-box-containing FNIP proteins might have similar effects on various host proteins. IMPORTANCE Megavirus baoshan encodes 82 FNIP proteins, more than any other reported mimiviruses. Their genetic and transcriptional features suggest that they are important for virus infection and adaption. Since most mimiviral FNIP proteins have the F-box domain, they were predicted to be involved in protein ubiquitylation. FNIP protein Mb0983 interacted with host SKP1 and cullin-1 through the F-box domain, supporting the idea that it is a part of the SCF-type ubiquitin E3 ligase complex. The substrates of Mb0983 for degradation were identified as the host small GTPases Rap1B and Rab7A. Combining the facts of the presence of a large number of FNIP genes in megavirus genomes, the extremely high expression level of the viral ubiquitin gene, and the reported observation that 35% of megavirus-infected amoeba cells died without productive infection, it is likely that megavirus actively explores the host ubiquitin-proteasome pathway in infection and that viral FNIP proteins play roles in the process.

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.

Lineage-Specific Genes and Family Expansions in Dictyostelid Genomes Display Expression Bias and Evolutionary Diversification during Development

Gene duplications generate new genes that can contribute to expression changes and the evolution of new functions. Genomes often consist of gene families that undergo expansions, some of which occur in specific lineages that reflect recent adaptive diversification. In this study, lineage-specific genes and gene family expansions were studied across five dictyostelid species to determine when and how they are expressed during multicellular development. Lineage-specific genes were found to be enriched among genes with biased expression (predominant expression in one developmental stage) in each species and at most developmental time points, suggesting independent functional innovations of new genes throughout the phylogeny. Biased duplicate genes had greater expression divergence than their orthologs and paralogs, consistent with subfunctionalization or neofunctionalization. Lineage-specific expansions in particular had biased genes with both molecular signals of positive selection and high expression, suggesting adaptive genetic and transcriptional diversification following duplication. Our results present insights into the potential contributions of lineage-specific genes and families in generating species-specific phenotypes during multicellular development in dictyostelids.

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.

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.

Glycomics, Glycoproteomics, and Glycogenomics: An Inter-Taxa Evolutionary Perspective

Glycosylation is a highly diverse set of co- and posttranslational modifications of proteins. For mammalian glycoproteins, glycosylation is often site-, tissue-, and species-specific and diversified by microheterogeneity. Multitudinous biochemical, cellular, physiological, and organismic effects of their glycans have been revealed, either intrinsic to the carrier proteins or mediated by endogenous reader proteins with carbohydrate recognition domains. Furthermore, glycans frequently form the first line of access by or defense from foreign invaders, and new roles for nucleocytoplasmic glycosylation are blossoming. We now know enough to conclude that the same general principles apply in invertebrate animals and unicellular eukaryotes-different branches of which spawned the plants or fungi and animals. The two major driving forces for exploring the glycomes of invertebrates and protists are (i) to understand the biochemical basis of glycan-driven biology in these organisms, especially of pathogens, and (ii) to uncover the evolutionary relationships between glycans, their biosynthetic enzyme genes, and biological functions for new glycobiological insights. With an emphasis on emerging areas of protist glycobiology, here we offer an overview of glycan diversity and evolution, to promote future access to this treasure trove of glycobiological processes.

Cytokinins in Dictyostelia - A Unique Model for Studying the Functions of Signaling Agents From Species to Kingdoms

Cytokinins (CKs) are a diverse group of evolutionarily significant growth-regulating molecules. While the CK biosynthesis and signal transduction pathways are the most well-understood in plant systems, these molecules have been identified in all kingdoms of life. This review follows the recent discovery of an expanded CK profile in the social amoeba, Dictyostelium discoideum. A comprehensive review on the present knowledge of CK biosynthesis, signal transduction, and CK-small molecule interactions within members of Dictyostelia will be summarized. In doing so, the utility of social amoebae will be highlighted as a model system for studying the evolution of these hormone-like signaling agents, which will set the stage for future research in this area.

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.

Description of Virulent Factors and Horizontal Gene Transfers of Keratitis-Associated Amoeba Acanthamoeba Triangularis by Genome Analysis

Acanthamoeba triangularis strain SH 621 is a free-living amoeba belonging to Acanthamoeba ribo-genotype T4. This ubiquitous protist is among the free-living amoebas responsible for Acanthamoeba keratitis, a severe infection of human cornea. Genome sequencing and genomic comparison were carried out to explore the biological functions and to better understand the virulence mechanism related to the pathogenicity of Acanthamoeba keratitis. The genome assembly harbored a length of 66.43 Mb encompassing 13,849 scaffolds. The analysis of predicted proteins reported the presence of 37,062 ORFs. A complete annotation revealed 33,168 and 16,605 genes that matched with NCBI non-redundant protein sequence (nr) and Cluster of Orthologous Group of proteins (COG) databases, respectively. The Kyoto Encyclopedia of Genes and Genomes Pathway (KEGG) annotation reported a great number of genes related to carbohydrate, amino acid and lipid metabolic pathways. The pangenome performed with 8 available amoeba genomes belonging to genus Acanthamoeba revealed a core genome containing 843 clusters of orthologous genes with a ratio core genome/pangenome of less than 0.02. We detected 48 genes related to virulent factors of Acanthamoeba keratitis. Best hit analyses in nr database identified 99 homologous genes shared with amoeba-resisting microorganisms. This study allows the deciphering the genome of a free-living amoeba with medical interest and provides genomic data to better understand virulence-related Acanthamoeba keratitis.

Investigation of potential pathogenicity of Willaertia magna by investigating the transfer of bacteria pathogenicity genes into its genome

Willaertia magna c2c maky is a thermophilic amoeba closely related to the genus Naegleria. This free-living amoeba has the ability to eliminate Legionella pneumophila, which is an amoeba-resisting bacterium living in an aquatic environment. To prevent the proliferation of L. pneumophila in cooling towers, the use of W. magna as natural biocide has been proposed. To provide a better understanding of the W. magna genome, whole-genome sequencing was performed through the study of virulence factors and lateral gene transfers. This amoeba harbors a genome of 36.5 megabases with 18,519 predicted genes. BLASTp analyses reported protein homology between 136 W. magna sequences and amoeba-resistant microorganisms. Horizontal gene transfers were observed based on the basis of the phylogenetic reconstruction hypothesis. We detected 15 homologs of N. fowleri genes related to virulence, although these latter were also found in the genome of N. gruberi, which is a non-pathogenic amoeba. Furthermore, the cytotoxicity test performed on human cells supports the hypothesis that the strain c2c maky is a non-pathogenic amoeba. This work explores the genomic repertory for the first draft genome of genus Willaertia and provides genomic data for further comparative studies on virulence of related pathogenic amoeba, N. fowleri.

Conditional expression explains molecular evolution of social genes in a microbe

Conflict is thought to play a critical role in the evolution of social interactions by promoting diversity or driving accelerated evolution. However, despite our sophisticated understanding of how conflict shapes social traits, we have limited knowledge of how it impacts molecular evolution across the underlying social genes. Here we address this problem by analyzing the genome-wide impact of social interactions using genome sequences from 67 Dictyostelium discoideum strains. We find that social genes tend to exhibit enhanced polymorphism and accelerated evolution. However, these patterns are not consistent with conflict driven processes, but instead reflect relaxed purifying selection. This pattern is most likely explained by the conditional nature of social interactions, whereby selection on genes expressed only in social interactions is diluted by generations of inactivity. This dilution of selection by inactivity enhances the role of drift, leading to increased polymorphism and accelerated evolution, which we call the Red King process.

Nucleocytoplasmic O-glycosylation in protists

O-Glycosylation is an increasingly recognized modification of intracellular proteins in all kingdoms of life, and its occurrence in protists has been investigated to understand its evolution and its roles in the virulence of unicellular pathogens. We focus here on two kinds of glycoregulation found in unicellular eukaryotes: one is a simple O-fucose modification of dozens if not hundreds of Ser/Thr-rich proteins, and the other a complex pentasaccharide devoted to a single protein associated with oxygen sensing and the assembly of polyubiquitin chains. These modifications are not required for life but contingently modulate biological processes in the social amoeba Dictyostelium and the human pathogen Toxoplasma gondii, and likely occur in diverse unicellular protists. O-Glycosylation that is co-localized in the cytoplasm allows for glycoregulation over the entire life of the protein, contrary to the secretory pathway where glycosylation usually occurs before its delivery to its site of function. Here, we interpret cellular roles of nucleocytoplasmic glycans in terms of current evidence for their effects on the conformation and dynamics of protist proteins, to serve as a guide for future studies to examine their broader significance.

The origins and evolution of macropinocytosis

In macropinocytosis, cells take up micrometre-sized droplets of medium into internal vesicles. These vesicles are acidified and fused to lysosomes, their contents digested and useful compounds extracted. Indigestible contents can be exocytosed. Macropinocytosis has been known for approaching 100 years and is described in both metazoa and amoebae, but not in plants or fungi. Its evolutionary origin goes back to at least the common ancestor of the amoebozoa and opisthokonts, with apparent secondary loss from fungi. The primary function of macropinocytosis in amoebae and some cancer cells is feeding, but the conserved processing pathway for macropinosomes, which involves shrinkage and the retrieval of membrane to the cell surface, has been adapted in immune cells for antigen presentation. Macropinocytic cups are large actin-driven processes, closely related to phagocytic cups and pseudopods and appear to be organized around a conserved signalling patch of PIP3, active Ras and active Rac that directs actin polymerization to its periphery. Patches can form spontaneously and must be sustained by excitable kinetics with strong cooperation from the actin cytoskeleton. Growth-factor signalling shares core components with macropinocytosis, based around phosphatidylinositol 3-kinase (PI3-kinase), and we suggest that it evolved to take control of ancient feeding structures through a coupled growth factor receptor. This article is part of the Theo Murphy meeting issue 'Macropinocytosis'.

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.

Cellular allorecognition and its roles in Dictyostelium development and social evolution

The social amoeba Dictyostelium discoideum is a tractable model organism to study cellular allorecognition, which is the ability of a cell to distinguish itself and its genetically similar relatives from more distantly related organisms. Cellular allorecognition is ubiquitous across the tree of life and affects many biological processes. Depending on the biological context, these versatile systems operate both within and between individual organisms, and both promote and constrain functional heterogeneity. Some of the most notable allorecognition systems mediate neural self-avoidance in flies and adaptive immunity in vertebrates. D. discoideum's allorecognition system shares several structures and functions with other allorecognition systems. Structurally, its key regulators reside at a single genomic locus that encodes two highly polymorphic proteins, a transmembrane ligand called TgrC1 and its receptor TgrB1. These proteins exhibit isoform-specific, heterophilic binding across cells. Functionally, this interaction determines the extent to which co-developing D. discoideum strains co-aggregate or segregate during the aggregation phase of multicellular development. The allorecognition system thus affects both development and social evolution, as available evidence suggests that the threat of developmental cheating represents a primary selective force acting on it. Other significant characteristics that may inform the study of allorecognition in general include that D. discoideum's allorecognition system is a continuous and inclusive trait, it is pleiotropic, and it is temporally regulated.

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.

Dictyostelium: An Important Source of Structural and Functional Diversity in Drug Discovery

The cellular slime mold Dictyostelium discoideum is an excellent model organism for the study of cell and developmental biology because of its simple life cycle and ease of use. Recent findings suggest that Dictyostelium and possibly other genera of cellular slime molds, are potential sources of novel lead compounds for pharmacological and medical research. In this review, we present supporting evidence that cellular slime molds are an untapped source of lead compounds by examining the discovery and functions of polyketide differentiation-inducing factor-1, a compound that was originally isolated as an inducer of stalk-cell differentiation in D. discoideum and, together with its derivatives, is now a promising lead compound for drug discovery in several areas. We also review other novel compounds, including secondary metabolites, that have been isolated from cellular slime molds.

Species recognition in social amoebae

Aggregative multicellularity requires the ability of cells to recognise conspecifics. Social amoebae are among the best studied of such organisms, but the mechanism and evolutionary background of species recognition remained to be investigated. Here we show that heterologous expression of a single Dictyostelium purpureum gene is sufficient for D. discoideum cells to efficiently make chimaeric fruiting bodies with D. purpureum cells. This gene forms a bidirectional pair with another gene on the D. purpureum genome, and they are both highly polymorphic among independent wild isolates of the same mating group that do not form chimaeric fruiting bodies with each other. These paired genes are both structurally similar to D. discoideum tgrB1/C1 pair, which is responsible for clonal discrimination within that species, suggesting that these tgr genes constitute the species recognition system that has attained a level of precision capable of discriminating between clones within a species. Analysis of the available genome sequences of social amoebae revealed that such gene pairs exist only within the clade composed of species that produce precursors of sterile stalk cells (prestalk cells), suggesting concurrent evolution of a precise allorecognition system and a new 'worker' cell-type dedicated to transporting and supporting the reproductive cells.

CRISPR/Cas9 and glycomics tools for Toxoplasma glycobiology

Infection with the protozoan parasite Toxoplasma gondii is a major health risk owing to birth defects, its chronic nature, ability to reactivate to cause blindness and encephalitis, and high prevalence in human populations. Unlike most eukaryotes, Toxoplasma propagates in intracellular parasitophorous vacuoles, but like nearly all other eukaryotes, Toxoplasma glycosylates many cellular proteins and lipids and assembles polysaccharides. Toxoplasma glycans resemble those of other eukaryotes, but species-specific variations have prohibited deeper investigations into their roles in parasite biology and virulence. The Toxoplasma genome encodes a suite of likely glycogenes expected to assemble N-glycans, O-glycans, a C-glycan, GPI-anchors, and polysaccharides, along with their precursors and membrane transporters. To investigate the roles of specific glycans in Toxoplasma, here we coupled genetic and glycomics approaches to map the connections between 67 glycogenes, their enzyme products, the glycans to which they contribute, and cellular functions. We applied a double-CRISPR/Cas9 strategy, in which two guide RNAs promote replacement of a candidate gene with a resistance gene; adapted MS-based glycomics workflows to test for effects on glycan formation; and infected fibroblast monolayers to assess cellular effects. By editing 17 glycogenes, we discovered novel Glc_0-2-Man₆-GlcNAc₂-type N-glycans, a novel HexNAc-GalNAc-mucin-type O-glycan, and Tn-antigen; identified the glycosyltransferases for assembling novel nuclear O-Fuc-type and cell surface Glc-Fuc-type O-glycans; and showed that they are important for in vitro growth. The guide sequences, editing constructs, and mutant strains are freely available to researchers to investigate the roles of glycans in their favorite biological processes.

Chemical compounds from Dictyostelium discoideum repel a plant-parasitic nematode and can protect roots

Slime mold species in the genus Dictyostelium are considered to have a close relationship with non-parasitic nematodes; they are sympatric in soils and can exhibit interspecific competition for food. We investigated whether this relationship extends to a plant-parasitic nematode that is active in the rhizosphere and has broad host specificity, damaging crops worldwide. Using a novel assay to examine the interaction between the cellular slime mold, Dictyostelium discoideum, and the plant-parasitic nematodes, Meloidogyne spp., we found that cellular slime molds can repel plant parasitic nematodes. Specifically, the repulsion activity was in response to chemical compounds released by cellular slime mold fruiting bodies. Under laboratory conditions, these soluble chemical extracts from fruiting bodies of D. discoideum showed repulsion activity strong enough to protect plant roots. The fruiting body cell extracts repelled but were not toxic to the plant-parasitic nematodes.

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.

Spo11-Independent Meiosis in Social Amoebae

Sex in social amoebae (or dictyostelids) has a number of striking features. Dictyostelid zygotes do not proliferate but grow to a large size by feeding on other cells of the same species, each zygote ultimately forming a walled structure called a macrocyst. The diploid macrocyst nucleus undergoes meiosis, after which a single meiotic product survives to restart haploid vegetative growth. Meiotic recombination is generally initiated by the Spo11 enzyme, which introduces DNA double-strand breaks. Uniquely, as far as is known among sexual eukaryotes, dictyostelids lack a SPO11 gene. Despite this, recombination occurs at high frequencies during meiosis in dictyostelids, through unknown mechanisms. The molecular processes underlying these events, and the evolutionary drivers that brought them into being, may shed light on the genetic conflicts that occur within and between genomes, and how they can be resolved.

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.

Eat Prey, Live: Dictyostelium discoideum As a Model for Cell-Autonomous Defenses

The soil-dwelling social amoeba Dictyostelium discoideum feeds on bacteria. Each meal is a potential infection because some bacteria have evolved mechanisms to resist predation. To survive such a hostile environment, D. discoideum has in turn evolved efficient antimicrobial responses that are intertwined with phagocytosis and autophagy, its nutrient acquisition pathways. The core machinery and antimicrobial functions of these pathways are conserved in the mononuclear phagocytes of mammals, which mediate the initial, innate-immune response to infection. In this review, we discuss the advantages and relevance of D. discoideum as a model phagocyte to study cell-autonomous defenses. We cover the antimicrobial functions of phagocytosis and autophagy and describe the processes that create a microbicidal phagosome: acidification and delivery of lytic enzymes, generation of reactive oxygen species, and the regulation of Zn²⁺, Cu²⁺, and Fe²⁺ availability. High concentrations of metals poison microbes while metal sequestration inhibits their metabolic activity. We also describe microbial interference with these defenses and highlight observations made first in D. discoideum. Finally, we discuss galectins, TNF receptor-associated factors, tripartite motif-containing proteins, and signal transducers and activators of transcription, microbial restriction factors initially characterized in mammalian phagocytes that have either homologs or functional analogs in D. discoideum.

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.

Natural Products from Social Amoebae

Natural products are invaluable sources of structural diversity and complexity ideally suited for the development of therapeutic agents. The search for novel bioactive molecules has prompted scientists to explore various ecological niches. Microorganisms have been shown to constitute such an important source. Despite their biosynthetic potential, social amoebae, that is, microorganisms with both a uni- and multicellular lifestyle, are underexplored regarding their secreted secondary metabolome. In this review, we present the structural diversity of amoebal natural products and discuss their biological functions as well as their total syntheses.

Isolation, Synthesis, and Biological Activity of Chlorinated Alkylresorcinols from Dictyostelium Cellular Slime Molds

Eight chlorinated alkylresorcinols, monochasiol A-H (1-8), were isolated from the fruiting bodies of Dictyostelium monochasioides. Compounds 1-8 were synthesized to confirm their structures and to obtain sufficient material for performing biological tests. Monochasiol A (1) selectively inhibited the concanavalin A-induced interleukin-2 production in Jurkat cells, a human T lymphocyte cell line. Monochasiols were biogenetically synthesized by the combination of biosynthetic enzymes relating to the principal polyketides, MPBD and DIF-1, produced by Dictyostelium discoideum.

The Origin of Animal Multicellularity and Cell Differentiation

Over 600 million years ago, animals evolved from a unicellular or colonial organism whose cell(s) captured bacteria with a collar complex, a flagellum surrounded by a microvillar collar. Using principles from evolutionary cell biology, we reason that the transition to multicellularity required modification of pre-existing mechanisms for extracellular matrix synthesis and cytokinesis. We discuss two hypotheses for the origin of animal cell types: division of labor from ancient plurifunctional cells and conversion of temporally alternating phenotypes into spatially juxtaposed cell types. Mechanistic studies in diverse animals and their relatives promise to deepen our understanding of animal origins and cell biology.

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.

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.

A core phylogeny of Dictyostelia inferred from genomes representative of the eight major and minor taxonomic divisions of the group

BACKGROUND

Dictyostelia are a well-studied group of organisms with colonial multicellularity, which are members of the mostly unicellular Amoebozoa. A phylogeny based on SSU rDNA data subdivided all Dictyostelia into four major groups, but left the position of the root and of six group-intermediate taxa unresolved. Recent phylogenies inferred from 30 or 213 proteins from sequenced genomes, positioned the root between two branches, each containing two major groups, but lacked data to position the group-intermediate taxa. Since the positions of these early diverging taxa are crucial for understanding the evolution of phenotypic complexity in Dictyostelia, we sequenced six representative genomes of early diverging taxa.

RESULTS

We retrieved orthologs of 47 housekeeping proteins with an average size of 890 amino acids from six newly sequenced and eight published genomes of Dictyostelia and unicellular Amoebozoa and inferred phylogenies from single and concatenated protein sequence alignments. Concatenated alignments of all 47 proteins, and four out of five subsets of nine concatenated proteins all produced the same consensus phylogeny with 100% statistical support. Trees inferred from just two out of the 47 proteins, individually reproduced the consensus phylogeny, highlighting that single gene phylogenies will rarely reflect correct species relationships. However, sets of two or three concatenated proteins again reproduced the consensus phylogeny, indicating that a small selection of genes suffices for low cost classification of as yet unincorporated or newly discovered dictyostelid and amoebozoan taxa by gene amplification.

CONCLUSIONS

The multi-locus consensus phylogeny shows that groups 1 and 2 are sister clades in branch I, with the group-intermediate taxon D. polycarpum positioned as outgroup to group 2. Branch II consists of groups 3 and 4, with the group-intermediate taxon Polysphondylium violaceum positioned as sister to group 4, and the group-intermediate taxon Dictyostelium polycephalum branching at the base of that whole clade. Given the data, the approximately unbiased test rejects all alternative topologies favoured by SSU rDNA and individual proteins with high statistical support. The test also rejects monophyletic origins for the genera Acytostelium, Polysphondylium and Dictyostelium. The current position of Acytostelium ellipticum in the consensus phylogeny indicates that somatic cells were lost twice in Dictyostelia.

A set of genes conserved in sequence and expression traces back the establishment of multicellularity in social amoebae

Background

The developmental cycle of Dictyostelid amoebae represents an early form of multicellularity with cell type differentiation. Mutant studies in the model Dictyostelium discoideum revealed that its developmental program integrates the actions of genes involved in signal transduction, adhesion, motility, autophagy and cell wall and matrix biosynthesis. However, due to functional redundancy and fail safe options not required in the laboratory, this single organism approach cannot capture all essential genes.

To understand how multicellular organisms evolved, it is essential to recognize both the conserved core features of their developmental programs and the gene modifications that instigated phenotypic innovation. For complex organisms, such as animals, this is not within easy reach, but it is feasible for less complex forms, such as the Dictyostelid social amoebas.

Results

We compared global profiles of gene expression during the development of four social amoebae species that represent 600 mya of Dictyostelia evolution, and identified orthologous conserved genes with similar developmental up-regulation of expression using three different methods. For validation, we disrupted five genes of this core set and examined the phenotypic consequences.

Conclusion

At least 71 of the developmentally regulated genes that were identified with all methods were likely to be already present in the last ancestor of all Dictyostelia. The lack of phenotypic changes in null mutants indicates that even highly conserved genes either participate in functionally redundant pathways or are necessary for developmental progression under adverse, non-standard laboratory conditions. Both mechanisms provide robustness to the developmental program, but impose a limit on the information that can be obtained from deleting single genes.

Electronic supplementary material

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

Evolutionary analyses of mitochondrial carrier family of dictyostelids

The transportation of solutes across the inner membrane of the mitochondria is catalyzed by a nuclear-coded family of transport proteins called mitochondrial carriers (MCs). Sequences from dictyostelid genome projects have facilitated analysis of the evolution of the dictyostelid mitochondrial carrier family (MCF). The average evolutionary distances between various regions in the MCF shows that the transmembrane region (TR) and conical pit region (CPR) are the only two conserved structural regions. A phylogenetic tree built using the concatenated orthologous TR and CPR sequences of 7 MCs showed that dictyostelids are similar to metazoans in this way. A close evolutionary relationship was observed between dictyostelids and metazoans in 4 MCs known to be related to ADP/ATP transport (MAA). This was further evidenced by the fact that dictyostelids have undergone gene expansion similar to that of metazoans during the evolution of MAA. Sequence logo analysis of CPR in MAA showed that dictyostelids have motifs similar to those of Metazoa. Combined with the conserved substrate binding site of 7 MCs in eukaryotes, it is postulated that dictyostelids are closely related to Metazoa with respect to the evolution of MAA.

Acanthamoeba and Dictyostelium Use Different Foraging Strategies

Amoeba often use cell movement as a mechanism to find food, such as bacteria, in their environment. The chemotactic movement of the soil amoeba Dictyostelium to folate or other pterin compounds released by bacteria is a well-documented foraging mechanism. Acanthamoeba can also feed on bacteria but relatively little is known about the mechanism(s) by which this amoeba locates bacteria. Acanthamoeba movement in the presence of folate or bacteria was analyzed in above agar assays and compared to that observed for Dictyostelium. The overall mobility of Acanthamoeba was robust like that of Dictyostelium but Acanthamoeba did not display a chemotactic response to folate. In the presence of bacteria, Acanthamoeba only showed a marginal bias in directed movement whereas Dictyostelium displayed a strong chemotactic response. A comparison of genomes revealed that Acanthamoeba and Dictyostelium share some similarities in G protein signaling components but that specific G proteins used in Dictyostelium chemotactic responses were not present in current Acanthamoeba genome sequence data. The results of this study suggest that Acanthamoeba does not use chemotaxis as the primary mechanism to find bacterial food sources and that the chemotactic responses of Dictyostelium to bacteria may have co-evolved with chemotactic responses that facilitate multicellular development.

Convergent evolution of tRNA gene targeting preferences in compact genomes

Background

In gene-dense genomes, mobile elements are confronted with highly selective pressure to amplify without causing excessive damage to the host. The targeting of tRNA genes as potentially safe integration sites has been developed by retrotransposons in various organisms such as the social amoeba Dictyostelium discoideum and the yeast Saccharomyces cerevisiae. In D. discoideum, tRNA gene-targeting retrotransposons have expanded to approximately 3 % of the genome. Recently obtained genome sequences of species representing the evolutionary history of social amoebae enabled us to determine whether the targeting of tRNA genes is a generally successful strategy for mobile elements to colonize compact genomes.

Results

During the evolution of dictyostelids, different retrotransposon types independently developed the targeting of tRNA genes at least six times. DGLT-A elements are long terminal repeat (LTR) retrotransposons that display integration preferences ~15 bp upstream of tRNA gene-coding regions reminiscent of the yeast Ty3 element. Skipper elements are chromoviruses that have developed two subgroups: one has canonical chromo domains that may favor integration in centromeric regions, whereas the other has diverged chromo domains and is found ~100 bp downstream of tRNA genes. The integration of D. discoideum non-LTR retrotransposons ~50 bp upstream (TRE5 elements) and ~100 bp downstream (TRE3 elements) of tRNA genes, respectively, likely emerged at the root of dictyostelid evolution. We identified two novel non-LTR retrotransposons unrelated to TREs: one with a TRE5-like integration behavior and the other with preference ~4 bp upstream of tRNA genes.

Conclusions

Dictyostelid retrotransposons demonstrate convergent evolution of tRNA gene targeting as a probable means to colonize the compact genomes of their hosts without being excessively mutagenic. However, high copy numbers of tRNA gene-associated retrotransposons, such as those observed in D. discoideum, are an exception, suggesting that the targeting of tRNA genes does not necessarily favor the amplification of position-specific integrating elements to high copy numbers under the repressive conditions that prevail in most host cells.

Electronic supplementary material

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

The multicellularity genes of dictyostelid social amoebas

The evolution of multicellularity enabled specialization of cells, but required novel signalling mechanisms for regulating cell differentiation. Early multicellular organisms are mostly extinct and the origins of these mechanisms are unknown. Here using comparative genome and transcriptome analysis across eight uni- and multicellular amoebozoan genomes, we find that 80% of proteins essential for the development of multicellular Dictyostelia are already present in their unicellular relatives. This set is enriched in cytosolic and nuclear proteins, and protein kinases. The remaining 20%, unique to Dictyostelia, mostly consists of extracellularly exposed and secreted proteins, with roles in sensing and recognition, while several genes for synthesis of signals that induce cell-type specialization were acquired by lateral gene transfer. Across Dictyostelia, changes in gene expression correspond more strongly with phenotypic innovation than changes in protein functional domains. We conclude that the transition to multicellularity required novel signals and sensors rather than novel signal processing mechanisms.

Unicellular social amoebae aggregate to form a multicellular life stage, making them a model system for the evolution of multicellularity. Here, Glöckner et al. use a comparative genomic and transcriptomic approach to determine the origin of the genes essential for multicellularity in the social amoebae.