Functional Constraints on Replacing an Essential Gene with Its Ancient and Modern Homologs

Fitness benefits of a synonymous substitution in an ancient EF-Tu gene depend on the genetic background

Synonymous mutations are changes to DNA sequence, which occur within translated genes but which do not affect the protein sequence. Although often referred to as silent mutations, evidence suggests that synonymous mutations can affect gene expression, mRNA stability, and even translation efficiency. A collection of both experimental and bioinformatic data has shown that synonymous mutations can impact cell phenotype, yet less is known about the molecular mechanisms and potential of beneficial or adaptive effects of such changes within evolved populations. Here, we report a beneficial synonymous mutation acquired via experimental evolution in an essential gene variant encoding the translation elongation factor protein EF-Tu. We demonstrate that this particular synonymous mutation increases EF-Tu mRNA and protein levels as well as global polysome abundance on RNA transcripts. Although presence of the synonymous mutation is clearly causative of such changes, we also demonstrate that fitness benefits are highly contingent on other potentiating mutations present within the genetic background in which the mutation arose. Our results underscore the importance of beneficial synonymous mutations, especially those that affect levels of proteins that are key for cellular processes.IMPORTANCEThis study explores the degree to which synonymous mutations in essential genes can influence adaptation in bacteria. An experimental system whereby an Escherichia coli strain harboring an engineered translation protein elongation factor-Tu (EF-Tu) was subjected to laboratory evolution. We find that a synonymous mutation acquired on the gene encoding for EF-Tu is conditionally beneficial for bacterial fitness. Our findings provide insight into the importance of the genetic background when a synonymous substitution is favored by natural selection and how such changes have the potential to impact evolution when critical cellular processes are involved.

A beneficial synonymous substitution in EF-Tu is contingent on genetic background

Synonymous mutations are changes to DNA sequence that occur within translated genes but which do not affect the protein sequence. Although often referred to as silent mutations, evidence suggests that synonymous mutations can affect gene expression, mRNA stability, and even translation efficiency. A collection of both experimental and bioinformatic data has shown that synonymous mutations can impact cell phenotype, yet less is known about the molecular mechanisms and potential of beneficial or adaptive effects of such changes within evolved populations. Here, we report a beneficial synonymous mutation acquired via experimental evolution in an essential gene variant encoding the translation Elongation Factor protein EF-Tu. We demonstrate that this particular synonymous mutation increases EF-Tu mRNA and protein levels, as well as the polysome abundance on global transcripts. Although presence of the synonymous mutation is clearly causative of such changes, we also demonstrate that fitness benefits are highly contingent on other potentiating mutations present within the genetic background in which the mutation arose. Our results underscore the importance of beneficial synonymous mutations, especially those that affect levels of proteins that are key for cellular processes.

Nitrogenase resurrection and the evolution of a singular enzymatic mechanism

The planetary biosphere is powered by a suite of key metabolic innovations that emerged early in the history of life. However, it is unknown whether life has always followed the same set of strategies for performing these critical tasks. Today, microbes access atmospheric sources of bioessential nitrogen through the activities of just one family of enzymes, nitrogenases. Here, we show that the only dinitrogen reduction mechanism known to date is an ancient feature conserved from nitrogenase ancestors. We designed a paleomolecular engineering approach wherein ancestral nitrogenase genes were phylogenetically reconstructed and inserted into the genome of the diazotrophic bacterial model, Azotobacter vinelandii, enabling an integrated assessment of both in vivo functionality and purified nitrogenase biochemistry. Nitrogenase ancestors are active and robust to variable incorporation of one or more ancestral protein subunits. Further, we find that all ancestors exhibit the reversible enzymatic mechanism for dinitrogen reduction, specifically evidenced by hydrogen inhibition, which is also exhibited by extant A. vinelandii nitrogenase isozymes. Our results suggest that life may have been constrained in its sampling of protein sequence space to catalyze one of the most energetically challenging biochemical reactions in nature. The experimental framework established here is essential for probing how nitrogenase functionality has been shaped within a dynamic, cellular context to sustain a globally consequential metabolism.

Incompatibility and Interchangeability in Molecular Evolution

There is remarkable variation in the rate at which genetic incompatibilities in molecular interactions accumulate. In some cases, minor changes-even single-nucleotide substitutions-create major incompatibilities when hybridization forces new variants to function in a novel genetic background from an isolated population. In other cases, genes or even entire functional pathways can be horizontally transferred between anciently divergent evolutionary lineages that span the tree of life with little evidence of incompatibilities. In this review, we explore whether there are general principles that can explain why certain genes are prone to incompatibilities while others maintain interchangeability. We summarize evidence pointing to four genetic features that may contribute to greater resistance to functional replacement: (1) function in multisubunit enzyme complexes and protein-protein interactions, (2) sensitivity to changes in gene dosage, (3) rapid rate of sequence evolution, and (4) overall importance to cell viability, which creates sensitivity to small perturbations in molecular function. We discuss the relative levels of support for these different hypotheses and lay out future directions that may help explain the striking contrasts in patterns of incompatibility and interchangeability throughout the history of molecular evolution.

