Repeated, selection-driven genome reduction of accessory genes in experimental populations

Ecological divergence of marine bacteria Alteromonas mediterranea

Alteromonas mediterranea, originally designated as A. macleodii, is a deep-sea ecotype that plays an important ecological role in the ocean. However, a comprehensive understanding of their biogeographic distribution and evolutionary histories remains limited. In this study, our analysis indicated that A. mediterranea members could adapt contrasting marine ecosystems and flourish in nutrient-rich habitats such as feces and coral reefs. No significant correlations between the relative abundance of A. mediterranea members and the environmental variables were identified. Phylogenetic analysis and geographic patterns of A. mediterranea strains suggested that they could be clustered into two clades (clade Ⅰ and clade Ⅱ). In contrast, many distinct genomic traits exist between these clades, such as the complete genes encoding cytochrome o ubiquinol oxidase only involved in clade Ⅱ. Genes were more likely to be lost in the evolutionary history of A. mediterranea relatives. Gene loss might be a major force in all phylogenetic groups driving the distinct clades. Adaptation to different biotopes resulted in the functional differentiation of A. mediterranea members, with the loss of genes encoding carbohydrate-active enzymes. Genes acquired horizontally from unclassified bacteria, and Proteobacteria represented by Gammaproteobacteria played key roles in the functional diversification of A. mediterranea in marine habitats. Given these data, these results are useful for information supplementation of A. mediterranea strains, particularly for making significant advances in understanding marine microbial ecology within different clonal frames using genome-wide recruitments.

Machine learning reveals the dynamic importance of accessory sequences for Salmonella outbreak clustering

Bacterial typing at whole-genome scales is now feasible owing to decreasing costs in high-throughput sequencing and the recent advances in computation. The unprecedented resolution of whole-genome typing is achieved by genotyping the variable segments of bacterial genomes that can fluctuate significantly in gene content. However, due to the transient and hypervariable nature of many accessory elements, the value of the added resolution in outbreak investigations remains disputed. To assess the analytical value of bacterial accessory genomes in clustering epidemiologically related cases, we trained classifiers on a set of genomes collected from 24 Salmonella enterica outbreaks of food, animal, or environmental origin. The models demonstrated high precision and recall on unseen test data with near-perfect accuracy in classifying clonal and short-term outbreaks. Annotating the genomic features important for cluster classification revealed functional enrichment of molecular fingerprints in genes involved in membrane transportation, trafficking, and carbohydrate metabolism. Importantly, we discovered polymorphisms in mobile genetic elements (MGEs) and gain/loss of MGEs to be informative in defining outbreak clusters. To quantify the ability of MGE variations to cluster outbreak clones, we devised a reference-free tree-building algorithm inspired by colored de Bruijn graphs, which enabled topological comparisons between MGE and standard typing methods. Systematic evaluation of clustering MGEs on an unseen dataset of 34 Salmonella outbreaks yielded mixed results that exemplified the power of accessory sequence variations when core genomes of unrelated cases are insufficiently discriminatory, as well as the distortion of outbreak signals by microevolution events or the incomplete assembly of MGEs.

IMPORTANCE

Gene-by-gene typing is widely used to detect clusters of foodborne illnesses that share a common origin. It remains actively debated whether the inclusion of accessory sequences in bacterial typing schema is informative or deleterious for cluster definitions in outbreak investigations due to the potential confounding effects of horizontal gene transfer. By training machine learning models on a curated set of historical Salmonella outbreaks, we revealed an enriched presence of outbreak distinguishing features in a wide range of mobile genetic elements. Systematic comparison of the efficacy of clustering different accessory elements against standard sequence typing methods led to our cataloging of scenarios where accessory sequence variations were beneficial and uninformative to resolving outbreak clusters. The presented work underscores the complexity of the molecular trends in enteric outbreaks and seeks to inspire novel computational ways to exploit whole-genome sequencing data in enteric disease surveillance and management.

Synthesis of a Truncated Microcystin Tetrapeptide Molecule from a Partial Mcy Gene Cluster in Microcystis Cultures and Blooms

Microcystis spp. threaten freshwater ecosystems through the proliferation of cyanobacterial harmful algal blooms (cyanoHABs) and production of the hepatotoxin, microcystin. While microcystin and its biosynthesis pathway, encoded by the mcy genes, have been well studied for over 50 years, a recent study found that Microcystis populations in western Lake Erie contain a transcriptionally active partial mcy operon, in which the A2 domain of mcyA and mcyB-C are present but the mcyD-J genes are absent. Here, we investigate the potential biosynthetic products and the evolutionary history of this partial operon. Our results reveal two candidate tetrapeptide constructs, with an X variable position, to be produced by strains with the partial operon. The partial operon appears necessary and sufficient for tetrapeptide biosynthesis and likely evolved from a single ancestor hundreds to tens of thousands of years ago. Bioactivity screens using Hep3B cells indicate a mild elevation of some markers of hepatotoxicity and inflammation, suggesting the need to further assess the effects of these novel secondary metabolites on freshwater ecosystems and public health. The need to assess these effects is even more pressing given the detection of tetrapeptides in both culture and western Lake Erie, which is a vital source of fresh water. Results from this study emphasize previous findings in which novel bacterial secondary metabolites may be derived from the molecular evolution of existing biosynthetic machinery under different environmental forcings.

Differential expression of the type III secretion system genes in Yersinia ruckeri: Preliminary investigations in different environmental conditions

Type III secretion system (T3SS) is an important virulence system in Gram-negative bacteria. In this investigation, different environmental conditions that regulate the expression of T3SS genes in Yersinia ruckeri were investigated aimed at obtaining a better understanding about its modulation after various environmental challenges. Four isolates of Y. ruckeri CSF007-82, ATCC29473, A7959-11 and YRNC10 were cultivated under the diverse in vitro challenges iron depletion, high salt, low pH and in the presence of fish serum or in the fish cell culture (Chinook Salmon Embryo - CHSE). The transcriptional modulation of the chromosomal genes ysaV, ysaC, ysaJ and prgH of ysa were investigated using quantitative real-time PCR. The expression of prgH, ysaV, ysaC and ysaJ was differentially expressed in all four strains under evaluation. The highest gene expression levels were observed for Y. ruckeri YRNC10 AN after addition of 0.3 M NaCl in Luria Bertani broth. The results obtained from this study provide initial insights into T3SS responses in Y. ruckeri, which pave the way for further studies aimed at expanding our knowledge on the functional roles of the T3SS genes in Y. ruckeri.

Domestication and evolutionary histories of specialized gut symbionts across cephalotine ants

The evolution of animals and their gut symbionts is a complex phenomenon, obscured by lability and diversity. In social organisms, transmission of symbionts among relatives may yield systems with more stable associations. Here, we study the history of a social insect symbiosis involving cephalotine ants and their extracellular gut bacteria, which come predominantly from host-specialized lineages. We perform multi-locus phylogenetics for symbionts from nine bacterial orders, and map prior amplicon sequence data to lineage-assigned symbiont genomes, studying distributions of rigorously defined symbionts across 20 host species. Based on monophyly and additional hypothesis testing, we estimate that these specialized gut bacteria belong to 18 distinct lineages, of which 15 have been successfully isolated and cultured. Several symbiont lineages showed evidence for domestication events that occurred later in cephalotine evolutionary history, and only one lineage was ubiquitously detected in all 20 host species and 48 colonies sampled with amplicon 16S rRNA sequencing. We found evidence for phylogenetically constrained distributions in four symbionts, suggesting historical or genetic impacts on community composition. Two lineages showed evidence for frequent intra-lineage co-infections, highlighting the potential for niche divergence after initial domestication. Nearly all symbionts showed evidence for occasional host switching, but four may, more often, co-diversify with their hosts. Through our further assessment of symbiont localization and genomic functional profiles, we demonstrate distinct niches for symbionts with shared evolutionary histories, prompting further questions on the forces underlying the evolution of hosts and their gut microbiomes.

