Keeping Pace with the Red Queen: Identifying the Genetic Basis of Susceptibility to Infectious Disease

Inferring strain-level mutational drivers of phage-bacteria interaction phenotypes arising during coevolutionary dynamics

The enormous diversity of bacteriophages and their bacterial hosts presents a significant challenge to predict which phages infect a focal set of bacteria. Infection is largely determined by complementary-and largely uncharacterized-genetics of adsorption, injection, cell take-over, and lysis. Here we present a machine learning approach to predict phage-bacteria interactions trained on genome sequences of and phenotypic interactions among 51 Escherichia coli strains and 45 phage λ strains that coevolved in laboratory conditions for 37 days. Leveraging multiple inference strategies and without a priori knowledge of driver mutations, this framework predicts both who infects whom and the quantitative levels of infections across a suite of 2,295 potential interactions. We found that the most effective approach inferred interaction phenotypes from independent contributions from phage and bacteria mutations, accurately predicting 86% of interactions while reducing the relative error in the estimated strength of the infection phenotype by 40%. Feature selection revealed key phage λ and Escherchia coli mutations that have a significant influence on the outcome of phage-bacteria interactions, corroborating sites previously known to affect phage λ infections, as well as identifying mutations in genes of unknown function not previously shown to influence bacterial resistance. The method's success in recapitulating strain-level infection outcomes arising during coevolutionary dynamics may also help inform generalized approaches for imputing genetic drivers of interaction phenotypes in complex communities of phage and bacteria.

Inferring strain-level mutational drivers of phage-bacteria interaction phenotypes arising during coevolutionary dynamics

The enormous diversity of bacteriophages and their bacterial hosts presents a significant challenge to predict which phages infect a focal set of bacteria. Infection is largely determined by complementary - and largely uncharacterized - genetics of adsorption, injection, cell take-over and lysis. Here we present a machine learning approach to predict phage-bacteria interactions trained on genome sequences of and phenotypic interactions amongst 51 Escherichia coli strains and 45 phage λ strains that coevolved in laboratory conditions for 37 days. Leveraging multiple inference strategies and without a priori knowledge of driver mutations, this framework predicts both who infects whom and the quantitative levels of infections across a suite of 2,295 potential interactions. We found that the most effective approach inferred interaction phenotypes from independent contributions from phage and bacteria mutations, accurately predicting 86 % of interactions while reducing the relative error in the estimated strength of the infection phenotype by 40 % . Feature selection revealed key phage λ and E. coli mutations that have a significant influence on the outcome of phage-bacteria interactions, corroborating sites previously known to affect phage λ infections, as well as identifying mutations in genes of unknown function not previously shown to influence bacterial resistance. The method's success in recapitulating strain-level infection outcomes arising during coevolutionary dynamics may also help inform generalized approaches for imputing genetic drivers of interaction phenotypes in complex communities of phage and bacteria.

Inference of Host-Pathogen Interaction Matrices from Genome-Wide Polymorphism Data

Host-pathogen coevolution is defined as the reciprocal evolutionary changes in both species due to genotype × genotype (G×G) interactions at the genetic level determining the outcome and severity of infection. While co-analyses of hosts and pathogen genomes (co-genome-wide association studies) allow us to pinpoint the interacting genes, these do not reveal which host genotype(s) is/are resistant to which pathogen genotype(s). The knowledge of this so-called infection matrix is important for agriculture and medicine. Building on established theories of host-pathogen interactions, we here derive four novel indices capturing the characteristics of the infection matrix. These indices can be computed from full genome polymorphism data of randomly sampled uninfected hosts, as well as infected hosts and their pathogen strains. We use these indices in an approximate Bayesian computation method to pinpoint loci with relevant G×G interactions and to infer their underlying interaction matrix. In a combined single nucleotide polymorphism dataset of 451 European humans and their infecting hepatitis C virus (HCV) strains and 503 uninfected individuals, we reveal a new human candidate gene for resistance to HCV and new virus mutations matching human genes. For two groups of significant human-HCV (G×G) associations, we infer a gene-for-gene infection matrix, which is commonly assumed to be typical of plant-pathogen interactions. Our model-based inference framework bridges theoretical models of G×G interactions with host and pathogen genomic data. It, therefore, paves the way for understanding the evolution of key G×G interactions underpinning HCV adaptation to the European human population after a recent expansion.