Very early evolution from the perspective of microbial ecology

The universal ancestor at the root of the species tree of life depicts a population of organisms with a surprising degree of complexity, posessing genomes and translation systems much like that of microbial life today. As the first life forms were most likely to have been simple replicators, considerable evolutionary change must have taken place prior to the last universal common ancestor. It is often assumed that the lack of earlier branches on the tree of life is due to a prevalence of random horizontal gene transfer that obscured the delineations between lineages and hindered their divergence. Therefore, principles of microbial evolution and ecology may give us some insight into these early stages in the history of life. Here, we synthesize the current understanding of organismal and genome evolution from the perspective of microbial ecology and apply these evolutionary principles to the earliest stages of life on Earth. We focus especially on broad evolutionary modes pertaining to horizontal gene transfer, pangenome structure, and microbial mat communities.

Low levels of tetracyclines select for a mutation that prevents the evolution of high-level resistance to tigecycline

In a collection of Escherichia coli isolates, we discovered a new mechanism leading to frequent and high-level tigecycline resistance involving tandem gene amplifications of an efflux pump encoded by the tet(A) determinant. Some isolates, despite carrying a functional tet(A), could not evolve high-level tigecycline resistance by amplification due to the presence of a deletion in the TetR(A) repressor. This mutation impaired induction of tetA(A) (encoding the TetA(A) efflux pump) in presence of tetracyclines, with the strongest effect observed for tigecycline, subsequently preventing the development of tet(A) amplification-dependent high-level tigecycline resistance. We found that this mutated tet(A) determinant was common among tet(A)-carrying E. coli isolates and analysed possible explanations for this high frequency. First, while the mutated tet(A) was found in several ST-groups, we found evidence of clonal spread among ST131 isolates, which increases its frequency within E. coli databases. Second, evolution and competition experiments revealed that the mutation in tetR(A) could be positively selected over the wild-type allele at sub-inhibitory concentrations of tetracyclines. Our work demonstrates how low concentrations of tetracyclines, such as those found in contaminated environments, can enrich and select for a mutation that generates an evolutionary dead-end that precludes the evolution towards high-level, clinically relevant tigecycline resistance.

A study on evolution of antimicrobial resistance reveals how sub-inhibitory concentrations of antibiotics enrich and select for a mutated allele that prevents evolution towards clinically significant levels of antibiotic resistance.

Gradients in gene essentiality reshape antibacterial research

Essential genes encode the processes that are necessary for life. Until recently, commonly applied binary classifications left no space between essential and non-essential genes. In this review, we frame bacterial gene essentiality in the context of genetic networks. We explore how the quantitative properties of gene essentiality are influenced by the nature of the encoded process, environmental conditions and genetic background, including a strain's distinct evolutionary history. The covered topics have important consequences for antibacterials, which inhibit essential processes. We argue that the quantitative properties of essentiality can thus be used to prioritize antibacterial cellular targets and desired spectrum of activity in specific infection settings. We summarize our points with a case study on the core essential genome of the cystic fibrosis pathobiome and highlight avenues for targeted antibacterial development.

The Minimal Translation Machinery: What We Can Learn From Naturally and Experimentally Reduced Genomes

The current theoretical proposals of minimal genomes have not attempted to outline the essential machinery for proper translation in cells. Here, we present a proposal of a minimal translation machinery based on (1) a comparative analysis of bacterial genomes of insects’ endosymbionts using a machine learning classification algorithm, (2) the empiric genomic information obtained from Mycoplasma mycoides JCVI-syn3.0 the first minimal bacterial genome obtained by design and synthesis, and (3) a detailed functional analysis of the candidate genes based on essentiality according to the DEG database (Escherichia coli and Bacillus subtilis) and the literature. This proposed minimal translational machinery is composed by 142 genes which must be present in any synthetic prokaryotic cell designed for biotechnological purposes, 76.8% of which are shared with JCVI-syn3.0. Eight additional genes were manually included in the proposal for a proper and efficient translation.

A positive correlation between GC content and growth temperature in prokaryotes

BACKGROUND

GC pairs are generally more stable than AT pairs; GC-rich genomes were proposed to be more adapted to high temperatures than AT-rich genomes. Previous studies consistently showed positive correlations between growth temperature and the GC contents of structural RNA genes. However, for the whole genome sequences and the silent sites of the codons in protein-coding genes, the relationship between GC content and growth temperature is in a long-lasting debate.

RESULTS

