Next Article in Journal
Emerging Functions of lncRNA Loci beyond the Transcript Itself
Next Article in Special Issue
Slaying (Yet Again) the Brain-Eating Zombie Called the “Isochore Theory”: A Segmentation Algorithm Used to “Confirm” the Existence of Isochores Creates “Isochores” Where None Exist
Previous Article in Journal
A Novel Synthetic Antibody Library with Complementarity-Determining Region Diversities Designed for an Improved Amplification Profile
Previous Article in Special Issue
On the Base Composition of Transposable Elements
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Review

Current Methods for Recombination Detection in Bacteria

by
Anton E. Shikov
1,2,
Yury V. Malovichko
1,2,
Anton A. Nizhnikov
1,2 and
Kirill S. Antonets
1,2,*
1
Laboratory for Proteomics of Supra-Organismal Systems, All-Russia Research Institute for Agricultural Microbiology (ARRIAM), 196608 St. Petersburg, Russia
2
Faculty of Biology, St. Petersburg State University (SPbSU), 199034 St. Petersburg, Russia
*
Author to whom correspondence should be addressed.
Int. J. Mol. Sci. 2022, 23(11), 6257; https://doi.org/10.3390/ijms23116257
Submission received: 4 March 2022 / Revised: 30 May 2022 / Accepted: 30 May 2022 / Published: 2 June 2022
(This article belongs to the Special Issue New Advances in Compositional Genome Evolution)

Abstract

:
The role of genetic exchanges, i.e., homologous recombination (HR) and horizontal gene transfer (HGT), in bacteria cannot be overestimated for it is a pivotal mechanism leading to their evolution and adaptation, thus, tracking the signs of recombination and HGT events is importance both for fundamental and applied science. To date, dozens of bioinformatics tools for revealing recombination signals are available, however, their pros and cons as well as the spectra of solvable tasks have not yet been systematically reviewed. Moreover, there are two major groups of software. One aims to infer evidence of HR, while the other only deals with horizontal gene transfer (HGT). However, despite seemingly different goals, all the methods use similar algorithmic approaches, and the processes are interconnected in terms of genomic evolution influencing each other. In this review, we propose a classification of novel instruments for both HR and HGT detection based on the genomic consequences of recombination. In this context, we summarize available methodologies paying particular attention to the type of traceable events for which a certain program has been designed.

1. Introduction

The bacterial genome is shaped by homologous recombination (HR) and horizontal or lateral gene transfer (HGT/LGT), with the latter represented by variable molecular mechanisms [1,2]. Recombination could be defined as an exchange of nucleotide sequences between different genomes or within a single genome [1]. If the donor sequence replaces the respective homologous (or homeologous, i.e., similar but not identical) region in the acceptor DNA molecule, then the process is called homologous recombination (HR) [3]. Broadly speaking, HGT could be defined as the incorporation of non-homologous genetic material into the donor genome which requires a long (>500 nucleotides) homologous region flanking the non-homologous segment [2,4]. During the incorporation, a direct RecA-dependent homologous recombination mediates the process, and it includes the excision of the transferred DNA fragment from the donor genome, and its integration into the recipient genome, implying two acts of homologous recombination. HR mostly affects core genes maintaining allelic diversity [5,6], while HGT induces the acquisition of accessory genes [7]. In bioinformatics literature, the term ”non-homologous recombination” (NHR) is sometimes used interchangeably with HGT [4,8], or NHR is seen as HGT-inducing machinery [9,10]; however, that is not always, if ever, true. In fact, DNA integration of mobile genetic elements into the recipient genome such as the integration of phages and genetic islands or conjugative transposons either by site-specific recombinases or by single-strand annealing proteins (SSAPs) requires micro-homologous and homologous sequences, respectively [11,12], that is, strictly speaking, this process could be treated as a type of homologous recombination. Nevertheless, it should be kept in mind that homologous recombination implies DNA strand exchange, whereas the integration processes mentioned do not include strand exchange. Therefore, in the current review by HR, we assume exchange between bacterial genomes and by HGT, we mean the incorporation of genetic material into the recipient genome driven by single-strand annealing (SSA) and/or site-specific recombination but not NHR. HR and HGT are interconnected with respect to the evolutionary dynamics of the bacterial genome. Horizontally transferred genes are often flanked by regions with a high HR rate [13] which could possibly maintain genome size by replacing/eliminating recently acquired genes [13,14]. Gene acquisition, loss, and replacement that are driven by HGT and HR often lead to the emergence of new pathogenic strains [15] and serotypes [16], including opportunistic pathogens [17], increased virulence [18], antibiotic resistance [19,20], immunity evasion [21,22], colonization of new hosts [23], and metabolic adaptations [24,25], thus, affecting public health.
Apart from practical implications, recombination exerts an effect on phylogenetic studies altering almost all trees’ parameters. Models applied in conventional phylogenetic analysis are based on the assumption that any parts of DNA or amino acid sequences determine the evolutionary history in the same way [26]. Nonetheless, if the data contain recombination events, the topologies of trees would differ depending on the part of the sequence, especially if the breakpoint is located in the middle of the sequence [1] which sometimes makes single locus-based phylogeny non-informative [27]. Furthermore, recombination exchange can result in terminal branches that are too long [28], loss of the molecular clock [28], non-uniform distribution of insertions and deletions [29], impossible to identify the common ancestor [30], and an erroneously high dN/dS ratio (the ratio of nonsynonymous to synonymous mutations) resulting in spurious signals of positive selection [31]. Using several housekeeping genes (5–20), namely, MLST (multilocus sequence typing) technique was proposed to overcome these issues; however, it cannot depict gene acquisition or replacement [5]. Progress in next-generation sequencing with high throughput has made it possible to use core genes in the genomes to reconstruct phylogenies, which is known as core genome MLST, or cgMLST. Unfortunately, it still cannot circumvent recombination-driven long terminal branches [32] or inaccurate topologies particularly when the selective pressure is high [33]. A prospective method to obtain trees with correct topology and branch lengths called the coarse-graining approach for phylogenetic reconstruction (CGP) has been devised recently, and it requires further studies to assess its effectiveness [34].
As stated above, HGT and HR are different, yet genomically connected processes. From a genomic perspective, it is virtually impossible to determine specific mechanisms and causes of a particular transfer and/or exchange event; therefore, researchers use indirect computational methods, namely, comparative genomics and phylogeny reconstruction. Here, we analyze state-of-the-art bioinformatics tools for detecting HGT and HR. We discuss conventional approaches as well as novel tools in the context of their pros and cons. We propose an integrated classification of the algorithms based on the ramifications of genetic exchanges, both HGT and HR. Finally, we examine major trends in modern tools’ designing new software and discuss the perspective of further developments.

2. A Brief Overview of Conventional Methods for Detecting Homologous Recombination and Horizontal Gene Transfer

Bioinformatics approaches for detecting genetic exchanges can be divided into several groups depending on the nature of the tasks set, applied algorithms, and genomic consequences that are analyzed. In the existing literature, researchers have separately discussed how to trace homologous recombination and HGT proposing distinct classifications. It is explainable as these two groups seem to have different goals: the former methods are aimed to calculate HR rates and detect chimeric loci in the closely related genomes [3,26], whereas the latter approaches reveal continuous genome regions, for example, genes or larger fragments, acquired from either related or evolutionary distinct species [2].
Considering the end goals of the analysis, methods for HR and HGT detection are divided depending on whether they accomplish: (i) revealing the evidence of exchanges/acquisitions, (ii) identifying mosaic sequences, (iii) finding breakpoint sites, or (iv) calculating recombination and HGT rates [3]. The first task is usually embedded into the latter ones; however, there are some algorithms designed only for revealing the fact of recombination in analyzed sequences applied mostly in HR studies. The second and the third goals are achieved by finding distinct local similarities among a subset of aligned sequences or via the identification of certain loci responsible for phylogenetic incongruences due to the exposure to recombination or horizontal transfer [26]. The last issue is mainly addressed by population genetics principles and phylogenetic analysis [35,36].
When describing the types of methods for HR analysis according to the statistical basis, it should be noted that they belong to so-called parametric and non-parametric methods. The former methods aims to calculate population parameters from a sample [3]. It implies revealing the average recombination frequency, which is achieved by population genetics methods based on a coalescent theory; therefore, these approaches assume the absence of selection and within-group subpopulations and constant population size [3]. The other methods rely on non-parametric statistics inferred directly from sequence alignments and/or tree topology [3]. A distinct methodology is reconstructing ancestral recombination graphs (ARGs) that include elements from all the aforementioned approaches and depict individual recombination events backed by population statistics. The non-parametric methods can be divided into five subclasses on the grounds of their algorithmic nature as follows:
  • Similarity methods are designed to reveal gene conversion by tracking anomalous identity in variable parts of the genome [37];
  • Distance methods find local dissimilarities between sequences using a sliding window technique [38];
  • Compatibility methods detect phylogenetic incongruence of individual sites from alignments and do not require the phylogeny itself [39,40];
  • Substitution distribution approaches group together sequences with similar patterns of integral substitution properties through comparison with the calculated model distribution [41];
  • Phylogenetic methods are based on topological differences between phylogenetic trees, and they represent the most frequently used class of methods in the current studies [42,43,44].
There are three groups of methods for revealing HGT, with two of them being similar to what is applied in HR detection. [2]. The first group is represented by parametric methods, that are aimed to find genetic loci with properties that differ from the genomic average, including GC content [45], oligonucleotide spectrum [46], DNA structure modeling [47], and genomic context [48]. The second group, namely, phylogenetic methods, falls into two subcategories: explicit and implicit phylogenetic methods [2] with the former comparing trees’ topologies and the latter analyzing distances between genomes [2]. The third group examines changes in synteny, i.e., the co-localization of genetic loci in the same regions [49].
As mentioned above, the interconnection between HGT and HR should not be ignored because simultaneous detection of these events can help to disentangle genome evolution. Moreover, the underlying algorithms in described methods are quite similar, and, furthermore, they actually deal with similar, but not opposite, goals, namely, finding loci subjected to recombination/transfer and calculating the frequency of such events. Different classifications do not contradict each other, thus, allowing us to unify them into a combined classification scheme based on the consequences of both HR and HGT (Figure 1). There are three possible scenarios leading to detectable signals in biological data. First, HR and HGT affect the relative positions of genes in the genome through loci gain/loss, repositioning, and duplication, thus, disrupting synteny which is especially conspicuous when comparing whole-genome sequences from diverse strains [49,50]. Second, phylogeny reconstruction based on different loci susceptible to HR or stemming from HGT would cause inconsistencies when collating different gene-based trees or comparing them to those representing species evolution [1,2]. Third, HR and HGT evoke traceable patterns of distributions of genomic properties, namely, single nucleotide polymorphisms (SNPs), alterations in GC-content, etc. [1,2,49]. While there are informative reviews discussing software coupled with guidelines to choose a particular method [1,2,3,26], recently, a lot of new tools have been devised which have not yet been systematically reviewed (Figure 1). Therefore, due to the progress in computational approaches and the occurrence of the novel tools, we discuss them in accordance with the proposed classification in the following section.

