Prokaryotic Evolution Research Articles

Long range segmentation of prokaryotic genomes by gene age and functionality.

Bacterial and archaeal genomes encompass numerous operons that typically consist of two to five genes. On larger scales, however, gene order is poorly conserved through the evolution of prokaryotes. Nevertheless, non-random localization of different classes of genes on prokaryotic chromosomes could reflect important functional and evolutionary constraints. We explored the patterns of genomic localization of evolutionarily conserved (ancient) and variable (young) genes across the diversity of bacteria and archaea. Nearly all bacterial and archaeal chromosomes were found to encompass large segments of 100-300 kb that were significantly enriched in either ancient or young genes. Similar clustering of genes with lethal knockout phenotype (essential genes) was observed as well. Mathematical modeling of genome evolution suggests that this long-range gene clustering in prokaryotic chromosomes reflects perpetual genome rearrangement driven by a combination of selective and neutral processes rather than evolutionary conservation.

Options and considerations for validation of prokaryotic names under the SeqCode

Stable taxon names for Bacteria and Archaea are essential for capturing and documenting prokaryotic diversity. They are also crucial for scientific communication, effective accumulation of biological data related to the taxon names and for developing a comprehensive understanding of prokaryotic evolution. However, after more than a hundred years, taxonomists have succeeded in valid publication of only around 30 000 species names, based mostly on pure cultures under the International Code of Nomenclature of Prokaryotes (ICNP), out of the millions estimated to reside in the biosphere. The vast majority of prokaryotic species have not been cultured and are becoming increasingly known to us via culture-independent sequence-based approaches. Until recently, such taxa could only be addressed nomenclaturally via provisional names such as Candidatus or alphanumeric identifiers. Here, we present options and considerations to facilitate validation of names for these taxa using the recently established Code of Nomenclature of Prokaryotes Described from Sequence Data (SeqCode). Community engagement and participation of relevant taxon specialists are critical and encouraged for the success of endeavours to formally name the uncultured majority.

We are doing it wrong: Putting homology before phylogeny in cyanobacterial taxonomy.

The rapid expansion of whole genome sequencing in bacterial taxonomy has revealed deep evolutionary relationships and speciation signals, but assembly methods often miss true nucleotide diversity in the ribosomal operons. Though it lacks sufficient phylogenetic signal at the species level, the 16S ribosomal RNA gene is still much used in bacterial taxonomy. In cyanobacterial taxonomy, comparisons of 16S-23S Internal Transcribed Spacer (ITS) regions are used to bridge this information gap. Although ITS rRNA region analyses are routinely being used to identify species, researchers often do not identify orthologous operons, which leads to improper comparisons. No method for delineating orthologous operon copies from paralogous ones has been established. A new method for recognizing orthologous ribosomal operons by quantifying the conserved paired nucleotides in a helical domain of the ITS, has been developed. The D1' Index quantifies differences in the ratio of pyrimidines to purines in paired nucleotide sequences of this helix. Comparing 111 operon sequences from 89 strains of Brasilonema, four orthologous operon types were identified. Plotting D1' Index values against the length of helices produced clear separation of orthologs. Most orthologous operons in this study were observed both with and without tRNA genes present. We hypothesize that genomic rearrangement, not gene duplication, is responsible for the variation among orthologs. This new method will allow cyanobacterial taxonomists to utilize ITS rRNA region data more correctly, preventing erroneous taxonomic hypotheses. Moreover, this work could assist genomicists in identifying and preserving evident sequence variability in ribosomal operons, which is an important proxy for evolution in prokaryotes.

Rapid functional activation of horizontally transferred eukaryotic intron-containing genes in the bacterial recipient.

Horizontal gene transfer has occurred across all domains of life and contributed substantially to the evolution of both prokaryotes and eukaryotes. Previous studies suggest that many horizontally transferred eukaryotic genes conferred selective advantages to bacterial recipients, but how these eukaryotic genes evolved into functional bacterial genes remained unclear, particularly how bacteria overcome the expressional barrier posed by eukaryotic introns. Here, we first confirmed that the presence of intron would inactivate the horizontally transferred gene in Escherichiacoli even if this gene could be efficiently transcribed. Subsequent large-scale genetic screens for activation of gene function revealed that activation events could rapidly occur within several days of selective cultivation. Molecular analysis of activation events uncovered two distinct mechanisms how bacteria overcome the intron barrier: (i) intron was partially deleted and the resulting stop codon-removed mutation led to one intact foreign protein or (ii) intron was intactly retained but it mediated the translation initiation and the interaction of two split small proteins (derived from coding sequences up- and downstream of intron, respectively) to restore gene function. Our findings underscore the likelihood that horizontally transferred eukaryotic intron-containing genes could rapidly acquire functionality if they confer a selective advantage to the prokaryotic recipient.