With a dataset much larger than previous studies (681 bacteria and 155 archaea with completely assembled genomes), our phylogenetic comparative analyses showed positive correlations between optimal growth temperature (Topt) and GC content both in bacterial and archaeal structural RNA genes and in bacterial whole genome sequences, chromosomal sequences, plasmid sequences, core genes, and accessory genes. However, in the 155 archaea, we did not observe a significant positive correlation of Topt with whole-genome GC content (GCw) or GC content at four-fold degenerate sites. We randomly drew 155 samples from the 681 bacteria for 1000 rounds. In most cases (> 95%), the positive correlations between Topt and genomic GC contents became statistically nonsignificant (P > 0.05). This result suggested that the small sample sizes might account for the lack of positive correlations between growth temperature and genomic GC content in the 155 archaea and the bacterial samples of previous studies. Comparing the GC content among four categories (psychrophiles/psychrotrophiles, mesophiles, thermophiles, and hyperthermophiles) also revealed a positive correlation between GCw and growth temperature in bacteria. By including the GCw of incompletely assembled genomes, we expanded the sample size of archaea to 303. Positive correlations between GCw and Topt appear especially after excluding the halophilic archaea whose GC contents might be strongly shaped by intense UV radiation.

CONCLUSIONS

This study explains the previous contradictory observations and ends a long debate. Prokaryotes growing in high temperatures have higher GC contents. Thermal adaptation is one possible explanation for the positive association. Meanwhile, we propose that the elevated efficiency of DNA repair in response to heat mutagenesis might have the by-product of increasing GC content like that happens in intracellular symbionts and marine bacterioplankton.

Kinetic Analysis Suggests Evolution of Ribosome Specificity in Modern Elongation Factor-Tus from "Generalist" Ancestors

It has been hypothesized that early enzymes are more promiscuous than their extant orthologs. Whether or not this hypothesis applies to the translation machinery, the oldest molecular machine of life, is not known. Efficient protein synthesis relies on a cascade of specific interactions between the ribosome and the translation factors. Here, using elongation factor-Tu (EF-Tu) as a model system, we have explored the evolution of ribosome specificity in translation factors. Employing presteady state fast kinetics using quench flow, we have quantitatively characterized the specificity of two sequence-reconstructed 1.3- to 3.3-Gy-old ancestral EF-Tus toward two unrelated bacterial ribosomes, mesophilic Escherichia coli and thermophilic Thermus thermophilus. Although the modern EF-Tus show clear preference for their respective ribosomes, the ancestral EF-Tus show similar specificity for diverse ribosomes. In addition, despite increase in the catalytic activity with temperature, the ribosome specificity of the thermophilic EF-Tus remains virtually unchanged. Our kinetic analysis thus suggests that EF-Tu proteins likely evolved from the catalytically promiscuous, “generalist” ancestors. Furthermore, compatibility of diverse ribosomes with the modern and ancestral EF-Tus suggests that the ribosomal core probably evolved before the diversification of the EF-Tus. This study thus provides important insights regarding the evolution of modern translation machinery.

Genetic Architecture and Fitness of Bacterial Interspecies Hybrids

Integration of a conjugative plasmid into a bacterial chromosome can promote the transfer of chromosomal DNA to other bacteria. Intraspecies chromosomal conjugation is believed responsible for creating the global pathogens Klebsiella pneumoniae ST258 and Escherichia coli ST1193. Interspecies conjugation is also possible but little is known about the genetic architecture or fitness of such hybrids. To study this, we generated by conjugation 14 hybrids of E. coli and Salmonella enterica. These species belong to different genera, diverged from a common ancestor >100 Ma, and share a conserved order of orthologous genes with ∼15% nucleotide divergence. Genomic analysis revealed that all but one hybrid had acquired a contiguous segment of donor E. coli DNA, replacing a homologous region of recipient Salmonella chromosome, and ranging in size from ∼100 to >4,000 kb. Recombination joints occurred in sequences with higher-than-average nucleotide identity. Most hybrid strains suffered a large reduction in growth rate, but the magnitude of this cost did not correlate with the length of foreign DNA. Compensatory evolution to ameliorate the cost of low-fitness hybrids pointed towards disruption of complex genetic networks as a cause. Most interestingly, 4 of the 14 hybrids, in which from 45% to 90% of the Salmonella chromosome was replaced with E. coli DNA, showed no significant reduction in growth fitness. These data suggest that the barriers to creating high-fitness interspecies hybrids may be significantly lower than generally appreciated with implications for the creation of novel species.

Ancestral Sequence Reconstructions of MotB Are Proton-Motile and Require MotA for Motility

