Comparative genome and transcriptome analyses of the social amoeba Acytostelium subglobosum that accomplishes multicellular development without germ-soma differentiation

MarFERReT, an open-source, version-controlled reference library of marine microbial eukaryote functional genes

Metatranscriptomics generates large volumes of sequence data about transcribed genes in natural environments. Taxonomic annotation of these datasets depends on availability of curated reference sequences. For marine microbial eukaryotes, current reference libraries are limited by gaps in sequenced organism diversity and barriers to updating libraries with new sequence data, resulting in taxonomic annotation of about half of eukaryotic environmental transcripts. Here, we introduce Marine Functional EukaRyotic Reference Taxa (MarFERReT), a marine microbial eukaryotic sequence library designed for use with taxonomic annotation of eukaryotic metatranscriptomes. We gathered 902 publicly accessible marine eukaryote genomes and transcriptomes and assessed their sequence quality and cross-contamination issues, selecting 800 validated entries for inclusion in MarFERReT. Version 1.1 of MarFERReT contains reference sequences from 800 marine eukaryotic genomes and transcriptomes, covering 453 species- and strain-level taxa, totaling nearly 28 million protein sequences with associated NCBI and PR2 Taxonomy identifiers and Pfam functional annotations. The MarFERReT project repository hosts containerized build scripts, documentation on installation and use case examples, and information on new versions of MarFERReT.

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.

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.

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.

Spliceosomal Prp8 intein at the crossroads of protein and RNA splicing

The spliceosome is a large ribonucleoprotein complex that removes introns from pre-mRNAs. At its functional core lies the essential pre-mRNA processing factor 8 (Prp8) protein. Across diverse eukaryotes, this protein cofactor of RNA catalysis harbors a self-splicing element called an intein. Inteins in Prp8 are extremely pervasive and are found at 7 different sites in various species. Here, we focus on the Prp8 intein from Cryptococcus neoformans (Cne), a human fungal pathogen. We solved the crystal structure of this intein, revealing structural homology among protein splicing sequences in eukaryotes, including the Hedgehog C terminus. Working with the Cne Prp8 intein in a reporter assay, we find that the biologically relevant divalent metals copper and zinc inhibit intein splicing, albeit by 2 different mechanisms. Copper likely stimulates reversible modifications on a catalytically important cysteine, whereas zinc binds at the terminal asparagine and the same critical cysteine. Importantly, we also show that copper treatment inhibits Prp8 protein splicing in Cne. Lastly, an intein-containing Prp8 precursor model is presented, suggesting that metal-induced protein splicing inhibition would disturb function of both Prp8 and the spliceosome. These results indicate that Prp8 protein splicing can be modulated, with potential functional implications for the spliceosome.

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.

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 genome Tree of Life for the Fungi kingdom

Fungi belong to one of the largest and most diverse kingdoms of living organisms. The evolutionary kinship within a fungal population has so far been inferred mostly from the gene-information-based trees ("gene trees"), constructed commonly based on the degree of differences of proteins or DNA sequences of a small number of highly conserved genes common among the population by a multiple sequence alignment (MSA) method. Since each gene evolves under different evolutionary pressure and time scale, it has been known that one gene tree for a population may differ from other gene trees for the same population depending on the subjective selection of the genes. Within the last decade, a large number of whole-genome sequences of fungi have become publicly available, which represent, at present, the most fundamental and complete information about each fungal organism. This presents an opportunity to infer kinship among fungi using a whole-genome information-based tree ("genome tree"). The method we used allows comparison of whole-genome information without MSA, and is a variation of a computational algorithm developed to find semantic similarities or plagiarism in two books, where we represent whole-genomic information of an organism as a book of words without spaces. The genome tree reveals several significant and notable differences from the gene trees, and these differences invoke new discussions about alternative narratives for the evolution of some of the currently accepted fungal groups.

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.

The Evolution of Aggregative Multicellularity and Cell-Cell Communication in the Dictyostelia

Aggregative multicellularity, resulting in formation of a spore-bearing fruiting body, evolved at least six times independently amongst both eukaryotes and prokaryotes. Amongst eukaryotes, this form of multicellularity is mainly studied in the social amoeba Dictyostelium discoideum. In this review, we summarise trends in the evolution of cell-type specialisation and behavioural complexity in the four major groups of Dictyostelia. We describe the cell-cell communication systems that control the developmental programme of D. discoideum, highlighting the central role of cAMP in the regulation of cell movement and cell differentiation. Comparative genomic studies showed that the proteins involved in cAMP signalling are deeply conserved across Dictyostelia and their unicellular amoebozoan ancestors. Comparative functional analysis revealed that cAMP signalling in D. discoideum originated from a second messenger role in amoebozoan encystation. We highlight some molecular changes in cAMP signalling genes that were responsible for the novel roles of cAMP in multicellular development.

Root of Dictyostelia based on 213 universal proteins

Dictyostelia are common soil microbes that can aggregate when starved to form multicellular fruiting bodies, a characteristic that has also led to their long history of study and widespread use as model systems. Ribosomal RNA phylogeny of Dictyostelia identified four major divisions (Groups 1-4), none of which correspond to traditional genera. Group 1 was also tentatively identified as sister lineage to the other three Groups, although not consistently or with strong support. We tested the dictyostelid root using universal protein-coding genes identified by exhaustive comparison of six completely sequenced dictyostelid genomes, which include representatives of all four major molecular Groups. A set of 213 genes are low-copy number in all genomes, present in at least one amoebozoan outgroup taxon (Acanthamoeba castellanii or Physarum polycephalum), and phylogenetically congruent. Phylogenetic analysis of a concatenation of the deduced protein sequences produces a single topology dividing Dictyostelia into two major divisions: Groups 1+2 and Groups 3+4. All clades in the tree are fully supported by maximum likelihood and Bayesian inference, and all alternative roots are unambiguously rejected by the approximately unbiased (AU) test. The 1+2, 3+4 root is also fully supported even after deleting clusters with strong individual support for this root, or concatenating all clusters with low support for alternative roots. The 213 putatively ancestral amoebozoan proteins encode a wide variety of functions including 21 KOG categories out of a total of 25. These comprehensive analyses and consistent results indicate that it is time for full taxonomic revision of Dictyostelia, which will also enable more effective exploitation of its unique potential as an evolutionary model system.