Whether coevolution determines distribution patterns of facultative symbionts among individuals of a host species is an open question. We used the Dictyostelium discoideum - Paraburkholderia system to investigate the role of host fitness and symbiont transmission in co-occurrence patterns of naturally occurring (native) host-symbiont strain pairs. P. bonniea is able to persistently infect D. discoideum in the lab, yet few wild D. discoideum isolates natively host this symbiont. If native host-symbiont combinations lead to higher fitness for host and symbiont compared to nonnative combinations, we assumed a stronger chance that coevolution had led to their co-occurrence. We approached the factors that influence the fitness consequences of host-symbiont interactions in terms of host defense mechanisms and symbiont virulence mechanisms. Given a specific combination of host and symbiont, variable host resistance and symbiont proliferation lead to variation in symbiont density, while variable host tolerance and symbiont benevolence lead to variation in host fitness relative to symbiont density. We used both vegetative and social stages of the host life cycle to approximate interactions between D. discoideum and P. bonniea , when hosts form social groups and symbiont horizontal transmission is facilitated by these social groups. We found no difference in host defense mechanisms between native and nonnative hosts of P. bonniea , but P. bonniea possesses variable phenotypes in proliferation and benevolence. More importantly, variation in symbiont benevolence led to increased horizontal transmission of benevolent strains compared to malevolent strains during the social stage of the host life cycle.

The "Amoeboid Predator-Fungal Animal Virulence Hypothesis" posits that interactions with environmental phagocytes shape the evolution of virulence traits in fungal pathogens. In this hypothesis, selection to avoid predation by amoeba inadvertently selects for traits that contribute to fungal escape from phagocytic immune cells. Here, we investigate this hypothesis in the human fungal pathogens Cryptococcus neoformans and Cryptococcus deneoformans. Applying quantitative trait locus (QTL) mapping and comparative genomics, we discovered a cross-species QTL region that is responsible for variation in resistance to amoeba predation. In C. neoformans, this same QTL was found to have pleiotropic effects on melanization, an established virulence factor. Through fine mapping and population genomic comparisons, we identified the gene encoding the transcription factor Bzp4 that underlies this pleiotropic QTL and we show that decreased expression of this gene reduces melanization and increases susceptibility to amoeba predation. Despite the joint effects of BZP4 on amoeba resistance and melanin production, we find no relationship between BZP4 genotype and escape from macrophages or virulence in murine models of disease. Our findings provide new perspectives on how microbial ecology shapes the genetic architecture of fungal virulence, and suggests the need for more nuanced models for the evolution of pathogenesis that account for the complexities of both microbe-microbe and microbe-host interactions.

A small subset of bacteria in soil interact directly with eukaryotes. Which ones do so can reveal what is important to a eukaryote and how eukaryote defenses might be breached. Soil amoebae are simple eukaryotic organisms and as such could be particularly good for understanding how eukaryote microbiomes originate and are maintained. One such amoeba, Dictyostelium discoideum, has both permanent and temporary associations with bacteria. Here we focus on culturable bacterial associates in order to interrogate their relationship with D. discoideum. To do this, we isolated over 250 D. discoideum fruiting body samples from soil and deer feces at Mountain Lake Biological Station. In one-third of the wild D. discoideum we tested, one to six bacterial species were found per fruiting body sorus (spore mass) for a total of 174 bacterial isolates. The remaining two-thirds of D. discoideum fruiting body samples did not contain culturable bacteria, as is thought to be the norm. A majority (71.4%) of the unique bacterial haplotypes are in Proteobacteria. The rest are in either Actinobacteria, Bacteriodetes, or Firmicutes. The highest bacterial diversity was found in D. discoideum fruiting bodies originating from deer feces (27 OTUs), greater than either of those originating in shallow (11 OTUs) or in deep soil (4 OTUs). Rarefaction curves and the Chao1 estimator for species richness indicated the diversity in any substrate was not fully sampled, but for soil it came close. A majority of the D. discoideum-associated bacteria were edible by D. discoideum and supported its growth (75.2% for feces and 81.8% for soil habitats). However, we found several bacteria genera were able to evade phagocytosis and persist in D. discoideum cells through one or more social cycles. This study focuses not on the entire D. discoideum microbiome, but on the culturable subset of bacteria that have important eukaryote interactions as prey, symbionts, or pathogens. These eukaryote and bacteria interactions may provide fertile ground for investigations of bacteria using amoebas to gain an initial foothold in eukaryotes and of the origins of symbiosis and simple microbiomes.