The bacterial flagellar motor (BFM) is a nanomachine that rotates the flagellum to propel many known bacteria. The BFM is powered by ion transit across the cell membrane through the stator complex, a membrane protein. Different bacteria use various ions to run their BFM, but the majority of BFMs are powered by either proton (H+) or sodium (Na+) ions. The transmembrane (TM) domain of the B-subunit of the stator complex is crucial for ion selectivity, as it forms the ion channel in complex with TM3 and TM4 of the A-subunit. In this study, we reconstructed and engineered thirteen ancestral sequences of the stator B-subunit to evaluate the functional properties and ionic power source of the stator proteins at reconstruction nodes to evaluate the potential of ancestral sequence reconstruction (ASR) methods for stator engineering and to test specific motifs previously hypothesized to be involved in ion-selectivity. We found that all thirteen of our reconstructed ancient B-subunit proteins could assemble into functional stator complexes in combination with the contemporary Escherichia coli MotA-subunit to restore motility in stator deleted E. coli strains. The flagellar rotation of the thirteen ancestral MotBs was found to be Na+ independent which suggested that the F30/Y30 residue was not significantly correlated with sodium/proton phenotype, in contrast to what we had reported previously. Additionally, four among the thirteen reconstructed B-subunits were compatible with the A-subunit of Aquifex aeolicus and able to function in a sodium-independent manner. Overall, this work demonstrates the use of ancestral reconstruction to generate novel stators and quantify which residues are correlated with which ionic power source.

A fine-scale map of genome-wide recombination in divergent Escherichia coli population

Recombination is one of the most important molecular mechanisms of prokaryotic genome evolution, but its exact roles are still in debate. Here we try to infer genome-wide recombination within a species, utilizing a dataset of 149 complete genomes of Escherichia coli from diverse animal hosts and geographic origins, including 45 in-house sequenced with the single-molecular real-time platform. Two major clades identified based on physiological, clinical and ecological characteristics form distinct genetic lineages based on scarcity of interclade gene exchanges. By defining gene-based syntenies for genomic segments within and between the two clades, we build a fine-scale recombination map for this representative global E. coli population. The map suggests extensive within-clade recombination that often breaks physical linkages among individual genes but seldom interrupts the structure of genome organizational frameworks as well as primary metabolic portfolios supported by the framework integrity, possibly due to strong natural selection for both physiological compatibility and ecological fitness. In contrast, the between-clade recombination declines drastically when phylogenetic distance increases to the extent where a 10-fold reduction can be observed, establishing a firm genetic barrier between clades. Our empirical data suggest a critical role for such recombination events in the early stage of speciation where recombination rate is associated with phylogenetic distance in addition to sequence and gene variations. The extensive intraclade recombination binds sister strains into a quasisexual group and optimizes genes or alleles to streamline physiological activities, whereas the sharply declined interclade recombination split the population into clades adaptive to divergent ecological niches.

Evolutionary History of Bioessential Elements Can Guide the Search for Life in the Universe

Our understanding of life in the universe comes from one sample, life on Earth. Current and next-generation space missions will target exoplanets as well as planets and moons in our own solar system with the primary goal of detecting, interpreting and characterizing indications of possible biological activity. Thus, understanding life's fundamental characteristics is increasingly critical for detecting and interpreting potential biological signatures elsewhere in the universe. Astrobiologists have outlined the essential roles of carbon and water for life, but we have yet to decipher the rules governing the evolution of how living organisms use bioessential elements. Does the suite of life's essential chemical elements on Earth constitute only one possible evolutionary outcome? Are some elements so essential for biological functions that evolution will select for them despite low availability? How would this play out on other worlds that have different relative element abundances? When we look for life in the universe, or the conditions that could give rise to life, we must learn how to recognize it in extremely different chemical and environmental conditions from those on Earth. We argue that by exposing self-organizing biotic chemistries to different combinations of abiotic materials, and by mapping the evolutionary history of metalloenzyme biochemistry onto geological availabilities of metals, alternative element choices that are very different from life's present-day molecular structure might result. A greater understanding of the paleomolecular evolutionary history of life on Earth will create a predictive capacity for detecting and assessing life's existence on worlds where alternate evolutionary paths might have been taken.

Synthetic cross-phyla gene replacement and evolutionary assimilation of major enzymes

The ability of DNA to produce a functional protein even after transfer to a foreign host is of fundamental importance in both evolutionary biology and biotechnology, enabling horizontal gene transfer in the wild and heterologous expression in the lab. However, the influence of genetic particulars on DNA functionality in a new host is poorly understood, as are the evolutionary mechanisms of assimilation and refinement. Here, we describe an automation-enabled large-scale experiment wherein Escherichia coli strains were evolved in parallel after replacement of the genes pgi or tpiA with orthologous DNA from donor species spanning all domains of life, from humans to hyperthermophilic archaea. Via analysis of hundreds of clones evolved for 50,000+ cumulative generations across dozens of independent lineages, we show that orthogene-upregulating mutations can completely mitigate fitness defects that result from initial non-functionality, with coding sequence changes unnecessary. Gene target, donor species and genomic location of the swap all influenced outcomes-both the nature of adaptive mutations (often synonymous) and the frequency with which strains successfully evolved to assimilate the foreign DNA. Additionally, time series DNA sequencing and replay evolution experiments revealed transient copy number expansions, the contingency of lineage outcome on first-step mutations and the ability for strains to escape from suboptimal local fitness maxima. Overall, this study establishes the influence of various DNA and protein features on cross-species genetic interchangeability and evolutionary outcomes, with implications for both horizontal gene transfer and rational strain design.