3. Current Bioinformatics Tools for Recombination Analysis

3.1. Synteny-Based Methods

Looking from the angle of genomic context, it is possible to find HGT signals in a synteny-aware way. Synteny has been defined as the degree of genomic conservation regarding the relative positions between genes [49]. Hence, changes in synteny can be traced to detect horizontally acquired genes by comparing the order of the loci in the defined genomic interval [49]. The so-called synteny index (SI) was proposed for such purposes and implemented in the Phylo SI software [51]. The synteny index denotes the number of shared gene pairs between most k genes both downstream and upstream of a selected shared ortholog. Then, the average values for all the genes within a pairwise comparison can be utilized to construct a synteny-aware phylogeny [51]. Later on, the SI was incorporated into the nearHGT tool together with constant relative mutability (CRM), another method of calculation that assumes mutation rates to remain constant for each gene within a genome [49]. For two orthologs in two species that exhibit increased similarity with other orthologs diverging in accordance with the mutability model, this approach reports a putative HGT event. Thus, in the beginning, possible HGT candidates are selected through SI calculation, and subsequently, patterns of gene divergence using CRM are defined. In the end, the chi-square test is performed to calculate the significance of the predicted events [49]. A further improvement considers the length of the transfer genes and also utilizes the Chernoff bound test instead of the chi-square test, thus, reducing the number of false-positive calls [50]. The nearHGT program has been applied to evaluate the HGT rate in Mycobacterium leprae, which displayed that pseudogenized loci were transferred with increased frequency in contrast to functional genes [9]. Unfortunately, the available nearHGT program only calculated the probability of HGT for a given set of sequences [49]. The prior steps of calculating the SI index and reporting possible HGT has not been provided as available scripts, thus, nearHGT is more of a conceptual method than a ready-to-use application.
Although other synteny-aware utilities do not report HGT events directly, they can indirectly point out candidates to explore. Lots of genomic browsers have been developed to visualize synteny, namely, BAGET for retrieving syntenic information for a certain gene [52], Synima to juxtapose loci between genomes [53], and SYN-View to investigate antibiotic resistance gene clusters [54]. Sibelia can obtain syntenic blocks in analyzed genomes [55], while SyntTax and SynTracker link them with taxonomical and strain-specific relationships [56,57]. Finally, current pan-genome analysis software now operates with synteny: PEPPAN enables one to retrieve putative HGT events from the accessory genes matrix through synteny-aware pangenome reconstruction [58] and Panaroo provides a graph with syntenic consecutive triplets of gene families, thus, detecting structural variations [59]. Finally, syntenic information could be obtained from gene-to-gene alignments with conventional tools [60,61].

3.2. Phylogenetic Methods

3.2.1. Phylogenetic Methods for HR Detection

One approach to finding present recombination events is called phylogenetic networks. In as much as recombination events lead to intermingling between evolutionally distant lineages, a conventional representation of the evolution as a tree does not reflect the actual phylogenesis. Given that phylogenetic networks pose a more suitable visualization for genetic exchange, there are two distinct types of phylogenetic networks, namely explicit and implicit [62]. The advantage of the former is their interpretability as phylogenetic trees because these networks possess information about parents and recombinants. Unfortunately, explicit networks are hardly obtainable in practical terms, in so far as many recombination events do not provide signals strong enough to distinguish them from mutations, in particular, when they affect conservative genes [26]. In contrast, implicit networks display the most conflicting clades where tree topology is disturbed, demonstrating alternative evolutionary scenarios to be verified with other techniques [62].
Once potential signals are found, it becomes possible to identify breakpoints and to find chimeric sequences. The combination of phylogenetic and distance approaches has revealed these regions that possibly transferred during recombination and the disentangling evolutionary relationships between analyzed sequences regarding these genetic exchanges [26]. Dividing sequences into parts can be carried out by a static procedure with constant borders [63] or dynamically by splitting into two chunks [38], applying a sliding window [41], or more complex heuristics [64]. Parental and recombinant sequences are usually determined by analyzing phylogenetic trees built on different parts of the sequences detected during the previous step. When a potential recombination event is identified, its statistical significance is evaluated, for example, by parametric bootstrap [65] or chi-square distribution [66].
At the moment, the most frequently applied novel programs to examine homologous recombination, as well as HGT, are based on phylogenetic methods. Among these, RDP4 [66] represents a user-friendly application implementing several algorithms with different partitioning schemes for identifying recombined sequences. Its advantages include utilizing a combination of phylogenetic and distance methods providing identification of parent–child relationships and breakpoints in recombined entries [26]. Its updated version, RDP5 [67] has incorporated extra statistical tests, namely, the Φw test [39], the four-gamete test [68], and adapted versions of the homoplasy test [43]. In RDP5, run time speed has been increased up to five times and the number of analyzed bacterial genomes up to 120 times [67]. Still, it cannot handle large batches of bacterial genomes, and therefore, it has been used to trace recombination predominantly in viral genomes, for example, in porcine reproductive and respiratory syndrome virus (PRRSV) [69], SARS-CoV-2 [70], human rhinovirus [71], and feline parvovirus [72]. However, it should be noted that the algorithm inherits limitations of phylogenetic algorithms, the most evident of which is its inability to reveal distant events [26]. Thus, this tool is more suitable for identifying recent events in sequences with moderate divergence and relatively small genomic datasets.
Another group of phylogenetic tools can apply the so-called clonal model [10,44,64,73]. This approach is aimed at scanning whole-genome sequences, in which conservative loci within housekeeping genes are used for phylogeny reconstruction. The chosen genes are considered to depict a clonal frame showing direct relationships between distinct clonal groups.
Gubbins starts with removing SNPs (single nucleotide polymorphisms) that do not fit the assumption of a constant per-site mutation rate, and then places these inconsistencies among the tree built on the remaining polymorphisms [44]. Among its applications, Gubbins has harnessed visualizing and characterizing recombination in Global Pneumococcal Sequence Clusters (GPSCs) [74] and pneumococcal capsular loci [75].
ClonalFrameML uses a pre-reconstructed starting tree and calculates the probability of engaging in recombination for each site using Bayesian maximum likelihood (ML) calculations [73]. ClonalFrameML has been widely used in bacterial genetics to evaluate within-population recombination rate in Prochlorococcus lineages [76], Staphylococcus aureus strains [77], and biosynthetic gene clusters in the Salinispora sp. [78].
Although BratNextGen and fastGEAR are not truly phylogenetic methods, they still operate with clonal relationships, hence, it is more appropriate to discuss them in the current section. However, they do not analyze single nucleotide polymorphisms (SNPs) directly but compare the distributions of variants within clonal lineages using hidden Markov model (HMM) approaches [10,64]. Notably, the latter represents an improvement of the former with higher statistical power. The ability of BratNextGen to reveal ancestral recombination has been applied in studies related to Streptomyces species [79], antibiotic-resistant Staphylococcus aureus strains [80], and differentiated Xylella fastidiosa isolates [81].
On the one hand, all the programs described provide a characterization of SNPs, revealing whether they originate from mutation or recombination, which allows calculating the r/m rate (the probability that a given site stems from recombination rather than mutation) as a proportion of recombination-derived variants. Moreover, these algorithms can handle large datasets due to their high computational capacities. On the other hand, all described tools cannot efficiently distinguish recombination from mutations in the presence of disruptive selection; they also lack statistical power when analyzing highly similar sequences [36]. Another limitation lies in the reliance on phylogenetic trees obtained by methods implying no recombination. Actually, such phylogenetic trees do not portray clonal relationships between ancestors and descendants, as the topology depicts different recombinational rates in diverging bacterial populations rather than sequential evolutionary development [82]. Keeping in mind the questionable feasibility of reflecting clonality even within conservative loci [82], the validity of matching recombination events to the overall phylogeny appears to be dubious. Therefore, it seems more valid to provide per-lineage recombination frequency instead of the overall rate. To sum up, the described tools allow examining large genomic datasets. Ancestral state reconstruction allows them to reveal possible ancestral events particularly optimized in the fastGEAR algorithm [10]. Moreover, due to single-lineage-based clonal relationships, ClonalFrameML [73], Gubbins [44], and BratNextGen [64] are tuned to analyze single bacterial linage with moderate diversity, while fastGEAR harnesses studying interspecies events in sequences with higher diversity [10].

3.2.2. Implicit Phylogenetic Methods to Reveal HGT