Genome streamlining in Parcubacteria transitioning from soil to groundwater

BACKGROUND

To better understand the influence of habitat on the genetic content of bacteria, with a focus on members of Candidate Phyla Radiation (CPR) bacteria, we studied the effects of transitioning from soil via seepage waters to groundwater on genomic composition of ultra-small Parcubacteria, the dominating CPR class in seepage waters, using genome resolved metagenomics.

RESULTS

Bacterial metagenome-assembled genomes (MAGs), (318 total, 32 of Parcubacteria) were generated from seepage waters and compared directly to groundwater counterparts. The estimated average genome sizes of members of major phyla Proteobacteria, Bacteroidota and Cand. Patescibacteria (Candidate Phyla Radiation - CPR bacteria) were significantly higher in soil-seepage water as compared to their groundwater counterparts. Seepage water Parcubacteria (Paceibacteria) exhibited 1.18-fold greater mean genome size and 2-fold lower mean proportion of pseudogenes than those in groundwater. Bacteroidota and Proteobacteria also showed a similar trend of reduced genomes in groundwater compared to seepage. While exploring gene loss and adaptive gains in closely related CPR lineages in groundwater, we identified a membrane protein, and a lipoglycopeptide resistance gene unique to a seepage Parcubacterium genome. A nitrite reductase gene was also identified and was unique to the groundwater Parcubacteria genomes, likely acquired from other planktonic microbes via horizontal gene transfer.

CONCLUSIONS

Overall, our data suggest that bacteria in seepage waters, including ultra-small Parcubacteria, have significantly larger genomes and higher metabolic enrichment than their groundwater counterparts, highlighting possible genome streamlining of the latter in response to habitat selection in an oligotrophic environment.

Reduced social function in experimentally evolved Dictyostelium discoideum implies selection for social conflict in nature

Many microbes interact with one another, but the difficulty of directly observing these interactions in nature makes interpreting their adaptive value complicated. The social amoeba Dictyostelium discoideum forms aggregates wherein some cells are sacrificed for the benefit of others. Within chimaeric aggregates containing multiple unrelated lineages, cheaters can gain an advantage by undercontributing, but the extent to which wild D. discoideum has adapted to cheat is not fully clear. In this study, we experimentally evolved D. discoideum in an environment where there were no selective pressures to cheat or resist cheating in chimaeras. Dictyostelium discoideum lines grown in this environment evolved reduced competitiveness within chimaeric aggregates and reduced ability to migrate during the slug stage. By contrast, we did not observe a reduction in cell number, a trait for which selection was not relaxed. The observed loss of traits that our laboratory conditions had made irrelevant suggests that these traits were adaptations driven and maintained by selective pressures D. discoideum faces in its natural environment. Our results suggest that D. discoideum faces social conflict in nature, and illustrate a general approach that could be applied to searching for social or non-social adaptations in other microbes.

The Multifaceted Role of Homologous Recombination in a Fastidious Bacterial Plant Pathogen

Homologous recombination plays a key function in the evolution of bacterial genomes. Within Xylella fastidiosa, an emerging plant pathogen with increasing host and geographic ranges, it has been suggested that homologous recombination facilitates host switching, speciation, and the development of virulence. We used 340 whole-genome sequences to study the relationship between inter- and intrasubspecific homologous recombination, random mutation, and natural selection across individual X. fastidiosa genes. Individual gene orthologs were identified and aligned, and a maximum likelihood (ML) gene tree was generated. Each gene alignment and tree pair were then used to calculate gene-wide and branch-specific r/m values (relative effect of recombination to mutation), gene-wide and branch-site nonsynonymous over synonymous substitution rates (dN/dS values; episodic selection), and branch length (as a proxy for mutation rate). The relationships between these variables were evaluated at the global level (i.e., for all genes among and within a subspecies), among specific functional classes (i.e., COGs), and between pangenome components (i.e., accessory versus core genes). Our analysis showed that r/m varied widely among genes as well as across X. fastidiosa subspecies. While r/m and dN/dS values were positively correlated in some instances (e.g., core genes in X. fastidiosa subsp. fastidiosa and both core and accessory genes in X. fastidiosa subsp. multiplex), low correlation coefficients suggested no clear biological significance. Overall, our results indicate that, in addition to its adaptive role in certain genes, homologous recombination acts as a homogenizing and a neutral force across phylogenetic clades, gene functional groups, and pangenome components. IMPORTANCE There is ample evidence that homologous recombination occurs frequently in the economically important plant pathogen Xylella fastidiosa. Homologous recombination has been known to occur among sympatric subspecies and is associated with host-switching events and virulence-linked genes. As a consequence, is it generally assumed that recombinant events in X. fastidiosa are adaptive. This mindset influences expectations of how homologous recombination acts as an evolutionary force as well as how management strategies for X. fastidiosa diseases are determined. Yet, homologous recombination plays roles beyond that of a source for diversification and adaptation. Homologous recombination can act as a DNA repair mechanism, as a means to facilitate nucleotide compositional change, as a homogenization mechanism within populations, or even as a neutral force. Here, we provide a first assessment of long-held beliefs regarding the general role of recombination in adaptation for X. fastidiosa. We evaluate gene-specific variations in homologous recombination rate across three X. fastidiosa subspecies and its relationship to other evolutionary forces (e.g., natural selection, mutation, etc.). These data were used to assess the role of homologous recombination in X. fastidiosa evolution.

On the evolution of natural product biosynthesis

Natural products are the raw material for drug discovery programmes. Bioactive natural products are used extensively in medicine and agriculture and have found utility as antibiotics, immunosuppressives, anti-cancer drugs and anthelminthics. Remarkably, the natural role and what mechanisms drive evolution of these molecules is relatively poorly understood. The exponential increase in genome and chemical data in recent years, coupled with technical advances in bioinformatics and genetics have enabled progress to be made in understanding the evolution of biosynthetic gene clusters and the products of their enzymatic machinery. Here we discuss the diversity of natural products, incorporating the mechanisms that govern evolution of metabolic pathways and how this can be applied to biosynthetic gene clusters. We build on the nomenclature of natural products in terms of primary, integrated, secondary and specialised metabolism and place this within an ecology-evolutionary-developmental biology framework. This eco-evo-devo framework we believe will help to clarify the nature and use of the term specialised metabolites in the future.

Mitochondrial genomes in the iconic reindeer lichens: Architecture, variation, and synteny across multiple evolutionary scales

Variation in mitochondrial genome composition across intraspecific, interspecific, and higher taxonomic scales has been little studied in lichen obligate symbioses. Cladonia is one of the most diverse and ecologically important lichen genera, with over 500 species representing an array of unique morphologies and chemical profiles. Here, we assess mitochondrial genome diversity and variation in this flagship genus, with focused sampling of two clades of the "true" reindeer lichens, Cladonia subgenus Cladina, and additional genomes from nine outgroup taxa. We describe composition and architecture at the gene and the genome scale, examining patterns in organellar genome size in larger taxonomic groups in Ascomycota. Mitochondrial genomes of Cladonia, Pilophorus, and Stereocaulon were consistently larger than those of Lepraria and contained more introns, suggesting a selective pressure in asexual morphology in Lepraria driving it toward genomic simplification. Collectively, lichen mitochondrial genomes were larger than most other fungal life strategies, reaffirming the notion that coevolutionary streamlining does not correlate to genome size reductions. Genomes from Cladonia ravenelii and Stereocaulon pileatum exhibited ATP9 duplication, bearing paralogs that may still be functional. Homing endonuclease genes (HEGs), though scarce in Lepraria, were diverse and abundant in Cladonia, exhibiting variable evolutionary histories that were sometimes independent of the mitochondrial evolutionary history. Intraspecific HEG diversity was also high, with C. rangiferina especially bearing a range of HEGs with one unique to the species. This study reveals a rich history of events that have transformed mitochondrial genomes of Cladonia and related genera, allowing future study alongside a wealth of assembled genomes.