Evolutionary stalling and a limit on the power of natural selection to improve a cellular module

Cells consist of molecular modules which perform vital biological functions. Cellular modules are key units of adaptive evolution because organismal fitness depends on their performance. Theory shows that in rapidly evolving populations, such as those of many microbes, adaptation is driven primarily by common beneficial mutations with large effects, while other mutations behave as if they are effectively neutral. As a consequence, if a module can be improved only by rare and/or weak beneficial mutations, its adaptive evolution would stall. However, such evolutionary stalling has not been empirically demonstrated, and it is unclear to what extent stalling may limit the power of natural selection to improve modules. Here we empirically characterize how natural selection improves the translation machinery (TM), an essential cellular module. We experimentally evolved populations of Escherichia coli with genetically perturbed TMs for 1,000 generations. Populations with severe TM defects initially adapted via mutations in the TM, but TM adaptation stalled within about 300 generations. We estimate that the genetic load in our populations incurred by residual TM defects ranges from 0.5 to 19%. Finally, we found evidence that both epistasis and the depletion of the pool of beneficial mutations contributed to evolutionary stalling. Our results suggest that cellular modules may not be fully optimized by natural selection despite the availability of adaptive mutations.

Emerging Frontiers in the Study of Molecular Evolution

A collection of the editors of Journal of Molecular Evolution have gotten together to pose a set of key challenges and future directions for the field of molecular evolution. Topics include challenges and new directions in prebiotic chemistry and the RNA world, reconstruction of early cellular genomes and proteins, macromolecular and functional evolution, evolutionary cell biology, genome evolution, molecular evolutionary ecology, viral phylodynamics, theoretical population genomics, somatic cell molecular evolution, and directed evolution. While our effort is not meant to be exhaustive, it reflects research questions and problems in the field of molecular evolution that are exciting to our editors.

The SNAP hypothesis: Chromosomal rearrangements could emerge from positive Selection during Niche Adaptation

The relative linear order of most genes on bacterial chromosomes is not conserved over evolutionary timescales. One explanation is that selection is weak, allowing recombination to randomize gene order by genetic drift. However, most chromosomal rearrangements are deleterious to fitness. In contrast, we propose the hypothesis that rearrangements in gene order are more likely the result of selection during niche adaptation (SNAP). Partial chromosomal duplications occur very frequently by recombination between direct repeat sequences. Duplicated regions may contain tens to hundreds of genes and segregate quickly unless maintained by selection. Bacteria exposed to non-lethal selections (for example, a requirement to grow on a poor nutrient) can adapt by maintaining a duplication that includes a gene that improves relative fitness. Further improvements in fitness result from the loss or inactivation of non-selected genes within each copy of the duplication. When genes that are essential in single copy are lost from different copies of the duplication, segregation is prevented even if the original selection is lifted. Functional gene loss continues until a new genetic equilibrium is reached. The outcome is a rearranged gene order. Mathematical modelling shows that this process of positive selection to adapt to a new niche can rapidly drive rearrangements in gene order to fixation. Signature features (duplication formation and divergence) of the SNAP model were identified in natural isolates from multiple species showing that the initial two steps in the SNAP process can occur with a remarkably high frequency. Further bioinformatic and experimental analyses are required to test if and to which extend the SNAP process acts on bacterial genomes.

All life on earth has evolved from a universal common ancestor with a specific order of genes on the chromosome. This order is not maintained in modern species and the standard hypothesis is that changes reflect a lack of strong selection on gene order. Here, we propose an alternative hypothesis, SNAP. The occupation of a novel environment by bacteria is generally a trade-off situation. For example, while the bacteria may not be adapted to grow well under the new conditions, they may benefit by not having to share available resources with other microorganisms. Bacterial populations frequently acquire duplications of chromosomal segments containing genes that can help them adapt to a new environment. Other genes that are also duplicated are not required in two copies so that over time a superfluous copy can be lost. Eventually, the process of duplication and gene loss can lead to the rearrangement of the gene order in the chromosomal segment. The major benefit of this model over the standard hypothesis is that the process is driven by positive selection and can reach fixation rapidly.

Comparative Genomics Screen Identifies Microbe-Associated Molecular Patterns from 'Candidatus Liberibacter' spp. That Elicit Immune Responses in Plants

Citrus huanglongbing (HLB), caused by phloem-limited 'Candidatus Liberibacter' bacteria, is a destructive disease threatening the worldwide citrus industry. The mechanisms of pathogenesis are poorly understood and no efficient strategy is available to control HLB. Here, we used a comparative genomics screen to identify candidate microbe-associated molecular patterns (MAMPs) from 'Ca. Liberibacter' spp. We identified the core genome from multiple 'Ca. Liberibacter' pathogens, and searched for core genes with signatures of positive selection. We hypothesized that genes encoding putative MAMPs would evolve to reduce recognition by the plant immune system, while retaining their essential functions. To efficiently screen candidate MAMP peptides, we established a high-throughput microtiter plate-based screening assay, particularly for citrus, that measured reactive oxygen species (ROS) production, which is a common immune response in plants. We found that two peptides could elicit ROS production in Arabidopsis and Nicotiana benthamiana. One of these peptides elicited ROS production and defense gene expression in HLB-tolerant citrus genotypes, and induced MAMP-triggered immunity against the bacterial pathogen Pseudomonas syringae. Our findings identify MAMPs that boost immunity in citrus and could help prevent or reduce HLB infection.