In revealing HGT events, explicit phylogenetic methods are presented by straightforward testing of topological similarity [83], decomposing trees’ initial partitions [84], pruning and regrafting subtrees [85], or selecting appropriate reconciliation models accounting for gene loss/duplication and homologous recombination events [86]. Implicit phylogenetic methods do not rely directly on juxtaposing species- and gene-based trees but summarize distances between genomes analyzed to reveal excessively related or different sequences by utilizing BLAST searches [87], disparities between species and gene distances [88], building so-called phylogenetic profiles characterizing patterns of gene presence/absence [89], and clustering polymorphisms [90]. Similar to homologous recombination, novel phylogenetic software to detect horizontal events has been devised recently. It should be noted, however, that the most current tools fall into an implicit category, therefore, these approaches are described here.
HGT-Finder implies a BLAST-based algorithm to provide a set of likely transferred sequences with a transfer index value and significance estimations [91]. The results of the BLAST search against the NCBI non-redundant protein (NCBI-nr) database are utilized to infer relative bit scores (R) calculated as a ratio of the observed bit score to the bit score of the same-sequence alignment. Simultaneously, taxonomic distance using the NCBI Taxonomy database (D) is evaluated as the number of taxonomic units in the query divided by the number of common units with the respective database hit [91]. Then, the transfer index is determined by the mean RD value for each hit genome divided by the number of genomes. Applying HGT-Finder has provided HGT screening in Burkholderia glumae [92] and Aspergillus sp. genomes [91].
HGTector is another tool depending on BLAST searches coupled with taxonomic inference. First, it categorizes genomic hits into three groups: self (the closest strains), close (the same genera or close family), and distal (other families, orders) [93]. The distributions of bit scores for the three categories are then followed by a gene-wise estimation of deviation from these distributions, indicating possible HGT-derived genes [93]. HGTector has been used to infer exchanges in Legionella sp. [94], Nocardia sp. [95], and Blautia sp. [92].
RecentHGT was developed to reveal HGT events between close species [96]. It performs global Needleman–Wunsch alignment of protein-coding sequences and builds the distribution accordingly. Next, particular hits are tested in terms of the inconsistency with the distribution [96]. The approach has successfully harnessed HGTs in Rhizobium strains [96,97].
HGT-Finder and HGTector are more sophisticated taxonomy-wise methods as compared with simple BLAST searches; however, it should be considered that they lack sensitivity as the success of detection depends on taxonomical distance [91,93]. Their design makes them more suitable for revealing HGT between distant bacterial lineages, for example, different taxonomic groups. Contrarily, RecentHGT, in its turn, is designed to detect genetic exchange in close lineages, and therefore can distinguish HGT events from highly conserved housekeeping genes with a reduced false-positive rate as compared with other tools [96].
Of the most current tools to mention, ShadowCaster represents a hybrid approach incorporating both composition-based support vector machines (SVMs) and implicit phylogenetic methods based on the phylogenetic shadow that is constructed on proteomes of species both closely related and distant to the analyzed ones [98]. ShadowCaster shows improved sensitivity as compared with other methods, and moreover, it can detect both close and distant events. For instance, it revealed the transfer of heavy metal resistance genes in Rhodanobacter denitrificans with high accuracy [98]. Nevertheless, while it looks promising, it does not reflect the direction of transfers [98]. As it was not benchmarked by comparing with RecentHGT, it is impossible to state which tool shows better performance, nevertheless, it could be proposed that due to a hybrid check implemented, ShadowCaster may be more sensitive and accurate.

3.3. Methods Based on Genetic Features

3.3.1. Compatibility Methods to Reveal HR

Being non-phylogenetic, compatibility methods now seem of great potential due to their ease and computational effectiveness. The basic approach of such evaluations is a so-called ”four-gamete test” [68]. If two sites provide a genealogy that should involve recurrent mutations to resolve evolutionary relationships, then, these sites are called phylogenetically incompatible, implying their occurrence through homoplasy or recombination [68]. In practice, it is almost impossible to tell recombination from homoplasy for highly similar sequences; nonetheless, one can summarize all homoplasic features and can compare results with the predictions of the model recombination-free distribution [3]. The most commonly used implementations of this approach are the homoplasy test [43] and its improvement, a Φw test [39], both depending on the frequency and distribution of incompatible sites.
The recently developed ptACR program identifies potential breakpoints with a sliding window followed by a permutational test to calculate the significance of found events [40]. Its architecture has ensured robustness to false-positive results checked on clinical isolates of Staphylococcus aureus [40] Nonetheless, ptACR’s disadvantage is the absence of strategies to handle gaps; thus, it is hard, if possible, to analyze divergent sequences with this utility [40], that is to say, this program is useful if the aim of the research is to reveal the most probable recombination events in sequences with moderate diversity.

3.3.2. Substitution Distribution-Based HR Detection Approaches

Similar to compatibility approaches, substitution distribution methods have regained attention due to their high speed as compared with phylogenetic approaches. HREfinder is a dynamic algorithm that divides the genome into blocks where each polymorphism is estimated to result from mutation, homologous recombination, or sequencing error [99]. The stepped validation guarantees obtaining events with high probability as tested in a Xanthomonas oryzae evolution study [100]. The sensitivity of HREfinder continuously grows with sequence diversity, while at the same time, a false-positive rate is coupled with it [99]. Hence, HREfinder just like ptACR, is suitable when dealing with moderately divergent sequences. Within the optimal diversion range, HREfinder detects mostly true events, however, it also tends to miss a lot of them because of detection thresholds [99].

3.3.3. Parametric Methods for HR Identification

Parametric methods are mostly aimed at evaluating the overall HR rate based on population genetics principles. [3]. Population recombination rate (p) is calculated as p = 4 N e r , where Ne is the effective population size and r stands for per-site recombination rate for one generation. Similarly, the population mutation rate is determined by the following equation: θ = 4 N e µ , where µ denotes per-site mutation rate. The p/θ ratio is considered to be an average quantitative variable characterizing recombination for a particular population [1].
One program implementing these methods is Mcorr [101]. This tool calculates the correlation of synonymous substitutions (correlation profiles), and the average recombination rate is delineated on the basis of these profiles [101]. The authors denoted a correlation profile as the probability of observing a difference at the i + l site for a randomly chosen site i, where l is the distance in nucleotides. The function P(l) is constant in the absence of recombination, whereas the presence of recombination causes a monotonic decrease of the P(l) function [101]. The method is highly useful in metagenomic studies, for example, subpopulations in soil metagenome [102] or multidrug-resistant Escherichia coli ST131 populations in the infant gut microbiome [101]. The presented statistic provides a vivid interpretable result reflecting the recombination rate, however, the congruity between this method and compatibility-based HR frequency calculation has not been assessed yet.

3.3.4. Ancestral Recombination Graphs

A distinct method combining phylogenetic incongruence detection, population genetics principles of coalescent theory, and phylogenetic networks is a reconstruction of so-called ancestral recombination graphs (ARGs) [63]. The ARG represents a directed graph in which the most probable site-to-site relationships are exhibited, thus, enabling lateral connections denoting horizontal events such as recombination, which is distinct from classic trees with acyclic topology determined by the average identity between sequences [63]. Being a hybrid approach, ARG construction can depict evolutionary histories that involve recombination coupled with the timed presentation of vertical inheritance, thus, providing a detailed evolution-wise report of recombination events [26].
Bacter, a Bayesian algorithm, has been applied to reconstruct ARGs based on the ClonalOrigin model and Markov chain Monte Carlo (MCMC) algorithm that are used jointly to infer genealogical relationships as well as homologous conversion events and the overall conversion rate [35]. This single-step procedure, instead of a stepped algorithm, improves detection and reduces uncertainty in the case of a poor phylogenetic signal [35]. Its application has accurately revealed previously undetected gene flow between pathogenic and nonpathogenic Escherichia coli serotype O157 representatives [35]. Still, the limitation of this tool is its dependence on a lot of parameters to be optimized for each study, poor throughput, and inefficiency if analyzing long genomes, especially big batches [34].
To handle the inference of ARGs on a large genomic scale, a computationally efficient alternative has been proposed. This approach is called topological data analysis (TDA) in which genomes are treated as points in a high dimensional space with pairwise distances delineated by genetic dissimilarities [103]. Loops in this space linking points occur in the presence of recombination, hence, summarizing loops generate a structure closely related to ARGs, namely, topological ARG (tARG) that depicts minimal recombination histories [103]. TARGet was designed in accordance with the aforementioned principles. While it was tested on eukaryotic organisms, it seems to be applicable for analyzing bacterial genomes, especially when examining large datasets [103]. Topological data analysis is a promising approach regarding its computational effectiveness, although tARG itself cannot depict the specific evolutionary histories behind the data [103]. Therefore, an available tool for recombination-wise bacterial evolution reconstruction, Bacter, is reasonable to apply when dealing with small genomes or parts of genomes, thus, it is necessary to develop computationally efficient tools possibly based on the principles of topological data analysis.

3.3.5. Parametric Methods for Finding HGT Events

Sample-based parametric methods in the context of an HGT analysis have been considered to be less accurate than phylogenetic methods which are dominant in the repertoire of HGT detection programs; however, recently, novel tools with better performance have been devised. They have been applied to obtain the most probable HGT-subjected parts of the genome and the overall transfer frequency. The respective HGT-rate computations rely on the calculation of the HGT-affected genome fraction [104], the ratio of gene gain to gene loss [105], or the total number of detected HGT events divided by the total number of compared genomes [106,107].
To reveal HGT-subjected parts, sequence clustering methods seem to be a perspective approach to deal with the constraints of current tools. The Clusterflock algorithm utilizes a model of self-organizing swarm intelligence originally proposed to imitate bird and insect behavior [108]. This model enables clustering based on a distance matrix with arbitrary distance metrics. The comparison of orthologous gene families’ (OGFs) clusters with obtained flocked clusters has revealed signals of HGT between sequences. Its application has disentangled a large-scale map of genetic exchanges in Staphylococcus aureus [108], still, the Clusterflock has not been benchmarked in the context of comparison with other tools or calculating accuracy and specificity.
The genome mosaic structure (gmos) algorithm was developed to overcome difficulties related to computational costs of full genome-comparison alignments [109]. This program performs local alignments for a given query sequence against subject genomes, refines the alignments according to the substitution models, and finally, overlaps the refined local alignments to gain the mosaic structure of the regions. The utility has been used to track mosaic sequences in the pathogenic Enterococcus faecium strain [109]. The advantage of such an approach is the ability to reveal both homologous recombination events and horizontally transferred genes. However, the latter is possible only if genomes possess sufficient similarity in transfer regions; moreover, the tool does not resolve the direction of transfer/exchange [109].
GeneMates is an R package to reveal co-transferred genes in bacterial genomes associated with mobile genetic elements [110]. In the package, the matrix of core genome SNPs coupled with allelic presence/absence matrix is analyzed using linear mixed models to generate a network of alleles that are most likely co-transferred together. This framework transcends simple co-occurrence tests according to a validation study of GeneMates on known antibiotic resistance genes in Escherichia coli and Salmonella Typhimurium; nonetheless, it is designed for a specific aim to identify intraspecies events, while its dependence on core SNPs may probably restrict the sensitivity of the analysis [110].
The abovementioned tools rely on completed and maximally annotated genomes. In contrast to it, Daisy is a reference-free method that processes short reads to detect HGT boundaries via split-read mapping and coverage information, and it leads to outperforming assembly-based approaches [111]. Its performance has been checked on a simulated H. pylori dataset and two real E. coli datasets [111]. While providing high sensitivity, Daisy relies on short reads only and requires genomes with explicitly defined suspected donor and acceptor, thus, it is not applicable to process long reads or it cannot compare bathes of genomes when donor and acceptor are unknown.

