Improving Phylogenies Based on Average Nucleotide Identity, Incorporating Saturation Correction and Nonparametric Bootstrap Support.

Whole-genome comparisons based on average nucleotide identities (ANI) and the genome-to-genome distance calculator have risen to prominence in rapidly classifying prokaryotic taxa using whole-genome sequences. Some implementations have even been proposed as a new standard in species classification and have become a common technique for papers describing newly sequenced genomes. However, attempts to apply whole-genome divergence data to the delineation of higher taxonomic units and to phylogenetic inference have had difficulty matching those produced by more complex phylogenetic methods. We present a novel method for generating statistically supported phylogenies of archaeal and bacterial groups using a combined ANI and alignment fraction-based metric. For the test cases to which we applied the developed approach, we obtained results comparable with other methodologies up to at least the family level. The developed method uses nonparametric bootstrapping to gauge support for inferred groups. This method offers the opportunity to make use of whole-genome comparison data, that is already being generated, to quickly produce phylogenies including support for inferred groups. Additionally, the developed ANI methodology can assist the classification of higher taxonomic groups.[Average nucleotide identity (ANI); genome evolution; prokaryotic species delineation; taxonomy.].

Systematic Detection of Large-Scale Multigene Horizontal Transfer in Prokaryotes.

Horizontal gene transfer (HGT) is central to prokaryotic evolution. However, little is known about the "scale" of individual HGT events. In this work, we introduce the first computational framework to help answer the following fundamental question: How often does more than one gene get horizontally transferred in a single HGT event? Our method, called HoMer, uses phylogenetic reconciliation to infer single-gene HGT events across a given set of species/strains, employs several techniques to account for inference error and uncertainty, combines that information with gene order information from extant genomes, and uses statistical analysis to identify candidate horizontal multigene transfers (HMGTs) in both extant and ancestral species/strains. HoMer is highly scalable and can be easily used to infer HMGTs across hundreds of genomes. We apply HoMer to a genome-scale data set of over 22,000 gene families from 103 Aeromonas genomes and identify a large number of plausible HMGTs of various scales at both small and large phylogenetic distances. Analysis of these HMGTs reveals interesting relationships between gene function, phylogenetic distance, and frequency of multigene transfer. Among other insights, we find that 1) the observed relative frequency of HMGT increases as divergence between genomes increases, 2) HMGTs often have conserved gene functions, and 3) rare genes are frequently acquired through HMGT. We also analyze in detail HMGTs involving the zonula occludens toxin and type III secretion systems. By enabling the systematic inference of HMGTs on a large scale, HoMer will facilitate a more accurate and more complete understanding of HGT and microbial evolution.

Horizontal gene transfer: building the web of life.

Horizontal gene transfer (HGT) is the sharing of genetic material between organisms that are not in a parent-offspring relationship. HGT is a widely recognized mechanism for adaptation in bacteria and archaea. Microbial antibiotic resistance and pathogenicity are often associated with HGT, but the scope of HGT extends far beyond disease-causing organisms. In this Review, we describe how HGT has shaped the web of life using examples of HGT among prokaryotes, between prokaryotes and eukaryotes, and even between multicellular eukaryotes. We discuss replacement and additive HGT, the proposed mechanisms of HGT, selective forces that influence HGT, and the evolutionary impact of HGT on ancestral populations and existing populations such as the human microbiome.

Impact of a homing intein on recombination frequency and organismal fitness.

Inteins are parasitic genetic elements that excise themselves at the protein level by self-splicing, allowing the formation of functional, nondisrupted proteins. Many inteins contain a homing endonuclease (HEN) domain and rely on its activity for horizontal propagation. However, successful invasion of an entire population will make this activity redundant, and the HEN domain is expected to degenerate quickly under these conditions. Several theories have been proposed for the continued existence of the both active HEN and noninvaded alleles within a population. However, to date, these models were not directly tested experimentally. Using the natural cell fusion ability of the halophilic archaeon Haloferax volcanii we were able to examine this question in vivo, by mating polB intein-positive [insertion site c in the gene encoding DNA polymerase B (polB-c)] and intein-negative cells and examining the dispersal efficiency of this intein in a natural, polyploid population. Through competition between otherwise isogenic intein-positive and intein-negative strains we determined a surprisingly high fitness cost of over 7% for the polB-c intein. Our laboratory culture experiments and samples taken from Israel's Mediterranean coastline show that the polB-c inteins do not efficiently take over an inteinless population through mating, even under ideal conditions. The presence of the HEN/intein promoted recombination when intein-positive and intein-negative cells were mated. Increased recombination due to HEN activity contributes not only to intein dissemination but also to variation at the population level because recombination tracts during repair extend substantially from the homing site.