Neutral evolution of cellular phenotypes

Eukaryotes exhibit a great diversity of cellular and subcellular morphologies, but their basic underlying architecture is fairly constant. All have a nucleus, Golgi, cytoskeleton, plasma membrane, vesicles, ribosomes, and all known lineages but one have mitochondrion-related organelles. Moreover, most eukaryotes undergo processes such as mitosis, meiosis, DNA recombination, and often perform feats such as phagocytosis, and amoeboid and flagellar movement. With all of these commonalities, it is obvious that eukaryotes evolved from a common ancestor, but it is not obvious how eukaryotes came to have their diverse structural phenotypes. Are these phenotypes adaptations to particular niches, their evolution dominated by positive natural selection? Or is eukaryotic cellular diversity substantially the product of neutral evolutionary processes, with adaptation either illusory or a secondary consequence? In this paper, we outline how a hierarchical view of phenotype can be used to articulate a neutral theory of phenotypic evolution, involving processes such as gene loss, gene replacement by homologues or analogues, gene duplication followed by subfunctionalization, and constructive neutral evolution. We suggest that neutral iterations of these processes followed by entrenchment of their products can explain much of the diversity of cellular, developmental, and biochemical phenotypes of unicellular eukaryotes and should be explored in addition to adaptive explanations.

Just the Two of Us? A Family of Pseudomonas Megaplasmids Offers a Rare Glimpse into the Evolution of Large Mobile Elements

Pseudomonads are ubiquitous group of environmental proteobacteria, well known for their roles in biogeochemical cycling, in the breakdown of xenobiotic materials, as plant growth promoters, and as pathogens of a variety of host organisms. We have previously identified a large megaplasmid present within one isolate of the plant pathogen Pseudomonas syringae, and here we report that a second member of this megaplasmid family is found within an environmental Pseudomonad isolate most closely related to Pseudomonas putida. Many of the shared genes are involved in critical cellular processes like replication, transcription, translation, and DNA repair. We argue that presence of these shared pathways sheds new light on discussions about the types of genes that undergo horizontal gene transfer (i.e., the complexity hypothesis) as well as the evolution of pangenomes. Furthermore, although both megaplasmids display a high level of synteny, genes that are shared differ by over 50% on average at the amino acid level. This combination of conservation in gene order despite divergence in gene sequence suggests that this Pseudomonad megaplasmid family is relatively old, that gene order is under strong selection within this family, and that there are likely many more members of this megaplasmid family waiting to be found in nature.

The Molecular Mechanisms Underlying Hidden Phenotypic Variation among Metallo-β-Lactamases

Genetic variation among orthologous genes has been largely formed through neutral genetic drift while maintaining the functional role of these genes. However, because the evolution of gene occurs in the context of each host organism, their sequence changes are also associated with adaptation to a specific environment. Thus, genetic variation can create critical phenotypic variation, particularly when genes are transferred to a new host by horizontal gene transfer. Unveiling "hidden phenotypic variation" is particularly important for genes that confer resistance to antibiotics. However, our understanding of the molecular mechanisms that underlie phenotypic variation remains limited. Here we sought to determine the extent of phenotypic variation in the B1 metallo-β-lactamase (MBL) family and its molecular basis by systematically characterizing eight MBL orthologs, including NDM-1 and VIM-2 and IMP-1. We found that these MBLs confer diverse levels of resistance. The phenotypic variation cannot be explained by variation in catalytic efficiency alone; rather, it is the combination of the catalytic efficiency and abundance of functional periplasmic enzyme that best predicts the observed variation in resistance. The level of functional periplasmic expression varied dramatically between MBL orthologs. This was the result of changes at multiple levels of each ortholog's: (1) quantity of mRNA, (2) amount of MBL expressed, and (3) efficacy of functional enzyme translocation to the periplasm. Overall, it is the interaction between each gene and the host's underlying cellular processes (transcription, translation, and translocation) that determines MBL genetic incompatibility through horizontal gene transfer. These host-specific processes may constrain the effective spread and deployment of MBLs to certain host species and could explain the current observed distribution bias.

Network hubs affect evolvability