Background and Objectives: Legionella is one of the most important water-related pathogens. Inside the water supply systems and the biofilms, Legionella interact with other bacteria and free-living amoeba (FLA). Several amoebas may serve as hosts for bacteria in aquatic systems. This study aimed to investigate the co-occurrence of Legionella spp. and FLA in drinking water supply systems. Materials and Methods: A total of 268 water samples were collected from apartment buildings, hotels, and public buildings. Detection of Legionella spp. was performed in accordance with ISO 11731:2017 standard. Three different polymerase chain reaction (PCR) protocols were used to identify FLA. Results: Occurrence of Legionella varied from an average of 12.5% in cold water samples with the most frequent occurrence observed in hot water, in areas receiving untreated groundwater, where 54.0% of the samples were Legionella positive. The occurrence of FLA was significantly higher. On average, 77.2% of samples contained at least one genus of FLA and, depending on the type of sample, the occurrence of FLA could reach 95%. In the samples collected during the study, Legionella was always isolated along with FLA, no samples containing Legionella in the absence of FLA were observed. Conclusions: The data obtained in our study can help to focus on the extensive distribution, close interaction, and long-term persistence of Legionella and FLA. Lack of Legionella risk management plans and control procedures may promote further spread of Legionella in water supply systems. In addition, the high incidence of Legionella-related FLA suggests that traditional monitoring methods may not be sufficient for Legionella control.

The evolution of eukaryotic life was predicated on the development of organelles such as mitochondria and plastids. During this complex process of organellogenesis, the host cell and the engulfed prokaryote became genetically codependent, with the integration of genes from the endosymbiont into the host nuclear genome and subsequent gene loss from the endosymbiont. This process required that horizontally transferred genes become active and properly regulated despite inherent differences in genetic features between donor (endosymbiont) and recipient (host). Although this genetic reorganization is considered critical for early stages of organellogenesis, we have little knowledge about the mechanisms governing this process. The photosynthetic amoeba Paulinella micropora offers a unique opportunity to study early evolutionary events associated with organellogenesis and primary endosymbiosis. This amoeba harbors a “chromatophore,” a nascent photosynthetic organelle derived from a relatively recent cyanobacterial association (∼120 million years ago) that is independent of the evolution of primary plastids in plants (initiated ∼1.5 billion years ago). Analysis of the genome and transcriptome of Paulinella revealed that retrotransposition of endosymbiont-derived nuclear genes was critical for their domestication in the host. These retrocopied genes involved in photoprotection in cyanobacteria became expanded gene families and were “rewired,” acquiring light-responsive regulatory elements that function in the host. The establishment of host control of endosymbiont-derived genes likely enabled the cell to withstand photo-oxidative stress generated by oxygenic photosynthesis in the nascent organelle. These results provide insights into the genetic mechanisms and evolutionary pressures that facilitated the metabolic integration of the host–endosymbiont association and sustained the evolution of a photosynthetic organelle.