Library-Assisted Evolution in Eukaryotic Cells Yield Adenine Base Editors with Enhanced Editing Specificity.

The current-generation adenine base editor (ABE) ABE8e, which has evolved from the prokaryotic evolution system, exhibits high efficiency in mediating A-to-Gconversion and is presumed to be promising for gene therapy. However, its much wider editing window and substantially higher off-target editing activity restricted its applications in precise base editing for therapeutic use. This study uses a library-assistedprotein evolution approach using eukaryotic cells to generate ABE variants with improved specificity and reduced off-target editing while maintaining high activity in human cells. The study generated an expanded set of ABEs with efficient editing activities and chose fourevolved variants thatoffered either similar or modestly higher efficiency within a narrower editing window of protospacer position ≈4-7 compared to that of ABE8e in human cells, which would enable minimized bystander editing. Moreover, these variants resulted in reduced off-target editing events when delivered as plasmid or mRNA into human cells. Finally, these variants can install both disease-suppressing mutations and disease-correcting mutations efficiently with minimal undesired bystander editing making them promising approaches for specific therapeutic edits. In summary, the work establishes a mutant-library-assisted protein evolution method in eukaryotic cells and generates alternative ABE variants as efficient tools for precise human genome editing.

Long range segmentation of prokaryotic genomes by gene age and functionality.

Bacterial and archaeal genomes encompass numerous operons that typically consist of two to five genes. On larger scales, however, gene order is poorly conserved through the evolution of prokaryotes. Nevertheless, non-random localization of different classes of genes on prokaryotic chromosomes could reflect important functional and evolutionary constraints. We explored the patterns of genomic localization of evolutionarily conserved (ancient) and variable (young) genes across the diversity of bacteria and archaea. Nearly all bacterial and archaeal chromosomes were found to encompass large segments of 100-300 kilobases that were significantly enriched in either ancient or young genes. Similar clustering of genes with lethal knockout phenotype (essential genes) was observed as well. Mathematical modeling of genome evolution suggests that this long-range gene clustering in prokaryotic chromosomes reflects perpetual genome rearrangement driven by a combination of selective and neutral processes rather than evolutionary conservation.

Moving toward the Inclusion of Epigenomics in Bacterial Genome Evolution: Perspectives and Challenges.

The universality of DNA methylation as an epigenetic regulatory mechanism belongs to all biological kingdoms. However, while eukaryotic systems have been the primary focus of DNA methylation studies, the molecular mechanisms in prokaryotes are less known. Nevertheless, DNA methylation in prokaryotes plays a pivotal role in many cellular processes such as defense systems against exogenous DNA, cell cycle dynamics, and gene expression, including virulence. Thanks to single-molecule DNA sequencing technologies, genome-wide identification of methylated DNA is becoming feasible on a large scale, providing the possibility to investigate more deeply the presence, variability, and roles of DNA methylation. Here, we present an overview of the multifaceted roles of DNA methylation in prokaryotes and suggest research directions and tools which can enable us to better understand the contribution of DNA methylation to prokaryotic genome evolution and adaptation. In particular, we emphasize the need to understand the presence and role of transgenerational inheritance, as well as the impact of epigenomic signatures on adaptation and genome evolution. Research directions and the importance of novel computational tools are underlined.

Defence systems and horizontal gene transfer in bacteria.

Horizontal gene transfer (HGT) is a fundamental process in prokaryotic evolution, contributing significantly to diversification and adaptation. HGT is typically facilitated by mobile genetic elements (MGEs), such as conjugative plasmids and phages, which often impose fitness costs on their hosts. However, a considerable number of bacterial genes are involved in defence mechanisms that limit the propagation of MGEs, suggesting they may actively restrict HGT. In our study, we investigated whether defence systems limit HGT by examining the relationship between the HGT rate and the presence of 73 defence systems across 12 bacterial species. We discovered that only six defence systems, three of which were different CRISPR-Cas subtypes, were associated with a reduced gene gain rate at the species evolution scale. Hosts of these defence systems tend to have a smaller pangenome size and fewer phage-related genes compared to genomes without these systems. This suggests that these defence mechanisms inhibit HGT by limiting prophage integration. We hypothesize that the restriction of HGT by defence systems is species-specific and depends on various ecological and genetic factors, including the burden of MGEs and the fitness effect of HGT in bacterial populations.