The impact of HGT on phylogenomic reconstruction methods.

Supermatrix and supertree analyses are frequently used to more accurately recover vertical evolutionary history but debate still exists over which method provides greater reliability. Traditional methods that resolve relationships among organisms from single genes are often unreliable because of the frequent lack of strong phylogenetic signal and the presence of systematic artifacts. Methods developed to reconstruct organismal history from multiple genes can be divided into supermatrix and supertree approaches. A supermatrix analysis consists of the concatenation of multiple genes into a single, possibly partitioned alignment, from which phylogenies are reconstructed using a variety of approaches. Supertrees build consensus trees from the topological information contained within individual gene trees. Both methods are now widely used and have been demonstrated to solve previously ambiguous or unresolved phylogenies with high statistical support. However, the amount of misleading signal needed to induce erroneous phylogenies for both strategies is still unknown. Using genome simulations, we test the accuracy of supertree and supermatrix approaches in recovering the true organismal phylogeny under increased amounts of horizontally transferred genes and changes in substitution rates. Our results show that overall, supermatrix approaches are preferable when a low amount of gene transfer is suspected to be present in the dataset, while supertrees have greater reliability in the presence of a moderate amount of misleading gene transfers. In the face of very high or very low substitution rates without horizontal gene transfers, supermatrix approaches outperform supertrees as individual gene trees remain unresolved and additional sequences contribute to a congruent phylogenetic signal.

Ancient horizontal gene transfer and the last common ancestors.

BACKGROUND: The genomic history of prokaryotic organismal lineages is marked by extensive horizontal gene transfer (HGT) between groups of organisms at all taxonomic levels. These HGT events have played an essential role in the origin and distribution of biological innovations. Analyses of ancient gene families show that HGT existed in the distant past, even at the time of the organismal last universal common ancestor (LUCA). Most gene transfers originated in lineages that have since gone extinct. Therefore, one cannot assume that the last common ancestors of each gene were all present in the same cell representing the cellular ancestor of all extant life. RESULTS: Organisms existing as part of a diverse ecosystem at the time of LUCA likely shared genetic material between lineages. If these other lineages persisted for some time, HGT with the descendants of LUCA could have continued into the bacterial and archaeal lineages. Phylogenetic analyses of aminoacyl-tRNA synthetase protein families support the hypothesis that the molecular common ancestors of the most ancient gene families did not all coincide in space and time. This is most apparent in the evolutionary histories of seryl-tRNA synthetase and threonyl-tRNA synthetase protein families, each containing highly divergent "rare" forms, as well as the sparse phylogenetic distributions of pyrrolysyl-tRNA synthetase, and the bacterial heterodimeric form of glycyl-tRNA synthetase. These topologies and phyletic distributions are consistent with horizontal transfers from ancient, likely extinct branches of the tree of life. CONCLUSIONS: Of all the organisms that may have existed at the time of LUCA, by definition only one lineage is survived by known progeny; however, this lineage retains a genomic record of heterogeneous genetic origins. The evolutionary histories of aminoacyl-tRNA synthetases (aaRS) are especially informative in detecting this signal, as they perform primordial biological functions, have undergone several ancient HGT events, and contain many sites with low substitution rates allowing deep phylogenetic reconstruction. We conclude that some aaRS families contain groups that diverge before LUCA. We propose that these ancient gene variants be described by the term "hypnologs", reflecting their ancient, reticulate origin from a time in life history that has been all but erased".

Reconstruction of ancestral 16S rRNA reveals mutation bias in the evolution of optimal growth temperature in the Thermotogae phylum.