4. Assessing the Effectiveness of Recombination Detection Software

To choose a particular algorithm to detect HR and HGT in biological data, it is useful to understand the expected rate of false-positive calls. Erroneous identification of recombination events may occur when analyzing extremely divergent sequences, given that in the tools applied, statistical power proportionally increases with sequence divergence [112]. However, handling substantially similar strains may also generate errors [113]. Some methods are also sensitive to asymmetric tree topology [112]. If linkage disequilibrium between nucleotide substitutions is used to predict recombination events, findings may actually represent evolutionary selection signals instead of genetic exchange [114]. A so-called “patchy-tachy” (PT) phenomenon describes sequences in which different partitions exhibit unequal evolution rates, which leads to an excess in false-positive results [115]. Tracking HGT can generate false-positive results as well. For instance, parametric methods based on codon usage are prone to a high rate of both false-positive and false-negative results [116]. In addition, similar to HR, false-positive HGT signals likely occur if comparing closely related strains [49]. Another essential source of misreported events relates to genomic data collection, namely, assembly procedures and PCR-gained chimeric sequences. For example, a comparative study of Mycobacterium tuberculosis genomes revealed that most of the recombination events described in the literature were artifacts [117]. They occurred due to inconsistencies in the genomic alignments in the case of reference-based genome assembly relying on the reference assembly already containing false-positive results; hence, in bacterial genomics, high-quality de novo assemblies should be preferred instead [117]. Sample preparation could provoke artificial recombination events both during PCR amplification and data analysis of sequencing data leading to the emergence of chimeric sequences [118,119]. These chimeric sequences are often presented in current databases, thus, making it difficult, if possible, to estimate the number of artefactual data possibly utilized as reference sequences in phylogenetic studies [26].
Given a great variety of cases in which correct detection of HGT and HR is hampered (Table 1), the limits of applications for the programs have to be quantitatively evaluated to ensure choosing the most accurate and sensitive algorithms. Therefore, it seems surprising that there is a lack of comparative analyses. In most cases, such studies include only a small number of algorithms to display the performance of the recently devised tool [10,44,98], whereas comprehensive examinations currently seem outdated [112,120]. Still, for such performance tests, one can apply genome evolution simulators under HR, such as SimBac [121] and Bacmeta [122]. Nevertheless, it should be borne in mind that these simulators are coalescent-based, implying a constant recombination rate and modeling neutral evolution. In contrast, cutting-edge technologies such as CoreSimul [123] include stochastic parameters imitating environmental changes accompanied by recombination. Similar to it, there are HGT simulators such as HgtSIM [124]. Finally, the most promising simulators capable of modeling both recombination and horizontal exchange such as SLiM [125] can be utilized to jointly analyze the detection of both HR and NHR, thus, providing a comprehensive evaluation of the genetic exchange map between bacterial populations.

5. Conclusions

Homologous recombination (HR) and horizontal gene transfer (HGT) in bacteria are fundamental mechanisms of their evolution, and these two processes are inextricably connected on a genomic scale. HR provides allelic diversity and causes genetic gain/loss [13]. It may well maintain genome stability by discarding unused HGT-obtained genes, and sometimes this gene loss intensity does not correlate to the overall HR rate [127]. HR and HGT are of importance for fundamental science and practical application. Therefore, genomic studies require special tools for the effective detection of these events. Recently, a host of programs have been devised, and the development is still going on. Having reviewed novel bioinformatics tools, we revealed that methods depend on the consequences of HR and HGT such as alterations in synteny, trees’ topologies incongruence, and altered distribution of genetic features (Figure 1). A great variety of available programs presents dozens of applications for studies with different goals and varying performances when used on diverse data. Programs such as Mcorr [101] or clonal frame model-based tools [44,64,73] can calculate overall HR rate, while nearHGT can evaluate HGT rate [49]. ARGs implemented in Bacter [35] are tuned to depict site-wise individual HR histories, thus, being computationally expensive, sensitive to divergence, and applicable for analyzing small sets of related genomes. Parent–child relationships for large blocks are also provided by RDP4/5 [66,67] in the case of HR, and similar donor-acceptor HGT directions could be identified with Daisy [111]. The tools also differentiate in preferred data to process. ClonalFrameML [73], Gubbins [44], and RDP4/5 [66,67] manage to detect recent HR events in moderately divergent sequences, while fastGEAR [10] is suitable for digging ancestral and recent recombination events in sequences with high divergence. If highly accurate detection of true recombination events is needed, ptACR [40] and HREfinder [99] seem to be useful, while, at the same time, they lack sensitivity. Similar to HR, RecentHGT [96] shows a lower false-positive rate being appropriately utilized to uncover recent transfers in similar sequences, whereas HGT-Finder [91] and HGTector [93] are tuned to trace events in distant genomes. Similar to fastGEAR [10], ShadowCaster [98] predicts both distant and close HGT events and potentially appears to be the most effective HGT-detecting tool by far. To sum up, state-of-the-art approaches for studying HR and HGT are characterized by different sensitivities and accuracies, and they find either recent or ancient events in similar, moderately different, or highly divergent sequences. We might conclude, that the tools reviewed show better performance when detecting some types of recombination events while being less effective to reveal others. Therefore, it looks promising to develop new software that incorporates hybrid approaches to improve recombination detection. Going further, given the genomic interrelation between HR and HGT affecting each other in terms of frequency and direction, a comprehensive framework equipped with both HR and HGT predictors would sufficiently broaden our understanding of the mechanisms driving the plasticity of bacterial genomes.

Author Contributions

Conceptualization, A.E.S. and K.S.A.; writing—original draft preparation, A.E.S., writing—review and editing, A.E.S., Y.V.M., A.A.N. and K.S.A.; visualization, A.E.S.; supervision, A.A.N. and K.S.A.; project administration, A.A.N. and K.S.A.; funding acquisition, K.S.A. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by the Russian Science Foundation (20-76-10044 to K.S.A.).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare that there is no conflict of interest regarding the publication of this paper.

Abbreviations

HRHomologous recombination
NHRNon-homologous recombination
HGTHorizontal gene transfer
LGTLateral gene transfer
SSAPsSingle-strand annealing proteins
SSASingle-strand annealing
ARGsAncestral recombination graphs
MLSTMultilocus sequence typing
CGPCoarse-graining approach for phylogenetic reconstruction
CRMConstant relative mutability
SISynteny index
PRRSVPorcine reproductive and respiratory syndrome virus
GPSCsGlobal pneumococcal sequence clusters
SNPsSingle nucleotide polymorphisms
HMMHidden Markov model
TDATopological data analysis
MCMCMarkov chain Monte Carlo
tARGTopological ARG
OGFsOrthologous gene families

