Abstract
A major goal in infectious disease research is to identify the human and pathogenic genetic variants that explain differences in microbial pathogenesis. However, neither pathogenic strain nor human genetic variation in isolation has proven adequate to explain the heterogeneity of disease pathology. We suggest that disrupted co-evolution between a pathogen and its human host can explain variation in disease outcomes, and that genome-by-genome interactions should therefore be incorporated into genetic models of disease caused by infectious agents. Genetic epidemiological studies that fail to take both the pathogen and host into account can lead to false and misleading conclusions about disease etiology. We discuss our model in the context of three pathogens, Helicobacter pylori, Mycobacterium tuberculosis and human papillomavirus, and generalize the conditions under which it may be applicable.
INTRODUCTION
Human response to infectious agents is known to be highly heritable, but identifying the genetic variants responsible for differences in disease susceptibility has proven difficult. Pathogenic variation has, in some cases, become a better predictor of disease outcome, but it too does not sufficiently predict whether a given individual or class of individuals will present with disease. Thus far, genetic epidemiological studies of infectious disease have typically sought to explain the inter-individual variation in disease phenotypes by assessing genetic factors in humans or pathogens alone, under the implicit assumption that these factors have effects that are essentially independent of each other. Here, we argue that genome-by-genome interactions between host and pathogen are likely to play a major role in infectious disease etiology, and as such, should be incorporated into genetic epidemiological models. In short, insofar as host and pathogen jointly determine disease phenotypes, no genetic variant in either should be considered harmful without taking the context of the other into account.
The term “interaction” has two related but distinct meanings in the context of infectious disease, one molecular, and one statistical. Here we refer mainly to the statistical meaning of the term. At the individual level, all aspects of pathogenesis involve molecular interactions of varying importance, e.g., between a pathogenic epitope and a host receptor. Such interactions can be detected statistically, however, only when multiple variants exist in a population and when specific pairings lead to different effects. In some cases, pathogenic variants may function independently of host variation, and vice versa. However, because many pathogens have co-existed with their human hosts for millennia and have likely co-evolved with them, we argue here that statistical interactions, where appropriately sought, will often be found, with profound biomedical implications.
Recent advances in genomics have provided both the impetus and the means to evaluate human–pathogen co-evolutionary hypotheses directly. Whole-genome sequencing of many pathogenic species has substantially improved the resolution with which we classify strains, and facilitated the detection of potentially virulent genetic variants. A clearer picture of microbial evolution has also emerged, marked by selective mechanisms such as rapid gene gain/loss and horizontal gene transfer (Pallen and Wren, 2007). Overlaying human genetic variation onto this emerging evolutionary picture of microbial diversity offers the potential to make the pathogenic process more transparent.
The past few decades have also seen an explosion in studies seeking to identify human susceptibility loci for infectious diseases (Rowell et al., 2012). Candidate gene and family based linkage studies have identified several common polymorphisms with clinical significance at the population level, such as the CCR5 deletion that protects against HIV (Samson et al., 1996; Picard et al., 2006; ). However, most human susceptibility is in fact polygenic, with individual polymorphisms conferring small marginal effects (Hill, 2001). Where infectious disease phenotypes deviate from the “one susceptibility locus – one infection” model, elucidating the genetic architecture underlying inter-individual variation has proven elusive.
While genome-wide association studies (GWAS) may be better designed to accommodate multifactorial phenotypes, those performed thus far on infectious diseases have typically been less informative than GWAS performed on complex non-communicable diseases (Jallow et al., 2009; Hill, 2012; Ko and Urban, 2013). A major challenge facing the GWAS of infectious disease has been the recruitment of a sufficient number of cases and matched controls to achieve adequate statistical power (Hill, 2012; Ko and Urban, 2013). Another potential drawback, and the one that concerns us here, is the fact that many infectious disease phenotypes depend on complex interactions between host and pathogen genomes. In such cases, the pooling together of human samples infected with even subtly different pathogenic strains can obscure genetic associations (Hill, 2012; Ko and Urban, 2013). A problem common to all GWAS is that the statistical effect sizes of biologically meaningful polymorphisms are often too small to pass significance thresholds after correction for multiple testing. This problem is exacerbated, however, when human polymorphisms (or networks of polymorphisms) (Wilfert and Schmid-Hempel, 2008) confer variable, or even opposite effects in the context of different pathogenic strains within the same study cohort. In this regard, it is perhaps telling that the most successful GWAS performed on infectious disease susceptibility to date have been on leprosy; the signal-to-noise ratios in these association studies may be higher because Mycobacterium leprae exhibits substantially less genetic heterogeneity than many other pathogens (Monot et al., 2009; Hill, 2012).
There is in fact strong empirical and theoretical justification for the hypothesis that the effects of susceptibility and virulence alleles in the respective gene pools of humans and pathogens are often contingent upon each other. The evolution of virulence is a dynamic process, easily perturbed by extrinsic variables over space and time, and therefore unlikely to follow the same trajectory in every population. For example, a spike in the density of hosts available for transmission can select for increased virulence, by reducing the cost of lethal harm (). If a pathogen is transmitted vertically (parent to child), the genetic factors that affect pathogenicity are “co-inherited” by host and pathogen, often promoting commensalism (Frank, 1996; Messenger et al., 1999). Even in these cases, the adventitious introduction of a microbial competitor can induce a commensal species to evolve a defensive toxin that harms the host, if only incidentally (; Frank and Schmid-Hempel, 2008). The evolution of defenses against pathogenic harm must also navigate fitness tradeoffs that vary with population, including tradeoffs pertaining to the correlated nature of complex traits (Lambrechts et al., 2006). As pathogens evolve rapidly, exerting strong selective pressures on different human populations, host phenotypes will respond in the ad hoc manner typical of evolution, limited by the available genetic variation at hand (Jacob, 1977). Whether the result is a steady-state equilibrium due to a perpetual “arms race” or a commensal detente, the same genes and pathways are unlikely to be involved in every population. As a consequence, when humans and pathogens migrate to new environments or admix, the ensuing disruption of co-evolutionary equilibria and loss of complementarity between host and pathogen genotypes may yield unpredictable and potentially deleterious biomedical consequences.
Our emphasis on the significance of mismatched traits is consistent with the genetic mosaic theory of co-evolution, which aims to account for why virtually all co-evolutionary interactions observed in natural populations show spatial variation in outcomes (Thompson et al., 2002; Thompson, 2014). The theory posits that co-evolution occurs in the context of geographically distinct “selection mosaics,” each characterized by a unique genetic and environmental profile, where environmental variables can include both biotic and abiotic factors. Every selection mosaic progresses toward its own co-evolutionary equilibrium, while gene flow between selection mosaics ensures that patterns of maladaptation will be common and detectable where properly studied (Thompson et al., 2002; Ridenhour and Nuismer, 2007).
Despite the likely etiological importance of human–pathogen co-evolution, attempts at empirical confirmation have been rare. Indeed, “proof” of co-evolution poses a formidable challenge, requiring a demonstration of increased reproductive fitness in each species driven by reciprocal changes in two genomes over time (Woolhouse et al., 2002). Although these criteria have been met in laboratory studies and in some natural populations (Lenski and Levin, 1985; Little, 2002; Little et al., 2006), a similarly rigorous assessment of human–pathogen co-evolution must accommodate long generation times and the genetic and phenotypic complexity of the human traits under selection. Nonetheless, substantial phenomenological evidence consistent with human–pathogen co-evolution now exists, including evidence of spatial patterns of parallel genetic variation between species, and of correlated functional changes at the molecular level (Kraaijeveld et al., 1998; Lively and Dybdahl, 2000; Funk et al., 2000; Woolhouse et al., 2002). The collection of high-density genomic data in paired human–pathogen samples and improvements in phenotypic data, as well as advances in pathogen genomics, should soon enable more explicit tests of the concept.
Our aim here is to summarize the growing body of evidence in favor of the hypothesis that genetic interactions driven by host and pathogen co-evolution can have significant implications for genetic epidemiological studies and biomedicine. While this is not a novel hypothesis, it remains understudied. We also underscore how recent advances in genomic technology provide new opportunities to test for genome-by-genome interactions, and offer suggestions on how to incorporate them into more accurate genetic models of disease.
HELICOBACTER PYLORI
Studies of Helicobacter pylori provide perhaps the best evidence in favor of human–pathogen co-evolution, and distinctly illustrate the power of the modern genetic toolkit to investigate it. H. pylori chronically infects the gastric epithelia of half the world’s population, causing peptic ulcers in 10–20% of those infected, and distal gastric carcinoma in ∼1% (Peek and Blaser, 2002; Jemal et al., 2011). The majority of individuals infected, however, suffer only from superficial gastritis in adulthood, while likely gaining protection against diseases such as esophageal cancer and reflux esophagitis, and more controversially, childhood asthma and diarrhea (Rothenbacher et al., 2000; Vaezi et al., 2000; ). That H. pylori should have a largely innocuous and potentially symbiotic relationship with its host follows from co-evolutionary theory, based on its vertical mode of transmission, its long-term colonization of a single host, and its ∼50,000 year association with Homo sapiens (Rothenbacher et al., 2002; Moodley et al., 2012). Why a fraction of individuals develop life-threatening clinical disease, on the other hand, requires explanation, with one possibility being the disruption of long-standing co-evolutionary relationships.
Although H. pylori-mediated diseases often advance to the clinical stage in late adulthood, their onset typically occurs during reproductive years (; Susser and Stein, 2002). Importantly, a disease need not have an especially large selection coefficient to shape allele frequency distributions in populations, especially over thousands of years (Ewald and Cochran, 2000). In fact, the historical fitness load of peptic ulcers, obtained by multiplying prevalence by selection coefficient, has been estimated to be similar to those for infectious diseases such as meningitis and rubella (). Also consistent with co-evolutionary theory is the fact that H. pylori-mediated gastric diseases occur disproportionately in men (Susser and Stein, 2002; Engel et al., 2003); H. pylori is usually, but not necessarily, transmitted by the mother, such that female fitness has likely exerted a stronger constraint against H. pylori virulence.
Some H. pylori virulence factors appear to increase the risk of serious clinical outcome regardless of host genotype. The cag pathogenicity island, present in some strains, encodes a type IV secretion system, and VacA encodes a pore-forming cytotoxin. Both have been implicated as carcinogenic risk factors, though neither is a necessary nor sufficient one (Wroblewski et al., 2010). Other virulence factors released by H. pylori include urease, which facilitates neutralization of the otherwise forbidding acidity of the gastric mucosa; NAP, which enables iron uptake; and arginase, which helps H. pylori subvert host macrophages. These, like most H. pylori virulence factors, operate to create a basal inflammatory state without generating an excessive immune response. Serious clinical disease reflects a disturbance of this balance (; ; Salama et al., 2013).
The maintenance of this balance also depends partly on human genetic factors (Lichtenstein et al., 2000; ; Mayerle et al., 2013a). Candidate gene studies on H. pylori-mediated diseases have implicated several gene polymorphisms that appear to affect risk, most notably in the interleukin-1 (IL-1) family of cytokines (Schneider et al., 2008). Recently, two GWAS assessing susceptibility to gastric cancer and H. pylori infection identified SNPs with odds ratios ranging from 1.3 to 1.4, mostly of uncertain biological function (Shi et al., 2011; El-Omar, 2013; Mayerle et al., 2013b, Table 1). These polymorphisms account for only a small proportion of the estimated heritability of disease phenotypes.
Table 1
| Disease/trait | Gene | SNP | Cases/controls | Population | p-value | OR1 | 95% CI2 | Reference |
|---|---|---|---|---|---|---|---|---|
| Gastric cancer | ZBTB20 | rs9841504 | 1006/2273 | Chinese | 1.7E-09 | 0.76 | [0.69–0.83] | Shi et al. (2011) |
| Gastric cancer | PRKAA1 | rs13361707 | 1006/2273 | Chinese | 7.6E-29 | 1.41 | [1.32–1.49] | Shi et al. (2011) |
| H. pylori serologic status | TLR10 | rs10004195 | 2623/7862 | European | 1.4E-18 | 0.70 | [0.65–0.76] | Mayerle et al. (2013b) |
| H. pylori serologic status | FCGR2A | rs368433 | 2623/7862 | European | 2.1E-08 | 0.73 | [0.65–0.85] | Mayerle et al. (2013b) |
| Tuberculosis | RCN1–WT1 | rs2057178 | 2127/5636 | African | 2.6E-09 | 0.77 | [0.71–0.84] | Thye et al. (2012) |
| Tuberculosis | RPS4XP18–UBE2CP2 | rs4331426 | 2237/3122 | African | 6.8E-09 | 1.19 | [1.13–1.27] | Thye et al. (2010) |
| Cervical cancer | EXOC1 | rs13117307 | 1364/3028 | Chinese | 9.7E-09 | 1.26 | [1.16–1.36] | Shi et al. (2013) |
| Cervical cancer | HLA-DPB2 | rs4282438 | 1364/3028 | Chinese | 4.5E-27 | 0.75 | [0.71–0.79] | Shi et al. (2013) |
| Cervical cancer | ZPBP2–GSDMB | rs8067378 | 1364/3028 | Chinese | 2.0E-08 | 1.18 | [1.11–1.25] | Shi et al. (2013) |
| Cervical cancer | – | rs9277952 | 1364/3028 | Chinese | 2.3E-09 | 0.85 | [0.81–0.90] | Shi et al. (2013) |
| Cervical cancer | MICA | rs2516448 | 2174/5002 | European | 1.6E-18 | 1.42 | [1.31–1.54] | |
| Cervical cancer | HLA-DRB1–HLA-DQA1 | rs9272143 | 2174/5006 | European | 9.3E-24 | 0.67 | [0.62–0.72] | |
| Cervical cancer | HLA-DPB2 | rs3117027 | 2171/4986 | European | 4.9E-08 | 1.25 | [1.15–1.35] |
Genetic variants identified by GWAS for phenotypes related to infection by H. pylori, M. tuberculosis, and human papillomavirus.
OR, odds ratio.
CI, confidence interval.
Studies of human or H. pylori genetics in isolation have generally failed to explain why populations with similar rates ofH. pylori infection exhibit strikingly different susceptibilities to gastric cancer. For example, in many African and South Asian countries, the low incidences of gastric cancer in the presence of almost universal rates of H. pylori infection remain a source of much speculation, and have been referred to collectively as the “African enigma” and the “Asian enigma” (Holcombe, 1992; ; Ghoshal et al., 2007). In Latin America, where H. pylori strains native to Amerindian populations have been largely displaced by European strains (Dominguez-Bello et al., 2008; ), the predominantly Amerindian populations living at high altitudes suffer disproportionately from gastric cancer relative to other populations with similar infection rates (; Torres et al., 2013). These and other points of evidence raise the possibility that the pathogenicity of a given H. pylori strain may vary with human genomic variation, and that some individuals may be better adapted to their infecting strains than others.
Modern genomic techniques have made the assessment of such hypotheses feasible. Over the past two decades, a comprehensive phylogeography of H. pylori has been constructed using multilocus sequence typing (MLST), a procedure by which polymorphisms in fragments from housekeeping genes are used to characterize bacterial isolates (Maiden et al., 1998). Analyses of samples from around the world have revealed a strong concordance between H. pylori phylogenetic clusters and the geographical locations from which they are derived (Falush et al., 2003; Moodley and Linz, 2009; Moodley et al., 2009). Ancestral H. pylori sequences inferred using MLST data also correspond to geographically defined human populations (Falush et al., 2003; Moodley et al., 2012). The typical modern H. pylori chromosome is now understood to be an amalgam of fragments from multiple ancestral sequences, a consequence of H. pylori’s high recombinogenicity (Suerbaum et al., 1998; Falush et al., 2003). The genome of an H. pylori isolate can thus be quantitatively resolved into ancestral proportions, which correlate with proportions of human ancestry in admixed populations (Kodaman et al., 2014). In some cases, the ancestries of H. pylori isolates outperform human mitochondria in differentiating ethnic groups (Wirth et al., 2004).
These shared patterns of ancestry are unlikely to have arisen merely from parallel divergence due to founder effects or neutral drift. Certainly, the well-documented evolvability of functional loci within H. pylori strains, even within single individuals over a 6 year span, argues for the importance of adaptive microevolution (Israel et al., 2001; Dorer et al., 2009). Furthermore, at least 25% of known genes, including genes involved in mucosal adherence and the evasion of host immunity, are absent in someH. pylori strains isolated from different ethnic groups (Salama et al., 2000; Gressmann et al., 2005). In at least one case, variants of an H. pylori gene (babA2) encode adhesion proteins that exhibit host-specific effects, a hallmark of co-evolution. BabA binds to blood group antigens, triggering the release of proinflammatory cytokines. Notably, Amerindians, who almost all carry blood group O, harbor strains with a BabA variant that has up to a 1500-fold greater binding affinity to blood group O ().
If we conclude from these patterns of genetic covariation that co-evolution between humans and H. pylori has occurred and that it has promoted commensalism, then we may ask whether individuals who develop serious clinical disease have inherited mutually ill-adapted sets of host and pathogen alleles. Under this hypothesis, we should expect to find significant interactions between specific pairs of host and pathogen loci in disease models. Toward this end, candidate pairs of loci can be tested based on biochemical evidence of protein–protein interactions, such as those between the adhesin BabA and the Lewis(b) antigen, its epithelial receptor (). However, the effect size of any single two-locus interaction may be relatively small, as gastric disease etiology is phenotypically heterogeneous, and likely to be influenced by a large number of human and H. pylori genes (El-Omar, 2013). Thus, characterizing the relevant loci in a biologically meaningful way will ultimately require a systems biological approach.
We recently took a broad-based view to assess the impact of human – H. pylori co-evolution on gastric disease, using ancestry estimates from both humans and their H. pylori isolates in the absence of knowledge of specific interacting loci (Kodaman et al., 2014). Our study participants were recruited from two Colombian populations with highly different rates of gastric cancer, despite a nearly universal prevalence of H. pylori infection in both. We found that the low-risk human, coastal population was of admixed African, European, and Amerindian ancestry, whereas the high-risk, Andean population was mainly of Amerindian ancestry, with a minority of European ancestry. Severity of gastric disease correlated with the proportion of African H. pylori ancestry in patients with primarily Amerindian ancestry. On the other hand, patients with a large proportion of African human ancestry infected by African H. pylori strains had the best prognoses, consistent with ancestral coadaptation, and likely pertinent to the “African enigma.” The interaction between Amerindian human ancestry and African H. pylori ancestry accounted for the difference in disease risk between mountain and coastal populations, whereas even the well-known virulence factor, CagA, did not. These findings are thus consistent with the idea that neither human nor H. pylori genetic variation confers susceptibility or virulence per se, but only in context (Figure 1).
FIGURE 1
These findings also bring to light how understanding co-evolutionary interactions can inform and improve public health measures. It has been suggested that because H. pylori dominates the gastric microbiome in infected persons and has been shown to confer some beneficial effects, large-scale antibiotic eradication programs may not be warranted (
MYCOBACTERIUM TUBERCULOSIS COMPLEX
Another interesting candidate to study from a co-evolutionary perspective is Mycobacterium tuberculosis (Mtb) and closely related species, believed to have co-existed with anatomically modern humans for ∼70,000 years (
In contrast to H. pylori, Mtb is transmitted horizontally, and must cause active disease to be transmitted (e.g., via coughing or sneezing). Because Mtb transmission increases with virulence, evolutionary theory predicts that strong selective pressures should favor increased virulence until the number of transmissions per infected host reaches a fitness-reducing limit (Knolle, 1989; Frank and Schmid-Hempel, 2008). Such a limit necessarily depends on population-specific parameters, of which host density is probably the most important (
Before advances in genotyping technology improved strain classification, the apparent genetic homogeneity of Mtb led investigators to believe that variation in disease outcome depended primarily on environmental and human genetic factors (Galagan, 2014). Twin and adoption studies provided compelling evidence for the involvement of human genetic variation as a risk modifier (
Since the advent of PCR-based genotyping techniques, it has become increasingly clear that Mtb genetic variation is non-trivial and clinically consequential (Malik and Godfrey-Faussett, 2005; Nicol and Wilkinson, 2008). Most notably, strains now recognized as part of the “Beijing family,” first genotyped in the 1990s following several drug-resistant outbreaks, have been found to exhibit greater efficiency of transmission and to cause more severe disease phenotypes in many animal models (Glynn et al., 2002; Reed et al., 2004; Parwati et al., 2010). Whole-genome sequencing of a large number of clinical Mtb isolates has since revealed over 30,000 Mtb SNPs, a large proportion of which are non-synonymous (
High-throughput sequence data have also enabled the construction of a robust phylogenetic tree, the major branches of which parallel human mitochondrial phylogeny (
It is thus likely that genetic factors in both Mtb and humans influence a wide range of TB phenotypes, including those pertaining to infectivity, progression from latent to active disease, and effectiveness of treatment (
Despite being an ancient strain with ample opportunity to spread beyond West Africa, M. africanum has not done so, possibly indicating host-specific adaptation (
While the above discussion has focused mainly on pulmonary TB, we note here that extra-pulmonary TB, a less common and more severe form of disease, may be especially amenable to analyses guided by co-evolutionary hypotheses. This form of the disease leads more quickly to fatality and results in fewer transmissions than the pulmonary form (Sharma and Mohan, 2004), which probably represents a non-optimal outcome in terms of Mtb fitness. However, data on extra-pulmonary TB to support co-evolutionary hypotheses – especially historical data pre-dating the antibiotic era and the HIV epidemic – are at present lacking (Tiemersma et al., 2011).
HUMAN PAPILLOMAVIRUS
Human papillomavirus (HPV) is the most common sexually transmitted infectious agent in the world, and the second most common infectious cause of cancer after H. pylori (
Papillomaviruses (PVs) are notable for their slow rate of evolution relative to other pathogens – only an order of magnitude higher than humans, in the case of HPV (Ong et al., 1993; Rector et al., 2007; Shah et al., 2010). This is commonly attributed to their use of high-fidelity host replication mechanisms (Van Doorslaer, 2013). A slow evolutionary rate precludes rapid adaptation to new hosts, and PV strains correspondingly show little evidence of inter-species transmission or related horizontal gene transfer (Herbst et al., 2009; Shah et al., 2010; Van Doorslaer, 2013). All carcinogenic types of HPV belong to a single genus of papillomaviruses that diverged from a common ancestor about 75 million years ago, predating the primate lineage (Rector et al., 2007; Van Doorslaer, 2013). By the emergence of H. sapiens, the common ancestor of HPV 16 and HPV 18 had diverged into separate species, and in fact HPV 16 and HPV 18 had already diverged from all other HPV types within their respective species clades (Lewin, 1993; Ong et al., 1993). Given this combination of early divergence, slow evolution, and strict host specialization, we would expect variants within HPV types independently to have similar phylogeographic patterns to that of H. sapiens. Global data collected for the two most frequently sexually transmitted types, HPV 16 and 18, reflect such a pattern (
Biochemical and bioinformatic analyses indicate that HPV evolution has not been entirely neutral. Viral genes expressed early during a PV infection, for example, appear to have evolved at different rates than those expressed late (Garcia-Vallve et al., 2005; Rector et al., 2007). Although most PV genes show signs of strong purifying selection, the exceptions appear to be important (
Because the integration of the HPV genome within the human genome is permanent, death of the host ends all possibility of viral multiplication and transmission. Even strains that damage the health of the host sufficiently to reduce human-to-human sexual contact can suffer a competitive disadvantage. Therefore, both host and pathogen should cooperate to prevent severe disease. As with H. pylori and MTB, there is some empirical evidence supporting the idea that humans and HPV types co-evolved to limit tumorigenesis, and that evolutionarily mismatched strains may be driving severe clinical outcomes. A study of high-grade cervical intraepithelial neoplasia (CIN) and invasive cervical cancer in an Italian cohort of Caucasian women demonstrated that non-European variants of HPV16, Af1 and AA, were found at an increased frequency in invasive lesions (Tornesello et al., 2004). A separate study of mostly Caucasian (81%) female university students in the United States showed that those infected with non-European HPV 16 variants were 6.5 times more likely to develop high-grade CIN than those with European variants (Xi et al., 1997). The same study demonstrated a similar HPV 16-related risk profile (4.5 relative risk) in a predominantly Caucasian (79%) population of women presenting at a sexually transmitted disease clinic (Xi et al., 1997). Finally, at the molecular level, there is some evidence that variants of the HPV 16 E6 protein, described above, may be better adapted for replication within specific hosts (
DISCUSSION
Taken together, the three examples above illustrate how co-evolution can promote a reduction in antagonism between pathogen and host, and in doing so leave discernible signatures on the genomes of both species. If, as we argue here, the disruption of historical co-evolutionary relationships can explain many differences in disease outcomes, knowledge of the conditions under which such relationships arise and dissolve will be helpful in defining genetic architecture of disease etiology. The applicability of this model depends, to a large extent, on the degree of integration between host and pathogen genomes, which can take many forms.
A long-standing association between humans and pathogens may be a necessary factor for cross-genomic integration, as with the three pathogens we have discussed. In contrast, many infectious diseases that occur epidemically are caused by zoonotic pathogens for which the human host is an evolutionary dead end, such as Salmonella enterica and Borrelia burgdorferi (Sokurenko et al., 2006; Falush, 2009). Other pathogens have had limited occasion to co-evolve with humans, because they cause disease primarily on an opportunistic basis (e.g., Streptococcus pneumonia or Clostridium difficile) or over a broad range of hosts (e.g., Toxoplasma gondii) (
For human-specific pathogens that cause endemic diseases and are not recent, the likelihood that severe disease is the outcome of a co-evolutionary mismatch should increase with the overlap between host and pathogen fitness. The pathogenicity of vertically transmitted pathogens, for example, should decrease over time, because such pathogens often depend on host survival (and possibly reproduction) for transmission. However, a strong overlap between host and pathogen fitness can also exist in the absence of vertical transmission. A horizontally transmitted pathogen, such as HPV, can evolve to be largely benign insofar as it depends on a healthy host for transmission.
When a pathogen’s fitness depends on its ability to cause damage to its human host, as with Mtb, attenuated antagonism becomes a special case, and its disruption becomes more difficult to detect and requires more evidence to confirm. While Mtb strains that increase the duration of a transmissible state will generally have a competitive advantage, the optimal duration can be expected to vary based on many population-level parameters, such as host density. This probably explains why modern Mtb lineages that are more common in high-density urban populations exhibit greater virulence. On the other hand, if horizontal transfer is confined to small, isolated populations, it may be considered effectively vertical. With such pathogens, a better understanding of the co-evolutionary history will be necessary to infer whether severe disease is caused by disrupted co-evolution or by another factor, such as infection by a universally more virulent strain or an opportunistic infection in an immunosuppressed patient.
The life history of the pathogen is also important in assessing the possibility and nature of co-evolution. A pathogen typically faces a tradeoff between fecundity and longevity. Increased fecundity within a host increases the probability (or rate) of transmission, but may negatively affect host lifespan or mobility (Frank and Schmid-Hempel, 2008). Therefore, a pathogen’s position on the continuum between greater fecundity and increased longevity will often reflect the degree to which its fitness depends on the health of the host. The case of HPV is somewhat of an exception in this regard. Host immune responses can induce diverse strategies, creating HPV types that are highly fecund, or less fecund with few virions per host. Whereas highly fecund types are more likely to transmit, they are also more likely to induce a vigorous immune response leading to clearance. Low fecundity types on the other hand, are more likely to persist as subclinical infections that can lead to prolonged inflammation and eventually cancer (
Another factor influencing the applicability of the model we propose is a pathogen’s recombinogenicity. In theory, a pathogen that recombines freely is more likely to be panmictic, and hence less likely to co-evolve with a particular human host population (
Recombination can also occur via horizontal gene transfer, as among species within the microbiome (Smillie et al., 2011; Ravel et al., 2011; Liu et al., 2012). This would suggest that co-evolution might be a relatively weak force in shaping microbiotal genetic variation. However, data possibly supporting human–microbiome co-evolution exist; for example, the strongest correlate of an individual’s microbiotal identity is ethnicity (
CONCLUSION
While the prospect of introducing co-evolutionary interactions into genetic epidemiology models may appear to add a new layer of complexity to an already difficult problem, a co-evolutionary perspective should help us construct more precise and accurate hypotheses, improving our ability to find real and reproducible results. Importantly, co-evolved genes will not be neutral in either species, which may make their identification easier. Although many methods exist to find loci that are candidates to have evolved under selection (
If the ultimate goal is to find interacting genes that have co-evolved to be benign and are subsequently disrupted in disease, we will need to identify differential patterns of concerted selection in paired human and pathogenic loci from different populations. The limiting factor to the development of appropriate methods toward this end has probably been the lack of prospectively collected paired genetic data for humans and pathogens. Once these data are available, existing methods to detect epistasis within a species can be adapted for cross-species analyses in the absence of a priori biological hypotheses. Where evidence for selection exists, genetic variants can be filtered prior to analyses to detect epistasis. Framing hypotheses in the context of biochemical and bioinformatic functional evidence or pre-existing evidence for association can hone study design even further. For example, using paired data and pathogenic genetic variation as the outcome variable, novel epitopes have been discovered in association studies (
Finally, we should note that the ultimate impact of this approach may extend beyond infectious diseases to what are traditionally considered non-communicable diseases. For example, we now recognize that both gastric and cervical cancers, as well as atherosclerosis, may have origins in infection (Libby et al., 2002; Porta et al., 2011). The number of such examples will certainly expand.
Statements
Acknowledgments
This study was supported by the National Center for Research Resources, grant UL1 RR024975-01, which is now at the National Center for Advancing Translational Sciences; National Cancer Institute Grant P01 CA28842, the Vanderbilt-Ingram Cancer Center, the Wendy Dio family and the TJ. Martell Foundation. Scott M. Williams was partially supported by P20 GM103534.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
REFERENCES
1
AbecasisG. R.AutonA.BrooksL. D.DePristoM. A.DurbinR. M.HandsakerR. E.et al (2012). An integrated map of genetic variation from 1,092 human genomes.Nature49156–65. 10.1038/nature11632
2
AchtmanM. (2008). Evolution, population structure, and phylogeography of genetically monomorphic bacterial pathogens.Annu. Rev. Microbiol.6253–70. 10.1146/annurev.micro.62.081307.162832
3
AguiletaG.RefregierG.YocktengR.FournierE.GiraudT. (2009). Rapidly evolving genes in pathogens: methods for detecting positive selection and examples among fungi, bacteria, viruses and protists.Infect. Genet. Evol.9656–670. 10.1016/j.meegid.2009.03.010
4
AjzenbergD.BanulsA. L.SuC.DumetreA.DemarM.CarmeB.et al (2004). Genetic diversity, clonality and sexuality in Toxoplasma gondii.Int. J. Parasitol.341185–1196. 10.1016/j.ijpara.2004.06.007
5
AndersonR. M.MayR. M. (1982). Coevolution of hosts and parasites.Parasitology85411–426. 10.1017/S0031182000055360
6
Aspholm-HurtigM.DailideG.LahmannM.KaliaA.IlverD.VikströmS.et al (2004). Functional adaptation of BabA, the H. pylori ABO blood group antigen binding adhesin.Science305519–522. 10.1126/science.1098801
7
BackstromA.LundbergC.KersulyteD.BergD. E.BorenT.AmqvistA. (2004). Metastability of Helicobacter pylori bab adhesin genes and dynamics in Lewis b antigen binding.Proc. Natl. Acad. Sci. U.S.A.10116923–16928. 10.1073/pnas.0404817101
8
BaldariC. T.LanzavecchiaA.TelfordJ. L. (2005). Immune subversion by Helicobacter pylori.Trends Immunol.26199–207. 10.1016/j.it.2005.01.007
9
BarreiroL. B.NeyrollesO.BabbC. L.TailleuxL.QuachH.McElreaveyK.et al (2006). Promoter variation in the DC-SIGN-encoding gene CD209 is associated with tuberculosis.PLoS Med.3:e20. 10.1371/journal.pmed.0030020
10
BarryC. E.IIIBoshoffH. I.DartoisV.DickT.EhrtS.FlynnJ.et al (2009). The spectrum of latent tuberculosis: rethinking the biology and intervention strategies.Nat. Rev. Microbiol.7845–855. 10.1038/nrmicro2236
11
BarthaI.CarlsonJ. M.BrummeC. J.McLarenP. J.BrummeZ. L.JohnM.et al (2013). A genome-to-genome analysis of associations between human genetic variation, HIV-1 sequence diversity, and viral control.Elife (Cambridge)2:e01123. 10.7554/eLife.01123
12
BensonA. K.KellyS. A.LeggeR.MaF.LowS. J.KimJ.et al (2010). Individuality in gut microbiota composition is a complex polygenic trait shaped by multiple environmental and host genetic factors.Proc. Natl. Acad. Sci. U.S.A.10718933–18938. 10.1073/pnas.1007028107
13
BernardH. U. (1994). Coevolution of papillomaviruses with human populations.Trends Microbiol.2140–143. 10.1016/0966-842X(94)90602-5
14
BernardH. U.BurkR. D.ChenZ.van DoorslaerK.zur HausenH.de VilliersE. M. (2010). Classification of papillomaviruses (PVs) based on 189 PV types and proposal of taxonomic amendments.Virology40170–79. 10.1016/j.virol.2010.02.002
15
BikE. M.EckburgP. B.GillS. R.NelsonK. E.PurdomE. A.FrancoisF.et al (2006). Molecular analysis of the bacterial microbiota in the human stomach.Proc. Natl. Acad. Sci. U.S.A.103732–737. 10.1073/pnas.0506655103
16
BlaserM. J.KirschnerD. (2007). The equilibria that allow bacterial persistence in human hosts.Nature449843–849. 10.1038/nature06198
17
BlaserM. J.ChenY.ReibmanJ. (2008). Does Helicobacter pylori protect against asthma and allergy?Gut57561–567. 10.1136/gut.2007.133462
18
BullJ. J.MolineuxI. J.RiceW. R. (1991). Selection of benevolence in a host–parasite system.Evolution45875–882. 10.2307/2409695
19
CampbellD. I.WarrenB. F.ThomasJ. E.FiguraN.TelfordJ. L.SullivanP. B. (2001). The African enigma: low prevalence of gastric atrophy, high prevalence of chronic inflammation in West African adults and children.Helicobacter6263–267. 10.1046/j.1083-4389.2001.00047.x
20
CarterJ. J.MadeleineM. M.SheraK.SchwartzS. M.Cushing-HaugenK. L.WipfG. C.et al (2001). Human papillomavirus 16 and 18 L1 serology compared across anogenital cancer sites.Cancer Res.611934–1940.
21
Carvajal-RodriguezA. (2008). Detecting recombination and diversifying selection in human alpha-papillomavirus.Infect Genet. Evol.8689–692. 10.1016/j.meegid.2008.07.002
22
CasanovaJ. L.AbelL. (2007). Human genetics of infectious diseases: a unified theory.EMBO J.26915–922. 10.1038/sj.emboj.7601558
23
Cavalli-SforzaL. L. (2005). The Human Genome Diversity Project: past, present and future.Nat. Rev. Genet.6333–340. 10.1038/nrg1596
24
CawsM.ThwaitesG.DunstanS.HawnT. R.LanN. T.ThuongN. T.et al (2008). The influence of host and bacterial genotype on the development of disseminated disease with Mycobacterium tuberculosis.PLoS Pathog.4:e1000034. 10.1371/journal.ppat.1000034
25
ChenD.Juko–PecirepI.HammerJ.IvanssonE.EnrothS.GustavssonI.et al (2013). Genome-wide association study of susceptibility loci for cervical cancer.J. Natl. Cancer Inst.105624–633. 10.1093/jnci/djt051
26
ChenZ.TeraiM.FuL.HerreroR.DeSalleR.BurkR. D. (2005). Diversifying selection in human papillomavirus type 16 lineages based on complete genome analyses.J. Virol.797014–7023. 10.1128/JVI.79.11.7014-7023.2005
27
ChibaT.SenoH.MarusawaH.WakatsukiY.OkazakiK. (2006). Host factors are important in determining clinical outcomes of Helicobacter pylori infection.J. Gastroenterol.411–9. 10.1007/s00535-005-1743-4
28
ChimusaE. R.ZaitlenN.DayaM.MollerM.van HeldenP. D.MulderN. J.et al (2014). Genome-wide association study of ancestry-specific TB risk in the South African coloured population.Hum. Mol. Genet.23796–809. 10.1093/hmg/ddt462
29
CochranG. M.EwaldP. W.CochranK. D. (2000). Infectious causation of disease: an evolutionary perspective.Perspect. Biol. Med.43406–448. 10.1353/pbm.2000.0016
30
ComasI.GagneuxS. (2011). A role for systems epidemiology in tuberculosis research.Trends Microbiol.19492–500. 10.1016/j.tim.2011.07.002
31
ComasI.CoscollaM.LuoT.BorrellS.HoltK. E.Kato-MaedaM.et al (2013). Out-of-Africa migration and Neolithic coexpansion of Mycobacterium tuberculosis with modern humans.Nat. Genet.451176–1182. 10.1038/ng.2744
32
ComstockG. W. (1978). Tuberculosis in twins: a re-analysis of the Prophit survey.Am. Rev. Respir. Dis.117621–624.
33
CorreaP.CuelloC.DuqueE.BurbanoL. C.GarciaF. T.BotanosO.et al (1976). Gastric cancer in Colombia. III. Natural history of precursor lesions.J. Natl. Cancer Inst.571027–1035. 10.1093/jnci/57.5.1027
34
CorreaP.PiazueloM. B. (2012). Evolutionary history of the Helicobacter pylori genome: implications for gastric carcinogenesis.Gut Liver621–28. 10.5009/gnl.2012.6.1.21
35
CoscollaM.GagneuxS. (2010). Does M. tuberculosis genomic diversity explain disease diversity?Drug Discov. Today Dis. Mech.7:e43–e59. 10.1016/j.ddmec.2010.09.004
36
DeFilippisV. R.AyalaF. J.VillarrealL. P. (2002). Evidence of diversifying selection in human papillomavirus type 16 E6 but not E7 oncogenes.J. Mol. Evol.55491–499. 10.1007/s00239-002-2344-y
37
de JongB. C.AntonioM.GagneuxS. (2010). Mycobacterium africanum – review of an important cause of human tuberculosis in West Africa.PLoS Negl. Trop. Dis.4:e744. 10.1371/journal.pntd.0000744
38
de JongB. C.HillP. C.AikenA.AwineT.AntonioM.AdetifaI. M.et al (2008). Progression to active tuberculosis, but not transmission, varies by Mycobacterium tuberculosis lineage in The Gambia.J. Infect. Dis.1981037–1043. 10.1086/591504
39
de MartelC.FerlayJ.FranceschiS.VignatJ.BrayF.FormanD.et al (2012). Global burden of cancers attributable to infections in 2008: a review and synthetic analysis.Lancet Oncol.13607–615. 10.1016/S1470-2045(12)70137-7
40
de SabletT.PiazueloM. B.ShafferC. L.SchneiderB. G.AsimM.ChaturvediR.et al (2011). Phylogeographic origin of Helicobacter pylori is a determinant of gastric cancer risk.Gut601189–1195. 10.1136/gut.2010.234468
41
Dominguez-BelloM. G.PerezM. E.BortoliniM. C.SalzanoF. M.PericchiL. R.Zambrano-GuzmánO.et al (2008). Amerindian Helicobacter pylori strains go extinct, as european strains expand their host range.PLoS ONE3:e3307. 10.1371/journal.pone.0003307
42
DoorbarJ. (2006). Molecular biology of human papillomavirus infection and cervical cancer.Clin. Sci. (Lond.)110525–541. 10.1042/CS20050369
43
DorerM. S.TalaricoS.SalamaN. R. (2009). Helicobacter pylori’s unconventional role in health and disease.PLoS Pathog.5:e1000544. 10.1371/journal.ppat.1000544
44
DyeC.WilliamsB. G. (2010). The population dynamics and control of tuberculosis.Science328856–861. 10.1126/science.1185449
45
DysonN.HowleyP. M.MungerK.HarlowE. (1989). The human papilloma virus-16 E7 oncoprotein is able to bind to the retinoblastoma gene product.Science243934–937. 10.1126/science.2537532
46
El-OmarE. M. (2013). Helicobacter pylori susceptibility in the GWAS era.JAMA3091939–1940. 10.1001/jama.2013.5590
47
EngelL. S.ChowW. H.VaughanT. L.GammonM. D.RischH. A.StanfordJ. L.et al (2003). Population attributable risks of esophageal and gastric cancers.J. Natl. Cancer Inst.951404–1413. 10.1093/jnci/djg047
48
EwaldP. W.CochranG. M. (2000). Chlamydia pneumoniae and cardiovascular disease: an evolutionary perspective on infectious causation and antibiotic treatment.J. Infect. Dis. 181(Suppl. 3), S394–S401. 10.1086/315602
49
FalushD. (2009). Toward the use of genomics to study microevolutionary change in bacteria.PLoS Genet.5:e1000627. 10.1371/journal.pgen.1000627
50
FalushD.WirthT.LinzB.PritchardJ. K.StephensM.KiddM.et al (2003). Traces of human migrations in Helicobacter pylori populations.Science2991582–1585. 10.1126/science.1080857
51
FennerL.EggerM.BodmerT.FurrerH.BallifM.BattegayM.et al (2013). HIV infection disrupts the sympatric host–pathogen relationship in human tuberculosis.PLoS Genet.9:e1003318. 10.1371/journal.pgen.1003318
52
FrankS. A. (1996). Models of parasite virulence.Q. Rev. Biol.7137–78. 10.1086/419267
53
FrankS. A.Schmid-HempelP. (2008). Mechanisms of pathogenesis and the evolution of parasite virulence.J. Evol. Biol.21396–404. 10.1111/j.1420-9101.2007.01480.x
54
FunkD. J.HelblingL.WernegreenJ. J.MoranN. A. (2000). Intraspecific phylogenetic congruence among multiple symbiont genomes.Proc. Biol. Sci.2672517–2521. 10.1098/rspb.2000.1314
55
GagneuxS. (2012). Host–pathogen coevolution in human tuberculosis.Philos. Trans. R. Soc. Lond. B Biol. Sci.367850–859. 10.1098/rstb.2011.0316
56
GagneuxS.DeRiemerK.VanT.Kato-MaedaM.de JongB. C.NarayananS.et al (2006). Variable host–pathogen compatibility in Mycobacterium tuberculosis.Proc. Natl. Acad. Sci. U.S.A.1032869–2873. 10.1073/pnas.0511240103
57
GalaganJ. E. (2014). Genomic insights into tuberculosis.Nat. Rev. Genet.15307–320. 10.1038/nrg3664
58
GaoL.TaoY.ZhangL.JinQ. (2010). Vitamin D receptor genetic polymorphisms and tuberculosis: updated systematic review and meta-analysis.Int. J. Tuberc. Lung Dis.1415–23.
59
Garcia-VallveS.AlonsoA.BravoI. G. (2005). Papillomaviruses: different genes have different histories.Trends Microbiol.13514–521. 10.1016/j.tim.2005.09.003
60
GhoshalU. C.TripathiS.GhoshalU. (2007). The Indian enigma of frequent H. pylori infection but infrequent gastric cancer: is the magic key in Indian diet, host’s genetic make up, or friendly bug?Am. J. Gastroenterol.1022113–2114. 10.1111/j.1572-0241.2007.01324_13.x
61
GlynnJ. R.WhiteleyJ.BifaniP. J.KremerK.van SoolingenD. (2002). Worldwide occurrence of Beijing/W strains of Mycobacterium tuberculosis: a systematic review.Emerg. Infect. Dis.8843–849. 10.3201/eid0805.020002
62
GressmannH.LinzB.GhaiR.PleissnerK. P.SchlapbachR.YamaokaY.et al (2005). Gain and loss of multiple genes during the evolution of Helicobacter pylori.PLoS Genet.1:e43. 10.1371/journal.pgen.0010043
63
GriggM. E.BonnefoyS.HehlA. B.SuzukiY.BoothroydJ. C. (2001). Success and virulence in toxoplasma as the result of sexual recombination between two distinct ancestries.Science294161–165. 10.1126/science.1061888
64
HeitmanJ. (2006). Sexual reproduction and the evolution of microbial pathogens.Curr. Biol.16R711–R725. 10.1016/j.cub.2006.07.064
65
HerbF.ThyeT.NiemannS.BrowneE. N.ChinbuahM. A.GyapongJ.et al (2008). ALOX5 variants associated with susceptibility to human pulmonary tuberculosis.Hum. Mol. Genet.171052–1060. 10.1093/hmg/ddm378
66
HerbstL. H.LenzJ.Van DoorslaerK.ChenZ.StacyB. A.WellenhanJ. F.Jr.et al (2009). Genomic characterization of two novel reptilian papillomaviruses, Chelonia mydas papillomavirus 1 and Caretta caretta papillomavirus 1.Virology383131–135. 10.1016/j.virol.2008.09.022
67
HershbergR.LipatovM.SmallP. M.ShefferH.NiemannS.HomolkaS.et al (2008). High functional diversity in Mycobacterium tuberculosis driven by genetic drift and human demography.PLoS Biol.6:e311. 10.1371/journal.pbio.0060311
68
HildesheimA.WangS. S. (2002). Host and viral genetics and risk of cervical cancer: a review.Virus Res.89229–240. 10.1016/S0168-1702(02)00191-0
69
HillA. V. (2001). The genomics and genetics of human infectious disease susceptibility.Annu. Rev. Genomics Hum. Genet.2373–400. 10.1146/annurev.genom.2.1.373
70
HillA. V. (2012). Evolution, revolution and heresy in the genetics of infectious disease susceptibility.Philos. Trans. R. Soc. Lond. B Biol. Sci.367840–849. 10.1098/rstb.2011.0275
71
HoalE. G.LewisL. A.JamiesonS. E.TanzerF.RossouwM.VictorT.et al (2004). SLC11A1 (NRAMP1) but not SLC11A2 (NRAMP2) polymorphisms are associated with susceptibility to tuberculosis in a high-incidence community in South Africa.Int. J. Tuberc. Lung Dis.81464–1471.
72
HoL.ChanS. Y.BurkR. D.DasB. C.FujinagaK.IcenogleJ. P.et al (1993). The genetic drift of human papillomavirus type 16 is a means of reconstructing prehistoric viral spread and the movement of ancient human populations.J. Virol.676413–6423.
73
HolcombeC. (1992). Helicobacter pylori: the African enigma.Gut33429–431. 10.1136/gut.33.4.429
74
HuibregtseJ. M.ScheffnerM.HowleyP. M. (1993a). Cloning and expression of the cDNA for E6-AP, a protein that mediates the interaction of the human papillomavirus E6 oncoprotein with p53.Mol. Cell. Biol.13775–784.
75
HuibregtseJ. M.ScheffnerM.HowleyP. M. (1993b). Localization of the E6-AP regions that direct human papillomavirus E6 binding, association with p53, and ubiquitination of associated proteins.Mol. Cell. Biol.134918–4927.
76
Human Microbiome Project Consortium. (2012). Structure, function and diversity of the healthy human microbiome.Nature486207–214. 10.1038/nature11234
77
HungI. F.WongB. C. (2009). Assessing the risks and benefits of treating Helicobacter pylori infection.Therap. Adv. Gastroenterol.2141–147. 10.1177/1756283X08100279
78
IntemannC. D.ThyeT.NiemannS.BrowneE. N.Amanua ChinbuahM.EnimilA.et al (2009). Autophagy gene variant IRGM-261T contributes to protection from tuberculosis caused by Mycobacterium tuberculosis but not by M. africanum strains.PLoS Pathog.5:e1000577. 10.1371/journal.ppat.1000577
79
IsraelD. A.SalamaN.KrishnaU.RiegerU. M.AthertonJ. C.FalkowS.et al (2001). Helicobacter pylori genetic diversity within the gastric niche of a single human host.Proc. Natl. Acad. Sci. U.S.A.9814625–14630. 10.1073/pnas.251551698
80
JacobF. (1977). Evolution and tinkering.Science1961161–1166. 10.1126/science.860134
81
JallowM.TeoY. Y.SmallK. S.RockettK. A.DeloukasP.ClarkT. G.et al (2009). Genome-wide and fine-resolution association analysis of malaria in West Africa.Nat. Genet.41657–665. 10.1038/ng.388
82
JemalA.BrayF.CenterM. M.FerlayJ.WardE.FormanD. (2011). Global cancer statistics.CA Cancer J. Clin.6169–90. 10.3322/caac.20107
83
KarlssonE. K.KwiatkowskiD. P.SabetiP. C. (2014). Natural selection and infectious disease in human populations.Nat. Rev. Genet.15379–393. 10.1038/nrg3734
84
KlingelhutzA. J.RomanA. (2012). Cellular transformation by human papillomaviruses: lessons learned by comparing high- and low-risk viruses.Virology42477–98. 10.1016/j.virol.2011.12.018
85
KnolleH. (1989). Host density and the evolution of parasite virulence.J. Theor. Biol.136199–207. 10.1016/S0022-5193(89)80226-7
86
KoD. C.UrbanT. J. (2013). Understanding human variation in infectious disease susceptibility through clinical and cellular GWAS.PLoS Pathog.9:e1003424. 10.1371/journal.ppat.1003424
87
KodamanN.PazosA.SchneiderB. G.PiazueloM. B.MeraR.SobotaR. S.et al (2014). Human and Helicobacter pylori coevolution shapes the risk of gastric disease.Proc. Natl. Acad. Sci. U.S.A.1111455–1460. 10.1073/pnas.1318093111
88
KraaijeveldA. R.Van AlphenJ. J.GodfrayH. C. (1998). The coevolution of host resistance and parasitoid virulence.Parasitology 116(Suppl. 1), S29–S45. 10.1017/S0031182000084924
89
LambrechtsL.FellousS.KoellaJ. C. (2006). Coevolutionary interactions between host and parasite genotypes.Trends Parasitol.2212–16. 10.1016/j.pt.2005.11.008
90
LenskiR. E.LevinB. R. (1985). Constraints on the coevolution of bacteria and virulent phage – a model, some experiments, and predictions for natural communities.Am. Nat.125585–602. 10.1086/284364
91
LewinR. (1993). Human Evolution3rd Edn.Boston: Blackwell Scientific Publications.
92
LewisS. J.BakerI.Davey SmithG. (2005). Meta-analysis of vitamin D receptor polymorphisms and pulmonary tuberculosis risk.Int. J. Tuberc. Lung Dis.91174–1177.
93
LibbyP.RidkerP. M.MaseriA. (2002). Inflammation and atherosclerosis.Circulation1051135–1143. 10.1161/hc0902.104353
94
LichtensteinP.HolmN. V.VerkasaloP. K.IliadouA.KaprioJ.KoskenvuoM.et al (2000). Environmental and heritable factors in the causation of cancer – analyses of cohorts of twins from Sweden, Denmark, and Finland.N. Engl. J. Med.34378–85. 10.1056/NEJM200007133430201
95
LittleT. J. (2002). The evolutionary significance of parasitism: do parasite-driven genetic dynamics occur ex silico?J. Evol. Biol.151–9. 10.1046/j.1420-9101.2002.00366.x
96
LittleT. J.WattK.EbertD. (2006). Parasite-host specificity: experimental studies on the basis of parasite adaptation.Evolution6031–38. 10.1111/j.0014-3820.2006.tb01079.x
97
LiuL.ChenX.SkogerboG.ZhangP.ChenR.HeS.et al (2012). The human microbiome: a hot spot of microbial horizontal gene transfer.Genomics100265–270. 10.1016/j.ygeno.2012.07.012
98
LivelyC. M.DybdahlM. F. (2000). Parasite adaptation to locally common host genotypes.Nature405679–681. 10.1038/35015069
99
MagnussonP. K.LichtensteinP.GyllenstenU. B. (2000). Heritability of cervical tumours.Int. J. Cancer88698–701. 10.1002/1097-0215(20001201)88:5<698::AID-IJC3>3.0.CO;2-J
100
MaidenM. C.BygravesJ. A.FeilE.MorelliG.RussellJ. E.UrwinR.et al (1998). Multilocus sequence typing: a portable approach to the identification of clones within populations of pathogenic microorganisms.Proc. Natl. Acad. Sci. U.S.A.953140–3145. 10.1073/pnas.95.6.3140
101
MalikA. N.Godfrey-FaussettP. (2005). Effects of genetic variability of Mycobacterium tuberculosis strains on the presentation of disease.Lancet Infect. Dis.5174–183. 10.1016/S1473-3099(05)01310-1
102
MayerleJ.den HoedC. M.SchurmannC.StolkL.HomuthG.PetersM. J.et al (2013a). Identification of genetic loci associated with Helicobacter pylori serologic status.JAMA3091912–1920. 10.1001/jama.2013.4350
103
MayerleJ.KuipersE. J.LerchM. M. (2013b). Genetic variants associated with susceptibility to Helicobacter pylori – reply.JAMA310976–977. 10.1001/jama.2013.194772
104
McFall-NgaiM.HadfieldM. G.BoschT. C.CareyH. V.Domazet-LosoT.DouglasA. E.et al (2013). Animals in a bacterial world, a new imperative for the life sciences.Proc. Natl. Acad. Sci. U.S.A.1103229–3236. 10.1073/pnas.1218525110
105
MessengerS. L.MolineuxI. J.BullJ. J. (1999). Virulence evolution in a virus obeys a trade-off.Proc. Biol. Sci.266397–404. 10.1098/rspb.1999.0651
106
MollerM.HoalE. G. (2010). Current findings, challenges and novel approaches in human genetic susceptibility to tuberculosis.Tuberculosis (Edinb.)9071–83. 10.1016/j.tube.2010.02.002
107
MonotM.HonoreN.GarnierT.ZidaneN.SherafiD.Paniz-MondolfiA.et al (2009). Comparative genomic and phylogeographic analysis of Mycobacterium leprae.Nat. Genet.411282–1289. 10.1038/ng.477
108
MoodleyY.LinzB. (2009). Helicobacter pylori sequences reflect past human migrations.Genome Dyn.662–74. 10.1159/000235763
109
MoodleyY.LinzB.BondR. P.NieuwoudtM.SoodyallH.SchlebuschC. M.et al (2012). Age of the association between Helicobacter pylori and man.PLoS Pathog.8:e1002693. 10.1371/journal.ppat.1002693
110
MoodleyY.LinzB.YamaokaY.WindsorH. M.BreurecS.WuJ. Y.et al (2009). The peopling of the Pacific from a bacterial perspective.Science323527–530. 10.1126/science.1166083
111
MungerK.BaldwinA.EdwardsK. M.HayakawaH.NguyenC. L.OwensM.et al (2004). Mechanisms of human papillomavirus-induced oncogenesis.J. Virol.7811451–11460. 10.1128/JVI.78.21.11451-11460.2004
112
NicolM. P.WilkinsonR. J. (2008). The clinical consequences of strain diversity in Mycobacterium tuberculosis.Trans. R. Soc. Trop. Med. Hyg.102955–965. 10.1016/j.trstmh.2008.03.025
113
OlesenR.WejseC.VelezD. R.BisseyeC.SodemannM.AabyP.et al (2007). DC-SIGN (CD209), pentraxin 3 and vitamin D receptor gene variants associate with pulmonary tuberculosis risk in West Africans.Genes Immun.8456–467. 10.1038/sj.gene.6364410
114
OngC. K.ChanS. Y.CampoM. S.FujinagaK.Mavromara-NazosP.LabropoulouV.et al (1993). Evolution of human papillomavirus type 18: an ancient phylogenetic root in Africa and intratype diversity reflect coevolution with human ethnic groups.J. Virol.676424–6431.
115
PallenM. J.WrenB. W. (2007). Bacterial pathogenomics.Nature449835–842. 10.1038/nature06248
116
ParwatiI.van CrevelR.van SoolingenD. (2010). Possible underlying mechanisms for successful emergence of the Mycobacterium tuberculosis Beijing genotype strains.Lancet Infect. Dis.10103–111. 10.1016/S1473-3099(09)70330-5
117
PeekR. M.Jr.BlaserM. J. (2002). Helicobacter pylori and gastrointestinal tract adenocarcinomas.Nat. Rev. Cancer228–37. 10.1038/nrc703
118
PicardC.CasanovaJ. L.AbelL (2006). Mendelian traits that confer predisposition or resistance to specific infections in humans.Curr. Opin. Immunol.18383–390. 10.1016/j.coi.2006.05.005
119
PlummerM.SchiffmanM.CastleP. E.Maucort-BoulchD.WheelerC. M.ALTS Group. (2007). A 2-year prospective study of human papillomavirus persistence among women with a cytological diagnosis of atypical squamous cells of undetermined significance or low-grade squamous intraepithelial lesion.J. Infect. Dis.1951582–1589. 10.1086/516784
120
PortaC.RiboldiE.SicaA. (2011). Mechanisms linking pathogens-associated inflammation and cancer.Cancer Lett.305250–262. 10.1016/j.canlet.2010.10.012
121
PortevinD.GagneuxS.ComasI.YoungD. (2011). Human macrophage responses to clinical isolates from the Mycobacterium tuberculosis complex discriminate between ancient and modern lineages.PLoS Pathog.7:e1001307. 10.1371/journal.ppat.1001307
122
RavelJ.GajerP.AbdoZ.SchneiderG. M.KoenigS. S.McCulleS. L.et al (2011). Vaginal microbiome of reproductive-age women.Proc. Natl. Acad. Sci. U.S.A. 108(Suppl. 1), 4680–4687. 10.1073/pnas.1002611107
123
RectorA.LemeyP.TachezyR.MostmansS.GhimS. J.Van DoorslaerK.et al (2007). Ancient papillomavirus–host co-speciation in Felidae.Genome Biol.8 R57. 10.1186/gb-2007-8-4-r57
124
ReedM. B.DomenechP.MancaC.SuH.BarczakA. K.KreiswirthB. N.et al (2004). A glycolipid of hypervirulent tuberculosis strains that inhibits the innate immune response.Nature43184–87. 10.1038/nature02837
125
ReilingN.HomolkaS.WalterK.BrandenburgJ.NiwinskiL.ErnstM.et al (2013). Clade-specific virulence patterns of Mycobacterium tuberculosis complex strains in human primary macrophages and aerogenically infected mice.MBio4pii:e00250-13. 10.1128/mBio.00250-13
126
RidenhourB. J.NuismerS. L. (2007). Polygenic traits and parasite local adaptation.Evolution61368–376. 10.1111/j.1558-5646.2007.00029.x
127
RothenbacherD.BlaserM. J.BodeG.BrennerH. (2000). Inverse relationship between gastric colonization of Helicobacter pylori and diarrheal illnesses in children: results of a population-based cross-sectional study.J. Infect. Dis.1821446–1449. 10.1086/315887
128
RothenbacherD.WinklerM.GonserT.AdlerG.BrennerH. (2002). Role of infected parents in transmission of Helicobacter pylori to their children.Pediatr. Infect. Dis. J.21674–679. 10.1097/00006454-200207000-00014
129
RowellJ. L.DowlingN. F.YuW.YesupriyaA.ZhangL.GwinnM. (2012). Trends in population-based studies of human genetics in infectious diseases.PLoS ONE7:e25431. 10.1371/journal.pone.0025431
130
SalamaN.GuilleminK.McDanielT. K.SherlockG.TompkinsL.FalkowS. (2000). A whole-genome microarray reveals genetic diversity among Helicobacter pylori strains.Proc. Natl. Acad. Sci. U.S.A.9714668–14673. 10.1073/pnas.97.26.14668
131
SalamaN. R.HartungM. L.MullerA. (2013). Life in the human stomach: persistence strategies of the bacterial pathogen Helicobacter pylori.Nat. Rev. Microbiol.11385–399. 10.1038/nrmicro3016
132
SamsonM.LibertF.DoranzB. J.RuckerJ.LiesnardC.FarberC. M.et al (1996). Resistance to HIV-1 infection in caucasian individuals bearing mutant alleles of the CCR-5 chemokine receptor gene.Nature382722–725. 10.1038/382722a0
133
SchiffmanM.HerreroR.DesalleR.HildesheimA.WacholderS.RodriguezA. C.et al (2005). The carcinogenicity of human papillomavirus types reflects viral evolution.Virology33776–84. 10.1016/j.virol.2005.04.002
134
SchneiderB. G.CamargoM. C.RyckmanK. K.SicinschiL. A.PiazueloM. B.ZabaletaJ.et al (2008). Cytokine polymorphisms and gastric cancer risk.Cancer Biol. Ther.7157–162. 10.4161/cbt.7.2.5270
135
ShahS. D.DoorbarJ.GoldsteinR. A. (2010). Analysis of host–parasite incongruence in papillomavirus evolution using importance sampling.Mol. Biol. Evol.271301–1314. 10.1093/molbev/msq015
136
SharmaS. K.MohanA. (2004). Extrapulmonary tuberculosis.Ind. J. Med. Res.120316–353.
137
ShiY.HuZ.WuC.DaiJ.LiH.DongJ.et al (2011). A genome-wide association study identifies new susceptibility loci for non-cardia gastric cancer at 3q13.31 and 5p13.1.Nat. Genet.431215–1218. 10.1038/ng.978
138
ShiY.LiL.HuZ.LiS.WangS.LiuJ.et al (2013). A genome-wide association study identifies two new cervical cancer susceptibility loci at 4q12 and 17q12.Nat. Genet.45918–922. 10.1038/ng.2687
139
SmillieC. S.SmithM. B.FriedmanJ.CorderoO. X.DavidL. A.Alm.E. J. (2011). Ecology drives a global network of gene exchange connecting the human microbiome.Nature480241–244. 10.1038/nature10571
140
SokurenkoE. V.GomulkiewiczR.DykhuizenD. E. (2006). Source-sink dynamics of virulence evolution.Nat. Rev. Microbiol.4548–555. 10.1038/nrmicro1446
141
StoreyA.ThomasM.KalitaA.HarwoodC.GardiolD.MantovaniF.et al (1998). Role of a p53 polymorphism in the development of human papillomavirus-associated cancer.Nature393229–234. 10.1038/30400
142
StuckiD.GagneuxS. (2013). Single nucleotide polymorphisms in Mycobacterium tuberculosis and the need for a curated database.Tuberculosis (Edinb.)9330–39. 10.1016/j.tube.2012.11.002
143
SuerbaumS.SmithJ. M.BapumiaK.MorelliG.SmithN. H.KunstmannE.et al (1998). Free recombination within Helicobacter pylori.Proc. Natl. Acad. Sci. U.S.A.9512619–12624. 10.1073/pnas.95.21.12619
144
SusserM.SteinZ. (2002). Civilization and peptic ulcer.Int. J. Epidemiol.3113–17. 10.1093/ije/31.1.13
145
The International HapMap Consortium. (2005). A haplotype map of the human genome.Nature4371299–1320. 10.1038/nature04226
146
ThompsonJ. N. (2014). Natural selection, coevolution, and the web of life.Am. Nat.183 iv–v. 10.1086/674238
147
ThompsonJ. N.NuismerS. L.GomulkiewiczR. (2002). Coevolution and maladaptation.Integr. Comp. Biol.42381–387. 10.1093/icb/42.2.381
148
ThyeT.NiemannS.WalterK.HomolkaS.IntemannC. D.ChinbuahA.et al (2011). Variant G57E of mannose binding lectin associated with protection against tuberculosis caused by Mycobacterium africanum but not by M. tuberculosis.PLoS ONE6:e20908. 10.1371/journal.pone.0020908
149
ThyeT.Owusu-DaboE.VannbergF. O.van CrevelR.CurtisJ.SahiratmadjaE.et al (2012). Common variants at 11p13 are associated with susceptibility to tuberculosis.Nat. Genet.44257–259. 10.1038/ng.1080
150
ThyeT.VannbergF. O.WongS. H.Owusu-DaboE.OseiI.GyapongJ.et al (2010). Genome-wide association analyses identifies a susceptibility locus for tuberculosis on chromosome 18q11.2.Nat. Genet.42739–741. 10.1038/ng.639
151
TiemersmaE. W.van der WerfM. J.BorgdorffM. W.WilliamsB. G.NagelkerkeN. J. (2011). Natural history of tuberculosis: duration and fatality of untreated pulmonary tuberculosis in HIV negative patients: a systematic review.PLoS ONE6:e17601. 10.1371/journal.pone.0017601
152
TorneselloM. L.DuraturoM. L.SalatielloI.BuonaguroL.LositoS.BottiG.et al (2004). Analysis of human papillomavirus type-16 variants in Italian women with cervical intraepithelial neoplasia and cervical cancer.J. Med. Virol.74117–126. 10.1002/jmv.20154
153
TorresJ.CorreaP.FerreccioC.Hernandez-SuarezG.HerreroR.Cavazza-PorroM.et al (2013). Gastric cancer incidence and mortality is associated with altitude in the mountainous regions of Pacific Latin America.Cancer Causes Control24249–256. 10.1007/s10552-012-0114-8
154
TurnbaughP. J.HamadyM.YatsunenkoT.Cantarel. B. L.Duncan. A.LeyR. E.et al (2009). A core gut microbiome in obese and lean twins.Nature457480–484. 10.1038/nature07540
155
VaeziM. F.FalkG. W.PeekR. M.VicariJ. J.GoldblumJ. R.Perez-PerezG. I.et al (2000). CagA-positive strains of Helicobacter pylori may protect against Barrett’s esophagus.Am. J. Gastroenterol.952206–2211. 10.1111/j.1572-0241.2000.02305.x
156
Van DoorslaerK. (2013). Evolution of the papillomaviridae.Virology44511–20. 10.1016/j.virol.2013.05.012
157
VelezD. R.HulmeW. F.MyersJ. L.StryjewskiM. E.AbbateE.EstevanR.et al (2009). Association of SLC11A1 with tuberculosis and interactions with NOS2A and TLR2 in African-Americans and Caucasians.Int. J. Tuberc. Lung Dis.131068–1076.
158
WheelerC. M. (2008). Natural history of human papillomavirus infections, cytologic and histologic abnormalities, and cancer.Obstet. Gynecol. Clin. North Am.35519–536; vii. 10.1016/j.ogc.2008.09.006
159
WhiteA. E.LivanosE. M.TlstyT. D. (1994). Differential disruption of genomic integrity and cell cycle regulation in normal human fibroblasts by the HPV oncoproteins.Genes Dev.8666–677. 10.1101/gad.8.6.666
160
WilfertL.Schmid-HempelP. (2008). The genetic architecture of susceptibility to parasites.BMC Evol. Biol.8:187. 10.1186/1471-2148-8-187
161
WirthT.WangX.LinzB.NovickR. P.LumJ. K.BlaserM.et al (2004). Distinguishing human ethnic groups by means of sequences from Helicobacter pylori: lessons from Ladakh.Proc. Natl. Acad. Sci. U.S.A.1014746–4751. 10.1073/pnas.0306629101
162
WoolhouseM. E.WebsterJ. P.DomingoE.CharlesworthB.LevinB. R. (2002). Biological and biomedical implications of the co-evolution of pathogens and their hosts.Nat. Genet.32569–577. 10.1038/ng1202-569
163
WroblewskiL. E.PeekR. M.Jr.WilsonK. T. (2010). Helicobacter pylori and gastric cancer: factors that modulate disease risk.Clin. Microbiol. Rev.23713–739. 10.1128/CMR.00011-10
164
XiL. F.KoutskyL. A.GallowayD. A.KuypersJ.HughesJ. P.WheelerM.et al (1997). Genomic variation of human papillomavirus type 16 and risk for high grade cervical intraepithelial neoplasia.J. Natl. Cancer Inst.89796–802. 10.1093/jnci/89.11.796
165
zur HausenH. (1989). Papillomavirus in anogenital cancer: the dilemma of epidemiologic approaches.J. Natl. Cancer Inst.811680–1682. 10.1093/jnci/81.22.1680
166
zur HausenH. (1991). Viruses in human cancers.Science2541167–1173. 10.1126/science.1659743
Summary
Keywords
host–pathogen co-evolution, human disease, Helicobacter pylori, Mycobacterium tuberculosis, human papillomavirus, genome–genome interactions
Citation
Kodaman N, Sobota RS, Mera R, Schneider BG and Williams SM (2014) Disrupted human–pathogen co-evolution: a model for disease. Front. Genet. 5:290. doi: 10.3389/fgene.2014.00290
Received
06 June 2014
Accepted
05 August 2014
Published
25 August 2014
Volume
5 - 2014
Edited by
José M. Álvarez-Castro, Universidade de Santiago de Compostela, Spain
Reviewed by
Dana Crawford, Vanderbilt University, USA; Sebastien Gagneux, Swiss Tropical and Public Health Institute, Switzerland
Copyright
© 2014 Kodaman, Sobota, Mera, Schneider and Williams.
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) or licensor are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Scott M. Williams, Department of Genetics, Geisel School of Medicine, Dartmouth College, HB-6044, Hanover, NH 03755, USA e-mail: scott.williams@dartmouth.edu
This article was submitted to Evolutionary and Population Genetics a section of the journal Frontiers in Genetics.
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.