A global survey of prokaryotic genomes reveals the eco-evolutionary pressures driving horizontal gene transfer

Horizontal gene transfer, the exchange of genetic material through means other than reproduction, is a fundamental force in prokaryotic genome evolution. Genomic persistence of horizontally transferred genes has been shown to be influenced by both ecological and evolutionary factors. However, there is limited availability of ecological information about species other than the habitats from which they were isolated, which has prevented a deeper exploration of ecological contributions to horizontal gene transfer. Here we focus on transfers detected through comparison of individual gene trees to the species tree, assessing the distribution of gene-exchanging prokaryotes across over a million environmental sequencing samples. By analysing detected horizontal gene transfer events, we show distinct functional profiles for recent versus old events. Although most genes transferred are part of the accessory genome, genes transferred earlier in evolution tend to be more ubiquitous within present-day species. We find that co-occurring, interacting and high-abundance species tend to exchange more genes. Finally, we show that host-associated specialist species are most likely to exchange genes with other host-associated specialist species, whereas species found across different habitats have similar gene exchange rates irrespective of their preferred habitat. Our study covers an unprecedented scale of integrated horizontal gene transfer and environmental information, highlighting broad eco-evolutionary trends.

Contingency, repeatability, and predictability in the evolution of a prokaryotic pangenome

Pangenomes exhibit remarkable variability in many prokaryotic species, much of which is maintained through the processes of horizontal gene transfer and gene loss. Repeated acquisitions of near-identical homologs can easily be observed across pangenomes, leading to the question of whether these parallel events potentiate similar evolutionary trajectories, or whether the remarkably different genetic backgrounds of the recipients mean that postacquisition evolutionary trajectories end up being quite different. In this study, we present a machine learning method that predicts the presence or absence of genes in the Escherichia coli pangenome based on complex patterns of the presence or absence of other accessory genes within a genome. Our analysis leverages the repeated transfer of genes through the E. coli pangenome to observe patterns of repeated evolution following similar events. We find that the presence or absence of a substantial set of genes is highly predictable from other genes alone, indicating that selection potentiates and maintains gene-gene co-occurrence and avoidance relationships deterministically over long-term bacterial evolution and is robust to differences in host evolutionary history. We propose that at least part of the pangenome can be understood as a set of genes with relationships that govern their likely cohabitants, analogous to an ecosystem's set of interacting organisms. Our findings indicate that intragenomic gene fitness effects may be key drivers of prokaryotic evolution, influencing the repeated emergence of complex gene-gene relationships across the pangenome.

Identification of recG genetic interactions in Escherichia coli by transposon sequencing.

DNA damage and subsequent DNA repair processes are mutagenic in nature and an important driver of evolution in prokaryotes, including antibiotic resistance development. Genetic screening approaches, such as transposon sequencing (Tn-seq), have provided important new insights into gene function and genetic relationships. Here, we employed Tn-seq to gain insight into the function of the recG gene, which renders Escherichia coli cells moderately sensitive to a variety of DNA-damaging agents when they are absent. The reported recG genetic interactions can be used in combination with future screens to aid in a more complete reconstruction of DNA repair pathways in bacteria.

Creating Polyploid Escherichia Coli and Its Application in Efficient L-Threonine Production.

Prokaryotic genomes are generally organized in haploid. In synthetic biological research, efficient chassis cells must be constructed to produce bio-based products. Here, the essential division of the ftsZ gene to create functional polyploid E. coli is regulated. The artificial polyploid E. coli containing 2-4 chromosomes is confirmed through PCR amplification, terminator localization, and flow cytometry. The polyploid E. coli exhibitsa larger cell size, and its low pH tolerance and acetate resistance are stronger than those of haploid E. coli. Transcriptome analysis showsthat the genes of the cell's main functional pathways are significantly upregulated in the polyploid E. coli. These advantages of the polyploid E. coli results in the highest reported L-threonine yield (160.3gL-1 ) in fed-batch fermentation to date. In summary, an easy and convenient method for constructing polyploid E. coli and demonstrated its application in L-threonine production is developed. This work provides a new approach for creating an excellent host strain for biochemical production and studying the evolution of prokaryotes and their chromosome functions.

