Epistasis, core-genome disharmony, and adaptation in recombining bacteria

Aidan J. Taylor,Koji Yahara,Ben Pascoe,Seungwon Ko,Leonardos Mageiros,Evangelos Mourkas,Jessica K. Calland,Santeri Puranen,Matthew D. Hitchings,Keith A. Jolley,Carolin M. Kobras,Sion Bayliss,Nicola J. Williams,Arnoud H. M. van Vliet,Julian Parkhill,Martin C. J. Maiden,Jukka Corander,Laurence D. Hurst,Daniel Falush,Paul Keim,Xavier Didelot,David J. Kelly,Samuel K. Sheppard
DOI: https://doi.org/10.1128/mbio.00581-24
IF: 6.4
2024-04-30
mBio
Abstract:ABSTRACT Recombination of short DNA fragments via horizontal gene transfer (HGT) can introduce beneficial alleles, create genomic disharmony through negative epistasis, and create adaptive gene combinations through positive epistasis. For non-core (accessory) genes, the negative epistatic cost is likely to be minimal because the incoming genes have not co-evolved with the recipient genome and are frequently observed as tightly linked cassettes with major effects. By contrast, interspecific recombination in the core genome is expected to be rare because disruptive allelic replacement is likely to introduce negative epistasis. Why then is homologous recombination common in the core of bacterial genomes? To understand this enigma, we take advantage of an exceptional model system, the common enteric pathogens Campylobacter jejuni and C. coli that are known for very high magnitude interspecies gene flow in the core genome. As expected, HGT does indeed disrupt co-adapted allele pairings, indirect evidence of negative epistasis. However, multiple HGT events enable recovery of the genome’s co-adaption between introgressing alleles, even in core metabolism genes (e.g., formate dehydrogenase). These findings demonstrate that, even for complex traits, genetic coalitions can be decoupled, transferred, and independently reinstated in a new genetic background—facilitating transition between fitness peaks. In this example, the two-step recombinational process is associated with C. coli that are adapted to the agricultural niche. IMPORTANCE Genetic exchange among bacteria shapes the microbial world. From the acquisition of antimicrobial resistance genes to fundamental questions about the nature of bacterial species, this powerful evolutionary force has preoccupied scientists for decades. However, the mixing of genes between species rests on a paradox: 0n one hand, promoting adaptation by conferring novel functionality; on the other, potentially introducing disharmonious gene combinations (negative epistasis) that will be selected against. Taking an interdisciplinary approach to analyze natural populations of the enteric bacteria Campylobacter , an ideal example of long-range admixture, we demonstrate that genes can independently transfer across species boundaries and rejoin in functional networks in a recipient genome. The positive impact of two-gene interactions appears to be adaptive by expanding metabolic capacity and facilitating niche shifts through interspecific hybridization. This challenges conventional ideas and highlights the possibility of multiple-step evolution of multi-gene traits by interspecific introgression.
microbiology
What problem does this paper attempt to address?
The core problem that this paper attempts to solve is: **How does horizontal gene transfer (HGT) introduce beneficial allele combinations in the bacterial core genome while avoiding genomic disharmony caused by negative epistasis?** Specifically, the author studied why homologous recombination is common in the bacterial core genome, although this recombination may disrupt highly adapted allele pairings. ### Problem Background 1. **The Role of Horizontal Gene Transfer (HGT)** - HGT can promote adaptive evolution by introducing beneficial alleles. - However, HGT may also introduce genomic disharmony through negative epistasis, especially when foreign genes have not co - evolved with the recipient genome. 2. **Differences between the Core Genome and the Non - core Genome** - For non - core genes (accessory genes), the cost of negative epistasis may be small because these genes usually exist in a tightly linked manner and have significant effects. - In contrast, inter - species recombination in the core genome is expected to be rare because disruptive allele replacements may introduce negative epistasis. 3. **Paradoxical Phenomenon** - Despite the above expectations, in actual observations, homologous recombination is common in the bacterial core genome. This has aroused the interest of researchers to try to explain this phenomenon. ### Research Objects The author selected **two important pathogenic bacteria in the genus Campylobacter - Campylobacter jejuni and Campylobacter coli** as research models. There is a large amount of inter - species gene flow between these two species, especially in the core genome. Specifically: - Some lineages of C. coli (such as CC - 828 and CC - 1150) contain a large number of gene fragments from C. jejuni, which provides an ideal model system for studying the impact of HGT on the core genome. - These lineages have adapted to the agricultural environment, indicating that HGT may have played an important role in their adaptation process. ### Main Findings 1. **HGT does disrupt co - adapted allele pairings**, which is indirect evidence of negative epistasis. 2. **Multiple HGT events can restore the co - adaptability of the genome**, even in core metabolic genes (such as formate dehydrogenase). 3. **Two - step Recombination Process**: The first recombination introduces negative selection pressure, but the subsequent second recombination restores the functional link, allowing the gene combination to re - establish harmony in a new genetic background. ### Conclusions This study shows that even for complex traits, gene alliances can be decoupled, transferred, and independently reconstructed in a new genetic background. This mechanism promotes adaptive evolution, especially in cases across species boundaries. In addition, this study challenges the traditional view of the evolution of polygenic traits and emphasizes the possibility of multi - step evolution through inter - species introgression. ### Key Formulas To describe the phenotypic interactions between genes, a continuous - time Markov chain (CTMC) model was used in the article to estimate the joint evolutionary parameters between biallelic loci. The specific formula is as follows: \[ P_{ij}(t)=\begin{pmatrix} \alpha_1 e^{-(\alpha_1+\beta_1)t}&(1 - e^{-(\alpha_1+\beta_1)t})\\ (1 - e^{-(\alpha_2+\beta_2)t})&\alpha_2 e^{-(\alpha_2+\beta_2)t} \end{pmatrix} \] where \(P_{ij}(t)\) represents the probability of transitioning from state \(i\) to state \(j\) within time \(t\), and \(\alpha_1,\alpha_2,\beta_1,\beta_2\) represent the independent substitution rates at different loci respectively. In this way, researchers can quantify and understand the phenotypic interactions between genes and their impact on adaptive evolution.