Intergenomic signatures of coevolution between Tasmanian devils and an infectious cancer

Coevolution is common and frequently governs host-pathogen interaction outcomes. Phenotypes underlying these interactions often manifest as the combined products of the genomes of interacting species, yet traditional quantitative trait mapping approaches ignore these intergenomic interactions. Devil facial tumor disease (DFTD), an infectious cancer afflicting Tasmanian devils (Sarcophilus harrisii), has decimated devil populations due to universal host susceptibility and a fatality rate approaching 100%. Here, we used a recently developed joint genome-wide association study (i.e., co-GWAS) approach, 15 y of mark-recapture data, and 960 genomes to identify intergenomic signatures of coevolution between devils and DFTD. Using a traditional GWA approach, we found that both devil and DFTD genomes explained a substantial proportion of variance in how quickly susceptible devils became infected, although genomic architectures differed across devils and DFTD; the devil genome had fewer loci of large effect whereas the DFTD genome had a more polygenic architecture. Using a co-GWA approach, devil-DFTD intergenomic interactions explained ~3× more variation in how quickly susceptible devils became infected than either genome alone, and the top genotype-by-genotype interactions were significantly enriched for cancer genes and signatures of selection. A devil regulatory mutation was associated with differential expression of a candidate cancer gene and showed putative allele matching effects with two DFTD coding sequence variants. Our results highlight the need to account for intergenomic interactions when investigating host-pathogen (co)evolution and emphasize the importance of such interactions when considering devil management strategies.

Integrating plant and fungal quantitative genetics to improve the ecological and agricultural applications of mycorrhizal symbioses

Finding and targeting genes that quantitatively contribute to agricultural and ecological processes progresses food production and conservation efforts. Typically, quantitative genetic approaches link variants in a single organism's genome with a trait of interest. Recently, genome-to-genome mapping has found genome variants interacting between species to produce the result of a multiorganism (including multikingdom) interaction. These were plant and bacterial pathogen genome interactions; plant-fungal coquantitative genetics have not yet been applied. Plant-mycorrhizae symbioses exist across most biomes, for a majority of land plants, including crop plants, and manipulate many traits from single organisms to ecosystems for which knowing the genetic basis would be useful. The availability of Rhizophagus irregularis mycorrhizal isolates, with genomic information, makes dual-genome methods with beneficial mutualists accessible and imminent.

Combining GWAS and population genomic analyses to characterize coevolution in a legume-rhizobia symbiosis

The mutualism between legumes and rhizobia is clearly the product of past coevolution. However, the nature of ongoing evolution between these partners is less clear. To characterize the nature of recent coevolution between legumes and rhizobia, we used population genomic analysis to characterize selection on functionally annotated symbiosis genes as well as on symbiosis gene candidates identified through a two-species association analysis. For the association analysis, we inoculated each of 202 accessions of the legume host Medicago truncatula with a community of 88 Sinorhizobia (Ensifer) meliloti strains. Multistrain inoculation, which better reflects the ecological reality of rhizobial selection in nature than single-strain inoculation, allows strains to compete for nodulation opportunities and host resources and for hosts to preferentially form nodules and provide resources to some strains. We found extensive host by symbiont, that is, genotype-by-genotype, effects on rhizobial fitness and some annotated rhizobial genes bear signatures of recent positive selection. However, neither genes responsible for this variation nor annotated host symbiosis genes are enriched for signatures of either positive or balancing selection. This result suggests that stabilizing selection dominates selection acting on symbiotic traits and that variation in these traits is under mutation-selection balance. Consistent with the lack of positive selection acting on host genes, we found that among-host variation in growth was similar whether plants were grown with rhizobia or N-fertilizer, suggesting that the symbiosis may not be a major driver of variation in plant growth in multistrain contexts.

Unraveling coevolutionary dynamics using ecological genomics