Seasonal variation of microbial community and methane metabolism in coalbed water in the Erlian Basin, China

Coalbed water is a semi-open system connecting underground coalbeds with the external environment. Microorganisms in coalbed water play an important role in coal biogasification and the carbon cycle. The community assemblages of microorganisms in such a dynamic system are not well understood. Here, we used high-throughput sequencing and metagenomic analysis to investigate microbial community structure and identify the potential functional microorganisms involved in methane metabolism in coalbed water in the Erlian Basin, a preferred low-rank coal bed methane (CBM) exploration and research area in China. The results showed that there were differences in the responses of bacteria and archaea to seasonal variation. Bacterial community structure was affected by seasonal variation but archaea was not. Methane oxidation metabolism dominated by Methylomonas and methanogenesis metabolism dominated by Methanobacterium may exist simultaneously in coalbed water.

Bacterial plasmid-associated and chromosomal proteins have fundamentally different properties in protein interaction networks

Plasmids facilitate horizontal gene transfer, which enables the diversification of pathogens into new anatomical and environmental niches, implying that plasmid-encoded genes can cooperate well with chromosomal genes. We hypothesise that such mobile genes are functionally different to chromosomal ones due to this ability to encode proteins performing non-essential functions like antimicrobial resistance and traverse distinct host cells. The effect of plasmid-driven gene gain on protein-protein interaction network topology is an important question in this area. Moreover, the extent to which these chromosomally- and plasmid-encoded proteins interact with proteins from their own groups compared to the levels with the other group remains unclear. Here, we examined the incidence and protein-protein interactions of all known plasmid-encoded proteins across representative specimens from most bacteria using all available plasmids. We found that plasmid-encoded genes constitute ~ 0.65% of the total number of genes per bacterial sample, and that plasmid genes are preferentially associated with different species but had limited taxonomical power beyond this. Surprisingly, plasmid-encoded proteins had both more protein-protein interactions compared to chromosomal proteins, countering the hypothesis that genes with higher mobility rates should have fewer protein-level interactions. Nonetheless, topological analysis and investigation of the protein-protein interaction networks' connectivity and change in the number of independent components demonstrated that the plasmid-encoded proteins had limited overall impact in > 96% of samples. This paper assembled extensive data on plasmid-encoded proteins, their interactions and associations with diverse bacterial specimens that is available for the community to investigate in more detail.

Bacterial genome reductions: Tools, applications, and challenges

Bacterial cells are widely used to produce value-added products due to their versatility, ease of manipulation, and the abundance of genome engineering tools. However, the efficiency of producing these desired biomolecules is often hindered by the cells’ own metabolism, genetic instability, and the toxicity of the product. To overcome these challenges, genome reductions have been performed, making strains with the potential of serving as chassis for downstream applications. Here we review the current technologies that enable the design and construction of such reduced-genome bacteria as well as the challenges that limit their assembly and applicability. While genomic reductions have shown improvement of many cellular characteristics, a major challenge still exists in constructing these cells efficiently and rapidly. Computational tools have been created in attempts at minimizing the time needed to design these organisms, but gaps still exist in modelling these reductions in silico. Genomic reductions are a promising avenue for improving the production of value-added products, constructing chassis cells, and for uncovering cellular function but are currently limited by their time-consuming construction methods. With improvements to and the creation of novel genome editing tools and in silico models, these approaches could be combined to expedite this process and create more streamlined and efficient cell factories.

Zooplankton as a Transitional Host for Escherichia coli in Freshwater

This study shows that Escherichia coli can be temporarily enriched in zooplankton under natural conditions and that these bacteria can belong to different phylogroups and sequence types (STs), including environmental, clinical, and animal isolates. We isolated 10 E. coli strains and sequenced the genomes of two of them. Phylogenetically, the two isolates were closer to strains isolated from poultry meat than to freshwater E. coli, albeit their genomes were smaller than those of the poultry isolates. After isolation and fluorescent protein tagging of strains ED1 and ED157, we show that Daphnia sp. can take up these strains and release them alive again, thus becoming a temporary host for E. coli. In a chemostat experiment, we show that this association does not prolong bacterial long-term survival, but at low abundances it also does not significantly reduce bacterial numbers. We demonstrate that E. coli does not belong to the core microbiota of Daphnia, suffers from competition by the natural Daphnia microbiota, but can profit from its carapax to survive in water. All in all, this study suggests that the association of E. coli with Daphnia is only temporary, but the cells are viable therein, and this might allow encounters with other bacteria for genetic exchange and potential genomic adaptation to the freshwater environment.

IMPORTANCE The contamination of freshwater with feces-derived bacteria is a major concern regarding drinking water acquisition and recreational activities. Ecological interactions promoting their persistence are still very scarcely studied. This study, which analyses the survival of E. coli in the presence of zooplankton, is thus of ecological and water safety relevance.

Aerobic Methoxydotrophy: Growth on Methoxylated Aromatic Compounds by Methylobacteriaceae

Pink-pigmented facultative methylotrophs have long been studied for their ability to grow on reduced single-carbon (C1) compounds. The C1 groups that support methylotrophic growth may come from a variety of sources. Here, we describe a group of Methylobacterium strains that can engage in methoxydotrophy: they can metabolize the methoxy groups from several aromatic compounds that are commonly the product of lignin depolymerization. Furthermore, these organisms can utilize the full aromatic ring as a growth substrate, a phenotype that has rarely been described in Methylobacterium. We demonstrated growth on p-hydroxybenzoate, protocatechuate, vanillate, and ferulate in laboratory culture conditions. We also used comparative genomics to explore the evolutionary history of this trait, finding that the capacity for aromatic catabolism is likely ancestral to two clades of Methylobacterium, but has also been acquired horizontally by closely related organisms. In addition, we surveyed the published metagenome data to find that the most abundant group of aromatic-degrading Methylobacterium in the environment is likely the group related to Methylobacterium nodulans, and they are especially common in soil and root environments. The demethoxylation of lignin-derived aromatic monomers in aerobic environments releases formaldehyde, a metabolite that is a potent cellular toxin but that is also a growth substrate for methylotrophs. We found that, whereas some known lignin-degrading organisms excrete formaldehyde as a byproduct during growth on vanillate, Methylobacterium do not. This observation is especially relevant to our understanding of the ecology and the bioengineering of lignin degradation.

Environmental stress leads to genome streamlining in a widely distributed species of soil bacteria

Bacteria have highly flexible pangenomes, which are thought to facilitate evolutionary responses to environmental change, but the impacts of environmental stress on pangenome evolution remain unclear. Using a landscape pangenomics approach, I demonstrate that environmental stress leads to consistent, continuous reduction in genome content along four environmental stress gradients (acidity, aridity, heat, salinity) in naturally occurring populations of Bradyrhizobium diazoefficiens (widespread soil-dwelling plant mutualists). Using gene-level network and duplication functional traits to predict accessory gene distributions across environments, genes predicted to be superfluous are more likely lost in high stress, while genes with multi-functional roles are more likely retained. Genes with higher probabilities of being lost with stress contain significantly higher proportions of codons under strong purifying and positive selection. Gene loss is widespread across the entire genome, with high gene-retention hotspots in close spatial proximity to core genes, suggesting Bradyrhizobium has evolved to cluster essential-function genes (accessory genes with multifunctional roles and core genes) in discrete genomic regions, which may stabilise viability during genomic decay. In conclusion, pangenome evolution through genome streamlining are important evolutionary responses to environmental change. This raises questions about impacts of genome streamlining on the adaptive capacity of bacterial populations facing rapid environmental change.

Squash root microbiome transplants and metagenomic inspection for in situ arid adaptations