Optimal growth temperature is a complex trait involving many cellular components, and its physiology is not yet fully understood. Evolution of continuous characters, such as optimal growth temperature, is often modeled as a one-dimensional random walk, but such a model may be an oversimplification given the complex processes underlying the evolution of continuous characters. Recent articles have used ancestral sequence reconstruction to infer the optimal growth temperature of ancient organisms from the guanine and cytosine content of the stem regions of ribosomal RNA, allowing inferences about the evolution of optimal growth temperature. Here, we investigate the optimal growth temperature of the bacterial phylum Thermotogae. Ancestral sequence reconstruction using a nonhomogeneous model was used to reconstruct the stem guanine and cytosine content of 16S rRNA sequences. We compare this sequence reconstruction method with other ancestral character reconstruction methods, and show that sequence reconstruction generates smaller confidence intervals and different ancestral values than other reconstruction methods. Unbiased random walk simulation indicates that the lower temperature members of the Thermotogales have been under directional selection; however, when a simulation is performed that takes possible mutations into account, it is the high temperature lineages that are, in fact, under directional selection. We find that the evolution of Thermotogales optimal growth temperatures is best fit by a biased random walk model. These findings suggest that it may be easier to evolve from a high optimal growth temperature to a lower one than vice versa.

Distribution and evolution of the mobile vma-1b intein.

Inteins are self-splicing parasitic genetic elements found in all domains of life. These genetic elements are found in highly conserved positions in conserved proteins. One protein family that has been invaded by inteins is the vacuolar and archaeal catalytic ATPase subunits (vma-1). There are two intein insertion sites in this protein, "a" and "b." The b site was previously thought to be only invaded in archaeal lineages. Here we survey the distribution and evolutionary histories of the b site inteins and show that the intein is present in more lineages than previously annotated, including a bacterial lineage, Mahella australiensis 50-1 BON. We present evidence, through ancestral character state reconstruction and substitution ratios between host genes and inteins, for several transfers of this intein between divergent species, including an interdomain transfer between the archaea and bacteria. Although inteins may persist within a single population or species for long periods of time, transfer of the vma-1b intein between divergent species contributed to the distribution of this intein.

Using the pan-genomic framework for the discovery of genomic islands in the haloarchaeon Halorubrum ezzemoulense.

In this study, we use pan-genomics to characterize the genomic variability of the widely dispersed halophilic archaeal species Halorubrum ezzemoulense (Hez). We include a multi-regional sampling of newly sequenced, high-quality draft genomes. The pan-genome graph of the species reveals 50 genomic islands that represent rare accessory genetic capabilities available to members. Most notably, we observe rearrangements that have led to the insertion/recombination/replacement of mutually exclusive genomic islands in equivalent genome positions ("homeocassettes"). These conflicting islands encode for similar functions, but homologs from islands located between the same core genes exhibit high divergence on the amino acid level, while the neighboring core genes are nearly identical. Both islands of a homeocassette often coexist in the same geographic location, suggesting that either island may be beyond the reach of selective sweeps and that these loci of divergence between Hez members are maintained and persist long term. This implies that subsections of the population have different niche preferences and rare metabolic capabilities. After an evaluation of the gene content in the homeocassettes, we speculate that these islands may play a role in the speciation, niche adaptability, and group selection dynamics in Hez. Though homeocassettes are first described in this study, similar replacements and divergence of genes on genomic islands have been previously reported in other Haloarchaea and distantly related Archaea, suggesting that homeocassettes may be a feature in a wide range of organisms outside of Hez.IMPORTANCEThis study catalogs the rare genes discovered in strains of the species Halorubrum ezzemoulense (Hez), an obligate halophilic archaeon, through the perspective of its pan-genome. These rare genes are often found to be arranged on islands that confer metabolic and transport functions and contain genes that have eluded previous studies. The discovery of divergent, but homologous islands occupying equivalent genome positions ("homeocassettes") in different genomes, reveals significant new information on genome evolution in Hez. Homeocassette pairs encode for similar functions, but their dissimilarity and distribution imply high rates of recombination, different specializations, and niche preferences in Hez. The coexistence of both islands of a homeocassette pair in multiple environments demonstrates that both islands are beyond the reach of selective sweeps and that these genome content differences between strains persist long term. The switch between islands through recombination under different environmental conditions may lead to a greater range of niche adaptability in Hez.

The effects of model choice and mitigating bias on the ribosomal tree of life.