References

  1. Posada, D.; Crandall, K.A.; Holmes, E.C. Recombination in evolutionary genomics. Annu. Rev. Genet. 2002, 36, 75–97. [Google Scholar] [CrossRef] [Green Version]
  2. Ravenhall, M.; Škunca, N.; Lassalle, F.; Dessimoz, C. Inferring horizontal gene transfer. PLoS Comput. Biol. 2015, 11, e1004095. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  3. Lemey, P.; Posada, D. Introduction to recombination detection. In The Phylogenetic Handbook; Vandamme, A.-M., Salemi, M., Lemey, P., Eds.; Cambridge University Press: Cambridge, UK, 2012; pp. 493–518. ISBN 9780511819049. [Google Scholar]
  4. Vos, M. Why do bacteria engage in homologous recombination? Trends Microbiol. 2009, 17, 226–232. [Google Scholar] [CrossRef] [PubMed]
  5. Didelot, X.; Maiden, M.C.J. Impact of recombination on bacterial evolution. Trends Microbiol. 2010, 18, 315–322. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  6. Cheng, K.; Rong, X.; Huang, Y. Widespread interspecies homologous recombination reveals reticulate evolution within the genus Streptomyces. Mol. Phylogenet. Evol. 2016, 102, 246–254. [Google Scholar] [CrossRef]
  7. Ochman, H.; Lawrence, J.G.; Groisman, E.A. Lateral gene transfer and the nature of bacterial innovation. Nature 2000, 405, 299–304. [Google Scholar] [CrossRef]
  8. Lassalle, F.; Didelot, X. Bacterial microevolution and the pangenome. In The Pangenome: Diversity, Dynamics and Evolution of Genomes; Springer: Cham, Switzerland, 2020; pp. 129–149. [Google Scholar] [CrossRef]
  9. Avni, E.; Montoya, D.; Lopez, D.; Modlin, R.; Pellegrini, M.; Snir, S. A phylogenomic study quantifies competing mechanisms for pseudogenization in prokaryotes. The Mycobacterium leprae case. PLoS ONE 2018, 13, e0204322. [Google Scholar] [CrossRef] [Green Version]
  10. Mostowy, R.; Croucher, N.J.; Andam, C.P.; Corander, J.; Hanage, W.P.; Marttinen, P. Efficient Inference of Recent and Ancestral Recombination within Bacterial Populations. Mol. Biol. Evol. 2017, 34, 1167–1182. [Google Scholar] [CrossRef] [Green Version]
  11. Steczkiewicz, K.; Prestel, E.; Bidnenko, E.; Szczepankowska, A.K. Expanding Diversity of Firmicutes Single-Strand Annealing Proteins: A Putative Role of Bacteriophage-Host Arms Race. Front. Microbiol. 2021, 12, 644622. [Google Scholar] [CrossRef]
  12. Subramaniam, S.; Erler, A.; Fu, J.; Kranz, A.; Tang, J.; Gopalswamy, M.; Ramakrishnan, S.; Keller, A.; Grundmeier, G.; Müller, D.; et al. DNA annealing by Redβ is insufficient for homologous recombination and the additional requirements involve intra- and intermolecular interactions. Sci. Rep. 2016, 6, 34525. [Google Scholar] [CrossRef]
  13. Iranzo, J.; Wolf, Y.I.; Koonin, E.V.; Sela, I. Gene gain and loss push prokaryotes beyond the homologous recombination barrier and accelerate genome sequence divergence. Nat. Commun. 2019, 10, 5376. [Google Scholar] [CrossRef] [Green Version]
  14. Ely, B. Recombination and gene loss occur simultaneously during bacterial horizontal gene transfer. PLoS ONE 2020, 15, 4–6. [Google Scholar] [CrossRef]
  15. Levin, B.R.; Cornejo, O.E. The population and evolutionary dynamics of homologous gene recombination in bacteria. PLoS Genet. 2009, 5, e1000601. [Google Scholar] [CrossRef] [Green Version]
  16. Gürtler, V.; Mayall, B.C. Genomic approaches to typing, taxonomy and evolution of bacterial isolates. Int. J. Syst. Evol. Microbiol. 2001, 51, 3–16. [Google Scholar] [CrossRef] [Green Version]
  17. Aujoulat, F.; Romano-Bertrand, S.; Masnou, A.; Marchandin, H.; Jumas-Bilak, E. Niches, population structure and genome reduction in Ochrobactrum intermedium: Clues to technology-driven emergence of pathogens. PLoS ONE 2014, 9, e0171448. [Google Scholar] [CrossRef] [Green Version]
  18. Hao, L.; Holden, M.T.G.; Wang, X.; Andrew, L.; Wellnitz, S.; Hu, F.; Whaley, M.; Sammons, S.; Knipe, K.; Frace, M.; et al. Distinct evolutionary patterns of Neisseria meningitidis serogroup B disease outbreaks at two universities in the USA. Microb. Genom. 2018, 4, 1–10. [Google Scholar] [CrossRef] [Green Version]
  19. Nudel, K.; Zhao, X.; Basu, S.; Dong, X.; Hoffmann, M.; Feldgarden, M.; Allard, M.; Klompas, M.; Bry, L. Genomics of Corynebacterium striatum, an emerging multidrug-resistant pathogen of immunocompromised patients. Clin. Microbiol. Infect. 2018, 24, 1016.e7–1016.e13. [Google Scholar] [CrossRef] [Green Version]
  20. Liu, L.; Cui, Y.; Zheng, B.; Jiang, S.; Yu, W.; Shen, P.; Ji, J.; Li, L.; Qin, N.; Xiao, Y. Analysis of tigecycline resistance development in clinical Acinetobacter baumannii isolates through a combined genomic and transcriptomic approach. Sci. Rep. 2016, 6, 1–12. [Google Scholar] [CrossRef] [Green Version]
  21. Štaudová, B.; Strouhal, M.; Zobaníková, M.; Čejková, D.; Fulton, L.L.; Chen, L.; Giacani, L.; Centurion-Lara, A.; Bruisten, S.M.; Sodergren, E.; et al. Whole Genome Sequence of the Treponema pallidum subsp. endemicum Strain Bosnia A: The Genome Is Related to Yaws Treponemes but Contains Few Loci Similar to Syphilis Treponemes. PLoS Negl. Trop. Dis. 2014, 8. [Google Scholar] [CrossRef]
  22. Guo, Q.; Mustapha, M.M.; Chen, M.; Qu, D.; Zhang, X.; Chen, M.; Doi, Y.; Wang, M.; Harrison, L.H. Evolution of sequence type 4821 clonal complex meningococcal strains in China from prequinolone to quinolone era, 1972–2013. Emerg. Infect. Dis. 2018, 24, 683–690. [Google Scholar] [CrossRef] [Green Version]
  23. Potnis, N.; Kandel, P.P.; Merfa, M.V.; Retchless, A.C.; Parker, J.K.; Stenger, D.C.; Almeida, R.P.P.; Bergsma-Vlami, M.; Westenberg, M.; Cobine, P.A.; et al. Patterns of inter- and intrasubspecific homologous recombination inform eco-evolutionary dynamics of Xylella fastidiosa. ISME J. 2019, 13, 2319–2333. [Google Scholar] [CrossRef]
  24. Rounge, T.B.; Rohrlack, T.; Kristensen, T.; Jakobsen, K.S. Recombination and selectional forces in cyanopeptolin NRPS operons from highly similar, but geographically remote Planktothrix strains. BMC Microbiol. 2008, 8, 1–10. [Google Scholar] [CrossRef] [Green Version]
  25. Bosch, R.; García-Valdés, E.; Moore, E.R.B. Complete nucleotide sequence and evolutionary significance of a chromosomally encoded naphthalene-degradation lower pathway from Pseudomonas stutzeri AN10. Gene 2000, 245, 65–74. [Google Scholar] [CrossRef]
  26. Martin, D.P.; Lemey, P.; Posada, D. Analysing recombination in nucleotide sequences. Mol. Ecol. Resour. 2011, 11, 943–955. [Google Scholar] [CrossRef]
  27. Archibald, J.M.; Roger, A.J. Gene duplication and gene conversion shape the evolution of archaeal chaperonins. J. Mol. Biol. 2002, 316, 1041–1050. [Google Scholar] [CrossRef]
  28. Schierup, M.H.; Hein, J. Consequences of recombination on traditional phylogenetic analysis. Genetics 2000, 156, 879–891. [Google Scholar] [CrossRef]
  29. Gribaldo, S.; Philippe, H. Ancient phylogenetic relationships. Theor. Popul. Biol. 2002, 61, 391–408. [Google Scholar] [CrossRef]
  30. Arenas, M.; Posada, D. The effect of recombination on the reconstruction of ancestral sequences. Genetics 2010, 184, 1133–1139. [Google Scholar] [CrossRef] [Green Version]
  31. Shriner, D.; Nickle, D.C.; Jensen, M.A.; Mullins, J.I. Potential impact of recombination on sitewise approaches for detecting positive natural selection. Genet. Res. 2003, 81, 115–121. [Google Scholar] [CrossRef] [Green Version]
  32. Hedge, J.; Wilson, D.J. Bacterial phylogenetic reconstruction from whole genomes is robust to recombination but demographic inference is not. mBio 2014, 5, 5–8. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  33. Stott, C.M.; Bobay, L.M. Impact of homologous recombination on core genome phylogenies. BMC Genom. 2020, 21, 1–10. [Google Scholar] [CrossRef] [PubMed]
  34. Pang, T.Y. A coarse-graining, ultrametric approach to resolve the phylogeny of prokaryotic strains with frequent homologous recombination. BMC Evol. Biol. 2020, 20, 1–13. [Google Scholar] [CrossRef] [PubMed]
  35. Vaughan, T.G.; Welch, D.; Drummond, A.J.; Biggs, P.J.; George, T.; French, N.P. Inferring ancestral recombination graphs from bacterial genomic data. Genetics 2017, 205, 857–870. [Google Scholar] [CrossRef] [Green Version]
  36. Hanage, W.P. Not so simple after all: Bacteria, their population genetics, and recombination. Cold Spring Harb. Perspect. Biol. 2016, 8, 1–18. [Google Scholar] [CrossRef] [Green Version]
  37. Ohta, T.; Basten, C.J. Gene conversion generates hypervariability at the variable regions of kallikreins and their inhibitors. Mol. Phylogenet. Evol. 1992, 1, 87–90. [Google Scholar] [CrossRef]
  38. Weiller, G.F. Phylogenetic profiles: A graphical method for detecting genetic recombinations in homologous sequences. Mol. Biol. Evol. 1998, 15, 326–335. [Google Scholar] [CrossRef]
  39. Bruen, T.C.; Philippe, H.; Bryant, D. A simple and robust statistical test for detecting the presence of recombination. Genetics 2006, 172, 2665–2681. [Google Scholar] [CrossRef] [Green Version]
  40. Lai, Y.P.; Ioerger, T.R. A statistical method to identify recombination in bacterial genomes based on SNP incompatibility. BMC Bioinform. 2018, 19, 450. [Google Scholar] [CrossRef]
  41. Gibbs, M.J.; Armstrong, J.S.; Gibbs, A.J. Sister-scanning: A Monte Carlo procedure for assessing signals in recombinant sequences. Bioinformatics 2000, 16, 573–582. [Google Scholar] [CrossRef]
  42. Taylor, J.C.; Martin, H.C.; Lise, S.; Broxholme, J.; Cazier, J.-B.; Rimmer, A.; Kanapin, A.; Lunter, G.; Fiddy, S.; Allan, C.; et al. Factors influencing success of clinical genome sequencing across a broad spectrum of disorders. Nat. Genet. 2015, 47, 717–726. [Google Scholar] [CrossRef]
  43. Maynard Smith, J.; Smith, N.H. Detecting recombination from gene trees. Mol. Biol. Evol. 1998, 15, 590–599. [Google Scholar] [CrossRef]
  44. Croucher, N.J.; Page, A.J.; Connor, T.R.; Delaney, A.J.; Keane, J.A.; Bentley, S.D.; Parkhill, J.; Harris, S.R. Rapid phylogenetic analysis of large samples of recombinant bacterial whole genome sequences using Gubbins. Nucleic Acids Res. 2015, 43, e15. [Google Scholar] [CrossRef] [Green Version]
  45. Daubin, V.; Lerat, E.; Perrière, G. The source of laterally transferred genes in bacterial genomes. Genome Biol. 2003, 4, R57. [Google Scholar] [CrossRef] [Green Version]
  46. Vernikos, G.S.; Parkhill, J. Interpolated variable order motifs for identification of horizontally acquired DNA: Revisiting the Salmonella pathogenicity islands. Bioinformatics 2006, 22, 2196–2203. [Google Scholar] [CrossRef] [Green Version]
  47. Worning, P.; Jensen, L.J.; Nelson, K.E.; Brunak, S.; Ussery, D.W. Structural analysis of DNA sequence: Evidence for lateral gene transfer in Thermotoga maritima. Nucleic Acids Res. 2000, 28, 706–709. [Google Scholar] [CrossRef] [Green Version]
  48. Vernikos, G.S.; Parkhill, J. Resolving the structural features of genomic islands: A machine learning approach. Genome Res. 2008, 18, 331–342. [Google Scholar] [CrossRef] [Green Version]
  49. Adato, O.; Ninyo, N.; Gophna, U.; Snir, S. Detecting horizontal gene transfer between closely related taxa. PLoS Comput. Biol. 2015, 11, e1004408. [Google Scholar] [CrossRef] [Green Version]
  50. Sevillya, G.; Adato, O.; Snir, S. Detecting horizontal gene transfer: A probabilistic approach. BMC Genom. 2020, 21, 106. [Google Scholar] [CrossRef]
  51. Shifman, A.; Ninyo, N.; Gophna, U.; Snir, S. Phylo SI: A new genome-wide approach for prokaryotic phylogeny. Nucleic Acids Res. 2014, 42, 2391–2404. [Google Scholar] [CrossRef]
  52. Hepp, B.; Da Cunha, V.; Lorieux, F.; Oberto, J. BAGET 2.0: An updated web tool for the effortless retrieval of prokaryotic gene context and sequence. Bioinformatics 2021, 37, 2750–2752. [Google Scholar] [CrossRef]
  53. Farrer, R.A. Synima: A Synteny imaging tool for annotated genome assemblies. BMC Bioinform. 2017, 18, 507. [Google Scholar] [CrossRef]
  54. Stahlecker, J.; Mingyar, E.; Ziemert, N.; Mungan, M.D. SYN-View: A Phylogeny-Based Synteny Exploration Tool for the Identification of Gene Clusters Linked to Antibiotic Resistance. Molecules 2021, 26, 144. [Google Scholar] [CrossRef]
  55. Minkin, I.; Patel, A.; Kolmogorov, M.; Vyahhi, N.; Pham, S. Sibelia: A Scalable and Comprehensive Synteny Block Generation Tool for Closely Related Microbial Genomes. In Proceedings of the Algorithms in Bioinformatics; Darling, A., Stoye, J., Eds.; Springer: Berlin/Heidelberg, Germany, 2013; pp. 215–229. [Google Scholar]
  56. Oberto, J. SyntTax: A web server linking synteny to prokaryotic taxonomy. BMC Bioinform. 2013, 14, 4. [Google Scholar] [CrossRef] [Green Version]
  57. Enav, H.; Ley, R.E. SynTracker: A synteny based tool for tracking microbial strains. bioRxiv 2021. [Google Scholar] [CrossRef]
  58. Zhou, Z.; Charlesworth, J.; Achtman, M. Accurate reconstruction of bacterial pan- and core genomes with PEPPAN. Genome Res. 2020, 30, 1667–1679. [Google Scholar] [CrossRef]
  59. Tonkin-Hill, G.; MacAlasdair, N.; Ruis, C.; Weimann, A.; Horesh, G.; Lees, J.A.; Gladstone, R.A.; Lo, S.; Beaudoin, C.; Floto, R.A.; et al. Producing polished prokaryotic pangenomes with the Panaroo pipeline. Genome Biol. 2020, 21, 180. [Google Scholar] [CrossRef]
  60. Altschul, S.F.; Madden, T.L.; Schäffer, A.A.; Zhang, J.; Zhang, Z.; Miller, W.; Lipman, D.J. Gapped BLAST and PSI-BLAST: A new generation of protein database search programs. Nucleic Acids Res. 1997, 25, 3389–3402. [Google Scholar] [CrossRef] [Green Version]
  61. Delcher, A.L.; Salzberg, S.L.; Phillippy, A.M. Using MUMmer to identify similar regions in large sequence sets. Curr. Protoc. Bioinform. 2003, 1, 10.3. [Google Scholar] [CrossRef]
  62. Tan, M.; Long, H.; Liao, B.; Cao, Z.; Yuan, D.; Tian, G.; Zhuang, J.; Yang, J. QS-Net: Reconstructing phylogenetic networks based on quartet and sextet. Front. Genet. 2019, 10, 1–9. [Google Scholar] [CrossRef] [Green Version]
  63. Bloomquist, E.W.; Suchard, M.A. Unifying vertical and nonvertical evolution: A stochastic arg-based framework. Syst. Biol. 2010, 59, 27–41. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  64. De Been, M.; Van Schaik, W.; Cheng, L.; Corander, J.; Willems, R.J. Recent recombination events in the core genome are associated with adaptive evolution in Enterococcus faecium. Genome Biol. Evol. 2013, 5, 1524–1535. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  65. Milne, I.; Lindner, D.; Bayer, M.; Husmeier, D.; McGuire, G.; Marshall, D.F.; Wright, F. TOPALi v2: A rich graphical interface for evolutionary analyses of multiple alignments on HPC clusters and multi-core desktops. Bioinformatics 2009, 25, 126–127. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  66. Martin, D.P.; Murrell, B.; Golden, M.; Khoosal, A.; Muhire, B. RDP4: Detection and analysis of recombination patterns in virus genomes. Virus Evol. 2015, 1, 1–5. [Google Scholar] [CrossRef] [Green Version]
  67. Martin, D.P.; Varsani, A.; Roumagnac, P.; Botha, G.; Maslamoney, S.; Schwab, T.; Kelz, Z.; Kumar, V.; Murrell, B. RDP5: A computer program for analyzing recombination in, and removing signals of recombination from, nucleotide sequence datasets. Virus Evol. 2021, 7, veaa087. [Google Scholar] [CrossRef]
  68. Hudson, R.R.; Kaplan, N.L. Statistical properties of the number of recombination events in the history of a sample of DNA sequences. Genetics 1985, 111, 147–164. [Google Scholar] [CrossRef]
  69. Vandenbussche, F.; Mathijs, E.; Tignon, M.; Vandersmissen, T.; Cay, A.B. WGS- versus ORF5-Based Typing of PRRSV: A Belgian Case Study. Viruses 2021, 13, 2419. [Google Scholar] [CrossRef]
  70. Islam, A.; Ferdous, J.; Sayeed, M.A.; Islam, S.; Kaisar Rahman, M.; Abedin, J.; Saha, O.; Hassan, M.M.; Shirin, T. Spatial epidemiology and genetic diversity of SARS-CoV-2 and related coronaviruses in domestic and wild animals. PLoS ONE 2021, 16, e0260635. [Google Scholar] [CrossRef]
  71. Luka, M.M.; Kamau, E.; de Laurent, Z.R.; Morobe, J.M.; Alii, L.K.; Nokes, D.J.; Agoti, C.N. Whole genome sequencing of two human rhinovirus A types (A101 and A15) detected in Kenya, 2016–2018. Wellcome Open Res. 2021, 6, 178. [Google Scholar] [CrossRef]
  72. Tucciarone, C.M.; Franzo, G.; Legnardi, M.; Lazzaro, E.; Zoia, A.; Petini, M.; Furlanello, T.; Caldin, M.; Cecchinato, M.; Drigo, M. Genetic Insights into Feline Parvovirus: Evaluation of Viral Evolutionary Patterns and Association between Phylogeny and Clinical Variables. Viruses 2021, 13, 1033. [Google Scholar] [CrossRef]
  73. Didelot, X.; Wilson, D.J. ClonalFrameML: Efficient Inference of Recombination in Whole Bacterial Genomes. PLoS Comput. Biol. 2015, 11, 1–18. [Google Scholar] [CrossRef] [Green Version]
  74. Gladstone, R.A.; Lo, S.W.; Goater, R.; Yeats, C.; Taylor, B.; Hadfield, J.; Lees, J.A.; Croucher, N.J.; van Tonder, A.J.; Bentley, L.J.; et al. Visualizing variation within Global Pneumococcal Sequence Clusters (GPSCs) and country population snapshots to contextualize pneumococcal isolates. Microb. Genom. 2020, 6, e000357. [Google Scholar] [CrossRef]
  75. Lo, S.W.; Gladstone, R.A.; van Tonder, A.J.; Du Plessis, M.; Cornick, J.E.; Hawkins, P.A.; Madhi, S.A.; Nzenze, S.A.; Kandasamy, R.; Ravikumar, K.L.; et al. A mosaic tetracycline resistance gene tet(S/M) detected in an MDR pneumococcal CC230 lineage that underwent capsular switching in South Africa. J. Antimicrob. Chemother. 2020, 75, 512–520. [Google Scholar] [CrossRef]
  76. Chen, Z.; Wang, X.; Song, Y.; Zeng, Q.; Zhang, Y.; Luo, H. Prochlorococcus have low global mutation rate and small effective population size. Nat. Ecol. Evol. 2021, 6, 183–194. [Google Scholar] [CrossRef]
  77. Gill, J.L.; Hedge, J.; Wilson, D.J.; MacLean, R.C. Evolutionary Processes Driving the Rise and Fall of Staphylococcus aureus ST239, a Dominant Hybrid Pathogen. mBio 2021, 12, e0216821. [Google Scholar] [CrossRef]
  78. Chase, A.B.; Sweeney, D.; Muskat, M.N.; Guillén-Matus, D.G.; Jensen, P.R. Vertical Inheritance Facilitates Interspecies Diversification in Biosynthetic Gene Clusters and Specialized Metabolites. mBio 2021, 12, e0270021. [Google Scholar] [CrossRef]
  79. Wang, J.; Li, Y.; Pinto-Tomás, A.A.; Cheng, K.; Huang, Y. Habitat Adaptation Drives Speciation of a Streptomyces Species with Distinct Habitats and Disparate Geographic Origins. mBio 2022, 13, e0278121. [Google Scholar] [CrossRef]
  80. Sawhney, S.S.; Ransom, E.M.; Wallace, M.A.; Reich, P.J.; Dantas, G.; Burnham, C.-A.D. Comparative Genomics of Borderline Oxacillin-Resistant Staphylococcus aureus Detected during a Pseudo-outbreak of Methicillin-Resistant S. aureus in a Neonatal Intensive Care Unit. mBio 2022, 13, e0319621. [Google Scholar] [CrossRef]
  81. Castillo, A.I.; Tsai, C.-W.; Su, C.-C.; Weng, L.-W.; Lin, Y.-C.; Cho, S.-T.; Almeida, R.P.P.; Kuo, C.-H. Genetic differentiation of Xylella fastidiosa following the introduction into Taiwan. Microb. Genom. 2021, 7, 727. [Google Scholar] [CrossRef]
  82. Sakoparnig, T.; Field, C.; van Nimwegen, E. Whole genome phylogenies reflect the distributions of recombination rates for many bacterial species. eLife 2021, 10, e65366. [Google Scholar] [CrossRef]
  83. Lerat, E.; Daubin, V.; Moran, N.A. From Gene Trees to Organismal Phylogeny in Prokaryotes: The Case of the γ-Proteobacteria. PLoS Biol. 2003, 1, e19. [Google Scholar] [CrossRef]
  84. Zhaxybayeva, O.; Gogarten, J.P.; Charlebois, R.L.; Doolittle, W.F.; Papke, R.T. Phylogenetic analyses of cyanobacterial genomes: Quantification of horizontal gene transfer events. Genome Res. 2006, 16, 1099–1108. [Google Scholar] [CrossRef] [Green Version]
  85. Baroni, M.; Grünewald, S.; Moulton, V.; Semple, C. Bounding the number of hybridisation events for a consistent evolutionary history. J. Math. Biol. 2005, 51, 171–182. [Google Scholar] [CrossRef] [Green Version]
  86. Szöllosi, G.J.; Boussau, B.; Abby, S.S.; Tannier, E.; Daubin, V. Phylogenetic modeling of lateral gene transfer reconstructs the pattern and relative timing of speciations. Proc. Natl. Acad. Sci. USA 2012, 109, 17513–17518. [Google Scholar] [CrossRef] [Green Version]
  87. Nelson, K.E.; Clayton, R.A.; Gill, S.R.; Gwinn, M.L.; Dodson, R.J.; Haft, D.H.; Hickey, E.K.; Peterson, J.D.; Nelson, W.C.; Ketchum, K.A.; et al. Evidence for lateral gene transfer between Archaea and bacteria from genome sequence of Thermotoga maritima. Nature 1999, 399, 323–329. [Google Scholar] [CrossRef] [PubMed]
  88. Clarke, G.D.P.; Beiko, R.G.; Ragan, M.A.; Charlebois, R.L. Inferring genome trees by using a filter to eliminate phylogenetically discordant sequences and a distance matrix based on mean normalized BLASTP scores. J. Bacteriol. 2002, 184, 2072–2080. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  89. Welch, R.A.; Burland, V.; Plunkett, G.; Redford, P.; Roesch, P.; Rasko, D.; Buckles, E.L.; Liou, S.R.; Boutin, A.; Hackett, J.; et al. Extensive mosaic structure revealed by the complete genome sequence of uropathogenic Escherichia coli. Proc. Natl. Acad. Sci. USA 2002, 99, 17020–17024. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  90. Didelot, X.; Falush, D. Inference of bacterial microevolution using multilocus sequence data. Genetics 2007, 175, 1251–1266. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  91. Nguyen, M.; Ekstrom, A.; Li, X.; Yin, Y. HGT-Finder: A New Tool for Horizontal Gene Transfer Finding and Application to Aspergillus genomes. Toxins 2015, 7, 4035–4053. [Google Scholar] [CrossRef] [Green Version]
  92. Cui, Z.; Wang, S.; Kakar, K.U.; Xie, G.; Li, B.; Chen, G.; Zhu, B. Genome Sequence and Adaptation Analysis of the Human and Rice Pathogenic Strain Burkholderia glumae AU6208. Pathogens 2021, 10, 87. [Google Scholar] [CrossRef]
  93. Zhu, Q.; Kosoy, M.; Dittmar, K. HGTector: An automated method facilitating genome-wide discovery of putative horizontal gene transfers. BMC Genom. 2014, 15, 717. [Google Scholar] [CrossRef] [Green Version]
  94. Shimada, S.; Nakai, R.; Aoki, K.; Kudoh, S.; Imura, S.; Shimoeda, N.; Ohno, G.; Watanabe, K.; Miyazaki, Y.; Ishii, Y.; et al. Characterization of the First Cultured Psychrotolerant Representative of Legionella from Antarctica Reveals Its Unique Genome Structure. Microbiol. Spectr. 2021, 9, e0042421. [Google Scholar] [CrossRef]
  95. Xu, S.; Li, Z.; Huang, Y.; Han, L.; Che, Y.; Hou, X.; Li, D.; Fan, S.; Li, Z. Whole genome sequencing reveals the genomic diversity, taxonomic classification, and evolutionary relationships of the genus Nocardia. PLoS Negl. Trop. Dis. 2021, 15, e0009665. [Google Scholar] [CrossRef]
  96. Li, X.; Tong, W.; Wang, L.; Rahman, S.U.; Wei, G.; Tao, S. A Novel Strategy for Detecting Recent Horizontal Gene Transfer and Its Application to Rhizobium Strains. Front. Microbiol. 2018, 9, 973. [Google Scholar] [CrossRef] [Green Version]
  97. Tong, W.; Li, X.; Wang, E.; Cao, Y.; Chen, W.; Tao, S.; Wei, G. Genomic insight into the origins and evolution of symbiosis genes in Phaseolus vulgaris microsymbionts. BMC Genom. 2020, 21, 186. [Google Scholar] [CrossRef]
  98. Sánchez-Soto, D.; Agüero-Chapin, G.; Armijos-Jaramillo, V.; Perez-Castillo, Y.; Tejera, E.; Antunes, A.; Sánchez-Rodríguez, A. ShadowCaster: Compositional Methods under the Shadow of Phylogenetic Models to Detect Horizontal Gene Transfers in Prokaryotes. Genes 2020, 11, 756. [Google Scholar] [CrossRef]
  99. Wang, W.B.; Jiang, T.; Gardner, S. Detection of Homologous Recombination Events in Bacterial Genomes. PLoS ONE 2013, 8, e75230. [Google Scholar] [CrossRef]
  100. Zhang, F.; Hu, Z.; Wu, Z.; Lu, J.; Shi, Y.; Xu, J.; Wang, X.; Wang, J.; Zhang, F.; Wang, M.; et al. Reciprocal adaptation of rice and Xanthomonas oryzae pv. oryzae: Cross-species 2D GWAS reveals the underlying genetics. Plant Cell 2021, 33, 2538–2561. [Google Scholar] [CrossRef]
  101. Lin, M.; Kussell, E. Inferring bacterial recombination rates from large-scale sequencing datasets. Nat. Methods 2019, 16, 199–204. [Google Scholar] [CrossRef]
  102. Crits-Christoph, A.; Olm, M.R.; Diamond, S.; Bouma-Gregson, K.; Banfield, J.F. Soil bacterial populations are shaped by recombination and gene-specific selection across a grassland meadow. ISME J. 2020, 14, 1834–1846. [Google Scholar] [CrossRef] [Green Version]
  103. Cámara, P.G.; Levine, A.J.; Rabadán, R. Inference of Ancestral Recombination Graphs through Topological Data Analysis. PLoS Comput. Biol. 2016, 12, e1005071. [Google Scholar] [CrossRef] [Green Version]
  104. Koonin, E.V.; Makarova, K.S.; Aravind, L. Horizontal gene transfer in prokaryotes: Quantification and classification. Annu. Rev. Microbiol. 2001, 55, 709–742. [Google Scholar] [CrossRef]
  105. Zamani-Dahaj, S.A.; Okasha, M.; Kosakowski, J.; Higgs, P.G. Estimating the Frequency of Horizontal Gene Transfer Using Phylogenetic Models of Gene Gain and Loss. Mol. Biol. Evol. 2016, 33, 1843–1857. [Google Scholar] [CrossRef] [Green Version]
  106. Jeong, H.; Nasir, A. A Preliminary List of Horizontally Transferred Genes in Prokaryotes Determined by Tree Reconstruction and Reconciliation. Front. Genet. 2017, 8, 112. [Google Scholar] [CrossRef] [Green Version]
  107. Vogan, A.A.; Higgs, P.G. The advantages and disadvantages of horizontal gene transfer and the emergence of the first species. Biol. Direct 2011, 6, 1. [Google Scholar] [CrossRef] [Green Version]
  108. Narechania, A.; Baker, R.; DeSalle, R.; Mathema, B.; Kolokotronis, S.O.; Kreiswirth, B.; Planet, P.J. Clusterflock: A flocking algorithm for isolating congruent phylogenomic datasets. Gigascience 2016, 5, s13742-016-0152-3. [Google Scholar] [CrossRef] [Green Version]
  109. Domazet-Lošo, M.; Domazet-Lošo, T. gmos: Rapid Detection of Genome Mosaicism over Short Evolutionary Distances. PLoS ONE 2016, 11, e0166602. [Google Scholar] [CrossRef] [Green Version]
  110. Wan, Y.; Wick, R.R.; Zobel, J.; Ingle, D.J.; Inouye, M.; Holt, K.E. GeneMates: An R package for detecting horizontal gene co-transfer between bacteria using gene-gene associations controlled for population structure. BMC Genom. 2020, 21, 658. [Google Scholar] [CrossRef] [PubMed]
  111. Trappe, K.; Marschall, T.; Renard, B.Y. Detecting horizontal gene transfer by mapping sequencing reads across species boundaries. Bioinformatics 2016, 32, i595–i604. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  112. Bay, R.A.; Bielawski, J.P. Recombination detection under evolutionary scenarios relevant to functional divergence. J. Mol. Evol. 2011, 73, 273–286. [Google Scholar] [CrossRef] [PubMed]
  113. Bertrand, Y.J.K.; Johansson, M.; Norberg, P. Revisiting Recombination Signal in the Tick-Borne Encephalitis Virus: A Simulation Approach. PLoS ONE 2016, 11, e0164435. [Google Scholar] [CrossRef]
  114. Reed, F.A.; Tishkoff, S.A. Positive selection can create false hotspots of recombination. Genetics 2006, 172, 2011–2014. [Google Scholar] [CrossRef] [Green Version]
  115. Sun, S.; Evans, B.J.; Golding, G.B. “Patchy-tachy” leads to false positives for recombination. Mol. Biol. Evol. 2011, 28, 2549–2559. [Google Scholar] [CrossRef] [Green Version]
  116. Friedman, R.; Ely, B. Codon usage methods for horizontal gene transfer detection generate an abundance of false positive and false negative results. Curr. Microbiol. 2012, 65, 639–642. [Google Scholar] [CrossRef]
  117. Godfroid, M.; Dagan, T.; Kupczok, A. Recombination Signal in Mycobacterium tuberculosis Stems from Reference-guided Assemblies and Alignment Artefacts. Genome Biol. Evol. 2018, 10, 1920–1926. [Google Scholar] [CrossRef] [Green Version]
  118. Meyerhans, A.; Vartanian, J.P.; Wain-Hobson, S. DNA recombination during PCR. Nucleic Acids Res. 1990, 18, 1687–1691. [Google Scholar] [CrossRef]
  119. Zagordi, O.; Klein, R.; Däumer, M.; Beerenwinkel, N. Error correction of next-generation sequencing data and reliable estimation of HIV quasispecies. Nucleic Acids Res. 2010, 38, 7400–7409. [Google Scholar] [CrossRef]
  120. Posada, D. Evaluation of Methods for Detecting Recombination from DNA Sequences: Empirical Data. Mol. Biol. Evol. 2002, 19, 708–717. [Google Scholar] [CrossRef] [Green Version]
  121. Brown, T.; Didelot, X.; Wilson, D.J.; Maio, N. De SimBac: Simulation of whole bacterial genomes with homologous recombination. Microb. Genom. 2016, 2, e000044. [Google Scholar] [CrossRef] [Green Version]
  122. Sipola, A.; Marttinen, P.; Corander, J. Bacmeta: Simulator for genomic evolution in bacterial metapopulations. Bioinformatics 2018, 34, 2308–2310. [Google Scholar] [CrossRef]
  123. Bobay, L.M. CoreSimul: A forward-in-time simulator of genome evolution for prokaryotes modeling homologous recombination. BMC Bioinform. 2020, 21, 1–7. [Google Scholar] [CrossRef]
  124. Song, W.; Steensen, K.; Thomas, T. HgtSIM: A simulator for horizontal gene transfer (HGT) in microbial communities. PeerJ 2017, 5, e4015. [Google Scholar] [CrossRef] [Green Version]
  125. Cury, J.; Haller, B.C.; Achaz, G.; Jay, F. Simulation of bacterial populations with SLiM. Peer Community J. 2022, 2, e7. [Google Scholar] [CrossRef]
  126. Marttinen, P.; Hanage, W.P.; Croucher, N.J.; Connor, T.R.; Harris, S.R.; Bentley, S.D.; Corander, J. Detection of recombination events in bacterial genomes from large population samples. Nucleic Acids Res. 2012, 40, 1–12. [Google Scholar] [CrossRef]
  127. Lehtinen, S.; Lehtinen, S.; Chewapreecha, C.; Chewapreecha, C.; Chewapreecha, C.; Lees, J.; Hanage, W.P.; Lipsitch, M.; Croucher, N.J.; Bentley, S.D.; et al. Horizontal gene transfer rate is not the primary determinant of observed antibiotic resistance frequencies in Streptococcus pneumonia. Sci. Adv. 2020, 6, eaaz6137. [Google Scholar] [CrossRef]