Arid zones contain a diverse set of microbes capable of survival under dry conditions, some of which can form relationships with plants under drought stress conditions to improve plant health. We studied squash (Cucurbita pepo L.) root microbiome under historically arid and humid sites, both in situ and performing a common garden experiment. Plants were grown in soils from sites with different drought levels, using in situ collected soils as the microbial source. We described and analyzed bacterial diversity by 16S rRNA gene sequencing (N = 48) from the soil, rhizosphere, and endosphere. Proteobacteria were the most abundant phylum present in humid and arid samples, while Actinobacteriota abundance was higher in arid ones. The β-diversity analyses showed split microbiomes between arid and humid microbiomes, and aridity and soil pH levels could explain it. These differences between humid and arid microbiomes were maintained in the common garden experiment, showing that it is possible to transplant in situ diversity to the greenhouse. We detected a total of 1009 bacterial genera; 199 exclusively associated with roots under arid conditions. By 16S and shotgun metagenomics, we identified dry-associated taxa such as Cellvibrio, Ensifer adhaerens, and Streptomyces flavovariabilis. With shotgun metagenomic sequencing of rhizospheres (N = 6), we identified 2969 protein families in the squash core metagenome and found an increased number of exclusively protein families from arid (924) than humid samples (158). We found arid conditions enriched genes involved in protein degradation and folding, oxidative stress, compatible solute synthesis, and ion pumps associated with osmotic regulation. Plant phenotyping allowed us to correlate bacterial communities with plant growth. Our study revealed that it is possible to evaluate microbiome diversity ex-situ and identify critical species and genes involved in plant-microbe interactions in historically arid locations.

Experimental evolution of the megaplasmid pMPPla107 in Pseudomonas stutzeri enables identification of genes contributing to sensitivity to an inhibitory agent

Horizontally transferred elements, such as plasmids, can burden host cells with various metabolic and fitness costs and may lead to other potentially detrimental phenotypic effects. Acquisition of the Pseudomonas syringae megaplasmid pMPPla107 by various Pseudomonads causes sensitivity to a growth-inhibiting substance that is produced in cultures by Pseudomonads during growth under standard laboratory conditions. After approximately 500 generations of laboratory passage of Pseudomonas stutzeri populations containing pMPPla107, strains from two out of six independent passage lines displayed resistance to this inhibitory agent. Resistance was transferable and is, therefore, associated with mutations occurring on pMPPla107. Resequencing experiments demonstrated that resistance is likely due to a large deletion on the megaplasmid in one line, and to a nonsynonymous change in an uncharacterized megaplasmid locus in the other strain. We further used allele exchange experiments to confirm that resistance is due to this single amino acid change in a previously uncharacterized megaplasmid protein, which we name SkaA. These results provide further evidence that costs and phenotypic changes associated with horizontal gene transfer can be compensated through single mutational events and emphasize the power of experimental evolution and resequencing to better understand the genetic basis of evolved phenotypes. This article is part of the theme issue 'The secret lives of microbial mobile genetic elements'.

What makes a megaplasmid?

Naturally occurring plasmids come in different sizes. The smallest are less than a kilobase of DNA, while the largest can be over three orders of magnitude larger. Historically, research has tended to focus on smaller plasmids that are usually easier to isolate, manipulate and sequence, but with improved genome assemblies made possible by long-read sequencing, there is increased appreciation that very large plasmids—known as megaplasmids—are widespread, diverse, complex, and often encode key traits in the biology of their host microorganisms. Why are megaplasmids so big? What other features come with large plasmid size that could affect bacterial ecology and evolution? Are megaplasmids 'just' big plasmids, or do they have distinct characteristics? In this perspective, we reflect on the distribution, diversity, biology, and gene content of megaplasmids, providing an overview to these large, yet often overlooked, mobile genetic elements.

This article is part of the theme issue ‘The secret lives of microbial mobile genetic elements’.

Ecological and evolutionary solutions to the plasmid paradox

The 'plasmid paradox' arises because, although plasmids are common features of bacterial genomes, theoretically they should not exist: rates of conjugation were believed insufficient to allow plasmids to persist by infectious transmission, whereas the costs of plasmid maintenance meant that plasmids should be purged by negative selection regardless of whether they encoded beneficial accessory traits because these traits should eventually be captured by the chromosome, enabling the loss of the redundant plasmid. In the decade since the plasmid paradox was described, new data and theory show that a range of ecological and evolutionary mechanisms operate in bacterial populations and communities to explain the widespread distribution and stable maintenance of plasmids. We conclude, therefore, that multiple solutions to the plasmid paradox are now well understood. The current challenge for the field, however, is to better understand how these solutions operate in natural bacterial communities to explain and predict the distribution of plasmids and the dynamics of the horizontal gene transfer that they mediate in bacterial (pan)genomes.

Microbial generalist or specialist: Intraspecific variation and dormancy potential matter

Microbial generalists and specialists coexist in the soil environment while having distinctive impacts on microbial community dynamics. In microbial ecology, the underlying mechanisms as to why a species is a generalist or a specialist remain ambiguous. Herein, we collected soils across a national scale and identified bacterial generalists and specialists according to niche breadth at the species level (OTU level), and the single-nucleotide differences in each species were measured to investigate intraspecific variation (at zero-radius OTU level). Compared with that of the specialists, the intraspecific variation of the generalists was much higher, which ensured their wider niche breadth and lower variability. The higher asynchrony and different niche preferences of conspecific individuals and the higher dormancy potential within the generalists further contributed to their stability in varying environments. Besides, generalists were less controlled by environmental filtering, which was indicated by the stronger signature of stochastic processes in their assembly, and had higher diversification and transition rates that allowed them to adapt to environmental changes to a greater extent than specialists. Overall, this study provides a new comprehensive understanding of the rules of assembly and the evolutionary roles of bacterial generalists and specialists. It also highlights the importance of intraspecific variation and the dormancy potential in the stability of species.

Plasmid fitness costs are caused by specific genetic conflicts enabling resolution by compensatory mutation

Plasmids play an important role in bacterial genome evolution by transferring genes between lineages. Fitness costs associated with plasmid carriage are expected to be a barrier to gene exchange, but the causes of plasmid fitness costs are poorly understood. Single compensatory mutations are often sufficient to completely ameliorate plasmid fitness costs, suggesting that such costs are caused by specific genetic conflicts rather than generic properties of plasmids, such as their size, metabolic burden, or gene expression level. By combining the results of experimental evolution with genetics and transcriptomics, we show here that fitness costs of 2 divergent large plasmids in Pseudomonas fluorescens are caused by inducing maladaptive expression of a chromosomal tailocin toxin operon. Mutations in single genes unrelated to the toxin operon, and located on either the chromosome or the plasmid, ameliorated the disruption associated with plasmid carriage. We identify one of these compensatory loci, the chromosomal gene PFLU4242, as the key mediator of the fitness costs of both plasmids, with the other compensatory loci either reducing expression of this gene or mitigating its deleterious effects by up-regulating a putative plasmid-borne ParAB operon. The chromosomal mobile genetic element Tn6291, which uses plasmids for transmission, remained up-regulated even in compensated strains, suggesting that mobile genetic elements communicate through pathways independent of general physiological disruption. Plasmid fitness costs caused by specific genetic conflicts are unlikely to act as a long-term barrier to horizontal gene transfer (HGT) due to their propensity for amelioration by single compensatory mutations, helping to explain why plasmids are so common in bacterial genomes.

The cost of evolved constitutive lac gene expression is usually, but not always, maintained during evolution of generalist populations

Beneficial mutations can become costly following an environmental change. Compensatory mutations can relieve these costs, while not affecting the selected function, so that the benefits are retained if the environment shifts back to be similar to the one in which the beneficial mutation was originally selected. Compensatory mutations have been extensively studied in the context of antibiotic resistance, responses to specific genetic perturbations, and in the determination of interacting gene network components. Few studies have focused on the role of compensatory mutations during more general adaptation, especially as the result of selection in fluctuating environments where adaptations to different environment components may often involve trade-offs. We examine whether costs of a mutation in lacI, which deregulated the expression of the lac operon in evolving populations of Escherichia coli bacteria, were compensated. This mutation occurred in multiple replicate populations selected in environments that fluctuated between growth on lactose, where the mutation was beneficial, and on glucose, where it was deleterious. We found that compensation for the cost of the lacI mutation was rare, but, when it did occur, it did not negatively affect the selected benefit. Compensation was not more likely to occur in a particular evolution environment. Compensation has the potential to remove pleiotropic costs of adaptation, but its rarity indicates that the circumstances to bring about the phenomenon may be peculiar to each individual or impeded by other selected mutations.