Coevolutionary interactions, from the delicate co-dependency in mutualistic interactions to the antagonistic relationship of hosts and parasites, are a ubiquitous driver of adaptation. Surprisingly, little is known about the genomic processes underlying coevolution in an ecological context. However, species comprise genetically differentiated populations that interact with temporally variable abiotic and biotic environments. We discuss the recent advances in coevolutionary theory and genomics as well as shortcomings, to identify coevolving genes that take into account this spatial and temporal variability of coevolution, and propose a practical guide to understand the dynamic of coevolution using an ecological genomics lens.

Negotiating mutualism: A locus for exploitation by rhizobia has a broad effect size distribution and context-dependent effects on legume hosts

In mutualisms, variation at genes determining partner fitness provides the raw material upon which coevolutionary selection acts, setting the dynamics and pace of coevolution. However, we know little about variation in the effects of genes that underlie symbiotic fitness in natural mutualist populations. In some species of legumes that form root nodule symbioses with nitrogen‐fixing rhizobial bacteria, hosts secrete nodule‐specific cysteine‐rich (NCR) peptides that cause rhizobia to differentiate in the nodule environment. However, rhizobia can cleave NCR peptides through the expression of genes like the plasmid‐borne Host range restriction peptidase (hrrP), whose product degrades specific NCR peptides. Although hrrP activity can confer host exploitation by depressing host fitness and enhancing symbiont fitness, the effects of hrrP on symbiosis phenotypes depend strongly on the genotypes of the interacting partners. However, the effects of hrrP have yet to be characterised in a natural population context, so its contribution to variation in wild mutualist populations is unknown. To understand the distribution of effects of hrrP in wild rhizobia, we measured mutualism phenotypes conferred by hrrP in 12 wild Ensifer medicae strains. To evaluate context dependency of hrrP effects, we compared hrrP effects across two Medicago polymorpha host genotypes and across two experimental years for five E. medicae strains. We show for the first time in a natural population context that hrrP has a wide distribution of effect sizes for many mutualism traits, ranging from strongly positive to strongly negative. Furthermore, we show that hrrP effect size varies across host genotypes and experiment years, suggesting that researchers should be cautious about extrapolating the role of genes in natural populations from controlled laboratory studies of single genetic variants.

Coevolutionary theory of hosts and parasites

Host and parasite evolution are closely intertwined, with selection for adaptations and counter-adaptations forming a coevolutionary feedback loop. Coevolutionary dynamics are often difficult to intuit due to these feedbacks and are hard to demonstrate empirically in most systems. Theoretical models have therefore played a crucial role in shaping our understanding of host-parasite coevolution. Theoretical models vary widely in their assumptions, approaches and aims, and such variety makes it difficult, especially for non-theoreticians and those new to the field, to: (1) understand how model approaches relate to one another; (2) identify key modelling assumptions; (3) determine how model assumptions relate to biological systems; and (4) reconcile the results of different models with contrasting assumptions. In this review, we identify important model features, highlight key results and predictions and describe how these pertain to model assumptions. We carry out a literature survey of theoretical studies published since the 1950s (n = 219 papers) to support our analysis. We identify two particularly important features of models that tend to have a significant qualitative impact on the outcome of host-parasite coevolution: population dynamics and the genetic basis of infection. We also highlight the importance of other modelling features, such as stochasticity and whether time proceeds continuously or in discrete steps, that have received less attention but can drastically alter coevolutionary dynamics. We finish by summarizing recent developments in the field, specifically the trend towards greater model complexity, and discuss likely future directions for research.

Challenges and Opportunities in Understanding Genetics of Fungal Diseases: Towards a Functional Genomics Approach

Infectious diseases are a leading cause of morbidity and mortality worldwide, and human pathogens have long been recognized as one of the main sources of evolutionary pressure, resulting in a high variable genetic background in immune-related genes. The study of the genetic contribution to infectious diseases has undergone tremendous advances over the last decades. Here, focusing on genetic predisposition to fungal diseases, we provide an overview of the available approaches for studying human genetic susceptibility to infections, reviewing current methodological and practical limitations. We describe how the classical methods available, such as family-based studies and candidate gene studies, have contributed to the discovery of crucial susceptibility factors for fungal infections. We will also discuss the contribution of novel unbiased approaches to the field, highlighting their success but also their limitations for the fungal immunology field. Finally, we show how a systems genomics approach can overcome those limitations and can lead to efficient prioritization and identification of genes and pathways with a critical role in susceptibility to fungal diseases. This knowledge will help to stratify at-risk patient groups and, subsequently, develop early appropriate prophylactic and treatment strategies.