The regulatory processes in cells are typically organized into complex genetic networks. However, it is still unclear how this network structure modulates the evolution of cellular regulation. One would expect that mutations in central and highly connected modules of a network (so-called hubs) would often result in a breakdown and therefore be an evolutionary dead end. However, a new study by Koubkova-Yu and colleagues finds that in some circumstances, altering a hub can offer a quick evolutionary advantage. Specifically, changes in a hub can induce significant phenotypic changes that allow organisms to move away from a local fitness peak, whereas the fitness defects caused by the perturbed hub can be mitigated by mutations in its interaction partners. Together, the results demonstrate how network architecture shapes and facilitates evolutionary adaptation.

Genes are organized into complex interaction networks, but it is unclear how network architecture affects evolution. This Primer explores a new study which uses experimental evolution to show how alterations in a gene central to a network affect evolutionary processes.

Heterologous Hsp90 promotes phenotypic diversity through network evolution

Biological processes in living cells are often carried out by gene networks in which signals and reactions are integrated through network hubs. Despite their functional importance, it remains unclear to what extent network hubs are evolvable and how alterations impact long-term evolution. We investigated these issues using heat shock protein 90 (Hsp90), a central hub of proteostasis networks. When native Hsp90 in Saccharomyces cerevisiae cells was replaced by the ortholog from hypersaline-tolerant Yarrowia lipolytica that diverged from S. cerevisiae about 270 million years ago, the cells exhibited improved growth in hypersaline environments but compromised growth in others, indicating functional divergence in Hsp90 between the two yeasts. Laboratory evolution shows that evolved Y. lipolytica-HSP90–carrying S. cerevisiae cells exhibit a wider range of phenotypic variation than cells carrying native Hsp90. Identified beneficial mutations are involved in multiple pathways and are often pleiotropic. Our results show that cells adapt to a heterologous Hsp90 by modifying different subnetworks, facilitating the evolution of phenotypic diversity inaccessible to wild-type cells.

Biological processes in living cells are often carried out by gene networks. Hubs are highly connected network components important for integrating signal inputs and generating responsive functional outputs. Heat shock protein 90 (Hsp90), a versatile hub in the protein homeostasis network, is a molecular chaperone essential for cell viability in all tested eukaryotic cells. In yeast, about a quarter of the expressed proteins are profoundly influenced when Hsp90 activity is reduced. Despite its pivotal role, we found that the function of Hsp90 has diverged between two yeast species, Yarrowia lipolytica and Saccharomyces cerevisiae, which split about 270 million years ago. To understand the impacts and adaptive strategies in cells with an altered network hub, we conducted laboratory evolution experiments using a S. cerevisiae strain in which native Hsp90 is replaced by its counterpart in Y. lipolytica. We observed different fitness gain or loss under various stress conditions in individual evolved clones, suggesting that cells adapted via different evolutionary paths. Genome sequencing and mutation reconstitution experiments show that beneficial mutations occurred in multiple Hsp90-related pathways that interact with each other. Our results show that a perturbed network allows cells to evolve a broader range of phenotypic diversity unavailable to wild-type cells.

Pooled CRISPR interference screening enables genome-scale functional genomics study in bacteria with superior performance

To fully exploit the microbial genome resources, a high-throughput experimental platform is needed to associate genes with phenotypes at the genome level. We present here a novel method that enables investigation of the cellular consequences of repressing individual transcripts based on the CRISPR interference (CRISPRi) pooled screening in bacteria. We identify rules for guide RNA library design to handle the unique structure of prokaryotic genomes by tiling screening and construct an E. coli genome-scale guide RNA library (~60,000 members) accordingly. We show that CRISPRi outperforms transposon sequencing, the benchmark method in the microbial functional genomics field, when similar library sizes are used or gene length is short. This tool is also effective for mapping phenotypes to non-coding RNAs (ncRNAs), as elucidated by a comprehensive tRNA-fitness map constructed here. Our results establish CRISPRi pooled screening as a powerful tool for mapping complex prokaryotic genetic networks in a precise and high-throughput manner.

Phylogenomic analysis demonstrates a pattern of rare and long-lasting concerted evolution in prokaryotes

Concerted evolution, where paralogs in the same species show higher sequence similarity to each other than to orthologs in other species, is widely found in many species. However, cases of concerted evolution that last for hundreds of millions of years are very rare. By genome-wide analysis of a broad selection of prokaryotes, we provide strong evidence of recurrent concerted evolution in 26 genes, most of which have lasted more than ~500 million years. We find that most concertedly evolving genes are key members of important pathways, and encode proteins from the same complexes and/or pathways, suggesting coevolution of genes via concerted evolution to maintain gene balance. We also present LRCE-DB, a comprehensive online repository of long-lasting concerted evolution. Collectively, our study reveals that although most duplicated genes may diverge in sequence over a long period, on rare occasions this constraint can be breached, leading to unexpected long-lasting concerted evolution in a recurrent manner.

Sishuo Wang and Youhua Chen present an analysis of concerted evolution in prokaryotes using a new computational pipeline, iSeeCE. They find evidence in 26 genes for recurrent concerted evolution, most of which last more than ~500 million years, and provide a database, LRCE-DB, for data exploration.