Forecasting of phenotypic and genetic outcomes of experimental evolution in Pseudomonas protegens

Experimental evolution with microbes is often highly repeatable under identical conditions, suggesting the possibility to predict short-term evolution. However, it is not clear to what degree evolutionary forecasts can be extended to related species in non-identical environments, which would allow testing of general predictive models and fundamental biological assumptions. To develop an extended model system for evolutionary forecasting, we used previous data and models of the genotype-to-phenotype map from the wrinkly spreader system in Pseudomonas fluorescens SBW25 to make predictions of evolutionary outcomes on different biological levels for Pseudomonas protegens Pf-5. In addition to sequence divergence (78% amino acid and 81% nucleotide identity) for the genes targeted by mutations, these species also differ in the inability of Pf-5 to make cellulose, which is the main structural basis for the adaptive phenotype in SBW25. The experimental conditions were changed compared to the SBW25 system to test if forecasts were extendable to a non-identical environment. Forty-three mutants with increased ability to colonize the air-liquid interface were isolated, and the majority had reduced motility and was partly dependent on the Pel exopolysaccharide as a structural component. Most (38/43) mutations are expected to disrupt negative regulation of the same three diguanylate cyclases as in SBW25, with a smaller number of mutations in promoter regions, including an uncharacterized polysaccharide synthase operon. A mathematical model developed for SBW25 predicted the order of the three main pathways and the genes targeted by mutations, but differences in fitness between mutants and mutational biases also appear to influence outcomes. Mutated regions in proteins could be predicted in most cases (16/22), but parallelism at the nucleotide level was low and mutational hot spot sites were not conserved. This study demonstrates the potential of short-term evolutionary forecasting in experimental populations and provides testable predictions for evolutionary outcomes in other Pseudomonas species.

A Metzincin and TIMP-Like Protein Pair of a Phage Origin Sensitize Listeria monocytogenes to Phage Lysins and Other Cell Wall Targeting Agents

Infection of mammalian cells by Listeria monocytogenes (Lm) was shown to be facilitated by its phage elements. In a search for additional phage remnants that play a role in Lm's lifecycle, we identified a conserved locus containing two XRE regulators and a pair of genes encoding a secreted metzincin protease and a lipoprotein structurally similar to a TIMP-family metzincin inhibitor. We found that the XRE regulators act as a classic CI/Cro regulatory switch that regulates the expression of the metzincin and TIMP-like genes under intracellular growth conditions. We established that when these genes are expressed, their products alter Lm morphology and increase its sensitivity to phage mediated lysis, thereby enhancing virion release. Expression of these proteins also sensitized the bacteria to cell wall targeting compounds, implying that they modulate the cell wall structure. Our data indicate that these effects are mediated by the cleavage of the TIMP-like protein by the metzincin, and its subsequent release to the extracellular milieu. While the importance of this locus to Lm pathogenicity remains unclear, the observation that this phage-associated protein pair act upon the bacterial cell wall may hold promise in the field of antibiotic potentiation to combat antibiotic resistant bacterial pathogens.

Genetic interaction network has a very limited impact on the evolutionary trajectories in continuous culture-grown populations of yeast

Background

The impact of genetic interaction networks on evolution is a fundamental issue. Previous studies have demonstrated that the topology of the network is determined by the properties of the cellular machinery. Functionally related genes frequently interact with one another, and they establish modules, e.g., modules of protein complexes and biochemical pathways. In this study, we experimentally tested the hypothesis that compensatory evolutionary modifications, such as mutations and transcriptional changes, occur frequently in genes from perturbed modules of interacting genes.

Results

Using Saccharomyces cerevisiae haploid deletion mutants as a model, we investigated two modules lacking COG7 or NUP133, which are evolutionarily conserved genes with many interactions. We performed laboratory evolution experiments with these strains in two genetic backgrounds (with or without additional deletion of MSH2), subjecting them to continuous culture in a non-limiting minimal medium. Next, the evolved yeast populations were characterized through whole-genome sequencing and transcriptome analyses. No obvious compensatory changes resulting from inactivation of genes already included in modules were identified. The supposedly compensatory inactivation of genes in the evolved strains was only rarely observed to be in accordance with the established fitness effect of the genetic interaction network. In fact, a substantial majority of the gene inactivations were predicted to be neutral in the experimental conditions used to determine the interaction network. Similarly, transcriptome changes during continuous culture mostly signified adaptation to growth conditions rather than compensation of the absence of the COG7, NUP133 or MSH2 genes. However, we noticed that for genes whose inactivation was deleterious an upregulation of transcription was more common than downregulation.

Conclusions

Our findings demonstrate that the genetic interactions and the modular structure of the network described by others have very limited effects on the evolutionary trajectory following gene deletion of module elements in our experimental conditions and has no significant impact on short-term compensatory evolution. However, we observed likely compensatory evolution in functionally related (albeit non-interacting) genes.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12862-021-01830-9.

Hidden paths to endless forms most wonderful: parasite-blind diversification of host quality

Evolutionary diversification can occur in allopatry or sympatry, can be driven by selection or unselected, and can be phenotypically manifested immediately or remain latent until manifested in a newly encountered environment. Diversification of host-parasite interactions is frequently studied in the context of intrinsically selective coevolution, but the potential for host-parasite interaction phenotypes to diversify latently during parasite-blind host evolution is rarely considered. Here, we use a social bacterium experimentally adapted to several environments in the absence of phage to analyse allopatric diversification of host quality-the degree to which a host population supports a viral epidemic. Phage-blind evolution reduced host quality overall, with some bacteria becoming completely resistant to growth suppression by phage. Selective-environment differences generated only mild divergence in host quality. However, selective environments nonetheless played a major role in shaping evolution by determining the degree of stochastic diversification among replicate populations within treatments. Ancestral motility genotype was also found to strongly shape patterns of latent host-quality evolution and diversification. These outcomes show that (i) adaptive landscapes can differ in how they constrain stochastic diversification of a latent phenotype and (ii) major effects of selection on biological diversification can be missed by focusing on trait means. Collectively, our findings suggest that latent-phenotype evolution should inform host-parasite evolution theory and that diversification should be conceived broadly to include latent phenotypes.

Accessory Genome Dynamics and Structural Variation of Shigella from Persistent Infections

Shigellosis is a diarrheal disease caused mainly by Shigella flexneri and Shigella sonnei Infection is thought to be largely self-limiting, with short- to medium-term and serotype-specific immunity provided following clearance. However, cases of men who have sex with men (MSM)-associated shigellosis have been reported where Shigella of the same serotype were serially sampled from individuals between 1 and 1,862 days apart, possibly due to persistent carriage or reinfection with the same serotype. Here, we investigate the accessory genome dynamics of MSM-associated S. flexneri and S. sonnei isolates serially sampled from individual patients at various days apart to shed light on the adaptation of these important pathogens during infection. We find that pairs likely associated with persistent infection/carriage and with a smaller single nucleotide polymorphism (SNP) distance, demonstrated significantly less variation in accessory genome content than pairs likely associated with reinfection, and with a greater SNP distance. We observed antimicrobial resistance acquisition during Shigella carriage, including the gain of an extended-spectrum beta-lactamase gene during carriage. Finally, we explored large chromosomal structural variations and rearrangements in seven (five chronic and two reinfection associated) pairs of S. flexneri 3a isolates from an MSM-associated epidemic sublineage, which revealed variations at several common regions across isolate pairs, mediated by insertion sequence elements and comprising a distinct predicted functional profile. This study provides insight on the variation of accessory genome dynamics and large structural genomic changes in Shigella during persistent infection/carriage. In addition, we have also created a complete reference genome and biobanked isolate of the globally important pathogen, S. flexneri 3a.IMPORTANCEShigella spp. are Gram-negative bacteria that are the etiological agent of shigellosis, the second most common cause of diarrheal illness among children under the age of five in low-income countries. In high-income countries, shigellosis is also a sexually transmissible disease among men who have sex with men. Within the latter setting, we have captured prolonged and/or recurrent infection with shigellae of the same serotype, challenging the belief that Shigella infection is short lived and providing an early opportunity to study the evolution of the pathogen over the course of infection. Using this recently emerged transmission scenario, we comprehensively characterize the genomic changes that occur over the course of individual infection with Shigella and uncover a distinct functional profile of variable genomic regions, findings that have relevance for other Enterobacteriaceae.