Novel genomic approaches to study antagonistic coevolution between hosts and parasites

Host-parasite coevolution is ubiquitous, shaping genetic and phenotypic diversity and the evolutionary trajectory of interacting species. With the advances of high throughput sequencing technologies applicable to model and non-model organisms alike, it is now feasible to study in greater detail (a) the genetic underpinnings of coevolution, (b) the speed and type of dynamics at coevolving loci, and (c) the genomic consequences of coevolution. This review focuses on three recently developed approaches that leverage information from host and parasite full genome data simultaneously to pinpoint coevolving loci and draw inference on the coevolutionary history. First, co-genome-wide association study (co-GWAS) methods allow pinpointing the loci underlying host-parasite interactions. These methods focus on detecting associations between genetic variants and the outcome of experimental infection tests or on correlations between genomes of naturally infected hosts and their infecting parasites. Second, extensions to population genomics methods can detect genes under coevolution and infer the coevolutionary history, such as fitness costs. Third, correlations between host and parasite population size in time are indicative of coevolution, and polymorphism levels across independent spatially distributed populations of hosts and parasites can reveal coevolutionary loci and infer coevolutionary history. We describe the principles of these three approaches and discuss their advantages and limitations based on coevolutionary theory. We present recommendations for their application to various host (prokaryotes, fungi, plants, and animals) and parasite (viruses, bacteria, fungi, and macroparasites) species. We conclude by pointing out methodological and theoretical gaps to be filled to extract maximum information from full genome data and thereby to shed light on the molecular underpinnings of coevolution.

Feedback between coevolution and epidemiology can help or hinder the maintenance of genetic variation in host-parasite models

Antagonistic coevolution has long been suggested to help maintain host genetic variation. Although ecological and epidemiological feedbacks are known to have important consequences on coevolutionary allele-frequency dynamics, their effects on the maintenance of genetic variation remains poorly understood. Here, we extend previous work on the maintenance of genetic variation in a classic matching alleles coevolutionary model by exploring the effects of ecological and epidemiological feedbacks, where both allele frequencies and population sizes are allowed to vary over time. We find that coevolution rarely maintains more host genetic variation than expected under neutral genetic drift alone. When and if coevolution maintains or depletes genetic variation relative to neutral drift is determined, predominantly, by two factors: the deterministic stability of the Red Queen allele-frequency cycles and the chance of allele fixation in the pathogen, as this results in directional selection and depletion of genetic variation in the host. Compared to purely coevolutionary models with constant host and pathogen population sizes, ecological and epidemiological feedbacks stabilize Red Queen cycles deterministically, but population fluctuations in the pathogen increase the rate of allele fixation in the pathogen, especially in epidemiological models. Our results illustrate the importance of considering the ecological and epidemiological context in which coevolution occurs when examining the impact of Red Queen cycles on genetic variation.

Tackling microbial threats in agriculture with integrative imaging and computational approaches

Pathogens and pests are one of the major threats to agricultural productivity worldwide. For decades, targeted resistance breeding was used to create crop cultivars that resist pathogens and environmental stress while retaining yields. The often decade-long process of crossing, selection, and field trials to create a new cultivar is challenged by the rapid rise of pathogens overcoming resistance. Similarly, antimicrobial compounds can rapidly lose efficacy due to resistance evolution. Here, we review three major areas where computational, imaging and experimental approaches are revolutionizing the management of pathogen damage on crops. Recognizing and scoring plant diseases have dramatically improved through high-throughput imaging techniques applicable both under well-controlled greenhouse conditions and directly in the field. However, computer vision of complex disease phenotypes will require significant improvements. In parallel, experimental setups similar to high-throughput drug discovery screens make it possible to screen thousands of pathogen strains for variation in resistance and other relevant phenotypic traits. Confocal microscopy and fluorescence can capture rich phenotypic information across pathogen genotypes. Through genome-wide association mapping approaches, phenotypic data helps to unravel the genetic architecture of stress- and virulence-related traits accelerating resistance breeding. Finally, joint, large-scale screenings of trait variation in crops and pathogens can yield fundamental insights into how pathogens face trade-offs in the adaptation to resistant crop varieties. We discuss how future implementations of such innovative approaches in breeding and pathogen screening can lead to more durable disease control.