Compensatory relationship between low-complexity regions and gene paralogy in the evolution of prokaryotes

The evolution of genomes in all life forms involves two distinct, dynamic types of genomic changes: gene duplication (and loss) that shape families of paralogous genes and extension (and contraction) of low-complexity regions (LCR), which occurs through dynamics of short repeats in protein-coding genes. Although the roles of each of these types of events in genome evolution have been studied, their co-evolutionary dynamics is not thoroughly understood. Here, by analyzing a wide range of genomes from diverse bacteria and archaea, we show that LCR and paralogy represent two distinct routes of evolution that are inversely correlated. The emergence of LCR is a prominent evolutionary mechanism in fast evolving, young protein families, whereas paralogy dominates the comparatively slow evolution of old protein families. The analysis of multiple prokaryotic genomes shows that the formation of LCR is likely a widespread, transient evolutionary mechanism that temporally and locally affects also ancestral functions, but apparently, fades away with time, under mutational and selective pressures, yielding to gene paralogy. We propose that compensatory relationships between short-term and longer-term evolutionary mechanisms are universal in the evolution of life.

Rhometa: Population recombination rate estimation from metagenomic read datasets.

Prokaryotic evolution is influenced by the exchange of genetic information between species through a process referred to as recombination. The rate of recombination is a useful measure for the adaptive capacity of a prokaryotic population. We introduce Rhometa (https://github.com/sid-krish/Rhometa), a new software package to determine recombination rates from shotgun sequencing reads of metagenomes. It extends the composite likelihood approach for population recombination rate estimation and enables the analysis of modern short-read datasets. We evaluated Rhometa over a broad range of sequencing depths and complexities, using simulated and real experimental short-read data aligned to external reference genomes. Rhometa offers a comprehensive solution for determining population recombination rates from contemporary metagenomic read datasets. Rhometa extends the capabilities of conventional sequence-based composite likelihood population recombination rate estimators to include modern aligned metagenomic read datasets with diverse sequencing depths, thereby enabling the effective application of these techniques and their high accuracy rates to the field of metagenomics. Using simulated datasets, we show that our method performs well, with its accuracy improving with increasing numbers of genomes. Rhometa was validated on a real S. pneumoniae transformation experiment, where we show that it obtains plausible estimates of the rate of recombination. Finally, the program was also run on ocean surface water metagenomic datasets, through which we demonstrate that the program works on uncultured metagenomic datasets.

Research progress on horizontal gene transfer and its functions in insects

Horizontal gene transfer (HGT) refers to the asexual exchange of genetic information between distantly related organisms. Although it is well acknowledged that HGT greatly contributes to the adaptive evolution in prokaryotes, its significance in shaping the orientation of eukaryote evolution remains obscure, especially in insect lineages. The massively expanded genomic data appears to be an excellent choice to uncover the mystery of HGTs in insects nowadays. Here we gather a body of evidence showing the HGT events from three broad donor origins, viruses, prokaryotes, and eukaryotes, in which most horizontally transferred genes are unlikely to be functional and will be eroded as a result of the difference of inheritable background between insect recipients and donor species. Nevertheless, particular interests in the prominent role of insect HGTs in maintaining and acquiring new functionalities have still been raised to underpin their adaptations. Among these, the previously investigated properties including reproduction regulation, detoxification of plant metabolites, formation of body color, and antimicrobic immunity are mainly included in different insect lineages. Albeit such cases are just the tip of the iceberg, we indeed demonstrate that HGT drives insect evolution, especially in coevolution with host plants, and additional explorations into its functions should be paid special attention to access the complex evolutionary history of insects in the near future.

How has horizontal gene transfer shaped the evolution of insect genomes?

As the most diverse group of animals on Earth, insects are key organisms in ecosystems. Horizontal gene transfer (HGT) refers to the transfer of genetic material between species by non-reproductive means. HGT is a major evolutionary force in prokaryotic genome evolution, but its importance in different eukaryotic groups, such as insects, has only recently begun to be understood. Genomic data from hundreds of insect species have enabled the detection of large numbers of HGT events and the elucidation of the functions of some of these foreign genes. Although quantification of the extent of HGT in insects broadens our understanding of its role in insect evolution, the scope of its influence and underlying mechanism(s) of its occurrence remain open questions for the field.

The squalene route to C30 carotenoid biosynthesis and the origins of carotenoid biosynthetic pathways

Carotenoids are isoprenoid lipids found across the tree of life with important implications in oxidative stress adaptations, photosynthetic metabolisms, as well as in membrane dynamics. The canonical view is that C40 carotenoids are synthesized from phytoene and C30 carotenoids from diapophytoene. Squalene is mostly associated with the biosynthesis of polycyclic triterpenes, although there have been suggestions that it could also be involved in the biosynthesis of C30 carotenoids. However, demonstration of the existence of this pathway in nature is lacking. Here, we demonstrate that C30 carotenoids are synthesized from squalene in the Planctomycetes bacteria and that this squalene route to C30 carotenoids is the most widespread in prokaryotes. Using the evolutionary history of carotenoid and squalene amino oxidases, we propose an evolutionary scenario to explain the origin and diversification of the different carotenoid and squalene-related pathways. Weshow that carotenoid biosynthetic pathways have been constantly transferred and neofunctionalized during prokaryotic evolution. One possible origin of the squalene pathway connects it with the one of C40 carotenoid synthesis of Cyanobacteria. The widespread occurrence of the squalene route to C30 carotenoids in Bacteria increases the functional repertoire of squalene, establishing it as a general hub of carotenoids and polycyclic triterpenes synthesis.

Discovery of integrons in Archaea: Platforms for cross-domain gene transfer.

Horizontal gene transfer between different domains of life is increasingly being recognized as an important evolutionary driver, with the potential to increase the pace of biochemical innovation and environmental adaptation. However, the mechanisms underlying the recruitment of exogenous genes from foreign domains are mostly unknown. Integrons are a family of genetic elements that facilitate this process within Bacteria. However, they have not been reported outside Bacteria, and thus their potential role in cross-domain gene transfer has not been investigated. Here, we discover that integrons are also present in 75 archaeal metagenome-assembled genomes from nine phyla, and are particularly enriched among Asgard archaea. Furthermore, we provide experimental evidence that integrons can facilitate the recruitment of archaeal genes by bacteria. Our findings establish a previously unknown mechanism of cross-domain gene transfer whereby bacteria can incorporate archaeal genes from their surrounding environment via integron activity. These findings have important implications for prokaryotic ecology and evolution.

Reconstructing horizontal gene flow network to understand prokaryotic evolution.

Horizontal gene transfer (HGT) is a major source of phenotypic innovation and a mechanism of niche adaptation in prokaryotes. Quantification of HGT is critical to decipher its myriad roles in microbial evolution and adaptation. Advances in genome sequencing and bioinformatics have augmented our ability to understand the microbial world, particularly the direct or indirect influence of HGT on diverse life forms. Methods for detecting HGT can be classified into phylogenetic-based and parametric or composition-based approaches. Here, we exploited the complementary strengths of both the approaches to construct a high confidence horizontal gene flow network. Our network is unique in its ability to detect the transfer of native genes of a genome to genomes from other taxa, thus establishing donor and recipient organisms (taxa), rather than through a post hoc analysis as is the practice with several other approaches. The scale-free horizontal gene flow network presented here provides new insights into modes of transfer for the exchange of genetic information and also illuminates differential gene flow across phyla.

Gene-influx-driven evolution.

Here we analyze the evolutionary process in the presence of continuous influx of genotypes with submaximum fitness from the outside to the given habitat with finite resources. We show that strong influx from the outside allows the low-fitness genotype to win the competition with the higher fitness genotype, and in a finite population, drive the latter to extinction. We analyze a mathematical model of this phenomenon and obtain the conditions for the transition from the high-fitness to the low-fitness genotype caused by the influx of the latter. We calculate the time to extinction of the high-fitness genotype in a finite population with two alleles and find the exact analytical dynamics of extinction for the case of many genes with epistasis. We solve a related quasispecies model for a single peak (random) fitness landscape as well as for a symmetric fitness landscape. In the symmetric landscape, a nonperturbative effect is observed such that even an extremely low influx of the low-fitness genotype drastically changes the steady state fitness distribution. A similar nonperturbative phenomenon is observed for the allele fixation time as well. The identified regime of influx-driven evolution appears to be relevant for a broad class of biological systems and could be central to the evolution of prokaryotes and viruses.