Deep-level relationships within Bacteria, Archaea, and Eukarya as well as the relationships of these three domains to each other require resolution. The ribosomal machinery, universal to all cellular life, represents a protein repertoire resistant to horizontal gene transfer, which provides a largely congruent signal necessary for reconstructing a tree suitable as a backbone for life's reticulate history. Here, we generate a ribosomal tree of life from a robust taxonomic sampling of Bacteria, Archaea, and Eukarya to elucidate deep-level intra-domain and inter-domain relationships. Lack of phylogenetic information and systematic errors caused by inadequate models (that cannot account for substitution rate or compositional heterogeneities) or improper model selection compound conflicting phylogenetic signals from HGT and/or paralogy. Thus, we tested several models of varying sophistication on three different datasets, performed removal of fast-evolving or long-branched Archaea and Eukarya, and employed three different strategies to remove compositional heterogeneity to examine their effects on the topological outcome. Our results support a two-domain topology for the tree of life, where Eukarya emerges from within Archaea as sister to a Korarchaeota/Thaumarchaeota (KT) or Crenarchaeota/KT clade for all models under all or at least one of the strategies employed. Taxonomic manipulation allows single-matrix and certain mixture models to vacillate between two-domain and three-domain phylogenies. We find that models vary in their ability to resolve different areas of the tree of life, which does not necessarily correlate with model complexity. For example, both single-matrix and some mixture models recover monophyletic Crenarchaeota and Euryarchaeota archaeal phyla. In contrast, the most sophisticated model recovers a paraphyletic Euryarchaeota but detects two large clades that comprise the Bacteria, which were recovered separately but never together in the other models. Overall, models recovered consistent topologies despite dataset modifications due to the removal of compositional bias, which reflects either ineffective bias reduction or robust datasets that allow models to overcome reconstruction artifacts. We recommend a comparative approach for evolutionary models to identify model weaknesses as well as consensus relationships.

The Evolutionary History of a DNA Methylase Reveals Frequent Horizontal Transfer and Within-Gene Recombination.

Inteins, often referred to as protein introns, are highly mobile genetic elements that invade conserved genes throughout the tree of life. Inteins have been found to invade a wide variety of key genes within actinophages. While in the process of conducting a survey of these inteins in actinophages, we discovered that one protein family of methylases contained a putative intein, and two other unique insertion elements. These methylases are known to occur commonly in phages as orphan methylases (possibly as a form of resistance to restriction-modification systems). We found that the methylase family is not conserved within phage clusters and has a disparate distribution across divergent phage groups. We determined that two of the three insertion elements have a patchy distribution within the methylase protein family. Additionally, we found that the third insertion element is likely a second homing endonuclease, and that all three elements (the intein, the homing endonuclease, and what we refer to as the ShiLan domain) have different insertion sites that are conserved in the methylase gene family. Furthermore, we find strong evidence that both the intein and ShiLan domain are partaking in long-distance horizontal gene transfer events between divergent methylases in disparate phage hosts within the already dispersed methylase distribution. The reticulate evolutionary history of methylases and their insertion elements reveals high rates of gene transfer and within-gene recombination in actinophages.

Neighboring inteins interfere with one another's homing capacity.

Inteins are mobile genetic elements that invade conserved genes across all domains of life and viruses. In some instances, a single gene will have several intein insertion sites. In Haloarchaea, the minichromosome maintenance (MCM) protein at the core of replicative DNA helicase contains four intein insertion sites within close proximity, where two of these sites (MCM-a and MCM-d) are more likely to be invaded. A haloarchaeon that harbors both MCM-a and MCM-d inteins, Haloferax mediterranei, was studied in vivo to determine intein invasion dynamics and the interactions between neighboring inteins. Additionally, invasion frequencies and the conservation of insertion site sequences in 129 Haloferacales mcm homologs were analyzed to assess intein distribution across the order. We show that the inteins at MCM-a and MCM-d recognize and cleave their respective target sites and, in the event that only one empty intein invasion site is present, readily initiate homing (i.e. single homing). However, when two inteins are present co-homing into an intein-free target sequence is much less effective. The two inteins are more effective when invading alleles that already contain an intein at one of the two sites. Our in vivo and computational studies also support that having a proline in place of a serine as the first C-terminal extein residue of the MCM-d insertion site prevents successful intein splicing, but does not stop recognition of the insertion site by the intein's homing endonuclease.