A need to consider the evolutionary genetics of host-symbiont mutualisms

Despite the ubiquity and importance of mutualistic interactions, we know little about the evolutionary genetics underlying their long-term persistence. As in antagonistic interactions, mutualistic symbioses are characterized by substantial levels of phenotypic and genetic diversity. In contrast to antagonistic interactions, however, we, by and large, do not understand how this variation arises, how it is maintained, nor its implications for future evolutionary change. Currently, we rely on phenotypic models to address the persistence of mutualistic symbioses, but the success of an interaction almost certainly depends heavily on genetic interactions. In this review, we argue that evolutionary genetic models could provide a framework for understanding the causes and consequences of diversity and why selection may favour processes that maintain variation in mutualistic interactions.

Detecting HLA-infectious disease associations for multi-strain pathogens

Human Leukocyte Antigen (HLA) molecules play a vital role helping our immune system to detect the presence of pathogens. Previous work to try and ascertain which HLA alleles offer advantages against particular pathogens has generated inconsistent results. We have constructed an epidemiological model to understand why this may occur. The model captures the epidemiology of a multi strain pathogen for which the host's ability to generate immunological memory responses to particular strains depends on that host's HLA genotype. We find that an HLA allele's ability to protect against infection, as measured in a case control study, depends on the population frequency of that HLA allele. Furthermore, our capability to detect associations between HLA alleles and infection with a multi strain pathogen may be affected by the properties of the pathogen itself (i.e R₀ and length of infectious period). Both host and pathogen genetics must be considered in order to identify true HLA associations. However, in the absence of detailed pathogen genetic information, a negative correlation between the frequency of an HLA type and its apparent protectiveness against disease caused by multi strain pathogen is a strong indication that the HLA type in question is well adapted to a subset of strains of that pathogen.

Inference of coevolutionary dynamics and parameters from host and parasite polymorphism data of repeated experiments

There is a long-standing interest in understanding host-parasite coevolutionary dynamics and associated fitness effects. Increasing amounts of genomic data for both interacting species offer a promising source to identify candidate loci and to infer the main parameters of the past coevolutionary history. However, so far no method exists to perform the latter. By coupling a gene-for-gene model with coalescent simulations, we first show that three types of biological costs, namely, resistance, infectivity and infection, define the allele frequencies at the internal equilibrium point of the coevolution model. These in return determine the strength of selective signatures at the coevolving host and parasite loci. We apply an Approximate Bayesian Computation (ABC) approach on simulated datasets to infer these costs by jointly integrating host and parasite polymorphism data at the coevolving loci. To control for the effect of genetic drift on coevolutionary dynamics, we assume that 10 or 30 repetitions are available from controlled experiments or several natural populations. We study two scenarios: 1) the cost of infection and population sizes (host and parasite) are unknown while costs of infectivity and resistance are known, and 2) all three costs are unknown while populations sizes are known. Using the ABC model choice procedure, we show that for both scenarios, we can distinguish with high accuracy pairs of coevolving host and parasite loci from pairs of neutrally evolving loci, though the statistical power decreases with higher cost of infection. The accuracy of parameter inference is high under both scenarios especially when using both host and parasite data because parasite polymorphism data do inform on costs applying to the host and vice-versa. As the false positive rate to detect pairs of genes under coevolution is small, we suggest that our method complements recently developed methods to identify host and parasite candidate loci for functional studies.

Genomics of sorghum local adaptation to a parasitic plant