Biochemical mechanisms determine the functional compatibility of heterologous genes

Elucidating the factors governing the functional compatibility of horizontally transferred genes is important to understand bacterial evolution, including the emergence and spread of antibiotic resistance, and to successfully engineer biological systems. In silico efforts and work using single-gene libraries have suggested that sequence composition is a strong barrier for the successful integration of heterologous genes. Here we sample 200 diverse genes, representing >80% of sequenced antibiotic resistance genes, to interrogate the factors governing genetic compatibility in new hosts. In contrast to previous work, we find that GC content, codon usage, and mRNA-folding energy are of minor importance for the compatibility of mechanistically diverse gene products at moderate expression. Instead, we identify the phylogenetic origin, and the dependence of a resistance mechanism on host physiology, as major factors governing the functionality and fitness of antibiotic resistance genes. These findings emphasize the importance of biochemical mechanism for heterologous gene compatibility, and suggest physiological constraints as a pivotal feature orienting the evolution of antibiotic resistance.

Sequence composition is thought to be a major factor governing the functionality of horizontally transferred genes. In contrast, Porse et al. show that phylogenetic origin, and the type of resistance mechanism, are major factors affecting the functionality of horizontally transferred antibiotic resistance genes.

Resurrecting ancestral genes in bacteria to interpret ancient biosignatures

Two datasets, the geologic record and the genetic content of extant organisms, provide complementary insights into the history of how key molecular components have shaped or driven global environmental and macroevolutionary trends. Changes in global physico-chemical modes over time are thought to be a consistent feature of this relationship between Earth and life, as life is thought to have been optimizing protein functions for the entirety of its approximately 3.8 billion years of history on the Earth. Organismal survival depends on how well critical genetic and metabolic components can adapt to their environments, reflecting an ability to optimize efficiently to changing conditions. The geologic record provides an array of biologically independent indicators of macroscale atmospheric and oceanic composition, but provides little in the way of the exact behaviour of the molecular components that influenced the compositions of these reservoirs. By reconstructing sequences of proteins that might have been present in ancient organisms, we can downselect to a subset of possible sequences that may have been optimized to these ancient environmental conditions. How can one use modern life to reconstruct ancestral behaviours? Configurations of ancient sequences can be inferred from the diversity of extant sequences, and then resurrected in the laboratory to ascertain their biochemical attributes. One way to augment sequence-based, single-gene methods to obtain a richer and more reliable picture of the deep past, is to resurrect inferred ancestral protein sequences in living organisms, where their phenotypes can be exposed in a complex molecular-systems context, and then to link consequences of those phenotypes to biosignatures that were preserved in the independent historical repository of the geological record. As a first step beyond single-molecule reconstruction to the study of functional molecular systems, we present here the ancestral sequence reconstruction of the beta-carbonic anhydrase protein. We assess how carbonic anhydrase proteins meet our selection criteria for reconstructing ancient biosignatures in the laboratory, which we term palaeophenotype reconstruction.This article is part of the themed issue 'Reconceptualizing the origins of life'.

Constraining the timing of the Great Oxidation Event within the Rubisco phylogenetic tree

Ribulose 1,5-bisphosphate (RuBP) carboxylase/oxygenase (RuBisCO, or Rubisco) catalyzes a key reaction by which inorganic carbon is converted into organic carbon in the metabolism of many aerobic and anaerobic organisms. Across the broader Rubisco protein family, homologs exhibit diverse biochemical characteristics and metabolic functions, but the evolutionary origins of this diversity are unclear. Evidence of the timing of Rubisco family emergence and diversification of its different forms has been obscured by a meager paleontological record of early Earth biota, their subcellular physiology and metabolic components. Here, we use computational models to reconstruct a Rubisco family phylogenetic tree, ancestral amino acid sequences at branching points on the tree, and protein structures for several key ancestors. Analysis of historic substitutions with respect to their structural locations shows that there were distinct periods of amino acid substitution enrichment above background levels near and within its oxygen-sensitive active site and subunit interfaces over the divergence between Form III (associated with anoxia) and Form I (associated with oxia) groups in its evolutionary history. One possible interpretation is that these periods of substitutional enrichment are coincident with oxidative stress exerted by the rise of oxygenic photosynthesis in the Precambrian era. Our interpretation implies that the periods of Rubisco substitutional enrichment inferred near the transition from anaerobic Form III to aerobic Form I ancestral sequences predate the acquisition of Rubisco by fully derived cyanobacterial (i.e., dual photosystem-bearing, oxygen-evolving) clades. The partitioning of extant lineages at high clade levels within our Rubisco phylogeny indicates that horizontal transfer of Rubisco is a relatively infrequent event. Therefore, it is possible that the mutational enrichment periods between the Form III and Form I common ancestral sequences correspond to the adaptation of key oxygen-sensitive components of Rubisco prior to, or coincident with, the Great Oxidation Event.