A structural dendrogram of the actinobacteriophage major capsid proteins provides important structural insights into the evolution of capsid stability.

Many double-stranded DNA viruses, including tailed bacteriophages (phages) and herpesviruses, use the HK97-fold in their major capsid protein to make the capsomers of the icosahedral viral capsid. After the genome packaging at near-crystalline densities, the capsid is subjected to a major expansion and stabilization step that allows it to withstand environmental stresses and internal high pressure. Several different mechanisms for stabilizing the capsid have been structurally characterized, but how these mechanisms have evolved is still not understood. Using cryo-EM structure determination of 10 capsids, structural comparisons, phylogenetic analyses, and Alphafold predictions, we have constructed a detailed structural dendrogram describing the evolution of capsid structural stability within the actinobacteriophages. We show that the actinobacteriophage major capsid proteins can be classified into 15 groups based upon their HK97-fold.

Home and away- the evolutionary dynamics of homing endonucleases.

BACKGROUND: Homing endonucleases (HEases) are a large and diverse group of site-specific DNAases. They reside within self-splicing introns and inteins, and promote their horizontal dissemination. In recent years, HEases have been the focus of extensive research due to their promising potential use in gene targeting procedures for the treatment of genetic diseases and for the genetic engineering of crop, animal models and cell lines. RESULTS: Using mathematical analysis and computational modeling, we present here a novel account for the evolution and population dynamics of HEase genes (HEGs). We describe HEGs as paradoxical selfish elements whose long-term persistence in a single population relies on low transmission rates and a positive correlation between transmission efficiency and toxicity. CONCLUSION: Plausible conditions allow HEGs to sustain at high frequency through long evolutionary periods, with the endonuclease frequency being either at equilibrium or periodically oscillating. The predictions of our model may prove important not only for evolutionary theory but also for gene therapy and bio-engineering applications of HEases.

An Efficient, Nonphylogenetic Method for Detecting Genes Sharing Evolutionary Signals in Phylogenomic Data Sets.

Assessing the compatibility between gene family phylogenies is a crucial and often computationally demanding step in many phylogenomic analyses. Here, we describe the Evolutionary Similarity Index (IES), a means to assess shared evolution between gene families using a weighted orthogonal distance regression model applied to sequence distances. The utilization of pairwise distance matrices circumvents comparisons between gene tree topologies, which are inherently uncertain and sensitive to evolutionary model choice, phylogenetic reconstruction artifacts, and other sources of error. Furthermore, IES enables the many-to-many pairing of multiple copies between similarly evolving gene families. This is done by selecting non-overlapping pairs of copies, one from each assessed family, and yielding the least sum of squared residuals. Analyses of simulated gene family data sets show that IES's accuracy is on par with popular tree-based methods while also less susceptible to noise introduced by sequence alignment and evolutionary model fitting. Applying IES to an empirical data set of 1,322 genes from 42 archaeal genomes identified eight major clusters of gene families with compatible evolutionary trends. The most cohesive cluster consisted of 62 genes with compatible evolutionary signal, which occur as both single-copy and multiple homologs per genome; phylogenetic analysis of concatenated alignments from this cluster produced a tree closely matching previously published species trees for Archaea. Four other clusters are mainly composed of accessory genes with limited distribution among Archaea and enriched toward specific metabolic functions. Pairwise evolutionary distances obtained from these accessory gene clusters suggest patterns of interphyla horizontal gene transfer. An IES implementation is available at https://github.com/lthiberiol/evolSimIndex.

The Evolutionary Origins of Extreme Halophilic Archaeal Lineages.

Interest and controversy surrounding the evolutionary origins of extremely halophilic Archaea has increased in recent years, due to the discovery and characterization of the Nanohaloarchaea and the Methanonatronarchaeia. Initial attempts in explaining the evolutionary placement of the two new lineages in relation to the classical Halobacteria (also referred to as Haloarchaea) resulted in hypotheses that imply the new groups share a common ancestor with the Haloarchaea. However, more recent analyses have led to a shift: the Nanohaloarchaea have been largely accepted as being a member of the DPANN superphylum, outside of the euryarchaeota; whereas the Methanonatronarchaeia have been placed near the base of the Methanotecta (composed of the class II methanogens, the Halobacteriales, and Archaeoglobales). These opposing hypotheses have far-reaching implications on the concepts of convergent evolution (distantly related groups evolve similar strategies for survival), genome reduction, and gene transfer. In this work, we attempt to resolve these conflicts with phylogenetic and phylogenomic data. We provide a robust taxonomic sampling of Archaeal genomes that spans the Asgardarchaea, TACK Group, euryarchaeota, and the DPANN superphylum. In addition, we assembled draft genomes from seven new representatives of the Nanohaloarchaea from distinct geographic locations. Phylogenies derived from these data imply that the highly conserved ATP synthase catalytic/noncatalytic subunits of Nanohaloarchaea share a sisterhood relationship with the Haloarchaea. We also employ a novel gene family distance clustering strategy which shows this sisterhood relationship is not likely the result of a recent gene transfer. In addition, we present and evaluate data that argue for and against the monophyly of the DPANN superphylum, in particular, the inclusion of the Nanohaloarchaea in DPANN.

Interaction of two strongly divergent archaellins stabilizes the structure of the Halorubrum archaellum.

Halophilic archaea from the genus Halorubrum possess two extraordinarily diverged archaellin genes, flaB1 and flaB2. To clarify roles for each archaellin, we compared two natural Halorubrum lacusprofundi strains: One of them contains both archaellin genes, and the other has the flaB2 gene only. Both strains synthesize functional archaella; however, the strain, where both archaellins are present, is more motile. In addition, we expressed these archaellins in a Haloferax volcanii strain from which the endogenous archaellin genes were deleted. Three Hfx. volcanii strains expressing Hrr. lacusprofundi archaellins produced functional filaments consisting of only one (FlaB1 or FlaB2) or both (FlaB1/FlaB2) archaellins. All three strains were motile, although there were profound differences in the efficiency of motility. Both native and recombinant FlaB1/FlaB2 filaments have greater thermal stability and resistance to low salinity stress than single-component filaments. Functional supercoiled Hrr. lacusprofundi archaella can be composed of either single archaellin: FlaB2 or FlaB1; however, the two divergent archaellin subunits provide additional stabilization to the archaellum structure and thus adaptation to a wider range of external conditions. Comparative genomic analysis suggests that the described combination of divergent archaellins is not restricted to Hrr. lacusprofundi, but is occurring also in organisms from other haloarchaeal genera.

Ancient horizontal gene transfer can benefit phylogenetic reconstruction.

Although horizontal gene transfer (HGT) is usually considered a disruptive force in recovering organismal phylogeny, it creates important phylogenetic information. In the 'net of life', the recipient of an ancient gene transfer can be the ancestor of a lineage that inherits the transferred gene; thus, the transferred gene marks the recipient and its descendants as a monophyletic group. Ancient gene transfer events can also reveal the order of emergence of donor and recipient lineages. In addition, these ancient events can significantly shape the genetic systems of the recipients and can play a part in their long-term evolution. In this article, we discuss the recent progress in phylogenetic application of ancient HGTs and describe two examples of transfer events to the ancestor of red algae and green plants that support a common origin of these two groups. We also address the potential pitfalls of this application.

Gene transfer: who benefits?

Horizontal gene and genome transfer forces us to recognize that life evolves by fusion as well as bifurcation of lineages, and necessitates the expansion of traditional views of evolution. This chapter reviews the role that horizontal gene transfer (HGT) may play in integrating selection at the gene, species, and community levels. Additionally, we provide an overview of the structure and content of this book, which reflects current thought in the dynamic field of HGT research.

Ancient gene transfer as a tool in phylogenetic reconstruction.

Although horizontal gene transfer (HGT) is often considered as a disruptive force in reconstructing organismal phylogeny, it can also be a valuable phylogenetic tool. A gene in the net of life is often horizontally transferred to the ancestor of a major lineage. If the gene is retained in the recipient and its descendants, it will constitute a shared derived character and mark the recipient and all descendants as a monophyletic group. Additionally, phylogenetically informative HGTs also provide information about the sequence of emergence of involved taxa, because the donor organism must have emerged at least as early as the recipient. Here we review the recent applications of ancient HGT events in reconstructing organismal phylogeny as well as the promise and potential pitfalls of this approach.