Figure 1. A combined classification of methods for detecting homologous recombination and horizontal gene transfer depending on the genomic consequences of the events. HR—homologous recombination, HGT—horizontal gene transfer, ARGs—ancestral recombination graphs.
Figure 1. A combined classification of methods for detecting homologous recombination and horizontal gene transfer depending on the genomic consequences of the events. HR—homologous recombination, HGT—horizontal gene transfer, ARGs—ancestral recombination graphs.
Ijms 23 06257 g001
Table 1. Current bioinformatics tools for detecting homologous recombination and horizontal gene transfer in genetic data. The table summarizes tools’ properties in terms of algorithms applied, input files and output results, type of detected events, advantages, and limitations.
Table 1. Current bioinformatics tools for detecting homologous recombination and horizontal gene transfer in genetic data. The table summarizes tools’ properties in terms of algorithms applied, input files and output results, type of detected events, advantages, and limitations.
ToolApplied ApproachMethod’s ClassInputOutputDetected EventsAdvantagesLimitationsReferences
Homologous Recombination (HR) Identification
RDP4/RDP5Combination of phylogenetic and distance methodsPhylogenetic and distance-basedAlignments in FASTA formatRecombination events with phylogenetic relationships and breakpoints coordinates for chimeric sequences in tabular formatRecentRobustness and providing the information on the direction of exchangesInability to reveal distant events and high computational costs[66,67]
GubbinsRevealing increased substitution rate among ML-tree branchesPhylogeneticAlignments in FASTA formatCoordinates of recombination events tabular format and their visualization on the genome alignmentRecent and ancestralPrecise reconstruction of ancestral stateHigh computational costs and possible false-positive results when analyzing trees with short branches (theoretically)[44]
ClonalFrameMLMaximal likelihood-based clonal modelPhylogeneticAlignments in FASTA format and guiding treePhylogeny regarding recombination and visualization of events’ coordinates on the genome alignment in tabular formatRecent and ancestralComputational effectivenessUnderestimation of recombination rate in datasets with intensive recombination[73]
BratNextGenBayesian modelingSubstitution distributionAlignments in FASTA formatCoordinates of the events in tabular format and visualization of transmitted regions on the genome alignmentRecent and ancestralComputational effectivenessFalse-negative results in the case of mosaic sequences with multiple recombination events[126]
fastGEARHMM algorithms coupled with Bayesian clusteringSubstitution distributionAlignments in FASTA formatCoordinates of ancestral and recent recombination events in tabular formatRecent and ancestralComputational effectiveness, high sensitivity, and handling of missing dataMissing events between closely related species[10]
ptACRGenome-wise average SNP compatibility calculationCompatibilityGap-free alignments in PHYLIP formatGenomic coordinates of recombination events in tabular formatRecentHigh accuracy and robustness to false-positive resultsInability to process alignments with gaps and high false-negative rate when processing divergent sequences[40]
HREfinderGenome partitioning into SNP-flanked blocksSubstitution distributionGenomes in FASTA format, tree in Newick format, and SNP list in tabular formatList of sequences subjected to recombination in tabular formatRecentHigh accuracyHigh false-negative rate when processing divergent sequences[99]
mcorrBuilding correlation profile of synonymous substitutionParametricAlignments in XMFA or BAM formatsTables and figures depicting the average recombination rateThe total rate of recent/ancient eventsThe ability to process raw reads and metagenomic dataHas not been compared to conventional r/m rate calculating tools[101]
BacterMarkov chain Monte Carlo (MCMC)ARGAlignments in FASTA formatAncestral recombination graph (ARG) in Newick formatRecentImproved detection of the events in the case of poor phylogenetic signalDependence on predetermined parameters and high computational costs[35]
TARGetTopological data analysis (TDA)ARGAlignments in FASTA format without gaps or segregating sites denoted by 1 and 0 Ancestral recombination graph (ARG) in XML format and positions of reticulate eventsRecentComputational effectivenessInability to process alignments with gaps[103]
Horizontal Gene Transfer (HGT) Detection
ClusterflockSelf-organizing flock algorithmParametricSequences and a distance matrixClusters of sequences in tabular formatRecentApplicability to any distance metrics and resilience to missing dataHas not been compared to the existing tools[108]
gmosPairwise local alignments with subsequent regions overlappingParametricQuery and subject genomes in FASTA formatStructural variants in FASTA formatRecentComputational effectiveness and the ability to reveal both HR and HGTDepends heavily on the high similarity between transferred regions[109]
GeneMatesAssociation tests with the linear-mixed model accounting for population structureParametricGenome assemblies in FASTA format and raw reads in FASTQ formatThe linkage network of horizontally co-transferred alleles in tabular formatRecentResolving co-occurred HGT eventsReduced sensitivity due to the dependence on core SNPs[110]
ShadowCasterSupport vector machine-based hybrid approachImplicit phylogenetic and parametricA query genome and proteome and list of related proteomes in FASTA formatThe list of HGT candidates with corresponding likelihood calculations in tabular formatRecent and ancestralHigh sensitivity when reveling both recent and ancient events and reduced false-positive rateDoes not determine the directions of transfers and processes only a single genome[98]
nearHGTCalculating synteny index (SI) followed by constant relative mutability (CRM) measurementSynteny-based and parametricReference and putatively transferred sequences in FASTA formatChi-square-based p-value denoting the probability of HGTRecentHigh sensitivityNo ready-made application is available[49]
HGT-FinderSimilarity ratio evaluation for proteins according to BLAST hits and taxonomic distance calculation based on the NCBI Taxonomy annotationImplicit phylogeneticThe BLAST search result and the NCBI Taxonomy databaseTabular format file with the transfer index value for a proteinRecentDetecting mostly true eventsHigh reliance on the taxonomic nomenclature and low sensitivity[91]
HGTectorAnalyzing BLAST hit distribution patterns according to predefined evolutionary categoriesImplicit phylogeneticFASTA files of amino acid sequences for each analyzed genomeList of candidate HGT-derived genes with the respective silhouette scores in tabular formatRecentInsensitive to gene loss, rate variations, and database errorsHigh reliance on the taxonomic nomenclature and low sensitivity[93]
RecentHGTThe expectation-maximization algorithm based on the sequence-similarity distribution of orthologous genesImplicit phylogeneticTabular file with strains information and RAST-annotated GenBank filePutative HGT events in chromosomal and plasmid regions in tabular formatRecentReduced false-positive rate when processing conserved genesMissing events when analyzing divergent sequences[96]
DaisyMapping-based detection relying on short read pairs and coverage informationParametricReads from the analyzed organism and poposed acceptor and donor genomes in FASTA formatA variant call format (VCF) file reporting HGT candidates meeting the predefined threshold and tabular format file with all potential eventsRecentOutperforms reference genome-based approaches if short reads are availableRequires short reads only and explicit specifying recipient and donor genomes[11]
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Shikov, A.E.; Malovichko, Y.V.; Nizhnikov, A.A.; Antonets, K.S. Current Methods for Recombination Detection in Bacteria. Int. J. Mol. Sci. 2022, 23, 6257. https://doi.org/10.3390/ijms23116257

AMA Style

Shikov AE, Malovichko YV, Nizhnikov AA, Antonets KS. Current Methods for Recombination Detection in Bacteria. International Journal of Molecular Sciences. 2022; 23(11):6257. https://doi.org/10.3390/ijms23116257

Chicago/Turabian Style

Shikov, Anton E., Yury V. Malovichko, Anton A. Nizhnikov, and Kirill S. Antonets. 2022. "Current Methods for Recombination Detection in Bacteria" International Journal of Molecular Sciences 23, no. 11: 6257. https://doi.org/10.3390/ijms23116257

APA Style

Shikov, A. E., Malovichko, Y. V., Nizhnikov, A. A., & Antonets, K. S. (2022). Current Methods for Recombination Detection in Bacteria. International Journal of Molecular Sciences, 23(11), 6257. https://doi.org/10.3390/ijms23116257

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop