The genome of the social amoeba Dictyostelium discoideum

Between a Pod and a Hard Test: The Deep Evolution of Amoebae

Amoebozoa is the eukaryotic supergroup sister to Obazoa, the lineage that contains the animals and Fungi, as well as their protistan relatives, and the breviate and apusomonad flagellates. Amoebozoa is extraordinarily diverse, encompassing important model organisms and significant pathogens. Although amoebozoans are integral to global nutrient cycles and present in nearly all environments, they remain vastly understudied. We present a robust phylogeny of Amoebozoa based on broad representative set of taxa in a phylogenomic framework (325 genes). By sampling 61 taxa using culture-based and single-cell transcriptomics, our analyses show two major clades of Amoebozoa, Discosea, and Tevosa. This phylogeny refutes previous studies in major respects. Our results support the hypothesis that the last common ancestor of Amoebozoa was sexual and flagellated, it also may have had the ability to disperse propagules from a sporocarp-type fruiting body. Overall, the main macroevolutionary patterns in Amoebozoa appear to result from the parallel losses of homologous characters of a multiphase life cycle that included flagella, sex, and sporocarps rather than independent acquisition of convergent features.

The La and related RNA-binding proteins (LARPs): structures, functions, and evolving perspectives

La was first identified as a polypeptide component of ribonucleic protein complexes targeted by antibodies in autoimmune patients and is now known to be a eukaryote cell-ubiquitous protein. Structure and function studies have shown that La binds to a common terminal motif, UUU-3'-OH, of nascent RNA polymerase III (RNAP III) transcripts and protects them from exonucleolytic decay. For precursor-tRNAs, the most diverse and abundant of these transcripts, La also functions as an RNA chaperone that helps to prevent their misfolding. Related to this, we review evidence that suggests that La and its link to RNAP III were significant in the great expansions of the tRNAomes that occurred in eukaryotes. Four families of La-related proteins (LARPs) emerged during eukaryotic evolution with specialized functions. We provide an overview of the high-resolution structural biology of La and LARPs. LARP7 family members most closely resemble La but function with a single RNAP III nuclear transcript, 7SK, or telomerase RNA. A cytoplasmic isoform of La protein as well as LARPs 6, 4, and 1 function in mRNA metabolism and translation in distinct but similar ways, sometimes with the poly(A)-binding protein, and in some cases by direct binding to poly(A)-RNA. New structures of LARP domains, some complexed with RNA, provide novel insights into the functional versatility of these proteins. We also consider LARPs in relation to ancestral La protein and potential retention of links to specific RNA-related pathways. One such link may be tRNA surveillance and codon usage by LARP-associated mRNAs. WIREs RNA 2017, 8:e1430. doi: 10.1002/wrna.1430 For further resources related to this article, please visit the WIREs website.

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.

G-Protein Dependent Signal Transduction and Ubiquitination in Dictyostelium

Signal transduction through G-protein-coupled receptors (GPCRs) is central for the regulation of virtually all cellular functions, and it has been widely implicated in human diseases. These receptors activate a common molecular switch that is represented by the heterotrimeric G-protein generating a number of second messengers (cAMP, cGMP, DAG, IP3, Ca²⁺ etc.), leading to a plethora of diverse cellular responses. Spatiotemporal regulation of signals generated by a given GPCR is crucial for proper signalling and is accomplished by a series of biochemical modifications. Over the past few years, it has become evident that many signalling proteins also undergo ubiquitination, a posttranslational modification that typically leads to protein degradation, but also mediates processes such as protein-protein interaction and protein subcellular localization. The social amoeba Dictyostelium discoideum has proven to be an excellent model to investigate signal transduction triggered by GPCR activation, as cAMP signalling via GPCR is a major regulator of chemotaxis, cell differentiation, and multicellular morphogenesis. Ubiquitin ligases have been recently involved in these processes. In the present review, we will summarize the most significant pathways activated upon GPCRs stimulation and discuss the role played by ubiquitination in Dictyostelium cells.

Optimization of trans-Splicing for Huntington's Disease RNA Therapy

Huntington's disease (HD) is a devastating neurodegenerative disorder caused by a polyglutamine (polyQ) expansion in exon 1 of the Huntingtin (HTT) gene. We have previously demonstrated that spliceosome-mediated trans-splicing is a viable molecular strategy to specifically reduce and repair mutant HTT (mtHTT). Here, the targeted tethering efficacy of the pre-mRNA trans-splicing modules (PTM) in HTT was optimized. Various PTMs that targeted the 3′ end of HTT intron 1 or the intron 1 branch point were shown trans-splice into an HTT mini-gene, as well as the endogenous HTT pre-mRNA. PTMs that specifically target the endogenous intron 1 branch point increased the trans-splicing efficacy from 1–5 to 10–15%. Furthermore, lentiviral expression of PTMs in a human HD patient iPSC-derived neural culture significantly reversed two previously established polyQ-length dependent phenotypes. These results suggest that pre-mRNA repair of mtHTT could hold therapeutic benefit and it demonstrates an alternative platform to correct the mRNA product produced by the mtHTT allele in the context of HD.

Retrotransposon Domestication and Control in Dictyostelium discoideum

Transposable elements, identified in all eukaryotes, are mobile genetic units that can change their genomic position. Transposons usually employ an excision and reintegration mechanism, by which they change position, but not copy number. In contrast, retrotransposons amplify via RNA intermediates, increasing their genomic copy number. Hence, they represent a particular threat to the structural and informational integrity of the invaded genome. The social amoeba Dictyostelium discoideum, model organism of the evolutionary Amoebozoa supergroup, features a haploid, gene-dense genome that offers limited space for damage-free transposition. Several of its contemporary retrotransposons display intrinsic integration preferences, for example by inserting next to transfer RNA genes or other retroelements. Likely, any retrotransposons that invaded the genome of the amoeba in a non-directed manner were lost during evolution, as this would result in decreased fitness of the organism. Thus, the positional preference of the Dictyostelium retroelements might represent a domestication of the selfish elements. Likewise, the reduced danger of such domesticated transposable elements led to their accumulation, and they represent about 10% of the current genome of D. discoideum. To prevent the uncontrolled spreading of retrotransposons, the amoeba employs control mechanisms including RNA interference and heterochromatization. Here, we review TRE5-A, DIRS-1 and Skipper-1, as representatives of the three retrotransposon classes in D. discoideum, which make up 5.7% of the Dictyostelium genome. We compile open questions with respect to their mobility and cellular regulation, and suggest strategies, how these questions might be addressed experimentally.

Lipopolysaccharide Detection across the Kingdoms of Life

Studies in recent years have uncovered a diverse set of eukaryotic receptors that recognize lipopolysaccharide (LPS), the major outer-membrane component of Gram-negative bacteria. Indeed, Toll-like receptors, G-protein-coupled receptors, integrins, receptor-like kinases, and caspases have emerged as important LPS-interacting proteins. In this review, the mammalian receptors that detect LPS are described. I highlight how no host protein is involved in all LPS responses, but a single lipid (phosphatidylinositol-4,5-bisphosphate) regulates many LPS responses, including endocytosis, phagocytosis, inflammation, and pyroptosis. I further describe LPS response systems that operate specifically in plants, and discuss potentially new LPS response systems that await discovery. This diversity of receptors for a single microbial product underscores the importance of host-microbe interactions in multiple kingdoms of life.

Fat-containing cells are eliminated during Dictyostelium development

Triacylglycerol is a universal storage molecule for metabolic energy in living organisms. However, Dictyostelium amoebae, that have accumulated storage fat from added fatty acids do not progress through the starvation period preceding the development of the durable spore. Mutants deficient in genes of fat metabolism, such as fcsA, encoding a fatty acid activating enzyme, or dgat1 and dgat2, specifying proteins that synthesize triacylglycerol, strongly increase their chances to contribute to the spore fraction of the developing fruiting body, but lose the ability to produce storage fat efficiently. Dictyostelium seipin, an orthologue of a human protein that in patients causes the complete loss of adipose tissue when mutated, does not quantitatively affect fat storage in the amoeba. Dictyostelium seiP knockout mutants have lipid droplets that are enlarged in size but reduced in number. These mutants are as vulnerable as the wild type when exposed to fatty acids during their vegetative growth phase, and do not efficiently enter the spore head in Dictyostelium development.

Summary: In contrast to many living organisms, storage fat is not beneficial for Dictyostelium cells when progressing through starvation and subsequent development of a dormant stage.

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.

Expanding Francisella models: Pairing up the soil amoeba Dictyostelium with aquatic Francisella

The bacterial genus Francisella comprises highly pathogenic species that infect mammals, arthropods, fish and protists. Understanding virulence and host defense mechanisms of Francisella infection relies on multiple animal and cellular model systems. In this review, we want to summarize the most commonly used Francisella host model platforms and highlight novel, alternative model systems using aquatic Francisella species. Established mouse and macrophage models contributed extensively to our understanding of Francisella infection. However, murine and human cells display significant differences in their response to Francisella infection. The zebrafish and the amoeba Dictyostelium are well-established model systems for host-pathogen interactions and open up opportunities to investigate bacterial virulence and host defense. Comparisons between model systems using human and fish pathogenic Francisella species revealed shared virulence strategies and pathology between them. Hence, zebrafish and Dictyostelium might complement current model systems to find new vaccine candidates and contribute to our understanding of Francisella infection.

Hydrophilic interaction anion exchange for separation of multiply modified neutral and anionic Dictyostelium N-glycans

The unusual nature of the N-glycans of the cellular slime mould Dictyostelium discoideum has been revealed by a number of studies, primarily based on examination of radiolabeled glycopeptides but more recently also by MS. The complexity of the N-glycomes of even glycosylation mutants is compounded by the occurrence of anionic modifications, which also present an analytical challenge. In this study, we have employed hydrophilic interaction anion exchange (HIAX) HPLC in combination with MALDI-TOF MS/MS to explore the anionic N-glycome of the M31 (modA) strain, which lacks endoplasmic reticulum α-glucosidase II, an enzyme conserved in most eukaryotes including Homo sapiens. Prefractionation with HIAX chromatography enabled the identification of N-glycans with unusual oligo-α1,2-mannose extensions as well as others with up to four anionic modifications. Due to the use of hydrofluoric acid treatment, we were able to discriminate isobaric glycans differing in the presence of sulphate or phosphate on intersected structures as opposed to those carrying GlcNAc-phosphodiesters. The latter represent biosynthetic intermediates during the pathway leading to formation of the methylphosphorylated mannose epitope, which may have a similar function in intracellular targeting of hydrolases as the mannose-6-phosphate modification of lysosomal enzymes in mammals. In conclusion, HIAX in combination with MS is a highly sensitive approach for both fine separation and definition of neutral and anionic N-glycan structures.

dictyExpress: a web-based platform for sequence data management and analytics in Dictyostelium and beyond

BACKGROUND

Dictyostelium discoideum, a soil-dwelling social amoeba, is a model for the study of numerous biological processes. Research in the field has benefited mightily from the adoption of next-generation sequencing for genomics and transcriptomics. Dictyostelium biologists now face the widespread challenges of analyzing and exploring high dimensional data sets to generate hypotheses and discovering novel insights.

RESULTS

We present dictyExpress (2.0), a web application designed for exploratory analysis of gene expression data, as well as data from related experiments such as Chromatin Immunoprecipitation sequencing (ChIP-Seq). The application features visualization modules that include time course expression profiles, clustering, gene ontology enrichment analysis, differential expression analysis and comparison of experiments. All visualizations are interactive and interconnected, such that the selection of genes in one module propagates instantly to visualizations in other modules. dictyExpress currently stores the data from over 800 Dictyostelium experiments and is embedded within a general-purpose software framework for management of next-generation sequencing data. dictyExpress allows users to explore their data in a broader context by reciprocal linking with dictyBase-a repository of Dictyostelium genomic data. In addition, we introduce a companion application called GenBoard, an intuitive graphic user interface for data management and bioinformatics analysis.

CONCLUSIONS

dictyExpress and GenBoard enable broad adoption of next generation sequencing based inquiries by the Dictyostelium research community. Labs without the means to undertake deep sequencing projects can mine the data available to the public. The entire information flow, from raw sequence data to hypothesis testing, can be accomplished in an efficient workspace. The software framework is generalizable and represents a useful approach for any research community. To encourage more wide usage, the backend is open-source, available for extension and further development by bioinformaticians and data scientists.

Adaptations to High Salt in a Halophilic Protist: Differential Expression and Gene Acquisitions through Duplications and Gene Transfers

The capacity of halophiles to thrive in extreme hypersaline habitats derives partly from the tight regulation of ion homeostasis, the salt-dependent adjustment of plasma membrane fluidity, and the increased capability to manage oxidative stress. Halophilic bacteria, and archaea have been intensively studied, and substantial research has been conducted on halophilic fungi, and the green alga Dunaliella. By contrast, there have been very few investigations of halophiles that are phagotrophic protists, i.e., protozoa. To gather fundamental knowledge about salt adaptation in these organisms, we studied the transcriptome-level response of Halocafeteria seosinensis (Stramenopiles) grown under contrasting salinities. We provided further evolutionary context to our analysis by identifying genes that underwent recent duplications. Genes that were highly responsive to salinity variations were involved in stress response (e.g., chaperones), ion homeostasis (e.g., Na+/H+ transporter), metabolism and transport of lipids (e.g., sterol biosynthetic genes), carbohydrate metabolism (e.g., glycosidases), and signal transduction pathways (e.g., transcription factors). A significantly high proportion (43%) of duplicated genes were also differentially expressed, accentuating the importance of gene expansion in adaptation by H. seosinensis to high salt environments. Furthermore, we found two genes that were lateral acquisitions from bacteria, and were also highly up-regulated and highly expressed at high salt, suggesting that this evolutionary mechanism could also have facilitated adaptation to high salt. We propose that a transition toward high-salt adaptation in the ancestors of H. seosinensis required the acquisition of new genes via duplication, and some lateral gene transfers (LGTs), as well as the alteration of transcriptional programs, leading to increased stress resistance, proper establishment of ion gradients, and modification of cell structure properties like membrane fluidity.

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.