The accelerated molecular dynamics (aMD) method has recently been shown to enhance the sampling of biomolecules in molecular dynamics (MD) simulations, often by several orders of magnitude. Here, we describe an implementation of the aMD method for the OpenMM application layer that takes full advantage of graphics processing units (GPUs) computing. The aMD method is shown to work in combination with the AMOEBA polarizable force field (AMOEBA-aMD), allowing the simulation of long time-scale events with a polarizable force field. Benchmarks are provided to show that the AMOEBA-aMD method is efficiently implemented and produces accurate results in its standard parametrization. For the BPTI protein, we demonstrate that the protein structure described with AMOEBA remains stable even on the extended time scales accessed at high levels of accelerations. For the DNA repair metalloenzyme endonuclease IV, we show that the use of the AMOEBA force field is a significant improvement over fixed charged models for describing the enzyme active-site. The new AMOEBA-aMD method is publicly available (http://wiki.simtk.org/openmm/VirtualRepository) and promises to be interesting for studying complex systems that can benefit from both the use of a polarizable force field and enhanced sampling.

Amoebae are found all around the world and play an essential role in the carbon cycle in the environment. Therefore, the behavior of amoebae is a crucial factor when considering the global environment. Amoebae change their distribution through amoeboid locomotion, which are classified into several modes. In the pressure-driven mode, intracellular hydrostatic pressure generated by the contraction of cellular cortex actomyosin causes the pseudopod to extend. During amoeboid locomotion, the cellular surface exhibits dynamic deformation. Therefore, to understand the mechanism of amoeboid locomotion, it is important to characterize cellular membrane dynamics. Here, to clarify membrane dynamics during pressure-driven amoeboid locomotion, we developed a polkadot membrane staining method and performed light-sheet microscopy in Amoeba proteus, which exhibits typical pressure-driven amoeboid locomotion. It was observed that the whole cell membrane moved in the direction of movement, and the dorsal cell membrane in the posterior part of the cell moved more slowly than the other membrane. In addition, membrane complexity varied depending on the focused characteristic size of the membrane structure, and in general, the dorsal side was more complex than the ventral side. In summary, the membrane dynamics of Amoeba proteus during pressure-driven locomotion are asymmetric between the dorsal and ventral sides. This article has an associated interview with the co-first authors of the paper.

Multinuclearity is a widespread phenomenon across the living world, yet how it is achieved, and the potential related advantages, are not systematically understood. In this study, we investigate multinuclearity in amoebae. We observe that non-adherent amoebae are giant multinucleate cells compared to adherent ones. The cells solve their multinuclearity by a stretchy cytokinesis process with cytosolic bridge formation when adherence resumes. After initial adhesion to a new substrate, the progeny of the multinucleate cells is more numerous than the sibling cells generated from uninucleate amoebae. Hence, multinucleate amoebae show an advantage for population growth when the number of cells is quantified over time. Multiple nuclei per cell are observed in different amoeba species, and the lack of adhesion induces multinuclearity in diverse protists such as Acanthamoeba castellanii, Vermamoeba vermiformis, Naegleria gruberi and Hartmannella rhysodes. In this study, we observe that agitation induces a cytokinesis delay, which promotes multinuclearity. Hence, we propose the hypothesis that multinuclearity represents a physiological adaptation under non-adherent conditions that can lead to biologically relevant advantages.

Symbiotic interactions exist within a parasitism to mutualism continuum that is influenced, among others, by genes and context. Dynamics of intracellular invasion, replication, and prevalence may underscore both host survivability and symbiont stability. More infectious symbionts might exert higher corresponding costs to hosts, which could ultimately disadvantage both partners. Here, we quantify infection patterns of diverse Paraburkholderia symbiont genotypes in their amoeba host Dictyostelium discoideum and probe the relationship between these patterns and host outcomes. We exposed D. discoideum to thirteen strains of Paraburkholderia each belonging to one of the three symbiont species found to naturally infect D. discoideum: Paraburkholderia agricolaris, Paraburkholderia hayleyella, and Paraburkholderia bonniea. We quantified the infection prevalence and intracellular density of fluorescently labeled symbionts along with the final host population size using flow cytometry and confocal microscopy. We find that infection phenotypes vary across symbiont strains. Symbionts belonging to the same species generally display similar infection patterns but are interestingly distinct when it comes to host outcomes. This results in final infection loads that do not strongly correlate to final host outcomes, suggesting other genetic factors that are not a direct cause or consequence of symbiont abundance impact host fitness.

Predatory protozoa play an essential role in shaping microbial populations. Among these protozoa, Acanthamoeba are ubiquitous in the soil and aqueous environments inhabited by Listeria monocytogenes. Observations of predator-prey interactions between these two microorganisms revealed a predation strategy in which Acanthamoeba castellanii assemble L. monocytogenes in aggregates, termed backpacks, on their posterior. The rapid formation and specific location of backpacks led to the assumption that A. castellanii may recruit L. monocytogenes by releasing an attractant. However, this hypothesis has not been validated, and the mechanisms driving this process remained unknown. Here, we combined video microscopy, microfluidics, single-cell image analyses, and theoretical modeling to characterize predator-prey interactions of A. castellanii and L. monocytogenes and determined whether bacterial chemotaxis contributes to the backpack formation. Our results indicate that L. monocytogenes captures are not driven by chemotaxis. Instead, random encounters of bacteria with amoebae initialize bacterial capture and aggregation. This is supported by the strong correlation between experimentally derived capture rates and theoretical encounter models at the single-cell level. Observations of the spatial rearrangement of L. monocytogenes trapped by A. castellanii revealed that bacterial aggregation into backpacks is mainly driven by amoeboid locomotion. Overall, we show that two nonspecific, independent mechanisms, namely random encounters enhanced by bacterial motility and predator surface-bound locomotion, drive backpack formation, resulting in a bacterial aggregate on the amoeba ready for phagocytosis. Due to the prevalence of these two processes in the environment, we expect this strategy to be widespread among amoebae, contributing to their effectiveness as predators.

Willaertia magna C2c maky is a thermophilic free-living amoeba strain that showed ability to eliminate Legionella pneumophila, a pathogenic bacterium living in the aquatic environment. The amoeba industry has proposed the use of Willaertia magna as a natural biocide to control L. pneumophila proliferation in cooling towers. Here, transcriptomic and proteomic studies were carried out in order to expand knowledge on W. magna produced in a bioreactor. Illumina RNA-seq generated 217 million raw reads. A total of 8,790 transcripts were identified, of which 6,179 and 5,341 were assigned a function through comparisons with National Center of Biotechnology Information (NCBI) reference sequence and the Clusters of Orthologous Groups of proteins (COG) databases, respectively. To corroborate these transcriptomic data, we analyzed the W. magna proteome using LC-MS/MS. A total of 3,561 proteins were identified. The results of transcriptome and proteome analyses were highly congruent. Metabolism study showed that W. magna preferentially consumed carbohydrates and fatty acids to grow. Finally, an in-depth analysis has shown that W. magna produces several enzymes that are probably involved in the metabolism of secondary metabolites. Overall, our multi-omic study of W. magna opens the way to a better understanding of the genetics and biology of this amoeba.

Dictyostelium cells are forest soil amoebae, which feed on bacteria and proliferate as solitary cells until bacteria are consumed. Starvation triggers a change in life style, forcing cells to gather into aggregates to form multicellular organisms capable of cell differentiation and morphogenesis. As a soil amoeba and a phagocyte that grazes on bacteria as the obligate source of food, Dictyostelium could be a natural host of pathogenic bacteria. Indeed, many pathogens that occasionally infect humans are hosted for most of their time in protozoa or free-living amoebae, where evolution of their virulence traits occurs. Due to these features and its amenability to genetic manipulation, Dictyostelium has become a valuable model organism for studying strategies of both the host to resist infection and the pathogen to escape the defense mechanisms. Similarly to higher eukaryotes, iron homeostasis is crucial for Dictyostelium resistance to invasive bacteria. Iron is essential for Dictyostelium, as both iron deficiency or overload inhibit cell growth. The Dictyostelium genome shares with mammals many genes regulating iron homeostasis. Iron transporters of the Nramp (Slc11A) family are represented with two genes, encoding Nramp1 and Nramp2. Like the mammalian ortholog, Nramp1 is recruited to phagosomes and macropinosomes, whereas Nramp2 is a membrane protein of the contractile vacuole network, which regulates osmolarity. Nramp1 and Nramp2 localization in distinct compartments suggests that both proteins synergistically regulate iron homeostasis. Rather than by absorption via membrane transporters, iron is likely gained by degradation of ingested bacteria and efflux via Nramp1 from phagosomes to the cytosol. Nramp gene disruption increases Dictyostelium sensitivity to infection, enhancing intracellular growth of Legionella or Mycobacteria. Generation of mutants in other "iron genes" will help identify genes essential for iron homeostasis and resistance to pathogens.

Altruism can be favored by high relatedness among interactants. We tested the effect of relatedness in experimental populations of the social amoeba Dictyostelium discoideum, where altruism occurs in a starvation-induced social stage when some amoebae die to form a stalk that lifts the fertile spores above the soil facilitating dispersal. The single cells that aggregate during the social stage can be genetically diverse, which can lead to conflict over spore and stalk allocation. We mixed eight genetically distinct wild isolates and maintained twelve replicated populations at a high and a low relatedness treatment. After one and ten social generations we assessed the strain composition of the populations. We expected that some strains would be out-competed in both treatments. In addition, we expected that low relatedness might allow the persistence of social cheaters as it provides opportunity to exploit other strains.

The effects of electronic polarization in biomolecular interactions will differ depending on the local dielectric constant of the environment, such as in solvent, DNA, proteins, and membranes. Here the performance of the AMOEBA polarizable force field is evaluated under nonaqueous conditions by calculating the solvation free energies of small molecules in four common organic solvents. Results are compared with experimental data and equivalent simulations performed with the GAFF pairwise-additive force field. Although AMOEBA results give mean errors close to "chemical accuracy," GAFF performs surprisingly well, with statistically significantly more accurate results than AMOEBA in some solvents. However, for both models, free energies calculated in chloroform show worst agreement to experiment and individual solutes are consistently poor performers, suggesting non-potential-specific errors also contribute to inaccuracy. Scope for the improvement of both potentials remains limited by the lack of high quality experimental data across multiple solvents, particularly those of high dielectric constant. © 2016 The Authors. Journal of Computational Chemistry Published by Wiley Periodicals, Inc.

Naegleria fowleri is a pathogenic, thermophilic, free-living amoeba which causes primary amebic meningoencephalitis (PAM). Penetrating the olfactory mucosa, the brain-eating amoeba travels along the olfactory nerves, burrowing through the cribriform plate to its destination: the brain's frontal lobes. The amoeba thrives in warm, freshwater environments, with peak infection rates in the summer months and has a mortality rate of approximately 97%. A major contributor to the pathogen's high mortality is the lack of sensitivity of N. fowleri to current drug therapies, even in the face of combination-drug therapy. To enable rational drug discovery and design efforts we have pursued protein production and crystallography-based structure determination efforts for likely drug targets from N. fowleri. The genes were selected if they had homology to drug targets listed in Drug Bank or were nominated by primary investigators engaged in N. fowleri research. In 2017, 178 N. fowleri protein targets were queued to the Seattle Structural Genomics Center of Infectious Disease (SSGCID) pipeline, and to date 89 soluble recombinant proteins and 19 unique target structures have been produced. Many of the new protein structures are potential drug targets and contain structural differences compared to their human homologs, which could allow for the development of pathogen-specific inhibitors. Five of the structures were analyzed in more detail, and four of five show promise that selective inhibitors of the active site could be found. The 19 solved crystal structures build a foundation for future work in combating this devastating disease by encouraging further investigation to stimulate drug discovery for this neglected pathogen.

Transposable elements amplify in genomes as selfish DNA elements and challenge host fitness because their intrinsic integration steps during mobilization can compromise genome integrity. In gene-dense genomes, transposable elements are notably under selection to avoid insertional mutagenesis of host protein-coding genes. We describe an example of convergent evolution in the distantly related amoebozoan Dictyostelium discoideum and the yeast Saccharomyces cerevisiae, in which the D. discoideum retrotransposon DGLT-A and the yeast Ty3 element developed different mechanisms to facilitate position-specific integration at similar sites upstream of tRNA genes. Transcription of tRNA genes by RNA polymerase III requires the transcription factor complexes TFIIIB and TFIIIC. Whereas Ty3 recognizes tRNA genes mainly through interactions of its integrase with TFIIIB subunits, the DGLT-A-encoded ribonuclease H contacts TFIIIC subunit Tfc4 at an interface that covers tetratricopeptide repeats (TPRs) 7 and 8. A major function of this interface is to connect TFIIIC subcomplexes τA and τB and to facilitate TFIIIB assembly. During the initiation of tRNA gene transcription τB is displaced from τA, which transiently exposes the TPR 7/8 surface of Tfc4 on τA. We propose that the DGLT-A intasome uses this binding site to obtain access to genomic DNA for integration during tRNA gene transcription.

G-quadruplexes (G4) are non-canonical DNA and/or RNA secondary structures formed in guanine-rich regions. Given their over-representation in specific regions in the genome such as promoters and telomeres, they are likely to play important roles in key processes such as transcription, replication or RNA maturation. Putative G4-forming sequences (G4FS) have been reported in humans, yeast, bacteria, viruses and many organisms. Here we present the first mapping of G-quadruplex sequences in Dictyostelium discoideum, the social amoeba. 'Dicty' is an ameboid protozoan with a small (34 Mb) and extremely AT rich genome (78%). As a consequence, very few G4-prone motifs are expected. An in silico analysis of the Dictyostelium genome with the G4Hunter software detected 249-1055 G4-prone motifs, depending on G4Hunter chosen threshold. Interestingly, despite an even lower GC content (as compared to the whole Dicty genome), the density of G4 motifs in Dictyostelium promoters and introns is significantly higher than in the rest of the genome. Fourteen selected sequences located in important genes were characterized by a combination of biophysical and biochemical techniques. Our data show that these sequences form highly stable G4 structures under physiological conditions. Five Dictyostelium genes containing G4-prone motifs in their promoters were studied for the effect of a new G4-binding porphyrin derivative on their expression. Our results demonstrated that the new ligand significantly decreased their expression. Overall, our results constitute the first step to adopt Dictyostelium discoideum as a 'G4-poor' model for studies on G-quadruplexes.

Naegleria fowleri is a parasitic unicellular free living eukaryotic amoeba. The parasite spreads through contaminated water and causes primary amoebic meningoencephalitis (PAM). Therefore, it is of interest to understand its molecular pathogenesis. Hence, we analyzed the parasite genome for miRNAs (microRNAs) that are non-coding, single stranded RNA molecules. We identified 245 miRNAs using computational methods in N. fowleri, of which five miRNAs are conserved. The predicted miRNA targets were analyzed by using miRanda (software) and further studied the functions by subsequently annotating using AmiGo (a gene ontology web tool).

Cryptococcus neoformans is an environmental yeast that can cause opportunistic infections in humans. As infecting animals does not form part of its normal life-cycle, it has been proposed that the virulence traits that allow cryptococci to resist immune cells were selected through interactions with environmental phagocytes such as amoebae. Here, we investigate the interactions between C. neoformans and the social amoeba Dictyostelium discoideum. We show that like macrophages, D. discoideum is unable to kill C. neoformans upon phagocytosis. Despite this, we find that the yeast pass through the amoebae with an apparently normal phagocytic transit and are released alive by constitutive exocytosis after ~80 min. This is the canonical pathway in amoebae, used to dispose of indigestible material after nutrient extraction. Surprisingly however, we show that upon either genetic or pharmacological blockage of constitutive exocytosis, C. neoformans still escape from D. discoideum by a secondary mechanism. We demonstrate that constitutive exocytosis-independent egress is stochastic and actin-independent. This strongly resembles the non-lytic release of cryptococci by vomocytosis from macrophages, which do not perform constitutive exocytosis and normally retain phagocytosed material. Our data indicate that vomocytosis is functionally redundant for escape from amoebae, which thus may not be the primary driver for its evolutionary selection. Nonetheless, we show that vomocytosis of C. neoformans is mechanistically conserved in hosts ranging from amoebae to man, providing new avenues to understand this poorly-understood but important virulence mechanism.

Dictyostelium discoideum is a eukaryotic microbial model system for multicellular development, cell-cell signaling, and social behavior. Key models of social evolution require an understanding of genetic relationships between individuals across the genome or possibly at specific genes, but the nature of variation within D. discoideum is largely unknown. We re-sequenced 137 gene fragments in wild North American strains of D. discoideum and examined the levels and patterns of nucleotide variation in this social microbial species. We observe surprisingly low levels of nucleotide variation in D. discoideum across these strains, with a mean nucleotide diversity (pi) of 0.08%, and no strong population stratification among North American strains. We also do not find any clear relationship between nucleotide divergence between strains and levels of social dominance and kin discrimination. Kin discrimination experiments, however, show that strains collected from the same location show greater ability to distinguish self from non-self than do strains from different geographic areas. This suggests that a greater ability to recognize self versus non-self may arise among strains that are more likely to encounter each other in nature, which would lead to preferential formation of fruiting bodies with clonemates and may prevent the evolution of cheating behaviors within D. discoideum populations. Finally, despite the fact that sex has rarely been observed in this species, we document a rapid decay of linkage disequilibrium between SNPs, the presence of recombinant genotypes among natural strains, and high estimates of the population recombination parameter rho. The SNP data indicate that recombination is widespread within D. discoideum and that sex as a form of social interaction is likely to be an important aspect of the life cycle.