Genetic Context Significantly Influences the Maintenance and Evolution of Degenerate Pathways

Understanding the evolution of novel physiological traits is highly relevant for expanding the characterization and manipulation of biological systems. Acquisition of new traits can be achieved through horizontal gene transfer (HGT). Here, we investigate drivers that promote or deter the maintenance of HGT-driven degeneracy, occurring when processes accomplish identical functions through nonidentical components. Subsequent evolution can optimize newly acquired functions; for example, beneficial alleles identified in an engineered Methylorubrum extorquens strain allowed it to utilize a “Foreign” formaldehyde oxidation pathway substituted for its Native pathway for methylotrophic growth. We examined the fitness consequences of interactions between these alleles when they were combined with the Native pathway or both (Dual) pathways. Unlike the Foreign pathway context where they evolved, these alleles were often neutral or deleterious when moved into these alternative genetic backgrounds. However, there were instances where combinations of multiple alleles resulted in higher fitness outcomes than individual allelic substitutions could provide. Importantly, the genetic context accompanying these allelic substitutions significantly altered the fitness landscape, shifting local fitness peaks and restricting the set of accessible evolutionary trajectories. These findings highlight how genetic context can negatively impact the probability of maintaining native and HGT-introduced functions together, making it difficult for degeneracy to evolve. However, in cases where the cost of maintaining degeneracy was mitigated by adding evolved alleles impacting the function of these pathways, we observed rare opportunities for pathway coevolution to occur. Together, our results highlight the importance of genetic context and resulting epistasis in retaining or losing HGT-acquired degenerate functions.

The enigmatic biology of rickettsiae: recent advances, open questions and outlook

Rickettsiae are obligate intracellular bacteria that can cause life-threatening illnesses and are among the oldest known vector-borne pathogens. Members of this genus are extraordinarily diverse and exhibit a broad host range. To establish intracellular infection, Rickettsia species undergo complex, multistep life cycles that are encoded by heavily streamlined genomes. As a result of reductive genome evolution, rickettsiae are exquisitely tailored to their host cell environment but cannot survive extracellularly. This host-cell dependence makes for a compelling system to uncover novel host-pathogen biology, but it has also hindered experimental progress. Consequently, the molecular details of rickettsial biology and pathogenesis remain poorly understood. With recent advances in molecular biology and genetics, the field is poised to start unraveling the molecular mechanisms of these host-pathogen interactions. Here, we review recent discoveries that have shed light on key aspects of rickettsial biology. These studies have revealed that rickettsiae subvert host cells using mechanisms that are distinct from other better-studied pathogens, underscoring the great potential of the Rickettsia genus for revealing novel biology. We also highlight several open questions as promising areas for future study and discuss the path toward solving the fundamental mysteries of this neglected and emerging human pathogen.

Evidence of gene nucleotide composition favoring replication and growth in a fastidious plant pathogen

Nucleotide composition (GC content) varies across bacteria species, genome regions, and specific genes. In Xylella fastidiosa, a vector-borne fastidious plant pathogen infecting multiple crops, GC content ranges between ∼51-52%; however, these values were gathered using limited genomic data. We evaluated GC content variations across X. fastidiosa subspecies fastidiosa (N = 194), subsp. pauca (N = 107), and subsp. multiplex (N = 39). Genomes were classified based on plant host and geographic origin; individual genes within each genome were classified based on gene function, strand, length, ortholog group, Core vs. Accessory, and Recombinant vs. Non-recombinant. GC content was calculated for each gene within each evaluated genome. The effects of genome and gene level variables were evaluated with a mixed effect ANOVA, and the marginal-GC content was calculated for each gene. Also, the correlation between gene-specific GC content vs. natural selection (dN/dS) and recombination/mutation (r/m) was estimated. Our analyses show that intra-genomic changes in nucleotide composition in X. fastidiosa are small and influenced by multiple variables. Higher AT-richness is observed in genes involved in replication and translation, and genes in the leading strand. In addition, we observed a negative correlation between high-AT and dN/dS in subsp. pauca. The relationship between recombination and GC content varied between core and accessory genes. We hypothesize that distinct evolutionary forces and energetic constraints both drive and limit these small variations in nucleotide composition.

Inference of symbiotic adaptations in nature using experimental evolution

Microbes must adapt to the presence of other species, but it can be difficult to recreate the natural context for these interactions in the laboratory. We describe a method for inferring the existence of symbiotic adaptations by experimentally evolving microbes that would normally interact in an artificial environment without access to other species. By looking for changes in the fitness effects microbes adapted to isolation have on their partners, we can infer the existence of ancestral adaptations that were lost during experimental evolution. The direction and magnitude of trait changes can offer useful insight as to whether the microbes have historically been selected to help or harm one another in nature. We apply our method to the complex symbiosis between the social amoeba Dictyostelium discoideum and two intracellular bacterial endosymbionts, Paraburkholderia agricolaris and Paraburkholderia hayleyella. Our results suggest P. hayleyella-but not P. agricolaris-has generally been selected to attenuate its virulence in nature, and that D. discoideum has evolved to antagonistically limit the growth of Paraburkholderia. The approach demonstrated here can be a powerful tool for studying adaptations in microbes, particularly when the specific natural context in which the adaptations evolved is unknown or hard to reproduce.

Mechanisms That Shape Microbial Pangenomes

Analyses of multiple whole-genome sequences from the same species have revealed that differences in gene content can be substantial, particularly in prokaryotes. Such variation has led to the recognition of pangenomes, the complete set of genes present in a species - consisting of core genes, present in all individuals, and accessory genes whose presence is variable. Questions now arise about how pangenomes originate and evolve. We describe how gene content variation can arise as a result of the combination of several processes, including random drift, selection, gain/loss balance, and the influence of ecological and epistatic interactions. We believe that identifying the contributions of these processes to pangenomes will need novel theoretical approaches and empirical data.

Reflections on the Predictability of Evolution: Toward a Conceptual Framework

Evolution is generally considered to be unpredictable because genetic variations are known to occur randomly. However, remarkable patterns of repeated convergent evolution are observed, for instance, loss of pigments by organisms living in caves. Analogous phenotypes appear in similar environments, sometimes in response to similar constraints. Alongside randomness, a certain evolutionary determinism also exists, for instance, the selection of particular phenotypes subjected to particular environmental constraints in the “evolutionary funnel.” We pursue the idea that eco-evolutionary specialization is in some way determinist. The conceptual framework of phenotypic changes entailing specialization presented in this essay explains how evolution can be predicted. We also discuss how the predictability of evolution could be tested using the case of metabolic specialization through gene losses. We also put forward that microorganisms could be key models to test and possibly make headway evolutionary predictions and knowledge about evolution.

CRISPR-CasΦ from huge phages is a hypercompact genome editor