La Deletion from Mouse Brain Alters Pre-tRNA Metabolism and Accumulation of Pre-5.8S rRNA, with Neuron Death and Reactive Astrocytosis

Human La antigen (Sjögren's syndrome antigen B [SSB]) is an abundant multifunctional RNA-binding protein. In the nucleoplasm, La binds to and protects from 3' exonucleases, the ends of precursor tRNAs, and other transcripts synthesized by RNA polymerase III and facilitates their maturation, while a nucleolar isoform has been implicated in rRNA biogenesis by multiple independent lines of evidence. We showed previously that conditional La knockout (La cKO) from mouse cortex neurons results in defective tRNA processing, although the pathway(s) involved in neuronal loss thereafter was unknown. Here, we demonstrate that La is stably associated with a spliced pre-tRNA intermediate. Microscopic evidence of aberrant nuclear accumulation of 5.8S rRNA in La cKO is supported by a 10-fold increase in a pre-5.8S rRNA intermediate. To identify pathways involved in subsequent neurodegeneration and loss of brain mass in the cKO cortex, we employed mRNA sequencing (mRNA-Seq), immunohistochemistry, and other approaches. This revealed robust enrichment of immune and astrocyte reactivity in La cKO cortex. Immunohistochemistry, including temporal analyses, demonstrated neurodegeneration, followed by astrocyte invasion associated with immune response and decreasing cKO cortex size over time. Thus, deletion of La from postmitotic neurons results in defective pre-tRNA and pre-rRNA processing and progressive neurodegeneration with loss of cortical brain mass.

Whole-Genome Restriction Mapping by "Subhaploid"-Based RAD Sequencing: An Efficient and Flexible Approach for Physical Mapping and Genome Scaffolding

Assembly of complex genomes using short reads remains a major challenge, which usually yields highly fragmented assemblies. Generation of ultradense linkage maps is promising for anchoring such assemblies, but traditional linkage mapping methods are hindered by the infrequency and unevenness of meiotic recombination that limit attainable map resolution. Here we develop a sequencing-based "in vitro" linkage mapping approach (called RadMap), where chromosome breakage and segregation are realized by generating hundreds of "subhaploid" fosmid/bacterial-artificial-chromosome clone pools, and by restriction site-associated DNA sequencing of these clone pools to produce an ultradense whole-genome restriction map to facilitate genome scaffolding. A bootstrap-based minimum spanning tree algorithm is developed for grouping and ordering of genome-wide markers and is implemented in a user-friendly, integrated software package (AMMO). We perform extensive analyses to validate the power and accuracy of our approach in the model plant Arabidopsis thaliana and human. We also demonstrate the utility of RadMap for enhancing the contiguity of a variety of whole-genome shotgun assemblies generated using either short Illumina reads (300 bp) or long PacBio reads (6-14 kb), with up to 15-fold improvement of N50 (∼816 kb-3.7 Mb) and high scaffolding accuracy (98.1-98.5%). RadMap outperforms BioNano and Hi-C when input assembly is highly fragmented (contig N50 = 54 kb). RadMap can capture wide-range contiguity information and provide an efficient and flexible tool for high-resolution physical mapping and scaffolding of highly fragmented assemblies.

TRE5-A retrotransposition profiling reveals putative RNA polymerase III transcription complex binding sites on the Dictyostelium extrachromosomal rDNA element

The amoeba Dictyostelium discoideum has a haploid genome in which two thirds of the DNA encodes proteins. Consequently, the space available for selfish mobile elements to expand without excess damage to the host genome is limited. The non-long terminal repeat retrotransposon TRE5-A maintains an active population in the D. discoideum genome and apparently adapted to this gene-dense environment by targeting positions ~47 bp upstream of tRNA genes that are devoid of protein-coding regions. Because only ~24% of tRNA genes are associated with a TRE5-A element in the reference genome, we evaluated whether TRE5-A retrotransposition is limited to this subset of tRNA genes. We determined that a tagged TRE5-A element (TRE5-A^bsr) integrated at 384 of 405 tRNA genes, suggesting that expansion of the current natural TRE5-A population is not limited by the availability of targets. We further observed that TRE5-A^bsr targets the ribosomal 5S gene on the multicopy extrachromosomal DNA element that carries the ribosomal RNA genes, indicating that TRE5-A integration may extend to the entire RNA polymerase III (Pol III) transcriptome. We determined that both natural TRE5-A and cloned TRE5-A^bsr retrotranspose to locations on the extrachromosomal rDNA element that contain tRNA gene-typical A/B box promoter motifs without displaying any other tRNA gene context. Based on previous data suggesting that TRE5-A targets tRNA genes by locating Pol III transcription complexes, we propose that A/B box loci reflect Pol III transcription complex assembly sites that possess a function in the biology of the extrachromosomal rDNA element.

Bordetella bronchiseptica exploits the complex life cycle of Dictyostelium discoideum as an amplifying transmission vector

Multiple lines of evidence suggest that Bordetella species have a significant life stage outside of the mammalian respiratory tract that has yet to be defined. The Bordetella virulence gene (BvgAS) two-component system, a paradigm for a global virulence regulon, controls the expression of many “virulence factors” expressed in the Bvg positive (Bvg⁺) phase that are necessary for successful respiratory tract infection. A similarly large set of highly conserved genes are expressed under Bvg negative (Bvg^-) phase growth conditions; however, these appear to be primarily expressed outside of the host and are thus hypothesized to be important in an undefined extrahost reservoir. Here, we show that Bvg^- phase genes are involved in the ability of Bordetella bronchiseptica to grow and disseminate via the complex life cycle of the amoeba Dictyostelium discoideum. Unlike bacteria that serve as an amoeba food source, B. bronchiseptica evades amoeba predation, survives within the amoeba for extended periods of time, incorporates itself into the amoeba sori, and disseminates along with the amoeba. Remarkably, B. bronchiseptica continues to be transferred with the amoeba for months, through multiple life cycles of amoebae grown on the lawns of other bacteria, thus demonstrating a stable relationship that allows B. bronchiseptica to expand and disperse geographically via the D. discoideum life cycle. Furthermore, B. bronchiseptica within the sori can efficiently infect mice, indicating that amoebae may represent an environmental vector within which pathogenic bordetellae expand and disseminate to encounter new mammalian hosts. These data identify amoebae as potential environmental reservoirs as well as amplifying and disseminating vectors for B. bronchiseptica and reveal an important role for the Bvg^- phase in these interactions.

Bordetella species are infectious bacterial respiratory pathogens of a range of animals, including humans. Bordetellae grow in two phenotypically distinct “phases,” each specifically expressing a large set of genes. The Bvg⁺ phase is primarily associated with respiratory tract infection (RTI) and has been well studied. The similarly large set of genes specifically expressed in the Bvg^- phase is poorly understood but has been proposed to be involved in some undefined environmental niche. Recently, we reported the presence of Bordetella species in many soil and water sources, indicating extensive exposure to predators. Herein, we show that the Bvg^- phase mediates B. bronchiseptica interactions with the common soil predator D. discoideum. Surprisingly, the bacterium not only can evade predation but can propagate and disseminate via the complex developmental process of D. discoideum. After multiple passages and over a million-fold expansion in association with D. discoideum, B. bronchiseptica retained the ability to efficiently colonize mice. The conservation of the genes involved in these two distinct phases raises the possibility of potential environmental sources for the frequently unexplained outbreaks of diseases caused by this and other Bordetella species.

Functional and evolutionary characterization of Ohr proteins in eukaryotes reveals many active homologs among pathogenic fungi

Ohr and OsmC proteins comprise two subfamilies within a large group of proteins that display Cys-based, thiol dependent peroxidase activity. These proteins were previously thought to be restricted to prokaryotes, but we show here, using iterated sequence searches, that Ohr/OsmC homologs are also present in 217 species of eukaryotes with a massive presence in Fungi (186 species). Many of these eukaryotic Ohr proteins possess an N-terminal extension that is predicted to target them to mitochondria. We obtained recombinant proteins for four eukaryotic members of the Ohr/OsmC family and three of them displayed lipoyl peroxidase activity. Further functional and biochemical characterization of the Ohr homologs from the ascomycete fungus Mycosphaerella fijiensis Mf_1 (MfOhr), the causative agent of Black Sigatoka disease in banana plants, was pursued. Similarly to what has been observed for the bacterial proteins, we found that: (i) the peroxidase activity of MfOhr was supported by DTT or dihydrolipoamide (dithiols), but not by β-mercaptoethanol or GSH (monothiols), even in large excess; (ii) MfOhr displayed preference for organic hydroperoxides (CuOOH and tBOOH) over hydrogen peroxide; (iii) MfOhr presented extraordinary reactivity towards linoleic acid hydroperoxides (k=3.18 (±2.13)×10⁸ M⁻¹ s⁻¹). Both Cys⁸⁷ and Cys¹⁵⁴ were essential to the peroxidase activity, since single mutants for each Cys residue presented no activity and no formation of intramolecular disulfide bond upon treatment with hydroperoxides. The pK_a value of the Cys_p residue was determined as 5.7±0.1 by a monobromobimane alkylation method. Therefore, eukaryotic Ohr peroxidases share several biochemical features with prokaryotic orthologues and are preferentially located in mitochondria.

•

Ohr/OsmC proteins are also present in lower eukaryotic organisms.

•

While Ohr proteins are massively present among Fungi, OsmC proteins are restricted to the cellular slime molds.

•

Eukaryotic Ohr and OsmC present a thiol dependent peroxidase activity similar to the bacterial counterparts.

•

Most of these eukaryotic enzymes are predominantly present in mitochondria.

Regulation of nucleosome positioning by a CHD Type III chromatin remodeler and its relationship to developmental gene expression in Dictyostelium

Nucleosome placement and repositioning can direct transcription of individual genes; however, the precise interactions of these events are complex and largely unresolved at the whole-genome level. The Chromodomain-Helicase-DNA binding (CHD) Type III proteins are a subfamily of SWI2/SNF2 proteins that control nucleosome positioning and are associated with several complex human disorders, including CHARGE syndrome and autism. Type III CHDs are required for multicellular development of animals and Dictyostelium but are absent in plants and yeast. These CHDs can mediate nucleosome translocation in vitro, but their in vivo mechanism is unknown. Here, we use genome-wide analysis of nucleosome positioning and transcription profiling to investigate the in vivo relationship between nucleosome positioning and gene expression during development of wild-type (WT) Dictyostelium and mutant cells lacking ChdC, a Type III CHD protein ortholog. We demonstrate major nucleosome positional changes associated with developmental gene regulation in WT. Loss of chdC caused an increase of intragenic nucleosome spacing and misregulation of gene expression, affecting ∼50% of the genes that are repositioned during WT development. These analyses demonstrate active nucleosome repositioning during Dictyostelium multicellular development, establish an in vivo function of CHD Type III chromatin remodeling proteins in this process, and reveal the detailed relationship between nucleosome positioning and gene regulation, as cells transition between developmental states.

Evolving Notch polyQ tracts reveal possible solenoid interference elements

Polyglutamine (polyQ) tracts in regulatory proteins are extremely polymorphic. As functional elements under selection for length, triplet repeats are prone to DNA replication slippage and indel mutations. Many polyQ tracts are also embedded within intrinsically disordered domains, which are less constrained, fast evolving, and difficult to characterize. To identify structural principles underlying polyQ tracts in disordered regulatory domains, here I analyze deep evolution of metazoan Notch polyQ tracts, which can generate alleles causing developmental and neurogenic defects. I show that Notch features polyQ tract turnover that is restricted to a discrete number of conserved “polyQ insertion slots”. Notch polyQ insertion slots are: (i) identifiable by an amphipathic “slot leader” motif; (ii) conserved as an intact C-terminal array in a 1-to-1 relationship with the N-terminal solenoid-forming ankyrin repeats (ARs); and (iii) enriched in carboxamide residues (Q/N), whose sidechains feature dual hydrogen bond donor and acceptor atoms. Correspondingly, the terminal loop and β-strand of each AR feature conserved carboxamide residues, which would be susceptible to folding interference by hydrogen bonding with residues outside the ARs. I thus suggest that Notch polyQ insertion slots constitute an array of AR interference elements (ARIEs). Notch ARIEs would dynamically compete with the delicate serial folding induced by adjacent ARs. Huntingtin, which harbors solenoid-forming HEAT repeats, also possesses a similar number of polyQ insertion slots. These results suggest that intrinsically disordered interference arrays featuring carboxamide and polyQ enrichment may constitute coupled proteodynamic modulators of solenoids.

Factors That Shape Eukaryotic tRNAomes: Processing, Modification and Anticodon-Codon Use

Transfer RNAs (tRNAs) contain sequence diversity beyond their anticodons and the large variety of nucleotide modifications found in all kingdoms of life. Some modifications stabilize structure and fit in the ribosome whereas those to the anticodon loop modulate messenger RNA (mRNA) decoding activity more directly. The identities of tRNAs with some universal anticodon loop modifications vary among distant and parallel species, likely to accommodate fine tuning for their translation systems. This plasticity in positions 34 (wobble) and 37 is reflected in codon use bias. Here, we review convergent evidence that suggest that expansion of the eukaryotic tRNAome was supported by its dedicated RNA polymerase III transcription system and coupling to the precursor-tRNA chaperone, La protein. We also review aspects of eukaryotic tRNAome evolution involving G34/A34 anticodon-sparing, relation to A34 modification to inosine, biased codon use and regulatory information in the redundancy (synonymous) component of the genetic code. We then review interdependent anticodon loop modifications involving position 37 in eukaryotes. This includes the eukaryote-specific tRNA modification, 3-methylcytidine-32 (m³C₃₂) and the responsible gene, TRM140 and homologs which were duplicated and subspecialized for isoacceptor-specific substrates and dependence on i⁶A₃₇ or t⁶A₃₇. The genetics of tRNA function is relevant to health directly and as disease modifiers.

Site-specific ADP-ribosylation of histone H2B in response to DNA double strand breaks

ADP-ribosyltransferases (ARTs) modify proteins with single units or polymers of ADP-ribose to regulate DNA repair. However, the substrates for these enzymes are ill-defined. For example, although histones are modified by ARTs, the sites on these proteins ADP-ribosylated following DNA damage and the ARTs that catalyse these events are unknown. This, in part, is due to the lack of a eukaryotic model that contains ARTs, in addition to histone genes that can be manipulated to assess ADP-ribosylation events in vivo. Here we exploit the model Dictyostelium to identify site-specific histone ADP-ribosylation events in vivo and define the ARTs that mediate these modifications. Dictyostelium histones are modified in response to DNA double strand breaks (DSBs) in vivo by the ARTs Adprt1a and Adprt2. Adprt1a is a mono-ART that modifies H2BE18 in vitro, although disruption of this site allows ADP-ribosylation at H2BE19. Although redundancy between H2BE18 and H2BE19 ADP-ribosylation is also apparent following DSBs in vivo, by generating a strain with mutations at E18/E19 in the h2b locus we demonstrate these are the principal sites modified by Adprt1a/Adprt2. This identifies DNA damage induced histone mono-ADP-ribosylation sites by specific ARTs in vivo, providing a unique platform to assess how histone ADP-ribosylation regulates DNA repair.

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.

Functional similarities between the dictyostelium protein AprA and the human protein dipeptidyl-peptidase IV

Autocrine proliferation repressor protein A (AprA) is a protein secreted by Dictyostelium discoideum cells. Although there is very little sequence similarity between AprA and any human protein, AprA has a predicted structural similarity to the human protein dipeptidyl peptidase IV (DPPIV). AprA is a chemorepellent for Dictyostelium cells, and DPPIV is a chemorepellent for neutrophils. This led us to investigate if AprA and DPPIV have additional functional similarities. We find that like AprA, DPPIV is a chemorepellent for, and inhibits the proliferation of, D. discoideum cells, and that AprA binds some DPPIV binding partners such as fibronectin. Conversely, rAprA has DPPIV-like protease activity. These results indicate a functional similarity between two eukaryotic chemorepellent proteins with very little sequence similarity, and emphasize the usefulness of using a predicted protein structure to search a protein structure database, in addition to searching for proteins with similar sequences.

Comparative Proteomics of Purified Pathogen Vacuoles Correlates Intracellular Replication of Legionella pneumophila with the Small GTPase Ras-related protein 1 (Rap1)

Legionella pneumophila is an opportunistic bacterial pathogen that causes a severe lung infection termed "Legionnaires' disease." The pathogen replicates in environmental protozoa as well as in macrophages within a unique membrane-bound compartment, the Legionella-containing-vacuole (LCV). LCV formation requires the bacterial Icm/Dot type IV secretion system, which translocates ca. 300 "effector proteins" into host cells, where they target distinct host factors. The L. pneumophila "pentuple" mutant (Δpentuple) lacks 5 gene clusters (31% of the effector proteins) and replicates in macrophages but not in Dictyostelium discoideum amoeba. To elucidate the host factors defining a replication-permissive compartment, we compare here the proteomes of intact LCVs isolated from D. discoideum or macrophages infected with Δpentuple or the parental strain Lp02. This analysis revealed that the majority of host proteins are shared in D. discoideum or macrophage LCVs containing the mutant or the parental strain, respectively, whereas some proteins preferentially localize to distinct LCVs. The small GTPase Rap1 was identified on D. discoideum LCVs containing strain Lp02 but not the Δpentuple mutant and on macrophage LCVs containing either strain. The localization pattern of active Rap1 on D. discoideum or macrophage LCVs was confirmed by fluorescence microscopy and imaging flow cytometry, and the depletion of Rap1 by RNA interference significantly reduced the intracellular growth of L. pneumophila Thus, comparative proteomics identified Rap1 as a novel LCV host component implicated in intracellular replication of L. pneumophila.

The Long Noncoding RNA Transcriptome of Dictyostelium discoideum Development

Dictyostelium discoideum live in the soil as single cells, engulfing bacteria and growing vegetatively. Upon starvation, tens of thousands of amoebae enter a developmental program that includes aggregation, multicellular differentiation, and sporulation. Major shifts across the protein-coding transcriptome accompany these developmental changes. However, no study has presented a global survey of long noncoding RNAs (ncRNAs) in D. discoideum To characterize the antisense and long intergenic noncoding RNA (lncRNA) transcriptome, we analyzed previously published developmental time course samples using an RNA-sequencing (RNA-seq) library preparation method that selectively depletes ribosomal RNAs (rRNAs). We detected the accumulation of transcripts for 9833 protein-coding messenger RNAs (mRNAs), 621 lncRNAs, and 162 putative antisense RNAs (asRNAs). The noncoding RNAs were interspersed throughout the genome, and were distinct in expression level, length, and nucleotide composition. The noncoding transcriptome displayed a temporal profile similar to the coding transcriptome, with stages of gradual change interspersed with larger leaps. The transcription profiles of some noncoding RNAs were strongly correlated with known differentially expressed coding RNAs, hinting at a functional role for these molecules during development. Examining the mitochondrial transcriptome, we modeled two novel antisense transcripts. We applied yet another ribosomal depletion method to a subset of the samples to better retain transfer RNA (tRNA) transcripts. We observed polymorphisms in tRNA anticodons that suggested a post-transcriptional means by which D. discoideum compensates for codons missing in the genomic complement of tRNAs. We concluded that the prevalence and characteristics of long ncRNAs indicate that these molecules are relevant to the progression of molecular and cellular phenotypes during development.

Context characterization of amino acid homorepeats using evolution, position, and order

Amino acid repeats, or homorepeats, are low complexity protein motifs consisting of tandem repetitions of a single amino acid. Their presence and relative number vary in different proteomes, and some studies have tried to address this variation, proteome by proteome. In this work, we present a full characterization of amino acid homorepeats across evolution. We studied the presence and differential usage of each possible homorepeat in proteomes from various taxonomic groups, using clusters of very similar proteins to eliminate redundancy. The position of each amino acid repeat within proteins, and the order of co-occurring amino acid repeats were also addressed. As a result, we present evidence about the unevenly evolution of homorepeats, as well as the functional implications of their relative position in proteins. We discuss some of these cases in their taxonomic context. Collectively, our results show evolutionary and positional signals that suggest that homorepeats have biological function, likely creating unspecific protein interactions or modulating specific interactions in a context dependent manner. In conclusion, our work supports the functional importance of homorepeats and establishes a basis for the study of other low complexity repeats. Proteins 2017; 85:709-719. © 2016 Wiley Periodicals, Inc.

Phylogenetic evidence for the modular evolution of metazoan signalling pathways

Communication among cells was paramount to the evolutionary increase in cell type diversity and, ultimately, the origin of large body size. Across the diversity of Metazoa, there are only few conserved cell signalling pathways known to orchestrate the complex cell and tissue interactions regulating development; thus, modification to these few pathways has been responsible for generating diversity during the evolution of animals. Here, we summarize evidence for the origin and putative function of the intracellular, membrane-bound and secreted components of seven metazoan cell signalling pathways with a special focus on early branching metazoans (ctenophores, poriferans, placozoans and cnidarians) and basal unikonts (amoebozoans, fungi, filastereans and choanoflagellates). We highlight the modular incorporation of intra- and extracellular components in each signalling pathway and suggest that increases in the complexity of the extracellular matrix may have further promoted the modulation of cell signalling during metazoan evolution. Most importantly, this updated view of metazoan signalling pathways highlights the need for explicit study of canonical signalling pathway components in taxa that do not operate a complete signalling pathway. Studies like these are critical for developing a deeper understanding of the evolution of cell signalling.This article is part of the themed issue 'Evo-devo in the genomics era, and the origins of morphological diversity'.

Effects of deletion of the receptor CrlA on Dictyostelium aggregation and MPBD-mediated responses are strain dependent and not evident in strain Ax2

The polyketide MPBD (4-methyl-5-pentylbenzene-1, 3-diol) is produced by the polyketide synthase SteelyA (StlA) in Dictyostelium discoideum. MPBD is required for appropriate expression of cAMP signalling genes involved in cell aggregation and additionally induces the spore maturation at the fruiting body stage. The MPBD signalling pathway for regulation of cell aggregation is unknown, but MPBD effects on sporulation were reported to be mediated by the G-protein coupled receptor CrlA in D. discoideum KAx3. In this study, we deleted the crlA gene from the same parental strain (Ax2) that was used to generate the MPBD-less mutant. We found that unlike the MPBD-less mutant, Ax2-derived crlA- mutants exhibited normal cell aggregation, indicating that in Ax2 MPBD effects on early development do not require CrlA. We also found that the Ax2/crlA- mutant formed normal spores in fruiting bodies. When transformed with PkaC, both Ax2 and Ax2/crlA- similarly responded to MPBD in vitro with spore encapsulation. Our data make it doubtful that CrlA acts as the receptor for MPBD signalling during the development of D. discoideum Ax2.

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

Background

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

Results

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

Conclusions

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

Electronic supplementary material

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

A polychromatic 'greenbeard' locus determines patterns of cooperation in a social amoeba

Cheaters disrupt cooperation by reaping the benefits without paying their fair share of associated costs. Cheater impact can be diminished if cooperators display a tag (‘greenbeard') and recognise and preferentially direct cooperation towards other tag carriers. Despite its popular appeal, the feasibility of such greenbeards has been questioned because the complex patterns of partner-specific cooperative behaviours seen in nature require greenbeards to come in different colours. Here we show that a locus (‘Tgr') of a social amoeba represents a polychromatic greenbeard. Patterns of natural Tgr locus sequence polymorphisms predict partner-specific patterns of cooperation by underlying variation in partner-specific protein–protein binding strength and recognition specificity. Finally, Tgr locus polymorphisms increase fitness because they help avoid potential costs of cooperating with incompatible partners. These results suggest that a polychromatic greenbeard can provide a key mechanism for the evolutionary maintenance of cooperation.

Cooperation can be stabilized against exploitation if cooperators can reliably recognize each other. Here, Gruenheit and colleagues show that different alleles of the Tgr locus of the social amoeba Dictyostelium discoideum underlie the ability of different strains to recognize and cooperate with socially compatible individuals.

Adenylate cyclase A acting on PKA mediates induction of stalk formation by cyclic diguanylate at the Dictyostelium organizer

Coordination of cell movement with cell differentiation is a major feat of embryonic development. The Dictyostelium stalk always forms at the organizing tip, by a mechanism that is not understood. We previously reported that cyclic diguanylate (c-di-GMP), synthesized by diguanylate cyclase A (DgcA), induces stalk formation. Here we used transcriptional profiling of dgca- structures to identify target genes for c-di-GMP, and used these genes to investigate the c-di-GMP signal transduction pathway. We found that knockdown of cAMP-dependent protein kinase (PKA) activity in prestalk cells reduced stalk gene induction by c-di-GMP, whereas PKA activation bypassed the c-di-GMP requirement for stalk gene expression. c-di-GMP caused a persistent increase in cAMP, which still occurred in mutants lacking the adenylate cyclases ACG or ACR, or the cAMP phosphodiesterase RegA. However, both inhibition of adenylate cyclase A (ACA) with SQ22536 and incubation of a temperature-sensitive ACA mutant at the restrictive temperature prevented c-di-GMP-induced cAMP synthesis as well as c-di-GMP-induced stalk gene transcription. ACA produces the cAMP pulses that coordinate Dictyostelium morphogenetic cell movement and is highly expressed at the organizing tip. The stalk-less dgca- mutant regained its stalk by expression of a light-activated adenylate cyclase from the ACA promoter and exposure to light, indicating that cAMP is also the intermediate for c-di-GMP in vivo. Our data show that the more widely expressed DgcA activates tip-expressed ACA, which then acts on PKA to induce stalk genes. These results explain why stalk formation in Dictyostelia always initiates at the site of the morphogenetic organizer.

Genome-wide identification and comparative expression analysis of LEA genes in watermelon and melon genomes

Late embryogenesis abundant (LEA) proteins are large and diverse group of polypeptides which were first identified during seed dehydration and then in vegetative plant tissues during different stress responses. Now, gene family members of LEA proteins have been detected in various organisms. However, there is no report for this protein family in watermelon and melon until this study. A total of 73 LEA genes from watermelon (ClLEA) and 61 LEA genes from melon (CmLEA) were identified in this comprehensive study. They were classified into four and three distinct clusters in watermelon and melon, respectively. There was a correlation between gene structure and motif composition among each LEA groups. Segmental duplication played an important role for LEA gene expansion in watermelon. Maximum gene ontology of LEA genes was observed with poplar LEA genes. For evaluation of tissue specific expression patterns of ClLEA and CmLEA genes, publicly available RNA-seq data were analyzed. The expression analysis of selected LEA genes in root and leaf tissues of drought-stressed watermelon and melon were examined using qRT-PCR. Among them, ClLEA-12-17-46 genes were quickly induced after drought application. Therefore, they might be considered as early response genes for water limitation conditions in watermelon. In addition, CmLEA-42-43 genes were found to be up-regulated in both tissues of melon under drought stress. Our results can open up new frontiers about understanding of functions of these important family members under normal developmental stages and stress conditions by bioinformatics and transcriptomic approaches.

EhNPC1 and EhNPC2 Proteins Participate in Trafficking of Exogenous Cholesterol in Entamoeba histolytica Trophozoites: Relevance for Phagocytosis

Entamoeba histolytica, the highly phagocytic protozoan causative of human amoebiasis lacks the machinery to synthesize cholesterol. Here, we investigated the presence of NPC1 and NPC2 proteins in this parasite, which are involved in cholesterol trafficking in mammals. Bioinformatics analysis revealed one Ehnpc1 and two Ehnpc2 genes. EhNPC1 appeared as a transmembrane protein and both EhNPC2 as peripheral membrane proteins. Molecular docking predicted that EhNPC1 and EhNPC2 bind cholesterol and interact with each other. Genes and proteins were identified in trophozoites. Serum pulse-chase and confocal microscopy assays unveiled that after trophozoites sensed the cholesterol source, EhNPC1 and EhNPC2 were organized around the plasma membrane in a punctuated pattern. Vesicles emerged and increased in number and size and some appeared full of cholesterol with EhNPC1 or EhNPC2 facing the extracellular space. Both proteins, but mostly EhNPC2, were found out of the cell associated with cholesterol. EhNPC1 and cholesterol formed networks from the plasma membrane to the nucleus. EhNPC2 appeared in erythrocytes that were being ingested by trophozoites, co-localizing with cholesterol of erythrocytes, whereas EhNPC1 surrounded the phagocytic cup. EhNPC1 and EhNPC2 co-localized with EhSERCA in the endoplasmic reticulum and with lysobisphosphatidic acid and EhADH (an Alix protein) in phagolysosomes. Immunoprecipitation assays confirmed the EhNPC1 and EhNPC2 association with cholesterol, EhRab7A and EhADH. Serum starved and blockage of cholesterol trafficking caused a low rate of phagocytosis and incapability of trophozoites to produce damage in the mouse colon. Ehnpc1 and Ehnpc2 knockdown provoked in trophozoites a lower intracellular cholesterol concentration and a diminished rate of phagocytosis; and Ehnpc1 silencing also produced a decrease of trophozoites movement. Trafficking of EhNPC1 and EhNPC2 during cholesterol uptake and phagocytosis as well as their association with molecules involved in endocytosis strongly suggest that these proteins play a key role in cholesterol uptake.

NPC1 and NPC2 proteins are involved in cholesterol trafficking in mammals. Using different approaches, we have detected the orthologues EhNPC1 and EhNPC2 proteins in Entamoeba histolytica. Trophozoites are particularly rich in membranes and vacuoles, but they do not possess the machinery to synthetize cholesterol. Thus, they are completely dependent on molecules able to “fish” cholesterol from the medium. The relevance of our findings lies in the fact that cholesterol is fundamental for endocytosis and motility; and, phagocytosis is an important nutritional and virulence factor for E. histolytica. In silico and experimental strategies, using U18666A to arrest cholesterol trafficking, as well as, knockdown mutants, showed that EhNPC1 and EhNPC2 participate in cholesterol uptake and trafficking in this parasite. They are secreted by trophozoites and directly involved in erythrophagocytosis and motility. Our findings revealed E. histolytica as one of the first protozoa in which these proteins are being characterized. Moreover, E. histolytica provides an excellent and less complicated model to elucidate the intricate event of cholesterol trafficking in eukaryotic cells. The relevance of cholesterol transport for the parasite virulence and the involvement of EhNPC1 and EhNPC2 in this process, make these proteins promising targets for therapy strategies development against the parasite.

The C-Terminal SynMuv/DdDUF926 Domain Regulates the Function of the N-Terminal Domain of DdNKAP

NKAP (NF-κB activating protein) is a highly conserved SR (serine/arginine-rich) protein involved in transcriptional control and splicing in mammals. We identified DdNKAP, the Dictyostelium discoideum ortholog of mammalian NKAP, as interacting partner of the nuclear envelope protein SUN-1. DdNKAP harbors a number of basic RDR/RDRS repeats in its N-terminal domain and the SynMuv/DUF926 domain at its C-terminus. We describe a novel and direct interaction between DdNKAP and Prp19 (Pre mRNA processing factor 19) which might be relevant for the observed DdNKAP ubiquitination. Genome wide analysis using cross-linking immunoprecipitation-high-throughput sequencing (CLIP-seq) revealed DdNKAP association with intergenic regions, exons, introns and non-coding RNAs. Ectopic expression of DdNKAP and its domains affects several developmental aspects like stream formation, aggregation, and chemotaxis. We conclude that DdNKAP is a multifunctional protein, which might influence Dictyostelium development through its interaction with RNA and RNA binding proteins. Mutants overexpressing full length DdNKAP and the N-terminal domain alone (DdN-NKAP) showed opposite phenotypes in development and opposite expression profiles of several genes and rRNAs. The observed interaction between DdN-NKAP and the DdDUF926 domain indicates that the DdDUF926 domain acts as negative regulator of the N-terminus.

Independent evolution of genomic characters during major metazoan transitions

Metazoan evolution encompasses a vast evolutionary time scale spanning over 600 million years. Our ability to infer ancestral metazoan characters, both morphological and functional, is limited by our understanding of the nature and evolutionary dynamics of the underlying regulatory networks. Increasing coverage of metazoan genomes enables us to identify the evolutionary changes of the relevant genomic characters such as the loss or gain of coding sequences, gene duplications, micro- and macro-synteny, and non-coding element evolution in different lineages. In this review we describe recent advances in our understanding of ancestral metazoan coding and non-coding features, as deduced from genomic comparisons. Some genomic changes such as innovations in gene and linkage content occur at different rates across metazoan clades, suggesting some level of independence among genomic characters. While their contribution to biological innovation remains largely unclear, we review recent literature about certain genomic changes that do correlate with changes to specific developmental pathways and metazoan innovations. In particular, we discuss the origins of the recently described pharyngeal cluster which is conserved across deuterostome genomes, and highlight different genomic features that have contributed to the evolution of this group. We also assess our current capacity to infer ancestral metazoan states from gene models and comparative genomics tools and elaborate on the future directions of metazoan comparative genomics relevant to evo-devo studies.

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.

Terpene Cyclases from Social Amoebae

Genome sequences of social amoebae reveal the presence of terpene cyclases (TCs) in these organisms. Two TCs from Dictyostelium discoideum converted farnesyl diphosphate into (2S,3R,6S,9S)-(-)-protoillud-7-ene and (3S)-(+)-asterisca-2(9),6-diene. The enzyme mechanisms and EI-MS fragmentations of the products were studied by labeling experiments.

A Diaphanous-related formin links Ras signaling directly to actin assembly in macropinocytosis and phagocytosis

Phagocytosis and macropinocytosis are Ras-regulated and actin-driven processes that depend on the dynamic rearrangements of the plasma membrane that protrudes and internalizes extracellular material by cup-shaped structures. However, the regulatory mechanisms underlying actin assembly in large-scale endocytosis remain elusive. Here, we show that the Diaphanous-related formin G (ForG) from the professional phagocyte Dictyostelium discoideum localizes to endocytic cups. Biochemical analyses revealed that ForG is a rather weak nucleator but efficiently elongates actin filaments in the presence of profilin. Notably, genetic inactivation of ForG is associated with a strongly impaired endocytosis and a markedly diminished F-actin content at the base of the cups. By contrast, ablation of the Arp2/3 (actin-related protein-2/3) complex activator SCAR (suppressor of cAMP receptor) diminishes F-actin mainly at the cup rim, being consistent with its known localization. These data therefore suggest that ForG acts as an actin polymerase of Arp2/3-nucleated filaments to allow for efficient membrane expansion and engulfment of extracellular material. Finally, we show that ForG is directly regulated in large-scale endocytosis by RasB and RasG, which are highly related to the human proto-oncogene KRas.

Dictyostelium discoideum as a Novel Host System to Study the Interaction between Phagocytes and Yeasts

The social amoeba Dictyostelium discoideum is a well-established model organism to study the interaction between bacteria and phagocytes. In contrast, research using D. discoideum as a host model for fungi is rare. We describe a comprehensive study, which uses D. discoideum as a host model system to investigate the interaction with apathogenic (Saccharomyces cerevisiae) and pathogenic (Candida sp.) yeast. We show that Dictyostelium can be co-cultivated with yeasts on solid media, offering a convenient test to study the interaction between fungi and phagocytes. We demonstrate that a number of D. discoideum mutants increase (atg1^-, kil1^-, kil2^-) or decrease (atg6^-) the ability of the amoebae to predate yeast cells. On the yeast side, growth characteristics, reduced phagocytosis rate, as well as known virulence factors of C. albicans (EFG1, CPH1, HGC1, ICL1) contribute to the resistance of yeast cells against predation by the amoebae. Investigating haploid C. albicans strains, we suggest using the amoebae plate test for screening purposes after random mutagenesis. Finally, we discuss the potential of our adapted amoebae plate test to use D. discoideum for risk assessment of yeast strains.

Huntingtin-associated protein 1: Eutherian adaptation from a TRAK-like protein, conserved gene promoter elements, and localization in the human intestine

Background

Huntingtin-associated Protein 1 (HAP1) is expressed in neurons and endocrine cells, and is critical for postnatal survival in mice. HAP1 shares a conserved “HAP1_N” domain with TRAfficking Kinesin proteins TRAK1 and TRAK2 (vertebrate), Milton (Drosophila) and T27A3.1 (C. elegans). HAP1, TRAK1 and TRAK2 have a degree of common function, particularly regarding intracellular receptor trafficking. However, TRAK1, TRAK2 and Milton (which have a “Milt/TRAK” domain that is absent in human and rodent HAP1) differ in function to HAP1 in that they are mitochondrial transport proteins, while HAP1 has emerging roles in starvation response. We have investigated HAP1 function by examining its evolution, and upstream gene promoter sequences. We performed phylogenetic analyses of the HAP1_N domain family of proteins, incorporating HAP1 orthologues (identified by genomic synteny) from 5 vertebrate classes, and also searched the Dictyostelium proteome for a common ancestor. Computational analyses of mammalian HAP1 gene promoters were performed to identify phylogenetically conserved regulatory motifs.

Results

We found that as recently as marsupials, HAP1 contained a Milt/TRAK domain and was more similar to TRAK1 and TRAK2 than to eutherian HAP1. The Milt/TRAK domain likely arose post multicellularity, as it was absent in the Dictyostelium proteome. It was lost from HAP1 in the eutherian lineage, and also from T27A3.1 in C. elegans. The HAP1 promoter from human, mouse, rat, rabbit, horse, dog, Tasmanian devil and opossum contained common sites for transcription factors involved in cell cycle, growth, differentiation, and stress response. A conserved arrangement of regulatory elements was identified, including sites for caudal-related homeobox transcription factors (CDX1 and CDX2), and myc-associated factor X (MAX) in the region of the TATA box. CDX1 and CDX2 are intestine-enriched factors, prompting investigation of HAP1 protein expression in the human duodenum. HAP1 was localized to singly dispersed mucosal cells, including a subset of serotonin-positive enterochromaffin cells.

Conclusion

We have identified eutherian HAP1 as an evolutionarily recent adaptation of a vertebrate TRAK protein-like ancestor, and found conserved CDX1/CDX2 and MAX transcription factor binding sites near the TATA box in mammalian HAP1 gene promoters. We also demonstrated that HAP1 is expressed in endocrine cells of the human gut.

Electronic supplementary material

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

Evolution of bacterial recombinase A (recA) in eukaryotes explained by addition of genomic data of key microbial lineages

Recombinase enzymes promote DNA repair by homologous recombination. The genes that encode them are ancestral to life, occurring in all known dominions: viruses, Eubacteria, Archaea and Eukaryota. Bacterial recombinases are also present in viruses and eukaryotic groups (supergroups), presumably via ancestral events of lateral gene transfer. The eukaryotic recA genes have two distinct origins (mitochondrial and plastidial), whose acquisition by eukaryotes was possible via primary (bacteria-eukaryote) and/or secondary (eukaryote-eukaryote) endosymbiotic gene transfers (EGTs). Here we present a comprehensive phylogenetic analysis of the recA genealogy, with substantially increased taxonomic sampling in the bacteria, viruses, eukaryotes and a special focus on the key eukaryotic supergroup Amoebozoa, earlier represented only by Dictyostelium We demonstrate that several major eukaryotic lineages have lost the bacterial recombinases (including Opisthokonta and Excavata), whereas others have retained them (Amoebozoa, Archaeplastida and the SAR-supergroups). When absent, the bacterial recA homologues may have been lost entirely (secondary loss of canonical mitochondria) or replaced by other eukaryotic recombinases. RecA proteins have a transit peptide for organellar import, where they act. The reconstruction of the RecA phylogeny with its EGT events presented here retells the intertwined evolutionary history of eukaryotes and bacteria, while further illuminating the events of endosymbiosis in eukaryotes by expanding the collection of widespread genes that provide insight to this deep history.

A pre-metazoan origin of the CRK gene family and co-opted signaling network

CRK and CRKL adapter proteins play essential roles in development and cancer through their SRC homology 2 and 3 (SH2 and SH3) domains. To gain insight into the origin of their shared functions, we have investigated their evolutionary history. We propose a term, crk/crkl ancestral (crka), for orthologs in invertebrates before the divergence of CRK and CRKL in the vertebrate ancestor. We have isolated two orthologs expressed in the choanoflagellate Monosiga brevicollis, a unicellular relative to the metazoans. Consistent with its highly-conserved three-dimensional structure, the SH2 domain of M. brevicollis crka1 can bind to the mammalian CRK/CRKL SH2 binding consensus phospho-YxxP, and to the SRC substrate/focal adhesion protein BCAR1 (p130^CAS) in the presence of activated SRC. These results demonstrate an ancient origin of the CRK/CRKL SH2-target recognition specificity. Although BCAR1 orthologs exist only in metazoans as identified by an N-terminal SH3 domain, YxxP motifs, and a C-terminal FAT-like domain, some pre-metazoan transmembrane proteins include several YxxP repeats in their cytosolic region, suggesting that they are remotely related to the BCAR1 substrate domain. Since the tyrosine kinase SRC also has a pre-metazoan origin, co-option of BCAR1-related sequences may have rewired the crka-dependent network to mediate adhesion signals in the metazoan ancestor.