Host-parasite coevolution can maintain high levels of genetic diversity in traits involved in species interactions. In many systems, host traits exploited by parasites are constrained by use in other functions, leading to complex selective pressures across space and time. Here, we study genome-wide variation in the staple crop Sorghum bicolor (L.) Moench and its association with the parasitic weed Striga hermonthica (Delile) Benth., a major constraint to food security in Africa. We hypothesize that geographic selection mosaics across gradients of parasite occurrence maintain genetic diversity in sorghum landrace resistance. Suggesting a role in local adaptation to parasite pressure, multiple independent loss-of-function alleles at sorghum LOW GERMINATION STIMULANT 1 (LGS1) are broadly distributed among African landraces and geographically associated with S. hermonthica occurrence. However, low frequency of these alleles within S. hermonthica-prone regions and their absence elsewhere implicate potential trade-offs restricting their fixation. LGS1 is thought to cause resistance by changing stereochemistry of strigolactones, hormones that control plant architecture and below-ground signaling to mycorrhizae and are required to stimulate parasite germination. Consistent with trade-offs, we find signatures of balancing selection surrounding LGS1 and other candidates from analysis of genome-wide associations with parasite distribution. Experiments with CRISPR-Cas9-edited sorghum further indicate that the benefit of LGS1-mediated resistance strongly depends on parasite genotype and abiotic environment and comes at the cost of reduced photosystem gene expression. Our study demonstrates long-term maintenance of diversity in host resistance genes across smallholder agroecosystems, providing a valuable comparison to both industrial farming systems and natural communities.

Dissecting the genetic architecture of a stepwise infection process

How a host fights infection depends on an ordered sequence of steps, beginning with attempts to prevent a pathogen from establishing an infection, through to steps that mitigate a pathogen's control of host resources or minimize the damage caused during infection. Yet empirically characterizing the genetic basis of these steps remains challenging. Although each step is likely to have a unique genetic and environmental signature, and may therefore respond to selection in different ways, events that occur earlier in the infection process can mask or overwhelm the contributions of subsequent steps. In this study, we dissect the genetic architecture of a stepwise infection process using a quantitative trait locus (QTL) mapping approach. We control for variation at the first line of defence against a bacterial pathogen and expose downstream genetic variability related to the host's ability to mitigate the damage pathogens cause. In our model, the water-flea Daphnia magna, we found a single major effect QTL, explaining 64% of the variance, that is linked to the host's ability to completely block pathogen entry by preventing their attachment to the host oesophagus; this is consistent with the detection of this locus in previous studies. In susceptible hosts allowing attachment, however, a further 23 QTLs, explaining between 5% and 16% of the variance, were mapped to traits related to the expression of disease. The general lack of pleiotropy and epistasis for traits related to the different stages of the infection process, together with the wide distribution of QTLs across the genome, highlights the modular nature of a host's defence portfolio, and the potential for each different step to evolve independently. We discuss how isolating the genetic basis of individual steps can help to resolve discussion over the genetic architecture of host resistance.

Ecological and Evolutionary Processes Shaping Viral Genetic Diversity

The contemporary genomic diversity of viruses is a result of the continuous and dynamic interaction of past ecological and evolutionary processes. Thus, genome sequences of viruses can be a valuable source of information about these processes. In this review, we first describe the relevant processes shaping viral genomic variation, with a focus on the role of host–virus coevolution and its potential to give rise to eco-evolutionary feedback loops. We further give a brief overview of available methodology designed to extract information about these processes from genomic data. Short generation times and small genomes make viruses ideal model systems to study the joint effect of complex coevolutionary and eco-evolutionary interactions on genetic evolution. This complexity, together with the diverse array of lifetime and reproductive strategies in viruses ask for extensions of existing inference methods, for example by integrating multiple information sources. Such integration can broaden the applicability of genetic inference methods and thus further improve our understanding of the role viruses play in biological communities.

Open questions: what are the genes underlying antagonistic coevolution?

Although the idea of coevolution was first presented 150 years ago, we still only vaguely understand the genetic basis of its workings. Identifying the genes responsible for coevolutionary interactions would enable us to distinguish between fundamentally different models of coevolution and would represent a milestone in population genetics and genomics.