CRISPR-Cas systems are found widely in prokaryotes, where they provide adaptive immunity against virus infection and plasmid transformation. We describe a minimal functional CRISPR-Cas system, comprising a single ~70-kilodalton protein, CasΦ, and a CRISPR array, encoded exclusively in the genomes of huge bacteriophages. CasΦ uses a single active site for both CRISPR RNA (crRNA) processing and crRNA-guided DNA cutting to target foreign nucleic acids. This hypercompact system is active in vitro and in human and plant cells with expanded target recognition capabilities relative to other CRISPR-Cas proteins. Useful for genome editing and DNA detection but with a molecular weight half that of Cas9 and Cas12a genome-editing enzymes, CasΦ offers advantages for cellular delivery that expand the genome editing toolbox.

Dynamics of Insect-Microbiome Interaction Influence Host and Microbial Symbiont

Insects share an intimate relationship with their gut microflora and this symbiotic association has developed into an essential evolutionary outcome intended for their survival through extreme environmental conditions. While it has been clearly established that insects, with very few exceptions, associate with several microbes during their life cycle, information regarding several aspects of these associations is yet to be fully unraveled. Acquisition of bacteria by insects marks the onset of microbial symbiosis, which is followed by the adaptation of these bacterial species to the gut environment for prolonged sustenance and successful transmission across generations. Although several insect-microbiome associations have been reported and each with their distinctive features, diversifications and specializations, it is still unclear as to what led to these diversifications. Recent studies have indicated the involvement of various evolutionary processes operating within an insect body that govern the transition of a free-living microbe to an obligate or facultative symbiont and eventually leading to the establishment and diversification of these symbiotic relationships. Data from various studies, summarized in this review, indicate that the symbiotic partners, i.e., the bacteria and the insect undergo several genetic, biochemical and physiological changes that have profound influence on their life cycle and biology. An interesting outcome of the insect-microbe interaction is the compliance of the microbial partner to its eventual genome reduction. Endosymbionts possess a smaller genome as compared to their free-living forms, and thus raising the question what is leading to reductive evolution in the microbial partner. This review attempts to highlight the fate of microbes within an insect body and its implications for both the bacteria and its insect host. While discussion on each specific association would be too voluminous and outside the scope of this review, we present an overview of some recent studies that contribute to a better understanding of the evolutionary trajectory and dynamics of the insect-microbe association and speculate that, in the future, a better understanding of the nature of this interaction could pave the path to a sustainable and environmentally safe way for controlling economically important pests of crop plants.

Large-scale and significant expression from pseudogenes in Sodalis glossinidius - a facultative bacterial endosymbiont

The majority of bacterial genomes have high coding efficiencies, but there are some genomes of intracellular bacteria that have low gene density. The genome of the endosymbiont Sodalis glossinidius contains almost 50 % pseudogenes containing mutations that putatively silence them at the genomic level. We have applied multiple 'omic' strategies, combining Illumina and Pacific Biosciences Single-Molecule Real-Time DNA sequencing and annotation, stranded RNA sequencing and proteome analysis to better understand the transcriptional and translational landscape of Sodalis pseudogenes, and potential mechanisms for their control. Between 53 and 74 % of the Sodalis transcriptome remains active in cell-free culture. The mean sense transcription from coding domain sequences (CDSs) is four times greater than that from pseudogenes. Comparative genomic analysis of six Illumina-sequenced Sodalis isolates from different host Glossina species shows pseudogenes make up ~40 % of the 2729 genes in the core genome, suggesting that they are stable and/or that Sodalis is a recent introduction across the genus Glossina as a facultative symbiont. These data shed further light on the importance of transcriptional and translational control in deciphering host-microbe interactions. The combination of genomics, transcriptomics and proteomics gives a multidimensional perspective for studying prokaryotic genomes with a view to elucidating evolutionary adaptation to novel environmental niches.

Destabilization of the Tumor-Inducing Plasmid from an Octopine-Type Agrobacterium tumefaciens Lineage Drives a Large Deletion in the Co-resident At Megaplasmid

Bacteria with multi-replicon genome organizations, including members of the family Rhizobiaceae, often carry a variety of niche-associated functions on large plasmids. While evidence exists for cross-replicon interactions and co-evolution between replicons in many of these systems, remarkable strain-to-strain variation is also observed for extrachromosomal elements, suggesting increased genetic plasticity. Here, we show that curing of the tumor-inducing virulence plasmid (pTi) of an octopine-type Agrobacterium tumefaciens lineage leads to a large deletion in the co-resident At megaplasmid (pAt). The deletion event is mediated by a repetitive IS-element, IS66, and results in a variety of environment-dependent fitness consequences, including loss of independent conjugal transfer of the plasmid. Interestingly, a related and otherwise wild-type A. tumefaciens strain is missing exactly the same large pAt segment as the pAt deletion derivatives, suggesting a similar event over its natural history. Overall, the findings presented here uncover a novel genetic interaction between the two large plasmids of A. tumefaciens and provide evidence for cross-replicon integration and co-evolution of these plasmids.

Microbial Experimental Evolution - a proving ground for evolutionary theory and a tool for discovery

Microbial experimental evolution uses controlled laboratory populations to study the mechanisms of evolution. The molecular analysis of evolved populations enables empirical tests that can confirm the predictions of evolutionary theory, but can also lead to surprising discoveries. As with other fields in the life sciences, microbial experimental evolution has become a tool, deployed as part of the suite of techniques available to the molecular biologist. Here, I provide a review of the general findings of microbial experimental evolution, especially those relevant to molecular microbiologists that are new to the field. I also relate these results to design considerations for an evolution experiment and suggest future directions for those working at the intersection of experimental evolution and molecular biology.

Pseudomonas aeruginosa MutL promotes large chromosomal deletions through non-homologous end joining to prevent bacteriophage predation

Pseudomonas aeruginosa is an opportunistic pathogen with a relatively large genome, and has been shown to routinely lose genomic fragments during environmental selection. However, the underlying molecular mechanisms that promote chromosomal deletion are still poorly understood. In a recent study, we showed that by deleting a large chromosomal fragment containing two closely situated genes, hmgA and galU, P. aeruginosa was able to form ‘brown mutants’, bacteriophage (phage) resistant mutants with a brown color phenotype. In this study, we show that the brown mutants occur at a frequency of 227 ± 87 × 10⁻⁸ and contain a deletion ranging from ∼200 to ∼620 kb. By screening P. aeruginosa transposon mutants, we identified mutL gene whose mutation constrained the emergence of phage-resistant brown mutants. Moreover, the P. aeruginosa MutL (PaMutL) nicking activity can result in DNA double strand break (DSB), which is then repaired by non-homologous end joining (NHEJ), leading to chromosomal deletions. Thus, we reported a noncanonical function of PaMutL that promotes chromosomal deletions through NHEJ to prevent phage predation.

Defining Division of Labor in Microbial Communities

In order to survive and reproduce, organisms must perform a multitude of tasks. However, trade-offs limit their ability to allocate energy and resources to all of these different processes. One strategy to solve this problem is to specialize in some traits and team up with other organisms that can help by providing additional, complementary functions. By reciprocally exchanging metabolites and/or services in this way, both parties benefit from the interaction. This phenomenon, which has been termed functional specialization or division of labor, is very common in nature and exists on all levels of biological organization. Also, microorganisms have evolved different types of synergistic interactions. However, very often, it remains unclear whether or not a given example represents a true case of division of labor. Here we aim at filling this gap by providing a list of criteria that clearly define division of labor in microbial communities. Furthermore, we propose a set of diagnostic experiments to verify whether a given interaction fulfills these conditions. In contrast to the common use of the term, our analysis reveals that both intraspecific and interspecific interactions meet the criteria defining division of labor. Moreover, our analysis identified non-cooperators of intraspecific public goods interactions as growth specialists that divide labor with conspecific producers, rather than being social parasites. By providing a conceptual toolkit, our work will help to unambiguously identify cases of division of labor and stimulate more detailed investigations of this important and widespread type of inter-microbial interaction.

Predictive genomic traits for bacterial growth in culture versus actual growth in soil

Relationships between microbial genes and performance are often evaluated in the laboratory in pure cultures, with little validation in nature. Here, we show that genomic traits related to laboratory measurements of maximum growth potential failed to predict the growth rates of bacteria in unamended soil, but successfully predicted growth responses to resource pulses: growth increased with 16S rRNA gene copy number and declined with genome size after substrate addition to soils, responses that were repeated in four different ecosystems. Genome size best predicted growth rate in response to addition of glucose alone; adding ammonium with glucose weakened the relationship, and the relationship was absent in nutrient-replete pure cultures, consistent with the idea that reduced genome size is a mechanism of nutrient conservation. Our findings demonstrate that genomic traits of soil bacteria can map to their ecological performance in nature, but the mapping is poor under native soil conditions, where genomic traits related to stress tolerance may prove more predictive. These results remind that phenotype depends on environmental context, underscoring the importance of verifying proposed schemes of trait-based strategies through direct measurement of performance in nature, an important and currently missing foundation for translating microbial processes from genes to ecosystems.

Predictive genomic traits for bacterial growth in culture versus actual growth in soil

Relationships between microbial genes and performance are often evaluated in the laboratory in pure cultures, with little validation in nature. Here, we show that genomic traits related to laboratory measurements of maximum growth potential failed to predict the growth rates of bacteria in unamended soil, but successfully predicted growth responses to resource pulses: growth increased with 16S rRNA gene copy number and declined with genome size after substrate addition to soils, responses that were repeated in four different ecosystems. Genome size best predicted growth rate in response to addition of glucose alone; adding ammonium with glucose weakened the relationship, and the relationship was absent in nutrient-replete pure cultures, consistent with the idea that reduced genome size is a mechanism of nutrient conservation. Our findings demonstrate that genomic traits of soil bacteria can map to their ecological performance in nature, but the mapping is poor under native soil conditions, where genomic traits related to stress tolerance may prove more predictive. These results remind that phenotype depends on environmental context, underscoring the importance of verifying proposed schemes of trait-based strategies through direct measurement of performance in nature, an important and currently missing foundation for translating microbial processes from genes to ecosystems.

Genomic and Transcriptomic Changes that Mediate Increased Platinum Resistance in Cupriavidus metallidurans

The extensive anthropogenic use of platinum, a rare element found in low natural abundance in the Earth's continental crust and one of the critical raw materials in the EU innovation partnership framework, has resulted in increased concentrations in surface environments. To minimize its spread and increase its recovery from the environment, biological recovery via different microbial systems is explored. In contrast, studies focusing on the effects of prolonged exposure to Pt are limited. In this study, we used the metal-resistant Cupriavidus metallidurans NA4 strain to explore the adaptation of environmental bacteria to platinum exposure. We used a combined Nanopore⁻Illumina sequencing approach to fully resolve all six replicons of the C. metallidurans NA4 genome, and compared them with the C. metallidurans CH34 genome, revealing an important role in metal resistance for its chromid rather than its megaplasmids. In addition, we identified the genomic and transcriptomic changes in a laboratory-evolved strain, displaying resistance to 160 µM Pt4+. The latter carried 20 mutations, including a large 69.9 kb deletion in its plasmid pNA4_D (89.6 kb in size), and 226 differentially-expressed genes compared to its parental strain. Many membrane-related processes were affected, including up-regulation of cytochrome c and a lytic transglycosylase, down-regulation of flagellar and pili-related genes, and loss of the pNA4_D conjugative machinery, pointing towards a significant role in the adaptation to platinum.

Predicting mutational routes to new adaptive phenotypes

Predicting evolutionary change poses numerous challenges. Here we take advantage of the model bacterium Pseudomonas fluorescens in which the genotype-to-phenotype map determining evolution of the adaptive 'wrinkly spreader' (WS) type is known. We present mathematical descriptions of three necessary regulatory pathways and use these to predict both the rate at which each mutational route is used and the expected mutational targets. To test predictions, mutation rates and targets were determined for each pathway. Unanticipated mutational hotspots caused experimental observations to depart from predictions but additional data led to refined models. A mismatch was observed between the spectra of WS-causing mutations obtained with and without selection due to low fitness of previously undetected WS-causing mutations. Our findings contribute toward the development of mechanistic models for forecasting evolution, highlight current limitations, and draw attention to challenges in predicting locus-specific mutational biases and fitness effects.

Metabolic adaptations underlying genome flexibility in prokaryotes

Even across genomes of the same species, prokaryotes exhibit remarkable flexibility in gene content. We do not know whether this flexible or "accessory" content is mostly neutral or adaptive, largely due to the lack of explicit analyses of accessory gene function. Here, across 96 diverse prokaryotic species, I show that a considerable fraction (~40%) of accessory genomes harbours beneficial metabolic functions. These functions take two forms: (1) they significantly expand the biosynthetic potential of individual strains, and (2) they help reduce strain-specific metabolic auxotrophies via intra-species metabolic exchanges. I find that the potential of both these functions increases with increasing genome flexibility. Together, these results are consistent with a significant adaptive role for prokaryotic pangenomes.

Experimental Design, Population Dynamics, and Diversity in Microbial Experimental Evolution

In experimental evolution, laboratory-controlled conditions select for the adaptation of species, which can be monitored in real time. Despite the current popularity of such experiments, nature's most pervasive biological force was long believed to be observable only on time scales that transcend a researcher's life-span, and studying evolution by natural selection was therefore carried out solely by comparative means. Eventually, microorganisms' propensity for fast evolutionary changes proved us wrong, displaying strong evolutionary adaptations over a limited time, nowadays massively exploited in laboratory evolution experiments. Here, we formulate a guide to experimental evolution with microorganisms, explaining experimental design and discussing evolutionary dynamics and outcomes and how it is used to assess ecoevolutionary theories, improve industrially important traits, and untangle complex phenotypes. Specifically, we give a comprehensive overview of the setups used in experimental evolution. Additionally, we address population dynamics and genetic or phenotypic diversity during evolution experiments and expand upon contributing factors, such as epistasis and the consequences of (a)sexual reproduction. Dynamics and outcomes of evolution are most profoundly affected by the spatiotemporal nature of the selective environment, where changing environments might lead to generalists and structured environments could foster diversity, aided by, for example, clonal interference and negative frequency-dependent selection. We conclude with future perspectives, with an emphasis on possibilities offered by fast-paced technological progress. This work is meant to serve as an introduction to those new to the field of experimental evolution, as a guide to the budding experimentalist, and as a reference work to the seasoned expert.

Genome analysis of the ubiquitous boxwood pathogen Pseudonectria foliicola

Boxwood (Buxus spp.) are broad-leaved, evergreen landscape plants valued for their longevity and ornamental qualities. Volutella leaf and stem blight, caused by the ascomycete fungi Pseudonectria foliicola and P. buxi, is one of the major diseases affecting the health and ornamental qualities of boxwood. Although this disease is less severe than boxwood blight caused by Calonectria pseudonaviculata and C. henricotiae, its widespread occurrence and disfiguring symptoms have caused substantial economic losses to the ornamental industry. In this study, we sequenced the genome of P. foliicola isolate ATCC13545 using Illumina technology and compared it to other publicly available fungal pathogen genomes to better understand the biology of this organism. A de novo assembly estimated the genome size of P. foliicola at 28.7 Mb (425 contigs; N50 = 184,987 bp; avg. coverage 188×), with just 9,272 protein-coding genes. To our knowledge, P. foliicola has the smallest known genome within the Nectriaceae. Consistent with the small size of the genome, the secretome, CAzyme and secondary metabolite profiles of this fungus are reduced relative to two other surveyed Nectriaceae fungal genomes: Dactylonectria macrodidyma JAC15-245 and Fusarium graminearum Ph-1. Interestingly, a large cohort of genes associated with reduced virulence and loss of pathogenicity was identified from the P. foliicola dataset. These data are consistent with the latest observations by plant pathologists that P. buxi and most likely P. foliicola, are opportunistic, latent pathogens that prey upon weak and stressed boxwood plants.